Abstract
As a massive number of nodes are expected to be connected to each other in next-generation wireless networks, the power supply of such multiple-antenna nodes is a challenge and requires new and sustainable solutions. Radio frequency energy harvesting (EH) is among the promising approaches to solve this issue. In this paper, two new mathematically tractable nonlinear EH (NL-EH) models are proposed, for which the system performance is evaluated in terms of average harvested power, throughput, and bit error probability. The system consists of one multiantenna power-constraint source harvesting its power from a dedicated power beacon and transmitting its signal to a destination equipped with multiple antennas. For a comprehensive analysis of the system, closed-form expressions are derived for Nakagami- fading channels and the special case of Rayleigh channels. For comparison purposes, the performance derivations in Nakagami- fading channels are provided for linear and piece-wise linear EH models given in the literature. Besides, the simulation results are obtained by applying the Monte-Carlo method for NL-EH models existing in the literature. The results provide a broadening view of EH systems and thoroughly compare the proposed NL-EH models with linear, piece-wise linear, and NL-EH models available in the literature. Hence, these provide better insights towards the analyses and design of EH systems.
1. Introduction
In next-generation wireless communications, a huge number of nodes with multiple antennas will be connected to each other, and power requirements will be the critical constraint for sustainability. Wireless energy-harvesting (EH) methods are considered among promising solutions to the energy shortage problem of these nodes. Considering the constraint energy resources and green communication, radio frequency (RF) EH provides a small step to solve the energy shortage because it replaces the usage of batteries and electronic power supplies at the sensor nodes. The EH node is equipped with an EH circuit which provides the ability to harvest the incoming signal energy. Wireless power transfer (WPT) and simultaneous wireless information and power transfer (SWIPT) are two EH methods considered in the design of battery-less nodes [1–3]. Based on the EH circuit feature, the power-constraint node can harvest the total incoming signal power or simultaneously harvest energy and process the information [1, 2]. Power splitting (PS) and time switching (TS) are examples of the SWIPT method. In the WPT, a dedicated power beacon (PB) or ambient power nodes are the sources from which the power-constraint nodes harvest their power [3]. In the above-mentioned type of EH circuits and EH sources, the important factors are the amount of harvested power and the characteristic of the EH circuit which directly affect the EH system performance. Linear, piece-wise linear, and nonlinear EH (NL-EH) models characterize the amount of harvested power [1, 4–8]. Most works in the literature consider the linear EH model which has a direct proportion between the input and output powers of the EH circuit [1, 9, 10]. On the other hand, the piece-wise linear EH model is considered in the literature [4, 5] which provides the same amount of energy as in the linear model for the small amount of input power. Different from the linear EH model, it saturates the output power to a fixed maximum value at high input powers.
With the measured data from the practical experimental tests given in [11–24], it is figured out that the input and output power relation of the EH circuit is nonlinear and the output power saturates to a maximum power even increasing the input power. Applying the curve fitting over the measurement data, some NL-EH models are proposed in the literature [6–8]. This characteristic of the EH node is different from the linear EH model in which the harvested power is directly proportional to the input power and the harvested power increases with rise in the input power of the EH circuit. Consequently, this results in an overestimation of the amount of harvested power, and this yields an inaccurate understanding of EH system performance. The system performance in terms of outage probability and bit error probability (BEP) is investigated well for the linear EH models [1, 2] and the piece-wise linear EH model [3] given in the literature. However, the system performance of the NL-EH models is mostly studied in terms of the amount of harvested energy and its statistical characteristics [6, 7, 25] since these NL-EH models are not mathematically tractable when investigating their bit error performances. Additionally, a small number of studies investigated the throughput [8]. This problem comes from the fact that for the NL-EH models given in the literature, it is impossible to calculate the probability density function (PDF) or cumulative density function (CDF) of the received signal-to-noise ratio (SNR), which prevents the derivations of BEP or throughput.
The ergodic capacity and throughput of dual-hop (DH) amplify-and-forward (AF) relaying system obtained from outage analysis are investigated in [1]. The results are provided for PS, TS, and ideal EH modes where a linear EH model is applied at the power-constraint node. The results are obtained for different parameters, and the authors show the discrepancy of the PS and TS EH modes in terms of the throughput from the ideal EH mode. The closed-form expressions for outage probability and BEP are derived in SWIPT-based DH AF and DH decode-and-forward (DF) relaying systems [2]. The channels are subject to Nakagami- fading, and the power-constraint relay node applies the linear EH model. The performance is provided for different system parameters, and generally, the PS EH mode outperforms the TS EH mode. System performance of the WPC DH AF relaying system is investigated in [26] where both S and R are power-constraint nodes and harvest their power from a dedicated PB. BER of the system is obtained for different parameters, and theoretical derivations are derived considering the linear EH model.
A DH AF relaying system with a power-constraint relay is considered in [4] where the relay applies the TS mode to harvest its power based on the piece-wise linear EH model. The outage capacity of the studied system is derived considering Nakagami- fading channels. A multirelay system where source (S) and destination (D) communicate with each other with the help of power-constraint relay nodes is investigated under the Rayleigh fading channel assumption [27]. Outage probability and throughput of the system are derived for the piece-wise linear EH model and imperfect channel state information (CSI). In the studied system, both transmit antenna and relay selections are considered for the transmission protocol. In [3], an overlay cognitive radio with a power-constraint secondary transmitter (ST) is investigated. The authors assumed that ST harvests its power from the primary transmitter and applies the piece-wise linear EH model. The closed-form expressions for the outage probabilities and BEPs are derived considering both linear and piece-wise linear EH models which provide a good insight for comparison purposes.
The throughput of the DH AF relaying system considering Rayleigh fading is studied in [5]. TS and PS EH modes with a segmented piece-wise linear EH model are applied in the battery-less relay node. The authors show that the PS mode outperforms the TS mode. SWIPT-based cooperative nonorthogonal multiple access (NOMA) system performance is studied in [28] where the near user is the power-constraint and applies a segmented piece-wise linear EH model. The theoretical derivations of outage probability and system throughput are provided for the Rayleigh fading channel.
In [7], a NL-EH model based on logistic (sigmoidal) function is proposed. Assuming a nonconvex optimization and simplifying it to a convex optimizing problem, the total harvested power is maximized and the average harvested power is calculated. The rate-energy trade-off region considering both linear and NL-EH models is obtained in [29]. The harvested power versus the number of receive antennas is investigated. In [30], the average harvested power of the point-to-point MIMO WPT system is calculated based on different system parameters. The authors study the maximization of harvested power jointly by optimizing the PS ratio factor and the transmit energy of PB. A MISO NOMA SWIPT system is investigated in [31] where the authors study the physical layer security by adopting the NL-EH model from [7]. The authors simplify the nonconvex optimization problem to a convex one.
The quadratic NL-EH model is proposed and studied in a multiuser orthogonal frequency-division multiple access (OFDMA) system [6]. Results are given compared to the linear EH model. Moreover, it is assumed that the parameters of the proposed function in [6] take values of and , and which is shown in Table 1. However, a high input power results in a negative value for harvested power. Hence, this proposed model is only valid for small input powers. A heuristic NL-EH model is proposed in [8] applying curve fitting. The throughput of the AF cooperative relaying system using PS and TS modes is investigated where the results are numerically obtained since the proposed model is not mathematically tractable.
The NL-EH models given above are based on the curve fitting of experimental measured data. The harvested powers given in [6–8] include the polynomial, exponential, or cubic functions of the input power. Consequently, all these NL-EH models given in the literature are not mathematically tractable, and the harvested power is not provided as PDF and CDF to further calculate the EH system performance. These are the main motivations behind our article.
In this paper, we propose two NL-EH models which are mathematically tractable in terms of the PDF, CDF of harvested power, and the system performance at the expense of a small gap in the accuracy between measured data points and those from our models. Apart from the NL-EH models given in the literature, our models depend directly on the system input power rather than the polynomial, exponential, or cubic functions of it. The considered system consists of a PB from where a power-constraint S harvests energy which is then used to transmit data to a D. The channels are assumed subject to Nakagami- and Rayleigh fading. The former is assumed when the channels between nodes are stronger and the distances are short. In the latter, the link between nodes is subjected to the deleterious effect of the environment in the lack of a line-of-sight component.
For both models, PDF and CDF expressions for the harvested power are obtained based on the curve fitting of the measurement data from [11, 15]. Moreover, the average harvested powers are calculated using the PDF of harvested power. Furthermore, the outage probability, throughput, and BEP of the considered system are derived in closed forms for the proposed two NL-EH models. Note that the NL-EH models considered in the literature do not provide closed-form expressions for the system performance since they are not mathematically tractable [6–8]. In the papers considering these models, only the amount of harvested power at the EH receiver is studied without investigating the system BEP or throughput performance.
Numerical results provide a comprehensive analysis for all the considered EH models in terms of system performance. The results show that the distance and number of antennas have crucial effects on the system performance which gives an insight into maximizing the system performance.
The main contributions of the paper are as follows: (i)Two mathematically tractable NL-EH models based on the curve fitting of measured data are proposed(ii)PDF and CDF of the harvested power for the proposed NL-EH models are provided(iii)Closed-form expressions for the average harvested power, throughput, and BEP are derived, and the results are supported by Monte-Carlo simulations(iv)Theoretical derivations of the system performance are provided for both Nakagami- and Rayleigh fading(v)Results are compared with existing linear, piece-wise linear, and NL-EH models from the literature
The paper is organized as follows. Section 2 describes the considered system model. The proposed NL-EH models are presented in Section 3. Expressions for the average harvested powers are obtained in Section 4. Section 5 deals with the derivation of throughput expressions. Bit error rate (BER) analyses are performed in Section 6. Theoretical and simulation results are presented in Section 7. Section 8 concludes the paper.
2. System Model
The considered half-duplex RF EH system where a power-constraint S harvests its power from a PB is given in Figure 1. S and D are equipped with and receive antennas, respectively, while PB is equipped with one antenna. S transmits its signal from one antenna to D using the harvested power. and denote the distances of links PB→S and S→D, respectively. Additionally, the path loss is represented as where , and is the path-loss coefficient. is S’s transmitted signal with . and represent the channel gain vectors of links PB→S and S→D, respectively. Furthermore, is the additive white Gaussian noise (AWGN) vector with entries as i.i.d complex Gaussian random variable for . Transmission of the S signal takes seconds where are allocated for harvesting energy and the remainder for IP. It is assumed that the CSI is perfectly known at node D and channels are assumed to be flat fading and remain constant during each transmission interval and vary independently from an interval to another. Line-of-sight (LoS) and distance are two determinant factors considered in the RF EH communication systems. Since the path loss causes a substantial decrease in the amount of harvested power, the EH communication is possible for short distances and in the presence of a LoS component. Generally, Rician fading is a proper channel model with its LoS component. However, mathematical derivations of the system performance considering the Rician PDF and CDF are not tractable. Meanwhile, Nakagami- fading is a good approximation of the Rician model [32, 33]. Hence, it is assumed that both links PB→S and S→D are subject to Nakagami- fading with parameters and , respectively. Here, we consider that the channel power is identical for each link. In other words, and for , and , respectively. In the first time slot, the input power of the EH circuit provided in Figure 1 is calculated as where is the energy-bearing signal transmitted from PB. In the second time slot, the received signal at D is given as where , , and is the transmit power of S. Considering (1), we have where is given in Figure 1. The received SNR after applying the maximum ratio combination (MRC) at node D is given using (2) as where . For simplicity, we define and in the sequel. The CDF and PDF of r.v.s and are given as [33]

