# Energy-Efficient Power Optimization with Spectrum Sensing Errors in OFDMA Cognitive Opportunistic Relay Links

Authors: Li Feng; Yujun Kuang; Xinchuan Fu; Binwei Wu
Journal: International Journal of Distributed Sensor Networks
 (2016)

---

## Abstract

We propose a novel algorithm to optimize the energy efficiency (EE) of OFDM-based cognitive opportunistic relaying links (CORL) under secondary users (SUs) incorrectly sensing the unlicensed spectrum. We formulate an optimization problem with imperfect sensing that satisfies a specified power budget for the secondary users (SUs), while restricting the interference to primary user (PU) in a statistical manner. Unlike all related works in the literature, we consider the effect of subcarrier transmission mode on the relaying links and we additionally consider the effect of limited sensing capabilities of the SUs. The optimization problem is nonconvex and it is transformed to an equivalent problem using the concept of fractional programming. With the aid of the fractional programming method, an EE-oriented power allocation policy with low complexity is proposed which adopts the bisection method to speed up the search of the optimum. Simulation results show that the EE deteriorates as the channel sensing error increases. Comparisons with relevant works from the literature show that the EE is slightly deteriorated if the SU does not account for spectrum sensing errors.


---

## Body

## 1. Introduction

Cognitive radio (CR) networks have emerged as an efficient solution to the problem of spectrum scarcity and its underutilization. This is achieved by granting SUs’ opportunistic access to the white spaces within PUs’ spectrum while controlling the interference to PUs. Orthogonal frequency-division multiplexing access (OFDMA) is recognized as an attractive technique for CR due to its spectrum shaping flexibility, adaptively in allocating vacant radio resources, and capability of analyzing the spectral activities of PUs [1]. Incorporating cooperation into cognitive radio networks results in substantial performance gains in terms of spectrum efficiency (SE) for both PUs and SUs [2]. Besides the SE, the EE becomes a key issue for future wireless networks since energy cost imposes both financial and ecological burden on its development. EE power allocation especially is of crucial importance for cognitive relaying network design [3]. This is because that high EE is a basic premise for SUs to achieve high utilization of the limited transmit power which is consumed to not only improve SE but also implement some additionally important functionalities, for example, spectrum sensing.As a result, energy-efficient resource management has attracted attention in both industry and academia recently, especially for the OFDMA system which is the most popular modulation technique for current wireless networks. Different from the throughput-oriented optimization targets, energy-efficient resource management aims at maximizing the energy efficiency of the wireless system. One of the energy-efficiency metrics is called “bits-per-Joule,” which is defined as the system throughput to unit-energy consumption. For instance, the EE-maximization problem in an OFDMA system under a maximum total power constraint in frequency-selective channels is addressed [4]. In [5], the authors studied the tradeoff between EE and SE in the downlink of OFDMA networks. They showed that the EE is quasi-concave in the SE. Then based on this observation, a power and subcarrier allocation algorithm is proposed. In the uplink of an OFDMA system, the EE is addressed in [6]. Furthermore, in [7], the EE of two-way relaying was compared with those of the one-way relaying, showing that two-way relay transmission is not always more energy-efficient than one-way relay transmission. Nevertheless, [4–7] aim at maximizing the EE of system without taking the interference by SUs in CRNs into account. In [8], a method named as water-filling factor aided search (WFAS) was proposed to maximize EE under multiple constraints with perfect channel state information (CSI) at CR source, but relaying was not considered. In [9], the energy efficiency of an OFDM-based system is maximized, where multiple radio access technologies are employed for parallel transmission. In [10], an energy-efficient resource management scheme is developed for a downlink multiuser OFDM system with distributed antennas while considering proportional fairness among users. In our previous work [11], we propose an optimal power allocation scheme to maximize the EE of OFDMA opportunistic relay which is first proposed in [12] to better exploit the frequency-selective channels. However, [11] has not considered the peak primary PU’s interference constraints and spectrum sensing errors. Note that [13] also studied the EE optimization problem in CR system with imperfect spectrum sensing, but they all focus on the frame design including optimal sensing duration and data-transmission duration, as well as the optimal transmission power instead of the power allocation among each subcarrier. Besides, the authors in [14–16] analyze the EE performance of CRNs with imperfect spectrum sensing while relay strategies are not applied. Although a solution for EE-maximization problem in relay-aided CRN is proposed in [17], the authors only consider the ideal situation, that is, ignoring sensing errors.Based on research in CR relaying system, the motivation of this paper is expressed as follows. In order to further improve the SUs’ performance in terms of EE metric, we introduce the opportunistic DF (Decode and Forward) relaying strategy [12] to CR relay-aided networks to better exploit the frequency-selective channels, unlike [5, 6] where always-relaying protocol was considered. On the other hand, we assume that the SUs can coexist with the PU in the presence of both idle and busy sensing decisions while adapting their transmission power according to the imperfect sensing results, which differs from [15, 18]. In fact, the perfect spectrum sensing results are unavailable in practice, which makes the past research too idealistic to achieve feasible schemes for real system. Therefore, it is of great importance to study the energy-efficient resource management scheme under imperfect channel sensing.In this paper, adaptive power allocation is investigated to maximize EE for CORL with spectrum sensing errors considered in the system model. To achieve the optimal solution, the EE-maximization problem is simplified and transformed into an equivalent concave form, and then we use the Lagrangian technique to transform the equivalent problem into a corresponding dual problem. Finally, optimal allocation algorithm is proposed. To the best of our knowledge, adaptive EE power allocation is investigated to maximize EE for CORL considering spectrum sensing errors which has not been discussed in the literature. Our main contributions of this paper are summarized as follows:(i)
An EE optimization problem for CORL system with imperfect sensing results is established, subject to the individual power budget, peak PU’s interference constraints, and circuit power consumption in the total power expenditure, as well as considering subcarrier transmission mode. Particularly, our model can be easily extended to many practical scenarios with necessary modifications.(ii)
We probe into the optimal power allocation scheme with incorrect sensing in CORL system. On the basis of CORL model, we proposed a novel EE-oriented optimal power allocation iterative algorithm by exploiting the fractional programming and bisection method to completely solve the optimization problem, which reduces computation complexity significantly and yields a good tradeoff between EE and computational complexity.(iii)
Finally, extensive numerical simulation results corroborate our theoretical analysis and demonstrate the effectiveness of the proposed method. We found out that opportunistic relay protocol as compared to always relay-aided transmission protocol in CR networks is able to achieve higher performance in terms of EE metric. Furthermore, comparisons with relevant works from the literature show that the EE is slightly deteriorated if the SU does not account for spectrum sensing errors.The rest of this paper is organized as follows. Section2 introduces the system model considering the imperfect sensing. The EE power allocation optimization problem is analyzed, and we outline the proposed algorithm for its solution in Sections 3 and 4. Finally, simulation results are presented in Section 5, and conclusions are drawn in Section 6.


## 2. System Model and Problem Formulation

### 2.1. System Model

