Physical Layer Security and Optimal Multi-Time-Slot Power Allocation of SWIPT System Powered by Hybrid Energy

In this paper, a new approach is proposed to solve the constrained optimization problem of saving grid energy and increasing safety in a simultaneous wireless information and power transfer (SWIPT) system. The traditional grid energy is combined with the renewable energy to form a hybrid energy, which provides power for the system to achieve green wireless transmission. The transfer process of SWIPT system is divided into multiple time slots. The renewable energy is harvested and stored in battery at each time slot. A multi-time-slot artificial noise-assisted transmission strategy is proposed to reduce the signal to noise ratio (SNR) of eavesdropping link. A power allocation algorithm based on multi-time-slot golden section is given, which performs one-dimensional search on the power ratio of artificial noise to determine the transmit power of source node. And then the allocation algorithm is utilized to dynamically configure the harvested renewable energy for each time slot. When the battery capacity is constant, the maximum renewable energy is being used to reduce the grid power consumption. Finally, the performances of proposed schemes are evaluated by simulations in terms of various tradeoffs.


Introduction
Exploiting new energy and transmission mode for wireless communication to reduce the grid power consumption and improve the physical layer security (PLS) is the main issue faced by the next generation (5G) of mobile communication [1,2], and is considered as a problem of wide concern.Based on the above analysis, we consider a simultaneous wireless information and power transfer (SWIPT) system powered by the hybrid energy to save grid resource, and utilize artificial noise (AN) to interference eavesdropper.
The SWIPT technology received extensive attention with separation of destination nodes and carrying energy transmission [3][4][5].But the energy receiver (ER) is easy to eavesdrop the information transmitted to the information receiver (IR) from source (S) in the SWIPT system [6,7].And then ER is seen as a potential eavesdropper that will affect physical layer security.Besides, according to the lower energy conversion efficiency at ER, the high energy consumption is caused in the SWIPT system.In terms of the security issue, the authors in [8] designed a secure beamforming scheme in a SWIPT system, which can minimize the whole transmit power while maintaining the secrecy rate and energy harvesting at ER.Meanwhile, one of the technologies for enhancing the secrecy performance is to transmit artificial noise by jammers to interfere eavesdroppers [9].The authors in [10,11] recently considered the PLS with the AN by treating energy harvesting (EH) receivers as potential eavesdroppers.A secure MISO SWIPT system was proposed in [12], where an artificial noise is embedded to the information bearing signal to interfere the eavesdroppers and to harvest power by all receivers.However, these security schemes are only suitable for single time slot.Besides, in order to save power consumption, the paper [13] proposed the use of green renewable energy to mobile receivers in a multiple-antenna SWIPT system.The receivers obtain the renewable energy, which cannot be controlled and are not always available [14,15].The authors in [16] proposed a using renewable energy scheme which divides time slots to allocate the dynamic renewable energy, and gave a power allocation algorithm based on two-stage filling water.In the literature [17][18][19], using the directional water-filling algorithm in different channel models is considered to allocate the harvesting energy, including the Gaussian fading channel [17], the broadcast channel [18], and the multiple-access channel [19] to maximize the throughput in finite time horizon.However, these algorithms cannot guarantee optimal resource allocation for multi-time-slot.Although many studies have investigated the performance in SWIPT, relatively less work has been done on multi-time-slot wireless transmission.
In this paper, a SWIPT model based on the renewable energy and multi-slot artificial noise-assisted transmission is proposed.According to the dynamic arrival characteristics of renewable energy, the power splitting ratio and the transmitted power are determined at source node.A part of the transmitted power is used to transmit the signal to IR, and the source served as a jammer utilizes the rest power to transmit AN to ER.Meanwhile, ER utilizes the AN to harvest energy to satisfy collection demand.Besides, we propose a power allocation algorithm which can fully utilize renewable energy and guarantee secure transmission under the minimum grid consumption.

