Research Article  Open Access
Efficient Numerical Methods for Analysis of Square Ratio of κμ and ημ Random Processes with Their Applications in Telecommunications
Abstract
We will provide statistical analysis of the square ratio of and random processes and its application in the signaltointerference ratio (SIR) based performance analysis of wireless transmission subjected to the influence of multipath fading, modelled by fading model, and undesired occurrence of cochannel interference (CCI), distributed as random process. First contribution of the paper is deriving exact closed expressions for the probability density function (PDF) and cumulative distribution function (CDF) of square ratio of and random processes. Further, a verification of accuracy of these PDF and CDF expressions was given by comparison with the corresponding approximations obtained by the highprecision quadrature formulas of Gaussian type with respect to the weight functions on . The computational procedure of such quadrature rules is provided by using the constructive theory of orthogonal polynomials and the MATHEMATICA package OrthogonalPolynomials created by Cvetković and Milovanović (2004). Capitalizing on obtained expression, important wireless performance criteria, namely, outage probability (OP), have been obtained, as functions of transmission parameters. Also, possible performance improvement is observed through a glance at SC (selection combining) reception employment based on obtained expressions.
1. Introduction
Let and be mutually independent Gaussian random processes, with variances . Modelling mean values of Gaussians, with and , and assuming initially to be an integer, the resulting random process is obtained as the following function of Gaussians: which follows distribution, whose PDF is given in the following form (cf. [1]): with and is the zerothe order modified Bessel function of the first kind [2, Eq. 8.445].
As a general probability distribution model, this model includes some classical distribution models as its particular cases, for example, Rician, and Nakagamim distribution models as special cases (as the OneSided Gaussian and the Rayleigh distributions since they also represent special cases of Nakagamim).
Let and , , be Gaussian random processes, with . Assuming initially to be an integer, resulting random process obtained as a function of Gaussians follows distribution, whose PDF is given in the following form (cf. [1]):with defined as , the ratio of arbitrary variances of independent Gaussians defined as , and distribution parameters and defined as , , . As a general probability distribution model, this model also includes some classical distribution models as its particular cases, for example, Nakagami (Hoyt), onesided Gaussian, Rayleigh, and Nakagami.
Multipath fading is physical phenomena that occurs as randomly delayed components of desired signal combine in constructive or destructive manner at the reception [1]. In the case when a line of sight (LOS) component between the transmitter and the receiver is present, the most general fading model which describes the shortterm signal variation is kappamu () fading model [3, 4]. Fading model, defined as the function of parameter , related to the dominant/scattered components powers quotient, and parameter , related to the number of propagating clusters, easily reduces to other fading models, by setting corresponding values for parameters and . When , the model transforms to the Rice model with arbitrary factor, [4]. By assigning , the observed model transforms to the Nakagami model [4]. Cochannel interference (CCI) signal which is transmitted at the same frequency as the carrier signal is also exposed to multipath fading phenomenon [5]. Due to complicated propagation of interfering signal, CCI at the reception is usually modelled with random process that describes small scale signal variations in general nonlineofsight condition. Recently proposed [6] etamu () fading distribution meets those conditions. As mentioned in [7], model is presented in function of the parameter , related to the inphase/inquadrature components scattered wave powers quotient of each cluster, and parameter , related to the number of multipath clusters in the environment. As a general distribution, this model spans through other some wellknown fading models as its particular cases, that is, Hoyt model and Nakagami model. Namely, Nakagami distribution can be obtained in an exact manner by assigning parameter values and (cf. [6]).
Observing interferencelimited fading environment, where corruptive effects of thermal noise can be ignored, it is necessary to determine properties of signaltointerference ratio (SIR) in order to carry out effective performance analysis of observed system. Indeed, SIR has been measured in base and mobile stations by using SIR estimators [8]. In order to determine behavior of instantaneous SIR random process, we must determine square ratio of and random processes, , and analyze its properties. Here, exact closed expressions for the probability density function (PDF) and cumulative distribution function (CDF) of square ratio of and random processes are shown. Further, a verification of accuracy of these PDF and CDF expressions was given by comparing with the corresponding approximations obtained by the highprecision quadrature formulas of Gaussian type with respect to the weight functions on . The computational procedure of such quadrature rules is provided by using Gautschi’s constructive theory of orthogonal polynomials (see [9, 10]) and the MATHEMATICA package OrthogonalPolynomials created by Cvetković and Milovanović (see [11, 12]). Based on obtained expressions, rapidly converging infiniteseries expression for wireless communication systems performance measure, outage probability (OP) is also presented, and OP is considered for different values of transmission parameters. At the end, possible performance improvement will be considered through a glance at SC (selection combining) reception employment.
The paper is organized as follows. System model and closedform of the PDF and CDF expressions are presented in Section 2. Section 3 is devoted to an alternative approach based on the constructive theory of orthogonal polynomials and the corresponding quadrature rules of Gaussian type. System performances and numerical results are presented in Sections 4 and 5, respectively. Finally, concluding comments are given in Section 6.
2. System Model and ClosedForm of the PDF and CDF Expressions
The desired signal follows distributed random process modelled as (cf. [3]) where with , denoting desired signal average power, while being the th order modified Bessel function of the first kind [2, Eq. 8.445]. Due to physical separation between the cochannel interferer and target receiver, which is very large in practice, it is very unlikely that direct LOS condition (s) exist, for CCI signal. The most general fading model, which does not include LOS condition is fading model, is chosen here for modelling random envelope fluctuations of CCI signal with corresponding PDF [6]: where stands for the average power of CCI, while denotes Gamma function [2, Eq. 8.310.1]. Here and are interfering signal parameters, written in the function of parameter as follows [6]:
Now, we will observe random envelope statistics in observed interferencelimited system. Here the novelty and significance of this model must be pointed out again. Namely, for the first time scenario is observed when desired signal envelope variations are caused by fading occurrence (most general fading model which includes LOS component existence), while CCI signal envelope variations are caused by fading occurrence (most general fading model that does not include LOS component). Previously, various other scenarios have been observed in [1, 13–15]. In [1, 14] LOS component existence in CCI signal has been assumed, which is not so common case, because of large physical distance between the interferer and receiver. In [13] nonLOS model was observed for desired signal, while Gamma approximation of random variables sum has been observed. Finally, neither in general nor in comprehensive study, such as [15], has this model been considered. In the case when CCI is much stronger than Gaussian noise in same channel, PDF for the instantaneous SIR, , can be determined according to [16] in the form
After substituting (6) and (8) into (10), we can derive PDF of SIR in the following closedform:where and stands for the average SIR.
Namely, with respect to wellknown series representation of modified Bessel function (cf. [2, Eq. 8.445]) substituting (6) and (8) in (10) results in By using the wellknown definition of the Gamma function [2, Eq. 8.310.1]: (10) can be reduced into (11). A similar procedure has been used in [7]. By taking finite numbers of terms in the last expression for we obtain its approximation in a closedform: The closedform expression (16) converges rapidly, since only (10–15) terms are needed to be summed in each sum in order to reach an accuracy of five significant decimal digits, for observed set of system parameters. Otherwise, in a case of slowly convergent series, we need some methods for accelerating their convergence (cf. [17–19]).
Following [16], the cumulative distribution function (CDF) can be obtained as
After substituting (11) into (17), by applying the same mathematical transformations given in Appendix of [7], we can obtain the following closedform expression of rapid convergence: where and is the incomplete Beta function [2, Eq. 8.38] (for some extended special functions see [20]). Detailed convergence analysis of this expression is provided in Table 1.

