Next Article in Journal
Theoretical Study on Design of Composite Lining Structure with Rigid–Flexible Combination for Trenchless Pipeline Rehabilitation
Previous Article in Journal
Chinese Named Entity Recognition Based on Knowledge Based Question Answering System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Sum of α-μ/Inverse Gamma Variates with Applications to Diversity Receivers

1
School of Information Engineering, Minzu University of China, Beijing 100081, China
2
Key Laboratory of Mining Disaster Prevention and Control, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(11), 5375; https://doi.org/10.3390/app12115375
Submission received: 20 April 2022 / Revised: 20 May 2022 / Accepted: 22 May 2022 / Published: 26 May 2022
(This article belongs to the Section Electrical, Electronics and Communications Engineering)

Abstract

:
How to model shadow fading by applying the inverse gamma (IGA) distribution has recently gained widespread attention in wireless transmissions. However, the sum of α-μ/IGA variates, and its applications under independent and/or correlated scenarios, have yet to be addressed in open research works. Hence, this paper provides a systematic investigation of the α-μ/IGA model. First, we derive the expressions of the fundamental statistics of the univariate and bivariate α-μ/IGA models including the probability density function, cumulative distribution function, and moment generating function, and propose a mixture α-μ model to approximate the α-μ/IGA model. Then, according to the above statistical expressions, the statistical properties of the sum of α-μ/IGA variates are obtained and employed in the maximal ratio combining receivers. Third, the novel exact and approximated expressions of some performance metrics of interest are also derived, for instance, the average bit/symbol error probability, the outage probability, the average channel capacity, and the effective rate. Moreover, to prove the asymptotic properties of the performance metrics at the high signal-to-noise regions, some examples are performed. Finally, we explore numerical analysis and simulations to demonstrate the accuracy of the theoretical expressions under the different channel and system parameters. These results provide some significant insights into the reliability design and deployment of some conventional and emerging wireless communication applications.

1. Introduction

In the past decades, how to model wireless fading channels accurately has been an important and inevitable topic in the research of conventional and emerging wireless communication systems, such as cellular networks, wearable communications, and vehicular wireless networks. Generally, the received signals of wireless transmission undergo either small-scale (multipath) fading or large-scale (shadowing) fading, or even the composite effects of both. For small-scale fading, some classical fading models have been reported in [1], namely the Rayleigh, the Rice, the Nakagami-m, the Weibull, and so on. To accurately characterize the envelope fluctuation of the received signals and better fit empirical data encountered in some emerging wireless applications, some more general and sophisticated versions of the above multipath fading models have also been proposed, such as η-μ, κ-μ, α-μ, α-η-μ, α-κ-μ and α-η-κ-μ in [2] and references therein. For large-scale fading, the lognormal (LN) model is a typical shadowing fading example [1]. Unfortunately, the main disadvantage of this model is that its probability density function (PDF) contains an integral form which limits its applications for performance analysis purposes. For this reason, two substituted versions of the LN model have been proposed to approximate it, namely the Gamma (GA) model [3] and the inverse Gaussian (IG) model [4]. However, in some practical wireless communication scenarios, the composite effects of both multipath and shadowing fading become a dominating fading phenomenon. To this end, the composite fading models have been extensively explored by combining multipath fading with shadowing, which comprise the LN-based composite model (such as Rayleigh/LN, Rice/LN, Nakagami/LN, Weibull/LN, α-η-μ/LN, and α-κ-μ/LN [5]), the GA-based composite model (such as Rayleigh/GA, Nakagami/GA, Weibull/GA, η-μ/GA, κ-μ/GA, and α-μ/GA [6]), the IG-based model (such as Rayleigh/IG, Nakagami/IG, η-μ/IG, κ-μ/IG, and α-μ/IG [7]), the line-of-sight (LOS) shadowing composite model (such as Rician shadowed [1] and κ-μ shadowed [8]), and the double shadowing composite model [9,10]. Furthermore, a double α-μ model in [11,12] has also been proposed to describe the composite multipath-shadowing models, where the first α-μ model denotes the multipath fading and the second one is the shadowing fading.
More recently, an inverse Gamma (IGA) model to characterize shadowing has been widely studied to achieve an appropriate balance between modeling accuracy and calculation complexity. For the empirical field measurements of the IGA model, the authors in [13] provided the experimental studies through wearable communication channels tests to prove that the underlying IGA model shadowing is very available compared with the LN and the GA models. At the same time, the authors in [14,15] also confirmed that describing light and moderate shadowing with the IGA model is a reasonable alternative for the LN and the GA models by considering more testing scenarios, including both indoor and outdoor with many different frequency bands. More importantly, [14,15] showed some mathematical expressions based on the IGA model are more tractable and simpler than those by using the LN or the IG model. By employing the IGA shadowing model, an unmanned aerial vehicle-to-ground communications channel model was presented in [16] which can provide a notable fit for empirical data.
For the application of the underlying IGA model, [17,18] first obtained the average symbol error rate of multiple amplitude modulation and the ergodic capacity over κ-μ/IGA and η-μ/IGA composite fading channels, respectively. After that, the same authors pointed out in [19] that the κ-μ/IGA and η-μ/IGA models not only provide a good fit for modeling the composite fading channels but also can be extended to the conventional and emerging wireless communication systems, such as cellular, wearable, and vehicular wireless applications. However, they only derived the statistical characteristics of the η-μ/IGA and κ-μ/IGA models, which include the PDF, the cumulative distribution function (CDF), the moment generating function (MGF), the higher-order moments, and the amount of fading. In [20], the effective capacity was analyzed over η-μ/IGA and κ-μ/IGA composite fading channels. In [21], the average symbol error probability (ASEP) and the channel capacity with adaptive schemes were obtained over η-μ/IGA and κ-μ/IGA composite fading channels. Moreover, some performance analysis was investigated over the inverse gamma pure shadowed fading channels in [22,23]. As a counterpart of η-μ/IGA and κ-μ/IGA, the α-μ/IGA composite model has also been involved in [24,25,26], where [24] obtained the performance analysis of the single-branch system with the aid of Fox’H functions, [25] developed the energy detection of single-branch and maximal-ratio combining (MRC) systems over independent and identical distributed (i.i.d.) conditions, and [26] provided an approximated performance analysis of MRC and equal gain combining (EGC) by using a mixed IGA distribution under i.i.d. conditions. In recent work, the authors in [27] presented the α-η-μ/IGA and α-κ-μ/IGA channel models and evaluated the performance of multi-hop links.
On the other hand, the Fisher-Snedecor F fading channel model in [28] can be considered as an alternative version of the IGA-based composite model based upon the fact that the square root of an IGA random variable (RV) follows an inverse Nakagami-m distribution. So far, this model has obtained extensive interest and a great number of works have been done in [29,30,31] and references therein. By considering the nonlinearity of the composite envelope of fading signals, the α- F , the α-η- F , the α-κ- F , and α-η-κ- F composite fading models were presented in [32,33,34], respectively. When α = 2, these three models can reduce to the Fisher-Snedecor F , the η-μ/IGA, and κ-μ/IGA, respectively. Although the α-η-μ/IGA, α-κ-μ/IGA, the α-η- F , the α-κ- F , and α-η-κ- F are versatile, the mathematical forms of their PDFs are still intractable, which renders them rather difficult to analyze the performance of wireless communication systems exactly.
Based on the above observation, most of the works involving the IGA shadowing model focused on the Fisher-Snedecor F , the η-μ/IGA, and κ-μ/IGA fading channels. Although the α-μ/IGA model has been considered in [24,25,26], to the best of the authors’ knowledge, a comprehensive investigation of wireless communication systems over α-μ/IGA composite fading channels has not yet been involved in the independent and non-identical distributed (i.n.i.d.) scenarios and/or the correlated scenarios. To this effect, we will further extend our work in [24] to provide a thorough investigation and some significant insights into the α-μ/IGA and its application in wireless communication systems. In addition, this model can also be potentially applied to describe the distribution of population densities and urban activities investigated in [35]. In this paper, we first drive the exact analytical formulations of the statistical characteristics of univariate and bivariate α-μ/IGA distributions and compare the differences between the α-μ/IGA and the α- F models, as well as a new mixture α-μ model and the mixture IGA model. Then, capitalizing on the obtained statistical expressions, the PDF, the CDF, and the MGF of the sum of α-μ/IGA variates are derived in the i.n.i.d. and the correlated scenarios in the light of the multivariate Fox’s H-function by utilizing the Mellin-Barnes type contour integral, respectively. At the same time, some exact approximated expressions are also deduced by using the mixture α-μ model. Thirdly, some performance metrics of interest for the single-link and MRC systems are investigated, respectively, including the outage probability (OP), the average bit/symbol error probability (ABEP/ASEP), and the average channel capacity and the effective rate. In addition, the asymptotic analysis of the ABEP of DPSK (Differential Phase Shift Keying) and the OP are considered at high signal-to-noise (SNR) regions. Finally, the numerical analysis and simulations are carried out to confirm the accuracy of our theoretical derivations in various scenarios.
The remainder of this paper is organized as follows: In Section 2, the statistical features of the univariate and bivariate α-μ/IGA distributions are investigated, and some comparisons are discussed. Some exact and approximated analytical expressions of the statistical properties of the sum of α-μ/IGA variates are obtained in Section 3. Section 4 derives the analytical expressions of performance metrics including the OP, the ABEP/ASEP, the channel capacity, and the effective capacity. In Section 5, we show some asymptotic analysis by using several examples in high SNR regions. Then, the numerical evaluations and simulations are performed and discussed in Section 6. Finally, we conclude the whole paper in Section 7.

2. Statistical Characteristics of α-μ/IGA Distribution

2.1. Univariate Case

