Estimation of the Reception Angle Distribution Based on the Power Delay Spectrum or Profile

The paper presents an estimation of the reception angle distribution based on temporal characteristics such as the power delay spectrum (PDS) or power delay profile (PDP). Here, we focus on such wireless environment, where the propagation phenomenon predominates in azimuth plane. As a basis to determine probability density function (PDF) of the angle of arrival (AOA), a geometrical channel model (GCM) in form of the multielliptical model for delayed scattering components and the von Mises’ PDF for local scattering components are used. Therefore, this estimator is called the distribution based on multielliptical model (DBMM). The parameters of GCM are defined on the basis of the PDS or PDP and the relative position of the transmitter and the receiver. In contrast to the previously known statistical models, DBMM ensures the estimation PDF of AOA by using the temporal characteristics of the channel for differing propagation conditions. Based on the results of measurements taken from the literature, DBMM verification, assessment of accuracy, and comparison with other models are shown. The results of comparison show that DBMM is the only model that provides the smallest least-squares error for different environments.


Introduction
In wireless communications, an angle distribution of the radio wave multipath components significantly affects the statistical properties of the received signals.Numerous measurements show that the environmental properties have a significant effect on the angle of arrival (AOA) of the receiving radio wave.Therefore, to perform statistical modeling of the impact of the propagation environment on deformations in time-spectrum of transmitted signals, knowledge of the probabilistic characteristics of angle is required.The search for new ways of transmitting information in wireless systems is primarily based on the results of simulation studies.Therefore, these characteristics are particularly important in these studies because their use provides the mapping of the effect of the object movement direction on the received signal.
One of the main ways of reception angle modeling is approximation of the measured data using standard probability density functions (PDFs).In this case, we talk about empirical models of angle distributions.In practice, Laplacian, Gaussian functions and Aulin's, Parson's models are the most commonly used distributions for azimuth [1] and elevation [2,3] angle, respectively.Unfortunately, there is no relationship that explicitly associates the parameters of these distributions with different types of propagation environments and the distance between the transmitter (Tx) and the receiver (Rx).Therefore, the use of geometric channel models (GCMs) gives the possibility of taking into account the impact of changes in the position of the objects (Tx, Rx) on the spatial properties of the received signals.These models reproduce the geometry of the spatial relationships among Tx, Rx, and the location scattering areas.The shape of the scattering areas and spatial density of the scatterers are the criteria that differentiate the individual models.In literature, models of PDF of AOA are based primarily on such scatterer regions as circle [4][5][6], ellipse [5][6][7], elliptical disc [8], hollowdisc [9,10], semispheroid [4,11], clipping semispheroid [12], ellipsoid [13], and bounded ellipsoid [14] for two-dimensional 2 International Journal of Antennas and Propagation (2D) and three-dimensional (3D) models, respectively.To the description of the spatial density of scatterers, such distributions are used as 2D uniform [4,5,9,15,16], 3D uniform [4,12], 2D Gaussian [17,18], 3D Gaussian [12], hyperbolic [19], Rayleigh and exponential [20], parabolic [6,14], "inverted parabolic" [21], and conical [22].GCMs provide the basis for theoretical analysis of PDFs of AOA.Each of these models is described by one, two, or more parameters.However, it is difficult to identify the physical premises that bind the values of these parameters with the propagation properties of environment.Therefore, there is no clear relationship between the PDF model parameters and the propagation conditions.For this reason, the verification of these PDFs in the literature is reduced to a measurement data approximation only.In this paper, the method of determining the structure of GCM is significantly different from existing solutions.In our proposition, the parameters of the temporal characteristics of channel are the basis for GCM.As a result, statistical model of AOA is PDF estimator closely related to the type of propagation environment.
In most land mobile access systems, the half power beamwidth (HPBW) of antenna pattern is 360 ∘ (omnidirectional antenna), or at least several tens of degrees in the azimuth plane.In elevation plane, HPBW has a few or at most a dozen degrees only.In this case, the phenomena modeling in the azimuth determines the mapping quality of the real signal.Therefore, these modeling conditions are decisive in simulation studies of wireless mobile channel.Unfortunately, none of the existing models does ensure the accuracy of the analytical description of the actual statistical properties of the AOA for different environmental conditions.Confirmed with the correctness of these observations are the results of a comparative analysis of the PDF of AOA for azimuth models that are presented in [23].This paper presents an evaluation of the accuracy of the mapping format of azimuth AOA by various 2D GCMs for different propagation scenarios.The results for selected measuring scenarios are reference data, whereas the accuracy measure of the models is the least-square error (LSE).The summary of [23] is that ". ..no one geometric model is best by all criteria and for all environments. ..." As an exception to this conclusion, we present an analytical model, estimator of the azimuth AOA distribution for multipath homogeneous environments.This PDF estimator is based on temporal characteristics of the channel such as the power delay spectrum (PDS) or power delay profile (PDP).For a given type of environment, the parameters of these characteristics for different positions of Tx/Rx are the estimation parameters.
The developed PDF estimator is based on the Parson-Bajwa multielliptical GCM [19].The geometry of this model results from the PDS or PDP.Literature analysis of the measurement results justifies such a methodology for determining PDF estimator.In [1], the correlation coefficient between delay and azimuth angle parameters is computed on the basis of the measurement data.Obtained linear regression line shows that these parameters are highly correlated [1, Figure 8].Therefore, the developed model ensures a clear reproduction of the angle statistical properties with respect to temporal characteristic that is associated with the type of environment.In simulation studies, this fact makes it possible to adopt the model to the specific environment that is described by temporal characteristics of radio channel.
The purpose of this study is to show that the developed PDF model of azimuth AOA for different environments provides the minimum value of the estimation error in comparison to the models previously presented in the literature.Our work concerns the uniform environments, which means that the comparative analysis is focused on the propagation scenarios with unimodal PDF of azimuth angle.
The paper is organized as follows.Section 2 presents the analysis principles of reception angle distribution for the developed model.Section 3 describes in detail AOA analysis for delayed components as well as methodology of von Mises' distribution application for local scattering components.Methodology of adaptation of model parameters is presented in Section 4. This adaptation is based on the channel characteristics and the parameters of measurement scenario.Based on the results of measurements taken from the literature, model verification, assessment of accuracy, and comparison with other PDFs of AOA models are shown in Section 5. Conclusion, which highlights the practical use of the model to predict the angle distribution of received signal at different environments, is presented in Section 6.

