Next Article in Journal
Impact on Reduction of Pollutant Emissions from Passenger Cars when Replacing Euro 4 with Euro 6d Diesel Engines Considering the Altitude Influence
Next Article in Special Issue
Realization of Licensed/Unlicensed Spectrum Sharing Using eICIC in Indoor Small Cells for High Spectral and Energy Efficiencies of 5G Networks
Previous Article in Journal
Pipe Sizing for Novel Heat Distribution Technology
Previous Article in Special Issue
Effect of Prediction Error of Machine Learning Schemes on Photovoltaic Power Trading Based on Energy Storage Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Performance Analysis of Amplify-and-Forward Bidirectional Relaying over Nakagami-m Fading Channels with Arbitrary Parameters

1
School of Information Engineering, Nanchang University, Nanchang 330031, China
2
School of Information Engineering, East China Jiaotong University, Nanchang 330013, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(7), 1277; https://doi.org/10.3390/en12071277
Submission received: 22 February 2019 / Revised: 25 March 2019 / Accepted: 28 March 2019 / Published: 3 April 2019

Abstract

:
The exact performance of amplify-and-forward (AF) bidirectional relay systems is studied in generalized and versatile Nakagami-m fading channels, where the parameter m is an arbitrary positive number. We consider three relaying modes: two, three, and four time slot bidirectional relaying. Closed form expressions of the moment generating function (MGF), higher order moments of signal-to-noise ratio (SNR), ergodic capacity, and average signal error probability (SEP) are derived, which are different from previous works. The obtained expressions are very concise, easy to calculate, and evaluated instantaneously without a complex summation operation, in contrast to the nested multifold numerical integrals and truncated infinite series expansions used in previous work, which lead to computational inefficiency, especially when the fading parameter m increases. Simulation results corroborate the correctness and tightness of the theoretical analysis.

1. Introduction

Relay technology is a promising solution for 5G communication in the future. The Nakagami-m environment is widely studied in relay networks. When it comes to the Nakagami-m distribution, most research institutions prefer an integer parameter m. For example, an approximate formula of the outage probability was studied in Reference [1] for bidirectional relaying systems with co-channel interferes. An exact outage expression was derived in Reference [2], when m is an integer. While multi-antenna technology was introduced in Reference [3]. The exact and asymptotic outage probability formulas were derived in Reference [4] when a bidirectional relaying system suffered from asymmetric traffic, which was quantified by a traffic pattern indicator and meant that the traffic transmitted and received by the two transceivers was different. Then, the author in Reference [4] took the asymptotic outage probability as an optimization objective, in order to find the optimal power allocation and relay location. Using the expectation property of the logarithmic function, the lower bound of the average rate was obtained in Reference [5] for bidirectional amplify-and-forward (AF) relaying networks that suffer from Nakagami-m fading channels. Resorting to polynomial theory, the lower bound involved many special functions, including the incomplete gamma function, factorial, Euler psi function, and hypergeometric function. By enlarging the original integral region, an approximation of the outage probability was provided in Reference [6]. Obviously, the authors in Reference [6] failed to get an exact formula in the original integral region, so they had to change the integral region. Some common average signal error probability (SEP) expressions were determined in Reference [7] for integer m, such as the 32-cross quadrature amplitude modulation, differentially encoded quadrature phase shift keying, and π /4 quadrature phase shift keying modulation schemes. The approximate probability density function and asymptotic outage probability in high signal-to-noise ratio (SNR) region were presented in Reference [8] for AF bidirectional relaying.
The papers published about non-integer m are obviously not as many as those about integer m. The outage probability and SEP were analyzed in Reference [9] for cascading Nakagami-m fading channels in AF bidirectional relaying systems. The diversity order of the outage probability was derived for an AF bidirectional relay network with beam-forming [10]. The authors in Reference [11] claimed the asymptotic formula was a gap in bidirectional relay networks. This claim indicated the difficulty for arbitrary m values, especially for non-integer m. When the SNR tended to be infinite, the asymptotic outage probability and average SEP were derived with the aid and knowledge of limit operations [11]. An exception appeared in Reference [12], where the exact outage probability, average SEP, and average sum rate were presented for non-integer m. However, this non-integer value only defined m as an integer plus one half. Clearly, this rule was too narrow for an arbitrary parameter and many other non-integers were not included in Reference [12].
Although the Nakagami-m fading channel has been extensively studied in bidirectional AF relaying networks, there are still some disadvantages. First, most previous works prefer an integer value for parameter m in the Nakagami-m environment for bidirectional relaying systems [1,2,3,4,5,6,7,8]. Obviously, it is easier to study an integer than any real number. When m is an integer, the complex gamma channel distribution is easily expanded into a finite summation of the exponential function and the power function, and so a closed form expression is easily yielded with integer m. But integer m often limits the application of the Nakagami-m fading channel, because the actual complex environment may not be well modeled by an integer parameter m. Also, the Nakagami-m fading channel does not specify that m is an integer. This brings the second shortcoming. When parameter m takes a non-integer value, a lower or upper bound, truncated infinite series expansion, and nested multifold numerical integrals have to be utilized to evaluate approximate performance [9,10,11,12]. This is because the gamma function is only expanded into an infinite series of exponential and power functions when m takes a non-integer value. Additionally, infinite series expansions appear in References [2,3] even for integer m, let alone non-integer m. Finally, most of the previous work links performance with the probability density function (PDF) or cumulative distribution function (CDF). This method is effective for simple PDF or CDF. However, for the complex gamma channel distribution, this method often fails in finding closed-form solutions when m is an arbitrary number. To our best knowledge, no investigation on the simple exact closed form formula has been done for arbitrary parameter m, due to the inherent complexity of this problem.
Stimulated by the above observations, this paper studies the accurate performance of AF bidirectional relay systems in the Nakagami-m fading channel with arbitrary parameter m. In our analysis, closed form expressions for the moment generating function (MGF), higher order moments of the SNR, ergodic capacitym and average SEP are derived, which differ from previous work. In addition, our formulas without summation operations are very easy and instantaneous to calculate, in sharp contrast to previous work, which had to resort to the evaluation of truncated infinite series or multifold numerical integrals. Finally, some mathematical notations are summarized in Table 1.

