## I Introduction

Cognitive radio (CR) is an important technology to maximize radio spectrum utilization efficiency[1]-[3]. In CR systems, the secondary network is allowed to share the spectrum allocated to the primary network provided that the interference caused by the secondary transmitter (ST) does not deteriorate the performance of the primary network. Consequently, the challenge is to maintain the interference at the primary receiver (PR) below a pre-determined threshold level. This can be achieved by adapting the ST transmit power that ensures satisfaction of the interference constraint at the PR [4].

Multiuser diversity is considered an important diversity technique to improve wireless communication systems performance [5]. Considering a multiuser network where the users experience independent fading conditions, the basic idea of multiuser diversity is to select the users with the best fading conditions for transmission or reception to obtain a specific performance gain. Multiuser diversity in CR systems has attracted much attention recently. Researchers analyze the performance of multiuser diversity techniques for uplink multiuser underlay CR systems without taking the interference from the primary network into consideration[6] -[11]. In particular, the ergodic capacity (throughput) of multiuser diversity gain of uplink multiuser underlay CR systems is investigated in [6]. In [7], the authors analyze the achievable capacity gain of uplink multiuser spectrum-sharing systems over dynamic fading environments. In [8], the outage probability and effective capacity are analyzed for opportunistic spectrum sharing in Rayleigh fading environment. In [9], the authors analyze the outage probability, average symbol error rate (SER) and ergodic capacity of an opportunistic multiuser cognitive network with multiple primary users assuming the channels in the secondary network are independent but not identical Nakagami- fading. In [10], [11] the authors analyze the outage probability and average capacity of multiuser diversity in single-input multiple-output (SIMO) spectrum sharing systems. The mentioned previous works did not consider the impact of interference from the primary network to the secondary network. However, in the practical CR systems the interference from the primary network will greatly affect the performance of the secondary network. The ergodic capacity of multiuser diversity in CR systems with interference from the primary network is investigated in [12]. Recently, the ergodic capacity of various multiuser scheduling schemes in downlink cognitive radio networks with interference from the primary network is analyzed in [13]; here, the authors analyze the ergodic capacity of multiuser diversity scheduling under the outage constraint of multiple primary user receivers and the secondary user (SU) maximum transmit power limit.

When the interference from the primary transmitter (PT) is much larger than the noise at the secondary receiver, the performance of the secondary network is limited by the interference from the primary transmitter and the quantity of interest is signal-to-interference ratio (SIR). Such CR system can be described as interference-limited underlay CR system [14]. For the interference-limited underlay CR systems considered in [14], the authors analyze the average bit error rate (BER) and outage probability for receive antenna selection schemes under discrete power adaptation at the ST.

### I-a Motivation

The mentioned previous works have only focused on conventional multiuser diversity in underlay CR systems where the SU with the best link quality is selected. Furthermore, no prior work has considered the performance of the conventional multiuser diversity for interference-limited underlay CR systems with continuous power adaptation at the ST. Accordingly, we focus in this paper on a generalized multiuser diversity scheme that features selection of the

-th best SU for interference-limited secondary multiuser network under continuous power adaptation. The -th best SU selection is of practical interest in underlay CR systems since the best SU may not be selected under given traffic conditions. This might happen when the best user is unavailable or occupied by other service requirements [15], in handoff situations [16] or due to scheduling delay [17]. Clearly, the -th best SU selection scheme includes the best SU selection at as a special case.In general, it is difficult to analyze the exact performance of the -th best SU selection scheme for arbitrary number of secondary users. Hence, we use the extreme value theorem (EVT) [18] to analyze the asymptotic performance (in the limit of large number of secondary users) of such selection scheme. As we will show later, EVT provides tractable and accurate asymptotic expressions for the average throughput, effective throughput, average bit error rate and outage probability. The derived asymptotic expressions are accurate for practical CR systems with not so large (realistic) number of secondary users. We validate the accuracy of the asymptotic expressions through numerical simulation.

It should be noted that the derived expressions are obtained by assuming perfect channel state information (CSI) of the secondary transmitter to primary receiver channel. The impact of imperfect CSI on the mathematical analysis is also investigated. As we will discuss later, the derived mathematical expressions under the perfect CSI can be used to deduce the system performance under the imperfect CSI case.

### I-B Related work and new contributions