3. Method of Quadrature Sums
In this section we give an alternative approach for getting very highprecision approximation of the probability density function (PDF), as well as one of the cumulative distribution function (CDF). Our goal is to obtain fast convergence of quadrature sums, , for integrals on halfinfinite ranges or finite intervals for functions of the form , where (≥0) is the socalled weight function and is a smooth factor if it is possible. We use here a computational procedure for constructing quadratures of Gaussian type with respect to the weight functions on , as well as ones on finite intervals, supposing the polynomial behavior of the function (for nonpolynomial cases see [21, 22], and for some nonstandard Gaussian quadratures see [23–25]). Our computational procedure provides the ability to verify the accuracy of expressions obtained in the previous section. As we will see, these quadrature processes give much better accuracy of the considered integrals.
As already mentioned PDF of this random process could be efficiently expressed by using Gaussian quadrature formulas, namely, starting from (10), which can be rewritten in the following form: with parameters , , and depending on (being functions of) normalized SIR expressed as
For calculating the weighted integral (20) we use the point generalized GaussLaguerre quadrature formula (cf. [26, pp. 324–328]) with respect to the generalized Laguerre weight function on , with a special selection of the parameter . The quadrature parameters, the nodes , and the weights (Christoffel numbers) can easily be constructed by the wellknown GolubWelsch procedure [27] implemented in the MATHEMATICA package OrthogonalPolynomials (see [11, 12]). This package is freely downloadable from the website http://www.mi.sanu.ac.rs/~gvm/.
In fact, the nodes are eigenvalues of the following symmetric tridiagonal Jacobi matrix (cf. [26, pp. 325–328]): and the weight coefficients are given by , , where is the first component of the eigenvector corresponding to the eigenvalue and normalized such that .
Remark 1. In some special cases for calculating the weight coefficients it is better to use the complete eigenvectors, instead of their first components. An analysis of such cases is given in [28] (see also [29, p. 339] for an application).
Elements in the Jacobi matrix (23) are recursion coefficients in the threeterm recurrence relation for polynomials orthogonal with respect to the generalized Laguerre weight function on : The coefficient can be arbitrary, but it is convenient to put in order to hold certain compact formulas, including the previously mentioned ones for Christoffel numbers , .
Remark 2. The recursion coefficients in the threeterm recurrence relation (24) are known in explicit form only for some narrow classes of orthogonal polynomials, among which are the most famous classical orthogonal polynomials (Jacobi, the generalized Laguerre, and Hermite polynomials) (cf. [26, pp. 121–146], [30], [31, pp. 19–24]). In the case of the socalled strongly nonclassical polynomials, the recursion coefficients must be constructed numerically using moment information on a given weight function (cf. [32–34]). For details on Gautschi’s constructive theory of orthogonal polynomials see [9, 10, 35], [26, pp. 160–166].
Knowing the first recursion coefficients in (24), we can obtain the nodes and the weights in (22) for any number of nodes .
It is known that if the integrand in the weighted integral (20) is smooth, then we can expect Gaussian quadrature applied to (20) to converge rapidly. According to the series representation (13), we see that is an entire function, and, therefore, the function is also entire. This means that, for calculating (20), this integral should be written in the following form: where and is defined by (27), and then quadrature formulas are applied with respect to the generalized GaussLaguerre weight with the parameter .
Then the sequence of quadrature sums , defined by shows a rapid convergence.
For example, for some typical parameters and , we present the graphic of in Figure 1 as an approximation of the function .
Relative errorsare presented in scale in Figure 2, for five different numbers of nodes in the quadrature sums (, , , , and ). As the exact PDF we take the approximation . For example, we can see that the quadrature sum with only six points gives an approximation of PDF for each , with the relative error less than ; that is, gives the result with at least six decimal digits, but the point quadrature sum approximates the PDF with at least eleven digits.
According to (17) the cumulative distribution function (CDF) can be transformed as For numerically calculating this integral we use the GaussLegendre quadrature rule (transformed to the interval ), whose parameters, the nodes and weights, and , , can be obtained from the eigenvalue problem for the corresponding Jacobi matrix of the form (23), with recursion coefficients In this symmetric case (with respect to ), it is clear that and , . For a sufficiently smooth function , the remaining term in (33) has the form [26, p. 333 & p. 97]
As before, for approximating , we use (29), so that we have that is, where is defined by (27).
The corresponding relative errors are presented in scale in Figure 3, for and in the quadrature sums (36), that is, (37). As the exact CDF we used the approximation .
As we can see, in the case when is large, this formula needs a bigger , and therefore the following approach is much more convenient. Namely, we use a decomposition of the interval into several subintervals for larger values of , for example, and then apply the GaussLegendre rule (33) to each of the subintervals. Thus, where the sums on the righthand side are empty for . The corresponding quadrature approximation, analogous to (36), becomes
Taking subintervals , the relative errors in the quadrature sums for are presented in Figure 4 for . We see that the last case, with four subintervals , , , and , gives the relative errors less than for each . The corresponding cumulative distribution function is presented in Figure 5.
4. System Performances
One of standard performance criteria, used for managing CCI level in faded channel, and for meeting gradeofservice and qualityofservice requirements, is OP [36]. When considering cellular systems, OP is defined as the probability that received SIR falls below a given outage threshold also known as a protection ratio
Diversity combining indicates a concept of upgrading reception system performances by combining replicas of information carrying signal [37]. By using separated antennas on reception diversity terminal, independent copies of desired signal are accomplished, without increase in signal bandwidth or transmitted signal power. Simplest space diversity combining method is selection combining (SC). In communication environments where the level of the CCI is significantly higher than the level of channel AWGN, SC outputs the antenna with the highest SIR. When reception system has uncorrelated branches (sufficient spacing between antennas), the PDF of the SIR at the output of the receiver is given by [38]: where denotes for the diversity reception order.
5. Numerical Results
At Figure 6 OP in the function of normalized SIR in dB (decibels, ) is presented. It is visible how smaller values of OP are obtained when parameter arises, since dominant component of desired signal then straightens. Also better performance (smaller values of OP) is obtained when fading, so that influences desired and CCI signal are less severe, namely, when parameters and take higher values. With SC reception performance improvement reached is also shown. Considerably lower OP values are obtained, comparing to the case with nodiversity reception. Finally it can be seen how balanced reception leads to achieving better quality of service. For example, it can be seen from the figure that when SC diversity has been used, for reaching the same level of OP (i.e., ), it is necessary to apply about 4 dB lower level power at the transmission, comparing to the level of transmitting power necessary to reach given OP level for the case of single branch reception when same values of propagation parameters are considered. Also, it can be seen from the figure that when the same level of transmitted signal power exists (i.e., 12 dB) higher level of outage (i.e., ) is obtained for single branch reception compared to case when SC diversity has been applied (i.e., ), when same values of propagation parameters are considered. Presented efficient numerical methods in such way could enable wireless system designers to perform tradeoff studies among the required values of system performances at the reception (quality of service) in order to determine the optimal choice of transmission parameters in the presence of their available constraints.
6. Conclusion and Future Work
In this paper the exact PDF and CDF expressions of the square ratio of and random processes were derived. Also, a verification of accuracy of these PDF and CDF expressions was given by comparing with the corresponding approximations obtained by the highprecision quadrature formulas of Gaussian type with respect to the weight functions on . The computational procedure of such quadrature rules is provided by using the contemporary constructive theory of orthogonal polynomials and the corresponding software. Our efficient numerical methods provide necessary control of level of error for the obtained approximation expressions. Also, the convergence of the proposed expressions is secured. Based on the provided analysis of this novel random process, interferencelimited wireless communication subjected to the influence of fading was observed when undesired occurrence of CCI has been modelled with random process. Rapid converging infiniteseries expressions are presented for received SIR’s PDF and CDF. Major communication performance measure, OP, has been calculated in terms of transmission parameters. OP upgrade accomplished by SC reception usage has also been analyzed.
Contribution presented in this paper provides potentials for applying proposed efficient numerical methods in carrying out analysis of other important wireless telecommunication performance measures for observed propagation scenario, such as channel capacity (CC) and average bit error probability (ABER). Closedform expressions for ABER values could be derived when both coherent and noncoherent modulation techniques are applied for wireless propagation over interferencelimited channels. In a similar manner CC values could be efficiently evaluated for proposed propagation scenario. By using provided highprecession approximation techniques that are based on presented quadrature formulas of Gaussian type, a verification of accuracy of these ABER and CC expressions could also be provided. Finally, presented approximations could find their application in many fields of telecommunication theory, since proposed solutions can be applied for performance analysis of large number of wireless transmission scenarios.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported in part by the Serbian Academy of Sciences and Arts (Project no. Φ96) and by the Serbian Ministry of Education, Science and Technological Development (Projects nos. OI 174015 and III 44006).
References
 S. Panić, M. Stefanović, J. Anastasov, and P. Spalević, Fading and Interference Mitigation in Wireless Communications, CRC Press, Taylor & Francis Group, USA, 2013.
 I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products Academic, New York, 1980. View at: MathSciNet
 J. C. S. Santos Filho and M. D. Yacoub, “Highly accurate κμ approximation to sum of M independent nonidentical Ricean variates,” IEEE Electronics Letters, vol. 41, no. 6, pp. 338339, 2005. View at: Publisher Site  Google Scholar
 S. R. Panic, D. M. Stefanović, I. M. Petrović, M. Č. Stefanović, J. A. Anastasov, and D. S. Krstić, “Secondorder statistics of selection macrodiversity system operating over Gamma shadowed κ–μ fading channels,” EURASIP Journal on Wireless Communications and Networking, vol. 2011, article 151, 2011. View at: Publisher Site  Google Scholar
 P. Stavroulakis, Interference Analysis and Reduction for Wireless Systems, Artech House, Inc., London, 2003.
 D. B. Da Costa and M. D. Yacoub, “The ημ Joint phaseenvelope distribution,” IEEE Antennas and Wireless Propagation Letters, vol. 6, pp. 195–198, 2007. View at: Publisher Site  Google Scholar
 P. Spalevic, M. Stefanovic, S. R. Panic, B. Milosevic, and Z. Trajcevski, “The cochannel interference effect on average error rates in Nakagamiq (Hoyt) fading channels,” Revue Roumaine des Sciences Techniques. Série Electrotechnique et Énergétique, vol. 56, no. 3, pp. 305–313, 2011. View at: Google Scholar
 M. Austin and G. Stuber, “Inservice signal quality estimation for TDMA cellular systems,” in Proceedings of the 6th International Symposium on Personal, Indoor and Mobile Radio Communications, pp. 836–840, Toronto, Ont., Canada. View at: Publisher Site  Google Scholar
 W. Gautschi, “On generating orthogonal polynomials,” SIAM Journal on Scientific Computing, vol. 3, no. 3, pp. 289–317, 1982. View at: Google Scholar  MathSciNet
 W. Gautschi, Orthogonal Polynomials, Oxford University Press, 2004. View at: MathSciNet
 A. S. Cvetković and G. V. Milovanović, “The Mathematica package OrthogonalPolynomials,” Facta Universitatis, Series: Mathematics and Informatics, no. 19, pp. 17–36, 2004. View at: Google Scholar  MathSciNet
 G. V. Milovanović and A. S. Cvetković, “Special classes of orthogonal polynomials and corresponding quadratures of Gaussian type,” Mathematica Balkanica, vol. 26, no. 12, pp. 169–184, 2012. View at: Google Scholar  MathSciNet
 K. Peppas, “Dualhop relaying communications with cochannel interference over ημ fading channels,” IEEE Transactions on Vehicular Technology, vol. 62, no. 8, pp. 4410–4115, 2013. View at: Google Scholar
 J. F. Paris, “Outage probability in ημ/ημ and κμ/ημ interferencelimited scenarios,” IEEE Transactions on Communications, vol. 61, no. 1, pp. 335–343, 2013. View at: Publisher Site  Google Scholar
 A. C. Moraes, D. B. Da Costa, and M. D. Yacoub, “An outage analysis of multibranch diversity receivers with cochannel interference in αμ, κμ, and ημ Fading scenarios,” Wireless Personal Communications, vol. 64, no. 1, pp. 3–19, 2012. View at: Publisher Site  Google Scholar
 M. Stefanovic, N. Sekulovic, M. Gligorijevic, M. Bandjur, P. Spalevic, and Z. Popovic, “Statistics of signal envelope in composite multipath fading/shadowing microcellular environment,” Technics Technologies Education Management, vol. 6, no. 4, pp. 1147–1151, 2011. View at: Google Scholar
 G. V. Milovanović, Numerical Analysis, University of Niš, Niš, 1979.
 G. V. Milovanović, “Summation of series and Gaussian quadratures,” in Approximation and Computation, R. V. M. Zahar, Ed., vol. 119, pp. 459–475, Birkhäuser Boston, Boston, MA, 1994. View at: Google Scholar  MathSciNet
 G. V. Milovanović, “On summation/integration methods for slowly convergent series,” Studia. Universitatis BabeşBolyai Mathematica, vol. 61, no. 3, pp. 359–375, 2016. View at: Google Scholar  MathSciNet
 M.J. Luo, G. V. Milovanovic, and P. Agarwal, “Some results on the extended beta and extended hypergeometric functions,” Applied Mathematics and Computation, vol. 248, pp. 631–651, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 G. V. Milovanović and A. S. Cvetković, “Gaussiantype quadrature rules for Müntz systems,” SIAM Journal on Scientific Computing, vol. 27, no. 3, pp. 893–913, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 G. V. Milovanović, “Müntz orthogonal polynomials and their numerical evaluation,” in Applications and computation of orthogonal polynomials, W. Gautschi, G. H. Golub, and G. Opfer, Eds., vol. 131, pp. 179–194, Birkhäuser, Basel, 1999. View at: Google Scholar  MathSciNet
 G. V. Milovanović and A. S. Cvetković, “Nonstandard Gaussian quadrature formulae based on operator values,” Advances in Computational Mathematics, vol. 32, no. 4, pp. 431–486, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 G. V. Milovanović and A. S. Cvetković, “Gaussian quadrature rules using function derivatives,” IMA Journal of Numerical Analysis (IMAJNA), vol. 31, no. 1, pp. 358–377, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 G. V. Milovanović, “Weighted quadrature formulas for semiinfinite range integrals,” Journal of Numerical Analysis and Approximation Theory, vol. 44, no. 1, pp. 69–80, 2015. View at: Google Scholar  MathSciNet
 G. Mastroianni and G. V. Milovanović, Interpolation Processes: Basic Theory and Applications, Springer, 2008. View at: MathSciNet
 G. H. Golub and J. H. Welsch, “Calculation of Gauss quadrature rules,” Mathematics of Computation, vol. 23, no. 106, article 221, 1969. View at: Publisher Site  Google Scholar  MathSciNet
 G. V. Milovanović and A. S. Cvetković, “Note on a construction of weights in Gausstype quadrature rule,” Facta Universitatis, Series: Mathematics and Informatics, no. 15, pp. 69–83, 2000. View at: Google Scholar  MathSciNet
 G. V. Milovanović, M. M. Spalević, and A. S. Cvetković, “Calculation of Gaussiantype quadratures with multiple nodes,” Mathematical and Computer Modelling, vol. 39, no. 23, pp. 325–347, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 R. P. Agarwal and G. V. Milovanović, “Extremal problems, inequalities, and classical orthogonal polynomials,” Applied Mathematics and Computation, vol. 128, no. 23, pp. 151–166, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 G. B. Djordjević and G. V. Milovanović, Special Classes of Polynomials, University of Niš, Faculty of Technology, Leskovac.
 G. V. Milovanović, “Construction and applications of Gaussian quadratures with nonclassical and exotic weight functions,” Studia. Universitatis BabeşBolyai Mathematica, vol. 60, no. 2, pp. 211–233, 2015. View at: Google Scholar  MathSciNet
 G. V. Milovanović, “Symbolicnumeric computation of orthogonal polynomials and Gaussian quadratures with respect to the cardinal Bspline,” Numerical Algorithms, vol. 76, no. 2, pp. 333–347, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 M. MasjedJamei and G. V. Milovanović, “Construction of Gaussian quadrature formulas for even weight functions,” Applicable Analysis and Discrete Mathematics, vol. 11, no. 1, pp. 177–198, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 W. Gautschi, Orthogonal Polynomials in Matlab: Exercises and Solutions, Software – Environments – Tools, SIAM, Philadelphia, PA, USA, 2016.
 M. Č. Stefanović, D. M. Milović, A. M. Mitić, and M. M. Jakovljević, “Performance analysis of system with selection combining over correlated Weibull fading channels in the presence of cochannel interference,” AEÜ  International Journal of Electronics and Communications, vol. 62, no. 9, pp. 695–700, 2008. View at: Publisher Site  Google Scholar
 D. Milic, S. Stanojcic, Z. Popovic, D. Stefanovic, and I. Petrovic, “Performance analysis of EGC combining over correlated Nakagamim fading channels,” Serbian Journal of Electrical Engineering, vol. 9, no. 2, pp. 293–300, 2012. View at: Publisher Site  Google Scholar
 M. K. Simon and M. S. Alouini, Digital Communication over Fading Channels, John Wiley & Sons, USA, 2000.
Copyright
Copyright © 2018 Gradimir V. Milovanović et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.