2. System Model

Consider an AF relaying system, where two transceivers S 1 and S 2 exchange messages via one relay R. We consider three bidirectional modes: two, three, and four time slot bidirectional relaying [13]. The transmit powers of the two transceivers and the relay are denoted as p 1 , p 2 , and p r , respectively. In three time slot bidirectional relaying, α 1 and α 2 denote the scaling factors of the total relay power p r assigned to S 1 and S 2 , respectively. Let h 1 and h 2 be the channel coefficients pertaining to the links S 1 R and S 2 R , respectively. N 0 represents the noise variance. All channel links suffer from independent Nakagami-m fading, where m is an arbitrary positive number. The equivalent SNR between two transceivers can be expressed by [13]:
γ = a γ 1 γ 2 γ 1 + γ 2 ,
where a, γ 1 , and γ 2 depend on the different relaying modes and are given in Table 2.

3. Performance Analysis

The PDF and CDF are often used to obtain the performance expression; most researchers write the performance expression as an integrand of the PDF or CDF. However, for an arbitrary parameter m, this method can only obtain approximate solutions by means of a lower or upper bound, nested multifold numerical integrals, or truncated infinite series. Even when the parameter m is an integer value, the performance expression is often the summation of multiple special functions. Obviously, the previous formulas are computationally inefficient and time consuming, especially when the parameter m increases [1,2,5,7,11]. To avoid this dilemma, we make use of the MGF for the inverse SNR. As will be seen later, our derived expressions allow us to evaluate and quantify performance at a low complexity.

3.1. Moment Generating Function

