Abstract

Optimal resource allocation for cooperative cognitive radio networks with opportunistic access to the licensed spectrum is studied. Resource allocation is based on minimizing the symbol error rate at the receiver. Both the cases of all-participate relaying and selective relaying are considered. The objective function is derived and the constraints are detailed for both scenarios. It is then shown that the objective functions and the constraints are nonlinear and nonconvex functions of the parameters of interest, that is, source and relay powers, symbol time, and sensing time. Therefore, it is difficult to obtain closed-form solutions for the optimal resource allocation. The optimization problem is then solved using numerical techniques. Numerical results show that the all-participate system provides better performance than its selection counterpart, at the cost of greater resources.

1. Introduction

The ever increasing wireless communication networks have put great stress on the already limited spectrum. Due to the fixed spectrum allocation policy, only the licensed users, otherwise known as primary users, are able to access the licensed spectrum. Additionally, the Federal Communications Commission (FCC) task force highlighted in their report the fact that at any given time only 2% of the spectrum is being used [1]. Therefore, ensuring better spectrum usage is of paramount importance.

Cognitive radios have been proposed to resolve this issue [2]. In cognitive radio, the unlicensed users, otherwise known as secondary users, first sense the licensed bands for spectrum holes (parts of the licensed spectrum which are not being employed by the primary users at some time in certain geographical location) [3]. Then, if a spectrum hole is found, the secondary users transmit data to the intended destination. However, the secondary user has to be careful so as not to cause interference to the primary user. The two stages of spectrum sensing and data transmission are related and for optimal performance must be optimized jointly. This is due to the probability of detection, 𝑃𝑑, and probability of false alarm, 𝑃𝑓, associated with spectrum sensing. If the secondary user with probability 𝑃𝑓 misses a spectrum hole, then it will keep silent and miss an opportunity to transmit, reducing throughput. However, if a transmission from the primary is missed, with probability 1βˆ’π‘ƒπ‘‘, then the secondary user transmits and causes interference to the primary user. Moreover, due to interference, the signal-to-noise-ratio (SNR) of the secondary user also decreases, decreasing the throughput and the symbol error rate (SER). Resource allocation that optimizes this sensing-throughput tradeoff has been discussed in [4]. Other optimal resource allocation algorithms for cognitive radio networks have been discussed in [5]. More specifically, in [5], the authors considered a multiband system and considered the two cases of sensing-based spectrum sharing and opportunistic spectrum access. However, both the above-mentioned works maximize the throughput.

In this paper, optimal resource allocation is discussed to minimize the SER. In order to achieve minimum SER, cooperation is introduced into the system as it decreases the SER due to diversity [6]. Hence, the transmitting secondary user now, upon sensing a spectrum hole, transmits to the relays as well as the destination. Power allocation for relay-assisted cognitive radio networks has been discussed in [7–15]. These works proposed strategies to maximize the throughput for a cognitive relay network that is allowed to share the frequency band with the primary user. Thus, they did not consider spectrum sensing for opportunistic access. Here, we consider an opportunistic system with amplify-and-forward (AF) relays. Full channel state information is assumed at the central controller which performs the resource allocation. Firstly, an all-participate (AP) system is discussed and it is shown that the optimization problem is nonconvex and hence cannot be solved using analytical means. It is then noted that the AP system is limited due to the systems resources being orthogonally distributed. To rectify this, a selection scheme is proposed and the optimal resource allocation, in this case, is discussed.

The rest of the paper is organized as follows. Section 2 gives the system model. The AP system is considered in Section 3. Section 4 details selective relaying. Numerical results are discussed in Section 5. Finally, Section 6 concludes the paper.

2. System Model

Consider a cognitive radio network in which the secondary source utilizes π‘š relays to send data to the secondary destination as shown in Figure 1. The secondary network only has opportunistic access to the licensed spectrum. Therefore, it needs to perform spectrum sensing. The source performs the spectrum sensing and then transmits information to the relays and destination if it finds a β€œspectral hole” in the first time slot. The relays then forward the received signal to the destination after amplification. In this paper narrowband channel is assumed. The source and the relays can use frequency orthogonal channels to avoid causing interference to each other in wideband channels. For ease of analysis, we consider time orthogonal channels here. Hence, a total of π‘š+1 time slots are used.

