This paper considers the scenario where multiple source nodes communicate with multiple destination nodes simultaneously with the aid of an amplify-and-forward relay equipped with massive antennas. In order to achieve optimal energy efficiency (EE) of the entire relay system, this paper investigates the power allocation problem for the multiple pairs of nodes at both the source nodes and the relay node, where the relay employs the backward and forward zero-forcing filters. Since the EE optimization problem cannot be solved analytically, we propose a two-phase power allocation method. Given power allocation of one phase, the optimal power allocation is derived for the other phase. Furthermore, two dual-iteration power allocation (DIPA) algorithms with performance approaching that of optimal EE are developed based on the instantaneous and statistic channel state information, respectively. Numerical results show that the proposed DIPA algorithms can greatly improve EE while guaranteeing spectrum efficiency (SE) when compared with the equal power allocation algorithm. Moreover, both algorithms suggest that deploying a rational number of antennas at the relay node and multiplexing a reasonable number of node pairs can improve on the EE and SE.
1. Introduction
With the rapid development of smart terminals and their new applications, the traffic of wireless communication networks grows exponentially [1]. Moreover, power consumption caused by the networks results in excessive carbon emissions and unaffordable operator expenditure [2]. Therefore, both spectrum efficiency (SE) and energy efficiency (EE) are considered as significant metrics for evaluating the performance of communication systems. In recent years, massive multiple-input multiple-output (MIMO) has been put forward to simultaneously serve tens of users when the base station (BS) employs hundreds of antennas [3]. Massive MIMO can not only improve SE, but also enhance system EE [3, 4]. Apart from these advantages, massive MIMO is able to enhance the link reliability due to its higher diversity gains [5], reduce interuser interference, and guarantee physical layer security owing to the extremely narrow beam [5, 6]. In a massive MIMO system, individual antenna failure of the antenna array will not be catastrophic to the system performance [5, 6]. The scheduling scheme can be made simply due to the channel harden phenomenon [7]. Certainly, since the maximum number of pilot sequences is bounded by the coherent time and bandwidth, pilot contamination due to pilot reuse in multicell scenarios imposes the ultimate limit on the performance of massive MIMO system [3, 5, 6].
Meanwhile, maximum ratio transmission (MRT)/zero forcing (ZF) and maximum ratio combination (MRC)/ZF techniques have been proven as practical precoders and detectors [5, 6]. The performance, including SE, EE, and reliability, of massive MIMO with different precoders/detectors has been investigated broadly for single cell or multicell [4, 8, 9]. However, relays with massive antennas serving multiple pairs of nodes are less studied. Such a relay system is a type of isolated system, which exchanges less necessary or no information with the networks. Taking the stricken area as an example, a temporary relay can be deployed, which is mainly used for mutual communication by rescue workers or local residents. As depicted in Figure 1, the ergodic SE of multiple pairs of nodes served by an amplify-and-forward (AF) relay with the ZF/MRT filter was analyzed when the number of antennas at the relay goes to infinity [10]. However, real-time SE with a finite number of antennas is more significant for a practical relay system. Moreover, realistic EE should be evaluated in consideration of both radiated and circuit powers.
Relay system serving multiple pairs of nodes with massive antennas.
This paper investigates real-time power allocation at both source nodes and relay node with a finite number of antennas in order to achieve optimal EE for the entire relay system. When the two-phase AF relay employs the backward and forward ZF filters, both the real-time SE and the lower bound of ergodic SE are first analyzed. Then, the EE optimization problem is formulated based on a realistic power consumption model, which consists of both radiated and circuit powers. Due to the intractable objective function, the optimization problem is decomposed into a two-phase power allocation problem. The optimal single-phase power allocation is derived when the power allocation of the other phase is given using convex optimization theory. Furthermore, considering both the real-time SE and ergodic SE, two dual-iteration power allocation (DIPA) algorithms capable of achieving near optimal EE performance are developed for the entire relay system, which can simultaneously improve EE and SE compared with the equal power allocation (EPA) algorithm.
The remainder of this paper is organized as follows. Section 2 introduces the system model and analyzes the system SE. The EE optimization problem is formulated in consideration of both radiated and circuit powers in Section 3. The optimal single-phase power allocation is derived and two DIPA algorithms for the entire relay system are developed in Section 4. Section 5 presents simulation results and Section 6 presents the conclusion of this paper.
Notation. Uppercase boldface letters and lowercase boldface letters denote matrices and vectors, respectively. (·)T, (·)H, and (·)-1 represent the transpose, conjugate transpose, and pseudoinverse of matrix/vector, respectively. diag{a} is a diagonal square matrix whose main diagonal is formed by vector a. ∥·∥ represents the Euclidian 2-norm. Consider [x]+=max{x,0} and [x]ab=min{max{a,x},b}. ℰ[·] represents mathematical expectation. 𝒩(μ,σ2) denotes the real Gaussian distribution with mean μ and variance σ2. 𝒞𝒩(μ,σ2) is a complex Gaussian distribution with mean μ and real/imaginary component variance σ2/2.
2. System Model and Spectrum Efficiency
A two-phase AF relay system model is illustrated in Figure 1. Multiple single-antenna source nodes S={S1,S2,…,SK} communicate with multiple single-antenna destination nodes D={D1,D2,…,DK} assisted by the relay node (R) equipped with N antennas, and Sk→Dk is a pair of nodes. In order to simplify the simulation in Section 5, the source and destination nodes are distributed in an area with minimum radius r0, maximum radius R0, and angle θ. But those assumptions have no effect on the theoretical deduction. This paper considers block Rayleigh fading channels and assumes that the perfect channel state information (CSI) can be acquired at R. We also assume that the timely power allocation information can be transmitted to nodes S or D free of errors after R implementing the power allocation algorithm.
In phase I, S transmits signals to R, and the received signals can be expressed as
(1)yR=GSR(diag{pS})1/2xS+nR,
where xS=[xS1,xS2,…,xSK]T is the modulated transmission symbol vector and E[xSxSH]=IK. pS=[pS1,pS2,…,pSK]T is the power allocation vector at S. The channel matrix from S to R is GSR=HSRDSR1/2. The large-scale fading matrix DSR=diag{βS1R,βS2R,…,βSKR} and βSkR=ϕdSkR-αξSkR, where ϕ is a constant related to the carrier frequency and antenna gain, dSkR is the distance between Sk and R, α∈[2,6] is the path loss exponent, and ξSkR is a log-normal shadow fading variable with distribution 10log10ξSkR~𝒞𝒩(0,σsh,SkR2). The elements of the fast fading matrix HSR=[hn,k]∈𝒞N×K are independent and identically distributed (i.i.d.) random variables with the distribution of 𝒞𝒩(0,1). nR~𝒞𝒩(0,σR2IN×N) is the additive white Gaussian noise (AWGN) at the relay node R.
At R, a backward ZF filter B=[b1T,b2T,…,bKT]T=(HSRHHSR)-1HSRH is employed, and its normalization matrix is QB=(diag{qB1,qB2,…,qBK})-1/2, where qBk=βSkRpSk+∥bk∥2σR2 (k=1,2,…,K). The signals after the backward filter can be written as
(2)xR=QBByR.
In phase II, R transmits signal xR to D, and the received signals at D can be expressed as
(3)yD=GRDFQF(diag{pR})1/2xR+nD,
where pR=[pR1,pR2,…,pRK]T is the power allocation vector at R for the multiple pairs of nodes. The channel from R to D is GRD=DRD1/2HRD. The large-scale fading matrix DRD=diag{βRD1,βRD2,…,βRDK} and βRDk=ϕdRDk-αξRDk, where ϕ and α are identical as in phase I, dRDk is the distance between R and Dk, and ξRDk is a log-normal shadow fading variable with the distribution of 10log10ξRDk~𝒞𝒩(0,σsh,RDk2). The fast fading matrix HRD=[hkn]∈𝒞K×N contains i.i.d. 𝒞𝒩(0,1) entries. Consider F=[f1,f2,…,fK]=HRDH(HRDHRDH)-1 is the forward ZF filter matrix at R, and QF=(diag{qF1,qF2,…,qFK})-1/2 is its normalization matrix with qFk=∥fk∥2. nD is the AWGN at D with the distribution of 𝒞𝒩(0,diag{σD12,σD22,…,σDK2}IK).
The received signal-to-noise ratio (SNR) at Dk can be readily shown as
(4)ρk=akpSkbkpRkakpSk+bkpRk+1,
where ak and bk are the channel-to-noise ratios (CNRs) of the links Sk→R and R→Dk, respectively, and can be written as:
(5)[ak,bk]=[βSkRσR2∥bk∥-2,βRDkσDk2∥fk∥-2].
Furthermore, denoting the SE vector of multiple node pairs by c=[c1,c2,…,cK]T, the overall SE of the entire relay system can be expressed as
(6)CΣ=∑k=1Kck=12∑k=1Klog2(1+ρk).
According to random matrix theory, both of the random variables ∥bk∥-2 and ∥fk∥-2 follow a gamma distribution with the parameters (N-K+1,1) [11, 12]. Then, the expectations of random variables ∥bk∥2 and ∥fk∥2 can be expressed as ℰ[∥bk∥2]=ℰ[∥fk∥2]=(N-K)-1. Owing to the convexity of function log2(1+1/x), the Jensen inequality ℰ[log2(1+1/x)]≥log2(1+1/ℰ[x]) holds. Therefore, the lower bound of the SE of the relay system can be given as
(7)C¯Σ=ℰ[CΣ]≥∑k=1Kc¯k=12∑k=1Klog2(1+ρ¯k),
where the equivalent SNRs and CNRs can be written as follows, respectively:
(8)ρ¯k=ℰ-1[ρk-1]=a¯kpSkb¯kpRka¯kpSk+b¯kpRk+1,(9)[a¯k,b¯k]=[(N-K)βSkRσR2,(N-K)βRDkσDk2].
The lower bound of ergodic SE in (7) is only dependent on the statistic CSI compared with the real-time SE in (6), which relates to the instantaneous CSI.
3. Energy Efficiency Formulation
According to [13], the power consumption model consists of both the radiated and circuit powers. The components of the source and destination nodes are similar to those in [13, 14]. The components of power consumption for the relay node are depicted in Figure 2. Then, the power consumption models of the source, relay, and destination nodes are given as follows, respectively:
(10)PSΣ=∑k=1K(ςSkηSkpSkB︸PA+pSC1ckB+︸coder&modulatorpSC2B︸OFDM+pSC3︸other),(11)PRΣ=ςRηR∑k=1KpRkB︸PA+pRC2NB︸OFDM+pRC3N︸filter,etc+pRC4︸constant+pRC1(7NK2+4NK-2K+7K3)B︸ZFdetector&precoder,(12)PDΣ=∑k=1K(pDC1ckB+︸demodulator&decoderpDC2B︸OFDM+pDC3︸other),
where ςSk and ςR represent the peak to average power ratios (PAPRs) of OFDM at the source nodes and relay node, respectively. ηSk and ηR denote the efficiencies of the power amplifiers (PAs) at the source nodes and relay node, respectively. B is the system bandwidth. pSCi (i=1,2,3), pRCi (i=1,2,3,4), and pDCi (i=1,2,3) are the constant coefficients of circuit power consumption at the source, relay, and destination nodes, respectively.
Components of power consumption at the relay node.
Therefore, the total power consumption per Hz of the whole relay system can be written as
(13)PΣ=1B(PSΣ+PRΣ+PDΣ)=∑k=1K(ςSkηSkpSk+ςRηRpRk+pC1ck)+g(N,K,B),
where pC1=pSC1+pDC1 and
(14)g(N,K,B)=pRC1(7NK2+4NK-2K+7K3)+pRC2N+(pSC2+pDC2)K+1B[pRC3N+(pSC3+pDC3)K+pRC4].
Based on the metric J/bit, the optimal EE problem of the overall relay system can be formulated as
(15)Ebo=minEb=min{PΣCΣ}s.t.{TS:0≤pSk≤PSk,max,k=1,2,…,K,TR:0≤∑k=1KpRk≤PR,max,
where PSk,max and PR,max are the maximum transmission powers at Sk and R, respectively.
4. Energy-Efficient Power Allocation4.1. DIPA Based on Instantaneous CSI
It is easy to see that (15) is a computationally intractable problem according to convex optimization theory. In this part, given the power allocation vector pS at the source nodes, the optimal power allocation vector pRo at the relay will be first derived. Then, taking advantage of the symmetry of the received SNR at D, the optimal power allocation vector pSo at the source nodes can be obtained if pR is given. Furthermore, we will propose an energy-efficient DIPA algorithm which is capable of achieving near-optimal EE performance for the entire relay system.
We present the following theorem to obtain pRo.
Theorem 1.
When the power allocation vector pS is given at the source nodes, the SE of the kth node pair achieving the optimal EE Ebo is given by
(16)cko={log2[(δRk+2-δRk)akpSk+12]}+,
where
(17)δRk=ςRakpSkln2ηRbk(Ebo-pC1).
Proof.
When pS is given in the feasible set, pRk can be given by plugging (4) into (6) as follows:
(18)pRk(ck)=(akpSk+1)/bk((akpSk/(22ck-1))-1).
Substituting (18) into (13), PΣ becomes a function of c. The first-order partial derivative of PΣ(c) can be shown as
(19)∂PΣ∂ck=akςR22ck+1bkηRlog2eak+pSk-1[ak-(22ck-1)pSk-1]2+pC1.
Moreover, the second-order partial derivative of PΣ(c) can be written as
(20)∂2PΣ∂ck2=2log2e[22ck+1pSk-1ak-(22ck-1)pSk-1+1]∂PΣ∂ck.
Thus, PΣ(c) is a quasiconvex function with respect to c, since the Hessian matrix of PΣ(c) is positive semidefinite; namely, ∇2PΣ(c)=diag{∂2PΣ/∂c12,∂2PΣ/∂c22,…,∂2PΣ/∂cK2}≥0. Besides, CΣ(c) is an affine function of c. Therefore, the lower-level set {c∣Eb(c)=PΣ(c)/CΣ(c)≤eb} is a convex set, and the objective function of (15) is quasiconvex.
Letting F(eb,c)=PΣ(c)-ebCΣ(c), F(eb,c) is a convex function for a given eb≥0 and decreases with eb. According to [15], when eb=Ebo, the objective function of the quasiconvex problem (15) can be transformed into
(21)Fo(Ebo,co)=minc≥0{PΣ(c)-EboCΣ(c)},
which can be solved by the following stationary condition owing to its convexity:
(22)∂F(Ebo,c)∂ck=∂PΣ(c)∂ck-Ebo=0.
Through simplification, (22) can be rewritten as
(23)22ck+2(akpSk+1)δRk2ck-(akpSk+1)=0.
Finally, solving the quadratic equation and saving the nonnegative root give rise to
(24)2ck=(δRk+2-δRk)akpSk+12.
By applying the logarithmic operation to both sides of (24) and considering the nonnegativity of SE, (16) in the theorem is thus proved.
Substituting (16) into (18) leads to the optimal power allocation pRko=[pRk(cko)]+=[pRk(Ebo,pSk)]+ for the kth node pair at R. Since both pRk(cko) and cko(Ebo) are monotonically increasing functions, the composite function pRko(Ebo) is also a monotonically increasing one. Therefore, there exists an upper limit EbU of the EE when considering the constraint 1TpR≤PR,max.
According to the definition of function F(eb,c), the following relationships hold [16]:
(25)F(eb,c){>0,eb<Ebo,=0,eb=Ebo,<0,eb>Ebo.
Then, given pS at S and considering the constraint TR, the optimal pRo at R can be achieved by the bisection algorithm, which will be described in the subalgorithm of Algorithm 1.
Algorithm 1: Energy-efficient dual-iteration power allocation.
Step 1. (Initialization): pSk=PSk,max(k=1,2,…,K), EbU=∞, EbL=EbL(pS) in (28), and ε>0.
Step 2. (Update for relay node): X=S and Y=R.
Sub-step 1: Executing the Sub-Algorithm.
Sub-step 2: pRko=pRk(pSk,cko), EbU=Ebo(pS) and EbL=EbL(pRo).
Step 3. (Update for source nodes): X=R and Y=S.
Sub-step 1: Executing the Sub-Algorithm.
Sub-step 2: pSko=pSk(pRk,cko), EbU=Ebo(pR) and EbL=Eb(pSo).
Step 4. (Judgment): if Ebo(pS)-Ebo(pR)≤ε, stop; else go back to Step 2.
Sub-Algorithm: Bi-section algorithm
Step 1. Let u=1/EbU and v=1/EbL.
Step 2. Compute κ=(u+v)/2, eb=1/κ and cko(eb,pX).
Step 3. Compute f=F(eb,co(eb,pX)). If f>0 and the constraint TY is valid, then u=κ; Else v=κ.
Step 4. If v-u<ε, Ebo(pX)=ebo; Else go back to Step 2.
Similarly, given the power allocation pR at R and taking advantage of the symmetry of the received SNR at Dk, the optimal SE of the kth node pair can be shown as
(26)cko={log2[(δSk+2-δSk)bkpRk+12]}+,
where
(27)δSk=bkςSkpRkln2akηSk(Ebo-pC1).
Furthermore, considering the constraint TS, one can obtain pSko=[pSk(cko)]0PSk,max=[pSk(Ebo,pRk)]0PSk,max. Therefore, the optimal pSo achieving the optimal EE at S is also attainable by the bisection algorithm when pR at R is given.
Considering the power allocation for both phases, we develop an energy-efficient DIPA algorithm for multiple node pairs at both S and R as detailed in Algorithm 1.
In the DIPA algorithm, given the transmission power pX at X(=SorR), ignoring the transmission power pY at Y(=RorS) and replacing the SE c with its upper bound cU(pX)=[c1U,c2U,…,cKU]T, a lower-bounded EE EbL can be obtained as follows:
(28)EbL(pX)=Pother(pX)1TcU(pX),
where
(29)Pother(pX)=PΣ-ςYηY1TpY
and a SE upper bound ckU for the kth node pair is
(30)ckU(pXk)=12log2(1+pXkγk),
where γk=ak if X=S; else γk=bk.
4.2. DIPA Based on Statistic CSI
As Algorithm 1 is dependent on the fast fading channels, namely, the instantaneous CSI, the different symbols on different subcarriers from the same source node should be calculated separately. Therefore, the complexity is overwhelming for practical systems. Thanks to the channel harden phenomenon of massive MIMO, the norm of each channel link can be approximated by the large-scale fading multiplied by the number of antennas [5–7]. Then, we can replace the instantaneous (ak,bk) in (5) with the statistic (a¯k,b¯k) in (9) in the DIAP algorithm. In other words, the instantaneous SE expression in (6) can be replaced with the ergodic lower bound in (7) during the derivations of Theorem 1. Thus, the proposed algorithm can be implemented efficiently owing to nearly the same large-scale fading for different subcarriers of each link. Moreover, the method of DIPA based on both the statistic CSI and instantaneous CSI can be easily extended to other precoders and detectors. And the stochastic analysis of relay with massive antennas under Markov channel is our future interest [17].
5. Simulation and Discussions
This section evaluates the EE and SE of the proposed energy-efficient DIPA algorithms compared with the EPA algorithm. The results of DIPA based on the instantaneous CSI and statistic CSI are denoted by “IC” and “SC”, respectively, in Figures 3–6. For the purpose of fair comparison, the total radiated power of the source nodes or relay node used in the EPA algorithm is made the same as the actual power consumed by the proposed DIPA algorithm. In order to simplify the simulation, without loss of generality, it is assumed that ηSk=ηR=η, ςSk=ςR=ς, σsh,SkR2=σsh,RDk2=σsh2, σR2=σDk2=σ2, and PSk,max=PS,max. Similar to [13, 14], the main parameters used in our simulations are listed in Table 1. During each instance of simulation, the source nodes or destination nodes are uniformly distributed in a given area with θ=60°, r0=35m, and different R0, as shown in Figure 1.
Simulation parameters for the relay system.
Parameters
Values
ϕ
1
η
0.35
ς
1.2705
α
3.7
σsh
8 dB
r0
35 m
σ2
−124 dBm/Hz
ε
1E-8
Δf
15 kHz
pDC1=2pSC1
0.2 mW/Δf
pRC1
0.5 μW/Δf
pRC2=2pSC2=2pDC2
0.8 mW/Δf
pRC3=2pSC3=2pDC3
2 W
pRC4
10 W
PS,max
1 μW/Hz
PR,max
1 mW/Hz
EE versus the number of antennas N at the relay.
EE versus the number of node pairs K.
SE versus the number of antennas N at the relay.
SE versus the number of node pairs K.
Figure 3 plots the average EE versus the number of antennas at the relay node with various radii R0, when the number of node pairs K is 10. Figure 4 shows the average EE versus the number of node pairs with different numbers of antennas N at the relay node when the radius R0 is 0.5 km. Both Figures 3 and 4 indicate that the proposed DIPA algorithms can improve the EE compared with the EPA algorithm, which validates the proposed DIPA algorithms. The results of DIPA based on the statistic CSI agree with the results of DIPA based on the instantaneous CSI. Thus, the DIPA based on the statistic CSI is more advisable for practical communication systems. Moreover, the EE first decreases and then increases with respect to both the number of antennas of the relay node and the number of multiple node pairs, since the rate of SE increase is first faster and then slower than the rate of power increase. Therefore, deploying a reasonable number of antennas at the relay node and multiplexing a rational number of node pairs can enhance the system EE. Furthermore, Figure 3 indicates that shrinking the radius of the serving area can enhance the EE, since less power is consumed to compensate the path loss.
The corresponding SEs to both Figures 3 and 4 are given in Figures 5 and 6, respectively. Because the algorithms consume the same amount of power, the DIPA algorithms perform better than the EPA algorithm in terms of SE according to the definition of EE. The DIPA based on the statistic CSI performs the same as the DIPA based on the instantaneous CSI. Both figures demonstrate that SE can be improved by deploying more antennas at the relay node. Multiplexing a rational number of node pairs and shrinking the radius of the severing area can improve system SE, as can be observed in Figures 6 and 5, respectively.
6. Conclusion
This paper studied energy-efficient power allocation at both source nodes and relay node for a two-phase AF relay system, where the relay is equipped with massive antennas and employs backward and forward ZF filters. Through decomposing the formulated EE optimization problem into two phases, two DIPA algorithms are proposed based on the instantaneous CSI and statistic CSI, respectively. Simulation results demonstrate the effectiveness of the proposed algorithms compared with the EPA algorithm. The DIPA algorithm based on the statistic CSI is more advisable compared with that based on the instantaneous CSI for practical systems in consideration of implementation complexity. Moreover, deploying a rational number of antennas at the relay and multiplexing a reasonable number of node pairs can improve the system performance. Shrinking the radius of the severing area is also an important measure for better SE and EE.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National 863 Program of China (2014AA01A705), National Key Technology R&D Program of China (2012ZX03004005), China Natural Science Funding (61271183), and Research Fund for the Doctoral Program of Higher Education (20130005120003).
HasanZ.BoostanimehrH.BhargavaV. K.Green cellular networks: a survey, some research issues and challenges20111345245402-s2.0-8485785746810.1109/SURV.2011.092311.00031FehskeA.FettweisG.MalmodinJ.BiczokG.The global footprint of mobile communications: the ecological and economic perspective201149855622-s2.0-8005170865610.1109/MCOM.2011.5978416MarzettaT. L.Noncooperative cellular wireless with unlimited numbers of base station antennas2010911359036002-s2.0-7844930160010.1109/TWC.2010.092810.091092NgoH. Q.LarssonE. G.MarzettaT. L.Energy and spectral efficiency of very large multiuser MIMO systems201361414361449RusekF.PerssonD.BuonK. L.LarssonE. G.MarzettaT. L.EdforsO.TufvessonF.Scaling up MIMO: opportunities and channenges with very large arrays20133014060LarssonE. G.TufvessonF.EdforsO.MarzettaT. L.Massive MIMO for next generation wireless systems2014522186195HochwaldB. M.MarzettaT. L.TarokhV.Multiple-antenna channel hardening and its implications for rate feedback and scheduling2004509189319092-s2.0-454423459310.1109/TIT.2004.833345HoydisJ.ten BrinkS.DebbahM.Massive MIMO in the UL/DL of cellular networks: how many antennas do we need?2013312160171YangH.MarzettaT. L.Performance of conjugate and zeroforcing beamforming in large-scale antenna system2013312172179HimalA. S.HienQ. N.TrungQ. D.ChauY.ErikG. L.Multipair amplify-and-forward relaying with very large antenna arraysProceedings of the IEEE International Conference on Communications (ICC ’13)June 2013Budapest, Hungary32283233CouilletR.DebbahM.2011New York, NY, USACambridge University PressZhaoL.ZhengK.LongH.ZhaoH.Performance analysis for downlink massive MIMO system with ZF precoding201310.1002/ett.2745KimH. S.DaneshradB.Energy-constrained link adaptation for MIMO OFDM wireless communication systems201099282028322-s2.0-7795692094610.1109/TWC.2010.062910.090983XuZ.YangC.LiG. Y.ZhangS.ChenY.XuS.Energyefficient configuration of spatial and frequency resources in MIMOOFDM sytems2013612564575IshedenC.ChongZ.JorswieckE.FettweisG.Framework for link-level energy efficiency optimization with informed transmitter201211829462957BoydS.VandenbergheL.2004New York, NY, USACambridge University PressZhengK.LiuF.LeiL.LinC.JiangY.Stochastic performance analysis of a wireless finite-state markov channel2013122782793