respectively, where , , , , , , , , and .
3. Linear and NL-EH Models
In the literature, different linear and NL-EH models are provided for the performance evaluation of EH communication systems [1, 4–8]. These EH models are listed in Table 1. For the linear EH model, the harvested power is linearly dependent on the input power of the EH circuit with a constant ratio parameter, namely, the conversion efficiency (). However, based on the data measurements from designed practical EH circuits, linear EH models overestimate the amount of the harvested energy especially for high input powers which consequently results in a misunderstanding of the system performance. Moreover, the EH circuits have nonlinear properties which can be modeled by a nonlinear function of the amount of harvested power with respect to the input power. Furthermore, the piece-wise linear EH model is provided in the literature where the harvested amount of energy is saturated to a threshold power as while acting the same as the linear EH model for the small input powers as seen from Table 1. In addition to these EH models, NL-EH models are also proposed, applying curve fitting from the measurement data experiments from different EH circuits (Table 1). Please note that the theoretical derivations of system performance are not mathematically tractable considering these known NL-EH models. However, the linear and piece-wise linear EH models are mathematically tractable and provide a general and simple analysis of the system. On the other hand, the proposed NL-EH models in the sequel enable mathematically tractable theoretical derivations of system performance. This provides an advantage on the design and evaluation of EH systems.
3.1. Proposed NL-EH Models
In this subsection, two EH models are proposed with their relations between input and output powers. The first proposed NL-EH model (NL-EH-1) is given as
Here, denotes the maximum saturation power of the EH circuit, and , , and are constants to be determined by curve fitting. Applying (7), the CDF of the harvested power of NL-EH-1 related to the input power is calculated as where , , and
Moreover, the PDF of NL-EH-1 is calculated using (8) as where is given in (6).
The second proposed NL-EH model (NL-EH-2) is given as
Similar to NL-EH-1, the CDF of harvested power is obtained as where . Taking the derivation of (12), the PDF of NL-EH-2 is calculated as
Figures 2(a) and 2(b) describe the curve fitting of the measured experimental data compared to the linear EH and NL-EH models existing in the literature and the two proposed models in (7) and (11). The results of the curve fitting are obtained based on the measured data from [11, 15]. The measured data are denoted by symbols while the fitted data are represented by lines and dashed lines.

