A Fair Resource Allocation Algorithm for Data and Energy Integrated Communication Networks

With the rapid advancement of wireless network technologies and the rapid increase in the number of mobile devices, mobile users (MUs) have an increasing high demand to access the Internet with guaranteed quality-of-service (QoS). Data and energy integrated communication networks (DEINs) are emerging as a new type of wireless networks that have the potential to simultaneously transfer wireless energy and information via the same base station (BS). This means that a physical BS is virtualized into two parts: one is transferring energy and the other is transferring information. The former is called virtual energy base station (eBS) and the latter is named as data base station (dBS). One important issue in such setting is dynamic resource allocation. Here the resource concerned includes both power and time. In this paper, we propose a fair data-and-energy resource allocation algorithm for DEINs by jointly designing the downlink energy beamforming and a power-and-time allocation scheme, with the consideration of finite capacity batteries at MUs and power sensitivity of radio frequency (RF) to direct current (DC) conversion circuits. Simulation results demonstrate that our proposed algorithm outperforms the existing algorithms in terms of fairness, beamforming design, sensitivity, and average throughput.


Introduction
The strength of network virtualization is upsurge, such as software defined networking (SDN) and collaborative radio access networks (C-RAN).For instance, SDN is expected to transform the way services are created, sourced, deployed, and supported [1][2][3].This paper discusses another type of virtualization, namely, a base station, being virtualized into providing not only information transferring but also energy transferring to charge mobile devices.This is largely driven by the fact that mobile devices, while getting more powerful in processing and networking, exhaust their battery more quickly.
To address this challenge, this paper utilizes energy harvesting into wireless communications, namely, the socalled DEINs (Data and energy integrated communication networks).With the development of energy harvesting (EH) technologies and wireless energy transfer (WET) techniques, the DEIN becomes an emerging trend focusing on the study of wireless power and information cooperation communications [4][5][6][7].Compared with simultaneous wireless information and power transfer (SWIPT) which mainly focuses on the physical layer, DEINs focus on the whole network system, resource allocation, and protocols design in different layers.The typical architecture of a DEIN is shown in Figure 1, which has three major components, that is, a virtual energy base station (eBS), a data base station (dBS), and massive mobile users/mobile devices.The wireless communication is divided into two parts.The virtual eBS first transfers energy to multiple mobile users (MUs) that do not have embedded energy sources via downlink (DL), and then MUs use the harvested energy to perform uplink (UL) wireless information transmission (WIT) to the dBS.
No existent works have taken into account the power sensitivity of RF-DC circuits when DL transfer energy in DEINs, which can lead to a falsely higher data rate when received RF signals cannot be converted into DC (i.e., energy transfer) if their power level is lower than the power sensitivity of an RF-DC circuit [8].Besides, none of these works have considered the possibility of energy overflow or the opportunities for users to optimize the use of harvested energy across UL WIT slots.It has been shown that a user using all available energy for WIT in each slot achieves a lower data rate than uniformly distributing energy between energy arrivals [9][10][11].Therefore, new dynamic time allocation schemes are needed since not all MUs can harvest energy in every slot, which take into consideration the policy that the energy harvesting of every user should not overflow in every DL WET phase and the energy harvested in the DL WET phase of former slot may be used in the UL WIT of the next.
In this paper, we devise a fair resource allocation algorithm for a DEIN based on dynamical time-slot division to guarantee the received power fairness among all MUs, by maximizing the minimum average DL WET power via optimizing the DL energy beamformer.In particular, considering the power sensitivity of the RF-DC conversion circuits, the virtual eBS transfers energy to these MUs whose estimated received power is larger than the corresponding certain threshold .The proposed algorithm benefits the following features.
(1) Fairness.Since there are more than one user in our scenario, fairness among all MUs should be taken into account naturally.While each MU is allocated the same slot for UL transmission by space-division-multiple-access (SDMA) leading to the fact that the increase or decrease of slot length would lead to similar trend of each MU's throughput, the fairness of the model is mainly reflected by the WET energy beamforming.To achieve the fairness, we aim to allocate the optimal beamformer to maximize the minimum received power among MUs, which can further lead to the fairness of throughput on the whole.
(2) Sensitivity.We study the UL transmission powered by WET in a DEIN and model the WET of every user as a Bernoulli process with different probability while taking into account the different sensitivity of RF-DC circuits and calculate the accurate WET probability of different users.
(3) Throughput Performance.Considering the sensitivity of RF-DC circuits, that is, the WET probability, we propose a new low complexity resource allocation algorithm to achieve near optimization throughput with finite capacity batteries, which employs the zero-forcing (ZF) based receive beamforming in the UL information transmission [12].
The rest of this paper is organized as follows.Section 2 illustrates some related works about DEINs in recent years.Section 3 presents a multiantenna DEIN model based on virtualization and formulates problems based on fairness.Section 4 presents a dynamic resource allocation to make a near throughput optimization for this problem.Section 5 provides simulation results to compare the performances of proposed solutions with exiting algorithms.Finally, Section 6 concludes the paper.
Notations.All lowercased and uppercased boldface letters represent vectors and matrices, respectively.Let tr(), det(),  −1 , and   denote the trace, determinant, inverse, and Hermitian of a symmetric matrix , respectively.C and R denote the set of complex and real matrices of size  × , and C and R denote the set of complex and real vectors of size  × 1, respectively.All the log(⋅) functions are of base 2 by default and ln(⋅) stands for the natural logarithm.All letters at the right bottom of different variables can be explained by the following:  shows the th slot and  is the different users.