In the α-μ/IGA composite fading environments, the amplitude of the multipath envelopes of the α-μ distribution suffers from random fluctuations induced by an IGA shadowing RV. Thus, the received signal envelope R can be expressed as
R 2 = Z R α u α = i = 1 μ Z ( X i 2 + Y i 2 ) ,
where R α μ denotes the envelope of α-μ distribution, Z follows the IGA distribution, α is the non-linear parameter and α > 0, μ represents the number of multipath clusters, and μ   1 / 2 , X i , and Y i are mutually independent Gaussian RVs with zero-mean and average power σ 2 , which represent the in-phase and quadrature components of cluster i. Hence, the PDF of R can be obtained by averaging the conditional PDF of the α-μ fading process over the statistics of the IGA shadowing RV. Mathematically,
f R ( r ) = 0 f R | Z ( r | z ) f Z ( z ) d z ,
where f R | Z ( r | z ) is the conditional PDF of the α-μ fading distribution in [36] with respect to RV z, as follows
f R | Z ( r | z ) = α β 0.5 α μ r α μ 1 Γ ( μ ) ( z Ω ) 0.5 α μ exp [ β 0.5 α r α ( z Ω ) 0.5 α ] ,
where E [ r 2 ] = Ω , β   = Γ ( μ + 2 α ) / Γ ( μ ) represents the power scale factor, E [ ]   denoting the statistical mean, and Γ ( · ) is the gamma function defined in ([37], Equation (8.310.1)). In (2), f Z ( z ) is the PDF of the IGA distribution to model shadowing fading in [15], and can be given by
f Z ( z ) = ( n 1 ) n Γ ( n ) z ( n + 1 ) exp ( n 1 z ) ,
where n denotes the shape parameter and n > 1 , E [ z ] = 1 . Note that n 1 denotes the received signals are subject to the heavy shadowing; on the contrary, n represents the absence of shadowing. Substituting (3) and (4) in (2), the PDF of the received signal envelope R can be expressed as
f R ( r ) = α β 0.5 α μ ( n 1 ) n r α μ 1 Γ ( μ ) Γ ( n ) Ω 0.5 α μ 0 z ( n + 0.5 α μ + 1 ) exp ( n 1 z ) exp [ β 0.5 α r α ( z Ω ) 0.5 α ] d y .
To obtain the closed-form expression of the PDF of the composite envelope R, we represent the exponential form exp ( x ) in terms of Fox’s H function H 0 , 1 1 , 0 [ x | ( 0 , 1 ) ]   by using the identity ([38], Equation (1.125)). After a simple variable transformation and based on the definition of Fox’s H function in [38] (Equation (1.2)), (5) can be rewritten as
f R ( r ) = α r 1 Γ ( μ ) Γ ( n ) H 1 , 1 1 , 1 [ β 0.5 α r α [ ( n 1 ) Ω ] 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) ] ,
where H p , q m , n [ · | · ] denotes a univariate Fox’s H-function. Similar to the η-μ/IGA and κ-μ/IGA models, the α-μ/IGA composite distribution is also a generalized composite distribution. As α = 2, this model is equivalent to the Fisher-Snedecor F composite fading model proposed in [28]. As μ = 1, the α-μ/IGA distribution reduces to Weibull/IGA. By setting α = 2 and μ = 1, this model results in Rayleigh/IGA, whilst n , the α-μ model becomes its special case. Here, we let the instantaneous SNR γ = r 2 E s / N 0 , where E s is the average energy per symbol and N 0 denotes the single-sided power spectral density, then the average SNR can be given as γ ¯ =   Ω E s / N 0 . Thus, by using (6) and applying the change of variables, the corresponding PDF of γ can be obtained as
f γ ( γ ) = α γ 1 2 Γ ( μ ) Γ ( n ) H 1 , 1 1 , 1 [ λ γ 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) ] ,
where λ = ( β / ( n 1 ) γ ¯ ) 0.5 α . Obviously, (7) can reduce to [39] (Equation (2)) as α = 2 by using the identity [38] (Equation (1.60)). Invoking the definition of the univariate Fox’s H-function and using the similar derived process of (6), the corresponding CDF and MGF of γ can be obtained, respectively, as
F γ ( γ ) = 0 γ f γ ( γ ) d γ = 1 Γ ( μ ) Γ ( n ) H 2 , 2 1 , 2 [ λ γ 0.5 α | ( μ , 1 ) , ( 0 , 1 ) ( 1 n , 0.5 α ) , ( 1 , 1 ) ] ,
M G F γ ( s ) = 0 exp ( s γ ) f γ ( γ ) d γ = α 2 Γ ( μ ) Γ ( n ) H 2 , 1 1 , 2 [ λ s 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) , ( 1 , 0.5 α ) ] .
It is noted here that (8) and (9) can also reduce to [40] (Equation (2)), [41] (Equation (2)) as α = 2, respectively. Interestingly, the Fisher-Snedecor F model is also a special case of the α- F distribution. In [32], the α- F model was deduced from the F model by using a power transformation, and its signal envelope R is expressed as R α = i = 1 μ Z ( X i 2 + Y i 2 ) . From the right-hand side of (1), the α-μ/IGA and α- F models have the identical mathematical form and physical parameters including α, μ, and n. However, they have different physical interpretations from the view of the physical processes of the signal envelopes. The former represents the mean of the non-linear multipath envelopes is subject to the shadowing effects following the IGA distribution, whereas the latter can be considered as the composite envelopes following the F distribution encounter the non-linear propagation in a nonhomogeneous environment.
For the sake of comparison, we evaluated the impacts of the physical parameters on the PDFs of both models in Figure 1. It can be observed from Figure 1 that as α = 2, both models reduce to the F one, while as α ≠ 2 and n = 5 (moderate shadowing), they diverge greatly from each other. With the increase in the value of shadowing parameters (n = 5→50), they tend to have similar results, which indicates the shadowing effects decrease, and they reduce to the α-μ model. Moreover, the parameter μ shows the same effect on the PDFs of both models when μ grows. This is because the nonlinear feature of the α- F model affects both μ and n, while only μ is affected by its nonlinear parameter in the α-μ/IGA model. As expected, the PDF curves of both models converge to around the average SNR ( γ ¯ = 1 ) when the values of α and μ get larger, which means good channel conditions.
Recently, several mixed PDF models were presented to provide some accurate approximated expressions via the Gaussian-Laguerre quadrature approximation and overcome certain mathematical difficulties, such as the mixture Gamma (MG) distribution [41] and the mixture inverse Gaussian distribution [7]. By applying similar approaches, the exact closed-form expression of (6) can also be approximated as a mixed form of the IGA or the α-μ distribution. In [26], the mixture IGA model has been considered to approximate the α-μ/IGA model just as the MG distribution was used to approximate the α-μ/GA model in [6]. However, it was reported in [26] that the number of the truncation term is large and goes up to 45 in order to achieve the error <10−5. To this end, we consider an alternative new form to approximate the α-μ/IGA model. By letting t = ( n 1 ) / z in (5) and performing some straightforward manipulations and several the changes of variables, (7) can be approximated as
f γ ( γ ) = i = 1 N ς i α γ 0.5 α μ 1 2 Γ ( μ ) ψ i 0.5 α μ exp [ γ 0.5 α ψ i 0.5 α ] ,
where ς i = ω i t i μ 1 / Γ ( n ) , ψ i = γ ¯ ( n 1 ) / ( β t i ) , N denotes the number of truncation terms, and ω i and t i are the abscissas and weight factors for the Gaussian-Laguerre integration, respectively. Here, we name the novel form in (10) as a mixed α-μ distribution model. In particular, this model is a generalized version of the MG distribution and can be referred to as the generalized MG (GMG) distribution. Most noteworthy is that the α-X multipath and α-μ-Y composited fading distributions can be simplified as this mixture sum form, where X denotes η-μ or κ-μ, and Y denotes LN, Gamma, or α-μ [11] distribution models.
Based on the definition in (8) and (9), and using (10), the corresponding approximated forms of the CDF and the MGF of γ can be obtained, respectively, as
F γ ( γ ) = i = 1 N ς i γ ( μ , ( γ / ψ i ) 0.5 α ) Γ ( μ ) ,
M G F γ ( s ) = i = 1 N ς i α 2 Γ ( μ ) H 1 , 1 1 , 1 [ ( s ψ i ) 0.5 α | ( μ , 1 ) ( 1 , 0.5 α ) ] ,
where γ ( · , · ) is the lower incomplete gamma function defined in [37] (Equation (8.310.1)). In order to prove the accuracy of the approximated distribution models, two classical criteria of error are often adopted to evaluate the error between the exact PDF and the approximated PDF in [41], namely the mean square error (MSE) and the Kullback-Leibler divergence (KLD). For the purposes of comparison, the approximated PDF of   γ based on the mixed IGA model can be expressed as [26]
f γ ( γ ) = i = 1 N ξ i η i n Γ ( n ) γ n 1 exp [ η i γ 1 ] ,
where ξ i = ω i t i μ 1 / Γ ( μ ) , η i = γ ¯ ( n 1 ) ( t i / β ) 2 / α .
In Figure 2a–d, we compare the exact PDF in (7), the mixed α-μ PDF in (10), and the mixed IGA PDF in (13) under different conditions, and present the impacts of various parameters on the PDFs. At the same time, the Monte Carlo simulations are also provided to validate the accuracy of the approximated distributions. From Figure 2a–d, it can be observed that the PDF curves of the mixed IGA model show the larger divergences when the values of the truncated terms (N), α, and/or μ are smaller (for example, N = 8, 12 in Figure 2a, α = 1.5 in Figure 2b, and μ = 1.5 in Figure 2d), or the values of the shadowing parameter, n, become larger (for example, n = 20 in Figure 2c), whereas it is interesting that the PDF curves of the mixed α-μ achieve good agreements with the exact PDF in (7) in most of the cases except the case where α = 5 and N = 8 in Figure 2b. Although increasing the value of N can reduce the divergences between the exact PDF and the mixed IGA, it renders more computation complexity and a longer running time. Furthermore, with the aid of the definitions of MSE and KLD in [42], the error analytical results are calculated under various fading scenarios in Table 1. From Table 1, we can also obtain the same conclusions as ones from Figure 2a–d. Despite the mixed α-μ distribution being used to simplify the α-κ-μ shadowed model in [43], its fading parameters (μ and m) are constrained to be positive integers and cannot be extended to other models, as well as having no application for it. To the best of the authors’ knowledge, the mixed α-μ PDF in (10) is novel and has not yet been reported in the open technical papers. Based on the above results, the mixed α-μ distribution will be discussed in the following sections.

2.2. Bivariate Case

In the performance evaluations of wireless digital communication systems, the bivariate distributions are often used to analyze the performance of diversity receivers over correlated fading environments, such as [1,29]. So far, the bivariate α-μ/IGA composited distribution has not been considered in the previous works. In this context, we will derive the analytical expressions of the statistical properties of the bivariate α-μ/IGA distribution.
Let Ri (i = 1, 2) represent the channel fading envelope of α-μ processes, then, the joint conditional PDF between R1 and R2 can be expressed from [36] (Equation (28)) as
f R 1 | W 1 , R 2 | W 2 ( r 1 | w 1 , r 2 | w 1 ) = α 1 α 2 r 1 0.5 α 1 ( μ + 1 ) 1 r 2 0.5 α 2 ( μ + 1 ) 1 ρ N ( μ 1 ) / 2 Γ ( μ ) ( ( β 1 / w 1 Ω 1 ) 0.5 α 1 ( β 2 / w 2 Ω 2 ) 0.5 α 2 ) ( μ + 1 ) ( 1 ρ N ) × exp [ 1 1 ρ N ( r 1 α 1 β 1 0.5 α 1 ( w 1 Ω 1 ) 0.5 α 1 + r 2 α 2 β 2 0.5 α 2 ( w 2 Ω 2 ) 0.5 α 2 ) ] I μ 1 [ 2 ( 1 ρ N ) ρ N r 1 α 1 r 2 α 2 β 1 0.5 α 1 β 2 0.5 α 2 ( w 1 Ω 1 ) 0.5 α 1 ( w 2 Ω 2 ) 0.5 α 2 ] ,
where α i is the nonlinear fading parameter, ρ N is the power correlation coefficient between R 1 2 and R 2 2 , and I μ 1 ( x ) represents the modified Bessel function of the first kind and order (μ − 1) defined in [36] (Equation (9.210/1)), β i =   Γ ( μ + 2 / α i ) / Γ ( μ ) . When multipath fading is superimposed on shadowing, W i slowly varies. Here, we assume W i follows the IGA distribution in (4) and E [ r i 2 ] = Ω i . By using a simple operation of Equation (8) given in [44], the joint PDF of the bivariate IGA distribution can be derived as
f W 1 , W 2 ( w 1 , w 2 ) = ( n 1 ) n + 1 ( w 1 w 2 ) 0.5 ( n + 3 ) Γ ( n ) ( 1 ρ G ) ρ G ( n 1 ) / 2 exp [ ( n 1 ) ( w 1 1 + w 2 1 ) 1 ρ G ] I n 1 [ 2 ( n 1 ) ρ G w 1 w 2 ( 1 ρ G ) ] ,
where n > 1   is the IGA shaping parameter, ρ G is the power correlation coefficient between w 1 2 and w 2 2 , and E [ w i ] = 1.
Thus, by using the same approach shown in [29] (Equation (5)), and performing some basic operations, the joint PDF of the bivariate α-μ/IGA composite distribution is written as
f R 1 , R 2 ( r 1 , r 2 ) = k = 0 l = 0 Φ ρ N k ρ G l r 1 1 r 2 1 k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 1 , 1 1 , 1 [ Ξ i r i α i | ( μ + k , 1 ) ( 1 n l , 0.5 α i ) ] ,
where Φ = α 1 α 2 ( 1 ρ G ) n ( 1 ρ N ) μ / [ Γ ( μ ) Γ ( n ) ] , Ξ i = ( ( 1 ρ G ) β i ) 0.5 α i / [ ( 1 ρ N ) ( Ω i ( n 1 ) ) 0.5 α i ] .
By setting γ i = r i 2 E s / N 0 and γ ¯ i = Ω i E s / N 0 , after a simple variable transformation, the joint PDF of γ 1 and γ 2 over the correlated α-μ/IGA channels can be obtained as
f γ 1 , γ 2 ( γ 1 , γ 2 ) = k = 0 l = 0 0.25 Φ ρ N k ρ G l γ 1 1 γ 2 1 k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 1 , 1 1 , 1 [ ϑ i γ i 0.5 α i | ( μ + k , 1 ) ( 1 n l , 0.5 α i ) ] ,
where ϑ i = ( ( 1 ρ G ) β i ) 0.5 α i / [ ( 1 ρ N ) ( ( n 1 ) γ ¯ i ) 0.5 α i ] . Hence, with the aid of [37] (Equation (3.194.1)), the corresponding joint CDF of   γ 1 and γ 2 can be yielded as
F γ 1 , γ 2 ( γ 1 , γ 2 ) = 0 γ 1 0 γ 2 f γ 1 , γ 2 ( γ 1 , γ 2 ) d γ 1 d γ 2   = k = 0 l = 0 0.25 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 2 , 2 1 , 2 [ ϑ i γ i 0.5 α i | ( μ + k , 1 ) , ( 0 , 0.5 α i ) ( 1 , 0.5 α i ) , ( 1 n l , 0.5 α i ) ] .
Moreover, based on the definition of the joint MGF in [1], the joint MGF of γ 1 and γ 2 can be given by
M γ 1 , γ 2 ( s 1 , s 2 ) = E [ exp ( s 1 γ 1 s 2 γ 2 ) ]   = k = 0 l = 0 0.25 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 2 , 1 1 , 2 [ ϑ i s i 0.5 α i | ( μ + k , 1 ) ( 1 , 0.5 α i ) , ( 1 n l , 0.5 α i ) ] .