We consider a scenario where a two-hop OFDM-based CR system coexists with a PU in the same geographical location. As shown in Figure1(a), the system comprises one PU, one SU-transmitter (ST), one SU-relay (SR), and one SU-destination (SD) in the system. Let us denote the set of the PU’s bands K=1,2,…,K including the occupied subcarriers set KO and spectrum holes (unoccupied subcarriers by PU) set KU. Thus, we can obtain K=KO∪KU,KU≠Φ. Each of PU’s bands has a fixed bandwidth of Δf Hz. The opportunistic DF protocol in [12] which is used assists ST transmission to SD. It can be seen in Figure 1(b) that data frame structure for the considered CORL is different from the always relay-aided transmission protocol which is always idle for ST in the second slot. Specifically, every data-transmission session takes two consecutive equal-duration time slots (TS1, TS2) and OFDM with k∈KU subcarrier is used. In the first time slot, the ST radiates OFDM symbols using Ps,k as the transmit power for subcarrier k while the SR and SD receive. The ST-to-SD and ST-to-SR channel coefficients for subcarrier k are hsd,k and hsr⁡,k, respectively. In the second time slot, we define the subcarriers transmission mode indicator θk, which is a binary integer variable; that is, θk∈{0,1}. θk=1 represents relay transmission mode (RTM) which means that the SR retransmit OFDM symbols using Pr,k as the transmit power. The SR-to-SD channel coefficient is hrd⁡,k for subcarrier k. In band θk=0 represents direct transmission mode (DTM) which means transmission is solely undertaken by the ST in two successive time slots, and the SR is inactive for subcarrier k. Here, we define θ≜θk∣k∈KU to facilitate further description. Based on the two signaling intervals, the SD exploits maximum ratio combining (MRC) to retrieve the message. We further assume noise variance within one OFDM subcarrier to be σr2 at SR and σd2 at SD. According to the Shannon capacity formula, the secondary achievable data rate for DTM and RTM over subcarrier k can be, respectively, expressed as(1)Rsd⁡=log2⁡1+Ps,kGsd,k,(2)Rco⁡sd=12min⁡log2⁡1+Ps,kGsr⁡,k︸Rsr⁡,log2⁡1+Ps,kGsd,k+Pr,kGrd⁡,k︸Rco⁡srd,where Gsd,k=hsd,k2/σd2, Gsr⁡,k=hsr⁡,k2/σr2, and Grd⁡,k=hrd⁡,k2/σd2 denote the channel gain to noise ratio on the kth subchannel. Consequently, the achievable sum data rate for CORL can be derived as(3)RTPS,PR=∑k∈KU1-θkRsd⁡PS+∑k∈KUθkRco⁡sd⁡PS,PR,where PS,PR≜Ps,k,Pr,k∣k∈KU denotes the power allocation on ST and SR.Figure 1
System model. (a) Transmission process of CORL. (b) A frame structure for one particular subcarrier in CORL.
(a)
(b)


### 2.2. Interference with Spectrum Sensing Errors

In CR system, PU can access the licensed spectrum at any time and the probability of PU using subcarrierj is denoted by POj. Let Qfj and Qdj denote the false alarm and detection probability of subcarrier j, respectively. Since this paper is focused on the EE problem, we consider Qfj, Qdj independent to reduce the complexity. Therefore, Qfj and Qdj are given by [19](4)Qfj=Qεσ2-1τfs,∀j=KO∪KU,Qdj=Qεσ2-1-γτfs2γ+1,∀j=KO∪KU,where Qx=2/π∫x∞e-t2/2dt is the Gaussian tail probability, ε denotes a common threshold used across all subcarriers, σ2 denotes the noise power received at each SU, and τ is the sensing time used by the SU when sensing the primary behavior in each frame. fs is the sampling frequency during the sensing time and γ is the PU’s average SNR. Via (4), Qdj can be rewritten as(5)QdjQfj=Q12τ+1Q-1Qfj-γτfs,∀j=KO∪KU.LetβOOj be the posterior probability that the SU detects subcarrier as being used by PU which is indeed occupied. Using Bayes’ theorem and the law of total probability [20], βOOj can be derived as(6)βOOj=PrHj,1∣H^j,1=PrH^j,1∣Hj,1PrHj,1∑Δ=0,1PrH^j,Δ∣Hj,ΔPrHj,Δ=QdjPOjQfj1-POj+QdjPOj,where Hj,1 and Hj,0 represent the events that PU is active and idle on subcarrier j and H^j,1, H^j,0 are the sensing results that subcarrier j is occupied or unoccupied by PU, respectively. βUUj is the posterior probability that the evidence subcarrier j is really unoccupied given that SU senses it to be unoccupied, which can be expressed as(7)βUUj=PrHj,0∣H^j,0=PrH^j,0∣Hj,0PrHj,0∑Δ=0,1PrH^j,Δ∣Hj,ΔPrHj,Δ=1-Qfj1-POj1-Qdj1-POj+1-QfjPOj.Note that for perfect sensing βOOj=1 and βUUj=1. There exist two cases in which subcarrier k may introduce interference to PU. One is that subcarrier k is sensed correctly to be occupied by PU and the other is that subcarrier k is sensed incorrectly to be unoccupied by PU. Taking the above into account, the average interference introduced into the PU over subcarrier k with unit transmission power [21] can be written as(8)Ik=∑j∈KOβOOjIk,jθk+∑j∈KU1-βUUjIk,jθk,k∈KU,where Ik,jθk indicates the interference introduced into PU on subcarrier j when ST or SR transmits on subcarrier k with unit transmission power [22], and it can be expressed as (9)Ik,jθk=∫j-1Δf-k-1/2ΔfjΔf-k-1/2Δf1-θkGsp,k+θkGrp,kϕfdf,where Gsp,k and Grp,k are, respectively, denoted as the channel gain from ST to PU and SR to PU over subcarrier k, respectively. ϕf=TssinπfTs/πfTs2 represents the power spectral density (PSD) of OFDM transmitted signal, and Ts represents the duration of OFDM symbol.


### 2.3. Problem Formulation in CORL

The overall transmission power consumption in a unit frame contains the transmit power on ST and SR, which is calculated by(10)Ptr⁡=12∑k∈KUPs,k+12∑k∈KU1-θkPs,k+θkPr,k.To transmit data, the energy consumption consists of two parts: the energy consumption of power amplifier related to transmit power and the circuit energy consumption incurred by signal processing and active circuit blocks. We further assume that the circuit power consumption of equipment has nothing to do with the state of transmission system, and its average value is constant [13, 23]. In conclusion, the system total power consumption consists of circuit consumptions and overall transmit power. Therefore, considering power amplifier efficiency ρ∈0,1, the total circuit power consumption can be expressed as(11)PTC=PCS+PCR+1ρ∑k∈KU122-θkPs,k+θkPr,k,where PCS, PCR are denoted as the ST and SR circuit consumption which can be a constant. Like [14–18, 23], the EE measured by the “throughput-per-Joule” metric is defined as ratio of total throughput and total power; that is, EE(θ)≜RTPS,PR/PTC. Hence, maximizing the average EE metric for the CORL system can be written as (12)PS∗,PR∗=argmaxPS,PR⁡EEθPS,PR,where PS∗,PR∗≜Ps,k∗,Pr,k∗∣k∈KU represents the optimal power allocation on ST and SR.


