1 Introduction

For multipath propagation environment, mapping the statistical properties of the angle of arrival (AOA) is important in the analytical and simulation studies of signals in wireless systems. Part I of this paper [1] is focused on a comparative analysis of empirical and geometrical (theoretical) models of the probability density function (PDF) of azimuth AOA (AAOA) for different environmental conditions. Empirical models are based on standard distributions such as Gaussian [25], Laplacian [2, 3, 612], logistic [1316], and von Mises [3, 1721], which are being modified due to limited support of PDF of AAOA. The theoretical models are defined by the geometrical structures that describe the spatial position of scattering elements. In [1], the comparative analysis of PDF models also takes into account the results from [22], which relate to the best geometrical models such as uniform elliptical (receiver (Rx) outside) [23], Gaussian [6, 24, 25], and Rayleigh circular (Rx outside) [26]. A comparison shows that the smallest approximation errors of the measured data are provided by modified Laplacian and uniform elliptical (Rx outside) PDF for empirical and geometrical models, respectively. For all the analyzed scenarios and measures of the approximation accuracy, the average error is smallest for modified Laplacian PDF. The simplicity of the analytical description of empirical models, whose detailed analysis is derived from mathematical statistics, is an additional asset in relation to other models. This feature of the empirical models plays a significant role in theoretical and simulation studies. Therefore, in analytical analysis, e.g. [3, 17], which take into account the impact of AOA on correlational and spectral properties of received signals, the empirical models are primarily used. Numerous measurements, e.g. [7, 27, 28], show that the statistical properties of AOA significantly depend on the type of the propagation environment and spatial parameters of the measurement scenario. In practice, this means that the use of empirical models in analytical studies requires their adaptation to research scenario. In COST 207 [29], the assessment of the propagation conditions and classification of environment type is based on the rms delay spread (DS), σ τ . This parameter is defined on the basis of the impulse response (IR), power delay profile (PDP), or power delay spectrum (PDS). Therefore, the fit of the model to the type of propagation environment consists of determining the relationship between the model parameter and DS of the environment. It is the fundamental problem of the practical use of empirical models in analytical and simulation studies. Despite the large number of publications devoted to PDF of AOA modeling, the method of fitting the model to any environment is not described in any paper. This paper fills this gap. Here, the method for determining the relationship between DS and the parameter of each PDF model is presented for all environment types. In analytical and simulation studies, this method makes it possible to fit each analyzed model to different propagation environments.

In Part I of this paper [1], the accuracy of AAOA measurement mapping results by empirical and geometrical PDF models is shown. Part II is focused on the adaptation method of the empirical models to the different types of propagation environments. In Sect. 2, the relationship between DS and the rms angle spread (AS) is defined on the basis of the measurement scenarios taken from the open literature. In Sect. 3, the relationships between DS and parameters of PDF models are determined by minimizing the mean square error. Section 4 includes conclusions, which highlights the practical relevance of this method for analytical and simulation studies.

2 Statistical Relationship Between DS and AS

Evaluation of propagation environment properties is based on the signal power distribution over time that is described by such characteristics as IR, PSD, or PDP. These characteristics are the basis to determine σ τ that defines the environment properties in the time domain. However, the propagation environment is described by parameters that are defined not only in time domain but also in the reception angle domain. In this case, the properties of the propagation environment are described by the standard deviation of reception angle, σ θ , which is defined based on PDF of AOA. In practice, the evaluation of the statistical properties of reception angle is based on the measurement of the angular power spread. When analyzing the reception angle in the azimuth plane, σ θ is defined on the basis of the power azimuth spectrum (PAS), P(θ), as

$$\sigma_{\theta } = \sqrt {\frac{{\int\limits_{{ - 180^{ \circ } }}^{{180^{ \circ } }} {\theta^{2} P\left( \theta \right){\text{d}}\theta } }}{{\int\limits_{{ - 180^{ \circ } }}^{{180^{ \circ } }} {P\left( \theta \right){\text{d}}\theta } }} - \left( {\frac{{\int\limits_{{ - 180^{ \circ } }}^{{180^{ \circ } }} {\theta P\left( \theta \right){\text{d}}\theta } }}{{\int\limits_{{ - 180^{ \circ } }}^{{180^{ \circ } }} {P\left( \theta \right){\text{d}}\theta } }}} \right)^{2} }$$
(1)