In our previous work [19], we used EVT to derive simple closed-form asymptotic expressions for the average throughput, effective throughput and average BER for the link with the -th highest signal-to-noise ratio (SNR) in traditional wireless communication systems with no spectrum sharing. We showed that the -th highest SNR converges uniformly in distribution to a Log-Gamma random variable ^{1}^{1}1 If is a Gamma random variable, then we say that is a Log-Gamma random variable whose support is the real line [20].. As a special case, if , the Log-Gamma reduces to the Gumbel random variable. The average throughput, effective throughput and average BER were derived for various channel models that are widely used to characterize fading in wireless communication systems such as Weibull, Gamma, and Gamma-Gamma.

Unlike[19], in this paper we consider a downlink interference-limited secondary multiuser network, where the noise at each secondary user receiver is negligible compared to the interference from the PT. We assume that the ST transmits information to the -th best secondary user (SU), namely, the SU with the -th highest SIR. Meanwhile, the ST adopts continuous limited power adaptation strategy to satisfy the instantaneous interference constraint at the PR. Our contribution is to utilize EVT to analyze the performance of the -th best SU for underlay CR systems. More specifically, we show that the -th highest SIR converges uniformly in distribution to an inverse gamma random variable for a fixed and large number of secondary users^{2}^{2}2Here, the -th highest SIR converges to an inverse gamma random variable. This is different from [19], where the -th highest SNR converges to a Log-Gamma random variable.. Then, we derive novel closed-form asymptotic expressions for the average and effective throughputs of the -th best SU employing continuous power adaptation at the ST with both limited and unlimited transmit power ^{3}^{3}3 A portion of this work was accepted for publication in 2018 IEEE Vehicular Technology Conference (VTC) [21]. In [21], we only studied the asymptotic average and effective throughputs of the -th best secondary user selection in uplink multiuser cognitive radio systems under unlimited transmit power. The VTC paper is available online at https://arxiv.org/pdf/1804.05257v2.pdf. Different from [21], in this paper we analyze the average and effective throughputs, average BER and outage probability of the -th best secondary user selection for interference limited CR systems under both limited and unlimited ST transmit power adaptation strategies. We also investigate the impact of imperfect CSI of the secondary transmitter to primary receiver channel. . Furthermore, novel closed-form asymptotic expressions for the average BER and outage probability with continuous power adaptation and unlimited ST power are derived. It should be noted that although we use EVT approach in this paper, the analysis and the derived expressions are completely different from what were obtained for traditional wireless communication systems with no spectrum sharing.

The rest of this paper is organized as follows. In Section II we discuss the system model. In Section III we analyze the asymptotic average throughput, effective throughput, average BER and outage probability of the -th best SU. Section IV includes numerical results and Section V concludes.

## Ii System Model

As shown in Fig. 1, we consider an underlay secondary network consisting of one ST equipped with a single antenna, and secondary users each equipped with a single antenna. The secondary network is sharing the spectrum of a primary network with one PT and one PR. The PT and PR are equipped with a single antenna each. Let denote the channel gain from the PT to the -th secondary user’s receiver (SU-Rx), where . Let and denote the channel gain from the ST to the PR and the -th SU-Rx, respectively. We assume that the primary network is far away from the secondary network and therefore and are assumed to be independent Rayleigh distributed random variables. This implies that the channel power gains, and

have probability density functions (PDFs)

and , respectively, where is the unit step function and the parameters and are the fading parameters. The channel power gains in the secondary network, , for , are assumed to be independent and identically distributed (i.i.d.) Gamma random variables with PDF^{4}

^{4}4 We assume that is a Nakagami-m random variable; hence, is a Gamma random variable. This a is a generalized fading model that includes many practical scenarios. First, Nakagami-m fading includes the Rayleigh fading as a special case when , then the distribution of becomes consistent with the distributions of and . Second, in the situation where line of-sight path exists between the secondary transmitter and secondary receivers the natural choice to use the Rician distribution to model the line of-sight effect. It is well known that the Rician distribution can be accurately approximated by the Nakagami-m distribution. Motivated by these reasons, we adopt the Nakagami-m fading to model the fading in the secondary network

(1) |

where the parameters and are positive reals and is the Gamma function.

Similar to [4], [7], [9], [14], [22] and [23], it is assumed that the ST has perfect CSI regarding the secondary transmitter to primary receiver channel, . The ST can be informed about through a mediate band manager between PR and ST [24] or by considering proper signaling [25]. However, the impact of imperfect CSI on the performance of the -th best SU will be discussed later in this paper.

With a perfect knowledge of , we consider a continuous power adaptation policy at the ST to control its interference to the PR such that the instantaneous transmit power of the ST is

