Performance Analysis of a User Selection Protocol in Cooperative Networks with Power Splitting Protocol-Based Energy Harvesting Over Nakagami-m / Rayleigh Channels

In this work, the system performance analysis of cooperative networks with power splitting protocol-based energy harvesting (EH) over Nakagami-m/Rayleigh channels is proposed. The exact-form expressions of the outage probability (OP) and ergodic capacity (EC) is demonstrated and derived. Using the proposed probabilistic models for wireless channels, we derive OP and EC as a research result. Finally, we conduct Monte Carlo simulations to verify a system performance analysis of the proposed system. The research results demonstrate the effectiveness of EH in the network over Nakagami-m/Rayleigh channels.


Introduction
Internet of Things (IoT) is considered a critical research area worldwide and has a significant impact on all activities in daily lives and industry [1][2][3].However, because of numerous disadvantages, some industrial and civil aspects of IoT are difficult to achieve.Energy limitations are the most significant problem in the long-term operation of wireless networks.Prolonging lifetimes and battery recharging/replacement are not practical solution for wireless networks.Wireless energy harvesting (EH) using radio frequency (RF) is considered a promising solution for prolonging the operating time of devices in wireless, energy-constrained cooperative networks [4][5][6].RF signals in EH wireless networks can transfer energy and information simultaneously.Cooperative networks are considered in connection with power sharing by helping the relay (R) between the source (S) and the destination (D) [6].In EH cooperative networks, an energy-constrained R, by simultaneously harvesting energy and transferring information, helps maintain constant operation and connection without the use of external energy sources.The authors in [7] considered a dual-hop cognitive, inter-vehicular, relay-assisted communication system where all communication links are not line of sight, and their fading is modeled by the binary Rayleigh fading distribution.An analytic performance evaluation of the bit error rate (BER) of underlay decode-and-forward cognitive networks with the best relay selection over Rayleigh multipath fading channels was introduced in [8].The authors in [9,10] presented a general framework for modeling and evaluating the performance of dual-hop, decode-and-forward (DF) relaying schemes over independent and not necessarily identically distributed (INID) Nakagami-m fading channels.The performance of dual-hop MIMO OSTBC transmission with multiple antennas employed at the transceiver and CSI-assisted AF relay over Nakagami-m fading is proposed in [11], and the performance of dual-hop, decode-and-forward relaying with relay selection (RS), analyzed over Nakagami fading channels, was studied in [12].In our previous papers, we consider the EH over Rician fading channels.The titles of these papers are as follows: A Performance Analysis for Half-Duplex Bidirectional Sensor Networks under Hardware Impairments [12], Adaptive EH Relaying Protocol for a Two-Way Half Duplex System Network over Rician Fading Channels [13], Energy Harvesting-based Spectrum Access with Incremental Cooperation, Relay Selection, and Hardware Noises [14], and Performance Enhancement for EH Based Two-way Relay Protocols in Wireless Ad-Hoc Networks with Partial and Full Relay Selection Methods [15].
The problem of a two-hop EH cooperative relaying network was studied in [16].In [17], a harvest-then-cooperate (HTC) protocol with AP points for S and R EH was investigated.The authors in [18] maximized the end-to-end achievable information transmission rate by assuming perfect channel state information (CSI) of the cooperative network with time switching (TS) and power splitting (PS) protocols.The authors in [19] investigated two-way relay beamforming optimization to maximize the achievable sum rate of simultaneous wireless information and power transfer (SWIPT) system with a full-duplex (FD) multiple-input multiple-output (MIMO), amplify-and-forward (AF) relay.To extend [19], the authors in [20] studied the joint optimization of two-way relay beamforming.Furthermore, the authors in [21] investigated SWIPT mechanisms in a relaying network with finite block length (FBL) codes.Moreover, the authors in [22] derived the theoretical symbol error probability (SEP) of cooperative systems with the best relay selection for Nakagami-m fading channels, and the authors in [23] proposed a novel DSTBC scheme with embedded adaptive DAF/amplify-and-forward (AAF) elements for data retransmission in multihop PLC networks with existing opportunistic listening (OL) protocols.The authors in [24] proposed a simple adaptive relaying protocol (ARP) for general relay networks, and in [25] we proposed a full rate DQOSTBC matrix with embedded adaptive DAF/AAF elements for four single-antenna relays.In spite of this research, cooperative relaying networks still need to be studied.This is the main aim of this paper.
There is very little research focus on both user selection protocol and EH with different Nakagami-m/Rayleigh channels.For this purpose, we analyzed the performance of user selection protocol cooperative networks with PS protocol-based EH over Nakagami-m/Rayleigh channels.Our article provides the following contributions.(i) We derive closed-form expressions of Outage probability (OP) and ergodic capacity (EC), (ii) propose user selection for choosing the best-received destination, (iii) investigate and compare the OP and EC with maximum and non-maximum capacity, and (iv) verify all analytical expressions by Monte Carlo simulations.The structure of the rest of the paper is as follows.Section 2 illustrates the system model, and Section 3 investigates the system performance.Section 4 proposes and discusses the research results, and conclusions are proposed in Section 5.

