Next Article in Journal
A Capacitor Voltage Balancing Approach Based on Mapping Strategy for MMC Applications
Previous Article in Journal
Rolling 3D Laplacian Pyramid Video Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Faculty of Electrical Engineering and Computer Science, VSB-Technical University of Ostrava 17, listopadu 15/2172, 708 33 Ostrava-Poruba, Czech Republic
2
Optoelectronics Research Group, Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City 70000, Vietnam
3
Center for Information Technology, Ho Chi Minh City University of Food Industry, Ho Chi Minh City 70000, Vietnam
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(4), 448; https://doi.org/10.3390/electronics8040448
Submission received: 27 February 2019 / Revised: 7 April 2019 / Accepted: 15 April 2019 / Published: 19 April 2019
(This article belongs to the Section Networks)

Abstract

:
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.

1. 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.

2. 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 (hRD) is a Nakagami-m fading channel, so between the relay R and the destinations, Di as hRDi 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 ρPs and receives information in (1−ρ)Ps 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:
y r = 1 d s r m ( 1 ρ ) h s r x s + n r .
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, x s 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, E { | x s | 2 } = P s , E { } 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
P r = E h d s r m ( T / 2 ) = η ρ P s | h s r | 2 ( T / 2 ) d s r m ( T / 2 ) = η ρ P s | h s r | 2 d s r m
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 nth destination at the second slot time can be expressed as
y d i = 1 d i m h r d i x r + n d i
where i ( 1 , 2 , , K ) , hrdi is the R to the ith D channel gain, di is the R to the D distance, ndi is the additive white Gaussian noise (AWGN) with variance N0, and E { | x r | 2 } = P r .
Here, the AF protocol is considered, and the amplifying factor can be given by
β = x r y r = P r ( 1 ρ ) P s | h s r | 2 d s r m + N 0 .
From Equations (1), (3) and (4), the received signal can be rewritten as
y d i = 1 d r d n m h r d i β y r + n d i = 1 d r d i m h r d i β [ 1 d s r m ( 1 ρ ) h s r x s + n r ] + n d i = ( 1 ρ ) h s r x s h r d i β d i m d s r m s i g n a l + 1 d i m h r d i β n r + n d i n o i s e
Therefore, the overall signal to noise ratio (SNR) from S to D can be given by
S N R = E { | s i g n a l | 2 } E { | n o i s e | 2 } = ( 1 ρ ) P s | h s r | 2 | h r d i | 2 β 2 d i m d s r m | h r d i | 2 β 2 N 0 d i m + N 0 = ( 1 ρ ) P s | h s r | 2 | h r d i | 2 | h r d i | 2 d s r m N 0 + N 0 d i m d s r m β 2 .
After doing some algebra, using the fact that N0 << Pr and denote γ 1 = | h s r | 2 , γ i = | h r d i | 2 , Equation (6) can be rewritten as
S N R = ( 1 ρ ) P s P r | h s r | 2 | h r d i | 2 | h r d i | 2 d s r m P r N 0 + ( 1 ρ ) P s | h s r | 2 N 0 d i m .
Combined with Equation (2), we have the SNR:
S N R = η ρ ( 1 ρ ) ψ | h s r | 2 | h r d i | 2 η ρ | h r d i | 2 d s r m + ( 1 ρ ) d i m
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-Di 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].