Distribution of Reception Angle: Principles of Analysis
As shown by numerous measurement data, for example, [24][25][26], due to multipath propagation of environment, a few to several replicas of the transmitted signal are received with various powers and time delays.This means that the received signal consists of the components.These components make up the time-clusters with a specific time delay,   ( = 1, . . .,  is number of clusters), relative to the direct path.The directivity, polarization, and gain from antenna patterns that are practical applications in wireless mobile communications are the reason for the dominance of the signal components, which are concentrated in the azimuth plane.It follows that the statistical properties of the received angle can be limited to the angle  in the azimuth plane.These comments indicate that the angular distribution of power can be represented as power azimuth spectrum (PAS), (), in the form where   () is PAS of the th time-cluster and  0 () is PAS for local scatterings, that is, for cluster with negligibly small delay relative to the direct path.However, the angle distribution and intensity of the direct path signal change with the change of Tx/Rx position.In this case, PDF of AOA, (), provides a basis for evaluation of the angular power dissipation and the Rice -factor describes the power distribution on direct path component and the components of local scatterings.The statistical description of PAS can be presented as [27]  () =  ⋅   () where  is the average power of the received signal,   is the power of all components delay by   (the th cluster),   () is PDF of AOA of multipath components delay by   ,  0 is the power of all components delay by zero,  0 () is PDF of AOA of multipath components from local scatterings (delay by zero), and (⋅) is the delta distribution.
Considering that  = ∫  − ()d = ∑  =0   , we have where   =   / is the normalized power of the th timecluster.Equation (3) shows that the statistical properties of the signal reception angle are determined by the components arriving with a delay to Rx and components from local scatterings.For delayed components, an analysis of reception angle is based on multielliptical channel model, whereas, for local scattering components, von Mises' PDF is used.Thus, the analytical form of PDF of AOA is called the distribution based on multielliptical model (DBMM).The basis for the DBMM parameters can be characteristics such as power delay profile (PDP) or power delay spectrum (PDS) that determine   values and the size of each ellipse.An additional input parameter is the height of the receiving antenna relative to the height of surrounding objects.This parameter is the basis for the assessment of local scattering intensity and determines the value of von Mises' PDF parameter.An analysis of the statistical properties of reception angle refers to homogeneous propagation environments.This means that such environments are dealt with for which the probability of a scatterer occurrence is the same for each direction when viewed from Tx.In practice, this means that the analysis covers those propagation scenarios that are characterized by unimodal PDF of the reception angle.In addition to the above considerations, to simplify the analyzed issues, the general assumptions that constitute the basis for most of GCMs are used [12,16,17]: (1) The radiation characteristics of the transmitting and receiving antennas are omnidirectional.
(2) Each propagation path from Tx to Rx consists of scatters on exactly one scattering element.
(3) For individual ellipses, each scatterer is a reradiating omnidirectional element with the same probability properties of scattering coefficient and uniform phase distribution.