System Model
We consider an EH relay network with one S, multiple Ds, and one R as shown in Figure 1 [13,14].In this system model, S, R, and D are working in a half-duplex (HD) mode with one antenna.The channel gain between node S and the relay R (h RD ) is a Nakagami-m fading channel, so between the relay R and the destinations, D i as h RDi is represented by Rayleigh fading channels.In this model, the direct link between S and D nodes is too weak without the help of a relay.The EH and information transmission (IT) for this proposed model system is presented in Figure 2. In this model, the transmission length time T is divided into two slots.In the time slot T/2, R harvests energy in ρP s and receives information in (1−ρ)P s from S. The remaining half-time slot T/2 is used for information transferring from R to D as in [15,16].In the first transmission phase, the received signal at R in period T/2 time can be formulated by the following:

Energy Harvesting Information transmission
In this equation, hsr is S to R channel gain, dsr is the distance between S and R, and m denotes the path loss exponent.Here, s x is the transmitted signal at S, nr is the additive white Gaussian noise (AWGN) with variance N0, and 0 1 ρ < < is the PS ratio at the relay R.Moreover, { } is the expectation operator, and Ps is the average transmit power at S.
After the first-time slot T/2, the harvested power at R could be obtained by where 0 1 η < ≤ is the energy conversion efficiency.
In the second time slot T/2, R retransmits the information from S to Di.The received signal at the n th destination at the second slot time can be expressed as where , hrdi is the R to the i th D channel gain, di is the R to the D distance, ndi is the additive white Gaussian noise (AWGN) with variance N0, and { } Here, the AF protocol is considered, and the amplifying factor can be given by   In the first transmission phase, the received signal at R in period T/2 time can be formulated by the following:

Energy Harvesting Information transmission
In this equation, hsr is S to R channel gain, dsr is the distance between S and R, and m denotes the path loss exponent.Here, s x is the transmitted signal at S, nr is the additive white Gaussian noise (AWGN) with variance N0, and 0 1 ρ < < is the PS ratio at the relay R.Moreover, { } is the expectation operator, and Ps is the average transmit power at S.
After the first-time slot T/2, the harvested power at R could be obtained by where 0 1 η < ≤ is the energy conversion efficiency.
In the second time slot T/2, R retransmits the information from S to Di.The received signal at the n th destination at the second slot time can be expressed as where , hrdi is the R to the i th D channel gain, di is the R to the D distance, ndi is the additive white Gaussian noise (AWGN) with variance N0, and { } Here, the AF protocol is considered, and the amplifying factor can be given by In the first transmission phase, the received signal at R in period T/2 time can be formulated by the following: In this equation, h sr is S to R channel gain, d sr is the distance between S and R, and m denotes the path loss exponent.Here, x s is the transmitted signal at S, n r is the additive white Gaussian noise (AWGN) with variance N 0 , and 0 < ρ < 1 is the PS ratio at the relay R.Moreover, E |x s | 2 = P s , E{•} is the expectation operator, and P s is the average transmit power at S.
After the first-time slot T/2, the harvested power at R could be obtained by where 0 < η ≤ 1 is the energy conversion efficiency.
In the second time slot T/2, R retransmits the information from S to D i .The received signal at the n th destination at the second slot time can be expressed as where i ∈ (1, 2, . . ., K), h rdi is the R to the i th D channel gain, d i is the R to the D distance, n di is the additive white Gaussian noise (AWGN) with variance N 0 , and E |x r | 2 = P r .
Here, the AF protocol is considered, and the amplifying factor can be given by From Equations ( 1), ( 3) and ( 4), the received signal can be rewritten as Therefore, the overall signal to noise ratio (SNR) from S to D can be given by After doing some algebra, using the fact that N 0 << P r and denote Equation ( 6) can be rewritten as Combined with Equation ( 2), we have the SNR: where ψ = P s /N 0 .
In our proposed model, we assume that the S-R link belongs to a Nakagami-m fading channel and the R-D i link belongs to a Rayleigh fading channel.In the next section, we derive the analytical expression of OP and EC as in [26][27][28][29][30][31].