(a)

(b)
It is seen from both Figures 2(a) and 2(b) that both the proposed NL-EH models fit better to data points than the linear and piece-wise linear models [1, 4]. However, as stated, the challenging issue of the existing NL-EH models are the difficulties in mathematical derivations of CDF and PDF of the harvested power and theoretical derivations of throughput as well as BER system performances. In addition, it is observed that apart from the NL-EH model in [6], the other NL-EH models provide similar curve-fitting functions. Moreover, the NL-EH model in [8] provides the best fit compared to all other NL-EH models. However, despite these similarities, these models given in the literature are not mathematically tractable which prevent theoretical investigation of the EH system performance. Meanwhile, the proposed NL-EH models which have an acceptable goodness in fitting are mathematically tractable.
Moreover, the goodness of the fits considering the NL-EH models are provided in Table 2 for [11, 15], respectively. (-squared) is the parameter that is statistically calculated from curve-fitting functions and data measurements [34]. The strength of the fit of a model is most commonly evaluated using which always takes values between and . This statistical parameter measures how well the data is close to being fit. The closer the value to 1, the fit is better. Considering both Figures 2 and 3 and Table 2, the proposed NL-EH model functions provide similar results as their counterparts given in the literature. However, the important factor in choosing the NL-EH model is having mathematical tractability while being well fit to the data measurements. In other words, a trade-off between accuracy and mathematical tractability is considered. Furthermore, the coefficient parameters of NL-EH-1 and NL-EH-2 considered in (7) and (11) are obtained with curve fitting and provided in Tables 3 and 4, respectively. As provided in both Tables 3 and 4, the parameters of these models are obtained for measurement results of [11–24].