## 2.1. System Model

We consider a scenario where a two-hop OFDM-based CR system coexists with a PU in the same geographical location. As shown in Figure1(a), the system comprises one PU, one SU-transmitter (ST), one SU-relay (SR), and one SU-destination (SD) in the system. Let us denote the set of the PU’s bands K=1,2,…,K including the occupied subcarriers set KO and spectrum holes (unoccupied subcarriers by PU) set KU. Thus, we can obtain K=KO∪KU,KU≠Φ. Each of PU’s bands has a fixed bandwidth of Δf Hz. The opportunistic DF protocol in [12] which is used assists ST transmission to SD. It can be seen in Figure 1(b) that data frame structure for the considered CORL is different from the always relay-aided transmission protocol which is always idle for ST in the second slot. Specifically, every data-transmission session takes two consecutive equal-duration time slots (TS1, TS2) and OFDM with k∈KU subcarrier is used. In the first time slot, the ST radiates OFDM symbols using Ps,k as the transmit power for subcarrier k while the SR and SD receive. The ST-to-SD and ST-to-SR channel coefficients for subcarrier k are hsd,k and hsr⁡,k, respectively. In the second time slot, we define the subcarriers transmission mode indicator θk, which is a binary integer variable; that is, θk∈{0,1}. θk=1 represents relay transmission mode (RTM) which means that the SR retransmit OFDM symbols using Pr,k as the transmit power. The SR-to-SD channel coefficient is hrd⁡,k for subcarrier k. In band θk=0 represents direct transmission mode (DTM) which means transmission is solely undertaken by the ST in two successive time slots, and the SR is inactive for subcarrier k. Here, we define θ≜θk∣k∈KU to facilitate further description. Based on the two signaling intervals, the SD exploits maximum ratio combining (MRC) to retrieve the message. We further assume noise variance within one OFDM subcarrier to be σr2 at SR and σd2 at SD. According to the Shannon capacity formula, the secondary achievable data rate for DTM and RTM over subcarrier k can be, respectively, expressed as(1)Rsd⁡=log2⁡1+Ps,kGsd,k,(2)Rco⁡sd=12min⁡log2⁡1+Ps,kGsr⁡,k︸Rsr⁡,log2⁡1+Ps,kGsd,k+Pr,kGrd⁡,k︸Rco⁡srd,where Gsd,k=hsd,k2/σd2, Gsr⁡,k=hsr⁡,k2/σr2, and Grd⁡,k=hrd⁡,k2/σd2 denote the channel gain to noise ratio on the kth subchannel. Consequently, the achievable sum data rate for CORL can be derived as(3)RTPS,PR=∑k∈KU1-θkRsd⁡PS+∑k∈KUθkRco⁡sd⁡PS,PR,where PS,PR≜Ps,k,Pr,k∣k∈KU denotes the power allocation on ST and SR.Figure 1
System model. (a) Transmission process of CORL. (b) A frame structure for one particular subcarrier in CORL.
(a)
(b)


## 2.2. Interference with Spectrum Sensing Errors

In CR system, PU can access the licensed spectrum at any time and the probability of PU using subcarrierj is denoted by POj. Let Qfj and Qdj denote the false alarm and detection probability of subcarrier j, respectively. Since this paper is focused on the EE problem, we consider Qfj, Qdj independent to reduce the complexity. Therefore, Qfj and Qdj are given by [19](4)Qfj=Qεσ2-1τfs,∀j=KO∪KU,Qdj=Qεσ2-1-γτfs2γ+1,∀j=KO∪KU,where Qx=2/π∫x∞e-t2/2dt is the Gaussian tail probability, ε denotes a common threshold used across all subcarriers, σ2 denotes the noise power received at each SU, and τ is the sensing time used by the SU when sensing the primary behavior in each frame. fs is the sampling frequency during the sensing time and γ is the PU’s average SNR. Via (4), Qdj can be rewritten as(5)QdjQfj=Q12τ+1Q-1Qfj-γτfs,∀j=KO∪KU.LetβOOj be the posterior probability that the SU detects subcarrier as being used by PU which is indeed occupied. Using Bayes’ theorem and the law of total probability [20], βOOj can be derived as(6)βOOj=PrHj,1∣H^j,1=PrH^j,1∣Hj,1PrHj,1∑Δ=0,1PrH^j,Δ∣Hj,ΔPrHj,Δ=QdjPOjQfj1-POj+QdjPOj,where Hj,1 and Hj,0 represent the events that PU is active and idle on subcarrier j and H^j,1, H^j,0 are the sensing results that subcarrier j is occupied or unoccupied by PU, respectively. βUUj is the posterior probability that the evidence subcarrier j is really unoccupied given that SU senses it to be unoccupied, which can be expressed as(7)βUUj=PrHj,0∣H^j,0=PrH^j,0∣Hj,0PrHj,0∑Δ=0,1PrH^j,Δ∣Hj,ΔPrHj,Δ=1-Qfj1-POj1-Qdj1-POj+1-QfjPOj.Note that for perfect sensing βOOj=1 and βUUj=1. There exist two cases in which subcarrier k may introduce interference to PU. One is that subcarrier k is sensed correctly to be occupied by PU and the other is that subcarrier k is sensed incorrectly to be unoccupied by PU. Taking the above into account, the average interference introduced into the PU over subcarrier k with unit transmission power [21] can be written as(8)Ik=∑j∈KOβOOjIk,jθk+∑j∈KU1-βUUjIk,jθk,k∈KU,where Ik,jθk indicates the interference introduced into PU on subcarrier j when ST or SR transmits on subcarrier k with unit transmission power [22], and it can be expressed as (9)Ik,jθk=∫j-1Δf-k-1/2ΔfjΔf-k-1/2Δf1-θkGsp,k+θkGrp,kϕfdf,where Gsp,k and Grp,k are, respectively, denoted as the channel gain from ST to PU and SR to PU over subcarrier k, respectively. ϕf=TssinπfTs/πfTs2 represents the power spectral density (PSD) of OFDM transmitted signal, and Ts represents the duration of OFDM symbol.


## 2.3. Problem Formulation in CORL

The overall transmission power consumption in a unit frame contains the transmit power on ST and SR, which is calculated by(10)Ptr⁡=12∑k∈KUPs,k+12∑k∈KU1-θkPs,k+θkPr,k.To transmit data, the energy consumption consists of two parts: the energy consumption of power amplifier related to transmit power and the circuit energy consumption incurred by signal processing and active circuit blocks. We further assume that the circuit power consumption of equipment has nothing to do with the state of transmission system, and its average value is constant [13, 23]. In conclusion, the system total power consumption consists of circuit consumptions and overall transmit power. Therefore, considering power amplifier efficiency ρ∈0,1, the total circuit power consumption can be expressed as(11)PTC=PCS+PCR+1ρ∑k∈KU122-θkPs,k+θkPr,k,where PCS, PCR are denoted as the ST and SR circuit consumption which can be a constant. Like [14–18, 23], the EE measured by the “throughput-per-Joule” metric is defined as ratio of total throughput and total power; that is, EE(θ)≜RTPS,PR/PTC. Hence, maximizing the average EE metric for the CORL system can be written as (12)PS∗,PR∗=argmaxPS,PR⁡EEθPS,PR,where PS∗,PR∗≜Ps,k∗,Pr,k∗∣k∈KU represents the optimal power allocation on ST and SR.