The System Performance
As in previous work [12,13,28,30], the closed-form expressions of OP and EC of this system model are derived in this section.At first, we will determine the probability density function (PDF) and the cumulative density function (CDF) of a random variable (RV) |h sr | 2 , h rd i 2 .As shown in [29], the PDF of RV γ 1 can be calculated by From Equation ( 9), the CDF of RV |h sr | 2 can be obtained with the help of Equation (8.353.4) in [32].
where Furthermore, we have the PDF and CDF of RV h rd i 2 as shown in [28] as follows: where is the mean of RV h rd i 2 .
In this system model, the best selection user US q is analyzed as follows: According to [19], the CDF of |h rd | 2 can be given by where we denote , and λ rd is the mean of RV |h rd | 2 .Finally, the PDF of random variable can be determined by From Equation ( 19), we have

The Outage Probability (OP)
Theorem 1 (OP-Closed Form).The closed-form expression of the OP of the proposed model system can be formulated as Proof of Theorem 1. See Appendix A.

Maximize Capacity
In this section, we will find the PS factor to maximize capacity.The overall capacity from source to destination can be given by Therefore, to maximize capacity, from Equation ( 18), we can see that max(C s,d ) ⇔ max(SNR) .
From Equation ( 8), SNR is a concave function of ρ, which, maximizing the SNR, can be given as follows: Theorem 2 (EC-Closed Form).The closed-form expression of the EC in the maximize mode of the system model can be expressed as Proof of Theorem 2. See Appendix B.