3. Statistical Characteristics of the Sum of α-μ/IGA RVs

In this section, the statistical properties of the sum of L α-μ/IGA RVs are investigated. The classical application of the sum of RVs in diversity systems is to evaluate the performance of MRC and EGC systems, and only the MRC system is considered here. To this effect, we first derive the exact closed-form expressions of PDF, CDF, and MGF of the instantaneous output SNR for the MRC system in the i.n.i.d. case. Then, their corresponding approximated expressions are obtained based on the mixture α-μ model, respectively. Finally, their correlated statistics are acquired when L = 2.

3.1. Exact Statistics of the Sum of the Independent RVs

Let γ M R C = l = 1 L γ l , where γ l is the i.n.i.d. α-μ/IGA RV with the PDF in (7), the average SNR γ l ¯ , the fading parameters α l , μ l , and n l , l = 1 , ,   L , and L is the number of the RVs, then the PDF of γ M R C is yielded by
f γ M R C ( γ ) = Θ 1 γ 1 H 0 , 1 : [ 2 , 1 ] l = 1 : L 0 , 0 : [ 1 , 2 ] l = 1 : L [ λ 1 γ 0.5 α 1 , , λ L γ 0.5 α L | ( 1 ; { 0.5 α l } l = 1 : L ) ( ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] l = 1 : L ] ,
where Θ 1 = l = 1 L α l / [ 2 Γ ( μ l ) Γ ( n l ) ] , H p , q : p 1 , q 1 : : p L , q L m , n : m 1 , n 1 : : m L , n L [ · ] denotes a multivariate Fox’s H-function defined in [38] (Equation (A.1)),   λ l = ( β l / ( n l 1 ) γ ¯ l ) 0.5 α l , β l = Γ ( μ l + 2 α l ) / Γ ( μ l ) .
Proof: Please refer to Appendix A.
Accordingly, based on the definition of the CDF and performing one definite integral of (20), the corresponding CDF of the sum of L α-μ/IGA RVs, γ M R C , can be expressed as
F γ M R C ( γ ) = Θ 1 H 0 , 1 : [ 2 , 1 ] l = 1 : L 0 , 0 : [ 1 , 2 ] l = 1 : L [ λ 1 γ 0.5 α 1 , , λ L γ 0.5 α L | ( 0 ; { 0.5 α l } l = 1 : L ) ( ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] l = 1 : L ] .
For the i.i.d. case, (20) and (21) can be simplified by omitting the subscript.

3.2. Approximated Statistics of the Sum of the Independent RVs

By using (A1) and (12), the approximated MGF of γ M R C can be obtained as
M G F γ M R C ( s ) = l = 1 L M G F γ l ( s ) = l = 1 L ( i l = 1 N l ς i l α l 2 Γ ( μ l ) H 1 , 1 1 , 1 [ ( s ψ i l ) 0.5 α l | ( μ l , 1 ) ( 1 , 0.5 α l ) ] ) .
For ease of calculation, we swap the order of the sum and the product in the right-hand side formula in (22); (22) can then be re-expressed in multiple summation forms as
M G F γ M R C ( s ) = i 1 = 1 N 1 i L = 1 N L l = 1 L ( ς i l α l 2 Γ ( μ l ) H 1 , 1 1 , 1 [ ( s ψ i l ) 0.5 α l | ( μ l , 1 ) ( 1 , 0.5 α l ) ] )
By utilizing (23) and the similar derivation procedure as Appendix A, the approximated PDF of γ M R C can be obtained as
f γ M R C ( γ ) = i 1 = 1 N 1 i L = 1 N L [ ( l = 1 L ς i l α l 2 Γ ( μ l ) ) γ 1 H 0 , 1 : [ 1 , 1 ] l = 1 : L 0 , 0 : [ 1 , 1 ] l = 1 : L [ ( γ / ψ i 1 ) 0.5 α 1 , , ( γ / ψ i L ) 0.5 α L | ( 1 ; { 0.5 α l } l = 1 : L ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 , 0.5 α l ) ] l = 1 : L ] ] .
From (21) and (24), we can see the summation forms in (24) bring the reduction of arguments in the multivariate Fox’s H-function. Importantly, (24) is a unified analytical expression and can be applied to more channel models where the PDF can be expressed as a mixture α-μ form in (10), such as α-η-μ and α-κ-μ, α-η-μ/lognormal and α-κ-μ/lognormal in [5], α-μ/gamma in [6], and α-μ/α-μ in [12]. Hence, the corresponding approximated CDF of γ M R C can be derived as
F γ M R C ( γ ) = i 1 = 1 N 1 i L = 1 N L [ ( l = 1 L ς i l α l 2 Γ ( μ l ) ) H 0 , 1 : [ 1 , 1 ] l = 1 : L 0 , 0 : [ 1 , 1 ] l = 1 : L [ ( γ / ψ i 1 ) 0.5 α 1 , , ( γ / ψ i L ) 0.5 α L | ( 0 ; { 0.5 α l } l = 1 : L ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 , 0.5 α l ) ] l = 1 : L ] ] .
In the i.i.d. case, (22) can be simplified as
M G F γ M R C ( s ) = ( i = 1 N ς i α 2 Γ ( μ ) H 1 , 1 1 , 1 [ ( s ψ i ) 0.5 α | ( μ , 1 ) ( 1 , 0.5 α ) ] ) L
By applying the multinomial theorem, (26) can be rewritten as
M G F γ M R C ( s ) = k 1 + + k N = L ( L k 1 , , k N ) 1 i N ( ς i α 2 Γ ( μ ) H 1 , 1 1 , 1 [ ( s ψ i ) 0.5 α | ( μ , 1 ) ( 1 , 0.5 α ) ] ) k i   I 1 ,
where k 1 + + k N = L = k 1 = 0 L k 2 = 0 L k 1 k N 1 = 0 L k 1 k N 2 , ( L k 1 , , k N ) = L ! k 1 ! k 2 ! k N ! .
To find the solution in (27), I1 in (27) can be expanded as
I 1 = ( α 2 Γ ( μ ) ) L i = 1 N ( ς i ) k i i = 1 N ( H 1 , 1 1 , 1 [ ( s ψ i ) 0.5 α | ( μ , 1 ) ( 1 , 0.5 α ) ] ) k i   I 2
With the help of the definition of the univariate Fox H-function, I2 in (28) can be expressed as
I 2 = ( 1 2 π j ) L 11 1 k 1 i 1 i k i N k N ( i = 1 N j = 1 k i Γ ( μ t i j ) Γ ( 0.5 α t i j ) ( s ψ i ) 0.5 α t i j ) d t 11 d t N k N = H 0 , 0 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N 0 , 0 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N [ ( s ψ 1 ) 0.5 α , , ( s ψ 1 ) 0.5 α   k 1 , , ( s ψ N ) 0.5 α | | [ ( μ , 1 ) ] l = 1 : k 1 [ ( 1 , 0.5 α ) ] l = 1 : k 1 | | [ ( μ , 1 ) ] l = 1 : k N [ ( 1 , 0.5 α ) ] l = 1 : k N ] .
By inserting (29) and (28) into (27), the closed-form expression of the MGF in (27) can be readily obtained in the i.i.d. case. Thus, by using (A3), after some mathematical manipulations, the approximated PDF of γ M R C in the i.i.d. case can be derived as
f γ M R C ( γ ) = k 1 + + k N = L ( L k 1 , , k N ) ( α 2 Γ ( μ ) ) L ( i = 1 N ( ς i ) k i ) γ 1   × H 0 , 1 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N 0 , 0 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N [ ( γ / ψ 1 ) 0.5 α , , ( γ / ψ 1 ) 0.5 α   k 1 , , ( γ / ψ N ) 0.5 α | ( 1 ; { 0.5 α } l = 1 : L ) | [ ( μ , 1 ) ] l = 1 : k 1 [ ( 1 , 0.5 α ) ] l = 1 : k 1 | | [ ( μ , 1 ) ] l = 1 : k N [ ( 1 , 0.5 α ) ] l = 1 : k N ] .
Analogous to (25), the approximated CDF of γ M R C in the i.i.d. case can be derived as
F γ M R C ( γ ) = k 1 + + k N = L ( L k 1 , , k N ) ( α 2 Γ ( μ ) ) L ( i = 1 N ( ς i ) k i )   × H 0 , 1 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N 0 , 0 : [ 1 , 1 ] l = 1 : k 1 [ 1 , 1 ] l = 1 : k N [ ( γ / ψ 1 ) 0.5 α , , ( γ / ψ 1 ) 0.5 α   k 1 , , ( γ / ψ N ) 0.5 α | ( 0 ; { 0.5 α } l = 1 : L ) | [ ( μ , 1 ) ] l = 1 : k 1 [ ( 1 , 0.5 α ) ] l = 1 : k 1 | | [ ( μ , 1 ) ] l = 1 : k N [ ( 1 , 0.5 α ) ] l = 1 : k N ] .

