Next Article in Journal
On Symmetrized Stochastic Harmonically Convexity and Hermite–Hadamard Type Inequalities
Next Article in Special Issue
Hilbert’s Double Series Theorem’s Extensions via the Mathieu Series Approach
Previous Article in Journal
Geometric Properties of Some Generalized Mathieu Power Series inside the Unit Disk
Previous Article in Special Issue
Fourier Transform of the Orthogonal Polynomials on the Unit Ball and Continuous Hahn Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Truncation Error in Statistical Description of Communication Signals over mm-Wave Channels

by
Zvezdan Marjanović
,
Dejan N. Milić
and
Goran T. Đorđević
*
Faculty of Electronic Engineering, University of Niš, A. Medvedeva 14, 18000 Niš, Serbia
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 569; https://doi.org/10.3390/axioms11100569
Submission received: 14 September 2022 / Revised: 7 October 2022 / Accepted: 12 October 2022 / Published: 19 October 2022
(This article belongs to the Special Issue Orthogonal Polynomials, Special Functions and Applications)

Abstract

:
This paper presents an illustration of how knowledge from the field of special functions, orthogonal polynomials and numerical series can be applied to solve a very important problem in the field of modern wireless communications. We present the formulas for the probability density function (PDF) and cumulative distribution function (CDF) of the composite signal envelope over an m m -Wave channel. The formulas for the PDF and CDF are expressed in the convergent infinity series form. The main contribution of the paper is in estimating the upper bounds for absolute truncation error in evaluating PDF and CDF of the signal envelope. We also derive the formulas for the required number of terms in the summation under the condition of achieving a given accuracy for typical values of channel parameters. In deriving these formulas, we use the alternating series estimation theorem, as well as some properties of orthogonal polynomials in order to derive upper bounds for hypergeometric functions. Based on the newly derived formulas, numerical results are presented and commented upon. The analytical results are verified by Monte Carlo simulations. The results are essential in the designing and performance estimating of the fifth-generation (5G) and beyond wireless networks.

1. Introduction

1.1. Motivation

In order to estimate accurately the performance of telecommunications systems, it is necessary to know the statistical characteristics of the signals that propagate through the telecommunication channel. The important metrics of a telecommunication system, such as error probability, ergodic capacity, and outage probability, can be determined by applying methods based on knowledge of the probability density function (PDF) and cumulative distribution function (CDF) of the detected signal envelope. It is desirable the PDF and CDF be expressed in a closed form. If it is not possible to derive these quantities in a closed form, then the aim is to express them in the form of convergent series. In the case when the PDF and CDF are known in the form of an infinite series, in order to calculate their numerical value, it is necessary to truncate the series. In those cases, it is important to be able to estimate the truncation error of these series [1,2]. For example, in wireless channels where the detected signal envelope is not constant, the error probability or the ergodic capacity is determined by integrating the PDF over all possible values of the signal envelope, so the accuracy of determining these quantities depends on the accuracy of the PDF calculation. In addition, the outage probability is determined directly on the basis of the CDF, so the calculation of the CDF is of exceptional importance for the determination of this quantity. Very often, both the PDF and CDF of the signal envelope over wireless channels contain different special functions and it is necessary to find some boundary values of these functions for different channel conditions, or it is necessary to perform a summation of these special functions when diversity techniques are applied at the receiver side [1,2].

1.2. Literature

High-frequency bands will be used in the fifth-generation (5G) and beyond wireless networks. In other words, signals will be transmitted over millimeter-wave ( m m -Wave) channels. Recently, it was shown experimentally that the statistical characteristics of signals transmitted in the band of 60 G Hz could be described by the two-wave diffuse-power (TWDP) model [3,4,5]. In addition, the authors of [6] illustrated that the fluctuating multiple-ray models had a good fit with experimental results in outdoor environments at 28 G Hz .
It is interesting that the TWDP model of fading over telecommunication channels was initially suggested in [7]. Namely, small-scale fading is characterized by the TWDP model whenever the received signal contains two strong (specular) multipath waves and many diffuse waves. An approximation for the PDF of the signal envelope variations was presented in [7], where the TWDP model was initially proposed. This model of fading has constantly been the subject of study by many researchers [8,9,10,11,12]. In Kim et al. [8], the authors emphasized some shortcomings of the approximation from [7], and they presented exact and approximate formulas for the bit-error rate when detecting a binary phase-shift keying signal transmitted over a TWDP channel for large values of the signal-to-noise ratio. In ref. [9], the authors derived some novel expressions for system performance metrics and presented interesting result showing that the TWDP fading model had a closed-form moment-generating function of the received signal envelope. In ref. [10], the joint estimation of the two parameters of the TWDP fading model was studied. A closed-form moment-based estimator was presented for fading parameters. The authors of [11,12] presented a novel way of parametrization of the TWDP channel model that was slightly different from the parametrization that had previously been utilized in [7,8,9,10]. One one hand, based on experimental results [3,4,5,6], this model of fading is suitable for m m -Wave channels. However, on the other hand, there is no closed-form solution for the PDF and CDF for the TWDP model. Consequently, any deeper theoretical analysis of these systems is hard because the PDF and CDF of the signal envelope are essential quantities for this analysis.

1.3. Contribution