## 3. Problem Analysis on EE Power Allocation

To maximize the EE of the CORL network while guaranteeing that the interference to the PU receiver is maintained below a predefined threshold, we formulate a transmit power allocation optimization problem under some practical constraints. For simplicity, let us collect all indicator and power variables inθ, PS, and PR, respectively, and define D≜θ,PS,PR. Mathematically, we can formulate the EE-maximization problem for CORL as follows:(13)OP1:maxD∑k∈KU1-θklog2⁡1+Ps,kGsd⁡,k+1/2∑k∈KUθkmin⁡Rsr⁡,Rco⁡srdPCS+PCR+1/2ρ∑k∈KU2-θkPs,k+θkPr,k,(14)s.t.∑k∈KUPs,k≤PSmax,(15)∑k∈KUPr,k≤PRmax,(16)∑k∈KUPs,kIkθk≤IPth,(17)∑k∈KU1-θkPs,k+θkPr,kIkθk≤IPth,(18)θk∈0,1,k=KU,(19)Ps,k≥0,Pr,k≥0,k∈KU,where PSmax and PRmax are denoted as the individual power limitations at ST and SR, respectively. IPth signifies the maximum interference power threshold prescribed by the PU. Constraints (16) and (17) in OP1 ensure that interference to PU should be restricted by a specified threshold to prevent the PU from severe performance degradation in TS1 and TS2. In its current form (13), it is obvious that the joint optimization problem OP1 is a nonconvex mixed-integer nonlinear program (MINLP) which is known to be NP-hard. However, the aim of this work is to maximize the EE metric of (13) subject to the individual power budget and peak PU’s interference constraints. According to the idea of subcarrier transmission mode indicator in [12], we introduced a straightforward method for CORL system for which the subcarrier k is selected RTM if Gsd,k≤Gsr⁡,k and Gsd,k≤Grd⁡,k. Otherwise, the DTM offers a better capacity. In what follows, we denoted two sets SDT and SRT to represent DTM and RTM, respectively, which are defined as follows:(20)SDT=k∣Gsd⁡,k>Gsr⁡,korGsd⁡,k<Gsr⁡,k,Gsd⁡,k>Grd⁡,k,SRT=k∣Gsd⁡,k≤Gsr⁡,k,Gsd⁡,k≤Grd⁡,k.⁡Theorem 1.
In a three-node DF relaying link network, when the system EE achieves maximum, each subcarrier rate for two hops should be equal. Namely, the most economical choice isRsdPs,k=RcosrdPs,k,Pr,k.Proof.
Suppose that the system EE achieves optimum; there existsk-, satisfying Rsd⁡Ps,k-≠Rco⁡srdPs,k-,Pr,k-. Because each subcarrier rate for two hops depends on the smaller one, this means one can reduce the transmit power of the larger rate hop in Rsd⁡Ps,k- and Rco⁡srdPs,k-,Pr,k-, leading to Rsd⁡Ps,k-=Rco⁡srdPs,k-,Pr,k-. By the definition of EE function, reducing transmission power of the larger rate hop will increase system EE under the condition that the system rate is a constant [4]. Obviously, the original hypothesis leads to contradiction. Therefore, the original proposition is true.According to Theorem1, we have the following relationship: (21)Ps,kGsr⁡,k=Ps,kGsd⁡,k+Pr,kGrd⁡,kifθk=1.Then Pr,k=χkPs,k, where χk=Gsr⁡,k-Gsd⁡,k/Grd⁡,k. Based on this classification, the equivalent problem of OP1 can be reformulated as(22)OP2:maxPSRTPSPCS+PCR+1/2ρ∑k∈SDT2Ps,k+∑k∈SRT1+χkPs,k,(23)s.t.∑k∈SDTPs,k+∑k∈SRTPs,k≤PSmax,(24)∑k∈SRTPs,k≤PRmax,(25)∑k∈SDTPs,kIk+∑k∈SRTPs,kIk≤IPth,(26)∑k∈SRT1+χkPs,kIk≤IPth,(27)Ps,k≥0,k∈SDT∪SRT,where RTPS represents the capacity for CORL system. It can be expressed as(28)RTPS=∑k∈SDTRsd⁡PS+∑k∈SRTRco⁡sd⁡PS.From OP2, we observe that constraints are either linear or convex, but objective function equation (22) is not a concave function. Actually, the problem of OP2 belongs to the quasi-concave programming, which has been proved in our previous work [11]. In the next section, we will show that we can obtain optimal solution of EE-maximization problem by exploiting special structure of the objective function. To this end, the monotonically increasing and strictly concave characteristic of the numerator RTPS in (22) is summarized in Theorem 2.Theorem 2.
Givenθ, RTPS for CORL is monotonically increasing and strictly concave with respect to PS.Proof.
From (3), RTPS can be expressed as(29)RTPS=1-θRsd⁡PS+θRco⁡sd⁡PS.According to [24], it is easy to know that Rsd⁡, Rsr⁡, and Rco⁡srd are monotonically increasing and strictly concave with respect to PS. On the other hand, we observe that the subcarriers transmission mode indicator θ=θk∈{0,1}∣k∈KU≥0 is defined as nonnegative integers, so we only need to prove that the second item in (29), that is, Rco⁡sd⁡PS, is monotonically increasing and strictly concave with respect to PS. Considering the relationships in (21), we can rewrite (3) as Rco⁡sd⁡PS=minRsr⁡PS,Rco⁡srd(χkPS). Then we have(30)Rco⁡sd⁡ξp1+1-ξp2=min⁡Rsr⁡ξp1+1-ξp2,Rco⁡srdξχ1p1+1-ξχ2p2>min⁡ξRsr⁡p1+1-ξRsr⁡p2,ξRco⁡srdχ1p1+1-ξRco⁡srdχ2p2≥ξmin⁡Rsr⁡p1,Rco⁡srdχ1p1+1-ξmin⁡Rsr⁡p2,Rco⁡srdχ2p2=ξRco⁡sd⁡p1+1-ξRco⁡sd⁡p2,where 0≤ξ≤1, ∀p1,p2∈dom⁡Rco⁡sd⁡(PS). Hence, RTPS is monotonically increasing and strictly concave with respect to PS.


## 4. Adaptive Power Allocation to Maximize EE

It is noticeable that objective function equation (22) in OP2 is not a concave function, and the solution for it will be of high complexity. In this section, we first use fraction programming [25] to transform the problem into an equivalent convex optimization problem and then use the Lagrangian technique to transform the equivalent problem into a corresponding dual problem. Subsequently, an optimal iterative algorithm is proposed.