Related Works
In this section, we introduce the related energy harvestingbased data-and-energy transmission techniques.
There are so many works focusing on designing resource allocation schemes for different networks.In [13], an energy-efficient context-aware resource allocation problem in caching-enabled ultradense small cells is investigated.In [14], a new analytical performance model is developed to evaluate the QoS of multihop cognitive radio networks.Moreover, we introduce some current researches on DEINs in the following.
Moreover, another problem, addressing the joint design of DL energy transfer and UL information transmission, is worth investigating [12,[23][24][25].UL WIT powered by DL WET in DEINs was studied in [23].Assuming perfect channel state information (CSI), the single-user DEIN scenario was studied in [23].The authors maximized the energy efficiency of uplink WIT by jointly optimizing the time duration and transmitted power for downlink WET.Besides, some work is focusing on resource allocation for optimization throughput in DEINs.In [24], a DEIN with single-antenna access point  (AP) and users has been studied for joint DL energy transfer and UL information transmission.A harvest-then-transmit protocol was proposed and orthogonal time allocations for the DL energy transfer and UL information transmissions of all users are jointly optimized to maximize the network throughput.In [12], a DEIN where one multiantenna AP coordinates energy transfer and information transfer to/from a set of single-antenna users was studied, and the author maximizes the minimum throughput among all users by a joint design of the DL-UL time allocation, the DL energy beamforming, and the UL transmit power allocation, as well as receive beamforming.A WET enabled massive MIMO system with imperfect CSI was studied in [25], where it is based on time-division-duplexing (TDD) protocol and each frame is divided into three phases: the UL channel estimation phase, the DL WET phase, and the UL WIT phase.
But all these works take the idea that the MUs transfer all the received power in DL, ignoring the power sensitivity of RF-DC circuits [9][10][11].Besides, no works take the advantage of the battery storage capacity to make a dynamic power allocation, which can have a higher throughput.Even though some works use the parameter (i.e., the battery storage capacity), they also overlook the fact that the capacity cannot overflow when making the resource allocation.In the following, we are going to discuss these questions.

System Model and Problem Formulation
In this section, we present a dynamical time slotted transmission scheme and analyze the DL WET phase and UL WIT phase.Finally, we formulate problems based on fairness.
We provide a DEIN model consisting of a BS with  antennas and  single-antenna MUs with finite battery capacity denoted by   ( = 1, . . ., ), as shown in Figure 2. It is assumed that  ≥ .It should be noted that the BS is virtualized into two parts, that is, a virtual eBS and a dBS, as shown in Figure 1.Each MU uses the harvested energy from DL WET phase via beamforming of the virtual eBS to power its UL information transmission via the dBS, under the assumption that two BSs and all MUs are perfectly synchronized and there is no other energy source of MUs.The total capacity of battery storage in every MU is  max .
A time slotted transmission scheme is considered as shown in Figure 3.We assume that the period of DL DIT in every slot can be neglected since there is little information to transmit in the DL phase.So, each slot has a constant period , consisting of two phases, namely, the DL WET phase of duration  ⋅  and the UL WIT phase of duration (1 − ) ⋅ , where 0 ≤  ≤ 1.The WET phase starts with several control frames, including the preamble, frame control header (FCH), DL map, and UL map.These frames define transmission parameters, such as coding schemes, available Mobile Information Systems resources, the duration of DL and UL transmission, and the WET probability (which will be defined in the following).Then, virtual energy BS transmits energy to   through wireless energy beamforming.The received power level and the energy harvested at   in slot  ( = 1, . . ., ) are denoted by  , and  , .It is worth noting that, due to the power sensitivity of RF energy harvesting circuits,   cannot harvest any RF energy if the received signal power  , is below a certain level.Thus, the received power at every MU can be used if the received power level  , is larger than the corresponding certain threshold (e.g., −10 dBm).This indicates that the WET of every MU follows a Bernoulli process with different probability   , where   stands for the probability of delivering energy from virtual eBS.Next, all MUs do UL WIT phase simultaneously via SDMA, which is powered by the energy stored in the batteries.For simplicity, we assume that  = 1 s and that the harvested energy is stored in the battery first and then used for UL information transmission.Note that since the length of control frames is much smaller than that of DL WET and UL WIT, we ignore the time duration of control frames in the following analysis.
We consider frame-based transmissions over flat-fading channels on a single frequency band [26] (i.e., which means the channel remains constant in each slot).Denote ℎ , ∈ C  * 1 as the UL channel of   in the th slot and we have where  0 denotes a constant determined by the RF propagation environment,   is the propagation distance,   is the path loss,   is Shadow fading, and  , ∈ C  * 1 is the matrix of Rayleigh fading coefficients and  , ∼ (0, 1).By exploiting the channel reciprocity, the DL transmission channel can be obtained as ℎ  , .For simplicity, we assume that   = 1 and that CSI is available at both two BSs and   .

DL WET Phase.
Assume that just one energy beam is to transmit energy from virtual energy BS to those users satisfying the probability   , since we just transmitted energy signals in the DL [27].Besides, we can allocate the optimal beamformer on the constraint that the power transferred at all of the antennas is equal to the eBS's transmitted power.Hence, we can assume the beamformer as a ratio of  which satisfies the fact that its norm is a unit value, where  is the transmitted power of eBS.Then, the practical power transferred can be denoted as the beamformer multiplied by eBS transmitted power .Moreover, ambient channel noise energy cannot be harvested.Thus, the DL received signal, received power, and harvested energy of   in the th slot are expressed as where  , ∼ (0,  2  ) is the receiver noise,   is the  × 1 beamforming and satisfies ‖  ‖ 2 = 1,  0 is the transmission signal and satisfies  2 0 ≤  max , where  max is the transmit power constraint, and   denotes the energy harvesting efficiency at   , which should satisfy 0 <   ≤ 1.For simplicity, we assume   = 1.

UL WIT Phase.
In the UL WIT phase of each slot, MUs use the harvested energy to power UL information transmission to the virtual dBS.For convenience, we assume the circuit energy consumption at   is 0. The received signal at the normal dBS in the th slot is given by where   ∈ C ×1 denotes the receiver additive white Gaussian noise (AWGN).It is assumed that   ∼ (0,  2  Γ).Besides, we assume that the normal information BS employs linear receivers to decode  , in the UL. , denotes the transmit signal of   and satisfies  2 , =   , , where   , is the transmit power of   .Specifically, let V , ∈ C ×1 denote the receive beamforming vector for decoding  , and define  = {V ,1 , . . ., V , }.In order to reduce complicity, we employ the ZF based receive beamforming in the normal information BS proposed by [12], which is not related to   and   .Define  −, = [ℎ ,1 , . . ., ℎ ,−1 , ℎ ,+1 , . . ., ℎ , ]  ,  = 1, . . ., , including all the UL channels except ℎ , .Then the singular value decomposition (SVD) of  −, is given as  −, =  , Λ    , =  , Λ , [ , Ỹ, ]  .Thus, the beamforming can be expressed as V ZF , = Ỹ, Ỹ , ℎ , /‖ Ỹ , ℎ , ‖.Then, throughput of   in bits/second/Hz (bps/Hz) can be expressed as where h, = ‖ Ỹ, ℎ , ‖ 2 .The energy consumption of UL WIT of   in the th slot is given by Let  , represent the amount of energy available in the battery of   at slot  and its updating function is as follows: There are two constraints in the UL WIT phase: the energy causality constraint and the battery storage constraint [9,11].Specifically, the energy causality constraint requires that the UL WIT can only use the energy harvested at the current and previous slots, and the battery storage constraint indicates that the energy available of   cannot exceed the maximum battery capacity at any time; that is, ) .(10) The expression (10) shows energy causality constraint, the battery storage constraint, and casualty CSI.We should note that (4) should satisfy the constraint  2 0 ≤  max .Thus, the optimization problem satisfying constraints (10) can be defined as max min lim where   = min(min  (( ,max − , )/ , ), 1) denotes the upper bound of   , preventing the battery overflow.