2.1. Received Signal Model

Based on the spectrum sensing result, there are two possible received signal models.

2.1.1. Without Interference from the Primary User

In this scenario, with probability 1βˆ’π‘ƒπ‘“, where 𝑃𝑓 is the probability of false alarm, the source correctly detects the presence of a β€œspectral hole” and transmits. The signals received at the destination and the relays areπ‘¦π‘Ÿπ‘‘=√𝐸STβ„Žπ‘Ÿπ‘‘π‘ +π‘›π‘Ÿπ‘‘,𝑦𝑖𝑑=√𝐸STβ„Žπ‘–π‘‘π‘ +𝑛𝑖𝑑,(1) where 𝑠 is the zero-mean and unit-energy transmitted symbol, 𝐸ST is the source energy, β„Žπ‘Ÿπ‘‘ is the channel response between the source and the destination, β„Žπ‘–π‘‘ is the channel response between the 𝑖th relay and the source, and π‘›π‘Ÿπ‘‘βˆΌCN(0,𝜎2π‘›π‘Ÿπ‘‘) and π‘›π‘–π‘‘βˆΌCN(0,𝜎2𝑛𝑖𝑑) are the complex Gaussian noise samples. The relay, after normalization and amplification, forwards the received signal to the destination. The signal after normalization becomes 𝑠𝑖=ξ„Άξ„΅ξ„΅βŽ·πΈST𝐸ST||β„Žπ‘–π‘‘||2+𝜎2π‘›π‘–π‘‘β„Žπ‘–π‘‘ξ„Άξ„΅ξ„΅βŽ·π‘ +1𝐸ST||β„Žπ‘–π‘‘||2+𝜎2𝑛𝑖𝑑𝑛𝑖𝑑.(2)

Therefore, the received signal at the destination from the 𝑖th relay isπ‘¦π‘Ÿπ‘–=βˆšπΈπ‘–β„Žπ‘Ÿπ‘–π‘ π‘–+π‘›π‘Ÿπ‘–,(3) where β„Žπ‘Ÿπ‘– is the known channel response between the receiver and the 𝑖th relay, 𝐸𝑖 is the 𝑖th relay's energy, and π‘›π‘Ÿπ‘–βˆΌCN(0,𝜎2π‘›π‘Ÿπ‘–) is the complex Gaussian noise. Substituting 𝑠𝑖 in (3) gives π‘¦π‘Ÿπ‘–=ξ„Άξ„΅ξ„΅βŽ·πΈST𝐸𝑖𝐸ST||β„Žπ‘–π‘‘||2+𝜎2π‘›π‘–π‘‘β„Žπ‘Ÿπ‘–β„Žπ‘–π‘‘π‘ +Μƒπ‘›π‘Ÿπ‘–,(4) where Μƒπ‘›π‘Ÿπ‘–βˆΌCN(0,ξ‚πœŽ2π‘›π‘Ÿπ‘–) and ξ‚πœŽ2π‘›π‘Ÿπ‘–=𝐸𝑖||β„Žπ‘Ÿπ‘–||2𝜎2𝑛𝑖𝑑𝐸ST||β„Žπ‘–π‘‘||2+𝜎2𝑛𝑖𝑑+𝜎2π‘›π‘Ÿπ‘–.(5)

Writing the π‘š+1 received signals in matrix form, one has𝐲=𝐑𝑠+𝐧,(6) where 1𝐲=πœŽπ‘›π‘Ÿπ‘‘π‘¦π‘Ÿπ‘‘1ξ‚πœŽπ‘›π‘Ÿ1π‘¦π‘Ÿ1β‹―1ξ‚πœŽπ‘›π‘Ÿπ‘šπ‘¦π‘Ÿπ‘šξƒ­π‘‡,βŽ‘βŽ’βŽ’βŽ£ξƒŽπ‘=𝐸ST𝜎2π‘›π‘Ÿπ‘‘β„Žπ‘Ÿπ‘‘ξ„Άξ„΅ξ„΅βŽ·πΈST𝐸1ξ‚πœŽ2π‘›π‘Ÿ1𝐸ST||β„Ž1𝑑||2+𝜎2𝑛1π‘‘ξ‚Γ—β„Žπ‘Ÿ1β„Ž1π‘‘β‹―ξ„Άξ„΅ξ„΅βŽ·πΈπ‘†π‘‡πΈπ‘šξ‚πœŽ2π‘›π‘Ÿπ‘šξ‚€πΈST||β„Žπ‘šπ‘‘||2+𝜎2π‘›π‘šπ‘‘ξ‚Γ—β„Žπ‘Ÿπ‘šβ„Žπ‘šπ‘‘ξƒ­,𝑇(7) and 𝐧 is a (π‘š+1) dimensional vector whose components are zero mean and unit variance complex Gaussian random variables. Therefore 𝐧 is also complex Gaussian with mean vector of all zeros and covariance matrix being the identity matrix, 𝐈(π‘š+1)Γ—(π‘š+1), that is 𝐧∼CN(𝟎,𝐈).

