Abstract
This paper focuses on radio resource allocation in OFDMA networks for maximizing the energy efficiency subject to the data rate requirements of users. We propose the energyefficient waterfilling structure to obtain the closedform optimal energyefficient power allocation for a given subcarrier assignment. Moreover, we establish a new sufficient condition for the optimal energyefficient subcarrier assignment. Based on the theoretical analysis, we develop a joint energyefficient resource allocation (JERA) algorithm to maximize the energy efficiency. Simulation results show that the JERA algorithm can yield optimal solution with significantly low computational complexity.
1. Introduction
As the high data rate applications are going to dominate the mobile services, energy efficiency (EE) is becoming more crucial in wireless communication networks. Energyefficient radio resource allocation is one of the effective ways to improve the EE of the OFDMA (orthogonal frequency division multiple access) networks [1]. Although radio resource allocation in OFDMA networks has been extensively studied, the major focus is on improving spectral efficiency (SE) which may not always coincide with EE [2].
Different from SEbased resource allocation schemes in which the total transmitting power is fixed, the EEbased schemes adjust the power level adaptively based on the channel conditions [3]. Accordingly, the classic waterfilling power allocation method cannot be applied directly due to the unknown total transmitting power. In order to determine the proper transmitting power, the perturbation functions of EE have been studied in [4, 5]. It has been shown that EE is strictly quasiconcave in SE [4] and in the total transmitting power [5]. However, the perturbation functions of EE cannot be expressed in an analytic expression; only the approximation algorithms have been proposed to find the nearoptimal solutions [4–6].
In this paper, we focus on joint subcarrier assignment and power allocation in OFDMA networks for maximizing the energy efficiency problem subject to the data rate requirements of users. The main contributions of our work are summarized as follows.(i)We prove that EE is strictly pseudoconcave with respect to power vector for a given subcarrier allocation, which guarantees that the solution satisfying the KKT conditions is also the global optimal. Using this property, we show that the optimal solution has a special EE waterfilling structure that is determined by only one variable. Based on this observation, we further provide the first closedform expression for the optimal energyefficient power allocation.(ii)According to the analysis, we propose an optimal energyefficient power allocation algorithm by sequentially searching within a finite number of waterlevel intervals. The computational complexity of the proposed algorithm is much lower than that of the algorithms approaching the optimal solution with iterative searching method.(iii)We provide a sufficient condition for the optimal energyefficient subcarrier assignment. This condition is the basis of the quick search method, because it allows us to easily determine whether a subcarrier assignment is optimal. As a result, by combining the sufficient condition and the EE waterfilling solution, we design a joint energyefficient resource allocation (JERA) algorithm, which not only achieves the optimal EE, but also outperforms the existing algorithms in terms of convergence rate.
The state of the art to solve the EE maximization problem is JIOO (joint inner and outerlayer optimization) algorithm [4, 5]. Relying on the strict quasiconcavity of the perturbation function of EE with respect to the total transmitting power, the JIOO algorithm solves the problem via twolayer iterations, where the inner layer aims to find the maximum EE and its derivative for a given total transmitting power, and the outer layer targets to search for the total transmitting power which results in the maximum EE by bisection search. The computational complexity is mainly determined by the outer layer iterations. Unfortunately, the number of the outer layer iterations becomes infinite when tolerance approaches zero. Different from this approach, the complexity of our algorithm depends on the finite number of the feasible subcarrier assignments. Moreover, the search space can be greatly reduced by using the sufficient condition provided by this paper. In addition, the EE waterfilling structure ensures that the results obtained by our algorithm are the exact optimal solutions other than approximation ones.
The rest of the paper is organized as follows. In Section 2, we describe the system model and formulate the EE maximization problem for the downlink OFDMA networks. The optimal energyefficient power allocation for a given subcarrier assignment and the sufficient condition for the optimal energyefficient subcarrier assignment are studied in Sections 3 and 4, respectively. An optimal energyefficient power allocation algorithm and a lowcomplexity joint resource allocation algorithm are elaborated in Section 5. Section 6 provides the performance evaluation of the proposed algorithms via simulations. Finally, we conclude this paper in Section 7.
2. System Model
Consider a downlink OFDMA network with one base station and users. Let and be the set of users and subcarriers, respectively. Define the subcarrier assignment matrix where means that subcarrier is allocated to user , and otherwise . The transmitting power allocation matrix is defined as , where represents the transmitting power allocated to user on subcarrier . Then, the maximum achievable data rate of user on subcarrier is given bywhere is the bandwidth of subcarrier and denotes the normalized channel power gain of user on subcarrier . Accordingly, the overall system data rate isand the total transmitting power is
In addition to transmitting power, the energy consumption also includes circuit power which is consumed by device electronics. The circuit power is modeled as a constant , which is independent of data transmission rate [7]. Accordingly, we define EE as the amount of bits transmitted per Joule of energy; that is, , where is the reciprocal of drain efficiency of power amplifier.
In our work, we consider maximizing EE under the minimum data rate requirements, , and the total transmitting power constraint, . Accordingly, this optimization problem can be formulated aswhere (4b) indicates the minimum data rate requirement of each user, (4c) is the total transmitting power constraint, and (4d) is the constraint on subcarrier assignment to ensure that each subcarrier is only assigned to one user.
Similar to the traditional spectralefficient resource model, P1 is a mixed integer nonlinear programming problem, and it is not trivial to obtain the global optimal solution to this problem. To solve the problem, we first decompose P1 into two subproblems, which include (1) the energyefficient power allocation for a fixed subcarrier assignment and (2) the energyefficient subcarrier assignment for a given total transmitting power . Then, based on the properties of the subproblems, we develop an algorithm to find the solution of joint energyefficient power allocation and subcarrier assignment to maximize EE.
3. Optimal EnergyEfficient Power Allocation
In this section, we analyze the optimal energyefficient power allocation (EPA) based on the EE waterfilling structure. All the major results are given by some theorems. In particular, Theorems 2 and 5 demonstrate that the global optimal solution to the energyefficient power optimization problem is with the EE waterfilling structure. Theorem 3 provides the corresponding closedform waterlevel, whose optimality is proved by Theorem 4.
3.1. EE WaterFilling Structure
Given the subcarrier assignment matrix , the set of subcarriers assigned to user can be denoted by and the power vector . Then P1 is reduced to the following EPA problem:where with , and denotes the index of the user assigned on subcarrier . , and . Since the numerator of is differentiable, positive, and strictly concave function in and the denominator is positive and affine in , is a strictly pseudoconcave function with respect to [8]. Besides, is differentiable and concave for all , and is positive and affine. Therefore, according to the KKT sufficient optimality theorem [9], any feasible solution satisfying the KKT conditions is also globally optimal for P2.
When the feasible solution set of P2 is nonempty, the minimum power vector to guarantee the minimum datarate requirement of each user must be a feasible solution to P2, which can be obtained by solving the following margin adaptive (MA) problem [10]:Resorting to the Lagrange dual theory, the optimal solution to P3 is given by , where is the root of the equationWe call the lowest power waterlevel of user .
More importantly, a series of feasible solutions to P2 can be constructed based on by raising the waterlevels of some users and maintaining that of the others. To be specific, by sorting all the users in ascending order of their lowest power waterlevels such that with , the region of the promotable waterlevel can be divided into intervals, that is, , named as waterlevel rise interval. For the th interval, given called EE waterlevel, if we raise the waterlevels of the first users to , while maintaining that of other users, a new feasible power solution can be obtained by the following EE waterfilling structure:where and . The corresponding data rate of each user satisfiesIt can be found that the data rate of each user with the same waterlevel is greater than the minimum requirement, while that of the other users is equal to the minimum requirement.
It is noteworthy that since in (8) can be obtained by solving (7) and hence can be expressed solely as a function of the waterlevel , accordingly, P2 can be transformed into a single variable problem. Furthermore, for any given total transmitting power , P2 is equivalent to P4 shown in the following. If the optimal solution to P4 has the EE waterfilling structure, we can deduce that P2 must be maximized at a power vector with the EE waterfilling structure. It can be proved by Theorem 1:where .
Theorem 1. Given , the optimal solution to P4 has the EE waterfilling structure (8).
Proof. Since P4 is a convex programming, the optimal solution must satisfy the KKT conditions; that is, there exist scalars and such thatAccording to (11), we haveBesides, based on the complementary slackness conditions (12), we can get that if . Let and . Therefore, can be further expressed as follows:where and . To maximize the overall data rate, the power allocated to the users in should be minimized. And then . Since , we get . On the other hand, in order to satisfy . Consequently, and . We can now conclude that has the EE waterfilling structure as (8).
According to the waterlevels of , the subcarrier set can be further divided into three subsets: , , and . The partial derivatives of have the following properties, which will be used in the proof of the following theorems.
Property 1. If is a feasible solution to P2, then (a) (b)if , where .
3.2. EE WaterLevel Interval
The minimum power vector is also with the EE waterfilling structure whose waterlevel is and the total transmitting power is . On the other hand, suppose is the optimal solution to P4 when ; according to Theorem 1, we have , where is the EE waterlevel. Since the feasible region of P1 is nonempty, the total transmitting power must satisfy . Hence, the corresponding power vector with the EE waterfilling structure in the feasible region of P1 should be subject to . Based on the strict pseudoconcavity of , we have the following theorem.
Theorem 2. Assume and . is the optimal solution to P2 if and only if , and is the optimal solution to P2 if and only if .
However, when and , whether there exists a EE waterlevel to make optimal is still not answered. We should study the relation between and the EE waterlevel .
3.3. Analytical Expression of
According to (8), is a piecewise function and is its discontinuity point. If we sort all discontinuity points in an ascending order such that , the interval can be divided into subintervals, that is, with . Moreover, is continuous when ( is named as the continuous power interval hereafter). To simplify our analysis and get the closedform , we further assume that the waterlevel rise interval is , and the continuous power interval is . Then, according to the definition of the subcarrier subset, we have , and let . In this case, can be transformed into a continuous function of the waterlevel ; that iswhere The domain of is . It is noteworthy that and are constant as long as the waterlevel rise interval and continuous power interval are determined.
Similar to , is also strictly pseudoconcave, and the firstorder derivative iswhereAccording to the firstorder optimality condition, a stationary point of is the root of the equation . The closedform expression of is given by Theorem 3.
Theorem 3. If there exists a stationary point in the domain of , its closedform expression is given bywhere represents the Lambert function.
The proof of the theorem can be found in [11].
3.4. EEOptimal WaterLevel
According to the strict pseudoconcavity, is maximized at the stationary point . However, whether the corresponding is the global optimal solution to P2 still needs to be verified.
Theorem 4. If is the stationary point of given by (21), then is the global optimal solution to P2.
Proof. Since , it can be verified that is a feasible solution to P2. Besides, because is a stationary point of , we have . According to Property 1, we can show that . Then it can be verified that satisfies the KKT conditions. Hence is the optimal solution to P2.
In addition, the existence of is proved by Theorem 5.
Theorem 5. If neither nor is the optimal solution to P2, there must exist such that is the optimal solution to P2.
Proof. According to the intermediate value theorem, to prove Theorem 5, we should show that there must exist a continuous power interval such that .
In fact, if neither nor is the optimal solution to P2, we can verify that and according to Property 1. Assume is divided into waterlevel rise intervals. It can be proved that there must exist a waterlevel rise interval such that . If there does not exist such an interval, it can be deduced that , which yields a contradiction.
Furthermore, assume that there are discontinuity points in such that . Similarly, there must exist an interval such that among continuous power intervals. Hence, there must be a satisfying . Based on Theorem 4, is the optimal solution to P2.
4. Optimal EnergyEfficient Subcarrier Assignment
In this section, we will provide a sufficient condition for the optimal energyefficient subcarrier assignment (ESA) based on the relation between EE and the total transmitting power . By utilizing this sufficient condition, a quick search method can be devised to obtain the optimal ESA, which will be described in the next section.
According to (4a)–(4e), a feasible ESA can be obtained by solving a rate adaptive (RA) problem for a given total transmitting power . Moreover, the maximum EE can only be achieved at one of three different total transmitting powers, including two boundary points ( and ) and a stationary point of the perturbation function of P1 [5]. To obtain the optimal ESA, it should first determine , which is an unknown value. Unfortunately, is difficult to be determined and only an approximation can be found by the iterative algorithms [5]. Therefore, only the suboptimal ESA can be obtained according to the approximate .
On the other hand, based on the EE waterfilling structure discussed in the previous section, the optimal ESA can be obtained by calculating the exact optimal EE for every feasible subcarrier assignment and then selecting the one with the maximum value. This exhaustive search is prohibitive for large and in a practical system. However, combining the EE waterfilling framework and the property of the perturbation function of P1, a sufficient condition for the optimal ESA can be established to greatly simplify the search.
Define as the perturbation function of P1, where represents the maximum overall data rate of the rate adaptive (RA) problem [10] with user data requirements for a given total transmitting power . Then, we have the following.
Theorem 6. For a feasible subcarrier assignment , the optimalEE waterlevel and the maximum EE are given by (21) and (17), respectively. Correspondingly, the total transmitting power and the overall data rate . If and , is the optimal ESA of P1.
Proof. To prove Theorem 6, we should show derived from the optimal EPA for the fixed subcarrier assignment is a stationary point of . In this case, consider the derivative of If , we can get . Since is the stationary point of , it has that based on (20). On the other hand, [11]. Then, ; that is, is a stationary point of . When , EE is maximized at the stationary point of . Therefore, is the optimal ESA.
Based on Theorem 6, the following proposition can be easily verified.
Proposition 7. If , then , where is the subcarrier assignment obtained by solving the RA problem with .
5. Joint EnergyEfficient Resource Allocation Algorithm
Based on the analysis in the previous sections, we develop an optimal energyefficient resource allocation algorithm with low complexity to solve P1, named as joint energyefficient resource allocation (JERA) algorithm. Different from the existing algorithms proposed in [4–6], the JERA algorithm consists of two layers to iteratively perform subcarrier assignment and power allocation so as to achieve the optimal solution. The aim of the outer layer is to find a feasible subcarrier assignment for a given total transmitting power, and the inner layer is in charge of energyefficient power allocation based on the obtained subcarrier assignment. Based on the EE waterfilling framework, the optimal EPA can be easily obtained in the inner layer. Meanwhile, the outer layer finds a series of subcarrier assignments such that the optimal EEs increase monotonously.
The JERA algorithm is shown in Algorithm 1. From Line (1) to Line (8), the algorithm first verifies whether the boundary point ( or ) is optimal by utilizing the strict quasiconcavity of shown in [5]. If neither is optimal, the algorithm will search the optimal ESA and calculate the corresponding EPA, as shown from Line (10) to Line (15). Specifically, starting from a feasible subcarrier assignment , JERA calculates the optimal EPA and hence a new , which is further used to find a new subcarrier assignment by solving a RA problem [10]. This procedure repeats until the sufficient condition in Theorem 6 is satisfied.