### 4.1. Adaptive Power Allocation

From Theorem2, we follow the fact that the numerator RTPS in (22) is a concave function, and the denominator of (22) is affine function of SUs’ power. Besides, all the constraints are convex set. Inspired by Dinkelbach’s algorithm in [25], we can transform this problem into a parameterized convex maximization problem. Primarily, a new objective function is defined as(31)TPS,q=∑k∈SDTRsd⁡PS+∑k∈SRTRco⁡sd⁡PS-q12ρ∑k∈SDT2Ps,k+∑k∈SRT1+χkPs,k+PCS+PCR,where q is a positive parameter and can be interpreted as a pricing factor for SUs’ power consumption. Hence, another problem is formulated as(32)OP3:maxPS,qTPS,q,s.t.25,26,27,28,29.Let S denote the feasible region of OP2 and OP3. Define Fq=maxPSTPS,q∣PS∈S as the maximum value of OP3 with each fixed q. Then, the optimal value and solution of OP3 can be defined as(33)fq=argmaxPSTPS,q∣PS∈S.The following lemma introduced by Dinkelbach’s algorithm [25] can relate OP2 and OP3, and the detailed proof of Lemma 3 can also be found in [25].Lemma 3.
The optimal solutionPS∗ achieves the optimal value q∗ of OP2, that is, q∗=maxPSTPS,q∣PS∈S=EE(θ)(PS∗), if and only if(34)Fq∗,PS∗=maxPS⁡TPS,q∗∣PS∈S=0,fq∗=PS∗.This lemma indicates that, at the optimal parameterq∗, the optimal solution to OP3 is also the optimal solution to OP2. Hence, solving OP2 can be realized by finding the optimal power allocation of OP3 for a given q and then update q until (31) is established. For a given q, the optimal power allocation can be obtained using convex theory [24] because of the convex characteristic of OP3. Hence, the existing water-filling power allocation approach gives the solution to it [26]. However, besides adapting the power allocation on all subcarriers, we need to consider subcarrier transmission mode. The Lagrange function for OP3 is constructed as(35)LPS,q,λ=-TPS,q+λ1∑k∈SDTPs,k+∑k∈SRTPs,k-PSmax+λ2∑k∈SRTPs,k-PRmax+λ3∑k∈SDTPs,kIk+∑k∈SRTPs,kIk-IPth+λ4∑k∈SRT1+χkPs,kIk-IPth,where λ1, λ2, λ3, and λ4 are the nonnegative Lagrangian dual variables for constraints (26)–(30), respectively. And the dual problem of OP3 is given by(36)D=minλii=14maxq,PS⁡Lq,Ps,k∣k∈KU,λii=14s.t.λii=14≥0.Using the Karush-Kuhn-Tucker (KKT) conditions [24], we can deduce the optimal power allocation of problem OP3, which can be written as(37)Ps,k∗=ρlog2⁡eq+ρξk-1Gsd⁡,k+,∀k∈SDT,(38)Ps,k∗=ρlog2⁡eq1+χk+ρϑk-1Gsr⁡,k+,∀k∈SRT,where ξk=λ1+λ3Ik and ϑk=2λ1+λ2+λ3Ik+λ4Ik1+χk and x+ denotes max0,x. In addition, subcarrier k is used for RTM communication; that is, k∈SRT. The corresponding relay transmission power at SR can be expressed as(39)Pr,k∗=Gsr⁡,k-Gsd⁡,kρlog2⁡eqGrd⁡,k+Gsr⁡,k-Gsd⁡,k+ρϑkGrd⁡,k-Gsr⁡,k-Gsd⁡,kGrd⁡,kGsr⁡,k+,∀k∈SRT.It can be found that the traditional water-filling method could not solve the primal problem directly, because the denominator of (37) and (38) contain a linear combination of the Lagrangian dual factors. Thus the incremental-update based subgradient method in [27] is introduced to derive the optimal dual factors λ∗=λi∗i=14 for power and interference constraints. For our problems, the corresponding iterative update λ is based on the following iteration procedure:(40)λ1l+1=λ1l-μlPSmax-∑k∈SDTPs,k-∑k∈SRTPs,k+λ2l+1=λ2l-μlPRmax-∑k∈SRTPs,k+λ3l+1=λ3l-μlIPth-∑k∈SDTPs,kIk-∑k∈SRTPs,kIk+λ4l+1=λ4l-μlIPth-∑k∈SRT1+χkPs,kIk+,where l refers to the iteration index and μl>0 denotes a sufficiently small positive step size for the lth iteration, and it is a sequence of step size which is defined in many types in [27]. It should be mentioned that small step size leads to slow convergence. Besides, each element of the gradient depends on the corresponding subcarrier’s channel gain, potentially differing from each other by orders of magnitude. Hence, a line search of the optimal step size needs to cover a large range to ensure global convergence on all subcarriers, which is computationally expensive. Therefore, in order to find the optimal step size, like [4], we define flu=λ[l]-u∇Lλ[l]+, which has also been proved to be concave in μl, and can quickly obtain the optimal step value of ul∗ by using bisection search algorithm summarized in [4]. After we obtain the optimal power assignment PS,PR for the given q, update q=EE(θ)PS,PR until Lemma 3 is satisfied; thus the equivalent optimal power assignment PS∗,PR∗ of OP3 can be achieved in the end. So, jointing the fractional programming and bisection method, an EE-oriented power allocation iterative optimization algorithm for CORL called ECORL is provided, which is described in Algorithm 1 to solve the power allocation of programming OP1.Algorithm 1:ECORL.
Input: the error tolerance ς>0, δ>0 and the maximum iteration number MaxIter;
Output: optimal EE power allocation policy sets of OP1;
Initialize maximum EE q0=q0, the iteration index n=0 and l=0, dual variables λ0=λinitial;
Compute  Gsr⁡,k,Gsd⁡,k,Grd⁡,k∣∀k∈KU, Then
Obtain subcarrier transmutation modeθ using (20);
While  (Fqn>ςandn≤MaxIter)  do
Updateλ via the subgradient method as follow:
Repeat
Step1 Compute Ps,kl, ∀k∈SDT and Ps,kl, ∀k∈SRT through (37) and (38) respectively;
Step2 Find the optimal step size ul∗ by using the bisection algorithm described in [27];
Step3 Update λ according to (40);
Step4 Calculate Δλ=Δλl+1-Δλl and Set l=l+1;
Until2-norm Δλ<δ (λ converges)
Setn=n+1 and Obtain PS through (37) and (38);
CalculatePR on SR by the corresponding (21);
Updateqn=EE(θ)PS,PR through (22);
End while
Return  q∗=qn, the ξ-optimal solution PS∗,PR∗ of OP1 and calculate the EE and capacity for CROL networks.Remark 4.
Note that, in the case ofq=0, EE-maximization problem is equivalent to SE-maximization. Consequently, for given maximum iteration number MaxIter and error tolerance ς and δ, the optimal EE and SE power allocation policy of OP1 can be easily obtained by ECORL, which will be validated by the simulation in Section 5.