(a)

(b)
Figures 3(a) and 3(b) show the histograms of the PDFs of NL-EH-1 and NL-EH-2 given by (10) and (13) for (6), respectively. Moreover, the histograms of both PDFs are provided for different channel powers and parameters where it is assumed that , and . It is seen from both Figures 3(a) and 3(b) that the histograms overlay the theoretical PDFs for the intended distributions. This verifies the accuracy of the derivations of both PDFs for NL-EH models. Please note that for all cases, the area of shaded regions is always equals to one which verifies the PDF calculations given in (10) and (13). Also, note that both histograms are bounded to the given boundaries defined in (10) and (13).
4. Average Harvested Power of the Proposed NL-EH Models
In this section, average harvested power of the two proposed NL-EH models are investigated. Additionally, for comparison purposes, the average harvested powers of the linear and piece-wise linear EH models are calculated.
4.1. Average Harvested Power of NL-EH-1 Model
The average harvested power of the NL-EH-1 model given in (7) is calculated as where
Here, is given in (6). Substituting (6) in (15) and (16) and applying ([35], 3.383-10), we have where is the upper incomplete Gamma function ([35], 8.350-2). Finally, the average harvested power for NL-EH-1 given in (7) is obtained by substituting (17) and (18) in (14).
4.2. Average Harvested Power of NL-EH-2 Model
The average harvested power of NL-EH-2 model given in (11) is calculated as where is given in (6). Applying ([35], 3.383-10), the average harvested power is obtained as
4.3. Average Harvested Power of Linear and Piece-Wise Linear EH Models
The average harvested powers of linear and piece-wise linear EH models given in Table 1 are written as
respectively, where . Substituting (6) in (21), simplifying and applying ([35], 3.326-2), the average harvested power of the linear EH model for the considered system is obtained as
Moreover, for the piece-wise linear EH model in (22), by substituting (6), we have where is given in (9). Here, (24) is calculated by ([35], 3.351-1).
Based on the obtained expressions given in equations (14), (19), (21), and (22), the average harvested power depends on the PDF of the input power as well as the channel path loss parameter, number of harvesting antennas , and the transmit power of PB, as seen from (6).
5. Outage and Throughput Analyses
In this section, closed-form expressions are derived for the outage probability and throughput of the overall system PB→S→D considering both proposed NL-EH models in Nakagami- fading channels. The throughput of the considered system is calculated as where and bits/sec/Hz is the target rate. Moreover, is the CDF of the received SNR given in (4).
5.1. Outage Probability and Throughput of NL-EH-1 Model
Considering (4) and (7), the conditional CDF of the received SNR is given as where is given in (5). Averaging over using (6), we have
(28) is rewritten by substituting (5) and (6) as where . By defining and after some mathematical simplification, (29) is rewritten as where . Using ([35], 1.111), we have
Applying ([35], 3.471-9), the CDF of the received SNR at D is calculated as where and
and is the modified Bessel function of the second kind ([35], 8.407-2). The throughput of the considered system for the NL-EH-1 model given in (7) is calculated by substituting from (32) into (26).
5.1.1. Special Case of Rayleigh Fading for NL-EH-1 Model
For the case where both links from PB→S and S→D are exposed to Rayleigh fading and both S and D are equipped with one antenna, the CDF and PDF of (5) and (6) are simplified as
respectively, where , , , and . Substituting (34) and (35) into (28), and taking and simplifying, we have
Using ([35], 3.471-9), the CDF is calculated as where . Then, substituting (37) into (26), the throughput for the special case of Rayleigh fading is obtained.
5.2. Outage Probability and Throughput of NL-EH-2 Model
Considering (4) and (11), the conditional CDF of the received SNR is given as
Averaging over (6), we have
Substituting (5) and (6) into (39), and simplifying, we have where . Using ([35], 1.111), we have
Applying ([35], 3.471-9) and simplifying, we have where
Then, substituting (42) into (26), the system throughput for NL-EH-2 is obtained.
5.2.1. Special Case of Rayleigh Fading for NL-EH-2 Model
The outage probability and throughput of the considered system for the NL-EH-2 model given in (11) is investigated under the assumption of Rayleigh fading for both links PB→S and S→D. Applying (34) and (35) into (39), we have
Using ([35], 3.471-9), the CDF of received SNR at node D is obtained as
Finally, the throughput of the considered system for NL-EH-2 is calculated by substituting (45) into (26).
5.3. Outage Probability and Throughput of Linear and Piece-Wise Linear Models
Considering Table 1, the conditional CDF of the received SNR at D given in (4) for both linear and piece-wise linear EH models is written as where is given in (5). For the linear EH model, the CDF is calculated as
Substituting (5) and (6) into (47) and applying ([35], 3.471-9), the CDF in (47) is rewritten as
For the piece-wise linear EH model given in Table 1, the CDF is expressed as
Considering (5) and (6), and simplifying, we have where . Here, is calculated by substituting (5) and (6) into (49) as
Moreover, applying ([35], 1.211-1), [36, 11], ([35], 9.31-2), and [36, 26], (51) is obtained as where and is the Meijer-G function ([35], 9.301). Finally, the throughputs of linear and piece-wise linear EH models are obtained by substituting (48) and (50), respectively, into (26).
Different from the NL-EH models previously given in the literature, we obtained closed-form outage probability expressions for the Nakagami- channels in equations (32) and (42) and for the special case of Rayleigh fading channels in equations (37) and (45) for NL-EH-1 and NL-EH-2 models, respectively.
6. BEP Analyses
In this section, closed-form expressions for the BEP of the considered system for NL-EH-1, NL-EH-2, linear, and piece-wise linear EH models are derived. The BEP is calculated as where is the CDF obtained in (32), (37), (42), and (45) for NL-EH-1 and NL-EH-2 in Nakagami- and Rayleigh fading channels, respectively. Here, constants and are modulation-specific parameters [33] and , being the modulation order. Under the assumption of Gray mapping, (53) provides a tight upper bound at high and medium SNR regions. Substituting (32), (37), (42), and (45) into (53), (53) is obtained in the general form of where and
(55) is calculated using ([35], 3.381-4). In (54), , , , and are obtained by substituting (32), (37), (42), and (45), respectively, as explained in the following subsections.
6.1. BEP Analysis for NL-EH-1 Model
The BEP for the NL-EH-1 model given in (7) is calculated for the Nakagami- fading channel, from (54), by taking for
Furthermore, assuming in (56), we have where and . Using ([35], 6.631-3) in (57) and simplifying, we have where and is the Whittaker function ([35], 9.22). Finally, the closed-form BEP for the NL-EH-1 model is calculated by substituting (58) into (54).
6.1.1. Special Case of Rayleigh Fading for NL-EH-1
The BEP for the Rayleigh fading channel is calculated replacing for
in (54). Moreover, using ([35], 6.614-4) in (59), we have where . Thus, the closed-form BEP expression for the NL-EH-1 model in Rayleigh fading is calculated by substituting (60) into (54).
6.2. BEP Analysis for NL-EH-2 Model
The BEP for the NL-EH-2 model given in (11) is calculated for the Nakagami- fading channel from (54) by replacing for
Assuming in (61), we have where and . Using ([35], 6.631-3) and simplifying, (62) is rewritten as
Finally, substituting (63) into (54), the closed-form BEP expression for NL-EH-2 is obtained.
6.2.1. Special Case of Rayleigh Fading for NL-EH-2
The BEP is calculated from (54) by taking and replacing
Using ([31], 6.614-4), we have where and . Hence, the BEP expression for the NL-EH-2 model under the Rayleigh fading channel is derived substituting (65) into (54).
6.3. BEP Analysis for Linear and Piece-Wise Linear EH Models
The conditional BEP of the system considering both linear and piece-wise linear EH models can be expressed as where is given in (46) and (5). Moreover, is obtained in (55), and
Applying ([35], 3.381-4) into (67) and after some mathematical simplifications, we have where . For the linear EH model, the unconditional BEP of the considered system is calculated as where where and are given in (67) and (6), respectively. Substituting, simplifying, and using Equation (10) of [36], we have
Using ([35], 7.813-1), we have
Finally, the BEP expression for the linear model is obtained substituting (70) and (73) into (69).
Now, considering the piece-wise linear EH model, the BEP of the considered system is written as where is calculated by substituting into (66), and is calculated using (9). Moreover, we have where where and are given in (67) and (6), respectively. Substituting into (77), simplifying, and applying , the upper bound of integration is changed from to , and (77) is rewritten as where where , , and are the integration boundaries for , , and , respectively. Considering and into (79) and applying [3, 36], we have where and and for and , respectively. Applying ([35], 07.20.07.0001.01) into (80), we have where is the Kummer confluent hypergeometric function given in ([35], 07.20.02.0001.01). Considering in (79) and applying Equation (3) in [36], we have where
Here, (84) and (85) are calculated using ([35], 3.351-2). Then, (83) is calculated by substituting (84) and (85). Moreover, (78) is calculated using (81), (82), and (83). Finally, substituting (75), , and into (66) and (9), respectively, the BEP expression for the piece-wise linear EH model is obtained.
Thanks to the mathematical tractability of our NL-EH models, we obtained closed-form expressions for the average BEP of the considered power-constraint system. Moreover, to gain a better insight through the new NL-EH models, we calculated BEP of the reference systems applying linear and piece-wise linear EH models. The detailed effects of the system parameters on the obtained BEP expressions are evaluated in the next section.
7. Performance Evaluation
In this section, simulation and theoretical results are obtained for different values of the system parameters. In all figures of this section, lines denote the theoretical results while symbols or symbols with dashed lines present the simulation results. In addition to the two proposed NL-EH models, the simulation results for linear, piece-wise linear, and existing NL-EH models are provided which gives an extra insight in analyzing the system performance. In all cases, the theoretical results are in perfect match with the Monte-Carlo simulation results. Theoretical curves for average harvested power are obtained from (14), (20), (21), and (22). Theoretical throughput curves are depicted by replacing (32), (42), (48), and (50) into (26). For the BER analysis, theoretical results are obtained by applying (54), (69), and (74). The fitting parameters for [11, 15] are listed in Tables 3 and 4, for the NL-EH-1 and NL-EH-2 models, respectively, while the coefficient of other models are provided in Table 5.
Unless stated differently, the fading parameters for each link are taken as , and the nodes S and D are both equipped with receive antennas. It is assumed that the modulation is binary phase-shift keying (BPSK) and the system throughput is obtained for bits/sec/Hz. Furthermore, it is assumed that and the path-loss coefficient is taken equal to [38]. Moreover, the channel power gains for both links are chosen as , and the noise power at node D is fixed to . Note that the optimization of the time slots dedicated for EH and IP is a difficult task due to the dependent parameters , , , , , , , , , , and .
Figure 4 describes the average harvested power at node S versus for data measurements taken from (14). The resulting curves are depicted for and . It is shown that the amount of average harvested power for is dBm which is approximately 25 dBm more compared to the case of . This verifies that the distance has a significant effect on the amount of harvested power. Moreover, it is seen that the average harvested powers of the NL-EH models for and are saturated at dBm and dBm, respectively. Please note that all NL-EH models are saturated at a predefined value. Furthermore, compared to dBm, more power is needed to saturate for .

