Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 4967613, 9 pages
https://doi.org/10.1155/2018/4967613

Research Article

Efficient Numerical Methods for Analysis of Square Ratio of κ-μ and η-μ Random Processes with Their Applications in Telecommunications

Gradimir V. Milovanović,1,2 Suad Suljović,³ Stefan R. Panić,⁴ Ismet Kalčo,⁵ and Mihajlo Č. Stefanović³

¹Serbian Academy of Sciences and Arts, Kneza Mihaila 35, 11000 Beograd, Serbia
²Faculty of Sciences and Mathematics, University of Niš, P.O. Box 224, 18000 Niš, Serbia
³Department of Telecommunications, Faculty of Electrical Engineering, University of Niš, 18000 Niš, Serbia
⁴Department of Informatics, Faculty of Natural Science and Mathematics, University of Priština, 383200 Kosovska Mitrovica, Serbia
⁵Faculty of Polytechnics, University of Zenica, Zenica, Bosnia and Herzegovina

Correspondence should be addressed to Stefan R. Panić; stefanpnc@yahoo.com

Received 12 October 2017; Accepted 30 January 2018; Published 11 March 2018

Academic Editor: Bo Shen

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.

Abstract.

We will provide statistical analysis of the square ratio of and random processes and its application in the signal-to-interference ratio (SIR) based performance analysis of wireless transmission subjected to the influence of multipath fading, modelled by fading model, and undesired occurrence of co-channel 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 high-precision 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 zero-the 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 Nakagami-m distribution models as special cases (as the One-Sided Gaussian and the Rayleigh distributions since they also represent special cases of Nakagami-m).

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), one-sided 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 short-term signal variation is kappa-mu () 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]. Co-channel 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 non-line-of-sight condition. Recently proposed [6] eta-mu () fading distribution meets those conditions. As mentioned in [7], model is presented in function of the parameter , related to the in-phase/in-quadrature 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 well-known 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 interference-limited fading environment, where corruptive effects of thermal noise can be ignored, it is necessary to determine properties of signal-to-interference 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 high-precision 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 infinite-series 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 closed-form 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 Closed-Form 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 co-channel 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 interference-limited 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] non-LOS 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 closed-form:

where

and stands for the average SIR.

Namely, with respect to well-known series representation of modified Bessel function (cf. [2, Eq. 8.445])

substituting (6) and (8) in (10) results in

By using the well-known 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 closed-form:

The closed-form 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 closed-form 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.

Table 1: Number of terms that should be added in expression (18) in order to reach accuracy at 5th significant digit.

 dB
 dB
 dB

3. Method of Quadrature Sums

In this section we give an alternative approach for getting very high-precision 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 half-infinite ranges or finite intervals for functions of the form , where (≥0) is the so-called 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 Gauss-Laguerre 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 well-known Golub-Welsch 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 three-term 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 three-term 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 so-called 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 Gauss-Laguerre 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 .

Figure 1: Approximation of PDF, , using nodes in (29).

Relative errors

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

Figure 2: Relative errors in obtained by the quadrature formula with nodes.

According to (17) the cumulative distribution function (CDF) can be transformed as

For numerically calculating this integral we use the Gauss-Legendre 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 .

Figure 3: Relative errors in the quadrature sums for and nodes.

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 Gauss-Legendre rule (33) to each of the subintervals. Thus,

where the sums on the right-hand 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.

Figure 4: Relative errors in the quadrature sums for , taking subintervals and .

Figure 5: Approximation of CDF, , for with four subintervals .

4. System Performances

One of standard performance criteria, used for managing CCI level in faded channel, and for meeting grade-of-service and quality-of-service 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 no-diversity 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 trade-off 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.

Figure 6: Outage probability versus normalized SIR.

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 high-precision 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, interference-limited wireless communication subjected to the influence of fading was observed when undesired occurrence of CCI has been modelled with random process. Rapid converging infinite-series 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). Closed-form expressions for ABER values could be derived when both coherent and noncoherent modulation techniques are applied for wireless propagation over interference-limited channels. In a similar manner CC values could be efficiently evaluated for proposed propagation scenario. By using provided high-precession 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.

J. C. S. Santos Filho and M. D. Yacoub, “Highly accurate κ-μ approximation to sum of M independent non-identical Ricean variates,” IEEE Electronics Letters, vol. 41, no. 6, pp. 338-339, 2005.