In this paper, we briefly show that the statistical characterization of the TWDP channel can be performed by applying the equivalence of this model with a situation where there are a useful signal, cochannel interference and an additive white Gaussian noise (AWGN). I. Kostic [13,14] derived the PDF and CDF for the envelope of the composite signal consisting of a narrowband useful signal, cochannel interference and an AWGN. Very recently, the equivalence between that scenario and the case where the composite signal envelope variations are described by the TWDP model was noticed in [11]. However, in the numerical evaluation of the PDF and CDF, there is still the important open problem of how many terms should be taken in the summation in both formulas. Because of that, our emphasis is on the derivation of upper bounds for the absolute truncation error in infinite series in formulas for the PDF and CDF. The alternating series estimation theorem [15] is used in determining the bound for the truncation error. Furthermore, in order to derive an explicit formula linking the number of terms in the series we are summing up and the truncation error, it is necessary to determine the boundary values of the Bessel function of the first kind and order ν defined by ([16], (8.431)), as well as some bounds for hypergeometric functions 2 F 1 m , m ; 1 ; · and 1 F 1 1 m ; 2 ; · that can be presented as generalized hypergeometric series defined by ([17], (1)). Based on these newly derived expressions of the upper bounds for the truncation error and newly derived expressions for the required number of terms to achieve a given truncation error, we present some numerical values of the required number of terms in the summation in order to achieve the given accuracy for representative values of channel parameters. In addition, we verify the analytical results by independent Monte Carlo simulations.

1.4. Structure

In Section 2, we briefly describe the fading model and present the parameters definition, including the PDF and CDF of the signal envelope. Section 3 presents a convergence analysis of the series in the formulas for the PDF and CDF. For the series in both PDF and CDF formulas, we estimate the absolute truncation error, and derive the relation between the number of terms in the summation and the given absolute accuracy. Numerical results followed by appropriate discussions are given in Section 4, while Section 5 presents basic concluding remarks.

2. Physical Background

Based on the measurement verifications from [3,4,5,6] and the analytical models presented in [7,8,9,10,11,12], the resulting received complex signal envelope is given by [7,8,9,10,11,12]
r = V 1 e j Ψ 1 + V 2 e j Ψ 2 + x F + j y F ,
The resulting signal envelope consists of two specular components and a diffuse part. Specular components have constant amplitudes V 1 and V 2 , and uniformly distributed phases ( Ψ 1 and Ψ 2 ) in the interval from 0 to 2 π . The diffuse component has a Rayleigh distribution, i.e., it consists of the in-phase and quadrature components having a Gaussian distribution with zero mean value and standard deviation denoted by σ F . These in-phase and quadrature components are denoted by x F and y F , respectively.
The complex fading can be presented in terms of the envelope r and argument θ as
V r = r e j θ ,
where the fading envelope is given by:
r = V 1 cos Ψ 1 + V 2 cos Ψ 2 + x F 2 + V 1 sin Ψ 1 + V 2 sin Ψ 2 + y F 2
This model of fading can be described in terms of two parameters denoted by K and Δ . The parameter K denotes the ratio of the power of the specular components to the power of the diffuse component. The parameter Δ is related to the ratio of the peak specular components power to the average specular components’ power. These two parameters are defined as [7,8,9,10]
K = Average specular power Average diffuse power = V 1 2 + V 2 2 2 σ F 2 ,
Δ = Peak specular power Average specular power 1 = 2 V 1 V 2 V 1 2 + V 2 2 .
When K 0 + , the case corresponds to the situation without specular components, i.e., the resulting PDF is the Rayleigh one. Based on [7,8,9,10], the typical values of parameter K for terrestrial mobile links have values as high as 20 dB. The larger the value of the parameter K, the larger the power of the specular components compared with the power of the diffuse component. In other words, the fading is shallower. The values of parameter Δ lie in the range from zero to one. When Δ is equal to one, the specular components have an equal amplitude, while when Δ is equal to zero, either specular component’s amplitude is equal to zero.
An alternative way for parametrization was recently suggested in [11,12]. The authors of [11,12] used the same definition of parameter K, but instead of parameter Δ , they introduced parameter Γ defined simply as
Γ = V 2 / V 1
The values of parameter Γ lie in the range from zero to one. When Γ is equal to one, specular components have an equal amplitude, while when Γ is equal to zero, either specular component’s amplitude is equal to zero.
For our analysis presented in this paper, both parametrizations can be applied equally. We use the second one for illustrative purposes. Without loss of generality, we assume that the mean squared value of the envelope is equal to one, i.e.,
r 2 ¯ = V 1 2 + V 2 2 + 2 σ F 2 = 1
From (4), (6) and (7), it follows that
V 1 = K / ( 1 + Γ 2 ) ( 1 + K ) , V 2 = V 1 Γ ,
σ F = 1 / 2 ( 1 + K ) .
On the basis of the analytical derivations presented in [13,14], the PDF and CDF of the signal envelope given by (3) can be presented, respectively, in the forms of
p ( r ) = r σ F e r 2 + V 1 2 + V 2 2 2 σ F 2 [ I 0 r V 1 σ F 2 I 0 r V 2 σ F 2 I 0 V 1 V 2 σ F 2 + 2 m = 1 + ( 1 ) m I m r V 1 σ F 2 I m r V 2 σ F 2 I m V 1 V 2 σ F 2 ]
F ( r ) = r e r 2 2 σ F 2 2 σ F 2 m = 1 + ( 1 ) m m ! V 1 2 2 σ F 2 m 2 F 1 m , m ; 1 ; V 2 2 V 1 2 1 F 1 1 m ; 2 ; r 2 2 σ F 2 ,
where I ν ( · ) , ν = 0 , 1 , denotes a modified Bessel function of the first kind and order ν ([16], (8.431)), 2 F 1 · , · ; · ; · , and 1 F 1 · ; · ; · denote the hypergeometric functions that can be presented as the generalized hypergeometric series defined by ([17], (1)).
Formulas (10) and (11) were originally derived in [13,14], and here they are adapted to describe the composite signal envelope in the case where a signal is transmitted over an m m -Wave channel, which is relevant for terrestrial mobile systems [3,4,5,6].