3. 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 s r | 2 , | h r d i | 2 . As shown in [29], the PDF of RV γ 1 can be calculated by
f | h s r | 2 ( x ) = x m | h s r | 2 1 ( m | h s r | 2 1 ) ! ( Ω | h s r | 2 ) m | h s r | 2 exp ( x Ω | h s r | 2 ) .
From Equation (9), the CDF of RV | h s r | 2 can be obtained with the help of Equation (8.353.4) in [32].
F | h s r | 2 ( x ) = 1 exp ( x Ω | h s r | 2 ) t = 0 m | h s r | 2 1 x t t ! ( Ω | h s r | 2 ) t
where Ω | h s r | 2 = λ s r m | h s r | 2 ; m | h s r | 2 is the Nakagami-m parameter and note that the case of m | h s r | 2 = 1 corresponds to Rayleigh fading; λ s r = E { | h s r | 2 } d s r m is the mean of RV | h s r | 2 .
Furthermore, we have the PDF and CDF of RV | h r d i | 2 as shown in [28] as follows:
f | h r d i | 2 ( x ) = 1 λ r d i e x / λ r d i
F | h r d i | 2 ( x ) = 1 e x / λ r d i
where λ r d i = E { | h r d i | 2 } d i m is the mean of RV | h r d i | 2 .
In this system model, the best selection user USq is analyzed as follows:
| h r d | 2 = max i = 1 , 2 , , K ( | h r d i | 2 ) .
According to [19], the CDF of | h r d | 2 can be given by
F | h r d | 2 ( y ) = p = 0 K ( 1 ) p C K p × e p y / λ r d ,
where we denote C K p = K ! p ! ( K p ) ! , and λ r d is the mean of RV | h r d | 2 .
Finally, the PDF of random variable can be determined by
f | h r d | 2 ( y ) = 1 λ r d p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) y / λ r d .
From Equation (19), we have
S N R = η ρ ( 1 ρ ) ψ | h s r | 2 | h r d | 2 η ρ | h r d | 2 d s r m + ( 1 ρ ) d i m .

3.1. 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
P o u t = 1 2 exp [ z d s r m Ω | h s r | 2 ( 1 ρ ) ψ ] t = 0 m γ 1 1 n = 0 t p = 0 K 1 ( 1 ) p K × C K 1 p z n + t + 1 2 ( λ r d ) n t 1 2 ( t n ) ! n ! ( Ω | h s r | 2 ψ ) n + t + 1 2 ( p + 1 ) n t + 1 2 ( d i m η ρ ) t n + 1 2 [ d s r m ( 1 ρ ) ] n × K n t + 1 ( 2 z d i m ( p + 1 ) Ω | h s r | 2 η ρ ψ λ r d )
Proof of Theorem 1.
See Appendix A. □

3.2. 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
C s , d = 1 2 log 2 ( 1 + S N R ) = 1 2 log 2 [ 1 + η ρ ( 1 ρ ) ψ | h s r | 2 | h r d | 2 η ρ | h r d | 2 d s r m + ( 1 ρ ) d i m ] .
Therefore, to maximize capacity, from Equation (18), we can see that max ( C s , d ) max ( S N R ) .
From Equation (8), SNR is a concave function of ρ, which, maximizing the SNR, can be given as follows:
ρ * = 1 1 + | h r d | η d s r m d i m .
Theorem 2 (EC—Closed Form).
The closed-form expression of the EC in the maximize mode of the system model can be expressed as
C A F = 1 ln 2 v = 1 N p ω v 1 F S N R max ( x v ) 1 + x v .
Proof of Theorem 2.
See Appendix B. □

4. 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 Ps/N0 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 Ps/N0 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 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.

5. 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.

Author Contributions

T.N.N. and M.T. created the main idea and developed the performance evaluation by analysis and simulation. M.T. wrote and revised the paper. T.-L.N. and D.-H.H. joined in the data analysis process. M.V. worked as an advisor to the group author of this project.

Funding