The MGFs of inverse γ 1 and γ 2 are found in (Equations (32) and (33), Reference [12]). Due to the statistical independence of γ 1 and γ 2 , the MGF of inverse γ is achieved via:
M 1 γ s = M 1 γ 1 s a M 1 γ 2 s a = 4 β s m s 2 β r m r 2 Γ m s Γ m r s a m s + m r 2 K m s 2 β s s a K m r 2 β r s a ,
where m s m r is the fading parameter between S 1 S 2 and R. β s and β r are scale parameters on the corresponding links. Γ · is the gamma function (Equation (8.311), [14]) and K ν · is the ν th order modified Bessel function of the second kind (Equation (8.432), [14]). According to the relationship between the MGF of a random variable and its inverse (Equation (42), [15]), the MGF of the SNR γ is given by:
M γ s = 1 - 2 s 0 J 1 2 s ξ M 1 γ ξ 2 d ξ = 1 - 2 s 0 J 1 2 s ξ H 0 , 2 2 , 0 β s ξ 2 a | 0 , 1 , m s , 1 H 0 , 2 2 , 0 β r ξ 2 a | 0 , 1 , m r , 1 d ξ = 1 - 1 Γ m s Γ m r H 2 , 0 : 0 , 0 , 2 : 2 1 , 0 , 0 , 2 , 2 β s a s β r a s | 1 , 1 , 0 , 1 - ; - - 0 , 1 , m s , 1 ; 0 , 1 , m r , 1 ,
where J ν · is the ν th order Bessel function of the first kind (Equation (8.402), [14]), H x and H x , y are the univariable and bivariable Fox’s functions (Equations (1.1.1) and (2.2.1), [16]), respectively. The second equality in Equation (3) is obtained by rewriting the Bessel function in terms of Fox’s function (Equation (2.9.19), [17]), i.e.,:
2 β s a m 2 K m 2 β s a = H 0 , 2 2 , 0 β s a | 0 , 1 , m , 1 .
From M γ s , the average SEP of some simple modulations can be immediately obtained. For example, the SEP of the differential phase shift keying (DPSK) modulation is given by P e = 1 2 M γ 1 .

3.2. Higher Order Moments

Dependent on (Equation (34), Reference [18]) and (Equation (6.576.4), [14]), the nth order moment of the SNR γ is given by:
E γ n = - 1 n ! 0 s n M 1 γ s s d s = 1 n - 1 ! 0 s n - 1 M 1 γ s d s = a n Γ m s + n Γ m r + n Γ m s + m r + n β r n Γ m s Γ m r Γ m s + m r + 2 n 2 F 1 n , m r + n ; m s + m r + 2 n ; 1 - β s β r = n - 1 ! a β r n B m s + m r + n , n B m s , n B m r , n 2 F 1 n , m r + n ; m s + m r + 2 n ; 1 - β s β r ,
where E · is the expectation operation, B · , · is the beta function (Equation (8.380.1), [14]), and 2 F 1 · is the Gaussian hypergeometric function (Equation (9.111), [14]).

3.3. Ergodic Capacity

Following a similar procedure as for higher order moments, according to the relationship between the ergodic capacity and the MGF of the inverse SNR (Equation (58), [18]), the ergodic capacity is given by:
E η log 2 1 + γ = η ln 2 0 1 - e - s s M 1 γ s d s = η ln 2 0 H 1 , 2 1 , 1 s | 0 , 1 0 , 1 , - 1 , 1 M 1 γ s d s = η a ln 2 β s Γ m s Γ m r H 2 , 1 : 0 , 0 , 2 : 2 2 , 1 , 0 , 1 , 2 a β s β r β s | 1 , 1 , m s + 1 , 1 0 , 1 ; - - 0 , 1 , - 1 , 1 ; 0 , 1 , m r , 1 ,
where η = 1 2 , 1 3 , and 1 4 correspond to two, three, and four time slot bidirectional relaying, respectively. The second equality in Equation (6) follows from (Equation (37), [19]):
1 - e - s s = H 1 , 2 1 , 1 s | 1 , 1 1 , 1 , 0 , 1 s = H 1 , 2 1 , 1 s | 0 , 1 0 , 1 , - 1 , 1 .

3.4. Average SEP