Statistical Model of Angle of Arrival
The purpose of this work is to determine the analytical form of PDF of AOA of radio waves for a homogeneous multipath environment.The presented analysis uses the geometric mapping and statistical description of propagation phenomena.Propagation of delayed components is mapped by multielliptical GCM, while PDF of AOA of local scattering components is described by von Mises' PDF.In [27], a similar method of modeling the effects of propagation phenomena is used in analysis of the power azimuth spectrum (PAS).
The empirical studies demonstrate that, in time domain, the signal structure at the output of channels is shaped by a few or several clusters of the multipath components.These time-clusters are differentiated according to the values of power and time delay, after which the components arrive to Rx.This means that, in the azimuth plane, the location of scatterers defines an ellipse for components that make up a single time-cluster.This fact is the basis for the use of the Parson and Bajwa model in the statistical analysis of AOA.In this case, GCM creates a set of  ellipses that represent time-clusters with various   .The foci are common to all the ellipses.Their location determines the position of Tx and Rx, which are located at a distance .Signal components, for which the arrival delay is negligibly small and is expressed only in the phase difference relative to the direct component, form a particular time-cluster.These components are the result of local scattering that occurs in proximity to Tx/Rx antenna.The above model reproduces statistical location of the scatterers and thus provides the basis for determining PDF of AOA.The geometric structure of the considered problem and the adopted notations on the example of a chosen ellipse are shown in Figure 1.In addition, the local scattering area is highlighted in this figure.
The selected scatterer on a coordinate  is marked as .The angle of departure (AOD) from Tx and AOA at Rx are denoted by  and , respectively.For each ellipse, the major axis 2  depends on both  and   relative to direct path; namely, where  1 is the path length between Tx and ,  2 is the path length between  and Rx, and  is the speed of light.

Statistical Properties of AOA for Delay Components.
In this section, for individual PDF components, a brief analysis from [27] is presented.Based on the geometric properties of an ellipse and taking into account coordinate system as in Figure 1, it results that where   = 0.5/  < 1 is the eccentricity of the th ellipse.
The coordinate  can be expressed by  and  as follows: On the basis of relations ( 4)-( 6), we get By transforming the above expression, we can express cos  as By taking into account the omnidirectional pattern of the transmitter antenna, for each of the paths, all directions of departure are equally likely and do not depend on the environment properties.Thus, the equal probability can be assumed for each AOD.It means that PDF of AOD can be expressed as   () = (2) −1 for  ∈ ⟨−, ⟩.Therefore, density function   (cos ) will take the following form: Let us note that Therefore, for paths arriving at Rx with   , PDF of  is but so, we have The definition of the ellipse eccentric indicates that   () depends on  and   .The first parameter describes the geometry of propagation scenario, and the second is clear from propagation properties of the environment.