2.1.2. With Interference from the Primary User

In this case, with probability 1βˆ’π‘ƒπ‘‘ where 𝑃𝑑 is the probability of detection, the source misses the transmission from the primary user and transmits, which causes interference. The signals at the destination from both the source and the relays now also include an interfering signal due to primary user activity isπ‘¦π‘Ÿπ‘‘=√𝐸STβ„Žπ‘Ÿπ‘‘π‘ +π‘›π‘Ÿπ‘‘+π‘¦πΌπ‘Ÿπ‘‘,𝑦𝑖𝑑=√𝐸STβ„Žπ‘–π‘‘π‘ +𝑛𝑖𝑑+𝑦𝐼𝑖𝑑,(8) where π‘¦πΌπ‘Ÿπ‘‘ and 𝑦𝐼𝑖𝑑 are the interference signals.

Taking into account the fact that the source and relays have no knowledge of the interfering signal and adopting the same approach as previously, one can writeπ‘¦π‘Ÿπ‘–=ξ„Άξ„΅ξ„΅βŽ·πΈST𝐸𝑖𝐸ST||β„Žπ‘–π‘‘||2+𝜎2π‘›π‘–π‘‘β„Žπ‘Ÿπ‘–β„Žπ‘–π‘‘π‘ +Μ‚π‘›π‘Ÿπ‘–+Μ‚π‘¦πΌπ‘Ÿπ‘–,(9) where Μ‚π‘›π‘Ÿπ‘–βˆΌCN(0,𝜎2π‘›π‘Ÿπ‘–) and 𝜎2π‘›π‘Ÿπ‘–=𝐸𝑖||β„Žπ‘Ÿπ‘–||2𝜎2𝑛𝑖𝑑𝐸ST||β„Žπ‘–π‘‘||2+𝜎2𝑛𝑖𝑑+𝜎2π‘›π‘Ÿπ‘–,̂𝑦𝐼𝑖=π‘¦πΌπ‘Ÿπ‘–+ξƒŽπΈπ‘–πΈST|β„Žπ‘–π‘‘|2+𝜎2π‘›π‘–π‘‘β„Žπ‘Ÿπ‘–π‘¦πΌπ‘–π‘‘.(10)

Again in matrix form one has𝐲𝐼=𝐑𝐼𝑠+𝐧𝐼+𝐘𝐼,(11) where 𝐲𝐼=1πœŽπ‘›π‘Ÿπ‘‘π‘¦π‘Ÿπ‘‘1ξπœŽπ‘›π‘Ÿ1π‘¦π‘Ÿ1β‹―1ξπœŽπ‘›π‘Ÿπ‘šπ‘¦π‘Ÿπ‘šξƒ­π‘‡,𝐑𝐼=βŽ‘βŽ’βŽ’βŽ£ξƒŽπΈST𝜎2π‘›π‘Ÿπ‘‘β„Žπ‘Ÿπ‘‘ξ„Άξ„΅ξ„΅βŽ·πΈST𝐸1𝜎2π‘›π‘Ÿ1𝐸ST||β„Ž1𝑑||2+𝜎2𝑛1π‘‘ξ‚Γ—β„Žπ‘Ÿ1β„Ž1π‘‘β‹―ξ„Άξ„΅ξ„΅βŽ·πΈSTπΈπ‘šξπœŽ2π‘›π‘Ÿπ‘šξ‚€πΈST||β„Žπ‘šπ‘‘||2+𝜎2π‘›π‘šπ‘‘ξ‚Γ—β„Žπ‘Ÿπ‘šβ„Žπ‘šπ‘‘ξƒ­π‘‡,𝐘𝐼=1πœŽπ‘›π‘Ÿπ‘‘π‘¦πΌπ‘Ÿπ‘‘1ξπœŽπ‘›π‘Ÿ1̂𝑦𝐼1β‹―1ξπœŽπ‘›π‘Ÿπ‘šΜ‚π‘¦πΌπ‘šξƒ­π‘‡,(12) and 𝐧𝐼∼CN(𝟎,𝐈).