Similarly, using (References [20], Equation (21)) and ([16], Equation (2.5.1)), the average SEP is given by:
E b e r f c c γ = - b π 0 G 1 , 3 2 , 0 c s | 1 1 2 , 0 , 0 M 1 γ s s d s = b - 0 b sin 2 c s π s M 1 γ s d s = b - b π Γ m s Γ m r H 2 , 0 : 0 , 0 , 2 : 2 1 , 0 , 0 , 2 , 2 β s a c β r a c | 1 2 , 1 , 0 , 1 - ; - - 0 , 1 , m s , 1 ; 0 , 1 , m r , 1 ,
where b and c are modulation specific constants, e r f c · is the complementary error function, and G · is Meijer’s G function (Equation (9.301), [14]). The third equality in Equation (8) follows from (Equation (2.9.7), [17]):
sin 2 c s π s = H 0 , 2 1 , 0 c s | - 1 2 , 1 , 0 , 1 s = c H 0 , 2 1 , 0 c s | - - 1 2 , 1 , - 1 , 1
Now, closed form expressions are found for the MGF, higher order moments, ergodic capacity, and average SEP in the AF bidirectional relaying system over a Nakagami-m fading environment, when m is an arbitrary positive number. These formulas are very compact and easy to calculate without complex numerical integration and summation operations. As far as we know, the formulas we provide are brand new and have not been publicly reported. However, in previous work, finding the closed form solution is a task far from non-intricate, and involved truncated infinite series expansions and nested multifold numerical integrals, which become verbose and slow even for small m values. Obviously, our derived formulas lower computational complexity.

4. Simulation Results

Simulation results are shown to examine the impacts of our compact formulas. The channel coefficients of all links are normalized to unity. For simplicity, we assume equal transmit power at the two transceivers and the relay node. The average SNR per hop is defined as p 1 / N 0 . As an example, we show the performance of the two time slot bidirectional relay system, while the three and four time slot bidirectional patterns can be presented similarly.
Figure 1 and Figure 2 show the average SEPs of DPSK and 4-ary pulse amplitude modulation (4PAM), respectively. Six sets of random non-integer fading parameters are selected. Other arbitrary parameters m can be simulated in a similar way. The reason for two figures here is that the error probability expressions of different modulation methods are different. The error probabilities of DPSK and 4PAM are derived from the MGF in Equation (3) and the complementary error function in Equation (8), respectively. Figure 1 and Figure 2 demonstrate the universality and correctness of the formulas. It is further observed that the two figures highlight the same decline speed, which implies the same slope and diversity order.
Figure 3 plots the first order moment E γ while other higher order moments can be drawn similarly using Equation (5). The theoretical values accurately coincide with the simulation curves. Notice that because of the symmetry of the SNR γ , Equation (5) still holds when the positions of β s and β r are swapped.
Finally, the ergodic capacity is shown in Figure 4, where the theoretical and simulated curves are indistinguishable from each other. The ergodic capacity increases as the fading parameters increase due to the improvement of the channel quality. It is observed that the capacity curves nearly keep the same slope in the high SNR region but shift to the right as m increases. This implies that m has no impact on the multiplexing gain, which always stays constant.

5. Conclusions

In this paper, we consider three types of bidirectional relay modes: two, three, and four time slot bidirectional relaying. In the Nakagami-m fading channel with arbitrary m value, we derive new, highly exact expressions for the MGF, higher order moments, ergodic capacity, and average SEP. All expressions are closed. In particular, the expressions are very concise and easy to calculate, as opposed to complex truncated infinite series expansions or nested multifold numerical integrals. Unlike most researchers who are used to starting a performance analysis with the PDF or CDF, we explore another effective tool, the MGF of the inverse SNR, which is often overlooked. Moreover, our analysis is based on an arbitrary value of m, in line with actual complex environments.

Author Contributions

Conceptualization, D.Q. and Y.W.; Methodology, D.Q. and T.Z.; Software, D.Q.; Writing—original draft preparation, D.Q.; Writing—review and editing, D.Q.

Funding