(2) |

where is the maximum instantaneous power available at the ST and is the maximum tolerable interference level at the PR. Assuming the noise at the -th SU-Rx is negligible compared to the interference from the PT, then the ST will select the -th best SU; namely, the SU with the -th highest SIR; i.e.,

(3) |

where , is the transmit power of the PT and is the PT interference power at the -th SU-Rx.

Let denote the instantaneous SIR at the -th best SU-Rx, where . According to [18], the PDF of can be expressed in terms of the PDF,

, and cumulative distribution function (CDF),

, of as(4) |

where the CDF and PDF of are given by [14]

(5) |

(6) |

respectively. Let denote the instantaneous throughput of the -th best SU, where and is the system bandwidth. Then, the average throughput of the -th best SU, , can be evaluated as

(7) |

in bit/s/Hz. The expectation in (7

) is taken over the joint distribution of random variables

and . Assuming a block fading channel, the effective throughput that can be supported by a wireless system under a statistical QoS constraint described by the delay QoS exponent is given by [26](8) |

where is a random variable which represents the instantaneous throughput during a single block and is the block length. implies there is no delay constraint and the effective throughput is then the ergodic (average) throughput of the corresponding wireless channel. Hence, the effective throughput of the -th best SU, , can be expressed as

(9) |

in bit/s/Hz, where and the expectation is taken over the joint distribution of and .

If we conisder a general class of modulation schemes whose conditional BER, , is given by

(10) |

where and are positive constants and is a random variable which represents the instantaneous received SIR, the average BER of the -th best SU can be expressed as

(11) |

where the expectation is taken over the joint distribution of and .

Due to the complicated nature of the distribution of the instantaneous SIR at the -th best SU-Rx, it is difficult to obtain exact expressions for , and . Therefore, in this paper we consider another approach based on EVT to analyze the performance of the -th best SU in terms of average throughput, effective throughput, outage probability and average BER.

## Iii Asymptotic Performance Analysis

In this section, we derive the limiting distribution of in Proposition 1 below. Then we use this result to analyze the average, effective throughputs, average BER and outage probability of the -th best SU.

### Iii-a The Limiting Distribution of

Proposition 1: Let denote the -th largest order statistic of i.i.d. random variables with a common CDF of , as expressed in (5), then for a fixed and , converges in distribution to a random variable with CDF

, which can be characterized by an inverse gamma distribution as

(12) |

where , and is the upper incomplete gamma function [27]. Furthermore, the PDF of , , can be obtained as

(13) |

Proof: We first investigate the limiting distribution of , which denotes the first largest order statistic of i.i.d. random variables. From Proposition 2 of [14], converges in distribution to a unit Fréchet distribution i.e.,

(14) |

where and .

Applying Proposition 1 of [19] with as in (14), it follows that for a fixed and , the sequence converges in distribution to a random variable with CDF of , which can be expressed in terms of as

(15) |

Using the fact that for an integer , can be finally expressed as in (12). By differentiating (12) we obtain (13).

It should be noted here that Proposition 1 of [19] can be applied for different CDF functions. In this paper we focus on the case when represents a Fréchet CDF. In this case has a limiting distribution of inverse gamma as shown in (12). This is different from what was obtained in [19], where Proposition 1 of [19] was applied for the case when represents Gumbel CDF and thus has a limiting distribution of Log-Gamma.

### Iii-B The Distribution of the ST Transmit Power

We consider a continuous power adaptation scheme in which the transmit power of the ST can be adapted with a power limit of ; therefore, the instantaneous transmit power of the ST is . Furthermore, we consider a continuous power adaptation scheme in which the transmit power of the ST can be adapted without any power limit, i.e., [22], [23]. In such case, the ST transmit power, , can be written as . We focus next on the PDF of the instantaneous transmit power of the ST, , then we use this PDF and Proposition 2 to evaluate the average and effective throughputs of the -th best SU.

Considering

is a continuous random variable, where

and is constant, then the CDF of the random variable , , can be given as(16) |

where is the CDF of the random variable and is the unit step function. Then it follows that the PDF of , , can be expressed as

(17) |

where is the PDF of the random variable and is the Dirac delta function, the derivative of .

Using the PDF of , and variable transformation then it follows that and . Finally we can write

(18) |

### Iii-C Average and Effective Throughputs

Proposition 2: The average and effective throughputs of the -th best SU for continuous power adaptation with limited ST power are respectively given by

(19) |

(20) |