2.2. Spectrum Sensing

Spectrum sensing is performed, by means of an energy detector, for the first 𝑑𝑠 seconds out of total time slot duration of 𝑇 seconds at the source node only. The remaining π‘‡βˆ’π‘‘π‘  is used for transmission, after detecting a β€œspectral hole”. The probabilities of detection and false alarm, according to [16], are given by π‘ƒπ‘‘βŽ›βŽœβŽœβŽœβŽ=π‘„πœ†βˆ’π‘βˆ’π›Ύπ‘‘ξ”2𝑁+2π›Ύπ‘‘ξ€ΈβŽžβŽŸβŽŸβŽŸβŽ ,𝑃𝑓=π‘„πœ†βˆ’π‘βˆšξƒͺ,2𝑁(13) respectively, where πœ† is the threshold of the energy detector, 𝑁=𝑑𝑠𝑓𝑠 is the number of samples, 𝑓𝑠 is the sampling frequency, 𝛾𝑑 equals 𝑁 times the SNR at the output of the detector and 𝑄(β‹―) is the Gaussian 𝑄-function.

3. All Participate System

In this section, an all-participate (AP) system is discussed. In such a system, all the relays forward the signal to the destination. Firstly, the optimization problem is formulated. Then the constraints on the objective function are derived. The SER at the destination is given by𝐻SER=𝑃0ξ€Έξ‚€π‘„ξ‚€βˆšπ‘˜π›Ύ01βˆ’π‘ƒπ‘“ξ€Έξ€·π»+𝑃1ξ€Έξ‚€π‘„ξ‚€βˆšπ‘˜π›ΎπΌξ€·ξ‚ξ‚1βˆ’π‘ƒπ‘‘ξ€Έ,(14) where 𝛾0 is the SNR after combining, 𝛾𝐼 is the signal-to-interference-plus-noise-ratio (SINR) after combining, π‘˜ is a constant which depends on the modulation scheme used, 𝑃(𝐻0) is the probability of no primary user transmission, and 𝑃(𝐻1)=1βˆ’π‘ƒ(𝐻0) is the probability of a primary user transmission. The SNR 𝛾0 can be found, assuming maximal ratio combining (MRC), as𝛾0=𝛼+π‘šξ“π‘–=1π›½π‘–ξ‚πœŽ2π‘›π‘Ÿπ‘–,(15) where 𝑝𝛼=ST𝑇𝑠||β„Žπ‘Ÿπ‘‘||2𝜎2π‘›π‘Ÿπ‘‘,𝛽𝑖=𝑝ST𝑝𝑖𝑇2𝑠||β„Žπ‘Ÿπ‘–||2||β„Žπ‘–π‘‘||2𝑝ST𝑇𝑠||β„Žπ‘–π‘‘||2+𝜎2𝑛𝑖𝑑,(16) where the source and relay energies have been replaced by 𝐸ST=𝑝ST𝑇𝑠,𝐸𝑖=𝑝𝑖𝑇𝑆,(17) where 𝑝ST and 𝑝𝑖s are the source and relay powers, respectively, and 𝑇𝑠 is the symbol time. Similarly, 𝛾𝐼 can be expressed as𝛾𝐼=ξ‚€βˆ‘π›Ό+π‘šπ‘–=1𝛽𝑖/𝜎2π‘›π‘Ÿπ‘–ξ‚ξ‚2βˆ‘π›Ό+π‘šπ‘–=1𝛽𝑖/𝜎2π‘›π‘Ÿπ‘–ξ‚+𝛼𝑐/𝜎2π‘›π‘Ÿπ‘‘ξ‚+βˆ‘π‘šπ‘–=1𝑑𝑖/𝜎2π‘›π‘Ÿπ‘–π›½ξ‚ξ‚€π‘–/𝜎2π‘›π‘Ÿπ‘–ξ‚,(18) where ξ‚Έ|||𝑦𝑐=πΈπΌπ‘Ÿπ‘‘|||2𝑑,(19)𝑖||=𝐸̂𝑦𝐼𝑖||2𝐸||,(20)̂𝑦𝐼𝑖||2ξ‚„ξ‚Έ|||𝑦=πΈπΌπ‘Ÿπ‘–|||2ξ‚Ήξ‚Έ|||𝑦+𝐸𝐼𝑖𝑑|||2𝐸𝑖𝐸𝑆𝑇|β„Žπ‘–π‘‘|2+𝜎2𝑛𝑖𝑑||β„Žπ‘Ÿπ‘–||2ξƒͺ.(21)