This research was funded by the Czech Ministry of Education, Youth and Sports within the grants No. SP2018/59 and SP2019/41 conducted at the VSB—Technical University of Ostrava.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Outage probability can be calculated as
P o u t = Pr ( S N R < z ) = Pr ( η ρ ( 1 ρ ) ψ | h s r | 2 | h r d | 2 η ρ | h r d | 2 d s r m + ( 1 ρ ) d i m < z )
where z = 2 2 R 1 is a threshold, and R is the source rate.
Equation (A1) can be rewritten as
P o u t = Pr { η ρ ( 1 ρ ) ψ | h s r | 2 | h r d | 2 < z η ρ | h r d | 2 d s r m + z ( 1 ρ ) d i m } = Pr [ | h s r | 2 < z η ρ | h r d | 2 d s r m + z ( 1 ρ ) d i m η ρ ( 1 ρ ) ψ | h r d | 2 ] = Pr [ | h s r | 2 < z d i m η ρ ψ | h r d | 2 + z d s r m ( 1 ρ ) ψ ]
P o u t = 0 F | h s r | 2 { z d i m η ρ ψ | h r d | 2 + z d s r m ( 1 ρ ) ψ | | h r d | 2 } f | h r d | 2 ( | h r d | 2 ) d ( | h r d | 2 ) .
From Equations (19) and (20), we have
P o u t = 1 0 exp [ z d i m η ρ ψ | h r d | 2 + z d s r m ( 1 ρ ) ψ Ω | h s r | 2 ] t = 0 m γ 1 1 [ z d i m η ρ ψ | h r d | 2 + z d s r m ( 1 ρ ) ψ ] t t ! ( Ω | h s r | 2 ) t × 1 λ q p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ q d ( | h r d | 2 )
P o u t = 1 exp [ z d s r m Ω | h s r | 2 ( 1 ρ ) ψ ] λ r d 0 t = 0 m γ 1 1 [ z d i m η ρ ψ | h r d | 2 + z d s r m ( 1 ρ ) ψ ] t t ! ( Ω | h s r | 2 ) t × exp [ z d i m Ω | h s r | 2 η ρ ψ | h r d | 2 ] × p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
Now by applying the equation ( x + y ) m = n = 0 m ( m n ) x m n y n to Equation (A5), the outage probability can demonstrate as follows:
P o u t = 1 exp [ z d s r m Ω | h s r | 2 ( 1 ρ ) ψ ] λ r d 0 t = 0 m γ 1 1 n = 0 t ( t n ) z t t ! ( Ω | h s r | 2 ) t ( d i m η ρ ψ | h r d | 2 ) t n [ d s r m ( 1 ρ ) ψ ] n × exp [ z d i m Ω | h s r | 2 η ρ ψ | h r d | 2 ] × p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
P o u t = 1 exp [ z d s r m Ω | h s r | 2 ( 1 ρ ) ψ ] λ r d t = 0 m γ 1 1 n = 0 t p = 0 K 1 ( 1 ) p K × C K 1 p z t ( t n ) ! n ! ( Ω | h s r | 2 ) t ψ t ( d i m η ρ ) t n [ d s r m ( 1 ρ ) ] n × 0 ( | h r d | 2 ) n t × exp [ z d i m Ω γ 1 η ρ ψ | h r d | 2 ] × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
Using the Table of Integral Equation (3.471,9) in [32], Equation (A7) can reformulated as
P o u t = 1 2 exp [ z d s r m Ω | h s r | 2 ( 1 ρ ) ψ ] λ r d t = 0 m γ 1 1 n = 0 t p = 0 K 1 ( 1 ) p K × C K 1 p z t ( t n ) ! n ! ( Ω | h s r | 2 ) t ψ t ( d i m η ρ ) t n [ d s r m ( 1 ρ ) ] n ( z d i m λ r d Ω | h s r | 2 η ρ ψ ( p + 1 ) ) n t + 1 2 × K n t + 1 ( 2 z d i m ( p + 1 ) Ω | h s r | 2 η ρ ψ λ r d )
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 S N R ρ 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 ρ*:
ρ * = 1 1 + | h r d | η d s r m d i m   or   ρ * = 1 1 | h r d | η d s r m d i m .
Since ρ * = 1 1 | h r d | η d s r m d i m results in a value of ρ* > 1 or ρ* < 0, we choose ρ * = 1 1 + | h r d | η d s r m d i m as the solution.
If Equation (A3) is placed into Equation (A1), SNRmax can be obtained as
S N R max = η ψ | h s r | 2 | h r d | 2 ( 1 + | h r d | η d s r m d i m ) 2 d i m .
EC analysis can be demonstrated as follows:
C A F = 0 f S N R max ( z ) log z ( 1 + z ) d z = 1 ln 2 0 1 F S N R max ( z ) 1 + z d z
F S N R max ( z ) = Pr ( S N R max < z ) = Pr [ η ψ | h s r | 2 | h r d | 2 ( 1 + | h r d | η d s r m d i m ) 2 d i m < z ] = Pr [ | h s r | 2 < ( 1 + | h r d | η d s r m d i m ) 2 z d i m η ψ | h r d | 2 ]
F S N R max ( z ) = 0 F | h s r | 2 [ ( 1 + | h r d | η d s r m d i m ) 2 z d i m η ψ | h r d | 2 | | h r d | 2 ] f | h r d | 2 ( | h r d | 2 ) d ( | h r d | 2 ) .
As in the previous section, we have
F S N R max ( z ) = 1 1 λ r d 0 exp [ ( 1 + | h r d | η d s r m d i m ) 2 z d i m Ω | h s r | 2 η ψ | h r d | 2 ] t = 0 m | h s r | 2 1 1 t ! ( Ω | h s r | 2 ) t × [ ( 1 + | h r d | η d s r m d i m ) 2 z d i m η ψ | h r d | 2 ] t × p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
In Equation (A14), we consider
A = exp [ ( 1 + | h r d | η d s r m d i m ) 2 z d i m Ω | h s r | 2 η ψ | h r d | 2 ] = exp [ z d i m + 2 z | h r d | η d s r m d i m + z | h r d | 2 η d s r m Ω | h s r | 2 η ψ | h r d | 2 ] = exp [ z d i m Ω | h s r | 2 η ψ | h r d | 2 ] exp [ 2 z d s r m d i m Ω | h s r | 2 ψ | h r d | η ] exp [ z d s r m Ω | h s r | 2 ψ ]
We then apply a Taylor series as follows:
exp [ 2 z d s r m d i m Ω | h s r | 2 ψ | h r d | η ] = q = 0 ( 2 z d s r m d i m Ω | h s r | 2 ψ | h r d | η ) q q ! = q = 0 ( 1 ) q ( 2 z d s r m d i m Ω | h s r | 2 ψ η ) q q ! | h r d | q .
Equation (A15) can then be rewritten as
A = q = 0 ( 1 ) q ( 2 z d s r m d i m Ω | h s r | 2 ψ η ) q q ! | h r d | q × exp [ z d s r m Ω | h s r | 2 ψ ] exp [ z d i m Ω | h s r | 2 η ψ | h r d | 2 ] .
If Equation (A14) is placed into Equation (A16) and the expression ( x + y ) m = n = 0 m ( m n ) x m n y n is applied, we have
F S N R max ( z ) = 1 1 λ r d 0 q = 0 ( 1 ) q ( 2 z d s r m d i m Ω | h s r | 2 ψ η ) q q ! | h r d | q × exp [ z d s r m Ω | h s r | 2 ψ ] exp [ z d i m Ω | h s r | 2 η ψ | h r d | 2 ] × t = 0 m | h s r | 2 1 ( z d i m ) t t ! ( Ω | h s r | 2 η ψ ) t | h r d | 2 t × n = 0 2 t ( 2 t n ) [ | h r d | η d s r m d i m ] n × p = 0 K 1 ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
F S N R max ( z ) = 1 1 λ r d t = 0 m | h s r | 2 1 n = 0 2 t p = 0 K 1 q = 0 0 exp [ z d s r m Ω | h s r | 2 ψ ] exp [ z d i m Ω | h s r | 2 η ψ | h r d | 2 ] × ( 2 ) q ( 2 t ) ! ( z ) t + q | h r d | n 2 t q ( d s r m ) q 2 ( d i m ) q 2 n ! t ! ( 2 t n ) ! q ! ( Ω | h s r | 2 ψ ) t + q η t + q / 2 × ( η d s r m d i m ) n 2 × ( 1 ) p C K 1 p K × e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
F S N R max ( z ) = 1 e z d s r m Ω | h s r | 2 ψ λ r d t = 0 m | h s r | 2 1 n = 0 2 t p = 0 K 1 q = 0 ( 1 ) p C K 1 p K × ( 2 ) q ( 2 t ) ! ( z ) t + q ( d s r m ) q + n 2 ( d i m ) q n 2 n ! t ! ( 2 t n ) ! q ! ( Ω | h s r | 2 ψ ) t + q η t + q / 2 n / 2 × 0 | h r d | n 2 t q × exp [ z d i m Ω | h s r | 2 η ψ | h r d | 2 ] e ( p + 1 ) | h r d | 2 / λ r d d ( | h r d | 2 )
Using the Table of Integral Equation (3.471,9) in [32], Equation (A18) can be reformulated as
F S N R max ( z ) = 1 2 e z d s r m Ω | h s r | 2 ψ λ r d t = 0 m | h s r | 2 1 n = 0 2 t p = 0 K 1 q = 0 ( 1 ) p C K 1 p K × ( 2 ) q ( 2 t ) ! ( z ) t + q ( d s r m ) q + n 2 ( d i m ) q n 2 n ! t ! ( 2 t n ) ! q ! ( Ω | h s r | 2 ψ ) t + q η t + q / 2 n / 2 × [ z d i m λ r d ( p + 1 ) Ω | h s r | 2 η ψ ] n 2 t q + 2 4 × K n 2 t q + 2 2 [ 2 z d i m ( p + 1 ) Ω | h s r | 2 η ψ λ r d ]
F S N R max ( z ) = 1 2 e z d s r m Ω | h s r | 2 ψ × t = 0 m | h s r | 2 1 n = 0 2 t p = 0 K 1 q = 0 ( 1 ) p C K 1 p K × ( 2 ) q ( 2 t ) ! ( d s r m ) q + n 2 ( d i m ) q n 2 t + 2 4 ( λ r d ) n 2 t q 2 4 n ! t ! ( 2 t n ) ! q ! η 2 t n + q + 2 4 ( p + 1 ) n 2 t q + 2 4 × ( z Ω | h s r | 2 ψ ) n + 2 t + 3 q + 2 4 × K n 2 t q + 2 2 [ 2 z d i m ( p + 1 ) Ω | h s r | 2 η ψ λ r d ]
C A F = 1 ln 2 0 1 F S N R max ( z ) 1 + z d z .
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
C A F = 1 ln 2 0 1 F S N R max ( z ) 1 + z d z = 1 ln 2 0 π / 2 1 F S N R max ( tan θ ) 1 + tan θ sec 2 θ d θ .
We can then apply an efficient NP-point Gauss–Chebychev quadrature (GCQ) formula ([19] Equation (25.4.39)) from [20] to numerically derive the Proof of Theorem 2.
Where x v = tan ( π 4 cos [ 2 v 1 2 N p π ] + π 4 ) , and ω v = π 2 sin ( 2 v 1 2 N p π ) 4 N p cos 2 ( π 4 cos [ 2 v 1 2 N p π ] + π 4 ) .