3. Convergence Analysis of Series in PDF and CDF

In order to evaluate the numerical value of the PDF and CDF given by (10) and (11), respectively, it is necessary to take a finite number of terms in the summation in (10) and (11). Because of that, our aim is to estimate the truncation errors appearing after truncating the series in (10) and (11), as well as to give a direct relation between the number of terms in the summation on one hand, and give the absolute accuracy in the summation and channel parameters on the other hand.

3.1. Convergence Analysis of Series in PDF

Our goal is to estimate firstly the upper bound for the absolute truncation error of the series in (10) and secondly to estimate in advance the number of terms required for a specified accuracy level.

3.1.1. Truncation Error of Series in PDF

The alternating series estimation theorem bounds the absolute truncation error on the upper side to the absolute value of the first term left out after stopping the summation [15]. In estimating the truncation error of the series in (10), we are actually focusing on the series
S 1 = m = 1 + ( 1 ) m I m r V 1 σ F 2 I m r V 2 σ F 2 I m V 1 V 2 σ F 2 .
The truncation error of series S 1 is defined as
ε N P D F = m = N + 1 + ( 1 ) m I m r V 1 σ F 2 I m r V 2 σ F 2 I m V 1 V 2 σ F 2 .
Proposition 1.
The upper bound for the absolute truncation error ε N P D F defined by (13) is given by
ε N P D F R N P D F = 1 ( N + 1 ) ! 3 r V 1 V 2 2 2 σ F 3 2 ( N + 1 ) ,
where V 1 , V 2 and σ F are defined by (8) and (9).
Proof of Proposition 1.
In order to estimate the upper bound for the truncation error ε N P D F , we start by estimating the upper bound for the function I m ( · ) . Using the inequality from Luke ([18], (6.25)), and [19]
1 I m ( x ) e x m ! x 2 m ,
it directly holds that
I m r V 1 σ F 2 I m r V 2 σ F 2 I m V 1 V 2 σ F 2 1 ( m ! ) 3 exp r ( V 1 + V 2 ) + V 1 V 2 σ F 2 r V 1 V 2 2 2 σ F 3 2 m
Consequently the absolute truncation error is bounded as it is presented by (14). □
Corollary 1.
On the basis of (14), it is evident that lim N + R N P D F = 0 , which means the series S 1 given by (12) is convergent.

3.1.2. Required Number of Terms in Evaluating PDF

Proposition 2.
The required number of terms in the summation in (12) to achieve the absolute truncation error ε N P D F can be expressed as N = m
m 1 2 log 1 d ( ε N P D F ) W 0 1 2 a e log 1 d ( ε N P D F ) , d ( ε N P D F ) < e 2 a , d ( ε N P D F ) 1 a e , d ( ε N P D F ) = 1 a , d ( ε N P D F ) = e 2 a 0 , otherwise
where a and d ( ε N P D F ) are defined by (20) and (26), respectively, and W 0 ( · ) is the principal branch of the Lambert W function.
Proof of Proposition 2.
In estimating the required number of terms of S 1 to achieve a given truncation error ε N P D F , we need to find N = min { m } so that
1 ( m ! ) 3 exp r ( V 1 + V 2 ) + V 1 V 2 σ F 2 r V 1 V 2 2 2 σ F 3 2 m ε N P D F
The inequality (18) can be written in the simpler form
a m m ! b ,
where a and b are given by
a = r V 1 V 2 2 2 σ F 3 2 / 3 ,
b = ε N P D F 3 exp r ( V 1 + V 2 ) + V 1 V 2 3 σ F 2 .
Using Stirling’s formula for factorial ([16], 8.327(2)) m ! > 2 π m m e m , we rewrite the inequality (19) as
a m m ! < a e m m 1 2 π m b .
Since m 1 , we can write
a e m m 1 2 π m a e m m 1 2 π b .
Formula (23) introduces a significant simplification 1 1 / m , but also enables the expression of the result in explicit form. Therefore, we expect the value of m to be overestimated, which is a trade-off for simplicity. The rightmost inequality
a e m m b 2 π
can be rearranged to show the explicit value for m expressed in terms of the Lambert W function [20]
m 1 2 log 1 d ( ε N P D F ) W 0 1 2 a e log 1 d ( ε N P D F ) , d ( ε N P D F ) < e 2 a ,
where
d ( ε N P D F ) = 2 π ε N P D F 2 / 3 exp 2 r ( V 1 + V 2 ) + V 1 V 2 3 σ F 2 ,
and W 0 ( · ) stands for the principal branch of the Lambert W function (also known as product logarithm).
For d ( ε N P D F ) e 2 a , the limiting value is m a . On the other hand, in the case of d ( ε N P D F ) > e 2 a , inequality (24) is fulfilled for any m 0 , so there is effectively no need for the summation. Furthermore, for the sake of completeness, we notice that
m a e , when d ( ε N P D F ) 1 .
Based on the previously mentioned equation, the required number of terms in the summation in (10) can be expressed as in (17). □
While there are recent results that enable the precise computation [21] of the Lambert W function, there have also been some attempts to produce analytical approximations for the principal branch, such as [22,23]. For our purposes, we found that lower bound (32) from [23]
W 0 ( x ) x 1 + log ( 1 + x ) 1 + 3 x + x log ( 1 + x ) 1 + log 1 + 2 x 1 + log ( 1 + x )
provided a reasonable compromise between accuracy and complexity. Therefore, in numerical computation we used this equation as a stand-in replacement for the Lambert W 0 function.
For illustrative purposes, we present the required number of terms (N) given by (17) versus the truncation error ε N P D F in Figure 1. The numerical results presented in Figure 1 were evaluated for r = 0.5 , K = 8 dB and Γ = 0.5 .
On the other hand, this figure presents the required number of terms in the summation when the termination criterion is an absolute value of the first term left after terminating the summation to be less than ε N P D F . In other words, if we denote the ( N + 1 ) th term of the series by b N + 1 = ( 1 ) N + 1 I m r V 1 / σ F 2 I m r V 2 / σ F 2 I m V 1 V 2 / σ F 2 , we perform the summation until the condition | b N + 1 | ε N PDF is satisfied. It is clear that the required number of terms estimated on the basis of (17) is larger than the number of terms determined under the condition that the absolute value of the first term left after terminating the summation is less than ε N P D F . That means we have derived a bound on the truncation error which gives us the required number of terms in the summation to be sure that a given accuracy will be satisfied.
Based on the results from Figure 1, it can be noted that the same number of terms in the summation is needed to reach the bound for the truncation error in the range of several orders of magnitude. For example, 16 terms are necessary to ensure the bound on the truncation error from 2 × 10 15 to 3 × 10 13 This is the consequence of the number of terms being an integer value, which leads to the step-like behavior in Figure 1.