S. R. Panic, D. M. Stefanović, I. M. Petrović, M. Č. Stefanović, J. A. Anastasov, and D. S. Krstić, “Second-order statistics of selection macro-diversity system operating over Gamma shadowed κ–μ fading channels,” EURASIP Journal on Wireless Communications and Networking, vol. 2011, article 151, 2011.

P. Stavroulakis, Interference Analysis and Reduction for Wireless Systems, Artech House, Inc., London, 2003.

D. B. Da Costa and M. D. Yacoub, “The η-μ Joint phase-envelope distribution,” IEEE Antennas and Wireless Propagation Letters, vol. 6, pp. 195–198, 2007.

P. Spalevic, M. Stefanovic, S. R. Panic, B. Milosevic, and Z. Trajcevski, “The co-channel interference effect on average error rates in Nakagami-q (Hoyt) fading channels,” Revue Roumaine des Sciences Techniques. Série Electrotechnique et Énergétique, vol. 56, no. 3, pp. 305–313, 2011.

M. Austin and G. Stuber, “In-service 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.

W. Gautschi, “On generating orthogonal polynomials,” SIAM Journal on Scientific Computing, vol. 3, no. 3, pp. 289–317, 1982.

W. Gautschi, Orthogonal Polynomials, Oxford University Press, 2004.

A. S. Cvetković and G. V. Milovanović, “The Mathematica package OrthogonalPolynomials,” Facta Universitatis, Series: Mathematics and Informatics, no. 19, pp. 17–36, 2004.

G. V. Milovanović and A. S. Cvetković, “Special classes of orthogonal polynomials and corresponding quadratures of Gaussian type,” Mathematica Balkanica, vol. 26, no. 1-2, pp. 169–184, 2012.

K. Peppas, “Dual-hop relaying communications with co-channel interference over η-μ fading channels,” IEEE Transactions on Vehicular Technology, vol. 62, no. 8, pp. 4410–4115, 2013.

J. F. Paris, “Outage probability in η-μ/η-μ and κ-μ/η-μ interference-limited scenarios,” IEEE Transactions on Communications, vol. 61, no. 1, pp. 335–343, 2013.

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.

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.

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.

G. V. Milovanović, “On summation/integration methods for slowly convergent series,” Studia. Universitatis Babeş-Bolyai Mathematica, vol. 61, no. 3, pp. 359–375, 2016.

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.

G. V. Milovanović and A. S. Cvetković, “Gaussian-type quadrature rules for Müntz systems,” SIAM Journal on Scientific Computing, vol. 27, no. 3, pp. 893–913, 2005.

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.

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.

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.

G. V. Milovanović, “Weighted quadrature formulas for semi-infinite range integrals,” Journal of Numerical Analysis and Approximation Theory, vol. 44, no. 1, pp. 69–80, 2015.

G. Mastroianni and G. V. Milovanović, Interpolation Processes: Basic Theory and Applications, Springer, 2008.

G. H. Golub and J. H. Welsch, “Calculation of Gauss quadrature rules,” Mathematics of Computation, vol. 23, no. 106, article 221, 1969.

G. V. Milovanović and A. S. Cvetković, “Note on a construction of weights in Gauss-type quadrature rule,” Facta Universitatis, Series: Mathematics and Informatics, no. 15, pp. 69–83, 2000.

G. V. Milovanović, M. M. Spalević, and A. S. Cvetković, “Calculation of Gaussian-type quadratures with multiple nodes,” Mathematical and Computer Modelling, vol. 39, no. 2-3, pp. 325–347, 2004.

R. P. Agarwal and G. V. Milovanović, “Extremal problems, inequalities, and classical orthogonal polynomials,” Applied Mathematics and Computation, vol. 128, no. 2-3, pp. 151–166, 2002.

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.

G. V. Milovanović, “Symbolic-numeric computation of orthogonal polynomials and Gaussian quadratures with respect to the cardinal B-spline,” Numerical Algorithms, vol. 76, no. 2, pp. 333–347, 2017.

M. Masjed-Jamei 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.

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.

D. Milic, S. Stanojcic, Z. Popovic, D. Stefanovic, and I. Petrovic, “Performance analysis of EGC combining over correlated Nakagami-m fading channels,” Serbian Journal of Electrical Engineering, vol. 9, no. 2, pp. 293–300, 2012.

M. K. Simon and M. S. Alouini, Digital Communication over Fading Channels, John Wiley & Sons, USA, 2000.