References

  1. Bi, S.; Ho, C.K.; Zhang, R. Wireless powered communication: Opportunities and challenges. IEEE Commun. Mag. 2015, 53, 117–125. [Google Scholar] [CrossRef]
  2. Niyato, D.; Kim, D.I.; Maso, M.; Han, Z. Wireless Powered Communication Networks: Research Directions and Technological Approaches. IEEE Wirel. Commun. 2017, 2–11. [Google Scholar] [CrossRef]
  3. Yu, H.; Lee, H.; Jeon, H. What is 5G? Emerging 5G Mobile Services and Network Requirements. Sustainability 2017, 9, 1848. [Google Scholar] [CrossRef]
  4. Zhou, X.; Zhang, R.; Ho, C.K. Wireless Information and Power Transfer: Architecture Design and Rate-Energy Tradeoff. IEEE Trans. Commun. 2013, 61, 4754–4767. [Google Scholar] [CrossRef]
  5. Lee, S.; Zhang, R.; Huang, K. Opportunistic Wireless Energy Harvesting in Cognitive Radio Networks. IEEE Trans. Commun. 2013, 12, 4788–4799. [Google Scholar] [CrossRef] [Green Version]
  6. Chen, H.; Zhai, C.; Li, Y.; Vucetic, B. Cooperative Strategies for Wireless-Powered Communications: An Overview. IEEE Wirel. Commun. 2018, 25, 112–119. [Google Scholar] [CrossRef] [Green Version]
  7. Duy, T.T.; Son, V.N.; Tung, V.T.; Alexandropoulos, G.C.; Duong, T.Q. Outage performance of cognitive cooperative networks with relay selection over double-Rayleigh fading channels. IET Commun. 2016, 10, 57–64. [Google Scholar] [CrossRef]
  8. Ho-Van, K.; Sofotasios, P.C.; Alexandropoulos, G.C.; Freear, S. Bit error rate of underlay decode-and-forward cognitive networks with best relay selection. J. Commun. Netw. 2015, 17, 162–171. [Google Scholar] [CrossRef] [Green Version]
  9. Alexandropoulos, G.C.; Papadogiannis, A.; Sofotasios, P.C. A Comparative Study of Relaying Schemes with Decode and Forward over Nakagami-Fading Channels. J. Comput. Netw. Commun. 2011, 2011, 1–14. [Google Scholar] [CrossRef]
  10. Nguyen, T.; Tran, M.; Nguyen, T.-L.; Ha, D.-H.; Voznak, M. Multisource Power Splitting Energy Harvesting Relaying Network in Half-Duplex System over Block Rayleigh Fading Channel: System Performance Analysis. Electronics 2019, 8, 67. [Google Scholar] [CrossRef]
  11. Yao, C.-H.; Pei, C.-X.; Guo, J. Performance analysis of two-way AF cooperative networks with relay selection over Nakagami-m fading channels. Inf. Technol. Appl. 2015, 9, 307–311. [Google Scholar]
  12. Nguyen, T.; Quang Minh, T.; Tran, P.; Vozňák, M. Energy Harvesting over Rician Fading Channel: A Performance Analysis for Half-Duplex Bidirectional Sensor Networks under Hardware Impairments. Sensors 2018, 18, 1781. [Google Scholar] [CrossRef]
  13. Nguyen, T.N.; Minh, T.H.; Tran, P.T.; Voznak, M. Adaptive Energy Harvesting Relaying Protocol for Two-Way Half Duplex System Network over Rician Fading Channel. Wirel. Commun. Mob. Comput. 2018, 2018, 1–10. [Google Scholar] [CrossRef]
  14. Nguyen, T.N.; Duy, T.T.; Luu, G.T.; Tran, P.T.; Vozňák, M. Energy Harvesting-based Spectrum Access with Incremental Cooperation, Relay Selection and Hardware Noises. Radioengineering 2017, 26, 240–250. [Google Scholar] [CrossRef]
  15. Nguyen, T.N.; Minh, T.H.; Tran, P.T.; Voznak, M.; Duy, T.T.; Nguyen, T.L.; Tin, P.T. Performance Enhancement for Energy Harvesting Based Two-way Relay Protocols in Wireless Ad-hoc Networks with Partial and Full Relay Selection Methods. Ad Hoc Netw. 2019, 84, 178–187. [Google Scholar] [CrossRef]
  16. Gündüz, D.; Devillers, B. Two-hop Communication with Energy Harvesting. In Proceedings of the 2011 4th IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP), San Juan, Puerto Rico, 13–16 December 2011. [Google Scholar] [CrossRef]
  17. Chen, H.; Li, Y.; Rebelatto, J.L.; Uchoa-Filho, B.F.; Vucetic, B. Harvest-Then-Cooperate: Wireless-Powered Cooperative Communications. IEEE Trans. Signal Process. 2015, 63, 1700–1711. [Google Scholar] [CrossRef] [Green Version]
  18. Xiong, K.; Fan, P.; Zhang, C.; Letaief, K.B. Wireless Information and Energy Transfer for Two-Hop Non-Regenerative MIMO-OFDM Relay Networks. IEEE J. Sel. Areas Commun. 2015, 33, 1595–1611. [Google Scholar] [CrossRef]
  19. Okandeji, A.A.; Khandaker, M.R.; Wong, K.K. Two-way Beamforming Optimization for Full-duplex SWIPT Systems. In Proceedings of the 2016 24th European Signal Processing Conference (EUSIPCO), Budapest, Hungary, 29 August–2 September 2016. [Google Scholar] [CrossRef]
  20. Okandeji, A.A.; Khandaker, M.R.; Wong, K.K.; Zheng, Z. Joint Transmit Power and Relay Two-Way Beamforming Optimization for Energy-Harvesting Full-Duplex Communications. In Proceedings of the 2016 IEEE Globecom Workshops (GC Wkshps), Washington, DC, USA, 4–8 December 2016. [Google Scholar] [CrossRef]
  21. Hu, Y.; Zhu, Y.; Schmeink, A. Simultaneous Wireless Information and Power Transfer in Relay Networks with Finite Blocklength Codes. In Proceedings of the 2017 23rd Asia-Pacific Conference on Communications (APCC), Perth, Australia, 11–13 December 2017. [Google Scholar] [CrossRef]
  22. Halima, N.B.; Boujemaa, H. Exact and Approximate Symbol Error Probability of cooperative systems with best relay selection and all participating relaying using Amplify and Forward or Decode and Forward Relaying over Nakagami-m fading channels. KSII Trans. Internet Inf. Syst. 2018, 12, 81–108. [Google Scholar]
  23. Tseng, S.-M.; Lee, T.-L.; Ho, Y.-C.; Tseng, D.-F. Distributed space-time block codes with embedded adaptive AAF/DAF elements and opportunistic listening for multihop power line communication networks. Int. J. Commun. Syst. 2015, 30, e2950. [Google Scholar] [CrossRef] [Green Version]
  24. Li, Y.; Vucetic, B. On the Performance of a Simple Adaptive Relaying Protocol for Wireless Relay Networks. In Proceedings of the VTC Spring 2008—IEEE Vehicular Technology Conference, Singapore, 11–14 May 2008. [Google Scholar]
  25. Tseng, S.-M.; Liao, C.-Y. Distributed Orthogonal and Quasi-Orthogonal Space-Time Block Code with Embedded AAF/DAF Matrix Elements in Wireless Relay Networks with Four Relays. Wirel. Pers. Commun. 2013, 75, 1187–1198. [Google Scholar] [CrossRef]
  26. Katti, S.; Gollakota, S.; Katabi, D. Embracing Wireless Interference. ACM SIGCOMM Comput. Commun. Rev. 2007, 37, 397–408. [Google Scholar] [CrossRef]
  27. Qin, J.; Zhu, Y.; Zhe, P. Broadband Analog Network Coding With Robust Processing for Two-Way Relay Networks. IEEE Commun. Lett. 2017, 21, 1115–1118. [Google Scholar] [CrossRef]
  28. Nasir, A.A.; Zhou, X.; Durrani, S.; Kennedy, R.A. Relaying Protocols for Wireless Energy Harvesting and Information Processing. IEEE Trans. Commun. 2013, 12, 3622–3636. [Google Scholar] [CrossRef] [Green Version]
  29. Deng, Y.; Wang, L.; Elkashlan, M.; Kim, K.J.; Duong, T.Q. Generalized Selection Combining for Cognitive Relay Networks Over Nakagami-m Fading. IEEE Trans. Signal Process. 2015, 63, 1993–2006. [Google Scholar] [CrossRef]
  30. Nguyen, T.N.; Tran, P.T.; Minh, T.H.Q.; Voznak, M.; Sevcik, L. Two-Way Half Duplex Decode and Forward Relaying Network with Hardware Impairment over Rician Fading Channel: System Performance Analysis. Electron. Electr. Eng. 2018, 24, 74–78. [Google Scholar] [CrossRef]
  31. Owen, D.B.; Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Technometrics 1965, 7, 78. [Google Scholar] [CrossRef] [Green Version]
  32. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 8th ed.; Daniel, Z., Victor, M., Eds.; Academic Press: Amsterdam, The Netherlands, 2015. [Google Scholar]
  33. Bhatnagar, M.R. On the Capacity of Decode-and-Forward Relaying over Rician Fading Channels. IEEE Commun. Lett. 2013, 17, 1100–1103. [Google Scholar] [CrossRef] [Green Version]
  34. Mouapi, A.; Hakem, N. A New Approach to Design Autonomous Wireless Sensor Node Based on RF Energy Harvesting System. Sensors 2018, 18, 133. [Google Scholar] [CrossRef]
  35. Wang, C.; Li, J.; Yang, Y.; Ye, F. Combining Solar Energy Harvesting with Wireless Charging for Hybrid Wireless Sensor Networks. IEEE Trans. Mob. Comput. 2018, 17, 560–576. [Google Scholar] [CrossRef]
  36. Chu, Z.; Zhou, F.; Zhu, Z.; Hu, R.Q.; Xiao, P. Wireless Powered Sensor Networks for Internet of Things: Maximum Throughput and Optimal Power Allocation. IEEE Internet Things J. 2018, 5, 310–321. [Google Scholar] [CrossRef] [Green Version]