3.2. Convergence Analysis of Series in CDF

As in the case of the PDF considered in the previous subsection, in this subsection, we derive the upper bound for the absolute truncation error when evaluating the CDF given by (11) and we derive the formula for the required number of terms to achieve a given accuracy in the summation in (11).

3.2.1. Truncation Error of Series in CDF

Here, we can also apply the alternating series estimation theorem bounding the absolute truncation error on the upper side to the absolute value of the first term left out after stopping the summation [15]. In estimating the truncation error in (11), we focus on series
S 2 = m = 1 + ( 1 ) m m ! V 1 2 2 σ F 2 m 2 F 1 m , m ; 1 ; V 2 2 V 1 2 1 F 1 1 m ; 2 ; r 2 2 σ F 2
The truncation error of series S 2 is
ε N C D F = m = N + 1 + ( 1 ) m m ! V 1 2 2 σ F 2 m 2 F 1 m , m ; 1 ; V 2 2 V 1 2 1 F 1 1 m ; 2 ; r 2 2 σ F 2
Proposition 3.
The upper bound for the truncation error ε N C D F defined by (30) is given by
ε N C D F R N C D F = N 1 / 2 N + 1 [ 2 a ( 1 + b ) ] N + 1 ( N + 1 ) ! ( ( ( N + 1 ) ! ) 1 / N + c ) N ( N + 1 ) ! .
where the parameters are defined as a = V 1 2 / ( 2 σ F 2 ) , b = V 2 2 / V 1 2 and c = r 2 / ( 2 σ F 2 ) .
Proof of Proposition 3.
In order to estimate the upper bound for the truncation error and the required number of terms in the series to achieve a given absolute accuracy, we derive the bound for the hypergeometric functions appearing in (29).
• Derivation of the upper bound for function 2 F 1 m , m ; 1 ; z .
Based on ([24], (07.23.03.0195.01)), it holds that
2 F 1 m , m ; 1 ; z = ( 1 z ) m P m 1 + z 1 z
where P m ( x ) denotes the Legendre polynomial defined in ([16], (8.910)). We start the derivation of the bound for the hypergeometric function in (29) by proving the equality
| P m ( x ) | < a m | x | m ,
where a m denotes the coefficient with x m in the Legendre polynomial of order m. The coefficient a m can be expressed as ([25], p. 40)
a m = ( 2 m ) ! 2 m ( m ! ) 2 = 2 m m 1 / 2 m
Since the zeros of Legendre polynomials are real and symmetrical, we write
P m ( x ) = a m x ( 1 + ( 1 ) m + 1 ) / 2 i = 0 m / 2 ( x 2 α i 2 ) ,
where α i are zeros of the Legendre polynomial of order m. The polynomials of odd order have an additional zero which is located at x = 0 . All zeros are located inside the interval ( 1 , 1 ) , implying 0 α i 2 < 1 and therefore | x 2 α i 2 | < x 2 when x > 1 , so we can directly write
| P m ( x ) | < a m | x | m , | x | > 1 , m > 1
Since 1 P m ( x ) 1 , for 1 x 1 , we can extend the inequality to arbitrary real x in the following way
| P m ( x ) | < 1 + a m | x | m , m > 1
Since ( 1 + x ) / ( 1 x ) 1 when 0 x < 1 , we use the first version (36) of the inequality, yielding
2 F 1 m , m ; 1 ; x a m ( 1 + x ) m , x 0 ,
where equality is achieved at x = 1. This derived inequality (38) is illustrated in Figure 2.
• Derivation of the upper bound for function 1 F 1 1 m ; 2 ; z .
A confluent hypergeometric function can be represented as ([24], (07.03.26.0002.01), (07.20.27.0001.01))
1 F 1 1 m ; 2 ; z = 1 m L m 1 1 ( z ) .
where L m α ( z ) represents the Laguerre function [16]. Love noted in [26] that there were only a few results published about Laguerre polynomials’ inequalities. Our result in this context is stated as:
1 m L m 1 1 ( z ) 1 m ! ( α ( m 1 ) + z ) m 1 = b m 1 ( 1 ) ( z ) m , α ( m 1 ) = ( m ! ) 1 / ( m 1 ) .
Using ([27], (5.17)), we get
L m 1 ( 0 ) = m + 1 .
By substituting z = 0 into the previous inequality, we get
L m 1 1 ( 0 ) = m 1 ( m 1 ) ! α ( m 1 ) m 1 ,
which yields an appropriate value for α ( m ) , as
α ( m ) = ( ( m + 1 ) ! ) 1 / m .
All the zeros of Laguerre polynomials are positive, so we can prove
L m 1 ( z ) < 1 m ! z m ,
when z is larger than the largest zero of L m 1 ( z ) . Therefore,
1 m L m 1 1 ( z ) 1 m ! ( α ( m 1 ) + z ) m 1 , z > z 1 ,
where z 1 is the largest zero of L m 1 1 ( z ) . Another bound can be proven for Laguerre polynomials, ([28], 10.18(14)), ([29], (18.14.8))
L m 1 ( x ) ( m + 1 ) exp x 2 = b m ( 2 ) ( x ) , z > 0 .
An illustration of the two bounds is shown in Figure 3. It can be seen that the bound b m ( 1 ) ( z ) is valid, but not as sharp as b m ( 2 ) ( z ) for smaller values of z that are below the largest zero. On the other hand, for large z, the bounds’ sharpness is reversed.
On the basis of the previously derived formulas and the alternative series truncation error theorem, Proposition 3 is proven. □
Corollary 2.
On the basis of (31), it is evident that lim x + R N C D F = 0 , which means the series S 2 given by (29) is convergent.
Regarding Figure 3, the crossing point x 0 between the two bounds can be used in the following statement:
bound = bound 2 , 0 x < x 0 bound 1 , x 0 x < +
To date, we have not seen in the literature this exact combination of bounds addressing the [ 0 , ) range. Additionally, the crossing point can be expressed in closed*form as
x 0 2 m = β m W 1 β m e β m ,
where β m = ( ( 1 + m ) ! ) 1 / m 2 m .

