Next Article in Journal
Digital Twins in Solar Farms: An Approach through Time Series and Deep Learning
Previous Article in Journal
Accelerating In-Transit Co-Processing for Scientific Simulations Using Region-Based Data-Driven Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy-Efficient Power Allocation in Non-Linear Energy Harvesting Multiple Relay Systems

1
Jiangsu Key Laboratory of Wireless Communications, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
2
Engineering Research Center of Health Service System Based on Ubiquitous Wireless Networks, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
*
Author to whom correspondence should be addressed.
Submission received: 18 April 2021 / Revised: 7 May 2021 / Accepted: 12 May 2021 / Published: 17 May 2021

Abstract

:
In this paper, to maximize the energy efficiency (EE) in the two-hop multi-relay cooperative decoding and forwarding (DF) system for simultaneous wireless information and power transmission (SWIPT), an optimal power allocation algorithm is proposed, in which the relay energy harvesting (EH) adopts a nonlinear model. Under the constraints, including energy causality, the minimum transmission quality of information and the total transmission power at the relays, an optimization problem is constructed to jointly optimize the transmit power and power-splitting (PS) ratios of multiple relays. Although this problem is a nonlinear fractional programming problem, an iterative algorithm is developed to obtain the optimal power allocation. In particular, the joint power allocation at multiple relays is first decoupled into a single relay power allocation, and then single-relay power allocation is performed by the Dinkelbach iteration algorithm, which can be proven that it is a convex programming problem. Its closed form solutions for different polylines of EH models are obtained by using mathematical methods, such as monotonicity, Lagrange multipliers, the KKT condition and the Cardan formula. The simulation results show the superiority of the power allocation algorithm proposed in this paper in terms of EE.

1. Introduction

According to statistics, the average annual energy consumption of the information and communications industry accounts for 3% of the global total energy consumption, and carbon emissions exceed 2% of the global total [1,2]. Collaborative relay technology and energy harvesting technology have emerged to meet the quality of communication while reducing energy consumption. In the post-5G era, because wireless networks are becoming more and more complex, how to design appropriate energy allocation strategies to improve system performance is a problem that needs further research [3].
Collaborative relay technology based on the energy harvesting (EH) model has been studied. In [4], the throughput maximization problem was solved via a successive convex approximation approach, in which the relay harvests energy in the natural environment. In [5], the authors analyzed the difference between an EH relay and a normal collaborative relay. The authors in [6] presented a power allocation algorithm in which the relay could switch between amplify and forward (AF) and decoding and forwarding (DF) modes to obtain the optimal transmit power and operational mode of the relay in order to maximize the energy efficiency (EE). The power allocation problem for a two-hop two-way single-relay system to maximize the EE was studied in [7]. Both source nodes and relay nodes use EH technology, and it was the aim of the authors to maximize the system capacity in [8]. The authors of [9] mainly analyzed and compared the advantages and disadvantages of the two schemes of simultaneous wireless information and power transmission (SWIPT). However, the above EH model assumes a linear EH model, in which the power conversion efficiency from the radio frequency (RF) to the direct current (DC) is a fixed constant. In fact, because the diodes, inductors and capacitors in the actual energy harvester [10] are nonlinear, the power conversion efficiency from the RF to the DC will vary with the input power of the energy harvester circuit.
There are also preliminary studies on the application of nonlinear EH models in relay systems. In [11], a power allocation algorithm for a nonlinear EH model that maximized the total harvesting power of the energy harvester was presented. The problem of maximizing the EE of a two-way DF single-relay network was studied in [12]. The simulation results in [13] show that the power allocation based on the nonlinear EH model has higher throughput than that based on the existing linear model. In [14], the authors investigated the interrupt performance of nonlinear EH relay networks with a PS scheme. However, the above works only studied power allocation for single-relay systems and not for multi-relay systems.
This paper presents a power allocation algorithm based on the nonlinear EH model for maximizing the EE in a two-hop multi-relay cooperative forwarding system. For clarity, we have listed the contributions as follows:
  • In this paper, a piecewise linear EH model is used to approximate the nonlinear EH model. Based on that, the EE maximization problem with quality of service (QoS) guaranteed is a complex, nondeterministic polynomial (NP) problem and can be decoupled into a single-relay power allocation problem. It is optimized based on the Dinkelbach iteration algorithm, and the optimization problem after iteration simplification is proven to be a convex programming problem.
  • Because the optimization problem is a maximum–minimum function problem, it is decomposed into two maximization problems to obtain the optimal solution by adding the corresponding constraints. For each maximization problem, simplified expressions are obtained according to the different polylines of EH models, whose closed-form solutions of the optimal relay transmit power and optimal power-splitting (PS) ratio are obtained by using mathematical methods such as monotonicity, Lagrange multipliers, the KKT condition and the Cardan formula.
  • The simulation results show that the piecewise linear EH model performs better than the traditional linear EH model, and the EE of the multi-relay system is better than that of the single-relay system.
The rest of this paper is organized as follows. In Section 2, the multi-relay cooperative DF system based on nonlinear EH is described, and the optimization problem is formulated to maximize the EE. In Section 3, the joint optimization algorithm of the relay transmit power and PS ratio is proposed. In Section 4, the algorithm is simulated, and its performance is analyzed. Finally, this paper is summarized in Section 5.

2. System Model and Optimization Problem

2.1. System Model