Dynamic Resource Allocation Algorithm
In this section, we present a dynamic resource allocation to make a near throughput optimization for the problem formulated in the former section, including DL beamformer design and allocation of WIT energy and duration.
As the optimization problem is formulated, we can divide it into two parts, one of which is the optimal beamformer designing of DL WET for fairness and the other is the optimal time and power allocation of UL WIT for throughput, since we employ the ZF based receive beamforming in the normal information BS proposed by [12], which is not related to   and   .

Design of WET Beamforming and Computation of EH
Probability.Considering fairness of the received power, we determine an optimal beamformer design to maximize the minimum received power for different MUs.To maximize the received power at   , we can set  2 0 =  max .Generally, the received power for   is denoted by where  , = ℎ , ℎ  , .Then problem on fairness could be formulated as max It is easy to relax the max.-min.process by introducing a slack variable   , and then problem ( 13) can be transformed into the following: Since the problem is convex, we can follow the approach of convex theory in [28] and consider its Lagrangian function is given by where  = [ 1 ,  2 , . . .,   ] ≥ 0 consists of Lagrange multipliers associated with received power constraints of MUs in problem (14).
The dual function of problem ( 14) is then given by The following rule can be used to determine whether problem ( 14) is feasible.For given   > 0, problem ( 14) is infeasible in the following interpretation if and only if there exists any  ≥ 0 such that () > 0.
Next, for given  ≥ 0, we would obtain () by the problem max With (12), the problem could be transformed as max where   = ∑  =1    , .Since  is a symmetric matrix, the optimal beamformer  *  could be obtained by the result of quadratic problem.After obtaining  *  for given   and , we can compute the corresponding  , ( *  ) and () in (15).If ( * ) ≤ 0, it indicates that problem ( 14) is infeasible.Therefore, we should decrease   and solve the feasibility problem in ( 14) again.On the other hand, if ( * ) ≤ 0, we can update  using ellipsoid method until  converges to  * , which denotes the maximizer of () or the optimal dual solution for problem (14).If ( * ) ≤ 0, it indicates that the problem is feasible.Then, we should increase   and solve the feasibility problem again.Finally,   * could be obtained numerically by iteratively updating   by a simple bisection search.Meanwhile, the corresponding  *  should be the optimal beamformer.Finally, we can compute the received power for each MU and determine whether they could harvest the energy, which can be successful only if  , ≥   .
To obtain the EH probability, we can iterate the above steps for  times, where  should be a large number and a different transmitting channel is generated under the rule of Rayleigh distribution for each time.The probability of   can be obtained by the ratio of successful harvested times and , which is denoted by where suc time  is the number of successful iterate times.
To summarize, the algorithm to solve the problem is given in Algorithm 1. Steps 5 to 20 give the solution to obtain an optimal beamformer in case that CSI is available.