for fixed and , where is the exponential integral function [28, Eq. (5.1.4)], is the upper incomplete gamma function [27, Eq. (8.350.2)] is the digamma function [27, Eq. (8.360.1)].

Proof: Average Throughput:

From Proposition 1, the CDF of approaches the CDF of for a fixed and , where the CDF of is as expressed in (12). Or equivalently, the PDF of can be approximated by the PDF of for a fixed and , where the PDF of is as expressed in (13). Then for a fixed , and conditioning on the ST transmit power , can be approximated as

(21) |

Noting that is an increasing function of , we have in (21) for large . Using this and variable transformation of , can be further approximated as

(22) |

where the above integral is evaluated with help of [27, Eq. (4.352. 1) ]. Averaging over the PDF of in (18) yields

(23) |

Using variable transformation of with help of [27, Eq. (4.331. 2) ] and after some basic algebraic manipulation, we have

(24) |

Combining (22), (23) and (24), it follows that is as expressed in (19).

Effective Throughput: Conditioning on the ST transmit power in (9) and by exploiting Lemma 2 of the Appendix, we infer that can be approximated as

(25) |

for fixed and . Making use as above of for large in (25) and variable transformation of , can be further approximated as

(26) |

where the above integral is evaluated with help of [27, Eq. (3.381.4)]. Averaging (26) over the PDF of in (18) yields

(27) |

Using variable transformation of and using the definition of the upper incomplete gamma function, , we have

(28) |

Combining (27), (28) and (9), it follows that is as expressed in (20).

While we focused in Proposition 2 on analyzing the average and effective throughputs under the limited ST power adaptation, i.e., , it should be noted that simpler expressions can be obtained under the unlimited ST power case. i.e., . These expressions are useful when [23], [22] and they serve as upper bounds on the average and effective throughput under the limited ST power case. Using the result from Proposition 2, we derive the average and effective throughputs of the -th best SU with unlimited ST power in the following corollary.

Corollary 1: The average and effective throughputs of the -th best SU for continuous power adaptation with unlimited ST power are respectively given by

(29) |

(30) |

for fixed and , where is the Euler constant.

Proof:

### Iii-D Average BER

We now derive the average BER for the limited and unlimited continuous ST power in the following proposition.

Proposition 3: The average BER of the -th best SU for continuous limited power adaptation scheme can be approximated as

(32) |

for fixed and , where is the modified Bessel function of the second kind and order [28, Eq. (8.407.1)] Furthermore, for the unlimited ST transmit power, the average BER of the -th best SU can be approximated as

(33) |

for fixed and , where is the Meijer G-function [29].

Proof: Conditioning on the ST transmit power in (11) and by exploiting Lemma 1 of the Appendix, we infer that can be approximated as

(34) |

for a fixed and , where the above integral is evaluated with help of [30, Eq. (2.11)]. It is hard to find an analytical expression for the average BER for the limited ST transmit power case. Therefore, averaging (34) over in (18) yields the average BER for limited ST power as in (32). For the unlimited ST transmit power, one can show that by letting in (32), we have

(35) |

The above integral can be expressed in terms of the Meijer G-function as in (33). It should be noted that the Meijer G-function can be easily and efficiently computed using the most standard software packages like MAPLE and MATHEMATICA.

### Iii-E Outage Probability

We now derive the outage probability for limited and unlimited continuous ST power in the following proposition.

Proposition 4: The outage probability of the -th best SU for continuous limited power adaptation scheme can be approximated as

(36) |

for fixed and . Furthermore, for the unlimited ST transmit power, the outage probability of the -th best SU can be approximated as

(37) |

for fixed and .

Proof: The outage probability of the -th best SU, , can be expressed as

(38) |

where as given in (18). From Proposition 1, the CDF of approaches the CDF of for fixed and , where the CDF of is as expressed in (12). Then, we have

(39) |

Making use of (18) in (39), the outage probability for the limited ST power is as in (36). For the unlimited ST transmit power, one can show that by letting in (36), we have

(40) |

where the integral above is evaluated with help of [27, Eq. (6.453)] after variable transformation of .

### Iii-F Effect of Imperfect CSI

In practical environments, the ST has only a partial channel knowledge of the ST to PR channel, . In this case, the CSI on provided to the ST is outdated due to the time-varying nature of the wireless link [31]. The outdated CSI can be described using the correlation model as [31].

(41) |

where is the outdated channel information available at the ST,

is a complex Gaussian random variable with zero mean and unit variance, and uncorrelated with