After substituting (15) and (18) in (14), the SER can be obtained as𝐻SER=𝑃0ξ€Έπ‘„βŽ›βŽœβŽœβŽξ„Άξ„΅ξ„΅ξ„΅βŽ·π‘˜βŽ›βŽœβŽœβŽπ›Ό+π‘šξ“π‘–=1π›½π‘–ξ‚πœŽ2π‘›π‘Ÿπ‘–βŽžβŽŸβŽŸβŽ βŽžβŽŸβŽŸβŽ Γ—ξƒ©ξƒ©1βˆ’π‘„πœ†βˆ’π‘βˆšξ€·π»2𝑁ξƒͺξƒͺ+𝑃1ξ€ΈβŽ›βŽœβŽœβŽœβŽβŽ›βŽœβŽœβŽœβŽ1βˆ’π‘„πœ†βˆ’π‘βˆ’π›Ύπ‘‘ξ”2𝑁+2π›Ύπ‘‘ξ€ΈβŽžβŽŸβŽŸβŽŸβŽ βŽžβŽŸβŽŸβŽŸβŽ π‘„Γ—βŽ›βŽœβŽœβŽœβŽξ„Άξ„΅ξ„΅ξ„΅ξ„΅βŽ·π‘˜ξ‚€βˆ‘π›Ό+π‘šπ‘–=1𝛽𝑖/𝜎2π‘›π‘Ÿπ‘–ξ‚ξ‚2βˆ‘π›Ό+π‘šπ‘–=1𝛽𝑖/𝜎2π‘›π‘Ÿπ‘–ξ‚+𝛼𝑐/𝜎2π‘›π‘Ÿπ‘‘ξ‚+βˆ‘π‘šπ‘–=1𝑑𝑖/𝜎2π‘›π‘Ÿπ‘–π›½ξ‚ξ‚€π‘–/𝜎2π‘›π‘Ÿπ‘–ξ‚βŽžβŽŸβŽŸβŽŸβŽ .(22)

Now we form the different constraints on the problem. First, we consider both individual power constraints at the source and the relay and a global power constraint on the whole system. Therefore, the constraints are given by0≀𝑝ST≀𝑝𝑇,0≀𝑝𝑖≀𝑝𝑖max,𝑃𝑓≀𝑃th𝑓,𝑝ST+π‘šξ“π‘–=1𝑝𝑖≀𝑝total,𝑅𝑇𝑠+𝑑𝑠≀𝑇,(23) where 𝑝ST is the power available at the source, 𝑝𝑖max is the power available at each relay, 𝑝total is the power available to the whole system, and 𝑃th𝑓 specifies the constraint on the probability of false alarm. The constraints on 𝑃𝑓, 𝑇𝑠, and 𝑑𝑠 are introduced to maintain an acceptable throughput. Next the two cases of global power constraint only and individual power constraints only are considered. For the case of global power constraint only, the constraints will be𝑃𝑓≀𝑃th𝑓,𝑝ST+π‘šξ“π‘–=1𝑝𝑖≀𝑝total,𝑅𝑇𝑠+𝑑𝑠≀𝑇.(24)

In the other scenario, the constraints are given by0≀𝑝ST≀𝑝𝑇,0≀𝑝𝑖≀𝑝𝑖max,𝑃𝑓≀𝑃th𝑓,𝑅𝑇𝑠+𝑑𝑠≀𝑇.(25)

