Energy-Efficiency Power Allocation for Cognitive Radio MIMO-OFDM Systems
Energy-Efficiency Power Allocation for Cognitive Radio MIMO-OFDM Systems
ETRI Journal. 2014. Aug, 36(4): 686-689
Copyright © 2014, Electronics and Telecommunications Research Institute(ETRI)
  • Received : August 27, 2013
  • Accepted : February 20, 2014
  • Published : August 01, 2014
Export by style
Cited by
About the Authors
Jiakuo, Zuo
Van Phuong, Dao
Yongqiang, Bao
Shiliang, Fang
Li, Zhao
Cairong, Zou

This paper studies energy-efficiency (EE) power allocation for cognitive radio MIMO-OFDM systems. Our aim is to minimize energy efficiency, measured by “Joule per bit” metric, while maintaining the minimal rate requirement of a secondary user under a total power constraint and mutual interference power constraints. However, since the formulated EE problem in this paper is non-convex, it is difficult to solve directly in general. To make it solvable, firstly we transform the original problem into an equivalent convex optimization problem via fractional programming. Then, the equivalent convex optimization problem is solved by a sequential quadratic programming algorithm. Finally, a new iterative energy- efficiency power allocation algorithm is presented. Numerical results show that the proposed method can obtain better EE performance than the maximizing capacity algorithm.
I. Introduction
Energy-efficient (EE) wireless communications have great ecological and economic benefits and are becoming more and more important in future wireless systems [1] .
Recently, an EE technique was introduced to cognitive radio (CR) systems for resource allocation, something that has become a hot topic of late. For CR OFDM systems, [2] proposed a search method making use of a water-filling factor to solve the energy-efficiency power allocation (EEPA) problem. In [3] , the EEPA problem is addressed via parametric programming, which is then followed by the presentation of an iterative algorithm. For CR MIMO systems, [4] studied the EE optimization problem of CR MIMO broadcast channels to improve system throughput for unit energy consumption. In [5] , the throughput and EE optimization under quality-of-service constraints for CR MIMO systems are studied. In [6] , a promising framework of spectrum sharing strategy selection, based on EE, is proposed for CR MIMO interference channels.
In this paper, we study the EEPA optimization for CR MIMO-OFDM systems. We try to minimize the energy efficiency of a secondary user (SU) and maintain the minimal rate requirements, without causing excessive interference to primary users (PUs). Since the EEPA problem is non-convex, to find the optimal solution the problem is first transformed into an equivalent convex problem via fractional programming (FP) [7] . Then, an optimal iterative algorithm based on sequential quadratic programming (SQP) [8] is presented.
II. Signal Model and Problem Statement
This paper considers the same CR network as described in [9] . In this CR network, there is a cognitive base station (CBS) and an SU coexisting with L PUs. The set of PUs is denoted by 𝗟 = {1, 2, ... , L }. Here, we consider the downlink scenario, whereby the CBS is equipped with Nt antennas and an SU has Nr antennas. Each PU is equipped with one receive antenna. The CR network adopts OFDM modulation for transmission. The available bands for an SU are divided into N subcarriers and the bandwidth for each subcarrier is Δ f Hz. According to [9] , the rate of the i th subcarrier is defined as
c i = j=1 n min Δf log 2 ( 1+ p ij a ij 2 N 0 ) ,
where n min =min( Nt , Nr ), N 0 is the additive white Gaussian noise (AWGN) variance, p ij is the transmit power at the i th subcarrier at the j th antenna, and aij are the singular values of the Hi ( Hi CNr× Nt is the channel matrix between CBS and an SU in the i th subcarrier). The mutual interference power (MIP) introduced by CBS to the l th PU is defined as
I l = i=1 N j=1 n min p ij | g ij l | 2 φ ij l ,
where the channel fading gains between the j th transmit antenna and the l th PU on the ith subcarrier are
g ij l
. The interference factor of the l th PU on the i th subcarrier is
φ ij l
[10] .
The EE of the CBS is defined as
E E CBS = P c + τ c ( i=1 N j=1 n min p ij ) i=1 N j=1 n min Δf log 2 ( 1+ p ij a ij 2 N 0 ) ,
where τc denotes power-amplifier efficiency and Pc denotes the power consumption of circuits and base-station facilities.
In this paper, our objective is to minimize EE CBS ; therefore, the EEPA problem is formulated as the following OP1
OP1    min p ij μ( { p ij } )  =  E E CBS , s.t.      C1: i=1 N j=1 n min p ij | g ij l | 2 φ ij l I th l ,  l,           C2: i=1 N j=1 n min p ij P max ,           C3: i=1 N j=1 n min Δf log 2 ( 1+ p ij a ij 2 N 0 ) R min ,           C4: p ij 0,  i=1,2,  ...  ,N,  j=1,2,  ...  , n min ,
where µ is the objective function value,
I th l
is the interference threshold of the l th PU, P max is the total power budget of CBS, and R min is the minimal rate requirement.
III. Optimal EEPA
It is difficult to solve OP1 directly, since it is non-convex. Here, we first transform the problem into a convex problem by using FP [7] and then propose an optimal power allocation method based on SQP [8] . We define the elements of the vectors p ∈ ℜ (N×nmin) , ϒl ∈ ℜ (N×nmin) , a ∈ ℜ (N×nmin) , and c ( p ) ∈ ℜ (L+2+N×nmin) , respectively as
p[ ( i−1 )× n min +j ]= p ij ,  i=1,2,  ...  ,N,  j=1,2,  ...  , n min ;
ϒ l [ ( i1 )× n min +j ]= | g ij l | 2 φ ij l ,  i=1,2,  ...  ,N,  j=1,2,  ...  , n min ;
a[ ( i1 )× n min +j ]= a ij 2 ,  i=1,2,  ...  ,N,  j=1,2,  ...  , n min ;
c( p )[m]={ n=1 N× n min p[ n ] ϒ l [ n ] I th l m=l,  l=1,2,  ...  ,L; n=1 N× n min p[ n ] P max   m=L+1; R min n=1 N× n min Δf log 2 ( 1+ a[n]p[n] N 0 ) m=L+2; p [ n ] m=L+2+n , n=1,2,  ...  ,N× n min .
Then, OP1 can be rewritten in vector form as OP2
OP2    min p μ( p )=E E CBS , s.t.       c( p )0.
Before solving OP2, we introduce a new problem OP3 as follows:
OP3   min p f( p,λ ),
f( p,λ )=( P c + τ c n=1 N× n min p[n] )λ n=1 N× n min Δf log 2 ( 1+ a[n]p[n] N 0 ) ,
λ is positive parameter, and denotes the feasible region of OP2.
θ(λ)= min p∈℘ f(p,λ)
p λ  = arg min p∈℘ f(p,λ)
the optimal value and solution of OP3, respectively. Then, OP2 and OP3 can be related to each other by the following lemma [7] .
Lemma 1: The optimal solution, p * , achieves the optimal value λ * = μ ( p * ) of OP2, if and only if, θ ( λ * ) = 0 and pλ* = p * .
By lemma 1, solving OP1 is equivalent to solving problem OP3 for a given λ and then updating λ until Lemma 1 is satisfied. In the following, SQP is used to solve OP3 for a given λ . The Jacobian matrix J ( p ) of c ( p ) is defined as
[ J( p ) ] m,n = c( p )[ m ] p[ n ] ={ ϒ l [ n ]       m=l( l=1,2,  ...  ,L ), 1               m=L+1, Δf ln2 a[ n ] N 0 +a[ n ]p[ n ]   m=L+2, 1​            m=L+2+n, 0​               mL+2+n.
The Lagrangian function is
L(p,β)=( P c + n=1 N× n min p[n] )λ n=1 N× n min Δf log 2 ( 1+ a[n]p[n] N 0 )                    + l=1 L β l ( n=1 N× n min p[n] ϒ l [n] I th l ) + β L+1 ( n=1 N× n min p[n] P max )                    + β L+2 ( R min n=1 N× n min Δf log 2 ( 1+ a[n]p[n] N 0 ) )                    n=1 N× n min β L+2+n p[n] ,
where β = [ β 1 , β 2 , ... , β L+1 , β L+2 , β L+3 , ... , β L+2+N×nmin ] T is the vector of Lagrangian multipliers. Next, we define
Q( p,β )=[ H( p,β )    J T ( p ) J( p )             0 ],
where H ( p , β ) is the Hessian matrix of L ( p , β ) ; the elements of which are defined as
[ H( p,β ) ] mn ={ Δf( λ+ β L+2 ) ln2 ( a[n] ) 2 ( N 0 +a[n]p[n] ) 2     m=n, 0                                                 mn.
The SQP method iteratively improves the estimates of ( pt , βt ) by finding the correction vectors
s t  =  [ ( s p t ) T ( s β t ) T ] T
at the t th iteration to construct the new estimates as follows:
[ p t+1 β t+1 ]=[ p t + η p t s p t β t + η β t s β t ],
η p t
η β t
are non-negative step sizes. The vector st is obtained by solving the following problem:
Q( p t , β t )× s t = [ g ( p t ) T c ( p t ) T ] T ,
where g ( pt ) is the gradient of the Lagrangian function (8); the elements of g ( p ) are
g( p )[ n ]=1 Δf( λ+ β L+2 )a[ n ] ln2( N 0 +a[ n ]p[ n ] ) + l=1 L β l ϒ l [ n ] + β L+1 β L+2+n .
The step sizes
η p t
η β t
also need to be addressed according to the solution point. To measure the progress, a merit function ϕ ( p ) is used. The step size
η p t
is selected so that
ϕ( p t + η p t s p t )
is sufficiently smaller than ϕ ( pt ), while
η β t
is set to be either
η p t
or 1. Here, we use the augmented Lagrangian penalty function to define the merit function
ϕ ω (p,ω)=f(p,λ)+c (p) T χ(p)+ ω 2 c(p) 2 2 ,
‖ • ‖ 2 2
represents the second-order norm, ω is a constant to be chosen, and
χ(p)= [ J (p) T J(p) ] 1 J (p) T ×( f 0 (p,λ)).
Here,∇( f 0 ( p , λ )) represents the gradient of function f 0 ( p , λ ). The multipliers χ ( p ) are the least-squares estimates of the optimal multipliers based on the first-order necessary conditions, that is, KKT conditions. Generally, ω is selected sufficiently large so that the merit function ϕω ( p , ω ) is bounded below. Now, we present the optimal EEPA algorithm based on FP and SQP, which is tabulated as follows:
Algorithm. EEPA. 1: Initialization: initial (p0 ,β0), p 0 * = p 0 , and ϕω(p0) 2: for    τ =0,1, ... , τmax 3:            update λ τ+1 =μ( p τ * ) 4:           fort =0,1, ... ,tmax 5:                     calculate Q( pt ,βt ) and [g(pt) c(pt)T]T 6:                     solve (12) to get st 7:                     select η p t such that ϕ( p t + η p t s p t )<ϕ( p t ) 8:                     update the estimates via (11) 9:              end for 10:            update p τ+1 * = p t max 11: end for
Note: The detailed proofs of the convergence of the inner loop (SQP) and outer loop (FP) can be found in [8] and [7] , respectively.
IV. Performance Simulation
In this section, some numerical results are presented to evaluate the performance of the proposed scheme. Here, we assume the existence of three PUs ( L =3) and that the bandwidths occupied by them are respectively 1 MHz, 2 MHz, and 5 MHz. In addition, there are sixteen OFDM subcarriers ( N = 16) for an SU, and the bandwidth for each subcarrier Δ f is set to 0.3125 MHz. Without loss of generality, the channel gains are assumed to be Rayleigh fading with an average power gain of 1 dB. We then set the following parameters: Nt = Nr = 4, R min = 10 bits/symbol, τc =1, Pc =10 –2 W, and N 0 =10 –13 W. For simplicity, let
I th l = I th .
All the results have been averaged over one thousand iterations.
Figure 1 shows the energy efficiency of EEPA versus parameter ω under different interference threshold Ith and total power budget P max . As shown in Fig. 1 , when ω is large enough, the performance of EEPA will be stable. Since it is impossible to present all cases of the effects of ω on the performance of EEPA, we will, in accordance with the results of Fig. 1 , set ω =10 7 in the following simulation. We compared EEPA with the capacity maximization scheme in [9] . Figure 2 depicts the energy efficiency versus P max , under Ith . The results show that the proposed algorithm can achieve lower energy efficiency than the scheme in [9] . The energy efficiency versus Ith under P max =10 –3 W is evaluated in Fig. 3 . These results depict similar properties to the results in Fig. 2 .
PPT Slide
Lager Image
Energy efficiency vs. parameter ω.
PPT Slide
Lager Image
Energy efficiency vs. total power budget.
PPT Slide
Lager Image
Energy efficiency vs. interference threshold.
V. Conclusion
In this letter, we investigated the EEPA problem in CR MIMO-OFDM systems. Since the objective function is nonconvex, we first construct an equivalent convex problem via FP, and then a new iterative power allocation algorithm is proposed based on SQP. From the simulation results, we observed that the proposed method can improve energy efficiency compared with the conventional capacity maximization scheme.
This work was supported by Natural Science Foundation of China (Grant No. 60872073, No. 51075068, No. 60975017, and No. 61301219) and Doctoral Fund of Ministry of Education of China (No. 20110092130004).
Feng D. 2013 “A Survey of Energy-Efficient Wireless Communications,” IEEE Commun. Surveys Tutorials 15 (1) 167 - 178    DOI : 10.1109/SURV.2012.020212.00049
Mao J.L. 2013 “Energy Efficiency Optimization for OFDM-Based Cognitive Radio Systems: A Water-Filling Factor Aided Search Method,” IEEE Trans. Wireless Commun. 12 (5) 2366 - 2375    DOI : 10.1109/TWC.2013.013013.121013
Yan W. 2012 “Optimal Energy-Efficient Power Allocation for OFDM-Based Cognitive Radio Networks,” IEEE Commun. Lett. 16 (9) 1420 - 1423    DOI : 10.1109/LCOMM.2012.070512.120662
Mao J.L. 2013 “Energy Efficiency Optimization for Cognitive Radio MIMO Broadcast Channels,” IEEE Commun. Lett. 17 (2) 337 - 340    DOI : 10.1109/LCOMM.2012.121912.122351
Akin S. , Gursoy M.C. 2013 “On the Throughput and Energy Efficiency of Cognitive MIMO Transmissions,” IEEE Trans. Veh. Technol. 62 (7) 3245 - 3260    DOI : 10.1109/TVT.2013.2254142
Zhong W. , Wang J. “Energy Efficient Spectrum Sharing Strategy Selection for Cognitive MIMO Interference Channels,” IEEE Trans. Signal Process. 61 (14) 3705 - 3717    DOI : 10.1109/TSP.2013.2259162
Dinkelbach W. 1967 “On Nonlinear Fractional Programming,” Manag. Sci. 13 (7) 492 - 498    DOI : 10.1287/mnsc.13.7.492
Boggs P.T. , Tolle J.W. 1996 “Sequential Quadratic Programming,”
Shahraki H.S. , Mohamed-pour K. 2011 “Power Allocation in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing-Based Cognitive Radio Networks,” IET Commun. 5 (3) 362 - 370    DOI : 10.1049/iet-com.2010.0262
Bansal G. , Hossain M.J. , Bhargava V.K. 2011 “Adaptive Power Loading for OFDM-Based Cognitive Radio Systems with Statistical Interference Constraint,” IEEE Trans. Wireless Commun. 10 (0) 2786 - 2791    DOI : 10.1109/TWC.2011.072011.100397