As shown in Figure 1, we considered the two-hop multi-relay cooperative forwarding system, including a source node S , a destination node D and N relay nodes R i for i = 1 , 2 , , N . All the R i can harvest energy, while S and D cannot harvest energy. Assuming that there is no direct link between S and D , the corresponding channel of each relay is independent of each other, where h i represents the channel gain from S to R i and g i represents the channel gain from R i to D .
Compared with PS, time switching (TS) cannot obtain a better trade-off between the transmission rate and the RF energy, and it was proven that the PS scheme achieves higher throughput than the TS scheme in [11]. Therefore, we assumed that R i could only harvest energy from the signal sent by S , and each transmission block period was divided into two time slots in the PS scheme of [13]. In the first time slot, R i receives the signal sent by S and divides the received signal power into two parts, according to the ratio of 1 ρ i : ρ i . One part enters the information processing module to decode the information, expressed as
y R i I D ( k ) = 1 ρ i ( P S h i s ( k ) + n i a ( k ) ) + n i I D ( k )
where E [ | s ( k ) | 2 ] = 1 ; P S is the transmit power of S , ρ i is the PS ratio of R i , n i a ( k ) C N ( 0 , σ a i 2 ) is the additive white Gaussian noise of the relay receiving antenna and n i I D ( k ) C N ( 0 , σ I D i 2 ) is the additive baseband Gaussian noise. Therefore, the SNR of the receiver of R i can be expressed as γ r i = ( 1 ρ i ) P S h i ( 1 ρ i ) σ a i 2 + σ I D i 2 .
The other part enters the EH module to provide energy for the information decoding of this time slot and the information forwarding of the next time slot. The received RF signal at R i used for energy harvesting is
y R i E H ( k ) = ρ i ( P S h i s ( k ) + n i a ( k ) )
Because of P S h i σ a i 2 , the power used by R i for EH is P R F i ( ρ i ) = ρ i P S h i .
In many cases, because of the nonlinearity of the EH model, the problems studied are often too complex to be solved, so the piecewise linear model can not only keep the nonlinearity of the original model to a maximum extent but also solve the problem well [15]. For EH, we used a piecewise linear model, which can be expressed as
P H i ( ρ i ) = { 0 , P R F i ( ρ i ) < P t h 1 a l P R F i ( ρ i ) + b l , P R F i ( ρ i ) [ P t h l , P t h l + 1 ] , l = 1 , 2 , , L 1 P m , P R F i ( ρ i ) > P t h L = a l ρ i P S h i + b l , P R F i ( ρ i ) [ P t h l , P t h l + 1 ] , l = 0 , 1 , 2 , , L
where P R F i is the receiving power of the EH circuit, P t h = { P t h l | 0 l L + 1 } is the dividing point separating P R F i into L + 1 segments, satisfying P t h 0 = 0 , P t h L + 1 = + , a l and b l are the slope and intercept of the linear function in the l t h ( 0 l L ) segment, respectively, satisfying a 0 = b 0 = a L = 0 , b L = P m , and P m denotes the maximum harvesting power.
When it has enough segments, it can be approximated to the nonlinear EH model, but this will also increase its complexity. Therefore, the real measurement data [16] of the energy harvester can be fitted by the piecewise linear function to replace the nonlinear EH model [10]. Figure 2 compares the harvested power in different EH models.
Because the transmission performance of DF is better than AF in most instances, the relay uses DF in the second time slot, which is to say that R i uses the harvested energy to forward the decoded information to D . Assuming that R i sends the signal x i ( k ) to D , the received RF signal at D is
y D i ( k ) = P i g i x i ( k ) + n D ( k )
where E [ | x i ( k ) | 2 ] = 1 , P i is the transmit power of R i and n D ( k ) C N ( 0 , σ D 2 ) is the additive white Gaussian noise caused by the antenna reception and radio frequency to baseband conversion at the destination node. Therefore, the signal-to-noise ratio (SNR) of the receiver of D can be expressed as γ d i = P i g i σ D 2 .
According to the maximum merge ratio criterion [17], the total SNR of all the links participating in the communication is γ = i = 1 N γ i = i = 1 N min ( γ r i , γ d i ) . Therefore, by using the Shannon formula, the total information transmission rate of S R D is
R S D = 1 2 log 2 ( 1 + γ ) = 1 2 log 2 ( 1 + i = 1 N γ i )
In addition, the total power consumption during the transmission of S R D is
P = P S + 2 P C + i = 1 N ( P i + P Q P H i ( ρ i ) )
where P C is the consumption of power of the S and D decoding circuits and P Q is the consumption of power of the R i circuit.

2.2. Optimization Problem

EE is defined as the ratio of the total information transmission rate to the total power consumption for all links, expressed as
η = R S D P = 1 2 log 2 ( 1 + i = 1 N γ i ) P S + 2 P C + i = 1 N ( P i + P Q P H i ( ρ i ) )
This paper studies how to jointly optimize the relay transmit power and the PS ratio to maximize the EE. The optimization problem can be expressed as
( P 1 ) max P i , ρ i η s . t . C 1 : 0 ρ i 1 C 2 : 0 P i + P Q P H i ( ρ i ) C 3 : γ i γ min C 4 : P t h l ρ i P S h i P t h l + 1 , l = 0 , 1 , 2 , , L C 5 : i = 1 N P i P max
where γ min = 2 R min 1 and R min is the minimum required rate for S R D , P max is the threshold of the relay nodes’ total transmit power, constraint C 2 means the energy causality, C 3 is the constraint of the quality of information transmission, C 4 is the constraint that the energy harvester works in the l t h linear regions and C 5 is the constraint of the total transmit power of all relays, and it is a requirement to consider interference across the network [18].

3. Joint Optimization Algorithm for the Relay Transmit Power and PS Ratio

The optimization problem (P1) is a multi-variable joint optimization problem and a complex NP problem. Because of the piecewise linear EH model, the optimization problem has different results for its different polyline segments.
In order to solve (P1), there are two main steps as follows. In the first step, the power allocation problem of multiple relays is decoupled into a single-relay power allocation problem by using the fixed-point iteration algorithm, which assumes that only the power of this relay is to be allocated and that the power allocation of other relays is known. For convenience, some iteration variables ψ i = w = 1 , w i N γ w , φ i = w = 1 , w i N P w and ϕ i = w = 1 , w i N ( P w + P Q P H w ( ρ w ) ) are introduced. In the second step, we assume ρ i = 1 and calculate the maximum number of segments which P R F i ( ρ i ) may belong to, expressed as S i ( s i { 0 , 1 , , S i } , i = 1 , 2 , , N ) . In a given number of segments s i , we transform this non-convex problem into a convex problem by subtraction and find a closed-form solution.
According to the Dinkelbach method [19], the maximum EE x * = max P i , ρ i η = 1 2 log 2 ( 1 + ψ i + γ i ( P i * , ρ i * ) ) P S + 2 P C + ϕ i + P i * + P Q P H i ( ρ i * ) is achieved if and only if the following equation is satisfied:
max P i , ρ i 1 2 log 2 ( 1 + ψ i + γ i ( P i , ρ i ) ) x * ( P S + 2 P C + ϕ i + P i + P Q P H i ( ρ i ) ) = 0
Therefore, we convert the fractional objective function to subtraction for a given segment s i , and (P1) can be rewritten as
( P 2 ) max P i , ρ i 1 2 log 2 ( 1 + ψ i + γ i ( P i , ρ i ) ) x * ( P S + 2 P C + ϕ i + P i + P Q P H i ( ρ i ) ) s . t . C 1 C 5
However, (P2) is still a non-convex problem. In order to convert this problem into a convex one, we let t i = 1 ρ i and introduce an auxiliary variable r i , expressed as
r i = 1 2 log 2 ( 1 + i = 1 N γ i ) = min ( 1 2 log 2 ( 1 + ψ i + γ r i ) , 1 2 log 2 ( 1 + ψ i + γ d i ) )
The optimization problem (P2) can be rewritten as
( P 3 ) max P i , t i f ( P i , t i , r i ) = r i x * ( P S + 2 P C + ϕ i + P i + P Q a s i ( 1 t i ) P S h i b s i ) s . t . C 1 : t min t i t max C 2 : 0 P i a s i ( 1 t i ) P S h i + b s i P Q C 3 : r i 1 2 log 2 ( 1 + γ min ) C 4 : 1 2 log 2 ( 1 + ψ i + t i P S h i t i σ a i 2 + σ I D i 2 ) r i C 5 : 1 2 log 2 ( 1 + ψ i + P i g i σ D 2 ) r i C 6 : P i P max φ i
where t min = 1 P t h s i + 1 P S h i and t max = 1 P t h s i P S h i .
Theorem 1.
The optimization problem (P3) is a convex programming problem.
Proof of Theorem 1.
The objective function of (P3) regarding P i , t i and r i is a linear function, and constraints C 1 , C 2 , C 3 and C 6 are linear constraints, while constraints C 4 and C 5 are nonlinear constraints, so (P2) is a nonlinear programming problem. Because linear functions are both convex and concave, the objective function f ( P i , t i , r i ) is convex, and constraints C 1 , C 2 , C 3 and C 6 are concave constraints. By letting f 1 ( t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i t i σ a i 2 + σ I D i 2 ) and f 2 ( P i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ D 2 ) , then we have 2 f 1 ( t i ) t i 2 = P S h i σ I D i 2 ( 2 a t t i + b t ) ( a t t i 2 + b t t i + c t ) 2 2 ln ( 2 ) < 0 , 2 f 2 ( P i ) P i 2 = g i 2 ( P i g i + ( 1 + ψ i ) σ D 2 ) 2 2 ln ( 2 ) < 0 , where a t = ( P S h i + ( 1 + ψ i ) σ a i 2 ) σ a i 2 , b t = ( P S h i + 2 ( 1 + ψ i ) σ a i 2 ) σ I D i 2 and c t = ( 1 + ψ i ) σ I D i 4 . Therefore, f 1 ( t i ) is a concave function for t i , and f 2 ( P i ) is also a concave function for P i ( i.e., the constraints C 4 and C 5 are concave constraints). From the definition of convex programming, the optimization problem (P3) is a convex programming problem. □
As Equation (11) is composed of two formulas, (P3) is a maximum–minimum function problem and can be decomposed into two maximization problems to obtain the optimal solution by adding the corresponding constraints, for generality, we assume that all noise variances are equal (i.e., σ a i 2 = σ I D i 2 = σ D 2 = σ 2 ).
1. If γ r i γ d i (i.e., P i t i P S h i ( t i + 1 ) g i ), we have r i = 1 2 log 2 ( 1 + ψ i + γ d i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) . The optimization problem can be expressed as
( P 4 ) max P i , t i y 1 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P S + 2 P C + ϕ i + P i + P Q a s i ( 1 t i ) P S h i b s i ) s . t . C 1 : t min t i t max C 2 : 0 P i a s i ( 1 t i ) P S h i + b s i P Q C 3 : P i g i σ 2 γ min ψ i C 4 : P i t i P S h i ( t i + 1 ) g i C 5 : P i P max φ i
2. If γ r i γ d i (i.e., P i t i P S h i ( t i + 1 ) g i ), we have r i = 1 2 log 2 ( 1 + ψ i + γ r i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) . The optimization problem can be expressed as
( P 5 ) max P i , t i y 2 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) x * ( P S + 2 P C + ϕ i + P i + P Q a s i ( 1 t i ) P S h i b s i ) s . t . C 1 : t min t i t max C 2 : 0 P i a s i ( 1 t i ) P S h i + b s i P Q C 3 : t i P S h i ( t i + 1 ) σ 2 γ min ψ i C 4 : P i t i P S h i ( t i + 1 ) g i C 5 : P i P max φ i
Due to the uncertainty of P R F i , according to the EH model in Equation (3), (P4) and (P5) are divided into three cases: Case 1, P R F i ( t i ) < P t h 1 ; Case 2, P R F i ( t i ) > P t h L ; and Case 3, P R F i ( t i ) [ P t h l , P t h l + 1 ] , l = 1 , 2 , , L 1 .