The individual power constraints are set to limit the interference suffered by the primary user in the case of missed detection. As there is no individual power constraint, the interference caused to the user in the global power constraint only case, where the primary user is only protected by spectrum sensing, is greater.

The problem with optimizing (22) is that it is a nonlinear and nonconvex function due to the Gaussian Q-function being nonlinear and, in general, nonconvex. Thus, the Lagrangian multiplier method [17] cannot be applied here to obtain closed form expressions of the optimal resource allocation. One has to resort to numerical techniques to obtain the optimal solution.

A special case of importance is the absence of the direct link between source and destination, because the relays take on a more prominent role in the presence of no direct link. In this case, the SER is can be obtained by setting 𝛼=0 in (22).

4. Selective Relaying

The drawback of the all-participate (AP) scheme discussed in the previous section is that to avoid causing interference, the source and the relay transmit on orthogonal channels. Hence, consuming a considerable amount of resources. In our discussion of a time orthogonal systems, π‘š+1 time slots are utilized for the transmission of one data frame. Additionally, as no sensing is performed at the relays, the primary may become active over any one of the π‘š time slots and cause interference.

To overcome these problems, a selection scheme is proposed in this section in which only one relay is selected to take part in forwarding the signal from the source. Now only 2 time slots are used in transmitting one frame of data and thus decreasing the likelihood of primary becoming active again during relay transmission. In the selection case, the SER is𝐻SER=𝑃0ξ€Έπ‘„βŽ›βŽœβŽœβŽξ„Άξ„΅ξ„΅ξ„΅βŽ·π‘˜βŽ›βŽœβŽœβŽπ›½π›Ό+π‘—ξ‚πœŽ2π‘›π‘Ÿπ‘—βŽžβŽŸβŽŸβŽ βŽžβŽŸβŽŸβŽ ξƒ©ξƒ©1βˆ’π‘„πœ†βˆ’π‘βˆšξ€·π»2𝑁ξƒͺξƒͺ+𝑃1ξ€Έπ‘„Γ—βŽ›βŽœβŽœβŽœβŽξ„Άξ„΅ξ„΅ξ„΅ξ„΅βŽ·π‘˜ξ‚€ξ‚€π›½π›Ό+𝑗/𝜎2π‘›π‘Ÿπ‘—ξ‚ξ‚2𝛽𝛼+𝑗/𝜎2π‘›π‘Ÿπ‘—ξ‚+𝛼𝑐/𝜎2π‘›π‘Ÿπ‘‘ξ‚+𝑑𝑗/𝜎2π‘›π‘Ÿπ‘—π›½ξ‚ξ‚€π‘—/𝜎2π‘›π‘Ÿπ‘—ξ‚βŽžβŽŸβŽŸβŽŸβŽ Γ—ξƒ©ξƒ©1βˆ’π‘„πœ†βˆ’π‘βˆ’π›Ύβˆš,2(𝑁+2𝛾)ξƒͺξƒͺ(26) where 𝑏𝑗, 𝑑𝑗, and 𝜎2π‘›π‘Ÿπ‘— correspond to the 𝑗th relay that is selected. The SER in (26) is first optimized for all the relays and the relay which gives the minimum optimum SER is selected. Again, all three cases given in (23), (24) and (25) of both global and individual constraints, global constraint only and individual constraint only are considered. The selection criteria of minimizing SER adds complexity. However, such a criteria provides results which can serve as a benchmark as minimizing SER is the optimal selection criteria.

It is again evident that, even in the selection case, the SER is still a nonlinear and nonconvex function. Therefore, one has to resort to numerical techniques to find the optimal solution. The special case of no direct link is again of particular interest and considered separately.

5. Numerical Results

In this section, numerical results are provided for the optimization problems discussed. First, the proposed AP system with optimal resource allocation is discussed and it is shown that the proposed AP schemes give better performance than the uniform power allocation (UPA) scheme. In UPA, the power is uniformly distributed among the source and the relays and the sensing time and the symbol time are set so that the inequality 𝑅𝑇𝑠+𝑑𝑠≀𝑇 is satisfied. The selection scheme is discussed next and its performance is compared with selection with UPA. To make it easy for the reader to follow the discussion, a glossary is included in Table 1.