The behavior of average harvested power versus is shown in Figure 5 considering the same system parameters for two different EH circuits given in [11, 15] for . It is observed that the EH circuit in [11] provides approximately 2 dBm more average harvested power compared to that in [15] considering the linear EH model. Moreover, it is shown that the average harvested power saturates at dBm and dBm for the EH circuits given in [11, 15], respectively, and results from the design of the EH system. This means that the EH circuit design has an important impact on the performance of the EH system since more power is harvested at the power-constraint node S. However, despite the EH circuit design, NL-EH models saturate at a predefined threshold value. Please note that all the NL-EH models along with the proposed NL-EH models have the same behavior on the amount of average harvested power.

In Figure 6, average harvested power is depicted versus . The results are provided for and dBm. It is concluded from Figure 6 that more antennas at node S result in higher average harvested power at the EH node which significantly affects the proposed system performance. However, it is shown that by increasing the number of receive antennas at node S, the average harvested power approximately saturates to the values between 38 dBm and 41 dBm for different NL-EH models when . In other words, node S with antennas provides an almost equal amount of harvested power compared to node S with antennas. Please note that the differences among NL-EH models are related to the curve fitting of the experimental measured data in a specified interval of input values. Moreover, it is seen from Figure 6 that the amount of average harvested power for the linear EH model is overestimated compared to the NL-EH and piece-wise linear EH models where for , the average harvested power for the linear EH model is 48.5 dBm. This means that the linear EH model overestimates approximately 8.5 dBm more average harvested power compared to NL-EH models.