3.1. Case 1

When P R F i ( t i ) < P t h 1 (i.e., s i = 0 ), we have P H i = 0 and P i = 0 , leading to a communication outage between R i and D . As a result, the EE is zero.

3.2. Case 2

When P R F i ( t i ) > P t h L (i.e., s i = L ), we have P H i ( t i ) = P m , t min = min ( 0 , 1 P t h L + 1 P S h i ) = 0 and t max = 1 P t h L P S h i .
Considering the constraints of (P4) and (P5), the optimization problem can be rewritten as three simplification scenarios:
  • Scenario I: When t i P S h i ( t i + 1 ) g i P 2 , we have t i t 2 , and (P4) can be rewritten as
    max P i , t i y 1 ( P i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q P m ) s . t . P 1 P i P 2 t min t i t max t i t 2
There is a solution to Equation (15) if and only if t 2 t max is satisfied. The objective function is independent of t i , so we usually let t i * = t max . Due to 2 y 1 P i 2 < 0 , y 1 ( P i ) is a convex function with respect to P i , and it can be solved by derivation. Let d y 1 d P i = 0 and have P 0 = 1 x * 2 ln 2 ( 1 + ψ i ) σ 2 g i . The solution I can be expressed as
{ P i I * = [ P 0 ] P 1 P 2 ρ i I * = 1 t i I * = 1 t max
where t 2 = P 2 g i P S h i P 2 g i , P 1 = max ( 0 , ( γ min ψ i ) σ 2 g i ) , P 2 = min ( P max φ i , P m P Q ) and [ x ] a b denotes max ( a , min ( x , b ) ) .
  • Scenario II: When t i P S h i ( t i + 1 ) g i P 2 , we have t i t 2 , and (P4) can be rewritten as
    max P i , t i y 1 ( P i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q P m ) s . t . P 1 P i t i P S h i ( t i + 1 ) g i t min t i t max t i t 2
There is a solution to Equation (17) if and only if t min t 2 is satisfied. The objective function is related to t i , so we use the Lagrange multiplier method to convert a constrained problem to an unconstrained problem max P i , t i L ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q P m ) + λ i ( t i P S h i ( t i + 1 ) g i P i ) , where λ i is the Lagrange multiplier. Since 2 L P i 2 < 0 and 2 L t i 2 < 0 , we use the KKT condition for the solution. Solution II can be expressed as
{ P i I I * = [ P 0 ] P 1 t i I I * P S h i ( t i I I * + 1 ) g i ρ i I I * = 1 t i I I * = 1 min ( t 2 , t max )
where t 2 = P 2 g i P S h i P 2 g i , P 0 = 1 x * 2 ln 2 ( 1 + ψ i ) σ 2 g i and P 1 = max ( 0 , ( γ min ψ i ) σ 2 g i ) .
  • Scenario III: Considering the constraints of (P5), the optimization problem can be rewritten as
    max P i , t i y 2 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) x * ( P i + P Q P m ) s . t . t i P S h i ( t i + 1 ) g i P i P 2 t min t i t max t 1 t i t 2