3.3. Statistics of the Two Correlated RVs

For a dual-branch MRC diversity receiver, the instantaneous output SNR per symbol is expressed as   γ M R C = γ 1 + γ 2   in [1]. Under a correlated fading environment, it is difficult to directly find a simple and closed-form expression of the PDF of   γ M R C . To this end, we considered the MGF-based approach to derive the PDF of   γ M R C   over correlated α-μ/IGA composite fading channels. Based on (19), the MGF of γ M R C can be written as
M γ M R C ( s ) = M γ 1 , γ 2 ( s , s )   = k = 0 l = 0 0.25 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 2 , 1 1 , 2 [ ϑ i s 0.5 α i | ( μ + k , 1 ) ( 1 , 0.5 α i ) , ( 1 n l , 0.5 α i ) ] .
By adopting similar steps in Appendix A and carrying out some mathematical manipulations, the PDF of γ M R C can be derived as
f γ M R C ( γ ) = k = 0 l = 0 0.25 Φ ρ N k ρ G l γ 1 k ! l ! Γ ( k + μ ) Γ ( l + n ) × H 0 , 1 : 2 , 1 : 2 , 1 0 , 0 : 1 , 2 : 1 , 2 [ ϑ 1 γ 0.5 α 1 , ϑ 2 γ 0.5 α 2 | ( 1 , 0.5 α 1 , 0.5 α 2 ) | ( μ + k , 1 ) ( 1 , 0.5 α 1 ) , ( 1 n l , 0.5 α 1 ) | ( μ + k , 1 ) ( 1 , 0.5 α 2 ) , ( 1 n l , 0.5 α 2 ) ] .
Then, the corresponding CDF of γ M R C can be yielded as
F γ M R C ( γ ) = k = 0 l = 0 0.25 Φ ρ N k ρ G l γ 1 k ! l ! Γ ( k + μ ) Γ ( l + n ) × H 0 , 1 : 2 , 1 : 2 , 1 0 , 0 : 1 , 2 : 1 , 2 [ ϑ 1 γ 0.5 α 1 , ϑ 2 γ 0.5 α 2 | ( 1 , 0.5 α 1 , 0.5 α 2 ) | ( μ + k , 1 ) ( 1 , 0.5 α 1 ) , ( 1 n l , 0.5 α 1 ) | ( μ + k , 1 ) ( 1 , 0.5 α 2 ) , ( 1 n l , 0.5 α 2 ) ] .

4. Performance Analysis

Capitalizing on the previous derived statistical characteristics of the single-link and MRC systems, some performance metrics of interest, such as the outage probability, the ABEP/ASEP, the ergodic capacity, and the effective capacity are evaluated and discussed in this section, respectively.

4.1. Outage Probability

The OP is an important performance metric of wireless communications systems operating over fading channels. It is defined as the probability that the instantaneous SNR at the receiver output, γ, falls below a predefined outage threshold, γth. Based on this definition, the outage probability of a wireless communication system over fading channels can be expressed as
P o u t = Pr ( 0 < γ < γ t h ) = 0 γ t h f γ ( γ ) d γ .
Therefore, the corresponding OP is readily deduced as follows
P o u t = F γ ( γ t h ) ,
where F γ ( γ t h ) can be obtained by using (8), (11), (21), (25), (31), and (34) replacing γ with γ𝑡 over different fading conditions, respectively.

4.2. Average BEP/SEP

ABEP/ASEP is a conventional important measure to reveal the wireless communications system behaviors over fading environments. In the previous technical literature, two classical methods are often used to evaluate the ABEP/ASEP of various digital modulation schemes. One is the PDF-based approach, and the other is the MGF-based approach. In this paper, the former is considered and the unified conditional BEP/SEP formulas for several classical modulation schemes, such as DPSK, BPSK (Binary Phase Shift Keying), MPSK (M-ary PSK), and MQAM (Multiple Quadrature Amplitude Modulation), are summarized in terms of the exponential function form and the complementary error function form in Table 2, where e r f c ( ) is the complementary error function defined in [37] (Equation (8.250.4)).
As is well-known, the ABEP/ASEP in fading environments can be obtained by averaging the conditional BEP/SEP ( P e ( ε | γ ) in Table 2) over the PDF ( f γ ( γ ) ) of the instantaneous SNR ( γ ), namely
P e ( ε ) = 0 P e ( ε | γ ) f γ ( γ ) d γ
Consequently, the ABEP/ASEP of the single-link and MRC systems in Section 2 and Section 3 can be deduced by substituting their corresponding PDFs into (37).

4.2.1. ABEP/ASEP in the Exponential Form

From Table 2, we can see the conditional BEP expressions of DPSK are written in the exponential form. For this case, the ABEP expressions for the single-link and MRC systems over independent and correlated fading channels can be derived by using the approaches to find the MGF of γ , such as (9). To this end, the exact ABEP expressions for various systems are obtained in Table 3.

4.2.2. ABEP/ASEP in the Erfc Function Form

A. Single-Link System

By inserting Pe ( ε | γ ) in the erfc function form and (7) into (37), the ABEP/ASEP of the single-link system over α-μ/IGA fading channels can be given as
P e ( ε ) = α C 2 Γ ( μ ) Γ ( n ) 0 γ 1 e r f c ( E γ ) H 1 , 1 1 , 1 [ λ γ 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) ] d γ   I 3   α D 2 Γ ( μ ) Γ ( n ) 0 γ 1 e r f c 2 ( E γ ) H 1 , 1 1 , 1 [ λ γ 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) ] d γ   I 4 .
In order to solve the integrals (I3 and I4) in (38), we employ the identities [45] (Equations (8.4.14.2) and (8.3.2.21)) to represent the erfc function in terms of Fox’s H function as
e r f c ( E γ ) = π 0.5 H 1 , 2 2 , 0 [ E γ | ( 0 , 1 ) , ( 0.5 , 1 ) ( 1 , 1 ) ] .
By inserting (39) into I3 and I4, with the help of [46] (Equation (2.8.4) and [47] (Equation (2.3)) and performing some algebraic operations, I3 and I4 can be obtained, respectively, as
I 3 = π 0.5 H 3 , 2 1 , 3 [ λ E 0.5 α | ( μ , 1 ) , ( 0 , 0.5 α ) ( 1 n , 0.5 α ) , ( 1 , 0.5 α ) , ( 0.5 , 0.5 α ) ] ,
I 4 = π 1 H 2 , 1 : 1 , 2 : 1 , 1 0 , 2 : 2 , 0 : 1 , 1 [ 1 , λ E 0.5 α | ( 0 , 1 , 0.5 α ) ( 1 , 1 , 0.5 α ) , ( 0.5 , 1 , 0.5 α ) | ( 0 , 1 ) , ( 0.5 , 1 ) ( 1 , 1 ) | ( μ , 1 ) ( 1 n , 0.5 α ) ] .
Consequently, if we substitute (40) into (38), the exact ABEP/ASEP expression of the single-link system over α-μ/IGA fading channels can be obtained. Following the similar procedure used in (40) and utilizing (10), the exact approximated ABEP/ASEP expressions of the single-link system by using the mixed α-μ model can be deduced as
P e ( ε ) = i = 1 N C α π 0.5 ς i 2 Γ ( μ ) H 2 , 2 1 , 2 [ ( ψ i E ) 0.5 α | ( μ , 1 ) , ( 0 , 0.5 α ) ( 1 , 0.5 α ) , ( 0.5 , 0.5 α ) ]   i = 1 N D α π 1 ς i 2 Γ ( μ ) H 2 , 1 : 1 , 2 : 0 , 1 0 , 2 : 2 , 0 : 1 , 0 [ 1 , ( ψ i E ) 0.5 α | ( 0 ; 1 , 0.5 α ) ( 1 ; 1 , 0.5 α ) , ( 0.5 ; 1 , 0.5 α ) | ( 0 , 1 ) , ( 0.5 , 1 ) ( 1 , 1 ) | ( μ , 1 ) ]

B. MRC System

In the i.n.i.d. case, by inserting (20) into (37), a unified ABEP/ASEP formula in the erfc function form for an L-branch MRC system over α-μ/IGA fading channels can be expressed as
P e ( ε ) = C π 0.5 Θ 1 H 1 , 1 : [ 2 , 1 ] l = 1 : L 0 , 1 : [ 1 , 2 ] l = 1 : L [ λ 1 E 0.5 α 1 , , λ L E 0.5 α L | ( 0 ; { 0.5 α l } l = 1 : L ) ( 0.5 ; { 0.5 α l } l = 1 : L ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] l = 1 : L ]   D π 1 Θ 1 H 2 , 2 : [ 2 , 1 ] l = 1 : L : 1 , 2 0 , 2 : [ 1 , 2 ] l = 1 : L : 2 , 0 [ λ 1 E 0.5 α 1 , , λ L E 0.5 α L , 1 | ( 1 ; { 0.5 α l } l = 1 : L , 0 ) , ( 0 ; { 0.5 α l } l = 1 : L , 1 ) ( 1 ; { 0.5 α l } l = 1 : L , 1 ) , ( 0.5 ; { 0.5 α l } l = 1 : L , 1 ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] l = 1 : L | ( 0 , 1 ) , ( 0.5 , 1 ) ( 1 , 1 ) ] .
Proof: Please refer to Appendix B.
Similarly, by inserting (24) into (37), a unified ABEP/ASEP formula in the erfc function form for an L-branch MRC system by using the mixed α-μ model can also be obtained. In the correlated case, by inserting (33) into (37), a unified expression of ABEP/ASEP in the erfc function form for a two-branch MRC system over α-μ/IGA fading channels can be expressed as
P e ( ε ) = k = 0 l = 0 0.25 C π 0.5 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n ) H 1 , 1 : 2 , 1 : 2 , 1 0 , 1 : 1 , 2 : 1 , 2 [ ϑ 1 E 0.5 α 1 , ϑ 2 E 0.5 α 2 | ( 0 ; 0.5 α 1 , 0.5 α 2 ) ( 0.5 ; 0.5 α 1 , 0.5 α 2 ) | ( μ + k , 1 ) ( 1 , 0.5 α 1 ) , ( 1 n l , 0.5 α 1 ) | ( μ + k , 1 ) ( 1 , 0.5 α 2 ) , ( 1 n l , 0.5 α 2 ) ]   k = 0 l = 0 0.25 D π 1 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n )   × H 2 , 2 : 2 , 1 : 2 , 1 : 1 , 2 0 , 2 : 1 , 2 : 1 , 2 : 2 , 0 [ ϑ 1 E 0.5 α 1 , ϑ 2 E 0.5 α 2 , 1 | ( 1 ; 0.5 α 1 , 0.5 α 2 , 0 ) , ( 0 ; 0.5 α 1 , 0.5 α 2 , 1 ) ( 1 ; 0.5 α 1 , 0.5 α 2 , 1 ) , ( 0.5 ; 0.5 α 1 , 0.5 α 2 , 1 ) | ( μ + k , 1 ) ( 1 , 0.5 α 1 ) , ( 1 n l , 0.5 α 1 ) | ( μ + k , 1 ) ( 1 , 0.5 α 2 ) , ( 1 n l , 0.5 α 2 ) | ( 0 , 1 ) , ( 0.5 , 1 ) ( 1 , 1 ) ] .

4.3. Average Channel Capacity