The System Model
We consider an energy cooperative network model as shown in Figure 1, which consists of a source, an information receiver, an energy receiver and a battery with a storage limit of E max at the source.The entire transmission process of the system is divided into N time slots and the size of each slot is t f .All channels are assumed to undergo blocking fading.f i , i ∈ 1, . . ., N represents the i-th slot channel gain of the S-IR link and g i , i ∈ 1, . . ., N represents the i-th slot channel gain of the S-ER link.The transmit power splitting ratio is defined as α i , i ∈ 1, . . . 1, N and the transmit power at S is defined as p i , i ∈ 1, . . ., N.
power by all receivers.However, these security schemes are only suitable for single time slot.Besides, in order to save power consumption, the paper [13] proposed the use of green renewable energy to mobile receivers in a multiple-antenna SWIPT system.The receivers obtain the renewable energy, which cannot be controlled and are not always available [14,15].The authors in [16] proposed a using renewable energy scheme which divides time slots to allocate the dynamic renewable energy, and gave a power allocation algorithm based on two-stage filling water.In the literature [17][18][19], using the directional water-filling algorithm in different channel models is considered to allocate the harvesting energy, including the Gaussian fading channel [17], the broadcast channel [18], and the multiple-access channel [19] to maximize the throughput in finite time horizon.However, these algorithms cannot guarantee optimal resource allocation for multi-time-slot.Although many studies have investigated the performance in SWIPT, relatively less work has been done on multi-time-slot wireless transmission.
In this paper, a SWIPT model based on the renewable energy and multi-slot artificial noise-assisted transmission is proposed.According to the dynamic arrival characteristics of renewable energy, the power splitting ratio and the transmitted power are determined at source node.A part of the transmitted power is used to transmit the signal to IR, and the source served as a jammer utilizes the rest power to transmit AN to ER.Meanwhile, ER utilizes the AN to harvest energy to satisfy collection demand.Besides, we propose a power allocation algorithm which can fully utilize renewable energy and guarantee secure transmission under the minimum grid consumption.

The System Model
We consider an energy cooperative network model as shown in Figure 1, which consists of a source, an information receiver, an energy receiver and a battery with a storage limit of max E at the source.The entire transmission process of the system is divided into N time slots and the size of each slot is f t .All channels are assumed to undergo blocking fading., 1,...,

The Mathematical Modeling
In this considered network, assuming that source is powered by hybrid energy consisted of the traditional power grid and the renewable energy.Therefore, the mixed transmit signal at source should be where s 0,i , s 1,i respectively represent the i-th information signal and AN signal at S, they are independent circular symmetric complex Gaussian random variables with zero mean and unit variance.
The received signal expression of IR and ER are formulated as where n ER,i and n IR,i respectively represent the additive white Gaussian noise (AWGN) at ER and IR with n ER,i ∼ CN(0, σ 2 2,i ) and n IR,i ∼ CN(0, σ 2 1,i ).Then the SNR at IR and ER are respectively given by According to Equations ( 4) and ( 5), the secrecy rate of i-th slot can be calculated as where [x] + = max(0, x).In this scheme, we assume that the system is based on timeslot changing.We introduce the following indicator function for the event of outage with respect to the minimum secrecy rate r 0 at each time slot: where the outage probability can be expressed as We consider a harvesting renewable energy plan at source as shown Figure 2, and assume that source is equipped a battery with initial power E 0 .When a time slot has just begun, the battery stores the harvested renewable energy, i.e., E i .During the wireless transmission process, the transmit power consists of two parts, one part of which is provided by the battery, i.e., p r i and the other obtains from grid, i.e., p g i .
The renewable energy of battery.
After all time slots, the collected energy of ER should satisfy the average collection requirement, i.e., where ξ ∈ is the energy collection efficiency depending on the rectification and energy harvesting circuitry, and ave p represents the need of average collection power.Besides, the collected renewable energy at each time slot must be stored in the battery before it is used, and the causal relationship is shown as Due to the limited battery capacity, the harvesting energy of source in each time slot cannot exceed the maximum limit of the battery, i.e., where max is the average harvesting energy of battery per second.And we define the maximum average energy consumption threshold of the grid in order to limit the actual grid consumption.Therefore, the constraint relation is expressed as where g ave p represents the average grid consumption.Besides, source is a power-constrained node and the maximum power allowed for transmission is peak p .Therefore, the transmit power of each slot must satisfy the following formula We aim at minimizing the secrecy outage probability for the IR by jointly optimizing i α and i p .Based on the definition in Equation ( 7), the optimization problem can be formulated as After all time slots, the collected energy of ER should satisfy the average collection requirement, i.e., 1 where ξ ∈ [0, 1] is the energy collection efficiency depending on the rectification and energy harvesting circuitry, and p ave represents the need of average collection power.Besides, the collected renewable energy at each time slot must be stored in the battery before it is used, and the causal relationship is shown as Due to the limited battery capacity, the harvesting energy of source in each time slot cannot exceed the maximum limit of the battery, i.e., where E max = MP r,ave t f with P r,ave = 20 dBm.M represents the battery storage capacity coefficient and P r,ave is the average harvesting energy of battery per second.And we define the maximum average energy consumption threshold of the grid in order to limit the actual grid consumption.Therefore, the constraint relation is expressed as where p g ave represents the average grid consumption.Besides, source is a power-constrained node and the maximum power allowed for transmission is p peak .Therefore, the transmit power of each slot must satisfy the following formula We aim at minimizing the secrecy outage probability for the IR by jointly optimizing α i and p i .Based on the definition in Equation ( 7), the optimization problem can be formulated as min (8), ( 9), ( 10), ( 11) and (12)