Statistical Properties of AOA for Local Scattering Components.
Uniform distribution is one of the basic models describing statistical properties of AOA for the local scattering components.This distribution is used to describe  0 () in Clarke's 2D isotropic scattering model [20].Examples of other distributions related to the local scatterings in the surroundings of the receiving antenna are presented in [15,18].As shown in [18], the von Mises PDF expressed by where  0 (⋅) is the zero-order modified Bessel function and  is a constant that controls the width of AOA of the scattering components, is a generalized distribution.
For  = 0, we obtain   () = (2) −1 while  → ∞ yields  0 () → ().This means that, with the increase in the scattering intensity of the propagation paths, the parameter value decreases.For small , this function approximates the cardioid while for large values it resembles the Gaussian PDF.In practice,  depends on the height and the density of objects in the surrounding of the receiving antenna.In [18], the presented measurement results show that if the antenna height is less than the height of the surrounding objects, then  ≤ 3.3.

Analytical PDF of AOA Model.
The dependences (3), (14), and ( 15) are the basis for the determination of the analytical form of PDF of AOA The above relation shows that   () is significantly dependent on   ,   , , and .This means that the analyzed function depends on the distance between Tx and Rx, the type of the environment, and the nature of the local scatterings.This means that   () changes with change of the position of Tx/Rx.Thus, for a given type of propagation environment, (16) describes PDF estimator of AOA, whose parameters are determined on the basis of changing temporal characteristics of channel.Also note that ( 16) is theoretical due to the presence of -distribution that describes PDF of AOA of the direct path.In practice, obtaining a result in the form of the distribution is not impossible.The reason for this is the limited measurement discreteness of angle that results from the limited width of the receiving antenna beam.Therefore, to evaluate the model approximation error with respect to measurement results, the distribution component should be replaced by a function that belongs to a set of delta-forming functions.These functions take nonzero values only for the interval that corresponds to the measurement discreteness of the angle, and their surface areas are equal to the normalized received signal power.

Example of Numerical Results
In this section, the methodology for determining PDF of AOA is described for selected test scenarios that the detailed characterization is presented in [1].For DBMM, the input data are identified on the basis of the measurement of PDS that is shown in the graph in this paper ([1, Figure 6]).Number of the time-clusters and the values of their parameters, that is, ,   , and   , are defined on the base of PDS selected points.These points are local extremes of the function, which is the difference between PDS and its trend.For the analyzed scenario, the function of trend is Laplacian function.The choice of points that identify the parameters of DBMM is shown in Figure 2, whereas the determined values are contained in Table 1.
The distance  = 1500 m between Tx and Rx is adopted based on the description of the scenario and analysis of measurement results presented in [4].To assure consistence in extracting numerical data from graphs the software WebPlotDigitizer [31] is used.The dataset, contained in Table 1, defined a first group of components that are presented in expression (16).This group describes the angular power dispersion of the signal components that arrive at the receiver with delays.Parameters of remaining components of   (), that is,  and , are determined by the propagation environment nature and parameters of the measuring system.
In [1], the described measurement campaign was conducted in urban area characterized by buildings ranging from four to six floors and placed on a slightly rolling terrain.The transmitting equipment is mounted in car.The transmitting antenna is omnidirectional whereas Rx is equipped with an eight-element uniform linear antenna array whose height corresponds to the average rooftop level, that is, 21 m above ground level.Thus, for the analyzed scenario, only the objects in immediate surroundings of Tx and the surfaces of roofs at the location of the receiving antenna are a source of local scatterings.This means that the local scatterings are targeted.In this case, it is the basis for the adoption of high-value of where   and   values are contained in Table 1.
In Figure 3, comparison of the empirical data with the theoretical analysis results that are obtained on the base of ( 17) is presented.
In subsequent sections,   is used to assess the accuracy of mapping the actual characteristics by the selected analytical models of PDF of AOA.

