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

Authors: Cezary Ziółkowski; Jan M. Kelner
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

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.


---

## Body

## 1. 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–6], ellipse [5–7], elliptical disc [8], hollow-disc [9, 10], semispheroid [4, 11], clipping semispheroid [12], ellipsoid [13], and bounded ellipsoid [14] for two-dimensional (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 is360° (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. Section2 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.


## 2. Distribution of Reception Angle: Principles of Analysis

As shown by numerous measurement data, for example, [24–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, τi (i=1,…,N 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), Pβ, in the form (1)Pβ=∑i=1NPiβ+P0β,where Piβ is PAS of the ith time-cluster and P0β 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,fβ, 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](2)Pβ=P·fββ=∑i=1NPifβiβ+1κ+1P0f0β+κκ+1P0δβ,where P is the average power of the received signal, Pi is the power of all components delay by τi (the ith cluster), fβiβ is PDF of AOA of multipath components delay by τi, P0 is the power of all components delay by zero, f0β is PDF of AOA of multipath components from local scatterings (delay by zero), and δ· is the delta distribution.Considering thatP=∫-ππPβdβ=∑i=0NPi, we have(3)fββ=∑i=1Npifβiβ+1κ+1p0f0β+κκ+1p0δβ,where pi=Pi/P is the normalized power of the ith time-cluster.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 pi 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.


## 3. 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 ofN ellipses that represent time-clusters with various τi. The foci are common to all the ellipses. Their location determines the position of Tx and Rx, which are located at a distance D. 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.Figure 1
Geometry of the analyzed problem.The selected scatterer on a coordinatex is marked as S. The angle of departure (AOD) from Tx and AOA at Rx are denoted by α and β, respectively. For each ellipse, the major axis 2ai depends on both D and τi relative to direct path; namely, (4)2ai=r1i+r2i=D+cτi,where r1i is the path length between Tx and S, r2i is the path length between S and Rx, and c is the speed of light.


### 3.1. 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 (5)r1i=ai-eix,r2i=ai+eix,where ei=0.5D/ai<1 is the eccentricity of the ith ellipse. The coordinate x can be expressed by α and β as follows: (6)x=r1icos⁡α+0.5D,x=r2icos⁡β-0.5D.On the basis of relations (4)–(6), we get(7)11+eicos⁡α+11-eicos⁡β=2aiai-0.5D.By transforming the above expression, we can express cos⁡α as(8)cos⁡α=aicos⁡β+0.5Deicos⁡β-Dai+0.5Dei-Dcos⁡β.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 asfαα=2π-1 for α∈-π,π. Therefore, density function fcαcos⁡α will take the following form: (9)fcαcos⁡α=fαcos⁡αdαdcos⁡α=12π11-cos2⁡α.Let us note that(10)dcos⁡αdcos⁡β=ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2=ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2.Thus, PDF of cos⁡β can be expressed as (11)fcβcos⁡β=fcαcos⁡βdcos⁡αdcos⁡β=12π11-aicos⁡β+0.5Deicos⁡β-D/ai+0.5Dei-Dcos⁡β2ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2=12π11-cos2⁡βai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β.Therefore, for paths arriving at Rx with τi, PDF of β is (12)fiββ=fcβcos⁡βdcos⁡βdβ=12πai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β,but (13)ai+0.5Dei2-D2=ai1+ei22-2ei2=ai1+ei2+2ei1+ei2-2ei=ai1+ei21-ei2=ai1+ei1-ei=ai1-ei2,so, we have (14)fiββ=12π1-ei21+ei2-2eicos⁡β,β∈-π,π.The definition of the ellipse eccentric indicates thatfiββ depends on D and τi. The first parameter describes the geometry of propagation scenario, and the second is clear from propagation properties of the environment.


### 3.2. 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 describef0β 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(15)f0β=exp⁡μcos⁡β2πI0μ,β∈-π,π,μ≥0,where I0· 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 fiββ=2π-1 while μ→∞ yields f0β→δβ. 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.


### 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 (16)fββ=12π∑i=1Npi1-ei21+ei2-2eicos⁡β+p01κ+1exp⁡μcos⁡β2πI0μ+p0κκ+1δβ.The above relation shows thatfββ is significantly dependent on ei, pi, κ, 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 fββ 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.