There is a solution to Equation (19) if and only if D = [ t 1 , t 2 ] [ t min , t max ] is satisfied and we let D = [ t a , t b ] . First, because Equation (19) is a linear programming problem for P i and y 2 P i < 0 , the optimal solution P i * = t i * P S h i ( t i * + 1 ) g i can be obtained according to monotonicity. Then, the binary objective function is transformed into a one-variable objective function y 2 ( t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) x * ( t i P S h i ( t i + 1 ) g i + P Q P m ) , which is only related to t i by the substitution of P i * . We let d y 2 d t i = 0 and have t 0 = g i x * 2 ln 2 ( 1 + ψ i ) σ 2 x * 2 ln 2 ( P S h i + ( 1 + ψ i ) σ 2 ) g i . Because d 2 y 2 d t i 2 can be positive or negative, we have t i * { t a , t b , [ t 0 ] t a t b } , and t i * is determined by the value of y 2 ( t i * ) . Solution III can be expressed as
{ P i I I I * = t i I I I * P S h i ( t i I I I * + 1 ) g i ρ i I I I * = 1 t i I I I * = 1 arg max t i * { t a , t b , [ t 0 ] t a t b } ( 1 2 log 2 ( 1 + ψ i + t i * P S h i ( t i * + 1 ) σ 2 ) x * t i * P S h i ( t i * + 1 ) g i )
where t 1 = ( γ min ψ i ) σ 2 P S h i ( γ min ψ i ) σ 2 , t 2 = P 2 g i P S h i P 2 g i and P 2 = min ( P max φ i , P m P Q ) .

3.3. Case 3

When P R F i ( t i ) [ P t h l , P t h l + 1 ] , l = 1 , 2 , , L 1 (i.e., s i = l ), we have P H i ( t i ) = a l · ( 1 t i ) P s h i + b l , t min = 1 P t h l + 1 P S h i and t max = 1 P t h l P S h i .
Considering the constraints of (P4) and (P5), the optimization problem can be rewritten as four simplification scenarios:
  • Scenario I: When a l ( 1 t i ) P S h i + b l P Q t i P S h i ( t i + 1 ) g i and a l ( 1 t i ) P S h i + b l P Q P max φ i , (P4) can be rewritten as
    max P i , t i y 1 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q a l ( 1 t i ) P S h i b l ) s . t . P 1 P i a l ( 1 t i ) P S h i + b l P Q t min t i t max t i max ( t 2 , t 3 )
