Advocated by the dual use of radio frequency signals, wireless energy transfer (WET) has attracted much attention for improving the system energy efficiency. In this context, simultaneous wireless information and power transfer (SWIPT) and wireless powered communication networks (WPCNs)[3, 4, 5] have been extensively studied in the literature. Moreover, since electromagnetic waves decay quickly over distance, energy beamforming is generally designed to achieve efficient WET. In WPCNs, a power station (PS) transfers wireless energy to some low-power users with a single antenna due to the hardware constraint. Afterwards, the users transmit information signals with the harvested energy. For the multiple users scenario, the signals are transmitted typically based on time division multiple access (TDMA) as in[3, 4, 5]. However, the spectrum efficiency can be greatly improved with appropriate interference management methods by allowing multiple users to transmit signals simultaneously.
In this paper, we consider the D2D communications underlaying WPCNs, where the low-power D2D transmitters (D2D-Txs) with one antenna must harvest energy from the PS equipped with multiple antennas before transmitting signals. In the considered scenario, all D2D-Txs transmit information signals simultaneously over the same spectrum resource. Our aim is to maximize the sum throughput via joint time scheduling and power control, while satisfying the energy causality constraints. For solving the formulated non-convex nonlinear problem, we develop a throughput maximization algorithm, where the problem is first transformed into a nonlinear fractional programming problem with a tactful reformulation. Then, an iterative algorithm is designed to address the equivalent problem by leveraging D.C. (difference of two convex functions) programming. In-depth simulations are conducted to evaluate the throughput performance under various system parameter configurations.
The rest of this paper is organized as follows. In Section II, we illustrate the system model and formulate the optimization problem. Then, we develop a throughput maximization algorithm in Section III. In Section IV, we present simulation results to verify the effectiveness of the proposed algorithm. Finally, we conclude the paper in Section V.
Ii System Model and Problem Formulation
Consider a WPCN with a PS equipped with antennas and low-power D2D pairs denoted by . The D2D pair carries a single antenna due to the size and cost constraints, such as the sensor node. With no embedded energy supply, each D2D-Tx first harvests energy from wireless signal transmitted by the PS (i.e., WET phase). Then, they utilize the harvested energy to transmit information signals to their intended receivers in the wireless information transmission (WIT) phase. According to the harvest-then-transmit protocol, in each block denoted by , the first amount of time, , is assigned to harvest energy for all D2D pairs, while the followed amount of time in the same block is assigned to transmit information signals. Followed by, we consider a normalized unit block time in the sequel without loss of generality. Then, there is . All the users considered in this paper operate on a single spectrum band.
In the WET phase, the transmitted signal is given by , where is the transmit power of the PS, and the beamformer is designed to improve the energy transfer efficiency and subject to . Let represent the
dimensional energy transfer channel vector between the PS and-th D2D-Tx. The energy harvested from the noise can be ignored since the noise power is usually much smaller than that of the PS. Therefore, the energy harvested at the -th D2D-Tx is given by with energy conversion efficiency . We apply the asymptotically optimal energy beamforming proposed in as
where controls the energy allocation among multiple D2D pairs and . In this paper, equal weight is designed for all D2D pairs, i.e., . Notice that other energy beamforming schemes as invetigated in can be employed and the proposed algorithm still works.
Denote as the channel power gain from the -th D2D-Tx to its receiver. The channel power gain of the interference link from the -th D2D-Tx to the -th D2D receiver (D2D-Rx) is denoted by . Since all D2D-Txs transmit information signals simultaneously over the same spectrum resource, the signal to interference plus noise ratio at the -th D2D-Rx is as follows:
where is the transmit power of -th D2D-Tx and is the noise power. The achievable throughput at the -th receiver in bits/second/Hz is thus given by
Intuitively, a D2D pair closer to the PS can harvest more energy in short time and vice versa, which potentially results in various energy constraints for different D2D pairs. To character this difference, the transmit power and time are jointly optimized here. The aim is to maximize the sum throughput of all D2D pairs via time scheduling and power control, while satisfying the energy causality constraints. Thus, the optimization problem can be formulated as the following:
where represents the non-ideal circuit power consumption (e.g., AC/DC converter, analog amplifier, and processor)[3, 4]. guarantees that the consuming energy by any D2D-Tx cannot exceed its harvested energy. and are the time and power control constraints. However, the feasible region of is non-convex, which means that the standard convex optimization methods cannot be used to efficiently solve this problem. Even if the time lengths of WET and WIT have been fixed, the investigated problem is still non-convex and hard to be addressed. In the next section, we propose an efficient throughput maximization algorithm by exploiting the problem structure.
Iii Throughput Maximization Algorithm for D2D Communications
The optimal solution to the problem (II) is generally difficult to be obtained since there are multiple local optima due to the non-convex nonlinear property. To this end, the formulated non-convex problem is first transformed into a nonlinear fractional programming problem with a tactful reformulation. Then, an iterative algorithm is designed to solve the equivalent problem by leveraging D.C. programming.
The time utilization is illustrated by the following Lemma.
Lemma 1: The optimal solution to the problem (II) is achieved if and only if all the time is used, i.e., .
To prove Lemma 1, we assume that is the optimal solution satisfying and , which means that there is remaining time available denoted by . If we can find a feasible solution to the optimization problem (II) in the remaining time, it demonstrates that the system throughput can also be improved. In other words, the solution is not the optimal solution.
The remaining time can also be divided into two parts according to the ratio . The first amount of time of the remaining time is used to harvest energy, while the remaining is to transmit information signals. Since is a feasible solution, there must be
Then, the following constraint also holds true:
This means that the solution is a feasible solution in the remaining time with , , and . So the remaining time can be used to improve the system throughput, which contradicts the assumption. The Lemma 1 has been proved. ∎
Based on Lemma 1, the constraint in problem (II) can be transformed as follows:
It is non-convex with respect to and , which hinders the application of standard convex optimization techniques. Although it can be transformed into a convex function by geometric programming, the resulting objective function will be much more complicated. To this end, replace and . The convenience of this replacement will be shown later. Thus, the constraint in problem (II) can be rewritten as
This alternation makes the constraint become linear. With this reformulation, the optimization problem (II) is equivalent to
It can be observed that the non-convex constraints are transformed into convex functions with a tactful reformulation. Moreover, the optimization problem (III) can be seen as a nonlinear fractional programming. Therefore, we try to search an optimal solution to the problem (II) by solving the equivalent problem (III).
Denote as the optimal solution of the considered problem (III), which is given by
The following Lemma 2 provides a guidance on designing an iterative approach to solve the problem (III).
Lemma 2: The optimal solution is achieved if and only if
Lemma 2 can be proven by following a similar approach as in. It is shown that the original fractional form in problem (III) can be transformed into a subtractive form with an equivalent solution, which indicates that an iterative algorithm can be designed to solve this problem. Specifically, initializing from a given , we should solve a sequence of following problems:
The overall procedure for solving the optimization problem (III) is presented in Algorithm 1. The convergence to the optimal solution is guaranteed.
Theorem 1: As long as the number of iterations is sufficiently large, the proposed algorithm will eventually approach the optimal solution.
See Appendix A for details. ∎
Although we have designed a framework to efficiently solve the problem (III), the optimization problem (III) is hard to be solved due to the non-convexity of the objective function. The optimal solution to this problem cannot be obtained so far. In the sequel, we develop an iterative algorithm to get a suboptimal solution to the problem (III).
Denote the following concave functions
The objective function in problem (III) can be expressed as
It can be observed that the objective function is the difference of two concave functions. A series of non-decreasing solutions can be obtained by iteratively optimizing the lower bound of objective function, which is given by the following Lemma 3.
Lemma 3: Given , the following function is a tight lower bound of the objective function in problem (III):
where the -th component of the is given by
Since is concave, based on the first-order condition of a concave function, we have . Thus, . Moreover, if , there is . So provides a tight lower bound of function . ∎
According to Lemma 3, an iterative algorithm can be developed to solve the optimization problem (III). In particular, initializing from a given , we can iteratively solve the following convex problem using standard convex optimization techniques:
There are many convex optimization techniques and they have been widely investigated. Due to the page limit, we omit the procedure of solving problem (III). Finally, the overall procedure for solving the optimization problem (III) is presented in Algorithm 2.
Lemma 4: The resulting values of Algorithm 2 at each iteration are non-decreasing, and the convergence is guaranteed.
Let be the solution at -th iteration. Since provides a lower bound for , there is
Moreover, because is the optimal solution at the th iteration, we have
Therefore, there is
It can be observed that the resulting values are non-decreasing at each iteration. Further, it must be upper bounded by the optimal value of (III). Thus, the convergence is guaranteed. ∎
Complexity Analysis: The proposed throughput maximization algorithm contains two-layers iterations. The outer-layer iteration is the update procedure defined by (III). The inner-layer iteration is to acquire a lower bound by solving the D.C. programming problem at given . The convex problem (III) can be solved with the complexity of by standard convex optimization techniques, such as interior point method, where is the number of D2D pairs. Assume that
are the computations in the two-layer iterations. The total complexity can be roughly estimated as.
Notably, we consider a centralized network in this paper, where the PS performs the proposed algorithm. The channel state information can be estimated by the PS according to the channel reciprocity.
Iv Simulations and Discussions
In this section, we perform in-depth simulations to evaluate the performance of the proposed algorithm in a m area, where multiple D2D pairs are randomly located and the maximum distance between D2D-Tx and D2D-Rx is m. The channel power gain is modeled as [4, 9], whereis the distance between the transmitter and receiver, and represents the path-loss exponent. Unless specified otherwise, the bandwidth is 1 MHz and noise power spectral density is dBm/Hz. The transmit power and number of antennas are 1 W and 10 for the PS, respectively. The energy conversion efficiency and circuit power consumption are 0.5 and 0.1 W. In all simulations, is set to start the algorithm and all results are averaged over 100 realizations.
The sum throughput versus the transmit power of the PS is shown in Fig. 1. For comparison, we also provide an energy transfer scheme without beamforming, namely the omnidirectional energy transfer (OET), and a TDMA-based algorithm where multiple users harvest energy and then transmit information signals based on TDMA. It can be observed that the proposed algorithm outperforms the OET and TDMA-based algorithm in all cases. Furthermore, we can observe that the growth rate gradually becomes slower as the transmit power increases. This is due to the fact that the mutual interference among D2D pairs dominates the system with sufficiently large transmit power. In addition, The final sum throughput would be better if the maximum distance between D2D-Tx and D2D-Rx is reduced. The reason is that smaller maximum distance results in better channel state.
In Fig. 2, the sum throughput is plotted against the circuit power consumption . It can be observed that the sum throughput decreases with an increasing circuit power consumption. Meanwhile, the throughput gain between the proposed algorithm and the TDMA-based algorithm is smaller. The reason is that D2D pairs have little energy for information transmission and some D2D pairs may even stop working since they have not enough energy. The impact of number of antennas is further investigated and simulation results are shown in Table I. The plot confirms the intuition that the sum throughput grows as more antennas are added at the PS since more antennas can make use of the spatial resource to improve diversity gain.
|Settings||Number of antennas ()|
Table II shows the WIT time for different settings, which characterizes the time split between energy harvesting and data transmission. The data transmission time grows with an increasing transmit power of the PS. This is due to the fact that the D2D-Txs can harvest more energy with high transmit power, which reduces the energy harvesting time. On the other hand, more D2D pairs not always results in an increasing WIT time since the mutual interference becomes serious, where some D2D pairs experienced serious interference do not work. It can be also observed that the value ( W) is smaller, which results from locally optimal points of non-convex optimization.
|Settings||Transmit power of the PS (W)|
In this paper, we investigated the resource allocation scheme for D2D communications underlaying WPCNs, where the harvest-then-transmit protocol was employed. We tried to maximize the sum throughput of all D2D pairs while satisfying the energy causality constraints. The considered joint time scheduling and power control problem was formulated as a non-convex optimization problem and then it was transformed into a nonlinear fractional programming problem with a tactful reformulation. By leveraging D.C. programming, a suboptimal solution of the non-convex problem can be obtained by iteratively solving a sequence of convex problems. In-depth simulations were conducted to validate the effectiveness of the proposed algorithm.
The proof of Theorem 1
The following proof follows the proof of algorithm convergence in paper. For notational convenience, denote and . For any feasible solution and , there is which means that always holds true for any feasible solution . Furthermore, let and as the optimal solutions for and , respectively. Then, there is If , So is a strictly monotonic decreasing function.
Denote as the optimal solution at the -th iteration and the according value is . From the iterative algorithm, we know . Then, we have Because , so there is . we can show that as long as the number of iterations is large enough, will eventually approach zero with the increasing since is a strictly decreasing function. ∎
-  X. Lu, P. Wang, D. Niyato, D. I. Kim, and Z. Han, “Wireless networks with RF energy harvesting: A contemporary survey,” IEEE Commun. Surv. & Tutor., vol. 17, no. 2, pp. 757-789, Nov. 2015.
-  K. Huang and E. Larsson, “Simultaneous information and power transfer for broadband wireless systems,” IEEE Trans. Signal Process., vol. 61, no. 23, pp. 5972-5986, Dec. 2013.
-  Z. Hadzi-Velkov, I. Nikoloska, G. K. Karagiannidis, and T. Q. Duong, “Wireless networks with energy harvesting and power transfer: Joint power and time allocation,” IEEE Signal Process. Lett., vol. 23, no. 1, pp. 50-54, Jan. 2016.
-  S. Pejoski, Z. Hadzi-Velkov, T. Q. Duong, and C. Zhong, “Wireless powered communication networks with non-ideal circuit power consumption,” IEEE Commun. Lett., vol. 21, no. 6, pp. 1429-1432, June 2017.
-  H. Ju and R. Zhang, “Throughput maximization in wireless powered communication networks,” IEEE Trans. Wireless Commun., vol. 13, no. 1, pp. 418-428, Jan. 2014.
-  S. Bi, C. K. Ho, and R. Zhang, “Wireless powered communication: Opportunities and challenges,” IEEE Commun. Mag., vol. 53, no. 4, pp. 117-125, Apr. 2015.
-  W. Huang, H. Chen, Y. Li, and B. Vucetic, “On the performance of multi-antenna wireless-powered communications with energy beamforming,” IEEE Trans. Veh. Technol., vol. 65, no. 3, pp. 1801-1808, Mar. 2016.
-  W. Zhao and S. Wang, “Resource sharing scheme for device-to-device communication underlaying cellular networks,” IEEE Trans. Commun., vol. 63, no. 12, pp. 4838-4848, Dec. 2015.
-  G. Yang, C. K. Ho, R. Zhang, and Y. L. Guan, “Throughput optimization for massive MIMO systems powered by wireless energy transfer,” IEEE J. Sel. Areas Commun., vol. 33, no. 8, pp. 1640-1650, Aug. 2015.
-  S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004.
-  W. Dinkelbach, “On Nonlinear Fractional Programming,” Manage. Sci., vol. 13, no. 7, pp. 492-498, Mar. 1967.
-  H. Wang, G. Ding, J. Wang, L. Wang, T. A. Tsiftsis, and P. K. Sharma, “Resource allocation for energy harvesting-powered D2D communications underlaying cellular networks,” in Proc. IEEE International Conference on Communications (ICC), May 2017, pp. 1-6.