Allocation of WIT Energy and Duration.
After obtaining the best WET beamformer, we have to determine an optimal allocation of WIT energy and duration.Considering multiple slots, the optimization of the problem would be approached if we keep the average transmitting energy in each slot.This is because uniformly distributing the energy between energy arrivals maximizes the data rate.As proposed in [29], according to the EH probability   , we take a fraction   of available energy of   to transmit information in UL duration.With this energy allocation, we only need to compute the optimal UL duration in each slot and can obtain a global optimization with multiple slots.
In the case of the system model, the total throughput of   in each slot is denoted by where h, ,   ,  , , and  , are considered as constant since they are calculated out before duration allocation.Like the method proposed in Section 4.1, the optimization problem can also be formulated as max where   = min(min  (( ,max − , )/ , ), 1) denotes the upper bound of   , preventing the battery from overflow.
Since the problem is convex, we can also solve the problem by the Lagrangian method, which is proposed in Section 4.1.What is different between two problems is how to obtain max   ∑  =1    , (  ).
Since ( 20) is concave, with  > 0, we can easily find that the following function is also concave: Consequently, one-dimensional search methods are considered to help us obtain the optimal   .
The algorithm to solve the problem is given in Algorithm 2.

Performance Results and Analysis
In this section, we focus on algorithm simulations and comparative analyses among algorithms by using MATLAB simulation tool.Moreover, we analyze simulation results based on the main problems which we have discussed in previous sections, including fairness, sensitivity, and throughput performance.