The Optimal Power Splitting Ratio and Transmit Power
The above mathematical model involves both α i and p i that need to be optimized.We utilize the two-stage optimization strategy in order to solve the above problem [20].At the first stage, one-dimensional search is performed for the power splitting ratio of each slot based on the golden section principle, and then the minimum required power is obtained.Secondly, we apply the Lagrange duality method to solve Equation ( 13) based on the ratio solution of the first stage.

The Power Splitting Ratio
Because of the transmission power of S is constrained, minimizing p i is also one of the goals [21].Given any 0 ≤ α i ≤ 1, let p(α i ) denotes the minimum required power to maintain the secrecy rate r 0 allowed for the system, i.e., R i (α i , p i ) = r 0 .It can be shown as By substituting ( 4) and ( 5) into (14), and then p i and α i should satisfy the condition One can show from (15) that it is a quadratic function of p i .Therefore, the solution of p i can be expressed as follows where α i is the independent variable, p(α i ) is the dependent variable.α * i is defined as the optimal power splitting ratio.∆ is given by One-dimensional search is performed for (15), and obtain α i which satisfies minimum p(α i ).When N is equal to 100, Figure 3 shows the optimal distribution of α i in i-th time slot with σ 2  1,i = σ 2 2,i = 0.01 and r 0 = 10 kbit/s.We see that the optimal distribution is in the range of 0.455 to 0.505, which changes according to the channel state.

The Optimal Transmit Power
In the first stage, we obtain the appropriate i α .In the second stage, the optimization problem is convex about i p .The Lagrangian function of ( 13) is expressed as where , , i i u λ ϕ and η are the dual variables associated with the average collected energy ave p at ER, the harvesting renewable energy i E at source, the battery capacity max E and the average grid consumption g ave p , respectively.Then the partial Lagrange dual function of ( 18) is expressed as , , , ) The dual problem of ( 13) is given by max ( , , , ) . .0, 0, , 0.
Substituting the optimal power splitting ratio from the first stage into Equation (18), we obtain the expression about the optimal transmit power of per time slot

The Optimal Transmit Power
In the first stage, we obtain the appropriate α i .In the second stage, the optimization problem is convex about p i .The Lagrangian function of ( 13) is expressed as where λ, u i , ϕ i and η are the dual variables associated with the average collected energy p ave at ER, the harvesting renewable energy E i at source, the battery capacity E max and the average grid consumption p g ave , respectively.Then the partial Lagrange dual function of ( 18) is expressed as The dual problem of ( 13) is given by Substituting the optimal power splitting ratio from the first stage into Equation (18), we obtain the expression about the optimal transmit power of per time slot where c = (η 21) is a cubic equation and is complex, we utilize MATLAB to solve the cubic equation and obtain the optimal power.