It is noteworthy that the sequence of EE generated by the JERA algorithm is monotonously increasing according to Proportion 1. Besides, the optimal EE is upper bounded, so that JERA algorithm must converge to the optimal solution.
The OEPA algorithm in Line (13) is devised to get the exactly optimal power allocation based on the EE waterfilling framework, which is summarized as Algorithm 2. According to Theorems 2 and 5, the OEPA algorithm consists of three phases. In the first phase, the minimum waterlevels ’s are obtained. If the power vector derived from ’s is not optimal, the algorithm proceeds to the second phase, where the maximum waterlevel under the total transmitting is found. If the maximum power vector with is optimal, the algorithm returns. Otherwise, it searches sequentially for the optimal waterlevel in the feasible region in the last phase. Since the number of the searched waterlevel intervals only depends on that of the users and the subcarriers, the average computational complexity of the OEPA algorithm is in the worst case, while that of the BPA (bisectionbased power adaptation) algorithm [5] is proportional to , where is the convergence tolerance. More importantly, the value of EE obtained by the OEPA algorithm is the exact optimum, whereas the BPA algorithm can only provide an approximation result.
6. Performance Evaluation
In this section, simulation results are given to verify the theoretical analysis and the performance of the proposed algorithms. In our simulation, the number of data subcarriers is set to be 72 and the bandwidth of each subcarrier is 15 kHz [5]. The block Rayleigh fading channel model is considered and the OkumuraHata path loss model is followed; that is, in decibels, where is the distance between transmitter and receiver in kilometers. The standard deviation of shadowing is 7 dB, and the thermal noise spectral density is dBm/Hz [4]. The circuit power is 20 W and the maximum transmitting power is 40 W for the base station [12]. The drain efficiency of power amplifier is assumed to be 38% [4]. Each user in the simulation has the same minimum rate requirement of 100 kbps.
First, we compare the performance of the OEPA algorithm with the other two algorithms: the BPA algorithm [5] and the MWF (multilevel water filling) algorithm [10] for a fixed subcarrier allocation. Although the MWF algorithm is a classical SEbased scheme rather than a EEbased scheme, it is used as a benchmark to measure the difference in the energy efficiency between the two classes of scheme. In this simulation example, the number of users is set to 30. The users are uniformly distributed in a circle centered at the BS with a variable radius. The results in Figures 1(a) and 1(b) show that the average EE and the system throughput of all the three algorithms decrease with the channel power gain. This is due to the fact that the average channelgaintonoise ratio (CNR) of each user decreases with the increase of the distance between the user and the BS, such that more power is required to combat with the severe channel fading to satisfy the data rate constraint of each user. Figure 1(c) illustrates that the SEbased scheme, that is, the MWF algorithm, always depletes all the power to maximize the system throughput. On the contrary, the power consumption of the EEbased scheme, that is, the OEPA algorithm and the BPA algorithm, is adaptive to the channel condition. It is worth noting that there exists a gap between the EE obtained by the OEPA algorithm and the BPA algorithm. The difference is especially 10.15 bit/joule on average when . This stems from the fact that the OEPA algorithm always finds the exact optimal solution, while the BPA algorithm only approaches the optimal value by iterative method. More importantly, the computational complexity of the OEPA algorithm is much lower than that of the BPA algorithm, which is demonstrated in the following simulation example.
(a) Evaluation and comparison of average EE
(b) Evaluation and comparison of average throughput
(c) Evaluation and comparison of average transmitting power
In order to compare the computational complexity, assume that there are 10 users evenly distributed in the network and 10000 independent experiments (with different user location) are conducted in the MATLAB environment. Figure 2 shows the CDF of the required CPU time for convergence in different algorithms, where represents the error tolerance of BPA algorithm. From Figure 2, it is obvious that the convergence speed of the OEPA algorithm is superior to that of the BPA algorithm. In addition, the complexity of the BPA algorithm increases evidently when the error tolerance becomes tighter. This is due to the fact that the OEPA algorithm can obtain the exact optimal solution by checking at most continuous power intervals in the worst case based on the closedform expression of EPA. On the contrary, the BPA algorithm is to search bidirectionally for the optimal transmitting power, which results in a higher computational complexity. On average, the CPU time for convergence of the OEPA algorithm is about 15.64% and 13.03% of that of the BPA algorithm with and , respectively.
Moreover, in order to verify the optimality of the proposed JERA algorithm, we compare the EE obtained by JERA with the global optimum obtained by exhaustive search. In this case, we consider a system with 9 subcarriers and 3 users to reduce the complexity of exhaustive search. As shown in Figure 3, the achieved optimal EE in both algorithms decrease with the distance between the BS and users. Moreover, the two curves match with each other very well. It demonstrates the proposed ESA search method based on Theorem 6 is effective and the JERA algorithm can obtain the global optimal ESA and EPA simultaneously.
In addition, we compare the performance of the JERA algorithm with that of the JIOO algorithm [5] with different number of users and subcarriers. The result is shown in Table 1; it can be observed that the EE of the JERA algorithm is superior to that of the JIOO algorithm due to the optimality of the solution obtained by the JERA algorithm. More importantly, the convergence rate of the JERA algorithm is significantly faster than that of the JIOO algorithm. Specifically, the number of iterations for convergence of the JERA algorithm is less than 5 in average, while the JIOO algorithm requires at least 29 iterations in average to approximate the optimum. It is worth noting that each iteration in both algorithms needs to solve a RA problem with rate requirements and total transmitting power constraint. Despite of the nonconvexity of this type RA problem, it has been proven that it can be solved efficiently by the Lagrange dual decomposition method with zero duality gap [13].
7. Conclusion
In this paper, we investigated the EE maximization problem under both the user rate requirements and the transmitting power constraint. Utilizing the EE waterfilling structure, we obtain the closedform of the optimal EPA. The sufficient condition for optimal ESA is also derived based on the relation between EE and power. Furthermore, we propose a lowcomplexity algorithm with joint ESA and EPA to address the energyefficient resource allocation in downlink OFDMAbased networks. Simulation results show that the proposed algorithm achieves the optimal energyefficient resource allocation with significantly reduced computational complexity compared with the iterative methods.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the National S & T Major Project of China under Grant no. 2014ZX03004003, the Fundamental Research Funds for the Central Universities (no. ZYGX2013J009), China, and EU FP7 Project CLIMBER (PIRSESGA2012318939).