## 3.1. 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 (5)r1i=ai-eix,r2i=ai+eix,where ei=0.5D/ai<1 is the eccentricity of the ith ellipse. The coordinate x can be expressed by α and β as follows: (6)x=r1icos⁡α+0.5D,x=r2icos⁡β-0.5D.On the basis of relations (4)–(6), we get(7)11+eicos⁡α+11-eicos⁡β=2aiai-0.5D.By transforming the above expression, we can express cos⁡α as(8)cos⁡α=aicos⁡β+0.5Deicos⁡β-Dai+0.5Dei-Dcos⁡β.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 asfαα=2π-1 for α∈-π,π. Therefore, density function fcαcos⁡α will take the following form: (9)fcαcos⁡α=fαcos⁡αdαdcos⁡α=12π11-cos2⁡α.Let us note that(10)dcos⁡αdcos⁡β=ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2=ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2.Thus, PDF of cos⁡β can be expressed as (11)fcβcos⁡β=fcαcos⁡βdcos⁡αdcos⁡β=12π11-aicos⁡β+0.5Deicos⁡β-D/ai+0.5Dei-Dcos⁡β2ai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β2=12π11-cos2⁡βai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β.Therefore, for paths arriving at Rx with τi, PDF of β is (12)fiββ=fcβcos⁡βdcos⁡βdβ=12πai+0.5Dei2-D2ai+0.5Dei-Dcos⁡β,but (13)ai+0.5Dei2-D2=ai1+ei22-2ei2=ai1+ei2+2ei1+ei2-2ei=ai1+ei21-ei2=ai1+ei1-ei=ai1-ei2,so, we have (14)fiββ=12π1-ei21+ei2-2eicos⁡β,β∈-π,π.The definition of the ellipse eccentric indicates thatfiββ depends on D and τi. The first parameter describes the geometry of propagation scenario, and the second is clear from propagation properties of the environment.


## 3.2. 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 describef0β 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(15)f0β=exp⁡μcos⁡β2πI0μ,β∈-π,π,μ≥0,where I0· 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 fiββ=2π-1 while μ→∞ yields f0β→δβ. 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.


## 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 (16)fββ=12π∑i=1Npi1-ei21+ei2-2eicos⁡β+p01κ+1exp⁡μcos⁡β2πI0μ+p0κκ+1δβ.The above relation shows thatfββ is significantly dependent on ei, pi, κ, 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 fββ 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.