Verification and Comparison of Models
The main aim of this section is to verify the accuracy and compare DBMM with other statistical AOA models for which comprehensive comparative analysis is presented in [23].According to the accepted assumptions in Section 2, comparative analysis refers to selected GCMs and measurement scenarios for which PDF of AOA is unimodal.On the base of the analysis results that are shown in [23] (Table 4), the comparison is reduced to a set of models that are best matched to the empirical data.Table 2 shows the selected test scenarios and corresponding models that provide the best fit to the measurement results, whereas Table 3 contains the analytical dependences that define PDFs of AOA for different GCMs based on [23].Additionally, in accordance with the suggestion made in the conclusion [23], the Gaussian model is included in comparison analysis.
In [23], LSE is used as an objective function with the minimizing criterion for determining the parameters of the analyzed models: The parameters of DBMM are determined otherwise.The bases for the determination of ,   , and   are the temporal channel characteristics such as PDS or PDP.For each scenario, the empirical data is not a uniform set.For [1,28,33] scenarios, PDS is the basis for DBMM parameters, so all of them take statistically constant values.PDP is used with respect to [26,29,30,34] scenarios.In these cases, (16) represents PDF estimator of AOA that is associated with the specific geometry of environment.The geometry of Rx antenna immediate surroundings determines the values of the other two parameters (, ).Due to the lack of data, the ranges of these parameters are adopted on the basis of the description of the measurement conditions.Similarly, as in [23], the criterion of minimizing LSE is used to select  and  values.In Section 4, the methodology for determining the parameters of DBMM and then PDF of AOA is shown in the example of results obtained from the Stockholm  measurement scenario [1].In the same way, this methodology is used for the second measurement scenario from [1] (Aarhus) and scenarios from [29,34].For measurement scenarios from Stockholm and Aarhus, different PDSs and the amount of the receiving antennas with respect to the environment objects are the causes of various parameters for DBMM.In measurement campaign that was conducted in Aarhus, the receiving antenna was mounted at higher altitude, that is, 12 m, than average rooftop level of the surrounding buildings.Thus, the scattering that occurs in the vicinity of the transmitting antenna especially determines the local scattering.This is the reason for reducing the scope of these scatterings, which is a prerequisite for increasing  from 60 for Stockholm to 120 for Aarhus.
Because of the nature of the propagation phenomena that are dominant in [26] measurement scenario, the methodology for determining PDF of AOA, on the basis of DBMM, is slightly modified.In this case, the amount of transmitting and receiving antennas was 2.7 m and 4.4 m, respectively.Mutual location of Tx and Rx and the average height of the buildings, that is, 8 m, to a significant extent determine the nature of the propagation environment.For this measurement scenario, an environment has the features of the canyon.This means that the propagation paths that run along the canyon have a dominant influence on the received signal.The geometry of the environment and the location of Tx show that the direct paths to the scattering places, located within the canyon, only concern the first three ellipses.Thus, for  1 = 0.02 s,  2 = 0.06 s, and  3 = 0.08 s,   () ( = 1, 2, 3) are different from zero for  ∈ ⟨−2.25 ∘ , 2.25 ∘ ⟩.Relative to the position of Rx, this range represents the canyon angular width that determines the location of the scatterers.Geometrically, these places are represented by three subsequent sections of ellipses that are presented in the surroundings of Tx.
As presented in [30], the measurement scenario does not contain characteristics that describe the propagation properties of the environment.This makes it difficult to define the parameters of DBMM for the analyzed scenario.However, measurement campaigns that are described in [29,30] relate to the same environment in Aalberg.Therefore, for these scenarios, the same parameter values of other models are adopted in [23].Also, in [28,33], the propagation properties of environments are not described.For these scenarios, DBMM parameters are defined on the basis of descriptions of the measurement conditions.They show that the measurement was carried out on frequencies of the standard GSM 1800 in rural and suburban environments for [28,33], respectively.These facts are a premise for the use of parameter values that describe the rural environment in accordance with COST 207 model [35].The descriptions of these scenarios also show that PAS measurements were conducted in line-of-sight (LOS) propagation conditions.The presence of the direct path determines AOA that theoretically is described by  0 [/( + 1)]().In this case, comparative analysis of analytical dependence (16) and empirical results is difficult due to () presence.In practice, the minimum beam width of the measurement antenna pattern is the basis to replace () by a triangular function.This function belongs to the set of -generating functions and maps the signal power arriving from the direction of the direct path.Hence, the span angle, Δ, for which the triangular function takes values different from zero, is associated with angular measurement discreteness that determines the beam width of the antenna pattern.Owing to the lack of accurate data, for measurement scenario [33], Δ = 2.5 ∘ is adopted.In [28], the presented graphs are averaged measurement results that have been obtained over the entire measurement route.Therefore, in this case, Δ is equal to the angular width of measurement route that is determined relative to the location of the receiving antenna.On the basis of the simplified scenario map, Δ = 15 ∘ is adopted.In both cases, the surface area of triangular functions depends on the level of direct path signal power and amounts  0 /( + 1).
Distance between Tx and Rx is one of the important parameters that determine the geometry of the models.For [26,34] scenarios,  values are strictly defined, whereas, for other measurement campaigns, they are adopted as the average Tx distance relative to Rx location.The set of all DBMM parameters that are determined on the basis of descriptions of analyzed measurement scenarios is contained    [23], the values of parameters that minimize LSE for selected models are presented in Table 5.
Assessment of the accuracy of individual models is based on an assessment of   (18) with respect to the measurement results for different environments.In Table 6, the results of numerical calculations are presented.
Additionally, the results obtained by Wong et al. [23] and the results for DBMM are shown in Table 7. Graphical illustrations of PDF of AOA models and empirical data for the analyzed scenarios are presented in Figures 4-11.
For different environments, the analysis of   and LSE show that, among the analytical models of PDF of AOA, DBMM is the model that accurately reproduces the results of measurements for all analyzed scenarios.The characteristics of the received signals and environmental conditions are the basis for DBMM parameters listed in Table 4.As a result, DBMM is closely related to the properties of each environment, which results in minimizing   and LSE in comparison with other models.Seemingly high degree of freedom that is in DBMM results from a larger number   of the parameters.However, these parameters are closely determined based on the transmission conditions and the characteristics of the received signals.
For analytical models of PDF of AOA that are presented in [23], the parameters are practically immeasurable.This fact significantly impedes their use to predict the angle