The channel capacity, in Shannon’s sense, is a core performance measure since it provides the maximum achievable transmission rate at which the errors are recoverable. The average channel capacity per unit bandwidth in fading channels can be expressed as
C ¯ γ = B ln 2 0 ln ( 1 + γ ) f γ ( γ ) d γ ,
where f γ ( γ ) is the PDF of the instantaneous SNR, γ .

4.3.1. Single-Link System

In order to find the closed-form expression of average channel capacity for the single-link system, we first use the identities in [45] (Equations (8.4.6.5) and (8.3.2.21)) to represent l n ( 1 + γ ) in terms of Fox’s H-function as H 2 , 2 1 , 2 [ γ | ( 1 , 1 ) , ( 0 , 1 ) ( 1 , 1 ) , ( 1 , 1 ) ] , then insert (7) into (44). After some mathematical manipulations, the exact expression of the average channel capacity for the single-link system over α-μ/IGA fading channels can be given by
C ¯ = 0.5 α ln 2 Γ ( μ ) Γ ( n ) H 3 , 3 3 , 2 [ λ | ( μ , 1 ) , ( 0 , 0.5 α ) , ( 0 , 0.5 α ) ( 1 n , 0.5 α ) , ( 0 , 0.5 α ) , ( 1 , 0.5 α ) ] .
Similarly, by inserting (10) into (44), the exact approximate expression of the average channel capacity for the single-link system over α-μ/IGA fading channels can be obtained as
C ¯ = i = 1 N 0.5 ς i α ln 2 Γ ( μ ) H 2 , 3 3 , 1 [ ψ i 0.5 α | ( μ , 1 ) , ( 0 , 0.5 α ) , ( 0 , 0.5 α ) ( 0 , 0.5 α ) , ( 1 , 0.5 α ) ] .

4.3.2. MRC System

In order to obtain the analytical expression of the average channel capacity of the MRC system and circumvent a mismatch with the definition of the multivariate H-function in [48,49], the final value theorem was adopted and an infinitesimal parameter close to zero was introduced. However, the multivariate H-function in the expression of Ergodic capacity in [48] (Equation (28)) did not exactly satisfy its definition in [38] (Equation (A.1)) because its parameter term was m = 1, not m = 0 as defined in [38] (Equation (A.1)). Moreover, the expression of Ergodic capacity in [49] (Equation (33)) avoided the above problem and caused a new case in which some associated coefficients α j ( i ) and β j ( i ) defined in [49] (Equation (34)) are negative real numbers (e.g., −1), not positive real numbers as defined in [38] (Equation (A.5)). Therefore, by using the final value theorem and inserting (20) into (44) along with some mathematical manipulations, an exact expression of average channel capacity for the MRC system over α-μ/IGA fading channels can be derived as
C ¯ γ M R C = ( Θ 1 / ln 2 ) H 1 , 1 : [ 2 , 1 ] l = 1 : L : 2 , 2 0 , 1 : [ 1 , 2 ] l = 1 : L : 1 , 2 [ λ 1 s 0.5 α 1 , , λ L s 0.5 α L , s 1 | ( 1 ; { 0.5 α l } l = 1 : L , 0 ) ( 1 ; { 0.5 α l } l = 1 : L , 1 ) | [ ( μ l , 1 ) ] l = 1 : L [ ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] l = 1 : L | ( 1 , 1 ) , ( 0 , 1 ) ( 1 , 1 ) , ( 1 , 1 ) ] ,
where s is a very small positive number that goes to zero, for instance, s = 10−6. Similarly, by inserting (24) into (44), a unified expression of average channel capacity by using the mixed α-μ model can also be obtained. Unfortunately, it will take too long to calculate (47) in order to obtain the exact results in numerical analysis. To this end, an MGF-based method is also considered to evaluate the average channel capacity in this section. Thus, the average channel capacity per unit bandwidth can be given as [50] (Equation (8))
C ¯ γ M R C = ( ln 2 ) 1 0 s 1 [ 1 M G F γ M R C ( s ) ] exp ( s ) d s .
where M G F γ M R C ( s ) is given in Equation (A1) in Appendix A. For the integral term in (48), it may be difficult to find an exact closed-form expression. By using the Gauss-Legendre quadrature method, the approximated form of (48) can be re-expressed as
C ¯ γ M R C = ( ln 2 ) 1 i = 1 N ω i t i 1 ( 1 M G F γ M R C ( t i ) ) .
For the correlated two-branch MRC system, based on the similar approach used in (47), by inserting (33) into (44), the average channel capacity per unit bandwidth over the correlated α-μ/IGA fading channels can be obtained as
C ¯ γ M R C = k = 0 l = 0 0.25 ( ln 2 ) 1 Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n )   × H 1 , 1 : 1 , 2 : 1 , 2 : 2 , 2 0 , 1 : 2 , 1 : 2 , 1 : 1 , 2 [ ϑ 1 s 0.5 α 1 , ϑ 2 s 0.5 α 2 , s 1 | ( 1 , 0.5 α 1 , 0.5 α 2 , 0 ) ( 1 , 0.5 α 1 , 0.5 α 2 , 1 ) | ( m + k , 1 ) ( 1 , 0.5 α 1 ) , ( 1 n l , 0.5 α 1 ) | ( m + k , 1 ) ( 1 , 0.5 α 2 ) , ( 1 n l , 0.5 α 2 ) | ( 1 , 1 ) , ( 0 , 1 ) ( 1 , 1 ) , ( 1 , 1 ) ] .
Similarly, the approximated form of the average channel capacity per unit bandwidth over the correlated α-μ/IGA fading channels is also obtained by inserting (32) into (49).

4.4. Effective Capacity

To evaluate the delay performance of some emerging wireless real-time applications under a quality of service (QoS) constraint, such as voice over IP and augmented reality/virtual reality, the effective capacity (EC), as an alternative performance metric of interest, has gained great attention. The EC is defined as the maximum constant arrival rate that a wireless fading channel can offer in a bid to ensure the QoS requirements. For wireless fading channels, the normalized EC can be expressed as
= 1 A log 2 ( E [ ( 1 + γ ) A ] ) = 1 A log 2 ( 0 ( 1 + γ ) A f γ ( γ ) d γ   I 5 ) ,
where A = θ T B ln 2 , θ denotes the QoS exponent characterizing the delay constraints, T is the block length, B represents the system bandwidth, and   f γ ( γ ) is the PDF of the instantaneous SNR (γ).

4.4.1. Single-Link System

By substituting (7) into (51), the integral term I5 in (51) can be yielded as
I 5 = α 2 Γ ( μ ) Γ ( n ) 0 γ 1 ( 1 + γ ) A H 1 , 1 1 , 1 [ λ γ 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) ] d γ .
To solve the integral in I5, we first represent ( 1 + γ ) A as H 1 , 1 1 , 1 [ γ | ( 0 , 1 ) ( 1 A , 1 ) ] / Γ ( A ) in terms of Fox’s H-function by using the identities in [45] (Equations (8.4.2.5) and (8.3.2.21)). Then with the aid of [46] (Equation (2.8.4)), after some mathematical manipulations, I5 can be calculated as
I 5 = α 2 Γ ( μ ) Γ ( n ) Γ ( A ) H 2 , 2 2 , 2 [ λ | ( μ , 1 ) , ( A , 0.5 α ) ( 1 n , 0.5 α ) , ( 1 , 0.5 α ) ] .
Then, plugging (53) into (51), the exact expression of the effective capacity for the single-link system over α-μ/IGA fading channels can be obtained. When we substitute (10) into (51), the approximated expression of the effective capacity for the single-link system over the α-μ/IGA fading channels can be deduced as
= 1 A log 2 ( i = 1 N 0.5 α ς i Γ ( μ ) Γ ( A ) H 1 , 2 2 , 1 [ ψ i 0.5 α | ( μ , 1 ) , ( A , 0.5 α ) ( 1 , 0.5 α ) ] ) .

4.4.2. MRC System

By inserting (20) into (51) and performing some similar manipulations as (52), the integral term I5 in (51) can be derived as
I 5 = ( Θ 1 / Γ ( A ) ) H 1 , 0 : [ 1 , 2 ] l = 1 : L 0 , 1 : [ 2 , 1 ] l = 1 : L [ 1 / λ 1 , , 1 / λ L | ( 1 A ; { 0.5 α l } l = 1 : L ) | [ ( n l , 0.5 α l ) , ( 0 , 0.5 α l ) ] l = 1 : L [ ( 1 μ l , 1 ) ] l = 1 : L ] .
For the two-branch correlated case, by inserting (33) into (51), the integral term I5 in (51) can be expressed as
I 5 = k = 0 l = 0 0.25 Φ ρ N k ρ G l [ Γ ( A ) ] 1 k ! l ! Γ ( k + μ ) Γ ( l + n )   × H 1 , 0 : 1 , 2 : 1 , 2 0 , 1 : 2 , 1 : 2 , 1 [ 1 / ϑ 1 , 1 / ϑ 2 | ( 1 A , 0.5 α 1 , 0.5 α 2 ) | ( 0 , 0.5 α 1 ) , ( n + l , 0.5 α 1 ) ( 1 μ k , 1 ) | ( 0 , 0.5 α 2 ) , ( n + l , 0.5 α 2 ) ( 1 μ k , 1 ) ] .
Hence, the effective capacity for the MRC system over independent and correlated α-μ/IGA fading channels can be obtained by plugging (55) and (56) into (51), respectively. Similarly, by inserting (24) into (51), a unified expression of the effective capacity by using the mixed α-μ model can also be obtained. As an alternative approach to finding the analytical expression of the effective capacity, an MGF-based approach has been adopted [51]. Here, we employ the Gauss-Legendre quadrature method again, an alternative approximated form of (51) can be expressed as
= 1 A log 2 ( 1 Γ ( A ) 0 x A 1 exp ( x ) M G F γ ( x ) d x ) = 1 A log 2 ( 1 Γ ( A ) i = 1 N ω i t i A 1 M G F γ ( t i ) ) ,
where M G F γ ( ) can be obtained by using (9), (22), (32), and (A1).

5. Asymptotic Analysis

In the performance evaluation of wireless communication systems, the asymptotic analysis at the high average SNR level has been a popular evaluation methodology to observe how the channel parameters affect the performance and obtain more significant insights for system deployments. Generally, the diversity order and the coding gain are two performance metrics of interest to compare different modulation systems/fading models in the previous works. In high average SNR regions, the asymptotic OP or ABEP (ASEP) can be defined as P x = ( G c · γ ¯ ) G d , where G c denotes the coding gain, G d is the diversity order, and P x represents the OP or ABEP/ASEP. In what follows, we will derive the asymptotic analytical expressions of the aforementioned communication systems over α-μ/IGA fading channels by taking the OP and the ABEP of DPSK as examples.

5.1. Single-Link System

As γ ¯ ,   based on (7) and (8) with the aid of [46] (Th(1.11)), the asymptotic forms of the CDF and MGF for the single-link system can be obtained, respectively, as
F γ ( γ ) Γ ( n + 0.5 α μ ) Γ ( 1 + μ ) Γ ( n ) ( λ γ 0.5 α ) μ ,
M G F ( s ) α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) 2 Γ ( μ ) Γ ( n ) ( λ s 0.5 α ) μ .
Thus, by using (58) and (36), the asymptotic form of the OP for the single-link system can be expressed as
P o u t ( ( n 1 ) β γ t h ( Γ ( 1 + μ ) Γ ( n ) Γ ( n + 0.5 α μ ) ) 2 / α μ · γ ¯ ) 0.5 α μ .
Similarly, from Table 3, the asymptotic form of the DPSK scheme for the single-link system can be derived as
P e ( ε ) ( ( n 1 ) B β ( 2 Γ ( μ ) Γ ( n ) A α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) ) 2 / α μ · γ ¯ ) 0.5 α μ .