## 4. 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, N, τi, and pi, 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.Table 1
Parameters for individual clusters (N=7).
Cluster numberi
P
i (dB)
P
i
p
i
τ
i (µs)
a
i (km)
e
i
0
0.00
1.0000
0.4003
0.00
—
—
1
−1.45
0.7161
0.2867
0.28
0.7920
0.9470
2
−4.91
0.3228
0.1292
0.78
0.8670
0.8651
3
−6.03
0.2495
0.0999
1.12
0.9180
0.8170
4
−10.64
0.0863
0.0345
1.71
1.0065
0.7452
5
−12.31
0.0587
0.0235
2.06
1.0590
0.7082
6
−13.72
0.0425
0.0170
2.43
1.1145
0.6729
7
−16.54
0.0222
0.0089
3.04
1.2060
0.6219Figure 2
Identification of parameters for individual clusters on the example of measurement scenario from [1] (Stockholm).The distanceD=1500m 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 softwareWebPlotDigitizer [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 fββ, 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 μ equal to 60. Because the measurement was carried out with non-line-of-sight (NLOS) between Tx and Rx, we have κ=0. Taking into account the above considerations, fββ reduces to the form (17)fββ=12π∑i=1Lpi1-ei21+ei2-2eicos⁡β+p0exp⁡μcos⁡β2πI0μ,where ei and pi values are contained in Table 1.In Figure3, comparison of the empirical data with the theoretical analysis results that are obtained on the base of (17) is presented.Figure 3
Comparing theoretical analysis to the empirical results in Pedersen et al. [1, Figure 4] (Stockholm).As a measure of the accuracy of measurement data by DBMM, the Kolmogorov-Smirnov statistic,Dm, is used [32]:(18)Dm=supβk⁡Fmβk-Fββk-β0,where Fββk=∫-πβkfββdβ is the cumulative distribution function (CDF), fββk represents PDF of AOA for analyzed model, Fmβk=(1/2)∑l=2kfmβl-1+fmβlβl-βl-1 and fmβk(k=1,…,K) denote the normalized CDF and PDF from empirical dataset, respectively, K refers to cardinality of the set of measurement data, and β0 is a parameter that calibrates PDF of AOA model to the direction of Tx-Rx for empirical data. For analyzed scenario, we have Dm=0.0264.In subsequent sections,Dm is used to assess the accuracy of mapping the actual characteristics by the selected analytical models of PDF of AOA.


## 5. 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.Table 2
Selected test scenarios and the best matched models.
Measurement scenario (empirical PDF of AOA)
PDF of AOA models
Matthews et al. [34, Figure 7]
(Wong et al. [23, Figure 4])
Gaussian
Rayleigh circular (Rx outside)
Pedersen et al. [33, Figure 1]
(Wong et al. [23, Figure 5])
Uniform elliptical (Rx outside)
Kuchar et al. [28, Figure 3]
(Wong et al. [23, Figure 6])
Gaussian
Rayleigh circular (Rx outside)
Takada et al. [26, Figure 4]
(Wong et al. [23, Figure 7])
Uniform elliptical (Rx outside)
Fleury et al. [29, Figure 16]
(Wong et al. [23, Figure 8])
Uniform elliptical (Rx outside)
Gaussian
Rayleigh circular (Rx outside)
Mogensen et al. [30, Figure 3]
(Wong et al. [23, Figure 9])
Uniform elliptical (Rx outside)
Pedersen et al. [1, Figure 4], Aarhus
(Wong et al. [23, Figure 10])
Uniform elliptical (Rx outside)
Pedersen et al. [1, Figure 4], Stockholm
(Wong et al. [23, Figure 11])
Uniform elliptical (Rx outside)
Uniform elliptical (Rx inside)Table 3
PDFs of AOA for selected GCMs.
Model
PDF of AOA,fββ (where fββ≠0)
Conditions
Uniform elliptical(Rx outside)[38]
2
π
b
a
3
D
a
2
cos
⁡
β
cos
2
⁡
β
-
cos
2
⁡
β
+
a
/
b
2
sin
2
⁡
β
1
-
D
/
a
2
1
/
2
b
/
a
2
cos
2
⁡
β
+
sin
2
⁡
β
2
β
∈
-
β
m
a
x
,
β
m
a
x, βmax=arctanbatan⁡arcsin⁡aD,2b<2a<D,2a, 2b: major and minor axis of elliptical disc for uniform density (centered at Tx but excluding Rx)
Uniform elliptical(Rx inside)[16, 17, 39]
1
2
π
L
/
D
L
/
D
2
-
1
L
/
D
2
-
1
L
/
D
-
cos
⁡
β
2
β
∈
-
π
,
π,L: major axis of elliptical disc for uniform density (support region focused on Tx and Rx)
Gaussian[4, 40, 41]
1
2
π
exp
⁡
-
1
2
D
σ
2
+
cos
⁡
β
2
2
π
D
σ
exp
⁡
cos
2
⁡
β
-
1
2
D
σ
2
erfc
⁡
-
cos
⁡
β
2
D
σ
β
∈
-
π
,
π,σ2: variance of Gaussian intensity (centered at Tx)
Rayleigh circular(Rx outside)[42]
1
2
π
exp
⁡
-
1
2
D
R
2
1
+
erf
⁡
cos
⁡
β
2
D
R
·
1
+
π
2
D
R
cos
⁡
β
exp
⁡
1
2
D
R
cos
⁡
β
2
β
∈
-
π
,
π,R≪D,R: parameter of Rayleigh intensity (centered at Tx)
D: distance between Tx and Rx.Table 4
Parameters of DBMM for analyzed measurement scenarios.
Measurement scenario
[34]
[33]
[28]
[26]
[29]
[30]
[1]Aarhus
[1]Stockholm
Location of measurement
Leeds, England
Bristol, England
Stuttgart, Germany
Yokosuka, Japan
Aalborg, Denmark
Aalborg, Denmark
Aarhus, Denmark
Stockholm, Sweden
D (m)
1800
5000
400
219
2100
2100
1500
1500
µ
60
240
90
200
140
140
120
60
κ
0
1
0.3
0
0
0
0
0
Δ
ϕ (°)
—
2.5
15
—
—
—
—
—
β
0 (°)
23.4
0
3.5
2.96
24.86
0.8
0.26
0.25
PDP/PDSbased on
[34, Figure 10]
COST 207 [35]rural area
COST 207 [35]rural area [modified]
[26, Figure 5]
[29, Figure 16]
[29, Figure 16]
[1, Figure 6]
[1, Figure 6]
N
6
5
5
7
4
4
5
7
p
0
0.8506
0.6024
0.8584
0.2115
0.4735
0.4735
0.4363
0.4003
p
i
0.1060, 0.0126, 0.0188, 0.0083, 0.0021, 0.0016.
0.2410, 0.0964, 0.0361, 0.0181, 0.0060.
0.0858, 0.0343, 0.0129, 0.0064, 0.0022.
0.2061, 0.0814, 0.3178, 0.1192, 0.0158, 0.0436, 0.0046.
0.4442, 0.0733, 0.0049, 0.0041.
0.4442, 0.0733, 0.0049, 0.0041.
0.3019, 0.2091, 0.0436, 0.0074, 0.0017.
0.2867, 0.1292, 0.0999, 0.0345, 0.0235, 0.0170, 0.0089.
τ
i (µs)
01.22, 01.86, 02.47, 06.52, 07.73, 10.84.
0.10, 0.20, 0.30, 0.40, 0.50.
0.10, 0.20, 0.30, 0.40, 0.50.
0.02, 0.06, 0.08, 0.14, 0.19, 0.32, 0.48.
0.95, 3.76, 5.61, 8.36.
0.95, 3.76, 5.61, 8.36.
0.14, 0.32, 0.83, 1.28, 1.95.
0.28, 0.78, 1.12, 1.71, 2.06, 2.43, 3.04.In [23], LSE is used as an objective function with the minimizing criterion for determining the parameters of the analyzed models:(19)LSE=1K∑k=1Kfmβk-fββk-β02.The parameters of DBMM are determined otherwise. The bases for the determination ofN, τi, and pi 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, fiββ(i=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 p0κ/κ+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 p0κ/κ+1.Distance between Tx and Rx is one of the important parameters that determine the geometry of the models. For [26, 34] scenarios, D 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 in Table 4. Based on [23], the values of parameters that minimize LSE for selected models are presented in Table 5.Table 5
Parameters of selected models that minimize LSE [23].
Measurement scenario
[34]
[33]
[28]
[26]
[29]
[30]
[1]Aarhus
[1]Stockholm
Uniform elliptical(Rx outside)
b/a
0.27
0.03
0.21
0.05
0.17
0.16
0.11
0.20
a/D
0.81
0.99
0.81
0.96
0.99
0.99
0.99
0.99
Uniform elliptical(Rx inside)
L/D
1.0204
1.0101
1.0101
1.0101
1.0101
1.0101
1.0101
1.0204
Gaussian
σ/D
0.14
0.02
0.11
0.03
0.10
0.10
0.07
0.13
Rayleigh circular(Rx outside)
R/D
0.139
0.027
0.110
0.032
0.104
0.098
0.073
0.134Assessment of the accuracy of individual models is based on an assessment ofDm (18) with respect to the measurement results for different environments. In Table 6, the results of numerical calculations are presented.Table 6
Kolmogorov-Smirnov statistic for selected models calibrated by empirical data.
Measurement scenario
[34]
[33]
[28]
[26]
[29]
[30]
[1]Aarhus
[1]Stockholm
Uniform elliptical (Rx outside)
0.0238
0.0654
0.0129
0.0733
0.0343
0.1001
0.0558
0.0309
Uniform elliptical (Rx inside)
0.0306
0.1890
0.0384
0.1584
0.0636
0.1148
0.0883
0.0314
Gaussian
0.0256
0.0909
0.0107
0.1017
0.0554
0.1004
0.0660
0.0429
Rayleigh circular (Rx outside)
0.0239
0.0684
0.0107
0.0937
0.0492
0.1029
0.0630
0.0395
DBMM
0.0196
0.0376
0.0083
0.0548
0.0271
0.0868
0.0473
0.0264Additionally, 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.Table 7
LSE for selected models calibrated by empirical data.
Measurement scenario
[34]
[33]
[28]
[26]
[29]
[30]
[1]Aarhus
[1]Stockholm
Uniform elliptical (Rx outside)
0.0054104
1.6483
0.0069964
0.12699
0.020099
0.014380
0.084363
0.043075
Uniform elliptical (Rx inside)
0.0098299
39.6960
0.0322370
1.15250
0.027642
0.018240
0.196290
0.040088
Gaussian
0.0041719
3.2328
0.0060807
0.14807
0.019722
0.019087
0.191920
0.083137
Rayleigh circular (Rx outside)
0.0041526
6.3259
0.0060807
0.18178
0.018841
0.018867
0.187050
0.081898
DBMM
0.0030141
1.0963
0.0054167
0.11389
0.012770
0.010482
0.050515
0.026164Figure 4
Fitting for various geometric models to the empirical data in Matthews et al. [34, Figure 7].Figure 5
Fitting for various geometric models to the empirical data in Pedersen et al. [33, Figure 1].Figure 6
Fitting for various geometric models to the empirical data in Kuchar et al. [28, Figure 3].Figure 7
Fitting for various geometric models to the empirical data in Takada et al. [26, Figure 4].Figure 8
Fitting for various geometric models to the empirical data in Fleury et al. [29, Figure 16].Figure 9
Fitting for various geometric models to the empirical data in Mogensen et al. [30, Figure 3].Figure 10
Fitting for various geometric models to the empirical data in Pedersen et al. [1, Figure 4], Aarhus.Figure 11
Fitting for various geometric models to the empirical data in Pedersen et al. [1, Figure 4], Stockholm.For different environments, the analysis ofDm 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 Dm 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 distribution of received signals for different environments. In the case of DBMM, it is different. The characteristics (IR, PDP, or PDS) that are the basis for this model define the propagation properties of the environment. As a result, DBMM has the ability to predict PDF of AOA for environments with different propagation properties.


## 6. Conclusion

In this paper, we presented analytical PDF model of AOA that takes into account the propagation properties of the environment on the basis of the temporal characteristics of channel. In the developed model (DBMM), PDFs of AOA for local and delayed scattering components are modeled by von Mises’ PDF and the distribution that results from geometry of the multielliptical propagation model. The parameters of each ellipse are defined on the basis of PDS or PDP. As a result, DBMM is closely related to the propagation conditions of environment and provides a flexible fit to the analyzed radio frequency range. For a given type of environment, DBMM is the spatial estimator of PDF of AOA, because the change of Tx/Rx position significantly affects the PDS and PDP. Therefore, the parameters of these characteristics are the estimation parameters of DBMM. Based on the measurement results presented in the literature, the mapping accuracy of the statistical properties of AOA was analyzed for varied propagation scenarios. To compare models, LSE andDm were adopted as a measure of the approximation and estimation accuracy, respectively. The obtained results show that DBMM provides the smallest LSE and Dm with respect to other models for all scenarios. This means that DBMM enables the estimation of PDF of AOA for different type of environments based on temporal channel characteristics and description of the propagation scenarios. In practice, the analytical simplicity of the model ensures its easy implementation for simulation studies that take into account the spatial propagation conditions in wireless channels. DBMM is the analytical description of the statistical properties of the reception angle in the Doppler multielliptical channel model (DMCM) [36, 37].The current trend in the development of channel models goes toward 3D modeling as it is shown in numerous papers, for example, [4, 11–13]. The presented model is limited to the azimuth plane, so in the near future, it will be used as a basis for the development of a 3D model.



---

*Source: 936406-2015-12-29.xml*