Figure 1. System model.
Figure 1. System model.
Electronics 08 00448 g001
Figure 2. EH and IT of the proposed system.
Figure 2. EH and IT of the proposed system.
Electronics 08 00448 g002
Figure 3. The outage probability (a) and ergodic capacity (b) versus Ps/N0 while K = 2.
Figure 3. The outage probability (a) and ergodic capacity (b) versus Ps/N0 while K = 2.
Electronics 08 00448 g003
Figure 4. The outage probability (a) and ergodic capacity (b) versus Ps/N0 while ρ = 0.5.
Figure 4. The outage probability (a) and ergodic capacity (b) versus Ps/N0 while ρ = 0.5.
Electronics 08 00448 g004
Figure 5. Comparison of the outage probability (a) and ergodic capacity (b) ersus Ps/N0.
Figure 5. Comparison of the outage probability (a) and ergodic capacity (b) ersus Ps/N0.
Electronics 08 00448 g005
Figure 6. The outage probability versus η.
Figure 6. The outage probability versus η.
Electronics 08 00448 g006
Figure 7. The outage probability (a) and ergodic capacity (b) versus the numbers of users.
Figure 7. The outage probability (a) and ergodic capacity (b) versus the numbers of users.
Electronics 08 00448 g007
Table 1. Simulation parameters.
Table 1. Simulation parameters.
SymbolNameValues
η Energy harvesting efficiency0.8
λ s r Mean of | h s r | 2 0.5
λ r d Mean of | h r d | 2 0.5
m γ 1 Nakagami m-factor3
zSNR threshold1
Ps/N0Source power to noise ratio0–20 dB
RSource rate0.5 bit/s/Hz
KNumber of users1–6
mPathloss exponent3
dsr = dithe distance of S-R link and R-D link, respectively 0.85