Parameter Settings.
In this section, several simulations are executed to testify the performance of the proposed algorithm in three aspects of fairness, sensitivity, and max.min.throughput.We make some comparison between the fair allocation algorithm of using fraction  of energy (UFPE) and the algorithm in [12], where MUs would use all the energy stored in their battery (UAE).Before the simulations, the following settings are used unless stated otherwise, as shown in Table 1.

Simulation Result and Analysis
5.2.1.Fairness.While each MU is allocated the same slot to transmit information uplink leading to the fact that the increase or decrease of slot length would lead to the same trend of throughput of each MU, the fairness of the model is mainly reflected on the WET energy beamforming.We compare our algorithm of allocating beamformer by channel in UFPE with the algorithm of allocating beamformer by distance weight in UAE.Since the fairness will be more obvious if received power difference among MUs is smaller, we can obtain the standard deviation of received power among all MUs.For simplicity, assuming that there are only two MUs, we can compute the results difference between two MUs rather than the standard deviation of them.Figure 4 shows the received power difference between two MUs versus the distance between  2 and the virtual eBS.It is easy to see that both curves achieve the minimum when the distance is around 10 m and increase as the distance increases or decreases from this value.It is because of the fact that the distance between  1 and BS is fixed to 10 m in this scenario; the factors influencing both channels would be much similar when the distance of  2 fluctuates around this value.We can also see that our algorithm UFPE shows a lower difference compared with UAE, indicating that fairness of our model is more highlighted.