There is a solution to Equation (21) if and only if max ( t 2 , t 3 ) t max is satisfied. First, because Equation (21) is a linear programming problem for t i and y 1 t i < 0 , t i * = max ( t 2 , t 3 , t min ) can be obtained according to monotonicity. Then, since 2 y 1 P i 2 < 0 , y 1 ( P i ) is a convex function with respect to P i , and it can be solved by derivation. Solution I can be expressed as
{ P i I * = [ P 0 ] P 1 a l ( 1 t i I * ) P S h i + b l P Q ρ i I * = 1 t i I * = 1 max ( t 2 , t 3 , t min )
where t 2 = b + b 2 4 a c 2 a is the positive root of the quadratic equation a t i 2 + b t i + c = 0 , a = a l g i P S h i > 0 , b = P S h i ( b l P Q ) g i , c = ( a l P S h i + b l P Q ) g i , t 3 = a l P S h i + b l P Q ( P max φ i ) a l P S h i , P 0 = 1 x * 2 ln 2 ( 1 + ψ i ) σ 2 g i and P 1 = max ( 0 , ( γ min ψ i ) σ 2 g i ) .
  • Scenario II: When P max φ i a l ( 1 t i ) P S h i + b l P Q and P max φ i t i P S h i ( t i + 1 ) g i , (P4) can be rewritten as
    max P i , t i y 1 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q a l ( 1 t i ) P S h i b l ) s . t . P 1 P i P max φ i t min t i t max t 4 t i t 3
There is a solution to Equation (23) if and only if D = [ t 4 , t 3 ] [ t min , t max ] is satisfied and we let D = [ t a , t b ] . Solved as in Equation (21), solution II can be expressed as
{ P i I I * = [ P 0 ] P 1 P max φ i ρ i I I * = 1 t i I I * = 1 t a
where P 0 , P 1 and t 3 are the same as in Scenario I: t 4 = ( P max φ i ) g i P S h i ( P max φ i ) g i .
  • Scenario III: When t i P S h i ( t i + 1 ) g i a l ( 1 t i ) P S h i + b l P Q and t i P S h i ( t i + 1 ) g i P max φ i , (P4) can be rewritten as
    max P i , t i y 1 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q a l ( 1 t i ) P S h i b l ) s . t . P 1 P i t i P S h i ( t i + 1 ) g i t min t i t max t i min ( t 2 , t 4 )
