Novel Method for 5 G Systems NLOS Channels Parameter Estimation

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 estimationwill be presented. Estimationwill 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.


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 lineof-sight (LOS) scenarios [1][2][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, laterarriving 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 (Peerto-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- distribution.Nakagami- 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][10][11][12][13][14].In [15][16][17] for the purpose of modelling observed 5G system propagation properties the Nakagami  parameter is directly computed from the measured data.As a general fading distribution, Nakagami- 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- fading model can also be transformed into Ricean fading model, by expressing  parameter in the function of Ricean  factor, as  = (1 + ) 2 /(1 + 2) [18].In that manner characterizing of the behavior of LOS fading envelope could be also performed.The Nakagami- fading exploits Nakagami probability density function (PDF) for the random envelope of received signal which is 2 International Journal of Antennas and Propagation written in the function of two parameters: scale parameter and shape parameter, called the fading severity parameter or -parameter.Determining  is a problem in Nakagami PDF estimation.For observed set of empirical fading signal data, the value of distribution parameter  should be estimated from it, in order to use the Nakagami- distribution to model a given set of values.Acquaintance of the  parameter is mandatory for the optimal reception of signals in Nakagami- fading environment [19], and in such case parameter  should be determined very accurately and very fast.Acquaintance of the  parameter is also necessary in the transmitter adaptation process, where its value could be feedback with respect to the channel information.Two most wellknown procedures used for the estimation of the Nakagami- fading parameter, , are (1) maximum likelihood (ML) estimation and (2) moment-based estimation.ML-based  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  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  parameter is the moment-based estimation.Observing the second and the fourth Nakagami- sample moments, estimation of parameter  was carried out in [22].An improvement of proposed method can be found in [23] where parameter  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- 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.

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 handoff 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.
Let  denote the random envelope of received Nakagami- faded signal, and let ż denote first derivative of  with respect to time.The joined probability density function (JPDF) of  and ż is then denoted with   ż( ż ).The level crossing rate (LCR) at the envelope  is defined as the rate at which a fading signal envelope crosses level  in a positive or a negative direction and is analytically expressed by formula [18]: JPDF can be expressed through the PDF of Nakagami- faded signal envelope,   (), and conditioned PDF of envelope time derivative,  ż | ( ż | ), as  ż  ( ż , ) =  ż | ( ż | ) ×   ().PDF expression for Nakagami- faded signal envelope, written in terms of two fading parameters,  and Ω, is given with where Γ() denotes the Gamma function [27].Common explanation of parameters  and Ω is that Ω defines average power of faded signal; that is, Ω = ( 2 ), while parameter  describes severity of fading process.As seen from [18] envelope time derivative,  ż | ( ż | ), is zero-mean Gaussian random process, defined with and variance of  2 ż =  2  2  Ω/.Here   is denoted by maximal Doppler shift frequency.
After substituting (2) and ( 3) into (1), LCR expression can be presented as In order to determine value of parameter  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,  1 and  2 and rates at which Nakagami- faded signal envelope has crossed those two levels, that is,   1 and   2 that can be determined as ) . ( International Journal of Antennas and Propagation 3 2.2.Real Time Parameter  Estimation for Unknown Ω.In order to determine value of parameter  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,  1 ,  2 , and  3 and rates at which Nakagami- faded signal envelope has crossed those three levels, that is,   1 ,   2 , and   3 .For envelope level  3 , LCR   3 can be expressed as follows: Now, in similar manner parameter  could be expressed in terms of  1 ,  2 ,  3 ,   1 ,   2 , and   3 as while now

Numerical Results
In Table 1 comparison between parameter  value obtained by Matlab simulation over  = 100000 samples of Nakagami- vector  and value of  parameter estimated by using (6) from this paper has been shown.Nakagami- random envelope vector  was obtained from Gaussian random processes by using its property that  = √∑  =1 ( 2  +  2  ), where   and   are zero-mean normally distributed Gaussian random processes, (0,  2 ) [18].LCR levels  1 and  2 can be selected in a way that  1 is level with highest number of crosses, while level  2 is selected in the near environment of  1 .Since it is known that envelope level  1 for which   1 has maximal value is probably the one for each stand  1 = mean(), level  1 can be selected in that way, while  2 level can be selected in the manner that  2 = mean() − var().We will observe  parameter values in the range of  ∈ [1,3], since in [15] it has been shown that values from this range of  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- fading processes, characterized with  ∈ [1,3], and defining LCR levels for each as explained, by using proposed method based on number of level crosses, we have obtained  estimated 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  value obtained by Matlab simulation over  = 1000000 samples of Nakagami- vector  and value of  parameter estimated using (6) from this paper has been shown.Absolute error obtained for  parameter estimation for both cases has been presented in Figure 1.
By using proposed method Nakagami  parameter can be efficiently determined from various measurement, provided for 5G channel characterization, as ones in [15][16][17], instead of parameter  computation by using moment-based model or the ML estimation.By using proposed method when expressing  parameter in the function of Ricean  factor, as  = (1 + ) 2 /(1 + 2), characterization of the behavior of LOS fading scenarios could be also performed.

Conclusion
Novel approach for Nakagami  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.

Figure 1 :
Figure 1: Absolute error of parameter estimation.