In this case, σ θ is called AS. For the empirical PAS, σ θ is marked as σ E .

The results of the numerous measurements show that σ τ and σ E increase with increasing urbanization of environment. In this paper, the relationship between σ τ and σ E is determined based on data from the seven measurement scenarios described in [7, 27, 28, 3032]. In [22], many measurement scenarios are used to assess the approximation errors of geometrical models of PDF of AAOA. For analysis of empirical models, a number of measurement data is limited to these scenarios that include the signal characteristics in both time and angle domain. For all selected measurement scenarios, σ τ and σ E are determined on the basis of PDS or PDP and PAS, respectively. The results of numerical calculations are included in Table 1.

Table 1 DSs and ASs—results of numerical calculation for selected measurement scenarios

The correlation coefficient, ρ, is the basis for assessing the degree of DS and AS correlation [7]

$$\rho = \frac{{\sum\nolimits_{n = 1}^{N} {\left( {\sigma_{E} \left[ n \right] - \bar{\sigma }_{E} } \right)\left( {\sigma_{\tau } \left[ n \right] - \bar{\sigma }_{\tau } } \right)} }}{{\sqrt {\sum\nolimits_{n = 1}^{N} {\left( {\sigma_{E} \left[ n \right] - \bar{\sigma }_{E} } \right)^{2} } \sum\nolimits_{n = 1}^{N} {\left( {\sigma_{\tau } \left[ n \right] - \bar{\sigma }_{\tau } } \right)^{2} } } }}$$
(2)

where N = 7 is the number of scenarios described in Part I [1], σ E [n] and σ τ [n] are the estimated AS and DS for the nth scenario, while \(\bar{\sigma }_{E} = 7.328^{ \circ }\) and \(\bar{\sigma }_{\tau } = 0. 7 2 0 3\,\upmu{\text{s}}\) are the sample means of the estimated ASs and DSs, respectively.

For empirical data included in Table 1, the correlation coefficient is computed to be ρ = 0.8537. This value indicates that AS and DS are highly correlated. This means that between σ τ and σ E is the functional relationship, which makes it possible to assess the reception angle spread on the basis of the spread of signal in the time domain. In [7, Fig. 6], the analysis results of measurement data are a prerequisite to the adaptation of a linear function for mapping the relationship between σ τ and σ E . The linear regression line between AS and DS is obtained by using the least square method [33, 34]. For the linear function

$$\tilde{\sigma }_{\theta } = a \cdot \sigma_{\tau } + b$$
(3)

where a and b are the coefficients described by [33]

$$\begin{aligned} a & = \frac{1}{\Delta }\left( {N\sum\limits_{n = 1}^{N} {\sigma_{\tau } \left[ n \right]\sigma_{\theta } \left[ n \right]} - \sum\limits_{n = 1}^{N} {\sigma_{\tau } \left[ n \right]} \sum\limits_{n = 1}^{N} {\sigma_{\theta } \left[ n \right]} } \right) \\ b & = \frac{1}{\Delta }\left( {\sum\limits_{n = 1}^{N} {\sigma_{\theta } \left[ n \right]} \sum\limits_{n = 1}^{N} {\sigma_{\tau }^{2} \left[ n \right]} - \sum\limits_{n = 1}^{N} {\sigma_{\tau } \left[ n \right]} \sum\limits_{n = 1}^{N} {\sigma_{\tau } \left[ n \right]\sigma_{\theta } \left[ n \right]} } \right) \\ \Delta & = N\sum\limits_{n = 1}^{N} {\sigma_{\tau }^{2} \left[ n \right]} - \left( {\sum\limits_{n = 1}^{N} {\sigma_{\tau } \left[ n \right]} } \right)^{2} \\ \end{aligned}$$
(4)

Based on Table 1, the linear relation between σ E and σ τ is

$$\tilde{\sigma }_{E} \left[ ^\circ \right] = 4.65 \cdot \sigma_{\tau } \left[ {\upmu{\text{s}}} \right] + 3.98$$
(5)

The empirical data and the regression line for DS and AS are presented in Fig. 1.

Fig. 1
figure 1

Empirical data and regression line for σ τ and σ E