5.2. MRC System

To obtain some useful insights into the MRC system as γ ¯ , by employing (59) and Equation (A1) in Appendix A, the asymptotic MGF of γ M R C in the i.i.d. case can be yielded as
M G F γ M R C ( s ) ( α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) 2 Γ ( μ ) Γ ( n ) λ μ ) L s 0.5 L α μ .
Then, with the aid of the methods used in Equation (A3) in Appendix A, the asymptotic PDF of γ M R C in the i.i.d. case can be obtained as
f γ M R C ( γ ) ( α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) 2 Γ ( μ ) Γ ( n ) λ μ ) L γ 0.5 L α μ 1 Γ ( 0.5 L α μ ) .
Thus, by integrating (63) with respect to γ and using (36), the corresponding asymptotic OP of the MRC system, when γ ¯ , can be easily derived as
P o u t ( ( 2 Γ ( μ ) Γ ( n ) [ Γ ( 1 + 0.5 L α μ ) ] 1 / L α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) ) 2 / α μ ( n 1 ) β γ t h · γ ¯ ) 0.5 L α μ .
Similarly, by using (62), the asymptotic ABEP of DPSK for the MRC system can be obtained as
P e ( ε ) ( B ( n 1 ) β A 2 / L α μ ( 2 Γ ( μ ) Γ ( n ) α Γ ( n + 0.5 α μ ) Γ ( 0.5 α μ ) ) 2 / α μ · γ ¯ ) 0.5 L α μ .
From the above asymptotic analysis, we can see that the diversity order of the single-link system is 0.5 α μ , and that of MRC systems is 0.5L α μ . As expected, the diversity systems can bring more benefits than the single-link system by the increase in the number of diversity branches. All the channel fading parameters show the impacts on the coding gain, while the diversity order depends on only the non-linear and the multipath parameters. Likewise, we can also derive the asymptotic expressions of the correlated cases and other modulation schemes and obtain some identical conclusions. In addition, the asymptotic results of the effective capacity can also be obtained at a high average SNR level by substituting (59) and (62) into (57).

6. Numerical and Simulation Results

To corroborate the derived analytical expressions in the previous sections, this section presents various numerical and simulation results under different α-μ/IGA fading channels. We adopted the Monte Carlo simulations to generate no less than 10e6 iterations. For the correlated scenarios, the simulation approaches described in [29] are considered by using a simple change of variables. In all figures, the simulation results provide a perfect matching with the numerical counterparts and confirm the availability of our theoretical derivations. In the numerical analysis, we provide some necessary methods to evaluate multivariate Fox’s H-functions and multiple-fold infinite/finite series sum in order to obtain accurate results. For the former, a Python code provided in [52] has been revised to compute them, such as (42) and (47). For the two-fold infinite series sum forms in the correlated cases, such as (43) and (50), we refer to the minimum numbers of truncated terms provided in [29] to meet the given target accuracy. Noted that kmin = lmin = 30 is considered for various values of the nonlinear parameter (α = 0.5, 1.5, 3.5) because it has similar effects on system performance as the multipath parameter (μ). Whereas for the finite series sum forms in (49), the number of terms (N) gradually increases as the average SNR ( γ ¯ ) grows, for example, N = 30 as γ ¯ = 10 dB and N = 150 as γ ¯ = 20 dB Even though the number of terms (N) is so large, it takes much less time to calculate (49) than (47) in order to get the same results. Finally, we consider some comparisons and discussions about the impacts of the non-linear parameter (α), the correlation coefficient ( ρ ), and the diversity branch (L) on the performance of different communication systems since the impacts of the multipath parameter (μ) and the shadowing parameter (n) on the performance have been profusely discussed in the previous works, such as [29,40,41,49].
In Figure 3 and Figure 4, we plot the OP of the single-link and MRC systems as a function of γ ¯ over α-μ/IGA fading channels with n = 5, μ = 2, and γ t h = 5 dB, respectively. As anticipated, it can be observed that the OP of all the systems is improved when γ ¯ increases. In Figure 3, some performance comparisons are given in the i.i.d. case. It is clear that the OP gets better with the increase in the number of diversity branches (L = 1, 2, and 3), and the OP of MRC outperforms that of the single-link. Furthermore, the asymptotic results at the high average SNR region and the numerical results by using the mixed α-μ model are also presented, respectively. We can see that the former keeps tight with the exact results at the high average SNR level while the latter matches well with the exact results in the whole average SNR region with the truncated term N = 12. These results demonstrate the validity of our derivations. In Figure 4, some i.n.i.d. and correlated examples are considered. As expected, we can get some same conclusions as those in Figure 3. However, for the correlated cases, the OP of MRC degrades with the increase in the correlation coefficients ( ρ = ρ N = ρ G = 0, 0.2, 0.5, and 0.7).
In Figure 5 and Figure 6, we illustrate the ABEP/ASEP of different modulation schemes for various communication systems as a function of the average SNR ( γ ¯ ) over α-μ/IGA fading channels. In Figure 5, the ABEP of DPSK for various communication systems are compared under i.i.d. fading scenarios. It can be seen that the ABEP is improved when α (α = 1.5 and 3), the number of diversity branches (L = 1, 2, and 3), and the average SNR ( γ ¯ = 0→20 dB) increases, respectively. Furthermore, the asymptotic results at the high average SNR region and the numerical results by using the mixed α-μ model also show similar behaviors as those in Figure 3. These results demonstrate the parameters α and L can provide the diversity gain, as shown in (64). In Figure 6, based on (42) and (43), we obtain the ASEP of MPSK for the MRC system in the i.n.i.d. and correlated cases. It is clear that the ASEP of MPSK degrades with the increase in the correlation coefficients ( ρ = ρ N = ρ G = 0 and 0.5) and/or the order of modulation (M = 2, 4, 8, and 16).
Figure 7 plots the average channel capacity per unit bandwidth of various communication systems versus the average SNR ( γ ¯ ) over several different communication conditions. It is evident that the average capacity improves with the growth of the value of α, μ, n, L, and/or γ ¯ , respectively. However, the increase in ρ   ( ρ = ρ N = ρ G ) leads to the reduction in the average capacity. Interestingly, the growing trends of average capacity become slight when the value of n gets larger (n = 5→50) and/or when the value of ρ gets smaller ( ρ = 0.5 0 ). These results suggest these parameters have less impact on the average capacity than other parameters (α, μ, L, and/or γ ¯ ). More importantly, the approximated analysis using (49) and the mixed α-μ model matches well with the exact numerical results. Moreover, Figure 8 shows the EC of various communication systems as a function of the average SNR ( γ ¯ ) over several different communication conditions with n = 5 and μ = 2. It can be seen from Figure 8 that the EC gets larger with the increase in the value of α, L, and/or γ ¯ , respectively. On the other hand, the EC gets smaller when the value of A and/or ρ   ( ρ = ρ N = ρ G ) grows, where the parameter A has a dominating impact on the EC. In addition, the approximated results by using (57) and the mixed α-μ model show good agreements with the exact numerical results.

7. Conclusions

In this paper, we presented a comprehensive investigation of the α-μ/IGA composite fading model and applied this model to the MRC system. First, we studied the fundamental statistics of the univariate and bivariate α-μ/IGA fading models including the PDF, CDF, and MGF, and proposed a mixture α-μ model to approximate the α-μ/IGA distribution. Then, as per the derived statistical expressions, the exact and approximated expressions of the statistical properties of the sum of α-μ/IGA variates by utilizing the multivariate Fox’s H-functions were obtained under the independent and the correlated fading environments. Third, some performance metrics of interest comprising the OP, ABEP/ ASEP, the average channel capacity, and the effective capacity were investigated in detail for the single-link and MRC systems, as well as their corresponding exact and approximated expressions, were derived in different channels conditions. Some asymptotic performance behaviors over α-μ/IGA fading were also provided in the high SNR regions. Finally, numerical analysis and Monte Carlo simulations were performed to prove the validity of the theoretical analysis under the different channel and system parameters. Our results will be helpful to improve the system reliability in the design and deployment of future communication applications, including 5G/6G wireless cellular systems, wearable communications, and vehicular wireless networks. Due to the simple and flexible form of the mixed α-μ model, we will apply it to other channel models and the complicated communication systems in future works, such as mmWave communications and Intelligent Reflecting Surface (IRS).

Author Contributions

Conceptualization, W.C. and G.W.; methodology, W.C. and T.M.; software, W.C. and T.M.; validation, W.C., T.M. and G.W.; investigation, T.M.; writing—original draft preparation, T.M. and W.C.; writing—review and editing, W.C. and G.W.; supervision, W.C.; project administration, W.C.; funding acquisition, G.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Open Project Fund of the Key Laboratory of Mine Disaster Prevention and Control under grant number SMDPC202102.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors appreciate and acknowledge the financial support from the Open Project Fund of the Key Laboratory of Mine Disaster Prevention and Control (Grant No. SMDPC202102).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of (20).
Since we assume the α-μ/IGA RVs are independent of each other, and the MGF of the lth RV can be given by using (9), then the MGF of γ M R C can be obtained as
M G F γ M R C ( s ) = l = 1 L M G F γ l ( s ) = l = 1 L A l H 2 , 1 1 , 2 [ λ l s 0.5 α l | ( μ l , 1 ) ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ] .
To get the desired results, with the aid of the definition of the univariate Fox H-function in [38], (A1) can be expanded as
M G F γ M R C ( s ) = Θ 1 ( 1 2 π j ) L ( 1 L l = 1 L ( Δ l ( λ l s 0.5 α l ) t l ) d t 1 d t L ) ,
where Θ 1 = l = 1 L α l / [ 2 Γ ( μ l ) Γ ( n l ) ] , Δ l =   Γ ( μ l t l ) Γ ( 0.5 α l t l ) Γ ( n l + 0.5 α l t l ) , j = 1 , and l denotes the lth appropriate contour which starts at the point   τ l j and goes to the point τ l j with τ l is a constant value, l { 1 , , L } .
Mathematically, the PDF of one RV can be found by applying the inverse Laplace transform of its MGF, then the PDF of γ M R C can be written as
f γ M R C ( γ ) = L 1 { M G F γ M R C ( s ) ; γ } = 1 2 π j M G F γ M R C ( s ) exp ( s γ ) d s ,
where 1 { · } represents the inverse Laplace transform. Thus, by substituting (A2) into (A3), the PDF of γ can be expressed as
f γ M R C ( γ ) = Θ 1 ( 1 2 π j ) L ( 1 L ( ϒ ) l = 1 L ( Δ l ( λ l ) t l ) d t 1 d t L ) ,
where ϒ = 1 2 π j s Δ 0 exp ( s γ ) d s , Δ 0 = l = 1 L 0.5 α l t l . By using [36] (Equation (8.315.1)), Υ can be solved as
ϒ = γ Δ 0 1 / Γ ( Δ 0 ) .
By plugging (A5) into (A4), we can re-express (A4) as
f γ M R C ( γ ) = Θ 1 γ 1 ( 1 2 π j ) L ( 1 L Γ 1 ( Δ 0 ) l = 1 L ( Δ l ( λ l γ 0.5 α l ) t l ) d t 1 d t L ) .
With the help of the definition of the multivariate H-function [37] (Equation (A.1)), the multiple Barnes-type contour integrals in (A6) can be written as the compact form in (20), which completes the proof. □

Appendix B