In Figure 7, the system throughput is depicted versus . The results are provided for different numbers of receive antennas and . It is seen that by increasing the number of antennas at both nodes S and D, the throughput is increased. Besides, it is observed from Figure 7 that all of the throughput curves start from approximately , reach their maximum values, and then get decreased by increasing the value of until they get toward zero. Moreover, it is observed that the linear EH model overestimates the system performance for all cases. Additionally, the more receive antennas at nodes S and D, the more overestimation of throughput is resulted. In other words, the differences between NL-EH models and linear EH model throughput curves are increased by increasing and . This is due to the fact that no saturation power is defined in the linear EH model.

Figure 8 presents BER versus when . The curves are obtained for . It is shown that at high SNR values, the system performance considering the NL-EH and piece-wise linear EH models go to the error floor while the linear EH model system performance is increased at high SNR values due to the overestimation of harvested power which causes an unreasonable expectation of the system performance. However, at low SNR values, the BER performances of linear/piece-wise linear/NL-EH models are equivalent. Additionally, based on the obtained parameters from curve fitting of [15], piece-wise linear and logistic EH models provide approximately the same performance. Similarly, the heuristic model, fraction model, and proposed NL-EH-1 and NL-EH-2 models provide similar system performances. This comes from the fact that the difference in these results is caused by curve-fitting functions which are calculated based on a specific interval of data measurements.