3.2.2. Required Number of Terms for Evaluating CDF

In this subsection, we derive two formulas giving an explicit relation between the required number of terms in the summation, and the value of the absolute truncation error. One equation is marked as “Upper bound 1”, and the other is marked as “Upper bound 2”.
Proposition 4.
The required number ( N 1 ) of terms in the summation in (29) under the condition of achieving absolute truncation error ε N C D F can be determined by N = m
m 2 e c + log 1 2 π 2 ε N C D F 2 2 W 1 4 a ( 1 + b ) e ( 2 e c + log 1 2 π 2 ε N C D F 2 ) ε N C D F < ε 1 = 1 π 1 2 e 4 a ( 1 + b ) + 2 e c 1 / 2 2 a ( 1 + b ) , ε N C D F = ε 1 0 , otherwise ,
where a = V 1 2 2 σ F 2 , b = V 2 2 V 1 2 and c = r 2 2 σ F 2 .
Proof of Proposition 4.
By regarding the generalized binomial coefficient as
m n = Γ ( m + 1 ) Γ ( n + 1 ) Γ ( m n + 1 ) m , n , m n Z 0 ,
and expanding m 1 / 2 m into a series around m + , we get
m 1 / 2 m = 1 m π + O ( m 3 / 2 ) .
It can be proved that the following holds
1 m π m 1 / 2 m
Unfortunately, after the previous simplification, (31) is still too complex for symbolic manipulation. In order to obtain results with a reasonable complexity, we use the limit
lim m + ( ( m ! ) 1 / ( m 1 ) + c ) m 1 m ! = exp ( e c ) .
Our upper bound is now
R m 1 C D F 2 a ( 1 + b ) m m ! exp ( e c ) m π < 2 a ( 1 + b ) m m ! exp ( e c ) π , m > 1 .
This leads to an equation for m given by
exp ( e c ) π 2 2 a ( 1 + b ) e m m ε .
The variable m can be expressed similarly as it was in the case of the PDF. The solution is given by (49). We note that the limiting value is
lim ε ε 1 2 e c + log 1 2 π 2 ε 2 2 W 1 4 a ( 1 + b ) e 2 e c + log 1 2 π 2 ε 2 = 2 a ( 1 + b ) ,
and for ε > ε 1 , the required m is m 0 . □
Proposition 5.
The required number ( N 2 ) of terms in the summation in (29) under the condition of achieving the truncation error can be determined by N = m
m c + log 1 2 π 2 ( ε N C D F ) 2 2 W 1 4 a ( 1 + b ) e ( c + log 1 2 π 2 ( ε N C D F ) 2 ) ε N C D F < ε 2 = 1 π 1 2 e 4 a ( 1 + b ) + c 1 / 2 2 a ( 1 + b ) , ε N C D F > ε 2 0 , otherwise ,
where a = V 1 2 2 σ F 2 , b = V 2 2 V 1 2 and c = r 2 2 σ F 2 . This solution is denoted by “Upper bound 2”.
Proof of Proposition 5.
By direct use of (46), it follows that
a m m ! 2 F 1 m , m ; 1 ; b 1 F 1 1 m ; 2 ; c < m 1 / 2 m [ 2 a ( 1 + b ) ] m m ! exp ( c / 2 ) ,
From inequality (58), (57) follows. □
By comparing (49) with (57), we can notice that the latter is obtained by substituting c with c / ( 2 e ) in the former.
Figure 4 presents the required number of terms (both N 1 and N 2 ) needed to achieve a given absolute truncation error. This figure also presents the required number of terms that is determined from the criterion that the summation be terminated when the absolute value of the next term is less than ε N C D F .
Figure 4 illustrates the required number of terms in the sum in (29). Firstly, we estimated the number of terms numerically in a way that we terminated the summation when the absolute value of the first term left after the summation was less than the accuracy level denoted by ε N C D F . Secondly, we determined the number of terms was the summation based on (49), and that value was denoted as “Upper bound 1”. Thirdly, we estimated the number of terms based on (57). That bound was denoted as “Upper bound 2”. It is clear that the values that were estimated numerically are closer to the values estimated based on (57) compared with those estimated based on (49).