Proof of (42).
By substituting (A6) into (57) where Pe ( ε | γ ) is in the complementary error function form, the ABEP/ASEP of an L-branch MRC system over α-μ/IGA fading channels can be obtained as
P e ( ε ) = C Θ 1 ( 1 2 π j ) L ( C 1 C 2 C L ( ϒ 1 Γ ( Δ 0 ) ) l = 1 L ( Δ l λ l t l ) d t 1 d t 2 d t L )   D Θ 1 ( 1 2 π j ) L ( C 1 C 2 C L ( ϒ 2 Γ ( Δ 0 ) ) l = 1 L ( Δ l λ l t l ) d t 1 d t 2 d t L ) ,
where ϒ 1 = 0 γ Δ 0 1 e r f c ( E γ ) d γ , ϒ 2 = 0 γ Δ 0 1 e r f c 2 ( E γ ) d γ .
To calculate the integral terms in ϒ 1 and ϒ 2 , we apply a similar method as in (40) and the Mellin-Barnes type contour integral of the univariate Fox H-function, ϒ 1 and ϒ 2 can be derived, respectively, as
ϒ 1 = Γ ( Δ 0 ) Γ ( 0.5 + Δ 0 ) π E Δ 0 Γ ( 1 + Δ 0 ) ,
ϒ 2 = E Δ 0 π ( 2 π j ) L + 1 Γ ( Δ 0 + t L + 1 ) Γ ( 0.5 + Δ 0 + t L + 1 ) Γ ( t L + 1 ) Γ ( 0.5 t L + 1 ) Γ ( 1 t L + 1 ) Γ ( 1 + Δ 0 + t L + 1 ) 1 t L + 1 d t L + 1 .
Finally, by inserting ϒ 1 and ϒ 2 into (A7), (42) can be obtained by using the compact form of the multivariate Fox H-function, which completes the proof. □

References

  1. Simon, M.K.; Alouini, M.S. Digital Communication over Fading Channels, 2nd ed.; Wiley: Hoboken, NJ, USA, 2005. [Google Scholar]
  2. Moualeu, J.M.; da Costa, D.B.; Lopez-Martinez, F.J.; de Souza, R.A.A. On the performance of α-η-κ-μ fading channels. IEEE Commun. Lett. 2019, 23, 967–970. [Google Scholar] [CrossRef]
  3. Abdi, A.; Kaveh, M. K distribution: An appropriate substitute for Rayleigh-lognormal distribution in fading-shadowing wireless channels. Electron. Lett. 1998, 34, 851–852. [Google Scholar] [CrossRef] [Green Version]
  4. Eltoft, T. The Rician inverse Gaussian distribution: A new model for non-Rayleigh signal amplitude statistics. IEEE Trans. Image Process. 2005, 14, 1722–1735. [Google Scholar] [CrossRef] [PubMed]
  5. Bhatt, M.; Soni, S.K. A unified performance analysis of energy detector over α-η-μ/lognormal and α-κ-μ/lognormal composite fading channels with diversity and cooperative spectrum sensing. Int. J. Electron. Commun. (AEÜ) 2018, 94, 367–376. [Google Scholar] [CrossRef]
  6. Al-Hmood, H.; Al-Raweshidy, H.S. Unified modeling of composite κ-μ/gamma, η-μ/gamma, and α-μ/gamma fading channels using a mixture gamma distribution with applications to energy detection. IEEE Antennas Wirel. Propag. Lett. 2017, 16, 104–108. [Google Scholar] [CrossRef]
  7. Chauhan, P.S.; Soni, S.K. Average SEP and channel capacity analysis over Generic/IG composite fading channels: A unified approach. Phys. Commun. 2019, 34, 9–18. [Google Scholar] [CrossRef]
  8. Paris, J.F. Statistical characterization of κ-μ shadowed fading. IEEE Trans. Veh. Technol. 2014, 63, 518–526. [Google Scholar] [CrossRef] [Green Version]
  9. Simmons, N.; da Silva, C.R.N.; Cotton, S.L.; Sofotasios, P.C.; Yacoub, M.D. Double shadowing the Rician fading model. IEEE Wirel. Commun. Lett. 2019, 8, 344–347. [Google Scholar] [CrossRef]
  10. Simmons, N.; Da Silva, C.R.N.; Cotton, S.L.; Sofotasios, P.C.; Yoo, S.K.; Yacoub, M.D. On shadowing the κ-μ fading model. IEEE Access 2020, 8, 120513–120536. [Google Scholar] [CrossRef]
  11. Leonardo, E.J.; Yacoub, M.D. The product of two α-μ variates and the composite α-μ multipathshadowing model. IEEE Trans. Veh. Technol 2015, 64, 2720–2725. [Google Scholar] [CrossRef]
  12. Badarneh, O.S. The α-μ/α-μ composite multipath-shadowing distribution and its connection with the extended generalized-K distribution. AEU-Int. J. Elect. Commun. 2016, 70, 1211–1218. [Google Scholar] [CrossRef]
  13. Yoo, S.; Cotton, S.; Zhang, L.; Sofotasios, P. The inverse Gamma distribution: A new shadowing model. In Proceedings of the 8th Asia-Pacific Conference on Antennas and Propagation, Incheon, Korea, 4–7 August 2019; pp. 475–476. [Google Scholar] [CrossRef]
  14. Ramırez-Espinosa, P.; Lopez-Martınez, F.J. On the utility of the inverse gamma distribution in modeling composite fading channels. In Proceedings of the IEEE Global Communications Conference, Waikoloa, HI, USA, 9–13 December 2019; pp. 1–6. [Google Scholar] [CrossRef] [Green Version]
  15. Ramírez-Espinosa, P.; López-Martínez, F.J. Composite fading models based on inverse Gamma shadowing: Theory and validation. IEEE Trans. Wirel. Commun. 2021, 20, 5034–5045. [Google Scholar] [CrossRef]
  16. Bithas, P.; Nikolaidis, V.; Kanatas, A.G.; Karagiannidis, G.K. UAV-to-ground communications: Channel modeling and UAV selection. IEEE Trans. Commun. 2020, 68, 5135–5144. [Google Scholar] [CrossRef]
  17. Yoo, S.K.; Cotton, S.L.; Sofotasios, P.C.; Matthaiou, M.; Valkama, M.; Karagiannidis, G.K. The κ-μ/inverse gamma fading model. In Proceedings of the IEEE PIMRC, Hong Kong, China, 30 August–2 September 2015; pp. 425–429. [Google Scholar]
  18. Yoo, S.K.; Sofotasios, P.C.; Cotton, S.L.; Matthaiou, M.; Valkama, M.; Karagiannidis, G.K. The η-μ/inverse gamma composite fading model. In Proceedings of the IEEE PIMRC, Hong Kong, China, 30 August–2 September 2015; pp. 166–170. [Google Scholar]
  19. Yoo, S.K.; Bhargav, N.; Cotton, S.L.; Sofotasios, P.C.; Matthaiou, M.; Valkama, M.; Karagiannidis, G.K. The κ-μ/inverse gamma and η-μ/inverse gamma composite fading models: Fundamental statistics and empirical validation. IEEE Trans. Commun. 2017, 69, 5514–5530. [Google Scholar] [CrossRef]
  20. Yoo, S.K.; Cotton, S.L.; Sofotasios, P.C.; Muhaidat, S.; Karagiannidis, G.K. Effective capacity analysis over generalized composite fading channels. IEEE Access 2020, 8, 123756–123764. [Google Scholar] [CrossRef]
  21. Pant, D.; Chauhan, P.S.; Soni, S.K.; Naithani, S. BER and channel capacity analysis of wireless system over κ-μ/inverse gamma and η-μ/inverse gamma composite fading model. Wirel. Netw. 2021, 27, 1251–1267. [Google Scholar] [CrossRef]
  22. Pant, D.; Chauhan, P.S.; Soni, S.K. Error probability and channel capacity analysis of wireless system over inverse gamma shadowed fading channel with selection diversity. Int. J. Commun. Syst. 2019, 32, e4083. [Google Scholar] [CrossRef]
  23. Upaddhyay, V.K.; Chauhan, P.S.; Soni, S.K. Effective capacity analysis of SIMO system with MRC and SC over Inverse-Gamma shadowing. Int. J. Electron. 2021, 109, 181–199. [Google Scholar] [CrossRef]
  24. Cheng, W.J.; Wang, X.T.; Xu, X.M. On the performance analysis of wireless transmission over α-μ/inverse gamma composite fading channels. In Proceedings of the IEEE/CIC International Conference on Communications in China (ICCC Workshops), Chongqing, China, 9–11 August 2020; pp. 64–69. [Google Scholar] [CrossRef]
  25. Maurya, R.; Chauhan, P.S.; Srivastava, S.; Soni, S.K.; Mishra, B. Energy detection investigation over composite α-μ/inverse-Gamma wireless channel. Int. J. Electron. Commun. 2020, 130, 153556. [Google Scholar] [CrossRef]
  26. Upaddhyay, V.K.; Soni, S.K.; Chauhan, P.S. An approximate statistical analysis of wireless channel over α-μ shadowed fading channel. Int. J. Commun. Syst. 2021, 34, e4884. [Google Scholar] [CrossRef]
  27. Goswami, A.; Kumar, A. Statistical Characterization and Performance Evaluation of α−η−μ/Inverse Gamma and α−κ−μ/Inverse Gamma Channels. Wirel. Pers. Commun. 2022, 124, 2313–2333. [Google Scholar] [CrossRef]
  28. Yoo, S.K.; Cotton, S.L.; Sofotasios, P.C.; Matthaiou, M.; Valkama, M.; Karagiannidis, G.K. The Fisher-Snedecor distribution: A simple and accurate composite fading model. IEEE Commun. Lett. 2017, 21, 1661–1664. [Google Scholar] [CrossRef] [Green Version]
  29. Cheng, W.J.; Wang, X.T. Bivariate Fisher–Snedecor distribution and its applications in wireless communication systems. IEEE Access 2020, 8, 146342–146360. [Google Scholar] [CrossRef]
  30. Cheng, W.J.; Wang, X.T.; Ma, T.F.; Wang, G. On the performance analysis of switched diversity combining receivers over Fisher–Snedecor composite fading channels. Sensors 2021, 21, 3014. [Google Scholar] [CrossRef]
  31. Zhang, J.Y.; Du, H.Y.; Sun, Q.; Ai, B.; Ng, D.W.K. Physical layer security enhancement with reconfigurable intelligent surface-aided networks. IEEE Trans. Inf. Forensics Secur. 2021, 16, 3480–3495. [Google Scholar] [CrossRef]
  32. Badarneh, O.S. The α- composite fading distribution: Statistical characterization and applications. IEEE Trans. Veh. Technol. 2020, 69, 8097–8106. [Google Scholar] [CrossRef]
  33. Badarneh, O.S. The α-η- and α-κ- composite fading distributions. IEEE Commun. Lett. 2020, 24, 1924–1928. [Google Scholar] [CrossRef]
  34. Badarneh, O.S.; Muhaidat, S.; Costa, D.B.D. The α-η-κ- composite fading distribution. IEEE Wirel. Commun. Lett. 2020, 9, 2182–2186. [Google Scholar] [CrossRef]
  35. Treanta, S. Optimization on the distribution of population densities and the arrangement of urban activities. Stat. Optim. Inf. Comput. 2018, 6, 208–218. [Google Scholar] [CrossRef]
  36. Yacoub, M.D. The α-μ distribution: A physical fading model for the stacy distribution. IEEE Trans. Veh. Technol. 2007, 56, 27–34. [Google Scholar] [CrossRef]
  37. Gradshteyn, I.; Ryzhik, I. Table of Integrals, Series, and Products, 7th ed.; Academic: New York, NY, USA, 2007. [Google Scholar]
  38. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theory and Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  39. Aldalgamouni, T.; Ilter, M.C.; Badarneh, O.S.; Yanikomeroglu, H. Performance analysis of fisher-snedecor composite fading channels. In Proceedings of the IEEE Middle East and North Africa Communications Conference, Jounieh, Lebanon, 18–20 April 2018; pp. 1–5. [Google Scholar] [CrossRef]
  40. Al-Hmood, H.; Al-Raweshidy, H.S. Selection combining scheme over non-identically distributed Fisher-Snedecor fading channels. IEEE Wirel. Commun. Lett. 2021, 10, 840–843. [Google Scholar] [CrossRef]
  41. Badarneh, O.S.; Costa, D.B.D.; Sofotasios, P.C.; Muhaidat, S.; Cotton, S.L. On the sum of Fisher-Snedecor variates and its application to maximal-ratio combining. IEEE Wirel. Commun. Lett. 2018, 7, 966–969. [Google Scholar] [CrossRef] [Green Version]
  42. Atapattu, S.; Tellambura, C.; Jiang, H. A mixture gamma distribution to model the SNR of wireless channels. IEEE Trans. Wirel. Commun. 2011, 10, 4193–4203. [Google Scholar] [CrossRef]
  43. Ramirez-Espinosa, P.; Moualeu, J.M.; da Costa, D.B.; Lopez-Martinez, F.J. The α-κ-μ shadowed fading distribution: Statistical characterization and applications. In Proceedings of the IEEE Global Communications Conference, Waikoloa, HI, USA, 9–13 December 2019; pp. 1–6. [Google Scholar] [CrossRef] [Green Version]
  44. Bithas, P.S.; Kanatas, A.G.; Matolak, D.W. Exploiting shadowing stationarity for antenna selection in V2V communications. IEEE Trans. Veh. Technol. 2018, 68, 1607–1615. [Google Scholar] [CrossRef]
  45. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series (More Special Functions); Gordon and Breach: Philadelphia, PA, USA, 1992; Volume 3. [Google Scholar]
  46. Kilbas, A.A.; Saigo, M. H-Transforms: Theory and Applications (Analytical Method and Special Function), 1st ed.; Chapman & Hall/CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar] [CrossRef]
  47. Mittal, P.K.; Gupta, K.C. An integral involving generalized function of two variables. Proc. Indian Acad. Sci.-Sect. A 1972, 75, 117–123. [Google Scholar] [CrossRef]
  48. Rahama, Y.A.; Ismail, M.H.; Hassan, M. On the sum of independent Fox’s H-Function variates with applications. IEEE Trans. Veh. Technol. 2018, 67, 6752–6760. [Google Scholar] [CrossRef]
  49. Du, H.Y.; Zhang, J.Y.; Cheng, J.L.; Ai, B. Sum of Fisher-Snedecor random variables and its applications. IEEE Open J. Commun. Soc. 2020, 1, 342–356. [Google Scholar] [CrossRef]
  50. Hamdi, K.A. Capacity of MRC on correlated rician fading channels. IEEE Trans. Commun. 2008, 56, 708–711. [Google Scholar] [CrossRef] [Green Version]
  51. Ji, Z.; Dong, C.; Wang, Y.; Lu, J. On the analysis of effective capacity over generalized fading channels. In Proceedings of the IEEE International Conference on Communications, Sydney, Australia, 10–14 June 2014; pp. 1977–1983. [Google Scholar] [CrossRef]
  52. Alhennawi, H.R.; El Ayadi, M.M.; Ismail, M.H.; Mourad, H.A.M. Closed-form exact and asymptotic expressions for the symbol error rate and capacity of the H-function fading channel. IEEE Trans. Veh. Technol. 2016, 65, 1957–1974. [Google Scholar] [CrossRef]