The BER performance curves versus are shown in Figure 9. Here, the results are provided for dBm, , , and . It is seen that by increasing , the BER performance gets better for both values. Also, the linear EH model outperforms the piece-wise linear and NL-EH models for and due to more harvested power. It comes from the fact that more receive antennas at node S provide more harvested energy which results in higher transmit power at node S and consequently high received SNR at node D. Hence, it ensures a better system performance.

In Figure 10, the BER performance versus is given for dBm, and . It is observed that by increasing the number of antennas from to , the system performance gets significantly better. Additionally, by increasing the number of antennas at S, more power is harvested which results in better system performance for the linear EH model. However, for the NL-EH and piece-wise linear EH models, at high input power, the harvested power is saturated at a threshold power value, and then, even by increasing the number of antennas at S, namely, , the system performance goes to the error floor. Obviously, for both and , the same performances are obtained compared to and , respectively. In other words, the S node with only antenna provides approximately equivalent performance compared to . Considering the results from both Figures 9 and 10, it is concluded that the performances of piece-wise linear and NL-EH models are mostly dependent on the value of rather than since provides only an increase in the harvested power which is saturated at some value. However, even with small harvested power values but more antennas at node D, the system performance is improved. However, for the linear EH model, since the system performance is directly related to , increasing results in better system performance.