A Power Allocation Algorithm
In this paper, we utilize the golden section principle to search the power splitting ratio, and then obtain the optimal transmit power at source based on the Lagrange duality method.Based on the above analysis, a power allocation algorithm structure is proposed.The steps of the algorithm are as follows: Step1: Initialize the optimization problem and algorithm parameters.We want to minimize the outage probability which includes the variable α i and p i .The renewable energy power harvested in the battery of each slot is p r i .Step2: Perform one-dimensional searching on α i and obtain the transmit power at source.Step3: The power allocation starts from the i-th slot with i = N, and the harvesting renewable energy cannot exceed E max .If p r i > E max , set p r i = E max .Step4: In this step, we compare the renewable energy p r i and the optimal power p i in i-th slot.When p i is greater then p r i , an allocation scheme is proposed referring to the model (a) in Figure 4.As shown in the model (a), all of the renewable energy of the i-th time slot is used, the rest of the transmit power is filled by grid, i.e., p i = p g i + p r i .And then set i = i − 1.In contrast, utilize the model (b) to allocate the hybrid power.In the model (b), find the j-th slot with j < i satisfying with p m > 0, and then set i = j.During this process, no grid power is used.
Step6: If there is still remaining renewable energy at rest of the time slot and p i < p peak , the remaining energy will be directly supplied to the source node of the corresponding time slot to transmit information for IR.
In this algorithm, we build an inverse sequence for the time slot to guarantee all the renewable energy be utilized.) ) 21) is a cubic equation and is complex, we utilize MATLAB to solve the cubic equation and obtain the optimal power.

A Power Allocation Algorithm
In this paper, we utilize the golden section principle to search the power splitting ratio, and then obtain the optimal transmit power at source based on the Lagrange duality method.Based on the above analysis, a power allocation algorithm structure is proposed.The steps of the algorithm are as follows: Step1: Initialize the optimization problem and algorithm parameters.We want to minimize the outage probability which includes the variable i α and i p .The renewable energy power harvested in the battery of each slot is Step2: Perform one-dimensional searching on i α and obtain the transmit power at source.

Step3:
The power allocation starts from the i-th slot with i = N, and the harvesting renewable Step4: In this step, we compare the renewable energy r i p and the optimal power i p in i-th slot.
When i p is greater then r i p , an allocation scheme is proposed referring to the model (a) in Figure 4.
As shown in the model (a), all of the renewable energy of the i-th time slot is used, the rest of the transmit power is filled by grid, i.e.,


, and then set i = j.During this process, no grid power is used.
Step6: If there is still remaining renewable energy at rest of the time slot and i p e a k p p < , the remaining energy will be directly supplied to the source node of the corresponding time slot to transmit information for IR.
In this algorithm, we build an inverse sequence for the time slot to guarantee all the renewable energy be utilized.

Simulation Results
In this section, we provide simulation results to validate our theoretical analysis and evaluate the performance of the proposed scheme.We run N = 100 frames in each simulation and set the maximum transmit power as 30 dBm, i.e., p peak = 30 dBm.The harvested energy of ER follows non-negative uniform distribution with mean p ave = 20 dBm.The energy conversion efficiency is set as ξ = 0.8 [22].Besides, in order to make fully use of the renewable energy, we set a little average grid consumption power as p g ave = 20 dBm. Figure 5 shows that the relationship between the power ratio, the outage probability and the transmit power.We define the power splitting ratio vector α = {α 1 , α 2 , . . . ,α N }.Three vector schemes are given: all the elements of vector are 0.1, all the elements of vector are 0.3 and all the elements of vector are the optimal search values.We can see that the transmit power gradually decreases and tends to be stable with the increasing outage probability.Besides, when the outage probability is fixed, the optimal α corresponds to the whole transmit power is minimum.When the difference between the power splitting ratio and the optimal α is larger, the difference of transmit power value will be larger.Therefore, in the case of satisfying a certain outage probability, utilizing the golden section to search the optimal α can effectively reduce the transmit power at source.non-negative uniform distribution with mean 20 dBm ave p = .The energy conversion efficiency is set as = 0.8 ξ [22].Besides, in order to make fully use of the renewable energy, we set a little average grid consumption power as 20 dBm g ave p = .
Figure 5 shows that the relationship between the power ratio, the outage probability and the transmit power.We define the power splitting ratio vector . Three vector schemes are given: all the elements of vector are 0.1, all the elements of vector are 0.3 and all the elements of vector are the optimal search values.We can see that the transmit power gradually decreases and tends to be stable with the increasing outage probability.Besides, when the outage probability is fixed, the optimal α corresponds to the whole transmit power is minimum.When the difference between the power splitting ratio and the optimal α is larger, the difference of transmit power value will be larger.Therefore, in the case of satisfying a certain outage probability, utilizing the golden section to search the optimal α can effectively reduce the transmit power at source.From Figure 6, we can see that the grid power consumption is gradually decreasing when the battery capacity increases.The power allocation algorithm which we proposed is better than the directly filling water algorithm and the two-stage filling water algorithm on reducing the consumption of grid.Therefore, fully utilizing renewable energy can improve the capacity, which reduces unnecessary power consumption of the grid as much as possible.
Figure 7 shows that the curve changing of the secrecy rate with the influence of M and i α .
According to the three-dimensional figure, we obtain that the secrecy rate gradually increases and tends to be stable with changing of M, and the curve of the secrecy rate has a peak at around 0.5 i α = . Therefore, choosing the appropriate battery capacity and the power splitting ratio can effectively improve the secrecy performance.From Figure 6, we can see that the grid power consumption is gradually decreasing when the battery capacity increases.The power allocation algorithm which we proposed is better than the directly filling water algorithm and the two-stage filling water algorithm on reducing the consumption of grid.Therefore, fully utilizing renewable energy can improve the capacity, which reduces unnecessary power consumption of the grid as much as possible.
Figure 7 shows that the curve changing of the secrecy rate with the influence of M and α i .According to the three-dimensional figure, we obtain that the secrecy rate gradually increases and tends to be stable with changing of M, and the curve of the secrecy rate has a peak at around α i = 0.5.Therefore, choosing the appropriate battery capacity and the power splitting ratio can effectively improve the secrecy performance.