### 4.2. Coverage and Complexity Analysis

The proposed algorithm ECORL will always converge to optimum provided by Theorem5, meaning that, for every ς, the power policy set PS that maximizes TPS,q is found. The algorithm execution stops if q is zero or less than the ς value.Theorem 5.
The ECORL will always converge to optimum.Proof.
To prove the convergence of ECORL, we letTPS,q=RTPS-qEPS, for ease of presentation, where EPS denotes the energy consumption in CORL system and can be written as(41)EPS=12ρ∑k∈SDT2Ps,k+∑k∈SRT1+χkPs,k+PCS+PCR.
Before going into the proof of convergence, we will first discuss an observation thatFq is decreasing in q. For the sake of generality, we assume q1<q2  q1,q2∈S; then(42)Fq2=maxPS⁡RTPS-q2EPS∣PS∈S=RTPS∗-q2EPS∗<RTPS∗-q1EPS∗<maxPS⁡RTPS-q1GPS∣PS∈S=Fq1.From (42), we know that Fq is decreasing in q. To prove the convergence of ECORL, we only need to show that Fq becomes smaller than ς with the number of iterations. We now show that q is nonincreasing in successive iterations of the algorithm. We have(43)TPS∗n-1,qn≤max⁡TPs∗n-1,qn∣PS∈S=Fqn=RTPS∗n-qnEPS∗n=qn+1EPS∗n-qnEPS∗n=qn+1-qnEPS∗n.Now, it follows that qn+1≥qn, because GPS∗n>0. From the above observation, we have that Fq is decreasing in q and q is nonincreasing in successive iterations of ECORL. Therefore, Fq is nonincreasing in successive iterations. Furthermore, Fqn does become zero from Lemma 3, which follows that fact that Fqn does become smaller than ς. Therefore, the ECORL algorithm will always converge to optimum.In addition, the ECORL has two iterations: one is for the total individual power constraint and the other is for the Lagrange multiplierλ for the total power constraint. The complexity of ECORL can be expressed as ONfNsNc, where Nf, Ns, respectively, denote the iteration numbers of the while-loop and repeat-loop in ECORL and Nc represent the number of subcarrier realizations. Nf and Ns especially rely on the choice of the corresponding error tolerance ς and δ, respectively. Hence, the EE optimization problem can be solved by the proposed algorithm with a considerable saving in terms of computational complexity, which is very applicable for the physical life.


## 4.1. Adaptive Power Allocation

From Theorem2, we follow the fact that the numerator RTPS in (22) is a concave function, and the denominator of (22) is affine function of SUs’ power. Besides, all the constraints are convex set. Inspired by Dinkelbach’s algorithm in [25], we can transform this problem into a parameterized convex maximization problem. Primarily, a new objective function is defined as(31)TPS,q=∑k∈SDTRsd⁡PS+∑k∈SRTRco⁡sd⁡PS-q12ρ∑k∈SDT2Ps,k+∑k∈SRT1+χkPs,k+PCS+PCR,where q is a positive parameter and can be interpreted as a pricing factor for SUs’ power consumption. Hence, another problem is formulated as(32)OP3:maxPS,qTPS,q,s.t.25,26,27,28,29.Let S denote the feasible region of OP2 and OP3. Define Fq=maxPSTPS,q∣PS∈S as the maximum value of OP3 with each fixed q. Then, the optimal value and solution of OP3 can be defined as(33)fq=argmaxPSTPS,q∣PS∈S.The following lemma introduced by Dinkelbach’s algorithm [25] can relate OP2 and OP3, and the detailed proof of Lemma 3 can also be found in [25].Lemma 3.
The optimal solutionPS∗ achieves the optimal value q∗ of OP2, that is, q∗=maxPSTPS,q∣PS∈S=EE(θ)(PS∗), if and only if(34)Fq∗,PS∗=maxPS⁡TPS,q∗∣PS∈S=0,fq∗=PS∗.This lemma indicates that, at the optimal parameterq∗, the optimal solution to OP3 is also the optimal solution to OP2. Hence, solving OP2 can be realized by finding the optimal power allocation of OP3 for a given q and then update q until (31) is established. For a given q, the optimal power allocation can be obtained using convex theory [24] because of the convex characteristic of OP3. Hence, the existing water-filling power allocation approach gives the solution to it [26]. However, besides adapting the power allocation on all subcarriers, we need to consider subcarrier transmission mode. The Lagrange function for OP3 is constructed as(35)LPS,q,λ=-TPS,q+λ1∑k∈SDTPs,k+∑k∈SRTPs,k-PSmax+λ2∑k∈SRTPs,k-PRmax+λ3∑k∈SDTPs,kIk+∑k∈SRTPs,kIk-IPth+λ4∑k∈SRT1+χkPs,kIk-IPth,where λ1, λ2, λ3, and λ4 are the nonnegative Lagrangian dual variables for constraints (26)–(30), respectively. And the dual problem of OP3 is given by(36)D=minλii=14maxq,PS⁡Lq,Ps,k∣k∈KU,λii=14s.t.λii=14≥0.Using the Karush-Kuhn-Tucker (KKT) conditions [24], we can deduce the optimal power allocation of problem OP3, which can be written as(37)Ps,k∗=ρlog2⁡eq+ρξk-1Gsd⁡,k+,∀k∈SDT,(38)Ps,k∗=ρlog2⁡eq1+χk+ρϑk-1Gsr⁡,k+,∀k∈SRT,where ξk=λ1+λ3Ik and ϑk=2λ1+λ2+λ3Ik+λ4Ik1+χk and x+ denotes max0,x. In addition, subcarrier k is used for RTM communication; that is, k∈SRT. The corresponding relay transmission power at SR can be expressed as(39)Pr,k∗=Gsr⁡,k-Gsd⁡,kρlog2⁡eqGrd⁡,k+Gsr⁡,k-Gsd⁡,k+ρϑkGrd⁡,k-Gsr⁡,k-Gsd⁡,kGrd⁡,kGsr⁡,k+,∀k∈SRT.It can be found that the traditional water-filling method could not solve the primal problem directly, because the denominator of (37) and (38) contain a linear combination of the Lagrangian dual factors. Thus the incremental-update based subgradient method in [27] is introduced to derive the optimal dual factors λ∗=λi∗i=14 for power and interference constraints. For our problems, the corresponding iterative update λ is based on the following iteration procedure:(40)λ1l+1=λ1l-μlPSmax-∑k∈SDTPs,k-∑k∈SRTPs,k+λ2l+1=λ2l-μlPRmax-∑k∈SRTPs,k+λ3l+1=λ3l-μlIPth-∑k∈SDTPs,kIk-∑k∈SRTPs,kIk+λ4l+1=λ4l-μlIPth-∑k∈SRT1+χkPs,kIk+,where l refers to the iteration index and μl>0 denotes a sufficiently small positive step size for the lth iteration, and it is a sequence of step size which is defined in many types in [27]. It should be mentioned that small step size leads to slow convergence. Besides, each element of the gradient depends on the corresponding subcarrier’s channel gain, potentially differing from each other by orders of magnitude. Hence, a line search of the optimal step size needs to cover a large range to ensure global convergence on all subcarriers, which is computationally expensive. Therefore, in order to find the optimal step size, like [4], we define flu=λ[l]-u∇Lλ[l]+, which has also been proved to be concave in μl, and can quickly obtain the optimal step value of ul∗ by using bisection search algorithm summarized in [4]. After we obtain the optimal power assignment PS,PR for the given q, update q=EE(θ)PS,PR until Lemma 3 is satisfied; thus the equivalent optimal power assignment PS∗,PR∗ of OP3 can be achieved in the end. So, jointing the fractional programming and bisection method, an EE-oriented power allocation iterative optimization algorithm for CORL called ECORL is provided, which is described in Algorithm 1 to solve the power allocation of programming OP1.Algorithm 1:ECORL.
Input: the error tolerance ς>0, δ>0 and the maximum iteration number MaxIter;
Output: optimal EE power allocation policy sets of OP1;
Initialize maximum EE q0=q0, the iteration index n=0 and l=0, dual variables λ0=λinitial;
Compute  Gsr⁡,k,Gsd⁡,k,Grd⁡,k∣∀k∈KU, Then
Obtain subcarrier transmutation modeθ using (20);
While  (Fqn>ςandn≤MaxIter)  do
Updateλ via the subgradient method as follow:
Repeat
Step1 Compute Ps,kl, ∀k∈SDT and Ps,kl, ∀k∈SRT through (37) and (38) respectively;
Step2 Find the optimal step size ul∗ by using the bisection algorithm described in [27];
Step3 Update λ according to (40);
Step4 Calculate Δλ=Δλl+1-Δλl and Set l=l+1;
Until2-norm Δλ<δ (λ converges)
Setn=n+1 and Obtain PS through (37) and (38);
CalculatePR on SR by the corresponding (21);
Updateqn=EE(θ)PS,PR through (22);
End while
Return  q∗=qn, the ξ-optimal solution PS∗,PR∗ of OP1 and calculate the EE and capacity for CROL networks.Remark 4.
Note that, in the case ofq=0, EE-maximization problem is equivalent to SE-maximization. Consequently, for given maximum iteration number MaxIter and error tolerance ς and δ, the optimal EE and SE power allocation policy of OP1 can be easily obtained by ECORL, which will be validated by the simulation in Section 5.


