Energy-Efficient Dual-Iteration Power Allocation for Two-Phase Relay System with Massive Antennas

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 themultiple 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.


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, realtime 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.
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.

System Model and Spectrum Efficiency
A two-phase AF relay system model is illustrated in Figure 1.Multiple single-antenna source nodes S = {S 1 , S 2 , . . ., S  } communicate with multiple single-antenna destination nodes D = {D 1 , D 2 , . . ., D  } assisted by the relay node (R) equipped with  antennas, and S  → D  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  0 , maximum radius  0 , 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 where x S = [ S 1 ,  S 2 , . . .,  S  ] T is the modulated transmission symbol vector and E[x S x H S ] = I  .p S = [ S 1 ,  S 2 , . . .,  S  ] T is the power allocation vector at S. The channel matrix from S to R is SR .The large-scale fading matrix D SR = diag{ S 1 R ,  S 2 R , . . .,  S  R } and  S  R =  − S  R  S  R , where  is a constant related to the carrier frequency and antenna gain,  S  R is the distance between S  and R,  ∈ [2,6] is the path loss exponent, and  S  R is a log-normal shadow fading variable with distribution 10log 10  S  R ∼ CN(0,  2 sh,S  R ).The elements of the fast fading matrix H SR = [ℎ , ] ∈ C × are independent and identically distributed (i.i.d.) random variables with the distribution of CN(0, The signals after the backward filter can be written as In phase II, R transmits signal x R to D, and the received signals at D can be expressed as where T is the power allocation vector at R for the multiple pairs of nodes.The channel  from R to D is where  and  are identical as in phase I,  RD  is the distance between R and D  , and  RD  is a lognormal shadow fading variable with the distribution of 10log 10  RD  ∼ CN(0,  2 sh,RD  ).The fast fading matrix −1 is the forward ZF filter matrix at R, and The received signal-to-noise ratio (SNR) at D  can be readily shown as where   and   are the channel-to-noise ratios (CNRs) of the links S  → R and R → D  , respectively, and can be written as: Furthermore, denoting the SE vector of multiple node pairs by c = [ 1 ,  2 , . . .,   ] T , the overall SE of the entire relay system can be expressed as According to random matrix theory, both of the random variables ‖b  ‖ −2 and ‖f  ‖ −2 follow a gamma distribution with the parameters ( −  + 1, 1) [11,12].Then, the expectations of random variables ‖b  ‖ 2 and ‖f  ‖ 2 can be expressed as ) holds.Therefore, the lower bound of the SE of the relay system can be given as where the equivalent SNRs and CNRs can be written as follows, respectively: 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.

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: where  S  and  R represent the peak to average power ratios (PAPRs) of OFDM at the source nodes and relay node, respectively. S  and  R denote the efficiencies of the power amplifiers (PAs) at the source nodes and relay node, respectively. is the system bandwidth. SC ( = 1, 2, 3),  RC ( = 1, 2, 3, 4), and  DC ( = 1, 2, 3) are the constant coefficients of circuit power consumption at the source, relay, and destination nodes, respectively.
International Journal of Antennas and Propagation Therefore, the total power consumption per Hz of the whole relay system can be written as where  C1 =  SC1 +  DC1 and Based on the metric J/bit, the optimal EE problem of the overall relay system can be formulated as where  S  ,max and  R,max are the maximum transmission powers at S  and R, respectively.

Energy-Efficient Power Allocation
4.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 p S at the source nodes, the optimal power allocation vector p o R at the relay will be first derived.Then, taking advantage of the symmetry of the received SNR at D, the optimal power allocation vector p o S at the source nodes can be obtained if p R 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 p o R .
Theorem 1.When the power allocation vector p  is given at the source nodes, the SE of the th node pair achieving the optimal EE  o  is given by where Proof.When p S is given in the feasible set,  R  can be given by plugging ( 4) into (6) as follows: Substituting ( 18) into ( 13),  Σ becomes a function of c.The first-order partial derivative of  Σ (c) can be shown as Moreover, the second-order partial derivative of  Σ (c) can be written as Thus,  Σ (c) is a quasiconvex function with respect to c, since the Hessian matrix of  Σ (c) is positive semidefinite; namely, and the objective function of ( 15) is quasiconvex.
Letting c) is a convex function for a given  b ≥ 0 and decreases with  b .According to [15], when  b =  o b , the objective function of the quasiconvex problem ( 15) can be transformed into which can be solved by the following stationary condition owing to its convexity: Through simplification, (22) can be rewritten as Finally, solving the quadratic equation and saving the nonnegative root give rise to 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 According to the definition of function ( b , c), the following relationships hold [16]: Then, given p S at S and considering the constraint T R , the optimal p o R at R can be achieved by the bisection algorithm, which will be described in the subalgorithm of Algorithm 1.
Similarly, given the power allocation p R at R and taking advantage of the symmetry of the received SNR at D  , the optimal SE of the th node pair can be shown as where Furthermore, considering the constraint T S , one can obtain . Therefore, the optimal p o S achieving the optimal EE at S is also attainable by the bisection algorithm when p R 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 p X at X(= S or R), ignoring the transmission power p Y at Y(= R or S) and replacing the SE c with its upper bound c U (p X ) = [ U 1 ,  U 2 , . . .,  U  ] T , a lower-bounded EE  L b can be obtained as follows: where and a SE upper bound  U  for the th node pair is where   =   if X = S; else   =   .

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 largescale fading multiplied by the number of antennas [5][6][7].
Then, we can replace the instantaneous (  ,   ) in ( 5) with the statistic (  ,   ) 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].

Simulation and Discussions
This section evaluates the EE and SE of the proposed energyefficient 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 , and  S  ,max =  S,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 ∘ ,  0 = 35m, and different  0 , as shown in Figure 1. Figure 3 plots the average EE versus the number of antennas at the relay node with various radii  0 , when the number of node pairs  is 10. Figure 4 shows the average EE versus the number of node pairs with different numbers of antennas  at the relay node when the radius  0 is 0.5 km.Both Figures 3 and 4   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.

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.

Figure 1 :
Figure 1: Relay system serving multiple pairs of nodes with massive antennas.

Figure 2 :
Figure 2: Components of power consumption at the relay node.

Figure 4 :
Figure 4: EE versus the number of node pairs .

Figure 5 :
Figure 5: SE versus the number of antennas  at the relay.

Figure 6 :
Figure 6: SE versus the number of node pairs .
o b ,  S  )] + for the th node pair at R. Since both  R  ( o  ) and  o  ( o b ) are monotonically increasing functions, the composite function  o R  ( o b) is also a monotonically increasing one.Therefore, there exists an upper limit  U b of the EE when considering the constraint 1 T p R ≤  R,max .Step 1. (Initialization):

Table 1 :
Simulation parameters for the relay system.EE versus the number of antennas  at the relay.
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