. The correlation coefficient () is a constant, which is used to evaluate the impact of channel estimation error and feedback delay on the CSI

[31]. It is assumed that the ST knows the outdated channel information and the correlation coefficient as well. In view ofbeing an exponentially distributed random variable with parameter

, the estimated channel power gain is also an exponentially distributed RV with parameter , where .As we discussed in Section II, when the ST has a perfect CSI of , it can access the spectrum if the peak interference power constraint can be satisfied. However, it is hard to satisfy the instantaneous interference constraint at the PR if only the outdated CSI is available at ST [31]. Therefore, a more flexible constraint based on a pre-selected interference outage probability is adopted [31], [32]. Considering the imperfect CSI effect, the transmit power of the ST in (2) can be rewritten as [31]

(42) |

where denotes the power margin factor which can be expressed as [31]

(43) |

where denotes the predetermined interference outage probability. As a special case, a power margin factor of , i.e., () indicates the perfect CSI of and therefore the ST transmit power in (42) reduces to (2).

For further practical considerations, we address the imperfect CSI of the channels in the secondary network, , for . The outdated CSI can be described as

(44) |

where is the outdated channel information of the -th secondary link available at the ST, is a complex Gaussian random variable with zero mean and unit variance, and uncorrelated with . The correlation coefficient () is a constant that describes the impact of outdated CSI. In view of being a Gamma distributed random variable with parameters and , the estimated channel power gain is also a Gamma distributed random variable with parameters and .

It should be noted that the expressions derived for average throughput, effective throughput, average BER and outage probability in previous subsections assuming a perfect CSI hold for the imperfect CSI case after replacing with and with , due to the imperfect CSI on . And replacing with due to the imperfect CSI on , where .

## Iv Numerial results

### Iv-a Perfect CSI

In this subsection, we numerically illustrate and verify the obtained asymptotic expressions in Section III under the perfect CSI condition, which refers to with ST transmit power is as in (2) as described at the end of Section III. F .

In Fig. 2, we plot the average throughput of the -th best SU versus the number of secondary users, , for unlimited ST power, , and limited ST power with , for . We validate the obtained asymptotic expressions for the average throughput using Monte Carlo simulations. We observe that the accuracy increases as increases. Furthermore, we observe that the asymptotic results are accurate for not so large (realistic) values of . For example, is considered sufficiently large to confirm the accuracy of the asymptotic results compared to simulations. This suggests that the EVT is a powerful approach that approximates the performance for realistic and large values of as well. In Fig. 3, we plot the average throughput of the best SU versus the ST power, , in dB. We observe that, compared to the simulations, the accuracy of the asymptotic average throughput increases as increases from 6 to 30. We also observe that the accuracy of the asymptotic average throughput increases as increases. Furthermore, for larger values of the asymptotic average throughput approaches the one with unlimited ST power, .

In Fig. 4, we plot the average throughput of the best SU versus the interference level, , in dB for unlimited ST power, , and limited ST power with . Some interesting observations can be made from this figure. First, we observe that, compared to the simulations, the accuracy of the asymptotic average throughputs increase as increases from 20 to 200. Second, as or increases, the accuracy of the asymptotic average throughputs also increases. Last, for the limited ST power, , the average throughput is saturated and it does not improve as . This is due to the fact that for higher values of , the ST will select with a higher probability.

In Fig. 5, we plot the effective throughput of the -th best SU versus the number of secondary users, , with unlimited ST power, and limited ST power with , for . We observe that the accuracy of the asymptotic effective throughput increases as increases. However, it is shown that for , the asymptotic effective throughputs is less accurate for small to moderate values of

. This is because the asymptotic analysis is more accurate for large

relative to a fixed . Consequently, if the value of is close enough to , it is expected that the asymptotic expression will be less accurate.In Fig. 6, we plot the effective throughput of the best SU versus the delay exponent, , for and different values of . We observe that the effective throughput significantly decreases for smaller values of . On the other hand, for reasonably large values of , the effective throughput does not significantly improve compared to the unlimited ST power case, . This is due to the fact that for higher values of the effective throughput is dominated by the interference level .

In Fig. 7, the outage probability of the -th best SU is plotted versus interference level, , in at with unlimited ST power, and limited ST power with , for at dB. The saturation in the outage probability for the limited ST power case is due to the fact that for higher values of , the ST will select for most of the time. In Fig. 8, we plot the asymptotic average BER as a function of the number of secondary users, , for the unlimited ST power with

Comments

There are no comments yet.