## 4.2. Coverage and Complexity Analysis

The proposed algorithm ECORL will always converge to optimum provided by Theorem5, meaning that, for every ς, the power policy set PS that maximizes TPS,q is found. The algorithm execution stops if q is zero or less than the ς value.Theorem 5.
The ECORL will always converge to optimum.Proof.
To prove the convergence of ECORL, we letTPS,q=RTPS-qEPS, for ease of presentation, where EPS denotes the energy consumption in CORL system and can be written as(41)EPS=12ρ∑k∈SDT2Ps,k+∑k∈SRT1+χkPs,k+PCS+PCR.
Before going into the proof of convergence, we will first discuss an observation thatFq is decreasing in q. For the sake of generality, we assume q1<q2  q1,q2∈S; then(42)Fq2=maxPS⁡RTPS-q2EPS∣PS∈S=RTPS∗-q2EPS∗<RTPS∗-q1EPS∗<maxPS⁡RTPS-q1GPS∣PS∈S=Fq1.From (42), we know that Fq is decreasing in q. To prove the convergence of ECORL, we only need to show that Fq becomes smaller than ς with the number of iterations. We now show that q is nonincreasing in successive iterations of the algorithm. We have(43)TPS∗n-1,qn≤max⁡TPs∗n-1,qn∣PS∈S=Fqn=RTPS∗n-qnEPS∗n=qn+1EPS∗n-qnEPS∗n=qn+1-qnEPS∗n.Now, it follows that qn+1≥qn, because GPS∗n>0. From the above observation, we have that Fq is decreasing in q and q is nonincreasing in successive iterations of ECORL. Therefore, Fq is nonincreasing in successive iterations. Furthermore, Fqn does become zero from Lemma 3, which follows that fact that Fqn does become smaller than ς. Therefore, the ECORL algorithm will always converge to optimum.In addition, the ECORL has two iterations: one is for the total individual power constraint and the other is for the Lagrange multiplierλ for the total power constraint. The complexity of ECORL can be expressed as ONfNsNc, where Nf, Ns, respectively, denote the iteration numbers of the while-loop and repeat-loop in ECORL and Nc represent the number of subcarrier realizations. Nf and Ns especially rely on the choice of the corresponding error tolerance ς and δ, respectively. Hence, the EE optimization problem can be solved by the proposed algorithm with a considerable saving in terms of computational complexity, which is very applicable for the physical life.


## 5. Simulation Results