8. Conclusion
The performance of a wireless-powered communication system equipped with multiple antennas at both nodes S and D has been investigated based on the practical NL-EH models. Two new mathematically tractable NL-EH models have been proposed for which the PDFs and CDFs of the harvested power and their histograms have been provided considering the data measurements taken from the literature. Different from the NL-EH models given in the literature which are not analytically tractable, these new NL-EH models make the performance analyses analytically tractable. Furthermore, considering the Nakagami- fading and as a special case the Rayleigh fading, closed-form expressions for the average harvested powers, outage probabilities, throughputs, and BER performances have been analytically derived. Moreover, for comparison purposes, the same performance measures have been analytically derived for linear and piece-wise linear EH models. In each case, theoretical derivations have been supported by Monte-Carlo simulation results and have been compared with the existent EH models given in the literature such as linear, piece-wise linear, and other NL-EH models. It has been shown that apart from the linear EH model which overestimates the realistic system performance, the NL-EH models and piece-wise linear EH model provide practical system performances. Finally, the results obtained from different NL-EH models provide a comprehensive insight into the performance of the considered system. We have considered the case of one transmit antenna in this paper; however, it is possible to have more antennas for transmission at the source and to even analyze the system performance based on transmit antenna selection (TAS). However, these assumptions are out of scope of our paper and left as a future work.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the İTÜ Vodafone Future Lab under Project ITUVF20180701P09. The work of Mohammadreza Babaei was supported by the Research Fund of the Istanbul Technical University under Project MDK-2019-42438.