Sensitivity.
Since the sensitivity of circuit of receiver exits, MUs may not always harvest the energy they received.The energy harvesting probability of two MUs versus different circuit sensitivity is shown in Figure 5.It is easy to see that the energy harvesting probabilities of MU1 and MU2 decrease with the circuit sensitivity since it will become more difficult for the receiver to harvest energy when the circuit sensitivity is much too high.What is more,  2 shows a higher probability than  1 at the same  since the distance between  2 and virtual eBS is shorter, causing more power to be received than  1 .Since  1 's channel is worse, its EH probability decreases rapidly when  is not too big.

Throughput Performance.
On the WIT phase, we aim to maximize the minimum throughput among MUs by obtaining an optimal power and time allocation.This subsection can be divided into two parts, one is for power allocation and the other is for time allocation.
(a) Power Allocation.In this part, we compare our algorithm UFPE, in which the power allocation is using fraction  of energy, with UAE, in which the power allocation is using up all the battery energy.The max.-min.throughput among MUs versus circuit sensitivity can be seen in Figure 6, which indicates that the max.-min.throughput of both algorithms decreases as the sensitivity grows.Moreover, UFPE shows a higher throughput than UAE when  is larger than 0.02 W, and the difference will enlarge as  grows.The throughput of UFPE is about 3.3 times of that of UAE when  achieves 0.09 W. This means that our algorithm performs better in the practical scenario where circuit sensitivity exits.
Figure 7 depicts the max.-min.throughput among MUs versus  max .It can be seen that the max.-min.throughput of both algorithms increases as  max grows.It is because of the fact that as  max grows, the constraints of battery would be more relaxed, which can help MUs allocate a more adaptive slot to transmit information.On the other hand, both curves increase more smoothly when  max is higher, since the energy harvested at MUs will be more difficult to make the battery fully charged, causing the battery constraint to work no more.What is more, UFPE also shows a better performance than UAE with each specific  max .(b) Slot Allocation.An optimal slot allocation is also necessary to obtain better throughput performance.Figure 8 shows the comparison of our algorithm (dynamically allocating  in each slot) and other models, where the slot  is fixed to several values during every period.Generally, we choose five values with  = 0.1, 0.3, 0.5, 0.7, 0.9, respectively.It can be seen that the algorithm of dynamic allocating  shows better performance than any other algorithm with fixed .The model of  = 0.5 performs better than any other model of fixed .On the contrary, the model of  = 0.1 shows the worst performance.

Conclusion
In this paper, we have studied a DEIN model based on virtualization, with a multiantenna virtual eBS, a multiantenna dBS, and  single-antenna MUs with finite battery capacity.We proposed a fair resource allocation algorithm by joint optimization of the DL-UL time allocation, DL energy beamforming, and UL transmit power allocation with ZF based receive beamforming.Considering fairness, we firstly design the DL energy beamforming to maximize the minimum received power.After calculating the probability of energy being transmitted from the virtual eBS to MUs, we propose a simple power allocation with a fraction  of available energy allocated for UL WIT and adaptive time allocation in each slot.The simulation results have shown that the proposed UFPE algorithm achieves good performance.

Figure 1 :
Figure 1: A typical architecture of a DEIN based on virtualization.

Figure 4 :Figure 5 :
Figure 4: Received power difference between two MUs versus the distance between  2 and virtual energy BS.

Figure 6 :
Figure 6: Maximum-minimum throughput between two MUs versus circuit sensitivity.

)
3.3.Problem Formulation.Considering fairness, we optimize maximum-minimum (max.-min.)average UL WIT throughput of all MUs by jointly optimizing the time allocation   , the DL energy beams   , and the UL transmit power allocation   , .Let  ZF ( , ) denote the UL data rate of   in slot  as a function of the allocated energy  , for UL transmission in slot .Notice that  , is a feasible energy allocation policy, which satisfies 0 ≤  , ≤  max ,  +1, = min ( , +  +1, −  , ,  max ) ,  , =  (, {  }