For the regression line, the root of the mean square error (RMSE), δ O is used to evaluate the approximation error

$$RMSE = \sqrt {\frac{1}{N}\sum\limits_{n = 1}^{N} {\left( {\sigma_{E} \left[ n \right] - \tilde{\sigma }_{\theta } \left[ n \right]} \right)^{2} } }$$
(6)

where σ E [n] and \(\tilde{\sigma }_{\theta } \left[ n \right]\) are empirical ASs from Table 1 and from regression line for the nth scenario, respectively. For the empirical data, \(\tilde{\sigma }_{\theta } \left[ n \right] = \tilde{\sigma }_{E} \left[ n \right]\).

On the basis of Eq. (6) and data contained in Table 1, the approximation error is computed to be RMSE = 1.4953°. The obtained result indicates that the approximation line yields an error of a value corresponding to the practical accuracy of the measurement angle. Later in this paper, RMSE is used to compare analyzed PDF models of AAOA.

3 Adaptation of the Empirical Models to Propagation Scenarios

The adaptation of the empirical model to the environment type involves selecting such value of model parameter that minimizes the approximation error with respect to measurement data. In Part I of this paper [1], the measure for the parameter selection and evaluation of the approximation error is the least-squares error (LSE). The optimal parameters of the models and corresponding LSEs are shown in [1, Table 4] and [1, Table 5], respectively. For all measurement scenarios and optimal parameters, each empirical model determines specified ASs that are shown in Table 2. To evaluate the accuracy of σ θ approximation by particular models, RMSEs relative to measurement results are also presented here.

Table 2 AS and RMSE for particular models and measurement scenarios

In time domain, each measurement scenario is identified by σ τ . Taking into account all measurement scenarios and their DSs, the data from Table 2 are the basis for the determination of the relationship between σ θ and σ τ for particular empirical models. The measurement results [7, Fig. 6] are a prerequisite to the adaptation of a linear function for the approximation of the relation between AS and DS for all models. Similarly as for (5), the regression lines for each model are obtained by using the method of least squares based on data from Tables 1 and 2. For particular models, the relationships are as follows:

  • modified Gaussian

    $$\tilde{\sigma }_{\theta } = 3.72 \cdot \sigma_{\tau } + 2.26$$
    (7)
  • modified Laplacian

    $$\tilde{\sigma }_{\theta } = 5.32 \cdot \sigma_{\tau } + 2.75$$
    (8)
  • modified logistic

    $$\tilde{\sigma }_{\theta } = 4.24 \cdot \sigma_{\tau } + 2.39$$
    (9)
  • von Mises

    $$\tilde{\sigma }_{\theta } = 3.73 \cdot \sigma_{\tau } + 2.26$$
    (10)

In Fig. 2, the graphical representations of these equations are presented.

Fig. 2
figure 2

Linear approximation of the relations between AS and DS for the empirical models

The graphs show that the approximation line for the modified Laplacian is the most convergent to the straight line for measurement data. The accuracy of mapping empirical relation of AS versus DS for modified Gaussian and von Mises PDF is practically the same. However, these models introduce a greater approximation error, especially for typical and bad urban environment.

The purpose of adapting the PDF model of AAOA is to provide such AS that corresponds to the propagation environment type. For each model, there is a close relationship between its parameter and σ θ . These relations for analyzed models are included in Table 3.

Table 3 AS versus parameters of particular PDF models

The above relationships are the basis for determining the parameters of models as a function of the environment type, that is, depending on σ τ . For the modified Gaussian and logistic models and taking (7), (9) into account, the relations between σ, s, and σ τ are, respectively

  • modified Gaussian

    $$\tilde{\sigma } = 3.73 \cdot \, \sigma_{\tau } + 2.23$$
    (11)
  • modified logistic

    $$\tilde{s} = 2.35 \cdot \, \sigma_{\tau } + 1.32$$
    (12)

The graphs of these equations are presented in Figs. 3 and 4, respectively.

Fig. 3
figure 3

σ versus σ τ for the modified Gaussian

Fig. 4
figure 4

s versus σ τ for the modified logistic model

In Table 3 for the modified Laplacian, the relationship between λ and σ θ is described as a hyperbolic function. Therefore, the relation λ versus σ τ is