Figure 1 :
Figure 1: Geometry of the analyzed problem.

Figure 2 :
Figure 2: Identification of parameters for individual clusters on the example of measurement scenario from [1] (Stockholm).

Figure 4 :
Figure 4: Fitting for various geometric models to the empirical data in Matthews et al. [34, Figure 7].

Figure 5 :
Figure 5: Fitting for various geometric models to the empirical data in Pedersen et al. [33, Figure 1].

Figure 6 :
Figure 6: Fitting for various geometric models to the empirical data in Kuchar et al. [28, Figure 3].

Figure 7 :
Figure 7: Fitting for various geometric models to the empirical data in Takada et al. [26, Figure 4].

Figure 8 :
Figure 8: Fitting for various geometric models to the empirical data in Fleury et al. [29, Figure 16].

Figure 9 :
Figure 9: Fitting for various geometric models to the empirical data in Mogensen et al. [30, Figure 3].

Table 2 :
Selected test scenarios and the best matched models.

Table 3 :
PDFs of AOA for selected GCMs.

Table 4 :
Parameters of DBMM for analyzed measurement scenarios.

Table 4 .
Based on

Table 6 :
Kolmogorov-Smirnov statistic for selected models calibrated by empirical data.

Table 7 :
LSE for selected models calibrated by empirical data.