4. Numerical Results and Discussion

In the previous section, we derived new closed-form expressions ((13) and (30)) for the upper bounds of the truncation errors when evaluating the PDF and CDF of the signal envelope. These expressions were derived in terms of elementary mathematical functions. In addition, we derived closed-form Formulas (17), (49) and (57) giving a straightforward relation between the number of terms in the summation on one hand, and the value of the truncation error and channel parameters on the other hand. In this section, we give some illustrations of applications of these formulas.
Figure 5 and Figure 6 present the numerical values of the upper bounds for the truncation error of the PDF under different channel conditions. Figure 5a presents the upper bound for the truncation error vs. the required number of terms in the summation to achieve this upper bound. We present several curves corresponding to different values of parameter K. The larger the value of parameter K, the larger the number of terms necessary to achieve the given bound for the truncation error. For example, to achieve the bound for a truncation error of 10 12 , the required number of terms in the summation increases from N = 8 to N = 21 when parameter K increases from K = 4 dB to K = 12 dB. In other words, a larger number of terms is required in channels with stronger specular components compared with the diffuse component. Figure 5b presents the dependence of the required number of terms in the summation vs. the truncation error. This dependence is given by (17). For values of the truncation error of 10 12 , the required numbers of terms according to (17) are 10 and 23. We can conclude that relation (17) gives numbers of terms larger than those numbers that can be noticed from the dependence in Figure 5a, which were obtained based on (14).
The influence of parameter Γ on the number of required terms for achieving a given upper bound for the truncation error is illustrated in Figure 6. The larger the value of parameter Γ , the larger the number of terms in the summation required to achieve a given upper bound for the truncation error. If we fix the bound for the truncation error to 10 12 , from Figure 6a, we can read that it is necessary to sum N = 10 and N = 15 terms if Γ = 0.2 and Γ = 1 , respectively. When applying (17), i.e., from the results presented in Figure 6b, we conclude that for a truncation error of 10 12 , relation (17) gives the estimation that it is necessary to have 11 and 16 terms in the summation. As it was highlighted during the derivation, (17) gives overestimated values compared with those values that can be read from the dependence defined by (14).
The effect of channel conditions on the required terms when evaluating a CDF value is illustrated in Figure 7. The high values of parameter K mean that a larger number of terms is required, as in the case of evaluating a PDF value. As in the case of evaluating PDF, the same conclusions are derived from the results presented in the figure. Relation (57), linking the required number of terms in the summation vs. the truncation error, gives overestimated values of the number of terms compared with those values that can be read from the dependence of the bounds for the truncation error vs. the number of terms, which are defined by (31). When comparing the results from Figure 5 and Figure 7, we can conclude that a larger number of terms is required when evaluating the CDF in comparison with the number of terms when evaluating the PDF, for the same conditions.
Figure 8 presents PDF curves that were evaluated based on (10), when the number of terms in the summation was estimated under the condition that the upper bound of the truncation error be less than 10 12 . The PDF curves obtained by simulations are also presented in Figure 8. These simulation results were estimated based on samples generated according to (3). The samples were generated in the software package Mathematica 13. Uniform and Gaussian random numbers were generated by using built-in subroutines in Mathematica 13. When showing the histograms of the Monte Carlo data, we used linearly distributed bins with a total bin count of 200, and a total number of simulation points of 10 6 . The standard deviation of the histogram relative to numerically computed values of the distribution was σ = 0.008 for the 10 dB case and σ = 0.0062 for 4 dB. It is visible that the analytical and simulation results are in agreement, which is a confirmation that the suggested estimation of the required number of terms in the summation presented in this paper is correct.

5. Conclusions

By applying the alternating series estimation theorem, we derived novel formulas for the upper bounds for the absolute truncation error when evaluating the PDF and CDF of a composite signal envelope transmitting over m m -Wave channels. In addition, we derived novel analytical expressions giving the required number of terms in the infinity sums in the PDF and CDF under the condition of achieving a given absolute accuracy. In deriving these main formulas, we first derived two formulas for the upper bounds for the special functions 2 F 1 m , m ; 1 ; · and 1 F 1 1 m ; 2 ; · , where some properties of orthogonal polynomials were used. All the formulas were derived in terms of elementary mathematical functions. Based on the results presented here, both series (in the PDF and CDF) converged fast and for all practical channel parameters, it was enough to evaluate up to 100 terms in order to achieve an absolute accuracy of 10 12 . A larger number of terms was required for evaluating the CDF in comparison with the number of terms required for evaluating the PDF, for the same channel conditions. To verify the analytical results, when evaluating PDF curves, we used in the summation the required number of terms predicted by the formulas derived here. These analytical curves were verified by independent Monte Carlo simulations. The results will have an important role in contemporary wireless communications theory for the design and analysis of beyond-5G networks.

Author Contributions

Conceptualization, G.T.Đ.; methodology, D.N.M., Z.M. and G.T.Đ.; software, D.N.M. and G.T.Đ.; validation, D.N.M., Z.M. and G.T.Đ.; formal analysis, D.N.M., Z.M. and G.T.Đ.; investigation, D.N.M., Z.M. and G.T.Đ.; resources, D.N.M. and G.T.Đ.; data curation, D.N.M. and G.T.Đ.; writing—original draft preparation, D.N.M., Z.M. and G.T.Đ.; writing—review and editing, D.N.M. and G.T.Đ.; visualization, D.N.M. and G.T.Đ.; supervision, G.T.Đ.; project administration, G.T.Đ.; funding acquisition, Z.M., D.N.M. and G.T.Đ. All authors have read and agreed to the published version of the manuscript.

Funding

The work of G. T. Đorđević was supported by the Science Fund of the Republic of Serbia under grant no. 7750284, Hybrid Integrated Satellite and Terrestrial Access Network—hi-STAR. The work of Đorđević, D. N. Milić and Z. Marjanović was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia.

Data Availability Statement

Not applicable.

Acknowledgments

G. T. Đorđević would like to thank to Ivo M. Kostic for many valuable discussions in the general field of wireless communications theory.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

5GFifth generation
AWGNAdditive white Gaussian noise
CDFCumulative distribution function
PDFProbability density function
TWDPTwo-wave diffuse power

References

  1. Simon, M.K.; Alouini, M.S. Digital Communication over Fading Channels, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2004. [Google Scholar] [CrossRef]
  2. Shankar, P.M. Fading and Shadowing in Wireless Systems, 2nd ed.; Springer: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  3. Mavridis, T.; Petrillo, L.; Sarrazin, J.; Benlarbi-Delai, A.; De Doncker, P. Near-body shadowing analysis at 60 GHz. IEEE Trans. Antennas Propag. 2015, 63, 4505–4511. [Google Scholar] [CrossRef]
  4. Kim, D.; Lee, H.; Kang, J. Comments on “Near-body shadowing analysis at 60 GHz”. IEEE Trans. Antennas Propag. 2017, 65, 3314. [Google Scholar] [CrossRef]
  5. Zöchmann, E.; Caban, S.; Mecklenbräuker, C.F.; Pratschner, S.; Lerch, M.; Schwarz, S.; Rupp, M. Better than Rician: Modelling millimetre wave channels as two-wave with diffuse power. EURASIP J. Wirel. Commun. Netw. 2019, 2019, 1–17. [Google Scholar] [CrossRef] [Green Version]
  6. Sánchez, J.D.V.; Urquiza-Aguiar, L.; Paredes Paredes, M.C. Fading channel models for mm-wave communications. Electronics 2021, 10, 798. [Google Scholar] [CrossRef]
  7. Durgin, G.D.; Rappaport, T.S.; De Wolf, D.A. New analytical models and probability density functions for fading in wireless communications. IEEE Trans. Commun. 2002, 50, 1005–1015. [Google Scholar] [CrossRef] [Green Version]
  8. Kim, D.; Lee, H.; Kang, J. Comprehensive analysis of the impact of TWDP fading on the achievable error rate performance of BPSK signaling. IEICE Trans. Commun. 2018, 101, 500–507. [Google Scholar] [CrossRef]
  9. Rao, M.; Lopez-Martinez, F.J.; Alouini, M.S.; Goldsmith, A. MGF approach to the analysis of generalized two-ray fading models. IEEE Trans. Wirel. Commun. 2015, 14, 2548–2561. [Google Scholar] [CrossRef] [Green Version]
  10. Lopez-Fernandez, J.; Moreno-Pozas, L.; Lopez-Martinez, F.J.; Martos-Naya, E. Joint parameter estimation for the two-wave with diffuse power fading model. Sensors 2016, 16, 1014. [Google Scholar] [CrossRef] [PubMed]
  11. Maric, A.; Kaljic, E.; Njemcevic, P. An alternative statistical characterization of TWDP fading model. Sensors 2021, 21, 7513. [Google Scholar] [CrossRef] [PubMed]
  12. Njemcevic, P.; Kaljic, E.; Maric, A. Moment-Based Parameter Estimation for the Γ-Parameterized TWDP Model. Sensors 2022, 22, 774. [Google Scholar] [CrossRef] [PubMed]
  13. Kostic, I. Envelope probability density function of the sum of signal, noise and interference. Electron. Lett. 1978, 14, 490–491. [Google Scholar] [CrossRef]
  14. Kostic, I. Cumulative distribution function of envelope of sum of signal, noise and interference. In Proceedings of the Telecommunication Forum (TELFOR), Belgrade, Yugoslavia, 26–28 November 1996; pp. 301–303. [Google Scholar]
  15. Milovanovic, G.V. Numerical Analysis, Part I; University of Nis: Nis, Yugoslavia, 1979. (In Serbian) [Google Scholar]
  16. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 6th ed.; Academic Press: New York, NY, USA, 2000. [Google Scholar]
  17. Choi, J.; Milovanović, G.V.; Rathie, A.K. Generalized summation formulas for the Kampé de Fériet function. Axioms 2021, 10, 318. [Google Scholar] [CrossRef]
  18. Luke, Y.L. Inequalities for generalized hypergeometric functions. J. Approx. Theory 1972, 5, 41–65. [Google Scholar] [CrossRef] [Green Version]
  19. Joshi, C.; Bissu, S. Some inequalities of Bessel and modified Bessel functions. J. Aust. Math. Soc. 1991, 50, 333–342. [Google Scholar] [CrossRef] [Green Version]
  20. Corless, R.M.; Gonnet, G.H.; Hare, D.E.; Jeffrey, D.J.; Knuth, D.E. On the LambertW function. Adv. Comput. Math. 1996, 5, 329–359. [Google Scholar] [CrossRef]
  21. Lóczi, L. Guaranteed-and high-precision evaluation of the Lambert W function. Appl. Math. Comput. 2022, 433, 127406. [Google Scholar] [CrossRef]
  22. Iacono, R.; Boyd, J.P. New approximations to the principal real-valued branch of the Lambert W-function. Adv. Comput. Math. 2017, 43, 1403–1436. [Google Scholar] [CrossRef]
  23. Howard, R.M. Analytical Approximations for the Principal Branch of the Lambert W Function. Eur. J. Math. Anal. 2022, 2, 14. [Google Scholar] [CrossRef]
  24. Wolfram Research, Inc. The Mathematical Functions Site. 1998–2022. Available online: http://functions.wolfram.com (accessed on 15 August 2022).
  25. Milovanovic, G.; Mitrinovic, D.; Rassias, T. Topics in Polynomials: Extremal Problems, Inequalities, Zeros; World Scientific: Singapore, 1994. [Google Scholar] [CrossRef]
  26. Love, E. Inequalities for Laguerre functions. J. Inequalities Appl. 1997, 1997, 936095. [Google Scholar] [CrossRef]
  27. Szegö, G. Orthogonal Polynomials, 4th ed.; American Mathematical Society, Colloquium Publications: Providence, RI, USA, 1975; Volume 23. [Google Scholar] [CrossRef]
  28. Erldelyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1955; Volume 2. [Google Scholar]
  29. NIST Digital Library of Mathematical Functions; Release 1.1.6 of 2022-06-30; Olver, F.W.J.; Daalhuis, A.B.O.; Lozier, D.W.; Schneider, B.I.; Boisvert, R.F.; Clark, C.W.; Miller, B.R.; Saunders, B.V.; Cohl, H.S.; McClain, M.A. (Eds.) . Available online: http://dlmf.nist.gov/ (accessed on 30 June 2022.).