$$\tilde{\lambda } = \frac{1}{{9.44 \cdot \sigma_{\tau } + 0.40}}$$
(13)

In Fig. 5, the graph of the equation is presented.

Fig. 5
figure 5

λ versus σ τ for the modified Laplacian

For the von Mises PDF, the complexity of relationship between κ and σ θ makes it difficult to use the least-squares method to determine the approximating function of κ versus σ τ . Therefore, in this case, the following approximation is adopted:

$${\sigma }_{\theta } = \sqrt {1 - \frac{{\mathrm{I}_{1} \left( \kappa \right)}}{{\mathrm{I}_{0} \left( \kappa \right)}}} \cong \frac{1}{{a_{M0} \cdot \kappa + b_{M0} }} + c_{M0} \cong \tilde{\sigma }_{\theta }$$
(14)

Taking (10) and (14) into account, the relation κ versus σ τ takes the form

$$\tilde{\kappa } = \frac{1}{{a_{M} \cdot \sigma_{\tau } + b_{M} }} + c_{M}$$
(15)

In (15), the parameters a M , b M , and c M are determined by the use of the least-squares method relative to optimal κ included in [1, Table 4] and data from Table 2. Obtained results are a M  = 0.0426, b M  = −0.0035, and c M  = 59.0287, so

$$\tilde{\kappa } = \frac{1}{{0.0426 \cdot \sigma_{\tau } - 0.0035}} + 59.0287$$
(16)

Graphical representation of (16) is shown in Fig. 6.

Fig. 6
figure 6

κ versus σ τ for the von Mises PDF

Based on the close relationship contained in Table 3 and for optimal κ from Table 2, RMSE is 2.687°, whereas for (14) is 2.852°. Therefore, the relative error of RMSE is about 1/16. It shows that (14) introduces about 16 times smaller error in relation to (10). This fact justifies the adaptation of (14) to approximate the relationship between σ θ and κ for the von Mises PDF.

For the particular empirical models, (11), (12), (13), and (16) provide the fit of models to the type environment defined based on σ τ . However, using these approximation increases the mapping error of AAOA statistical properties when compared with the optimal parameters. As in Part I [1], LSE is used to access the PDF mapping error of AAOA by empirical models. The basis for assessing the mapping error are the average LSE, for the optimal parameters and the average LSE R , for values from (11), (12), (13), and (16). For the particular models, the results of the numerical calculations of these measures and their difference are contained in Table 4.

Table 4 Average LSEs for the particular models

The obtained results make it possible to assess the mapping error of PDF of AAOA for each model. Comparison of LSE R and LSE shows that the adaptation method provides the smallest increment of mapping error for the von Mises PDF. A slightly greater increase in error occurs for the modified Laplacian. However, this model provides the smallest error in the evaluation of PDF of AAOA because it gives the smallest LSE and LSE R . Relatively worse mapping of real AAOA properties by the modified Gaussian, modified logistic, and von Mises PDFs do not exclude the possibility of their use in theoretical analysis. In many cases, the simplified description of the statistical properties of AAOA is ensured by PDF model that introduces a greater approximation error. In this paper, the presented adaptation method provides the opportunity to evaluate this error.

4 Conclusion

The paper discusses the adaptation of the PDF models of AAOA to study the propagation scenarios that take into account the diversity of environments. In this publication, the simplicity of the analytical description and a better fit to the measurement results are reasons to focus attention on the empirical models. Especially, these properties are important in analytical and simulation studies of signals in wireless systems. For diverse propagation environments, such PDF models of AAOA as the modified Gaussian, modified Laplacian, modified logistic, and von Mises have been considered in the presented comparative analysis. Similarly as in COST 207 [29], DS is adopted for the classification of the propagation environment type. The adaptation method of empirical model boils down to determining the relationship between DS and the model parameter. The obtained relationships make it possible to fit each PDF model to the real properties of AAOA that correspond to each propagation environment. For seven propagation scenarios, the measurement results are used as a reference data. A comparative analysis shows that the modified Laplacian model provides the smallest fitting error to measurement data. However, in theoretical and simulation studies, other empirical models can also be used due to the relatively small error and simple analytical description in relation to the geometric models. In this case, the presented adaptation method enables the assessment of the approximation error that results from the application of the empirical model.