An interior-point algorithm was used to perform the optimization. The MATLAB function fmincon, which performs constrained optimization, is used to run the interior-point algorithm. To ensure that the algorithm converged to the optimal solution, the algorithm was run for a large number of initial values. All the variances are set as equal, 𝜎2π‘›π‘Ÿπ‘‘=𝜎2𝑛𝑖𝑑=𝜎2π‘›π‘Ÿπ‘–=𝜎2. The constraint of 𝑃𝑓 is set at 𝑃th𝑓=0.1. The total time duration is taken to be 100 ms. Binary phase shift keying (BPSK) is the modulation scheme employed. Due to the the number of samples and the sensing time being linearly related, the results are plotted against the number of samples.

The relationship between the number of samples (𝑁) and SER is shown in Figure 2. As one can clearly observe from Figure 2, there is an optimal value of the number of samples, hence for the sensing time, which minimizes the SER. This is because of the tradeoff between symbol time, 𝑇𝑠 and sensing time, 𝑑𝑠. Increasing sensing time means higher probability of detection which leads to a lower SER. However, an increase in sensing time comes at the cost of a decrease in symbol time which leads to lower 𝛾0 and 𝛾𝐼. Therefore, the SER increases. Similarly, decreasing sensing time implies a lower probability of detection and, in turn, higher SER. However, it also means high values of 𝛾0 and 𝛾𝐼 due to increase in symbol time, which in turn leads to lower SER. Hence, there exists an optimal value. This optimal value is affected by the primary user's SNR. The higher the primary user's SNR, the lower the optimal value of the sensing time will be as it takes shorter time to reach the same value of 𝑃𝑑 as it would take for a lower SNR of the primary user.

Figure 3 shows the relationship between the symbol time, 𝑇𝑠, and the SER. The relationship follows the opposite pattern as the sensing time. This is due to the constraint relating sensing time and the symbol time. Therefore, when the optimal value of the sensing time is low, the optimal value of the symbol time is high.

Figure 4 shows the SER performance of the different AP schemes plotted against 𝛾𝑠, where 𝛾𝑠=1/𝜎2. As expected, for the case with a direct link, the three optimal resource allocation (ORA) scenarios, global constraint only (GL), individual constraints only (Ind), and both global and individual constraints, outperform the uniform power allocation (UPA) and the direct link only for all values of 𝛾𝑠 and the gap in performance becomes wider with increasing 𝛾𝑠.

In Figure 4, in the case where there is no direct link (NDL) between source and destination, the performance is a little different. In this case, all three ORA schemes, AP-ORA-NDL, AP-ORA-GL-NDL, and AP-ORA-Ind-NDL, outperform the uniform power allocation scheme (UPA-NDL). However, for 𝛾𝑠 less than 0 dB, the direct link only case provides better SER performance than the three ORA cases with no direct link. The three ORA schemes with no direct link start to catch up to the direct link only scenario after 0 dB and completely outperform it after 5 dB. This phenomena coupled with the fact that AP-ORA-NDL, AP-ORA-GL-NDL, and AP-ORA-Ind-NDL, are handily outperformed by AP-ORA, AP-ORA-GL, and AP-ORA-Ind, respectively, demonstrate the significance of the presence of a link between the source and destination.

Comparing the different constraints, AP-ORA gives the worst performance in both scenarios of direct and no direct link. This is due to the fact that AP-ORA is constrained both globally and individually. Thus, even if one relay has more favourable conditions, the power allocated to it cannot exceed 𝑝𝑖max which is not the case in global constraint only case. In global constraint only case more power can be allocated to the source and relay which has more favourable conditions. The comparison between global constraint only and individual constraints only requires further elaboration.

First, the global constraints only and individual only scenarios are compared in the no direct link case. Here, AP-ORA-Ind-NDL provides lower SER than AP-ORA-GL-Ind for all values of 𝛾𝑠. AP-ORA-GL-NDL has the advantage of allocating more power to relays with better channel conditions. However, AP-ORA-Ind-NDL makes up for this advantage by having more total power in the system as it is not constrained by a total power constraint.

