The distributions of random variables are of interest in many areas of science. In this paper, the probability density function (PDF) and cumulative distribution function (CDF) of ratio of products of two random variables and random variable are derived. Random variables are described with Rayleigh, Nakagami-m, Weibull, and α-μ distributions. An application of obtained results in performance analysis of multihop wireless communication systems in different transmission environments described in detail. The proposed mathematical analysis is also complemented by various graphically presented numerical results.
1. Introduction
The distribution of the ratio of random variables is of interest in statistical analysis in biological and physical sciences, econometrics, and ranking and selection [1]. It has been studied by several authors [2–5]. The ratio of random variables is also of interest in analysing wireless communication systems [6–9]. In wireless communication systems, one of problems is fading caused by multipath propagation. When a received signal experiences fading during transmission, signal envelope fluctuates over time [10]. There is a very wide range of statistical models for describing the statistical behaviour of signal envelope where accuracy and veracity depend on communication scenario and propagation environment. The most frequently applied models in the open technical literature are Rayleigh, Nakagami-m, Weibull, and α-μ. In addition to fading, the wireless transmission can be also subjected to shadowing which is the result of large obstacles and deviations in terrain profile between transmitter and receiver. In such composite environments, signal envelope can be modelled with product of two random variables [11, 12]. A fundamental requirement in planning wireless communication systems, besides providing specified quality of service (QoS), is to ensure high system capacity. Cochannel interference (CCI) is a result of frequency reuse which is essential in increasing system capacity [13].
In the areas of deployment, connectivity, and capacity, multihop systems have a number of advantages over traditional communication networks. Relaying techniques enable network connectivity where traditional architectures are impractical due to location constraints and can be applied to cellular, wireless local area networks (WLANs) and hybrid networks. The performance analysis of multihop wireless communication systems operating in fading channels has been an important field of research in the past few years [14–17].
The main contribution of this paper is the derivation of useful closed-form expressions for the probability density function (PDF) and cumulative distribution function (CDF) of ratio of product of two random variables and random variable. Among the rest, these expressions can be used as general theoretical tool in the performance analysis of wireless communication systems since the product of random variables in nominator can represent signal envelope which suffers from fading and shadowing while the random variable in denominator can represent CCI envelope. Therefore, the ratio of such random variables is signal-to-interference ratio (SIR). Rayleigh, Nakagami-m, Weibull, and α-μ distributions are included in our analysis to enable wide range of use of results presented in the paper, that is, to model different transmission environments. Capitalizing on closed-form expressions derived in the paper in the terms of Meijer G functions, the outage probability of multihop system, as an important and widely accepted performance measure, is determined and discussed.
2. Statistics of Ratio of Products of Two Random Variables and Random Variable
In this section, statistical analysis of ratio of products of two random variables and random variable, λ=xz/y, is presented. For derivation of the corresponding PDF of λ, we have to obtain the PDF of products of two random variables, t=xz, using the following equation:pt(t)=∫0∞|J|px(tz)pz(z)dz,
where |J| is the Jacobian transformation given by |J|=|dx/dt|=1/z. Using similar procedure, the PDF of λ can be derived aspλ(λ)=∫0∞|J|pt(λy)py(y)dy,
where |J|=|dt/dλ|=y.
Cumulative distribution function can be obtained by definition asFλ(λ)=∫0λpλ(λ)dλ.
First, we present different distributions describing statistical behaviour of variables x,y, and z.
(a) Rayleigh Distribution
In the case of Rayleigh distribution, the PDFs of random variables are given by
px(x)=2Ωxxe-x2/Ωx,py(y)=2Ωyye-y2/Ωy,pz(z)=2Ωzze-z2/Ωz,
where Ωx=ε〈x2〉,Ωy=ε〈y2〉,Ωz=ε〈z2〉, and ε〈·〉 denotes expectation.
(b) Weibull Distribution
The Weibull random variables are distributed according to
px(x)=βΩxxβ-1e-xβ/Ωx,py(y)=βΩyyβ-1e-yβ/Ωy,pz(z)=βΩzzβ-1e-zβ/Ωz,
whereΩx=ε〈xβ〉,Ωy=ε〈yβ〉,Ωz=ε〈zβ〉, and β is Weibull parameter (for the special case of β=2 Weibull PDFs reduces to Rayleigh PDFs).
(c) Nakagami-m Distribution
The Nakagami-m PDFs are
px(x)=2mxmxΓ(mx)Ωxmxx2mx-1e-mxx2/Ωx,py(y)=2mymyΓ(my)Ωymyy2my-1e-myy2/Ωy,pz(z)=2mzmzΓ(mz)Ωzmzz2mz-1e-mzz2/Ωz,
where Ωx=ε〈x2〉,Ωy=ε〈y2〉,Ωz=ε〈z2〉, Γ(·) is gamma function, and mx,my, and mz are Nakagami-m parameters that range from 0.5 to ∞. Nakagami-m distribution spans via m parameter through wide range of distributions which includes the one-sided Gaussian distribution (m=0.5) and the Rayleigh distribution (m=1) as special cases.
(d) α-μ Distribution
The α-μ distributed random variables are described with following equations:
px(x)=α(μxΩx)μxxαμx-1Γ(μx)e-μxxα/Ωx,py(y)=α(μyΩy)μyyαμy-1Γ(μy)e-μyyα/Ωy,pz(z)=α(μzΩz)μzzαμz-1Γ(μz)e-μzzα/Ωz,
where Ωx=ε〈xα〉,Ωy=ε〈yα〉,Ωz=ε〈zα〉, α is parameter related to the nonlinearity, and μx,μy, and μz are the inverse of the normalized variance of xα,yα and zα, respectively (μx,μy,μz≥0.5). The α-μ distribution is a general distribution that includes as special cases Nakagami-m distribution for α=2 and Weibull distribution for μ=1.
Applying the procedure described at the beginning of this section with the aid of [18, Equations (3.461) and (6.631)], [19], [20, Equation (26)], the PDFs and CDFs of λ in different scenarios can be expressed in terms of Meijer G functions Gp,qm,n(γ|a1,a2,…,apb1,b2,…,bq) as
The PDFs and CDFs of λ in all considered scenarios are presented in Figures 1, 2, 3 and 4.
PDF and CDF of λ in Rayleigh scenario.
PDF and CDF of λ in Weibull scenario.
PDF and CDF of λ in Nakagami-m scenario.
PDF and CDF of λ in α-μ scenario.
3. An Application of Presented Theoretical Results in Wireless Communication Systems: The Performance Analysis of Multihop Systems
Multihop (N-hop) transmission is a technique by which the channel from the source (S) to the destination (D) is split into several, possibly shorter, links using relays (Ri,i=1,N-1̅), as shown in Figure 5.
Multihop transmission.
Depending on the nature of the radio propagation environment, the statistical behaviour of signal envelope can be described with different distributions. This is reason why our analysis includes so many statistical models. The Rayleigh model is frequently used to describe multipath fading with no direct line-of-sight (LOS) path [10]. The Weibull distribution exhibits an excellent fit to experimental fading channel measurements, for both indoor [21] and outdoor [22] environments. Also good results are provided in urban environments. The Nakagami-m distribution has gained widespread application in the modelling of physical radio channels since it shows great agreement with experimentally obtained results [23]. Fading severity in Weibull and Nakagami-m environments is described by Weibull parameter β and Nakagami parameter m, respectively. As these parameters increase fading severity decreases. The α-μ distribution is a general fading distribution that can be used to better represent the small-scale variation of the fading signal [24]. Parameter α is related to the nonlinearity of the environment, while parameter μ is associated to the number of multipath [25].
The product of random variables, ti=xizi, represents signal envelope which suffers from fading and shadowing while the random variable yi represents CCI envelope at the input of ith terminal (i=1,N̅). The random variable λi=xizi/yi presents SIR value at the input of ith terminal. Therefore, the mathematical results presented in the previous section can be efficiently used for performance analysis of relayed communication systems in fading environment. Highly important and widely accepted system performance indicator is outage probability defined as probability of having SIR value lower than predetermined threshold λ0 which defines required QoS. The system failure can occur in sections S-R1,R1-R2,R2-R3,…,Rj-1-Rj,…,RN-1-D when some of the values of λ1,λ2,λ3,…,λj,…,λN are below the predetermined threshold λ0. The PDF of minimum of λi, λ=min{λ1,λ2,…,λN}, is important for analyzing multihop relayed communication systems in which the source terminal communicates with the destination terminal through a number of relay terminals. It can be obtained based on previous results aspλ(λ)=∑n=1Npλn(λ)∏k=1k≠nN(1-Fλk(λ)).
The outage probability of multihop system is defined asPout=∫0λ0pλ(λ)dλ.
As an illustrative example, the outage probability of dual-hop and triple-hop communication systems can be obtained asPout=Fλ1(λ)(1-Fλ2(λ))+Fλ2(λ)(1-Fλ1(λ))+Fλ1(λ)Fλ2(λ),Pout=Fλ1(λ)(1-Fλ2(λ))(1-Fλ3(λ))+Fλ2(λ)(1-Fλ1(λ))(1-Fλ3(λ))+Fλ3(λ)(1-Fλ1(λ))(1-Fλ2(λ))+Fλ1(λ)Fλ2(λ)(1-Fλ3(λ))+Fλ1(λ)Fλ3(λ)(1-Fλ2(λ))+Fλ2(λ)Fλ3(λ)(1-Fλ1(λ))+Fλ1(λ)Fλ2(λ)Fλ3(λ),
respectively. The curves representing PDF of λ and outage probability of multihop (dual-hop and triple-hop) communication system in Weibull fading environment are presented in Figure 6. Without loss of generality, we assumed that the ratios of average powers are equal on all terminals inputs, that is, γi=ΩxiΩzi/Ωyi=γ,i=1,N¯. Namely, signal is amplified in terminal so that the ratio of average powers at the input of the next terminal is equal to the ratio of average powers at the input of previous one. It is noticeable that the outage probability is higher for lower values of Weibull parameter, that is, in environment with higher fading severity. Also, having in mind that γi=γ, i=1,N¯, the higher values of N imply larger distance between the source and destination terminals and higher values of the outage probability.
PDF of λ = min {λ1,λ2,…,λN} and outage probability of dual-hop (N=2) and triple-hop (N=3) communication system for different values of γ in Weibull fading environment.
4. Conclusion
The PDF and CDF of ratio of product of two random variables and random variable λ=xz/y=t/y have been derived. Rayleigh, Weibull, Nakagami-m, and α-μ statistical models are included in the paper, so that other researchers and engineers could use our results in wide range of scenarios in many areas of science. An application of these results for the wireless communications community has been described. Namely, presented results can help the designers of wireless communication systems to simulate different wireless environments where fading and shadowing affect desired signal and fading affects CCI and readjust system parameters in order to meet the QoS demands. In our future work, we hope to obtain analytical expressions for PDF and CDF of ratio of products of random variables which will have application in analysis of wireless communication systems where both desired signal and CCI are affected by fading and shadowing simultaneously which is realistic scenario in modern urban areas.
Acknowledgment
This work was supported by the Ministry of Science of Serbia within the Project TR-32023.
NadarajahS.KotzS.On the ratio of Pearson type VII and Bessel random variables2005200541911992-s2.0-3454759026010.1155/JAMDS.2005.191BasuA.P.LochnerR. H.On the distribution of the ratio of two random variables having generalized life distributions197113281287Pham-GiaT.Distributions of the ratios of independent beta variables and applications20002912269327152-s2.0-0012457557KorhonenP. J.NarulaS. C.The probability distribution of the ratio of the absolute values of two normal variables198933173182MarsagliaG.Ratios of normal variables and ratios of sums of uniform variables196560193204OkuiS.Effects of CIR selection diversity with two correlated branches in the m-fading channel20004810163116332-s2.0-003429107610.1109/26.871387StefanovičM. C.DraćaD. L.PanajotovičA. S.SekulovičN. M.sekulani@bankerinter.netPerformance analysis of system with L-branch selection combining over correlated Weibull fading channels in the presence of cochannel interference201023213915010.1002/dac.1050KaragiannidisG. K.Performance analysis of SIR-based dual selection diversity over correlated Nakagami-m fading channels2003525120712162-s2.0-014195398410.1109/TVT.2003.816612PanajotovićA. S.StefanovićM. Č.DračaD. L.Performance analysis of system with selection combining over correlated Rician fading channels in the presence of cochannel interference20096312106110662-s2.0-7174908844610.1016/j.aeue.2008.08.001SimonM. K.AlouiniM.-S.20001stNew York, NY, USAWileyCorazzaG. E.VatalaroF.Statistical model for land mobile satellite channels and its application to nongeostationary orbit systems19944337387422-s2.0-002849403910.1109/25.312773BabalisP. G.CapsalisC. N.Impact of the combined slow and fast fading channel characteristics on the symbol error probability for multipath dispersionless channel characterized by a small number of dominant paths19994756536572-s2.0-003267244510.1109/26.768753ParsonsJ. D.20002ndNew York, NY, USAWileyKaragiannidisG. K.TsiftsisT. A.MallikR. K.Bounds for multihop relayed communications in Nakagami-m fading200654118222-s2.0-3134447595310.1109/TCOMM.2005.861679HasnaM. O.AlouiniM. S.End-to-end performance of transmission systems with relays over Rayleigh-fading channels200326112611312-s2.0-004261045110.1109/TWC.2003.819030HasnaM. O.AlouiniM. S.Outage probability of multihop transmission over Nakagami fading channels2003752162182-s2.0-003760450610.1109/LCOMM.2003.812178ZhongC.JinS.WongK. K.Outage probability of dual-hop relay channels in the presence of interference69th IEEE Vehicular Technology Conference (VTC '09)April 2009152-s2.0-7034968059610.1109/VETECS.2009.5073759GradshteynI. S.RyzhikI. M.19945thNew York, NY, USAAcademichttp://functions.wolfram.com/07.45.26.0005.01AdamchikV. S.MarichevO. I.Algorithm for calculating integrals of hypergeometric type functions and its realization in reduce systemProceedings of International Symposium on Symbolic and Algebraic Computation (ISSAC '90)August 1990212224BabichF.LombardiG.Statistical analysis and characterization of the indoor propagation channel20004834554642-s2.0-003389819610.1109/26.837048TzeremesG.ChristodoulouC. G.Use of Weibull distribution for describing outdoor multipath fadingIEEE Antennas and Propagation Society International SymposiumJune 20022322352-s2.0-0036078715NakagamiM.The m-distribution—a general formula of intensity distribution of rapid fading1960Oxford, UKPergamon Press336YacoubM. D.michel@decom.fee.unicamp.brThe α-μ distribution: a general fading distribution213th IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC '02)September 200262963310.1109/PIMRC.2002.1047298Da CostaD. B.YacoubM. D.FraidenraichG.Second-order statistics of equal-gain and maximal-ratio combining for the α-μ (generalized gamma) fading distribution9th IEEE International Symposium on Spread Spectrum Symposium on Spread Spectrum Techniques and Applications (ISSSTA '06)August 20063423462-s2.0-4614911970110.1109/ISSSTA.2006.311791