Aggregate Interference Modeling in Cognitive Radio Networks with Power and Contention Control
In this paper, we present an interference model for cognitive radio (CR) networks employing power control, contention control or hybrid power/contention control schemes. For the first case, a power control scheme is proposed to govern the transmission power of a CR node.
For the second one, a contention control scheme at the media access control (MAC) layer, based on carrier sense multiple access with collision avoidance (CSMA/CA), is proposed to coordinate the operation of CR nodes with transmission requests.
The probability density functions of the interference received at a primary receiver from a CR network are first derived numerically for these two cases. For the hybrid case, where power and contention controls are jointly adopted by a CR node to govern its transmission, the interference is analyzed and compared with that of the first two schemes by simulations. Then, the interference distributions under the first two control schemes are fitted by log-normal distributions with greatly reduced complexity. Moreover, the effect of a hidden primary receiver on the interference experienced at the receiver is investigated. It is demonstrated that both power and contention controls are effective approaches to alleviate the interference caused by CR networks. Some in-depth analysis of the impact of key parameters on the interference of CR networks is given via numerical studies as well.
Index Terms – Cognitive radio, interference modeling, hidden primary receiver.
With the requirement to improve spectrum utilization, the newly emerging cognitive radio (CR) technology – has attracted increasing attention. A CR network is envisioned to be capable of reusing the unused or underutilized spectra of incumbent systems (also known as primary networks) by sensing its surrounding environment and adapting its operational parameters autonomously. A CR system may coexist with a primary network on either an interference-free or interference-tolerant basis . For the former case, the CR system only exploits the unused spectra of the primary network, which consequently guarantees no interference to primary users. For the latter case, the CR system is allowed to share the spectra assigned to the primary network, under the condition that the CR network must not impose detrimental interference on the primary network. Therefore, modeling and analyzing the interference caused by CR networks is of great importance to reveal how the service of a primary network is deteriorated and how CR networks may be deployed.
In the literature, the existing research on interference modeling for CR networks mainly falls into three categories: spatial, frequency-domain and accumulated interference modeling. For spatial interference modeling, the fraction of white spaces available for CR networks was investigated in  and . In , the region of interference for CR receivers and region of communication for CR transmitters were studied for the case where a CR network coexists with a cellular network. The interference from CR devices to wireless microphones operating in TV bands was analyzed in , where the loss of reliable communication area of a wireless microphone due to the existence of CR devices was examined. CR interference in the frequency domain was also researched in the literature, e.g., the interference due to out-of-band emission of a wireless regional area network (WRAN) was analyzed in .
As for accumulated interference modeling, in , the aggregate interference power from a sea of CR transmitters surrounding a primary receiver was derived. Also, the accumulated CR transmission power perceived at a primary receiver was given by integrating over the “CR sea” with a certain power density. The performance of a primary system was evaluated in  in terms of outage probability caused by the interference from CR networks. The outage probability was derived for both underlay and overlay spectrum sharing cases. In  the aggregate interference from multiple CR transmitters following a Poisson point process was approximated by a Gamma distribution and the probability of interference at a primary receiver was also given. It is worth noting that only pathloss was assumed for the interfering channel in –. Their work was extended by taking both shadowing and fading into account in  and . Moreover, the probability density function (PDF) for accumulated interference and outage probability due to the aggregate interference from CR nodes were also derived in  and , respectively.
However, in all the previous works –, the CR transmitters were assumed to transmit at a fixed power level, i.e., no power control for CR transmitters was considered. Moreover, the CR nodes were all assumed to communicate with each other simultaneously. Thus, no contention control scheme was employed at the cognitive media access control (MAC) layer. Some preliminary results on CR interference modeling were obtained in  by incorporating either power or contention control scheme. In this paper, we extend the aggregate interference modeling in the following aspects. Firstly, a more realistic power control scheme than that in  is proposed, and a hybrid power/contention control scheme is introduced. Secondly, the PDFs of interference perceived at a primary network from a CR network are derived numerically for the cases of power or contention control. The interference distribution of the hybrid control scheme is also analyzed and compared with that of the pure power control and pure contention control schemes by simulations. Furthermore, for the power and contention control schemes, their interference distributions are fitted by log-normal distributions, which greatly reduces computational complexity compared to a numerical approach to obtain PDFs. Finally, the impact of a hidden primary receiver on the aggregate interference is investigated for all the three schemes. The impact of several key parameters on the resulting interference is evaluated as well, which provides some insights for the deployment of CR networks.
The remainder of this paper is organized as follows. The system model is elaborated in Section II. The detailed interference modeling is presented in Section III. In Section IV, the interference distributions are approximated by log-normal distributions. We incorporate the hidden primary receiver problem in Section V. The impact of several key parameters on the interference is analyzed via numerical studies in Section VI. Finally, Section VII concludes the paper.
Ii System Model
The system model is illustrated in Fig. 1. It consists of a CR network coexisting with a primary transmitter-receiver pair. The active CR transmitters are distributed in a 2-dimensional plane outside the interference region (IR) of the primary receiver as shown in Fig. 1. The IR is a disk centered at the primary receiver with a radius . CR transmission is forbidden within this circular IR in order to protect the primary receiver against co-channel interference from the surrounding CR transmitters, since it is assumed that all the CR transmitters reside in the same frequency spectrum as the primary transmitter. We model the aggregate interference received at the primary receiver due to the existence of a CR network and investigate the impact of CR network deployment parameters on the resulting aggregate interference.
The underlying interference channels from CR transmitters to the primary receiver experience pathloss, shadowing and fading. The pathloss function is
where is the distance between the th active CR transmitter and the primary receiver and is the pathloss exponent. The composite model for shadowing and fading can be expressed as the product of the long term shadowing and the short term multipath fading. In this paper, log-normal shadowing and Nakagami fading are considered. Let denote the channel gain for the composite shadowing and fading of the interference channel from the th active CR transmitter to the primary receiver. The PDF of the composite channel gain can be approximated by the following log-normal distribution 
where the mean and variance can be expressed as
with standing for the Nakagami shape factor and and denoting the standard mean and variance of the log-normal distribution, respectively.
Let denote the transmission power of the th active CR transmitter. The accumulated power of the instantaneous interference received at the primary receiver can be expressed as
In this paper, we investigate the characteristics of the aggregate interference from all CR transmitters employing the following three different schemes: (i) power control, (ii) contention control, and (iii) hybrid power/contention control.
Ii-a Power Control
In this scenario, the distribution of active CR transmitters follows a Poisson point process with a density parameter for the density of CR transmitters on the plane.
The transmission power of a CR transmitter is governed by the following power control law
where is the distance from the th active CR transmitter to its nearest neighbouring active CR transmitter, is the power control exponent, is the maximum transmission power for CR transmitters, and is the power control range, which determines the minimum leading to maximum CR transmission power . Compared to the power control law in , a new parameter is introduced here to adjust the range of the power control. We assume that the power control exponent is equal to the pathloss exponent in (1) throughout the paper. The above proposed power control scheme is designed in such a manner that the interference caused by the th active CR transmitter to its nearest active CR transmitter due to pathloss is . It is clear that within the power control range , this interference is equal to a constant . But beyond the power control range, the interference is less than that constant. In other words, at any CR transmitter the interference from the nearest neighbouring CR transmitter is capped and independent of the nearest neighbour distance within the power control range. It is worth noting that for each CR transmitter the information of its nearest neighbour distance is indispensable to determine its transmission power. Therefore, to facilitate the abovementioned power control scheme, either a central console having the global position information of all active CR transmitters or a distributed sensing scheme for CR transmitters like pilot sensing  is required, the detail of which is, however, beyond the scope of this work. When CR transmitters follow a Poisson point distribution with a density , the PDF of can be given as 
Ii-B Contention Control
Unlike the previous power control scheme, for the case of contention control every active CR transmitter has fixed transmission power , but their transmission is governed by contention control to determine which CR transmitters can transmit at a given time. We assume that the multiple access protocol carrier sense multiple access with collision avoidance (CSMA/CA) is employed, like in IEEE 802.11 networks. Every CR transmitter senses the medium before transmission. If the medium is busy, namely, the CR transmitter detects transmission from other CR transmitters within its contention region, it defers its transmission. Otherwise, the CR transmitter starts its transmission. As a result of the contention control shown, all the active CR transmitters are separated from each other by at least the contention distance, which is the minimum distance between two concurrent CR transmitters.
The distribution of the active CR transmitters under the contention control can be modeled as a Matern-hardcore (MH) point process , which can be considered as a thinned process from a Poisson point process . The thinning operation deletes some points from the original Poisson process under certain criteria. The MH process is the result of dependent thinning from a Poisson point process , i.e., deleting or retaining a point depends on previous deletion operations. The mathematical expression of the MH process is given by 
Each point in the original Poisson point process is marked with a random variable uniformly distributed in (0,1), while is a disk centered at point with the radius . The retaining probability for the MH process, which is the probability of a point from a Poisson point process with a density surviving the thinning process, is given by 
Ii-C Hybrid Power/Contention Control
The aforementioned power control scheme regulates the transmission power of each CR transmitter according to its nearest neighbouring transmitter distance, while the contention control determines which CR transmitter can transmit at a time instant with fixed transmission power. A natural extension of the above two interference management schemes is to implement both schemes in the same system. This is termed hybrid power/contention control and it works in the following manner. The contention control scheme is first applied, resulting in a set of active CR transmitters following an MH point process. Then, a power control scheme similar to (6) is employed to adjust the transmission power of each active CR transmitter according to the distance to the nearest neighbouring active transmitter. The following power control law is adopted in the hybrid control scheme
where is the distance from an active CR transmitter to its nearest neighbouring active CR transmitter, is the power control exponent as in (6), and is the power control range similar to in (6) except that it also determines the maximum transmission power, i.e., . It is obvious that a larger leads to a larger maximum CR transmission power and, consequently, longer communication range for CR transmitters. The above power control law (10) guarantees that when a pathloss channel is considered for each active CR transmitter, the perceived interference caused by its nearest neighbouring CR transmitter is , which is (i) a constant within the power control range and (ii) less than the constant when the distance is larger than the power control range.
Iii Interference Modeling
We intend to model the aggregate interference from CR transmitters employing the three different interference management schemes introduced in Section II by finding their corresponding PDFs. We apply the methodology used, for example, in  and  to derive the PDFs. First, the characteristic functions of the interference under different system models are derived. Then, the PDFs of the aggregate interference are obtained by performing an inverse Fourier transform on their characteristic functions.
Iii-a Power Control
When all the CR transmitters follow a Poisson point process distribution and employ the power control scheme proposed in (6), we can adopt the characteristic function-based method as in , - and obtain the following characteristic function of the aggregate interference at a primary receiver from all CR transmitters
where is the PDF of the transmission power of a CR transmitter defined in (6) and
In (12), denotes the inverse function of in (1). For the derivation of (11), the following fact is used: the distances from the th CR transmitter to the primary receiver have independent and identical uniform distributions for a given number of CR transmitters . Their PDFs have the following form 
when CR transmitters are distributed within an annular ring with inner radius and outer radius . In (11), is a function of as shown in (6), so the expectation of over equals that of over . Using the PDF of given in (7), (11) can be rewritten as
Moreover, (14) can be written as (see Appendix A for the detailed derivation procedure)
Finally, we obtain the PDF of the interference by performing the inverse Fourier transform on as
Equations (III-A) and (16) serve as general expressions for the characteristic function and PDF, respectively, of the interference under the power control scheme. As a special case, when the pathloss exponent and the radius of the interference region , the PDF can be further simplified through similar steps to that used in  and obtained as
The detailed derivation procedure for can be found in Appendix B.
Iii-B Contention Control
As mentioned in Section II.B, the distribution of CR transmitters can be modeled as an MH point process when the contention control is adopted. The MH process is a dependent thinning process from the original Poisson point process, which means that the positions of CR transmitters are correlated to each other. However, it is very difficult to obtain the distribution function like (13) for an MH point process in order to model the distance from an active CR transmitter to the primary receiver. Instead, we approximate the MH point process as an independent thinned Poisson point process with retaining probability given by (9). Then, the transmission power for the th CR transmitter is , which is a random variable taking values or with probabilities and , respectively. To this end, the contention control scheme can be interpreted as follows: all the CR transmitters still follow the original Poisson point process with intensity , but the th CR transmitter has probability to transmit at power level . The characteristic function of the accumulated interference can be found as
The detailed derivation of (19) is presented in Appendix C.
It is worth noting that the approximation for the MH point process actually ignores the dependence among the CR transmitters and treats an MH point process as a result of independent thinning process from an original Poisson point process. The accuracy of this approximation is evaluated in Section IV.
Iii-C Hybrid Power/Contention Control
So far, the PDFs of the interferences received at a primary receiver from a CR network employing power control and contention control schemes have been derived. In order to model the aggregate interference under the hybrid control scheme, the nearest neighbouring distance distribution function analogous to (7) for active CR transmitters is indispensable to evaluate the transmission power designated in (10). Unfortunately, there is no closed-form expression for the nearest neighbour distance distribution function for an MH point process . Alternatively, several estimators have been used to statistically estimate the nearest neighbour distance distribution function in practice . However, statistical estimation is not practical for deriving the characteristic function in our case. Thus, we approach this problem numerically.
The PDF for the aggregate interference under the hybrid control scheme is simulated in Fig. 2, where the interference PDFs for power and contention control are given as well for the purpose of comparison. It can be seen from this figure that both the mean and variance of the aggregate interference increase for the hybrid control scheme compared to either power or contention control schemes. However, the boosted interference is paid off by the increased CR communication area (coverage) for the hybrid control scheme. We define the coverage of each CR transmitter as a circular disk centered at a CR transmitter with radii being , and for power control, contention control and hybrid power/contention control schemes, respectively. Then, the received signal power at cell edge of a CR transmitter due to pathloss is , and for the above three aforementioned schemes. For the sake of comparison, let and , which guarantees that the strength of the received signal power at cell edge of a CR transmitter is the same for all the three schemes. The overall coverage of the CR netwrok under different control schemes can be investigated numerically. With this setup, the overall coverage ratio for the power control, contention control and hybrid power/contention control is , respectively. Two interesting facts are unveiled from this experiment. Firstly, the power control scheme leads to slightly smaller interference and slightly lager coverage compared to the contention control scheme, which suggests that power control is preferable to contention control in terms of lower resulting interference and larger coverage if the CR system can afford the complexity introduced by implementing the power control scheme. Secondly, the hybrid scheme tends to cause higher interference, but it greatly enlarges the coverage compared to power and contention control schemes.
Iv Analytical Approximation
In the previous section, to derive the PDFs for aggregate interference, the characteristic function-based method has been used which consists of two steps. Namely, characteristic function computation and Fourier transformation. This interference modeling approach is extremely computation-intensive, since generally closed-form expressions are not admitted for either step and the computations in both steps have to be performed numerically. It is desirable to model the aggregate interference with less complexity. An alternative approach to model the interference, which greatly reduces complexity, is to approximate interference PDFs as certain known distributions. Observations from Fig. 2 suggest that the interference distribution for either power or contention control is positively skewed and heavy-tailed, which suggests a log-normal distribution. Thus, in this section, we fit the aggregate interference under power and contention control schemes to log-normal distributions. The theory behind the log-normal fitting is based on the following two facts. It has been shown that the sum of interference from uniformly distributed interferers in a circular area is asymptotically log-normal , . This ensures that the aggregate interference in these two schemes can be approximated as log-normal distributed. Meanwhile, the sum of randomly weighted log-normal variables can be modelled as a log-normal distribution as well , which guarantees that the aggregate interference is still log-normal distributed even if the effect of shadow fading (2) is taken into account. In what follows, the log-normal fitting is performed using a cumulant-matching approach , where the first two order cumulants of the aggregate interference in (5) are used to estimate the mean and variance of the log-normal distribution function. Therefore, the exact PDFs of interference can be obtained. Fortunately, these cumulants have closed-form expressions for both control schemes. Consequently, it significantly reduces the complexity compared to the interference modeling carried out in Section III.
For the PDF of a log-normal variable
its mean and variance can be estimated using its first two order cumulants and as follows :
In the context of interference distribution fitting, the th cumulant of the aggregate interference can be obtained from its characteristic function via the following equation
Iv-a Power Control
To evaluate the accuracy of the log-normal approximation for the power control case, some comparisons are performed in Fig. 3(a). It can be seen from Fig. 3(a) that there is fairly good agreement between the interference PDFs derived in Section III and the approximated counterparts. This approximation approach can be applied to both the pathloss-only and shadow fading channels.
Iv-B Contention Control
Following the similar steps as in Appendix D and given the characteristic function (19) for the aggregate interference under contention control and also using (24), we can find the th cumulant of aggregate interference as
The accuracy evaluation of log-normal approximation under the contention control scheme is also performed and shown in Fig. 3(b). It can be seen from this figure that the log-normal approximation is fairly accurate compared to the simulated interference PDFs for either pathloss-only or shadow fading channels. Moreover, the derived interference PDF obtained from (16) and (19) is validated against the simulated counterpart in Fig. 3(b) as well, which suggests that the approximation for the MH point process in the derivation is reasonable.
V Imperfect Primary System Knowledge
In practice, some information about the primary system may not be perfectly known. One prominent example is the location of the primary receivers, which is usually required by CR networks in order to protect primary receivers from interfering CR transmitters. However, this information is not always available, especially in the case of passive primary receivers, i.e., when the primary receivers are hidden from CR networks. It is widely accepted that passive receiver detection techniques can be used or developed in the context of CR networks. For example, one of such primary receiver detection techniques is reported in . Nevertheless, its applicability is still not convincingly viable since it requires deploying sensor nodes close to primary receivers and much coordination is involved between these sensors and CR networks as well. The most commonly used and also the simplest approach to protect the primary receiver is to regulate the transmission of the CR network based on primary transmitter sensing, assuming that primary receivers are in close proximity to the primary transmitter. In this section, we evaluate the effect of a hidden primary receiver on the resulting interference to primary receivers.
Consider a primary and CR coexisting systems depicted in Fig. 4, where an IR with radius centered at the primary transmitter is introduced. All CR transmitters are distributed in the shaded concentric ring with inner radius and outer radius . Let be the angle between the line joining the primary receiver and a CR transmitter and the line joining the primary transmitter-receiver pair. The distance from the CR transmitter to the primary transmitter is and the distance between the primary transmitter-receiver pair is . Then, the distance between the CR transmitter and the primary receiver can be expressed as
where is distributed as in (13) and is uniformly distributed in if a Poisson point process is assumed for the CR transmitter distribution.
V-a Power Control
Under the power control scheme proposed in Section II.A and the system model given in Fig. 4, the characteristic function of aggregate interference can be written as follows (see Appendix E for the detailed derivation):
Applying the log-normal approximation method used in Section IV, we obtain the th cumulant of interference as
As can be seen from (V-A), unlike (IV-A), the th cumulant does not have a closed-form expression. However, the complexity of obtaining the exact interference PDF from (V-A) is still smaller than that of the numerical method in Section III.
An experiment is performed in Fig. 5(a) to examine the effect of hidden primary receiver on the resulting interference compared to the interference for the case of perfect knowledge of primary receiver location. We consider a pathloss-only channel in this figure. It can be seen from the figure that the hidden primary receiver problem boosts the interference in terms of increased interference mean and variance. This figure also shows that the log-normal approximation still fits well the interference distribution in this scenario.
V-B Contention Control
Under the contention control scheme proposed in Section II.B and the system model given in Fig. 4, the characteristic function of aggregate interference can be expressed as
Using the same log-normal approximation method as in Section IV, the th cumulant of interference can be written as
The effect of hidden primary receiver under contention control is evaluated in Fig. 5(b), where a pathloss-only channel is assumed. As we can see from this figure, the uncertainty about the primary receiver location leads to interference with larger mean and variance as compared to that in the case with perfect knowledge of primary receiver location. Moreover, it can be seen from this figure that the log-normal fitting for the interference is fairly accurate and the approximation approach is still applicable in this scenario.
For the case of hybrid power/contention control, the effect of hidden priamry receiver cannot be examined analytically because the closed-form interference PDF is not available. Therefore, it is analyzed numerically in Fig. 6, whose initial setup is the same as the one used in Fig. 5(b) except that the power control range is m. It can bee seen from Fig. 6 that the uncertainty about the primary receiver location boosts the interference in terms of increased mean, variance, and heavier tails for the hybrid control scheme as well. More interestingly, another two facts can be found by comparing Figs. 5 and 6: (i) the hidden primary receiver phenomonon has similar impact on the pure power and pure contention control schemes; (ii) the hybrid power/contention control scheme is less sensible to the phenomenon of hidden primary receiver than any of the other two schemes.
Vi Numerical Studies
The aggregate interference power from CR transmitters employing power control or contention control is investigated numerically in this section. For the power control scheme, Fig. 7(a) shows the effect of different power control parameters on their resulting aggregate interference. The detailed setup for the initial power control scheme is as follows: the maximum transmission power for each CR transmitter W, the density of CR transmitter user/10m, the IR radius m, the power control range m, the pathloss exponent and the power control exponent . From the two rightmost PDFs in this figure, it can seen that introducing power control scheme actually shifts the interference distribution leftwards compared to the distribution without power control. It means that the power control scheme can reduce the interference experienced at the primary receiver in terms of reducing its mean and slightly decreasing its variance. When deploying a CR network under the power control scheme, its resulting interference can be controlled by manipulating several parameters including , , , and . It can be seen in Fig. 7(a) that the interference can be reduced by either decreasing the maximum transmission power and/or CR density, or increasing the power control range and/or IR radius. Interestingly, it also suggests that adjusting the IR radius is an effective way to control the interference, since the interference is more sensitive to the IR radius than to any other parameter as demonstrated in Fig. 7(a). Meanwhile, the interference is least sensitive to the CR user density in the sense that halving leads to higher interference compared to doubling , halving or doubling .
For the contention control scheme, the impact of contention control parameters on the resulting interference is depicted in Fig. 7(b), whose initial setup is the same as that of Fig. 7(a) except that the transmission power for each CR transmitter is W and the contention control range is m. It can be seen from the two rightmost PDFs in Fig. 7(b) that the contention control scheme results in an interference distribution with reduced mean like the power control scheme in Fig. 7(a). Meanwhile, the interference can be reduced by decreasing , , and/or increasing or . It can be observed by comparing Fig. 7(b) with Fig. 7(a) that (i) increasing the IR radius is an effective approach to reduce the interference for both the power and contention control schemes. However, the power control scheme is more sensitive to the IR radius than the contention control one; (ii) reducing the transmission power and/or CR transmitter density affects the interference in the very similar manner for these two control schemes.
Finally, the impact of shadow fading on the aggregate interference is investigated for different values of the Nakagami shaping factor under power and contention control schemes, respectively, in Figs. 8(a) and 8(b). The initial setup in this example is the same as the one used for Figs. 7(a) and 7(b), except that the standard variance is dB. When the interfering channel becomes a Rayleigh channel, which is dominated by the log-normal shadowing. Whereas, when the fluctuations of the channel are reduced significantly compared to the Rayleigh fading channel. One fact observed in Fig. 8 is that the interference distributions have larger variance and heavier tails when shadow fading is incorporated for both control schemes. Interestingly, fading tends to make the interference distribution more heavy-tailed than shadowing, i.e., the interference under shadowing has better outage property than that under fading. Moreover, the shadow fading has the similar effect for both control schemes.
Interference at a primary receiver caused by CR transmitters with power control, contention control, and hybrid power/contention control schemes has been characterized. The PDFs of interference in the first two cases have been evaluated analytically while, the interference distribution under the hybrid power/contention control has been studied numerically. It has been found that the proposed power control and contention control schemes are two effective approaches to alleviate interference caused by CR transmitters. The hybrid control scheme causes higher interference to a primary receiver, but leads to larger CR coverage as compared to either power or contention control schemes. Then, the interference distributions for power and contention control schemes have been approximated by log-normal distributions using the cumulant-matching approach where the interference PDFs have been obtained with reduced complexity. Furthermore, the effect of a hidden primary receiver on the perceived interference has also been investigated for the primary receiver. Numerical studies have demonstrated the impact of some CR deployment parameters on the resulting aggregate interference under power and contention control schemes. It has been shown that increasing the IR radius is an effective way to reduce the interference. Moreover, the power control scheme is more sensitive to the IR radius than the contention control counterpart. Finally, the impact of shadow fading on the aggregate interference has been analyzed as well.
Z. Chen, C.-X. Wang, X. Hong, and J. Thompson acknowledge the support from the Scottish Funding Council for the Joint Research Institute in Signal and Image Processing between the University of Edinburgh and Heriot-Watt University, as a part of the Edinburgh Research Partnership in Engineering and Mathematics (ERPem). S. A. Vorobyov acknowledges the support in part from the Natural Sciences and Engineering Research Council (NSERC) of Canada and in part from the Alberta Ingenuity Foundation, Alberta, Canada. X. Ge acknowledges the support from National Natural Science Foundation of China (NSFC) (Grant No.: 60872007), National 863 High Technology Program of China (Grant No.: 2009AA01Z239), and the Ministry of Science and Technology (MOST) of China, International Science and Technology Collaboration Program (Grant No.: 0903). F. Zhao acknowledges the support from the NSFC (Grant No.: 60872022). C.-X. Wang and F. Zhao acknowledge the support of the Key Laboratory of Cognitive Radio and Information Processing (Guilin University of Electronic Technology), Ministry of Education, China. The authors acknowledge the support from the RCUK for the UK-China Science Bridges Project: R&D on (B)4G Wireless Mobile Communications.
A. Derivation of (III-A)
B. Derivation of (18)
C. Derivation of (19)
Following similar steps as in , the characteristic function of the aggregate interference can be expressed as
The integral in the last equality of (VII) can be written as
D. Derivation of (IV-A)
The first equality of (VII) is obtained based on the following fact
-  S. Haykin, “Cognitive radio: brain-empowered wireless communications,” IEEE J. Sel. Areas Commun., vol. 23, no. 2, pp. 201–220, Sept. 2005.
-  I. F. Akyildiz, W. Y. Lee, M. C. Vuran, and S. Mohanty, “NeXt generation/dynamic spectrum access/cognitive radio wireless networks: A survey,” Computer Networks, vol. 50, no. 13, pp. 2127–2159, Sept. 2006.
-  Q. Zhao and B. M. Sadler, “A survey of dynamic spectrum access,” IEEE Signal Process. Mag., vol. 24, no. 3, pp. 79–89, May 2007.
-  C.-X. Wang, H.-H. Chen, X. Hong, and M. Guizani, “Cognitive radio network management: tuning in to real-time conditions,” IEEE Veh. Technol. Mag., vol. 3, no. 1, pp. 28–35, Mar. 2008.
-  C.-X. Wang, X. Hong, H.-H. Chen, and J. S. Thompson, “On capacity of cognitive radio networks with average interference power constraints,” IEEE Trans. Wireless Commun., vol. 8, no. 4, pp. 1620–1625, Apr. 2009.
-  X. Hong, C.-X. Wang, J. Thompson, and Y. Zhang, “Demystifying white spaces,” in Proc. IEEE ICCCAS’08, Xiamen, China, May 2008, pp. 350–354.
-  X. Hong, C.-X. Wang, H.-H. Chen, and Y. Zhang, “Secondary spectrum access networks: recent development on the spatial models,” IEEE Veh. Technol. Mag., vol. 4, no. 2, pp. 36–43, June 2009.
-  T. Kamakaris, D. Kivanc-Tureli, and U. Tureli, “Interference model for cognitive coexistence in cellular systems,” in Proc. IEEE GLOBECOM’07, Washington, DC, USA, Nov. 2007, pp. 4175–4179.
-  R. S. Dhillon and T. X. Brown, “Models for analyzing cognitive radio interference to wireless microphones in TV bands,” in Proc. IEEE DySPAN’08, Chicago, USA, Oct. 2008, pp. 1–10.
-  G. L. Stuber, S. Almalfouh, and D. Sale, “Interference analysis of TV band whitespace,” in Proc. IEEE, vol. 97, no. 4, pp. 741–754, Apr. 2009.
-  N. Hoven and A. Sahai, “Power scaling for cognitive radio,” in Proc. IEEE WNCMC’05, Hawaii, USA, June 2005, pp. 250–255.
-  R. Menon, R. M. Buehrer and J. Reed, “Outage probability based comparison of underlay and overlay spectrum sharing techniques,” in Proc. IEEE DySPAN’05, Baltimore, USA, Nov. 2005, pp. 101–109.
-  M. Timmers, S. Pollin, A. Dejonghe, A. Bahai, L. Van der Perre, and F. Catthoor, “Accumulative interference modeling for cognitive radios with distributed channel access,” in Proc. IEEE CrownCom’08, Singapore, May 2008.
-  X. Hong, C.-X. Wang, and J. S. Thompson, “Interference modeling of cognitive radio networks,” in Proc. IEEE VTC’08-Spring, Singapore, May 2008, pp. 1851–1855.
-  R. Menon, R. Buehrer, and J. Reed, “On the impact of dynamic spectrum sharing techniques on legacy radio systems,” IEEE Trans. Wireless Commun., vol. 7, no. 11, pp. 4198–4207, Nov. 2008.
-  Z. Chen, C.-X. Wang, X. Hong, J. Thompson, S. A. Vorobyov and X. Ge, “Interference modeling for cognitive radio networks with power or contention control,” in Proc. IEEE WCNC 2010, Sydney, Australia, Apr. 2010.
-  G. L. Stuber, Principles of Mobile Communication, 2nd Edition, Boston: Kluwer Academic Publishers, 2001.
-  H. Q. Nguyen, F. Baccelli and D. Kofman, “A stochastic geometry analysis of dense IEEE 802.11 networks,” in Proc. IEEE INFOCOM’07, Anchorage, USA, May 2007, pp. 1199–1207.
-  M. Ghosh, “Text on FFT-based pilot sensing,” IEEE 802.22, doc. no. 22-07-0298-01-0000, July 2007.
-  D. Stoyan, W. S. Kendall, and J. Mecke, Stochastic Geometry and Its Applications, Chichester: John Wiley & Sons, 1986.
-  E. S. Sousa and J. A. Silvester, “Optimum transmission range in a direct-sequence spread-spectrum multihop pack radio network,” IEEE J. Sel. Areas Commun., vol. 8, no. 5, pp. 762–771, June 1990.
-  X. Yang and A. P. Pertropulu, “Co-channel interference modeling and analysis in a Poisson field of interferers in wireless communications,” IEEE Trans. Signal Process., vol. 51, no. 1, pp. 63–76, Jan. 2003.
-  P. C. Pinto and M. Z. Win, “Communication in a poisson field of interferers,” in Proc. IEEE 40th Annual Conf. Inform. Sciences and Systems, Princeton, USA, Mar. 2006, pp. 432–437.
-  J. E. Paloheimo, “On a theory of search,” Biometrika, vol. 58, no. 1, pp. 61–75, Apr. 1971.
-  D. Stoyan, “On estimators of the nearest neighbour distance distribution function for stationary point processes,” Metrika, vol. 64, no. 2, pp. 139-150, Feb. 2006.
-  J. Salo, L. Vuokko, H. M. El-Sallabi, and P. Vainikainen, “An additive model as a physical basis for shadow fading,” IEEE Trans. Veh. Technol., vol. 56, no. 1, pp. 13–26, Jan. 2007.
-  M. Pratesi, F. Santucci, and F. Graziosi, “Generalized moment matching for the linear combination of log-normal RVs: application to outage analysis in wireless systems,” IEEE Trans. Wireless Commun., vol. 5, no. 5, pp. 1122–1132, May 2006.
-  C. C. Chan and S. V. Hanly, “Calculating the outage probability in a CDMA network with spatial poisson traffic,” IEEE Trans. Veh. Technol., vol. 50, no. 1, pp. 183–204, Jan. 2001.
-  R. Menon, R. M. Buehrer, and J. H. Reed, “Impact of exclusion region and spreading in spectrum-sharing ad hoc networks,” in Proc. 1st Int. Workshop on Technology and Policy for Accessing Spectrum, TAPAS06, Aug. 2006.
-  J. Aitchison and J. A. C. Brown, The Lognormal Distribution, Cambridge University Press, Cambridge UK, 1957.
-  B. Wild and K. Ramchandran, “Detecting primary receivers for cognitive radio applications,” in Proc. IEEE DySPAN 2005, Nov. 2005, pp. 124–130.
Fig. 1. System model for CR networks coexisting with a primary network ( m).
Fig. 2. Comparison of interference distributions for power, contention and hybrid power/contention control schemes (100 m, 3 user/m, 4, m, , W, W, m and ).
Fig. 3. Log-normal approximation for interference distribution under (a) power control (100 m, 3 user/m, 4, m, , W, and dB) or (b) contention control (100 m, 3 user/m, 4, m, W, and dB).
Fig. 4. Imperfect knowledge of primary receiver location - the primary receiver is hidden from all CR transmitters distributed in the shaded region.
Fig. 5. Log-normal approximation for interference distribution with a hidden primary receiver under (a) power control (200 m, 3 user/m, 4, m, , W and ) or (b) contention control (200 m, 3 user/m, 4, m, W and ).
Fig. 6. Impact of hidden primary receiver on interference distribution for CR networks under hybrid power/contention control scheme (200 m, 3 user/m, 4, , , W, and ).
Fig. 7. Impact of various CR deployment parameters on the aggregated interference for CR networks with (a) power control (100 m, 3 user/m, 4, m, and W) or (b) contention control (100 m, 3 user/m, 4, 20 m, and W).
Fig. 8. Impact of shadow fading on the aggregated interference for CR networks with (a) power control (100 m, 3 user/m, 4,