Now consider the direct link case. Here, AP-ORA-GL outperforms AP-ORA-Ind at low values of 𝛾𝑠 due to the presence of the direct link. As discussed before, the direct link is quite important, hence, as AP-ORA-GL is not limited by individual constraints, more power can be allocated to the direct link. This is not the case for AP-ORA-Ind. Therefore, AP-ORA-GL gives lower SER. However, with an increase in 𝛾𝑠, the noise decreases and the greater total power in the case of AP-ORA-Ind comes more into play. Thus, AP-ORA-Ind starts to outperform AP-ORA-GL at higher values of 𝛾𝑠. However, we must keep in mind that in the global power constraint only case, the interference to the primary is greater than those in the other two cases. Hence, the advantage in performance at low 𝛾𝑠 comes at the cost of greater interference to the primary.

Figure 5 shows the SER performance of the AP system as a function of the number of relays, π‘š. A similar pattern to Figure 4 is observed. The ORA schemes outperform the UPA schemes in both cases of direct and no direct link. Among the proposed ORA schemes, AP-ORA-Ind-NDL provides lower SER than AP-ORA-NDL and AP-ORA-GL-NDL in the no direct link scenario while in the direct link scenario AP-ORA is outperformed by both AP-ORA-Ind and AP-ORA-GL. In addition, AP-ORA-GL has better performance than AP-ORA-Ind for a small number of relays. However, as the number of relays increases AP-ORA-Ind surpasses AP-ORA-GL in terms of performance due to greater total power. Moreover, AP-ORA-NDL, AP-ORA-GL-NDL, and AP-ORA-Ind-NDL even start to outperform UPA for a large number of relays which shows the gain in performance due to ORA.

Figure 6 shows the performance of the various selection schemes as a function of 𝛾𝑠. The comparison in performance follows a similar pattern as in the AP case with the proposed selection ORA schemes outperforming their UPA counterparts and direct link only scenario. However, there is one major difference. In the presence of a direct link, Sel-ORA-Ind gives poorer performance than Sel-ORA-GL even for high values of 𝛾𝑠. Only at around 15 dB does Sel-ORA-Ind starts to catch up to Sel-ORA-GL. This is due to the fact that as pointed out in the AP system, in the case of global constraints only more power can be allocated to the source. However, unlike AP, in Sel there is only one additional relay which implies less total power for Sel-ORA-Ind and, therefore, requires a high value of 𝛾𝑠 to make the difference in total power count.

SER performance for selective relaying as a function of the number of relays is shown in Figure 7. Again, the main difference from the AP case is that Sel-ORA-GL outperforms Sel-ORA-Ind even for a large number of relays. This is due to the fact that even though the number of relays increase, the total power for Sel-ORA-Ind remains constant as only one relay in addition to the source takes part in data transmission. An interesting point to note here is that the there seems to be a minimum threshold for the SER for the selective system.

Figures 8 and 9 show the performance comparison between the AP and Sel system as a function of 𝛾𝑠 and π‘š, respectively. The comparison is presented separately for clarity, as if it was included in the previous figures, they would have become cluttered. From Figure 8, one can see that AP scheme outperforms the selection scheme in all scenarios, however, the gap in performance is not too big. This is due to the fact that the total number of relays is 3. If π‘š is increased, the performance gap will also increase. Still, one has to keep in mind the extra cost and spectral inefficiency associated with the AP scheme. This becomes more clear when the Figure 9 is examined.

As one can see, the difference in performance between the respective AP schemes and Sel schemes increases with increase in number of relays. As discussed earlier, the Sel schemes look to be bounded by a minimum threshold. Due to this, Sel with direct link scenarios even fall below the AP with no direct link scenarios for a large number of relays.

6. Conclusions

In this paper, ORA for a cognitive relay network has been discussed. It has been shown that for an AP system that ORA improves SER performance and the discussed schemes outperform the UPA schemes. The importance of the direct link between the source and the destination has also been demonstrated. Among the different constraints on the system, the case of both individual and global constraints gives the worst performance while global constraints only is the best for low 𝛾𝑠. However, this comes at the cost of greater interference to the primary user. The individual constraints only case takes over as the best scheme as 𝛾𝑠 increases.

It was then noted that the AP scheme consumes considerable resources and is spectrally inefficient. Therefore, a simple relay selection scheme has been proposed. Optimal resource allocation was then discussed for the selection scheme. The performance comparison of the AP and Sel shows that while AP provides better SER performance, it comes at the cost of considerable resources.

Acknowledgment

This work was supported by King Abdullah University of Science and technology (KAUST).