Numerical Results and Discussion
Like in previous studies [33][34][35][36], we conducted Monte Carlo simulations to verify analytical expressions of the OP and EC of the proposed system in the above section.In addition, we investigated the effect of the primary system parameter on the system performance in terms of OP and EC.All other simulation parameters are listed in Table 1. Figure 3a,b plots the influence of P s /N 0 on the OP and EC of the proposed system.In Figure 3, the main parameters are as follows: K = 2, R = 0.5, and ρ = 0.2 and 0.6.In this figure, the OP of the case ρ = 0.2 and 06 and the maximum capacity are also proposed for comparison.It is observed that the simulation values of the OP match the values from the mathematical analysis.In connection with the effect of ρ, the OP decreases, and EC increases as ρ varies from 0.2 to 0.6.When P s /N 0 increases from 0 to 20 dB, the OP decreases and EC significantly increases.Furthermore, the higher the value of ρ is, the faster the OP decreases and the EC increases.In addition, we can see that the OP and EC of the model system in the maximum capacity case are better in comparison with the other cases, with other values of ρ.This can be observed based on the mathematical analysis in Equations ( 17) and (20).Figure 4a,b depicts the effect of Ps/N0 on the OP and the EC.We set R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6 in Figure 4a and R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6, respectively.From Figure 4a, the OP decreases when Ps/N0 increases from 0 to 20 dB, and OP decreases faster with a higher K. On the other hand, the EC increases significantly, while Ps/N0 rises from 0 to 20 dB.Furthermore, the EC is higher with the higher K value.In all research results, the simulation and analytical results are the same.Moreover, Figure 5a,b plots the influence of the ratio Ps/N0 on the OP and the EC of the model system.Here, the cases Ray-Ray, Naka-Ray in the non-maximize and maximize modes are compared with each other in the same system condition.In the simulation, we set the ratio Ps/N0 increased from 0 to 20 dB, ρ = 0.2, and K = 2 for the OP and ρ = 0.4 and K = 3 for the EC, respectively.Figure 5a shows that the OP decreases faster in the Naka-Ray case with maximum capacity compared with other cases.In the same way, the EC increases faster in the Ray-Ray case with non-maximum capacity in Figure 5b.Here we can see that the system performance in the maximum capacity case is better than in the non-maximum capacity case.Furthermore, the simulation results agreed with the mathematical analysis of the above section.Figure 4a,b depicts the effect of P s /N 0 on the OP and the EC.We set R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6 in Figure 4a and R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6, respectively.From Figure 4a, the OP decreases when P s /N 0 increases from 0 to 20 dB, and OP decreases faster with a higher K. On the other hand, the EC increases significantly, while P s /N 0 rises from 0 to 20 dB.Furthermore, the EC is higher with the higher K value.In all research results, the simulation and analytical results are the same.Figure 4a,b depicts the effect of Ps/N0 on the OP and the EC.We set R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6 in Figure 4a and R = 0.5 bps, ρ = 0.5, and K = 1, 3, and 6, respectively.From Figure 4a, the OP decreases when Ps/N0 increases from 0 to 20 dB, and OP decreases faster with a higher K. On the other hand, the EC increases significantly, while Ps/N0 rises from 0 to 20 dB.Furthermore, the EC is higher with the higher K value.In all research results, the simulation and analytical results are the same.Moreover, Figure 5a,b plots the influence of the ratio Ps/N0 on the OP and the EC of the model system.Here, the cases Ray-Ray, Naka-Ray in the non-maximize and maximize modes are compared with each other in the same system condition.In the simulation, we set the ratio Ps/N0 increased from 0 to 20 dB, ρ = 0.2, and K = 2 for the OP and ρ = 0.4 and K = 3 for the EC, respectively.Figure 5a shows that the OP decreases faster in the Naka-Ray case with maximum capacity compared with other cases.In the same way, the EC increases faster in the Ray-Ray case with non-maximum capacity in Figure 5b.Here we can see that the system performance in the maximum capacity case is better than in the non-maximum capacity case.Furthermore, the simulation results agreed with the mathematical analysis of the above section.Moreover, Figure 5a,b plots the influence of the ratio P s /N 0 on the OP and the EC of the model system.Here, the cases Ray-Ray, Naka-Ray in the non-maximize and maximize modes are compared with each other in the same system condition.In the simulation, we set the ratio P s /N 0 increased from 0 to 20 dB, ρ = 0.2, and K = 2 for the OP and ρ = 0.4 and K = 3 for the EC, respectively.Figure 5a shows that the OP decreases faster in the Naka-Ray case with maximum capacity compared with other cases.In the same way, the EC increases faster in the Ray-Ray case with non-maximum capacity in Figure 5b.Here we can see that the system performance in the maximum capacity case is better than in the non-maximum capacity case.Furthermore, the simulation results agreed with the mathematical analysis of the above section.Figure 6 illustrates the OP versus the system factors η.In this case, we set Ps/N0 = 10 dB, ρ = 0.3, and K = 1 and 4 in each case.From Figure 6, the OP decreases when η increases from 0 to 1, and the OP in the maximum capacity case is better than the non-maximum capacity case.Here, the simulation results agreed with the mathematical analysis of the above section.Furthermore, the OP and the EC of the proposed system in connection with the number of users are presented in Figure 7a,b.Similarly, the OP decreases and the EC increases, while the number of users varies from 0 to 10 and have better values in the maximum capacity case.In all of them, the simulation and analytical mathematical results agreed well with each other.Figure 6 illustrates the OP versus the system factors η.In this case, we set P s /N 0 = 10 dB, ρ = 0.3, and K = 1 and 4 in each case.From Figure 6, the OP decreases when η increases from 0 to 1, and the OP in the maximum capacity case is better than the non-maximum capacity case.Here, the simulation results agreed with the mathematical analysis of the above section.Furthermore, the OP and the EC of the proposed system in connection with the number of users are presented in Figure 7a,b.Similarly, the OP decreases and the EC increases, while the number of users varies from 0 to 10 and have better values in the maximum capacity case.In all of them, the simulation and analytical mathematical results agreed well with each other.Figure 6 illustrates the OP versus the system factors η.In this case, we set Ps/N0 = 10 dB, ρ = 0.3, and K = 1 and 4 in each case.From Figure 6, the OP decreases when η increases from 0 to 1, and the OP in the maximum capacity case is better than the non-maximum capacity case.Here, the simulation results agreed with the mathematical analysis of the above section.Furthermore, the OP and the EC of the proposed system in connection with the number of users are presented in Figure 7a,b.Similarly, the OP decreases and the EC increases, while the number of users varies from 0 to 10 and have better values in the maximum capacity case.In all of them, the simulation and analytical mathematical results agreed well with each other.