In this section, we present numerical results to illustrate the performance of the proposed EE-maximizing power adaptation methods. Consider the scenario as the following: the location of ST is(0,0)m and the coordinates of PU, SR, and SD are fixed on (500,288)m, (500,-288)m, and (1000,0)m, respectively. We consider the CR relay-aided system consists of K=64 subcarriers; that is, K=1,2,…,64, the probability of spectrum holes ≥50%. The bandwidth of each subcarrier Δf and the duration of the OFDM symbol Ts in TS1, TS2 are assumed as 0.3125 MHz and 4 us, respectively. Further define circuit power consumption PCS=3 W, PCR=2 W, fs=100 MHz, τ=2 ms, γ=-5 dB, ρ=0.5, PSmax=PRmax=Ptr⁡th, and the noise variance σd2=σr2=σp2=10-6W. In the primary network, the central location of the PU is defined as c=16. Assume the activity probability POj equals 0.8·exp-j-c2/50 [14]. For ease of analysis, the detection Qdj is uniformly distributed over 0.90,0.99. Channel complex gains hk are picked from a Rayleigh fading channel with the following distribution: hk=CN0,1/L1+dα [12], where the path loss exponent α is 3, distance is dm, and the number of taps L is 16. Finally, given hk2, σd2, σr2, and σp2 the channel gains-to-noise-power ratios Gsd⁡,k, Gsr⁡,k, Grd⁡,k, Gsp,k, and Grp,k are calculated as described in Section 2.Figures2 and 3, respectively, present the power allocation for maximizing EE (EE-Max) and maximizing SE (SE-Max) versus each subcarrier link when running ECORL algorithm. As shown in Figures 2 and 3, the notation (×) indicates the subcarrier is occupied by PU, and the notations (+) and (□) at the top of the figure signify the opportunistic relay link transmission mode; that is, θk=1 and θk=0, respectively. It is needed to mention that the corresponding power on SR has to be used over two successive time slots (value shown by the solid curve) when notation (□) is active. It can be shown that the individual power budget is split among ST power (solid line) and SR power (dotted line) appropriately and effectively under the peak PU’s interference constraints, which demonstrates that ECORL algorithm has excellent performance on power allocation both for EE-Max and for SE-Max scheme. It is needed to mention that the EE-Max problem is equivalent to SE-Max when we set q=0 in ECORL algorithm, which confirms the conclusions in Remark 4. Moreover, in contrast to optimal power allocation between EE-Max and SE-Max as shown in Figures 2 and 3, it is clear that the EE-Max scheme on each subcarrier link is far smaller than SE-Max scheme, for the reason that the SE-Max scheme greedily grows with the interference tolerance all the time until the transmit power is used up. We also found that we should use low power transmission to guarantee high EE rather than augment transmission power budget. Consequently, it validates that the proposed method is better for EE not for SE, which is valuable in physical life.Figure 2
Power allocation on SU’s each subcarrier link with maximizing EE considering sensing errors.Figure 3
Power allocation on SU’s each subcarrier link for maximizing SE considering sensing errors.Figure4 demonstrates the EE-Max and SE-Max solved by the proposed ECORL algorithm with/without considering sensing errors (cse/wcse) with respect to Ptr⁡th, where the interference threshold IPth is 1×10-3W. Here, the maximizing EE without considering sensing errors is obtained by the proposed method in [18] to facilitate comparison. In Figure 4, as Ptr⁡th increases, the EE is first increasing and then decreasing, because when Ptr⁡th becomes larger, the EE performance is subject not to the individual power constraint but to the interference constraint. Besides, the interference constraint is first gradually bound and then strictly bound. From the figure, we observe that the EE by the proposed method first is the same as that by maximizing SE, while it is larger than the latter when transmitted power goes larger. Because system rate R is a logarithmic function of transmission power Ptr⁡ and it will consume much more transmission power Ptr⁡ to further improve R when the capacity is greater than a certain degree, thus the EE decreases. It also conveys that the EE without sensing errors is better than that with sensing errors. More importantly, we found that the opportunistic relay protocol as compared to always relay-aided transmission protocol is able to effectively improve performance in terms of EE metric. This is because introducing the opportunistic DF relaying strategy into CR relay-aided networks can better exploit the frequency-selective channels.Figure 4
System EE versus individual power budgetPtr⁡th.In Figure5, we illustrate the EE versus the interference threshold IPth with/without considering sensing errors (cse/wcse) under different individual total transmit power. It can be observed that the greater the total individual transmit power, the higher the EE. When the transmit power is relatively high, for example, Ptr⁡th=10 W, the EE performance is mainly decided by the interference threshold. However, when the transmit power is low enough, for example, Ptr⁡th=3 W, the EE is constrained by the total transmit power and will be constant as the interference threshold increases. Also, the figure shows that the performance of ECORL which has taken spectrum sensing errors into consideration has a reasonable loss compared to that without considering sensing errors depending on the value of the total transmit power. When the interference IPth constraints are relatively small, the EE achieved without considering sensing errors [18] is about 6.5% larger than that gained considering sensing errors. This is due to the fact that the strategy proposed in [18] is EE-Max with the total power and interference constraints, and it does not consider sensing errors. Besides, we found the error gap caused by imperfect sensing will become smaller as IPth increases.Figure 5
System EE versus the interference thresholdIPth with perfect and imperfect sensing.Figure6 displays the system EE with respect to βUUj with/without considering sensing errors (cse/wcse). For simplicity, we assume a common reliability of decision βUU instead of βUUj. Clearly, the system EE increases as βUU increases, for the reason that SU can obtain more chance to access spectrum as βUU increases. This is due to the fact that SU gains more access probability, and the dominated interference to the primary system decreases. Hence, the improved throughput, to some extent, leads to larger energy efficiency. Besides, the system EE considering sensing errors is not better than that without considering case because the more interference resulted from PU. It also explains that the sensing errors can weaken the performance of the secondary system, and it is important to improve the accuracy of spectrum sensing. Figure 7 demonstrates the relationship between the normalized system EE and the number of iterations. Here we assume the channel gain of each subcarrier has Rayleigh distribution with a unit average. The system EE is normalized by the optimal value. It can be shown that the sequence of iterations produces a monotonically increasing objective and always converges. After at most six iterations, the achieved normalized system EE converges to the final optimal value. And the normalized system EE has been up to 90% of the final optimal value on the average after only four iterations; meanwhile, the values of all the samples are more than 80% of the final optimal value. So we can take four as the maximum number of iterations in the practical design to greatly reduce the complexity of the algorithm. Besides, we can see that ECORL algorithm converges very fast to the global optimum.Figure 6
System EE versusβUU with perfect and imperfect sensing.Figure 7
Convergence rate of ECORL: relationship between normalized system EE and iterations.Finally, we consider the scenario that the coordinate of SR moves from point(0,-288) m to point (1000,-288) m with step 200 m; that is, PosxSR∈0,1000 m. In Figure 8, we depict the system EE versus different PosxSR for ECORL under the three special transmission mode cases (such as CORL, always relay, and no relay) with Ptr⁡th=10W and IPth=1×10-3W. We first find out that the EE curve of CORL system is relatively flat and always higher than that of other cases in CR networks. This accounts for the fact CORL system uses opportunistic DF relaying, which better exploits the flexibility of transmission mode selection for the power reduction. The second case is the always cooperative transmission; that is, RTM is used at each subcarrier. We can observe that the curve first increases and then decreases with PosxSR increasing. Meanwhile, the system EE will increase as PosxSR increases and achieve their maximum EE when SR lies in the middle range both (such as PosxSR∈400,600m) in CORL and in always relay cases. The third one is the noncooperative transmission; that is, DTM is used at each subcarrier. The EE curve which is first increasing and then decreasing is quite different from the previous two cases. It is clear that DTM is more energy-efficient than RTM when the SR lies in the point of (0,-288)m or (1000,-288) m, for the reason that it is more likely to have one of Gsr⁡,k or Grd⁡,k be much smaller than Gsd,k, which can explain the observation. Additionally, we also see that the sensing errors can weaken the performance of the secondary transmission, and it is important to improve the accuracy of spectrum sensing.Figure 8
System EE versus differentPosxSR.


## 6. Conclusion

In this paper, we have studied the resource allocation problem for EE power allocation in OFDM-based cognitive opportunistic relay links with spectrum sensing errors considered. To maximize the EE of the SUs under joint individual transmit power and interference constraints, we proposed an optimal power allocation algorithm using equivalent conversion and transform the equivalent problem into a corresponding Lagrangian dual problem. The simulation results show that when imperfect spectrum sensing is not taken into account, excessive interference will be introduced to PU; however, the EE is about6.5% larger than that obtained by ECORL scheme. Meanwhile, the proposed strategy can improve EE significantly compared to the always relay-aided scheme in CR networks, and the proposed algorithm exhibits a good convergence performance both theoretically and in simulation analysis. In the future, the energy-efficient resource allocation problems for more complicated green cognitive radio networks (e.g., multiuser scenario with imperfect channel state information) should be considered.



---

*Source: 2716042-2016-02-14.xml*