This work was supported by the National Natural Science Foundation of China under Grant 61761030 and by the China Postdoctoral Science Foundation under Grant 2017M622103.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. da Costa, D.B.; Ding, H.; Yacoub, M.D.; Ge, J. Two-way relaying in interference limited AF cooperative networks over Nakagami-m fading. IEEE Trans. Veh. Technol. 2012, 61, 3766–3771. [Google Scholar] [CrossRef]
  2. Upadhyay, P.K.; Prakriya, S. Performance of two-way opportunistic relaying with analog network coding over Nakagami-m fading. IEEE Trans. Veh. Technol. 2011, 60, 1965–1971. [Google Scholar] [CrossRef]
  3. Yang, K.; Yang, N.; Xing, C.; Wu, J. Relay antenna selection in MIMO two-way relay networks over Nakagami-m fading channels. IEEE Trans. Veh. Technol. 2014, 63, 2349–2362. [Google Scholar] [CrossRef]
  4. Shukla, M.K.; Yadav, S.; Purohit, N. Performance evaluation and optimization of traffic-aware cellular multiuser two-way relay networks over Nakagami-m fading. IEEE Syst. J. 2018, 12, 1933–1944. [Google Scholar] [CrossRef]
  5. Yang, J.; Fan, P.; Ding, Z. Capacity of AF two-way relaying with multiuser scheduling in Nakagami-m fading. Electron. Lett. 2012, 48, 1432–1434. [Google Scholar] [CrossRef]
  6. Guo, H.; Ge, J. Performance analysis of two-way opportunistic relaying over Nakagami-m fading channels. Electron. Lett. 2011, 47, 150–152. [Google Scholar] [CrossRef]
  7. Kumar, N.; Bhatia, V. Exact ASER analysis of rectangular QAM in two-way relaying networks over Nakagami-m fading channels. IEEE Wirel. Lett. 2016, 5, 548–551. [Google Scholar] [CrossRef]
  8. Lei, X.; Fan, L.; Michalopoulos, D.S.; Fan, P.; Hu, R.Q. Outage probability of TDBC protocol in multiuser two-way relay systems with Nakagami-m fading. IEEE Commun. Lett. 2013, 17, 487–490. [Google Scholar] [CrossRef]
  9. İlhan, H. Performance analysis of two-way AF relaying systems over cascaded Nakagami-m fading channels. IEEE Signal Process. Lett. 2012, 19, 332–335. [Google Scholar] [CrossRef]
  10. Guo, H.; Ge, J. Diversity order of multiuser two-way relay networks with beamforming over Nakagami-m fading channels. Sci. China Inf. Sci. 2011, 54, 1986–1990. [Google Scholar] [CrossRef]
  11. Yang, J.; Fan, P. Asymptotic analysis of multi-branch two-way amplify-and-forward relaying in Nakagami-m fading with arbitrary fading parameter. J. Mod. Transp. 2012, 20, 249–254. [Google Scholar] [CrossRef]
  12. Yang, J.; Fan, P.; Duong, T.Q.; Lei, X. Exact performance of two-way AF relaying in Nakagami-m fading environment. IEEE Trans. Wirel. Commun. 2011, 10, 980–987. [Google Scholar] [CrossRef]
  13. Louie, R.H.Y.; Li, Y.; Vucetic, B. Practical physical layer network coding for two-way relay channels: Performance analysis and comparison. IEEE Trans. Wirel. Commun. 2010, 9, 764–777. [Google Scholar] [CrossRef]
  14. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products, 7th ed.; Academic Press: Cambridge, MA, USA, 2007. [Google Scholar]
  15. Trigui, I.; Affes, S.; Stéphenne, A. On the ergodic capacity of amplify-and-forward relay channels with interference in Nakagami-m fading. IEEE Trans. Commun. 2013, 61, 3136–3145. [Google Scholar] [CrossRef]
  16. Mathai, A.M.; Saxena, R.K. The H-Function with Applications in Statistics and Other Disciplines; John Wiley & Sons: Hoboken, NJ, USA, 1978. [Google Scholar]
  17. Kilbas, A.A.; Saigo, M. H-Transforms: Theory and Applications; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  18. Aalo, V.A.; Peppas, K.P.; Efthymoglou, G.P.; Alwakeel, M.M.; Alwakeel, S.S. Serial amplify-and-forward relay transmission systems in Nakagami-m fading channels with a poisson interference field. IEEE Trans. Veh. Technol. 2014, 63, 2183–2196. [Google Scholar] [CrossRef]
  19. Aalo, V.A.; Efthymoglou, G.P.; Soithong, T.; Alwakeel, M.; Alwakeel, S. Performance analysis of multi-hop amplify-and-forward relaying systems in Rayleigh fading channels with a poisson interference field. IEEE Trans. Wirel. Commun. 2014, 13, 24–35. [Google Scholar] [CrossRef]
  20. Fredj, K.B.; Aïssa, S. Performance of amplify-and-forward systems with partial relay selection under spectrum-sharing constraints. IEEE Trans. Wirel. Commun. 2012, 11, 500–504. [Google Scholar] [CrossRef]