Conclusions
This study proposes a performance analysis of a user selection protocol in cooperative networks with PS protocol-based EH over Nakagami-m/Rayleigh channels.The analytical closed-form expressions for the OP and the EC is derived.Moreover, the closed-form expression of the EC is derived in cases of maximum and non-maximum capacity.Finally, we use Monte Carlo simulations to derive analytical expressions in connection with all possible system parameters.From the research results, we can see that the simulation and the analytical analysis agree well.

Appendix A
Outage probability can be calculated as where − is a threshold, and R is the source rate.Equation (A1) can be rewritten as ( ) From Equations ( 19) and ( 20), we have

Conclusions
This study proposes a performance analysis of a user selection protocol in cooperative networks with PS protocol-based EH over Nakagami-m/Rayleigh channels.The analytical closed-form expressions for the OP and the EC is derived.Moreover, the closed-form expression of the EC is derived in cases of maximum and non-maximum capacity.Finally, we use Monte Carlo simulations to derive analytical expressions in connection with all possible system parameters.From the research results, we can see that the simulation and the analytical analysis agree well.

Appendix A
Outage probability can be calculated as where z = 2 2R − 1 is a threshold, and R is the source rate.Equation (A1) can be rewritten as From Equations ( 19) and ( 20), we have (A6) Using the Table of Integral Equation (3.471,9) in [32], Equation (A7) can reformulated as where K v (•) is the modified Bessel function of the second kind and the v-th order.

Appendix B
It is easy to observe that ∂ 2 SNR ∂ρ 2 is negative for 0 < ρ < 1.Hence, we conclude that SNR is a concave function of ρ (0 < ρ < 1).We can find the value of ρ that maximizes SNR by differentiating the SNR concerning ρ, and we can then equate it to zero.After some algebraic calculations, we have the following possible solutions for ρ * : If Equation (A3) is placed into Equation (A1), SNR max can be obtained as EC analysis can be demonstrated as follows: As in the previous section, we have In Equation (A14), we consider We then apply a Taylor series as follows: exp We can observe that the involving integral in Equation (A22) is difficult to solve in a closed form.However, by changing the variable of the integration in Equation (B13) as z = tan θ, we have 1 + tan θ sec 2 θdθ.(A23)

Figure 2 .
Figure 2. EH and IT of the proposed system.

Figure 2 .
Figure 2. EH and IT of the proposed system.

Figure 2 .
Figure 2. EH and IT of the proposed system.

Figure 3 .
Figure 3.The outage probability (a) and ergodic capacity (b) versus P s /N 0 while K = 2.

Figure 5 .
Figure 5.Comparison of the outage probability (a) and ergodic capacity (b) ersus P s /N 0 .

Figure 7 .
Figure 7.The outage probability (a) and ergodic capacity (b) versus the numbers of users.
by applying the equation (x + y) m = m n=0 m n x m−n y n to Equation (A5), the outage probability can demonstrate as follows: value of ρ * > 1 or ρ * < 0, we choose ρ * =