Share and Cite

MDPI and ACS Style

Nguyen, T.N.; Tran, M.; Nguyen, T.-L.; Ha, D.-H.; Voznak, M. Performance Analysis of a User Selection Protocol in Cooperative Networks with Power Splitting Protocol-Based Energy Harvesting Over Nakagami-m/Rayleigh Channels. Electronics 2019, 8, 448. https://doi.org/10.3390/electronics8040448

AMA Style

Nguyen TN, Tran M, Nguyen T-L, Ha D-H, Voznak M. Performance Analysis of a User Selection Protocol in Cooperative Networks with Power Splitting Protocol-Based Energy Harvesting Over Nakagami-m/Rayleigh Channels. Electronics. 2019; 8(4):448. https://doi.org/10.3390/electronics8040448

Chicago/Turabian Style

Nguyen, Tan N., Minh Tran, Thanh-Long Nguyen, Duy-Hung Ha, and Miroslav Voznak. 2019. "Performance Analysis of a User Selection Protocol in Cooperative Networks with Power Splitting Protocol-Based Energy Harvesting Over Nakagami-m/Rayleigh Channels" Electronics 8, no. 4: 448. https://doi.org/10.3390/electronics8040448

APA Style

Nguyen, T. N., Tran, M., Nguyen, T. -L., Ha, D. -H., & Voznak, M. (2019). Performance Analysis of a User Selection Protocol in Cooperative Networks with Power Splitting Protocol-Based Energy Harvesting Over Nakagami-m/Rayleigh Channels. Electronics, 8(4), 448. https://doi.org/10.3390/electronics8040448

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