Figure 1. Average signal error probability (SEP) of the differential phase shift keying (DPSK) modulation. SNR: signal-to-noise ratio.
Figure 1. Average signal error probability (SEP) of the differential phase shift keying (DPSK) modulation. SNR: signal-to-noise ratio.
Energies 12 01277 g001
Figure 2. Average SEP of 4-ary pulse amplitude modulation (4PAM).
Figure 2. Average SEP of 4-ary pulse amplitude modulation (4PAM).
Energies 12 01277 g002
Figure 3. Comparison of the average SNR.
Figure 3. Comparison of the average SNR.
Energies 12 01277 g003
Figure 4. Comparison of the ergodic capacity.
Figure 4. Comparison of the ergodic capacity.
Energies 12 01277 g004
Table 1. Mathematical notations.
Table 1. Mathematical notations.
NotationDescription
α 1 Scaling factor of the total relay power assigned to S 1
α 2 Scaling factor of the total relay power assigned to S 2
B · , · Beta function
e r f c · Complementary error function
2 F 1 · Gaussian hypergeometric function
Γ · Gamma function
G · Meijer’s G function
H · Fox’s function
K ν · The ν th order modified Bessel function of the second kind
Table 2. Parameter values for the three relay modes.
Table 2. Parameter values for the three relay modes.
Relay ModeTransmission Directiona γ 1 γ 2
Two time slot bidirectional relay S 2 R S 1 p r p 1 + p r p 2 h 2 2 N 0 p r + p 1 h 1 2 N 0
S 1 R S 2 p r p r + p 2 p 1 h 1 2 N 0 p r + p 2 h 2 2 N 0
Three time slot bidirectional relay S 2 R S 1 p r p r + α 1 p 1 α 2 p 2 h 2 2 N 0 p r + α 1 p 1 h 1 2 N 0
S 1 R S 2 p r p r + α 2 p 2 α 1 p 1 h 1 2 N 0 p r + α 2 p 2 h 2 2 N 0
Four time slot bidirectional relay S 1 R S 2 1 p 1 h 1 2 N 0 p r h 2 2 N 0
S 2 R S 1 1 p 2 h 2 2 N 0 p r h 1 2 N 0

Share and Cite

MDPI and ACS Style

Qin, D.; Wang, Y.; Zhou, T. Exact Performance Analysis of Amplify-and-Forward Bidirectional Relaying over Nakagami-m Fading Channels with Arbitrary Parameters. Energies 2019, 12, 1277. https://doi.org/10.3390/en12071277

AMA Style

Qin D, Wang Y, Zhou T. Exact Performance Analysis of Amplify-and-Forward Bidirectional Relaying over Nakagami-m Fading Channels with Arbitrary Parameters. Energies. 2019; 12(7):1277. https://doi.org/10.3390/en12071277

Chicago/Turabian Style

Qin, Dong, Yuhao Wang, and Tianqing Zhou. 2019. "Exact Performance Analysis of Amplify-and-Forward Bidirectional Relaying over Nakagami-m Fading Channels with Arbitrary Parameters" Energies 12, no. 7: 1277. https://doi.org/10.3390/en12071277

APA Style

Qin, D., Wang, Y., & Zhou, T. (2019). Exact Performance Analysis of Amplify-and-Forward Bidirectional Relaying over Nakagami-m Fading Channels with Arbitrary Parameters. Energies, 12(7), 1277. https://doi.org/10.3390/en12071277

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