Conclusions and Future Works
This paper considers the problems of safety and energy in a SWIPT system, in presence of potential eavesdropper, i.e., ER.The paper shows that secure transmission can be achieved by adding AN based on multi-slot.And we have optimized the power splitting ratio and the transmit power at source based on the time domain.Also, we propose a power allocation algorithm for the system to allocate the renewable energy at each time slot to guarantee the minimum grid power consumption.
The simulation examples have validated the proposed methods and confirmed the effect on the secrecy performance and grid energy.A system with a full duplex node or IR served as a jammer may be interesting future topic.

Conclusions and Future Works
This paper considers the problems of safety and energy in a SWIPT system, in presence of potential eavesdropper, i.e., ER.The paper shows that secure transmission can be achieved by adding AN based on multi-slot.And we have optimized the power splitting ratio and the transmit power at source based on the time domain.Also, we propose a power allocation algorithm for the system to allocate the renewable energy at each time slot to guarantee the minimum grid power consumption.
The simulation examples have validated the proposed methods and confirmed the effect on the secrecy performance and grid energy.A system with a full duplex node or IR served as a jammer may be interesting future topic.

Conclusions and Future Works
This paper considers the problems of safety and energy in a SWIPT system, in presence of potential eavesdropper, i.e., ER.The paper shows that secure transmission can be achieved by adding AN based on multi-slot.And we have optimized the power splitting ratio and the transmit power at source based on the time domain.Also, we propose a power allocation algorithm for the system to allocate the renewable energy at each time slot to guarantee the minimum grid power consumption.
The simulation examples have validated the proposed methods and confirmed the effect on the secrecy performance and grid energy.A system with a full duplex node or IR served as a jammer may be interesting future topic.
i-th slot channel gain of the S-ER link.The transmit power splitting ratio is defined as , 1,...1,

Figure 1 .
Figure 1.The system model.Figure 1.The system model.

Figure 1 .
Figure 1.The system model.Figure 1.The system model.

,
. M represents the battery storage capacity coefficient and , r ave P

Figure 2 .
Figure 2. The renewable energy of battery.
then set i = i − 1.In contrast, utilize the model (b) to allocate the hybrid power.In the model (b), find the j-th slot with j < i satisfying with

Figure 4 .
Figure 4. Two types of power allocation model.(a) The less renewable energy allocation model; (b) The enough renewable energy allocation model.

Figure 4 .
Figure 4. Two types of power allocation model.(a) The less renewable energy allocation model; (b) The enough renewable energy allocation model.

Figure 5 .
Figure 5.The average transmit power with different i α , M = 3,

Figure 6 .
Figure 6.The grid power consumption with different power allocation algorithm.

Figure 7 .
Figure 7.The secrecy rate with the influence of M and i α .
M the grid energy consumption (mW)the directly filling water algorithm the two-stage filling water algorithm the power allocation algorithm proposed in this paper

Figure 7 .
Figure 7.The secrecy rate with the influence of M and i α .
M the grid energy consumption (mW)the directly filling water algorithm the two-stage filling water algorithm the power allocation algorithm proposed in this paper 0