There is a solution to Equation (25) if and only if t min min ( t 2 , t 4 ) is satisfied. Using the Lagrange multiplier method to convert a constrained problem to an unconstrained problem max P i , t i L ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + P i g i σ 2 ) x * ( P i + P Q a l ( 1 t i ) P S h i b l ) + λ i ( t i P S h i ( t i + 1 ) g i P i ) , where λ i is the Lagrange multiplier, since 2 L P i 2 < 0 and 2 L t i 2 < 0 , the optimal solution can be obtained by the KKT condition as follows
{ P i I I I * = max ( P 1 , t i I I I * P S h i ( t i I I I * + 1 ) g i ) ρ i I I I * = 1 t i I I I * = 1 [ t 0 ] t min min ( t 2 , t 4 , t max )
where t 0 is the real root of the univariate cubic equation a ( t 0 ) 3 + b ( t 0 ) 2 + c t 0 + d = 0 , a = x * 2 ln 2 a l g i ( P S h i + ( 1 + ψ i ) σ 2 ) , b = x * 2 ln 2 a l g i ( 2 P S h i + 3 ( 1 + ψ i ) σ 2 ) , c = x * 2 ln 2 ( a l g i + 1 ) ( P S h i + ( 1 + ψ i ) σ 2 ) + 4 x * ln 2 a l g i ( 1 + ψ i ) σ 2 g i , d = x * 2 ln 2 ( a l g i + 1 ) ( 1 + ψ i ) σ 2 g i , t 2 and P 1 are the same as in Scenario I and t 4 is the same as in Scenario II.
In order to solve this univariate cubic equation, we use Girolamo Cardan’s formula. We let t 0 = Y b 3 a and change it into a standard univariate cubic equation Y 3 + p Y + q = 0 , where p = c a b 2 3 a 2 and q = d a b c 3 a 2 + 2 b 3 27 a 3 . The discriminant is Δ = ( q 2 ) 2 + ( p 3 ) 3 and can be divided into three cases: when Δ > 0 , there is a real root Y 1 = u + v ; when Δ = 0 , there are two real roots Y 1 = u + v and Y 2 = w · u + w 2 and when Δ < 0 , there are three real roots Y 1 = u + v , Y 2 = w u + w 2 v and Y 3 = w 2 u + w v , where u = q 2 Δ 3 , v = q 2 + Δ 3 and w = 1 + 3 i 2 .
  • Scenario IV: Considering the constraints of (P5), the optimization problem can be rewritten as
    max P i , t i y 2 ( P i , t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) x * ( P i + P Q a l ( 1 t i ) P S h i b l ) s . t . t i P S h i ( t i + 1 ) g i P i min ( a l ( 1 t i ) P S h i + b l P Q , P max φ i )   t min t i t max     t 1 t i min ( t 2 , t 4 )
There is a solution to Equation (27) if and only if D = [ t 1 , min ( t 2 , t 4 ) ] [ t min , t max ] is satisfied and we let D = [ t a , t b ] . Solved as in Equation (19), solution IV can be expressed as
{ P i I V * = t i I V * P S h i ( t i I V * + 1 ) g i ρ i I V * = 1 t i I V * = 1 arg max t i * { t a , t b , [ t 0 ] t a t b } y 2 ( t i * )
where y 2 ( t i ) = 1 2 log 2 ( 1 + ψ i + t i P S h i ( t i + 1 ) σ 2 ) x * ( t i P S h i ( t i + 1 ) g i + P Q a l ( 1 t i ) P S h i b l ) , t 1 = ( γ min ψ i ) σ 2 P S h i ( γ min ψ i ) σ 2 and t 0 , t 2 and t 4 are the same as in Scenario III.
The power allocation algorithm in this paper consists of an outer loop algorithm and an inner loop algorithm. Since only single-relay power is assumed to be allocated and other relay power allocations are known in the solution, a fixed-point iteration algorithm is needed to obtain joint power optimization for all relays (i.e., an outer loop algorithm based on a fixed-point iteration, given by Algorithm 1). In the single-relay power allocation, given the initial optimal EE x * and the segment s i of the EH model, the Dinkelbach iteration algorithm is used to obtain the optimal relay transmit power, the optimal PS ratio and the optimal EE of R i (i.e., the inner loop algorithm based on the Dinkelbach iteration, given by Algorithm 2). The workflow of this power allocation algorithm is shown in Figure 3.
Algorithm 1. Outer loop algorithm based on a fixed-point iteration.
1: Initialize the maximum iterations I max 1 and the maximum error toleranceδ. Set the iteration indexm = 0, the optimal relay transmission power P i 0 and the optimal PS ratio ρ i 0 , ( i = 1 , ... , N ) .
2:  repeat
3:   For fixed P m w and ρ m w , solve (P3) by Algorithm 2 to get P m + 1 w , ρ m + 1 w and η * , i = 1 , ... , N
4:   if | P i m + 1 P i m | < δ and | ρ i m + 1 ρ i m | < δ then
5:      Convergence = true
6:      return P i * = P i m + 1 , ρ i * = ρ i m + 1
7:   else
8:      m = m + 1 and Convergence = false
9:   end if
10:  until Convergence = true or m = I max 1
Algorithm2. Inner loop algorithm based on the Dinkelbach iteration.
1:We set ρ i = 1 and Calculate maximum number of segments Si.
2:for s i = 1 : S i .
3:   Initialize the maximum iterations I max 2 , the maximum error tolerance δ and the optimal energy efficiency x * . Set iteration index k = 0.
4:   repeat
5:       ψ i = w = 1 , w i N γ w , φ i = w = 1 , w i N P w , ϕ i = w = 1 , w i N ( P w + P Q P H w ( ρ w ) )
6:       if si = L then
7:          Solve the power allocation by using (16), (18), (20) to get Pi, ρi.
8:       else
9:          Solve the power allocation by using (22), (24), (26), (28) to get Pi, ρi.
10:       end if
11:       Update x k + 1 * .
12:       if | x k + 1 * x k * | < δ then
13:             Convergence = true
14:             return P i s i = P i , ρ i s i = ρ i , η s i = x k + 1 *
15:       else
16:          k = k + 1 and Convergence = false
17:       end if
18:    until Convergence = true or k = I max 2
19:end for
20: [ η * , s i ] = max ( η s i )
21:return P i m + 1 = P i s i ( s i ) , P i m + 1 = P i s i ( s i ) , η *
In the inner loop algorithm, in the worst case, the Dinkelbach iteration is I max 2 , and it needs to be performed ( L + 1 ) times. The outer loop algorithm contains the I max 1 inner loop algorithm, where the fixed-point iteration is N in the worst case. Therefore, the complexity is O ( ( L + 1 ) N I max 1 I max 2 ) .

4. Simulations

In this section, the simulation results will be analyzed, and the effectiveness of the proposed algorithm will be proven. The simulation scenario is shown in Figure 1, where the distance between S and D is d = 10   m and the relay is randomly distributed in the circle with the midpoint of S and D as the center and radius, expressed as r = 2   m . Like [12], all the channels are Rayleigh fading channels, where the path loss exponent α is α = 3 and the small-scale fading is Rayleigh fading with a mean value of 0 and variance of 1. In accordance with [12,20,21], we set the simulation parameters as follows: P S = 30   dBm , P C = 10   mW , P Q = 0.1   mW , σ 2 = 70   dBm and R min = 3   bps/Hz . The parameters of the EH model in Figure 2 are shown in Table 1.
Figure 4 shows the variation of the EE in different EH models under the same relay total transmit power threshold. It can be seen that the EE increased as the number of relays increased. This is because, compared with a single relay, multiple relays can improve the information transmission rate. Since the relay uses the harvested energy for DF without consuming fixed energy, the increase in the number of relays can also improve the EE of the system. In addition, it can also be observed that the performance of the piecewise linear EH model was better than that of the linear EH model in [7]. This is because the piecewise linear EH model improves the accuracy of the traditional linear EH one. The more segments, the better the performance, but its growth rate will decrease with the increase in the segment number. Therefore, the piecewise linear EH model with L = 5 was selected for simulation analysis in this paper.
Figure 5 shows that the EE changed with the total transmit power threshold of relays P max , and the number of relays was N = 4 . It can be seen that the EE increased with the increase of P max , and the increase gradually decreased until it no longer increased. This was because the energy harvested by the EH model was limited. When P max exceeds the total power that can be used after EH, the EE will not increase with the increase of P max . In this simulation scenario, when P max > 40   mW , the EE would not change significantly with the increase in P max .
Figure 6 shows that the EE varied with the different transmit powers of source P S . It can be seen that the EE first increased with the increase of P S and then decreased with the increase of P S . This is because when P S is very small, the energy harvested by the relay may not be enough to complete the information transmission, which makes the EE relatively low. With the increase of P S , the relay can harvest more sufficient energy to complete the information transmission, so the EE correspondingly increases. However, when P S exceeds a certain value, the energy consumption is too high, and the EE will decrease with the increase of P S . In this simulation scenario, when P S = 18   dBm , the EE is the highest. In addition, compared with the single-relay in [13], multi-relay cooperation can effectively improve performance.
Figure 7 shows the change in EE with the noise power σ 2 . The lower the noise power, the higher the energy efficiency. This is because the reduction of noise power can effectively improve the information transmission rate. Figure 8 shows the change in EE with a different path loss exponent α . The larger the path loss exponent α , the smaller the channel gain, the smaller the information transmission rate, the smaller the harvested energy and the smaller the energy efficiency. Figure 9 shows the change in EE with the distance d between S and D. The larger the distance d , the smaller the information transmission rate and the smaller the energy efficiency. Figure 10 shows the change in EE with the position of relays. The closer the relay is to the source node, the higher the energy efficiency.

5. Conclusions

In this paper, a power allocation algorithm based on a nonlinear EH model is proposed for a two-hop half-duplex multiple cooperative forwarding relays system. The fixed-point iterative algorithm and Dinkelbach iterative algorithm are used to solve the optimization problem by jointly optimizing the transmit power and PS ratios of multiple relays to maximize the EE. The simulation results show that, compared with the conventional linear EH model, the designed EH model had better performance in terms of EE. Meanwhile, the superior EE performance of the multi-relay system was also verified by comparison with different numbers of relays. Compared with the linear model, when the piecewise linear model (n = 5) was used, the energy efficiency increased by 2%. Compared with a single relay, when four relays were used for cooperation, the energy efficiency increased by 26%.

Author Contributions

Conceptualization, H.P. and Q.Z.; methodology, H.P.; software: H.P.; validation, H.P. and Q.Z.; formal analysis, H.P.; investigation, H.P.; resources, Q.Z.; data curation, H.P.; writing, H.P.; visualization, H.P.; supervision, Q.Z.; project administration, Q.Z.; funding acquisition, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (61971239).

Data Availability Statement

The data and codes presented in this study are available from the corresponding authors by request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fettweis, G.; Zimmermann, E. ICT energy consumption-trends and challenges. In Proceedings of the 11th International Symposium on Wireless Personal Multimedia Communications, Lisbon, Portugal, 27–29 November 2008; pp. 2–6. [Google Scholar]
  2. Bianzino, A.; Chaudet, C.; Rossi, D.; Rougier, J.-L. A Survey of Green Networking Research. IEEE Commun. Surv. Tutor. 2010, 99, 1–18. [Google Scholar] [CrossRef] [Green Version]
  3. Want, R.; Farkas, K.I.; Narayanaswami, C. Energy Harvesting and Conservation. IEEE Pervasive Comput. 2005, 4, 14–17. [Google Scholar] [CrossRef]
  4. Singh, K.; Ku, M.; Lin, J.; Rathanarajah, T. Toward Optimal Power Control and Transfer for Energy Harvesting Amplify-and-Forward Relay Networks. IEEE Trans. Wirel. Commun. 2018, 17, 4971–4986. [Google Scholar] [CrossRef] [Green Version]
  5. Ozel, O.; Tutuncuoglu, K.; Yang, J.; Ulukus, S.; Yener, A. Resource management for fading wireless channels with energy harvesting nodes. In Proceedings of the IEEE Infocom, Shanghai, China, 10–15 April 2011. [Google Scholar]
  6. Gurrala, K.K.; Navya, K.; Sravya, M.; Ramakrishna, E. Maximized Energy Efficiency Based Power Allocation Strategy in Wireless Cooperative Network. In Proceedings of the 2019 TEQIP III Sponsored International Conference on Microwave Integrated Circuits, Photonics and Wireless Networks (IMICPW), Tiruchirappalli, Tamil Nadu, India, 22–24 May 2019; pp. 64–68. [Google Scholar]
  7. Zhang, C.; Du, H.; Ge, J. Energy-Efficient Power Allocation in Energy Harvesting Two-Way AF Relay Systems. IEEE Access 2017, 5, 3640–3645. [Google Scholar] [CrossRef]
  8. Orhan, O.; Erkip, E. Optimal Transmission Policies for Energy Harvesting Two-hop Networks. In Proceedings of the Information Sciences & Systems, Princeton, NJ, USA, 21–23 March 2012. [Google Scholar]
  9. Chen, Y.; Shi, R.; Feng, W.; Ning, G. AF Relaying with Energy Harvesting Source and Relay. IEEE Trans. Veh. Technol. 2017, 66, 874–879. [Google Scholar] [CrossRef] [Green Version]
  10. Ng, D.W.K.; Lo, E.S.; Schober, R. Wireless information and power transfer: Energy efficiency optimization in OFDMA systems. IEEE Trans. Wireless Commun. 2013, 12, 6352–6370. [Google Scholar] [CrossRef] [Green Version]
  11. Boshkovska, E.; Ng, D.W.K.; Zlatanov, N.; Schober, R. Practical Non-Linear Energy Harvesting Model and Resource Allocation for SWIPT Systems. IEEE Commun. Lett. 2015, 19, 2082–2085. [Google Scholar] [CrossRef] [Green Version]
  12. Shi, L.; Ye, Y.; Hu, R.Q.; Zhang, H. Energy Efficiency Maximization for SWIPT Enabled Two-Way DF Relaying. IEEE Signal Process. Lett. 2019, 26, 755–759. [Google Scholar] [CrossRef]
  13. Lu, G.; Shi, L.; Ye, Y. Maximum Throughput of TS/PS Scheme in an AF Relaying Network with Non-Linear Energy Harvester. IEEE Access 2018, 6, 26617–26625. [Google Scholar] [CrossRef]
  14. Li, T.; Fan, P.; Letaief, K.B. Outage Probability of Energy Harvesting Relay-Aided Cooperative Networks Over Rayleigh Fading Channel. IEEE Trans. Veh. Technol. 2016, 65, 972–978. [Google Scholar] [CrossRef] [Green Version]
  15. Shi, L.; Zhao, L.; Liang, K.; Chu, X.; Wu, G.; Chen, H. Profit maximization in wireless powered communications with improved non-linear energy conversion and storage efficiencies. In Proceedings of the 2017 IEEE International Conference on Communications (ICC), Paris, France, 21–25 May 2017; pp. 1–6. [Google Scholar]
  16. Guo, J.; Zhu, X. An improved analytical model for RF-DC conversion efficiency in microwave rectifiers. In Proceedings of the 2012 IEEE/MTT-S International Microwave Symposium Digest, Montreal, QC, Canada, 17–22 June 2012; pp. 1–3. [Google Scholar]
  17. Grover, P.; Sahai, A. Shannon meets Tesla: Wireless information and power transfer. In Proceedings of the IEEE International Symposium on Information Theory Proceedings, Austin, TX, USA, 13–18 June 2010; pp. 2363–2367. [Google Scholar]
  18. Doosti-Aref, A.; Ebrahimzadeh, A. Adaptive Relay Selection and Power Allocation for OFDM Cooperative Underwater Acoustic Systems. IEEE Trans. Mobile Comput. 2017, 99, 1. [Google Scholar] [CrossRef]
  19. Dinkelbach, W. On nonlinear fractional programming. Manag. Sci. 1967, 13, 492–498. [Google Scholar] [CrossRef]
  20. Gupta, A.; Singh, K.; Sellathurai, M. Time-Switching EH-Based Joint Relay Selection and Resource Allocation Algorithms for Multi-User Multi-Carrier AF Relay Networks. IEEE Trans. Green Commun. Netw. 2019, 3, 505–522. [Google Scholar] [CrossRef]
  21. Li, C.; Wu, P.; Xia, M. Energy Efficiency Maximization of AF Relaying SWIPT Systems with Energy Recycling. In Proceedings of the 2019 IEEE 89th Vehicular Technology Conference (VTC2019-Spring), Kuala Lumpur, Malaysia, 28 April–1 May 2019; pp. 1–5. [Google Scholar]
Figure 1. The system model.
Figure 1. The system model.
Algorithms 14 00155 g001
Figure 2. Comparison of energy harvesting in different EH models [7,10,16].
Figure 2. Comparison of energy harvesting in different EH models [7,10,16].
Algorithms 14 00155 g002
Figure 3. The workflow of the power allocation algorithm.
Figure 3. The workflow of the power allocation algorithm.
Algorithms 14 00155 g003
Figure 4. Energy efficiency versus number of relays under different EH models.
Figure 4. Energy efficiency versus number of relays under different EH models.
Algorithms 14 00155 g004
Figure 5. Energy efficiency versus the threshold of the relays’ total transmit power P max .
Figure 5. Energy efficiency versus the threshold of the relays’ total transmit power P max .
Algorithms 14 00155 g005
Figure 6. Energy efficiency versus the transmit power of source P S .
Figure 6. Energy efficiency versus the transmit power of source P S .
Algorithms 14 00155 g006
Figure 7. Energy efficiency versus number of relays with different noise power σ 2 .
Figure 7. Energy efficiency versus number of relays with different noise power σ 2 .
Algorithms 14 00155 g007
Figure 8. Energy efficiency versus number of relays with different path loss exponent α .
Figure 8. Energy efficiency versus number of relays with different path loss exponent α .
Algorithms 14 00155 g008
Figure 9. Energy efficiency versus number of relays with different distance d between S and D.
Figure 9. Energy efficiency versus number of relays with different distance d between S and D.
Algorithms 14 00155 g009
Figure 10. Energy efficiency versus the positions of the relays.
Figure 10. Energy efficiency versus the positions of the relays.
Algorithms 14 00155 g010
Table 1. Parameters of the EH model.
Table 1. Parameters of the EH model.
Linear EH model P t h = [ 0 , 60 ]   mW , a = [ 0 . 3957 , 0 ] , b = [ 0 , 23 . 742 ]   mW
Piecewise linear EH model (L = 4) P t h = [ 0 , 5 , 15 , 60 ]   mW , a = [ 0.4 , 0.88 , 0.2877 , 0 ] , b = [ 0 , 2.4 , 6.4845 , 23.7465 ]   mW
Piecewise linear EH model (L = 5) P t h = [ 0 , 5 , 15 , 25 , 60 ]   mW , a = [ 0.4 , 0.88 , 0.896 , 0.1139 , 0 ] , b = [ 0 , 2.4 , 2.64 , 16.9125 , 23.7465 ]   mW
Piecewise linear EH model (L = 6) P t h = [ 0 , 5 , 15 , 25 , 40 , 60 ]   mW , a = [ 0 . 4 , 0 . 88 , 0 . 896 , 0 . 2406 , 0 . 018625 , 0 ] , b = [ 0 , - 2 . 4 , - 2 . 64 , 13 . 745 , 22 . 629 , 23 . 7465 ]   mW
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, H.; Zhu, Q. Energy-Efficient Power Allocation in Non-Linear Energy Harvesting Multiple Relay Systems. Algorithms 2021, 14, 155. https://0-doi-org.brum.beds.ac.uk/10.3390/a14050155

AMA Style

Pan H, Zhu Q. Energy-Efficient Power Allocation in Non-Linear Energy Harvesting Multiple Relay Systems. Algorithms. 2021; 14(5):155. https://0-doi-org.brum.beds.ac.uk/10.3390/a14050155

Chicago/Turabian Style

Pan, Huifang, and Qi Zhu. 2021. "Energy-Efficient Power Allocation in Non-Linear Energy Harvesting Multiple Relay Systems" Algorithms 14, no. 5: 155. https://0-doi-org.brum.beds.ac.uk/10.3390/a14050155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop