# Novel Method for 5G Systems NLOS Channels Parameter Estimation

Authors: Vladeta Milenkovic; Stefan Panic; Dragan Denic; Dragan Radenkovic
Journal: International Journal of Antennas and Propagation
 (2017)

---

## Abstract

For the development of new 5G systems to operate in mm bands, there is a need for accurate radio propagation modelling at these bands. In this paper novel approach for NLOS channels parameter estimation will be presented. Estimation will be performed based on LCR performance measure, which will enable us to estimate propagation parameters in real time and to avoid weaknesses of ML and moment method estimation approaches.


---

## Body

## 1. Introduction

Exploitation of unused mm wave spectrum (spectrum between 6 and 300 GHz) is an efficient solution for meeting the standards for 5G networks enormous data demand growth explosion. Because of that characterization and modelling of such channel propagation in urban environments is one of most important tasks in developing novel 5G mobile access networks. Many propagation studies, performed at these bands, for these types of applications, consider line-of-sight (LOS) scenarios [1–3]. However, non-line-of-sight (NLOS) scenarios occur more often. They not only occur in cases when transmitting and receiving antenna are separated by obstructions and there is no clear optical path between the antennas, but also occur in cases when there is indeed clear optical path between the antennas, but antennas are not aligned or boresight [4]. Various outdoor propagation measurements conducted in urban environments have shown that, for reasonable level of signal-to-noise ratio (SNR), that is, higher than 5 dB, NLOS signals have stable first arriving signal levels, even if they are weaker than stronger, later-arriving multipath components, counter to LOS signals with strongest first arriving component at the reception. Measurements provided at 38 GHz (base station-to-mobile access scenario [5] and Peer-to-Peer Scenario [6]), 60 GHz (Peer-to-Peer Scenario and vehicular scenario [7]), and 73 GHz [8] have clearly identified existence of NLOS conditions.One of the most intensively used statistical models for characterizing the complex behavior and random nature of NLOS fading envelope is the Nakagami-m distribution. Nakagami-m fading model closely approximates data values that are obtained by providing real-time measurements in indoor and outdoor wireless environments and is often used in analyzing propagation performances in 5G systems [9–14]. In [15–17] for the purpose of modelling observed 5G system propagation properties the Nakagami m parameter is directly computed from the measured data. As a general fading distribution, Nakagami-m fading model includes in itself other distributions and its simple family form allows obtaining closed-form analytical results for wireless communication link standard performance criterion measures. Namely, Nakagami-m fading model can also be transformed into Ricean fading model, by expressing m parameter in the function of Ricean K factor, as m=(1+K)2/1+2K [18]. In that manner characterizing of the behavior of LOS fading envelope could be also performed. The Nakagami-m fading exploits Nakagami probability density function (PDF) for the random envelope of received signal which is written in the function of two parameters: scale parameter and shape parameter, called the fading severity parameter or m-parameter. Determining m is a problem in Nakagami PDF estimation. For observed set of empirical fading signal data, the value of distribution parameter m should be estimated from it, in order to use the Nakagami-m distribution to model a given set of values. Acquaintance of the m parameter is mandatory for the optimal reception of signals in Nakagami-m fading environment [19], and in such case parameter m should be determined very accurately and very fast. Acquaintance of the m parameter is also necessary in the transmitter adaptation process, where its value could be feedback with respect to the channel information. Two most well-known procedures used for the estimation of the Nakagami-m fading parameter, m, are (1) maximum likelihood (ML) estimation and (2) moment-based estimation. ML-based m parameter estimation problem reduces to the problem of solving some transcendental equations written in the form of logarithmic and digamma functions. During the process of Gamma shape parameter estimation authors of [20] have developed most famous ML-based estimator of m parameter. Recently, another two approximate ML-based estimators have been proposed in [21], where estimation is carried out by observing approximations of digamma function, that is, first-order approximation and second-order approximation. Opposite approach to the ML estimation of m parameter is the moment-based estimation. Observing the second and the fourth Nakagami-m sample moments, estimation of parameter m was carried out in [22]. An improvement of proposed method can be found in [23] where parameter m estimation is carried on capitalizing on the first and the third sample moments. Group of new moment estimators form based on noninteger sample moments (estimators based on integer sample moment are their special cases), along with simulation study, was proposed in [24]. The generalized method of moments (GMM) is introduced for the estimation of the Nakagami-m fading parameter in [25]. However, it is known that sample moments are often subjected to the effects of outliers (even a small portion of extreme values, outliers, can affect the Gaussian parameters, especially the higher order moments). Moreover, occurrence of outliers is especially problematic when higher order sample moments are used for estimation, since estimation inaccuracy arises in such cases. Providing the best moment-based estimator is still major issue that should be addressed. Specifically, the level crossing rate (LCR), which provides us with a measure of the average number of crossings per second at which the envelope crosses a specified signal level in positive or negative direction, is an important second-order statistical quantity that characterizes the rate of occurrence of fade [18]. However, analytical solution for the LCR often depends mainly on the envelope distribution of the considered process.


## 2. Estimation Methods

### 2.1. Estimation for Known Power ParameterΩ

Various power estimation techniques have been implemented over years with some advantages and disadvantages. Accurate estimation of the average power of received fading signal is crucial for many reasons. Namely, power control and hand-off decisions in wireless communications are mainly based on the accurate estimation of the average signal power. Wireless communication link quality is also indicated through some system criterion measures, such as channel access, hand-off, and power control, that can be determined mainly based on the local mean signal levels. For example, the fading signal power could be estimated by using (38) from [26] without requiring any other parameter to be estimated.Letz denote the random envelope of received Nakagami-m faded signal, and let z˙ denote first derivative of z with respect to time. The joined probability density function (JPDF) of z and z˙ is then denoted with pzz˙(zz˙). The level crossing rate (LCR) at the envelope z is defined as the rate at which a fading signal envelope crosses level z in a positive or a negative direction and is analytically expressed by formula [18]:(1)Nzz=∫0∞z˙pz˙zz˙,zdz˙.JPDF can be expressed through the PDF of Nakagami-m faded signal envelope, pz(z), and conditioned PDF of envelope time derivative, pz˙∣z(z˙∣z), as pz˙z(z˙,z)=pz˙∣z(z˙∣z)×pz(z). PDF expression for Nakagami-m faded signal envelope, written in terms of two fading parameters, m and Ω, is given with(2)pzz=2mmz2m-1ΩmΓmexp⁡-mz2Ω,where Γ(a) denotes the Gamma function [27]. Common explanation of parameters m and Ω is that Ω defines average power of faded signal; that is, Ω=E(z2), while parameter m describes severity of fading process. As seen from [18] envelope time derivative, pz˙∣z(z˙∣z), is zero-mean Gaussian random process, defined with(3)pz˙∣zz˙∣z=12πσz˙exp⁡-z˙22σz˙2 and variance of σz˙2=π2fm2Ω/m. Here fm is denoted by maximal Doppler shift frequency.After substituting (2) and (3) into (1), LCR expression can be presented as(4)Nzz=mm-1/22πfmΩm-1/2Γmz2m-1exp⁡-mz2Ω. In order to determine value of parameter m as a function of received signal level and the number of times at which level is crossed, for known value of Ω, we only need to observe two signal levels, that is, z1 and z2 and rates at which Nakagami-m faded signal envelope has crossed those two levels, that is, Nz1 and Nz2 that can be determined as(5)Nz1=mm-1/22πfmΩm-1/2Γmz12m-1exp⁡-mz12Ω,Nz2=mm-1/22πfmΩm-1/2Γmz22m-1exp⁡-mz22Ω.Now, from (5) parameter m could be expressed in terms of z1, z2, Nz1, and Nz2 as(6)m=ln⁡Nz1/Nz2+ln⁡z1/z22ln⁡z1/z2+z22-z12/Ω.


### 2.2. Real Time Parameterm Estimation for Unknown Ω

In order to determine value of parameterm as a function of received signal level and the number of times at which level is crossed, for known value of Ω, we now need to observe three signal levels, that is, z1,z2, and z3 and rates at which Nakagami-m faded signal envelope has crossed those three levels, that is, Nz1, Nz2, and Nz3. For envelope level z3, LCR Nz3 can be expressed as follows:(7)Nz3=mm-1/22πfmΩm-1/2Γmz32m-1exp⁡-mz32Ω. Now, in similar manner parameter m could be expressed in terms of z1,z2, z3, Nz1, Nz2, and Nz3 as(8)m=121+z22-z12lnNz1/Nz3-z32-z12lnNz1/Nz2z22-z12lnz1/z3-z32-z12lnz1/z2while now(9)Ω=mz12-z22lnNz2/Nz1+2m-1lnz1/z2.


## 2.1. Estimation for Known Power ParameterΩ

Various power estimation techniques have been implemented over years with some advantages and disadvantages. Accurate estimation of the average power of received fading signal is crucial for many reasons. Namely, power control and hand-off decisions in wireless communications are mainly based on the accurate estimation of the average signal power. Wireless communication link quality is also indicated through some system criterion measures, such as channel access, hand-off, and power control, that can be determined mainly based on the local mean signal levels. For example, the fading signal power could be estimated by using (38) from [26] without requiring any other parameter to be estimated.Letz denote the random envelope of received Nakagami-m faded signal, and let z˙ denote first derivative of z with respect to time. The joined probability density function (JPDF) of z and z˙ is then denoted with pzz˙(zz˙). The level crossing rate (LCR) at the envelope z is defined as the rate at which a fading signal envelope crosses level z in a positive or a negative direction and is analytically expressed by formula [18]:(1)Nzz=∫0∞z˙pz˙zz˙,zdz˙.JPDF can be expressed through the PDF of Nakagami-m faded signal envelope, pz(z), and conditioned PDF of envelope time derivative, pz˙∣z(z˙∣z), as pz˙z(z˙,z)=pz˙∣z(z˙∣z)×pz(z). PDF expression for Nakagami-m faded signal envelope, written in terms of two fading parameters, m and Ω, is given with(2)pzz=2mmz2m-1ΩmΓmexp⁡-mz2Ω,where Γ(a) denotes the Gamma function [27]. Common explanation of parameters m and Ω is that Ω defines average power of faded signal; that is, Ω=E(z2), while parameter m describes severity of fading process. As seen from [18] envelope time derivative, pz˙∣z(z˙∣z), is zero-mean Gaussian random process, defined with(3)pz˙∣zz˙∣z=12πσz˙exp⁡-z˙22σz˙2 and variance of σz˙2=π2fm2Ω/m. Here fm is denoted by maximal Doppler shift frequency.After substituting (2) and (3) into (1), LCR expression can be presented as(4)Nzz=mm-1/22πfmΩm-1/2Γmz2m-1exp⁡-mz2Ω. In order to determine value of parameter m as a function of received signal level and the number of times at which level is crossed, for known value of Ω, we only need to observe two signal levels, that is, z1 and z2 and rates at which Nakagami-m faded signal envelope has crossed those two levels, that is, Nz1 and Nz2 that can be determined as(5)Nz1=mm-1/22πfmΩm-1/2Γmz12m-1exp⁡-mz12Ω,Nz2=mm-1/22πfmΩm-1/2Γmz22m-1exp⁡-mz22Ω.Now, from (5) parameter m could be expressed in terms of z1, z2, Nz1, and Nz2 as(6)m=ln⁡Nz1/Nz2+ln⁡z1/z22ln⁡z1/z2+z22-z12/Ω.


## 2.2. Real Time Parameterm Estimation for Unknown Ω

In order to determine value of parameterm as a function of received signal level and the number of times at which level is crossed, for known value of Ω, we now need to observe three signal levels, that is, z1,z2, and z3 and rates at which Nakagami-m faded signal envelope has crossed those three levels, that is, Nz1, Nz2, and Nz3. For envelope level z3, LCR Nz3 can be expressed as follows:(7)Nz3=mm-1/22πfmΩm-1/2Γmz32m-1exp⁡-mz32Ω. Now, in similar manner parameter m could be expressed in terms of z1,z2, z3, Nz1, Nz2, and Nz3 as(8)m=121+z22-z12lnNz1/Nz3-z32-z12lnNz1/Nz2z22-z12lnz1/z3-z32-z12lnz1/z2while now(9)Ω=mz12-z22lnNz2/Nz1+2m-1lnz1/z2.


## 3. Numerical Results

In Table1 comparison between parameter m value obtained by Matlab simulation over N=100000 samples of Nakagami-m vector z and value of m parameter estimated by using (6) from this paper has been shown. Nakagami-m random envelope vector z was obtained from Gaussian random processes by using its property that z=∑i=1m‍(Xi2+Yi2), where Xi and Yi are zero-mean normally distributed Gaussian random processes, N(0,σ2)  [18]. LCR levels z1 and z2 can be selected in a way that z1 is level with highest number of crosses, while level z2 is selected in the near environment of z1. Since it is known that envelope level z1 for which Nz1 has maximal value is probably the one for each stand z1=mean(z), level z1 can be selected in that way, while z2 level can be selected in the manner that z2=mean(z)-var(z). We will observe m parameter values in the range of m∈[1,3], since in [15] it has been shown that values from this range of m parameter correspond to the case of measurements performed at 60 GHz over NLOS channel conditions, when the user was mobile in a range of different indoor and outdoor small cell scenarios. After generating Nakagami-m fading processes, characterized with m∈[1,3], and defining LCR levels for each as explained, by using proposed method based on number of level crosses, we have obtained mestimated for each of the generated process. As it can be seen from Table 1, shown results provide very good match and proposed method has good accuracy. In Table 2 comparison between parameter m value obtained by Matlab simulation over N=1000000 samples of Nakagami-m vector z and value of m parameter estimated using (6) from this paper has been shown. Absolute error obtained for m parameter estimation for both cases has been presented in Figure 1.Table 1
Comparison between parameterm value obtained by Matlab simulation over N=100000 samples of Nakagami-m vector z and value of m parameter estimated by using (6) from this letter (z1=mean(z), z2=mean(z)-var(z), and Ω=E(z2)).
m
1
1.25
1.5
1.75
2
2.25
2.5
2.75
3
mestimated
1.071
1.2403
1.506
1.7507
2.009
2.2509
2.518
2.7712
3.0121Table 2
Comparison between parameterm value obtained by Matlab simulation over N=1000000 samples of Nakagami-m vector z and value of m parameter estimated by using (6) from this letter (z1=mean(z), z2=mean(z)-var(z), and Ω=E(z2)).
m
1
1.25
1.5
1.75
2
2.25
2.5
2.75
3
mestimated
1.038
1.2462
1.5042
1.7503
2.003
2.2509
2.512
2.7683
3.0117Figure 1
Absolute error of parameter estimation.By using proposed method Nakagamim parameter can be efficiently determined from various measurement, provided for 5G channel characterization, as ones in [15–17], instead of parameter m computation by using moment-based model or the ML estimation. By using proposed method when expressing m parameter in the function of Ricean K factor, as m=(1+K)2/1+2K, characterization of the behavior of LOS fading scenarios could be also performed.


## 4. Conclusion

Novel approach for Nakagamim parameter estimation based on LCR has been proposed in this paper. It has been shown that the proposed method not only avoids weaknesses of ML and moment method estimation approaches, but also is quite accurate, simple for realization, and can be used in real time. This interesting approach could also be used for fading parameters estimation of various NLOS channels in real time by using channel emulators for short, subwavelength distances in 5G multielement antenna realizations.



---

*Source: 5236246-2017-05-11.xml*