Figure 1. Required number of terms computed numerically and using the bound given by (17).
Figure 1. Required number of terms computed numerically and using the bound given by (17).
Axioms 11 00569 g001
Figure 2. Simple upper bound for 2 F 1 ( m , m ; 1 ; x ) , for m = 20 .
Figure 2. Simple upper bound for 2 F 1 ( m , m ; 1 ; x ) , for m = 20 .
Axioms 11 00569 g002
Figure 3. Two upper bounds for | 1 F 1 ( 1 m ; 2 ; x ) | , for m = 20 .
Figure 3. Two upper bounds for | 1 F 1 ( 1 m ; 2 ; x ) | , for m = 20 .
Axioms 11 00569 g003
Figure 4. Illustration of the required number of terms in the summation in the CDF to achieve a specified absolute truncation error.
Figure 4. Illustration of the required number of terms in the summation in the CDF to achieve a specified absolute truncation error.
Axioms 11 00569 g004
Figure 5. Effect of parameter K on accuracy of evaluating PDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Figure 5. Effect of parameter K on accuracy of evaluating PDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Axioms 11 00569 g005
Figure 6. Effect of parameter Γ on accuracy of evaluating PDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Figure 6. Effect of parameter Γ on accuracy of evaluating PDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Axioms 11 00569 g006
Figure 7. Effect of parameter K on accuracy of evaluating CDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Figure 7. Effect of parameter K on accuracy of evaluating CDF: (a) upper bound for truncation error vs. number of terms in summation; (b) required number of terms vs. truncation error.
Axioms 11 00569 g007
Figure 8. Analytical and Monte Carlo simulation results.
Figure 8. Analytical and Monte Carlo simulation results.
Axioms 11 00569 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marjanović, Z.; Milić, D.N.; Đorđević, G.T. Estimation of Truncation Error in Statistical Description of Communication Signals over mm-Wave Channels. Axioms 2022, 11, 569. https://doi.org/10.3390/axioms11100569

AMA Style

Marjanović Z, Milić DN, Đorđević GT. Estimation of Truncation Error in Statistical Description of Communication Signals over mm-Wave Channels. Axioms. 2022; 11(10):569. https://doi.org/10.3390/axioms11100569

Chicago/Turabian Style

Marjanović, Zvezdan, Dejan N. Milić, and Goran T. Đorđević. 2022. "Estimation of Truncation Error in Statistical Description of Communication Signals over mm-Wave Channels" Axioms 11, no. 10: 569. https://doi.org/10.3390/axioms11100569

APA Style

Marjanović, Z., Milić, D. N., & Đorđević, G. T. (2022). Estimation of Truncation Error in Statistical Description of Communication Signals over mm-Wave Channels. Axioms, 11(10), 569. https://doi.org/10.3390/axioms11100569

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