Figure 1. PDFs of the α-μ/IGA and α- F models for different channel parameters with γ ¯ = 1 .
Figure 1. PDFs of the α-μ/IGA and α- F models for different channel parameters with γ ¯ = 1 .
Applsci 12 05375 g001
Figure 2. Impacts of different system and channel parameters on the exact PDF in (7), the mixed α-μ PDF in (10), and the mixed IGA PDF in (13).
Figure 2. Impacts of different system and channel parameters on the exact PDF in (7), the mixed α-μ PDF in (10), and the mixed IGA PDF in (13).
Applsci 12 05375 g002
Figure 3. Outage probability of various communication systems as a function of the average SNR over i.i.d. α-μ/IGA fading channels with n = 5, μ = 2, and γ t h = 5 dB.
Figure 3. Outage probability of various communication systems as a function of the average SNR over i.i.d. α-μ/IGA fading channels with n = 5, μ = 2, and γ t h = 5 dB.
Applsci 12 05375 g003
Figure 4. Outage probability of various communication systems as a function of the average SNR ( γ ¯ ) over i.n.i.d. and correlated α-μ/IGA fading channels with n = 5, μ = 2, and γ t h = 5 dB.
Figure 4. Outage probability of various communication systems as a function of the average SNR ( γ ¯ ) over i.n.i.d. and correlated α-μ/IGA fading channels with n = 5, μ = 2, and γ t h = 5 dB.
Applsci 12 05375 g004
Figure 5. ABEP of DPSK for various communication systems as a function of the average SNR ( γ ¯ ) over i.i.d. α-μ/IGA fading channels with n = 5 and μ = 2.
Figure 5. ABEP of DPSK for various communication systems as a function of the average SNR ( γ ¯ ) over i.i.d. α-μ/IGA fading channels with n = 5 and μ = 2.
Applsci 12 05375 g005
Figure 6. ASEP of MPSK for the MRC system as a function of the average SNR ( γ ¯ ) over i.n.i.d. and correlated α-μ/IGA fading channels with n = 5 and μ = 2.
Figure 6. ASEP of MPSK for the MRC system as a function of the average SNR ( γ ¯ ) over i.n.i.d. and correlated α-μ/IGA fading channels with n = 5 and μ = 2.
Applsci 12 05375 g006
Figure 7. Average channel capacity per unit bandwidth of various communication systems as a function of the average SNR ( γ ¯ ) over different α-μ/IGA fading channels.
Figure 7. Average channel capacity per unit bandwidth of various communication systems as a function of the average SNR ( γ ¯ ) over different α-μ/IGA fading channels.
Applsci 12 05375 g007
Figure 8. Effective capacity of various communication systems as a function of the average SNR ( γ ¯ ) over different communication conditions with n = 5 and μ = 2.
Figure 8. Effective capacity of various communication systems as a function of the average SNR ( γ ¯ ) over different communication conditions with n = 5 and μ = 2.
Applsci 12 05375 g008
Table 1. Comparisons of the MSE and the KLD for the exact PDF and the mixed PDFs with different channel parameters (α, μ, and n) and the truncated term (N).
Table 1. Comparisons of the MSE and the KLD for the exact PDF and the mixed PDFs with different channel parameters (α, μ, and n) and the truncated term (N).
Channel
Parameters
N = 8N = 20
Mixed IGAMixed α-μMixed IGAMixed α-μ
MSEKLDMSEKLDMSEKLDMSEKLD
α = 2, μ = 2n = 20.0011580.0072851.148 × 10−50.0093369.947 × 10−50.0010563.008 × 10−86.149 × 10−4
n = 200.039360.40762.766 × 10−50.0095830.0084850.091452.933 × 10−191.816 × 10−9
α = 2, n = 5μ = 10.051021.53894.759 × 10−112.399 × 10−50.024870.52769.694 × 10−208.243 × 10−10
μ = 51.491 × 10−52.3875 × 10−41.258 × 10−64.8473 × 10−41.886 × 10−83.608 × 10−61.047 × 10−101.048 × 10−5
α = 2, n = 5α = 10.14490.31542.345 × 10−125.906 × 10−60.035690.080062.536 × 10−172.433 × 10−8
α = 52.526 × 10−58.651 × 10−42.366 × 10−40.0076681.012 × 10−61.378 × 10−48.094 × 10−74.4059 × 10−4
Table 2. Theoretical conditional BEP/SEP for several classical modulation schemes.
Table 2. Theoretical conditional BEP/SEP for several classical modulation schemes.
BEP / SEP   ( P e ( ε | γ ) ) Modulation SchemeParameter Values
A e x p ( B γ ) DPSK A = 0.5 ,   B = 1
C e r f c ( E γ ) D e r f c 2 ( E γ ) Coherent BPSK C = 0.5 ,   D = 0 ,   E = 1
MPSK C = 1 ,   D = 0 ,   E = sin 2 ( π / M ) , M > 4
Square MQAM C = 2 ( 1 1 / M ) ,   D = ( 1 1 / M ) 2 , E = 1.5 / ( M 1 )
Table 3. ABEP expressions of DPSK for various systems.
Table 3. ABEP expressions of DPSK for various systems.
Communication Systems ABEP (Pe)
Single-link system by using (9) A α 2 Γ ( μ ) Γ ( n ) H 2 , 1 1 , 2 [ λ B 0.5 α | ( μ , 1 ) ( 1 n , 0.5 α ) , ( 1 , 0.5 α ) ]
Single-link system by using (12) i = 1 N A ς i α 2 Γ ( μ ) H 1 , 1 1 , 1 [ ( B ψ i ) 0.5 α | ( μ , 1 ) ( 1 , 0.5 α ) ]
L-branch MRC system
by using (A1)
A l = 1 L A l H 2 , 1 1 , 2 [ λ l B 0.5 α l | ( μ l , 1 ) ( 1 n l , 0.5 α l ) , ( 1 , 0.5 α l ) ]
L-branch MRC system
by using (23)
i 1 = 1 N 1 i L = 1 N L A l = 1 L ( ς i l α l 2 Γ ( μ l ) H 1 , 1 1 , 1 [ ( B ψ i l ) 0.5 α l | ( μ l , 1 ) ( 1 , 0.5 α l ) ] )
Two-branch MRC system by using (32) k = 0 l = 0 0.25 A Φ ρ N k ρ G l k ! l ! Γ ( k + μ ) Γ ( l + n ) i = 1 2 H 2 , 1 1 , 2 [ ϑ i B 0.5 α i | ( μ + k , 1 ) ( 1 n l , 0.5 α i ) , ( 1 , 0.5 α i ) ]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, W.; Ma, T.; Wang, G. On the Sum of α-μ/Inverse Gamma Variates with Applications to Diversity Receivers. Appl. Sci. 2022, 12, 5375. https://doi.org/10.3390/app12115375

AMA Style

Cheng W, Ma T, Wang G. On the Sum of α-μ/Inverse Gamma Variates with Applications to Diversity Receivers. Applied Sciences. 2022; 12(11):5375. https://doi.org/10.3390/app12115375

Chicago/Turabian Style

Cheng, Weijun, Tengfei Ma, and Gang Wang. 2022. "On the Sum of α-μ/Inverse Gamma Variates with Applications to Diversity Receivers" Applied Sciences 12, no. 11: 5375. https://doi.org/10.3390/app12115375

APA Style

Cheng, W., Ma, T., & Wang, G. (2022). On the Sum of α-μ/Inverse Gamma Variates with Applications to Diversity Receivers. Applied Sciences, 12(11), 5375. https://doi.org/10.3390/app12115375

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