International Journal of Antennas and Propagation

Volume 2013, Article ID 176704, 17 pages

http://dx.doi.org/10.1155/2013/176704

## Quantification of Scenario Distance within Generic WINNER Channel Model

^{1}Department for Power, Electronic and Communication Engineering, Faculty of Technical Sciences, Trg D. Obradovića 6, 2100 Novi Sad, Serbia^{2}Electronic Measurement Research Lab, Ilmenau University of Technology, PSF 100565, 98684 Ilmenau, Germany^{3}Digital Broadcasting Research Laboratory, Ilmenau University of Technology, PSF 100565, 98684 Ilmenau, Germany

Received 27 July 2012; Revised 29 October 2012; Accepted 1 November 2012

Academic Editor: Ai Bo

Copyright © 2013 Milan Narandžić et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Starting from the premise that stochastic properties of a radio environment can be abstracted by defining scenarios, a generic MIMO channel model is built by the WINNER project. The parameter space of the WINNER model is, among others, described by normal probability distributions and correlation coefficients that provide a suitable space for scenario comparison. The possibility to quantify the distance between reference scenarios and measurements enables objective comparison and classification of measurements into scenario classes. In this paper we approximate the WINNER scenarios with multivariate normal distributions and then use the mean Kullback-Leibler divergence to quantify their divergence. The results show that the WINNER scenario groups (A, B, C, and D) or propagation classes (LoS, OLoS, and NLoS) do not necessarily ensure minimum separation within the groups/classes. Instead, the following grouping minimizes intragroup distances: (i) indoor-to-outdoor and outdoor-to-indoor scenarios (A2, B4, and C4), (ii) macrocell configurations for suburban, urban, and rural scenarios (C1, C2, and D1), and (iii) indoor/hotspot/microcellular scenarios (A1, B3, and B1). The computation of the divergence between Ilmenau and Dresden measurements and WINNER scenarios confirms that the parameters of the C2 scenario are a proper reference for a large variety of urban macrocell environments.

#### 1. Introduction

In order to maximize transmission efficiency, wireless communication systems are forced to exploit the spatial and temporal dimensions of the radio channel to the full. The design and performance analysis of such system requires the channel model to reflect all relevant propagation aspects, which imposes serious constraints on the minimal complexity of the model.

This paper concentrates on the class of geometry-based stochastic channel models (GSCMs) that offer good trade-off between complexity and performance (realism). These models deal with physical ray propagation and therefore implicitly or explicitly include the geometry of the propagation environment. The flexible structure of GSCMs enables the representation of different propagation environments by simple adjustment of model parameters, which is referred to as *generic* property. The generic models introduce abstract classes called *scenarios* that act as stochastic equivalents for many similar radio environments. As discussed in [1] these scenarios are not necessarily distinguished by the quantification of parametric space, but they represent a convenient terminology to designate typical deployment and propagation conditions. From history, it was COST207 [2] that started classifying environments based on the type of dispersion (delay spread and delay window), so some intuitive consideration of the (rather limited) parameter space was involved. But, these types of differences have (almost) never been sought later on, when defining new scenarios, especially following deployment schemes.

Nowadays we can distinguish between two major classes of generic models: COST 259/273/2100 ([3–5], resp.) and 3GPP SCM [6]/WINNER [7, 8]. The first one defines spatial regions where interacting objects become “visible,” that is, contributing to the total received field. The second class offers an abstraction of the environment in parametric space by using delay and angular spreads, cross-polarization, shadowing, -factor, and so forth. These generic models have been made by joint effort of many institutions; otherwise provision of parameters for different scenarios would be unattainable. Despite different model structures, significant overlap of propagation scenario definitions exists between SCM/WINNER and COST models [9].

The need for generic models follows from the ever growing concept of heterogeneous networks, requiring simultaneous representation of multiple scenarios or transitions between scenarios. For this purpose scenarios of generic models provide a uniform modeling approach and decrease the perceived complexity of handling different environments.

A reduction of the number of scenarios in generic models also reduces the necessary time and effort for design and performance evaluation of communication systems. Since every environment is specific, the classification of propagation environments into the different (reference) scenarios is not a simple task: how many classes suffice and how much divergence within a class should be tolerated? Obviously, a meaningful answer can be provided only if we have a metric to quantify the similarity between propagation environments. Providing such a metric is the main goal of this paper.

In absence of a scenario distance metric, reference scenarios are typically formed as a combination of system deployment schemes, mobility assumptions, and narrative description of environments, as illustrated in Section 2 for the WINNER reference propagation scenarios. The structure of the WINNER channel model and an approximation of its parametric space with multivariate normal distribution are given in Section 3. Section 4 describes measurement experiments used for the validation of the proposed distance metric and WINNER scenario parameters. The mean Kullback-Leibler divergence is introduced in Section 5 and is exploited to quantify the similarity between the approximated WINNER scenarios. Necessary modifications of the WINNER correlation coefficients are explained in the same section. Section 6 presents the results of measurement classification based on the introduced divergence measure, and Section 7 concludes the paper.

#### 2. WINNER Reference Propagation Scenarios

The WINNER (Wireless World Initiative New Radio) project [10] was conducted in three phases (I, II, and +) from 2004 until 2010, with the aim to define a single ubiquitous radio access system concept, scalable and adaptable to different short range and wide area scenarios. The effects of the radiopropagation on the overall system design are abstracted by the introduction of Reference Propagation Scenarios (RPSs). RPSs are related to WINNER system-deployment schemes, being suitably selected to represent different coverage ranges: wide area (WA), metropolitan area (MA), and local area (LA), and each deployment scheme was described by as few propagation scenarios as possible. The outcome is that the WINNER scenarios cover some typical cases, without the intent to encounter all possible propagation environments.

The WINNER reference propagation scenarios [7] are determined by the aspects that have immediate impact on the radio-signal propagation: (i)propagation environment, (a)LoS/NLoS condition, (b)limited distance range, (ii)terminal positions (heights) with respect to environment, (iii)mobility model (terminal speed), (iv)carrier frequency range/bandwidth.

Due to different propagation mechanisms under LoS and NLoS conditions, they are distinguished and separately characterized in all applicable physical environments.

All WINNER reference propagation scenarios are represented by generic channel model. This model, called WINNER channel Model (WIM), has been developed within the 3GPP Spatial Channel Model (SCM) framework. By its nature, these models are representing the wideband MIMO channels in static environments for nonstationary users. The MATLAB implementations of SCM and WIM are publicly available through the project website [10]. At the end of the phase II, WIM was parameterized for 12 different scenarios, being listed in Table 1. The full set of WIM RPS parameters can be found in Sections 4.3 and 4.4 of the WINNER deliverable D1.1.2 [7]. Relations between WINNER reference propagation scenarios and WIM parameters are illustrated in Figure 1.

The characterization of the reference propagation scenarios and parametrization of the generic model are based on channel sounding results. In order to collect relevant data, a large number of measurement campaigns have been carried out during the project. However, the realization of large-scale campaigns and the subsequent processing of the results are both complex and time consuming. As a consequence, the WINNER “scenario” is formed on the basis of measurement results that are gathered by different institutions and are individually projected on the parameter set of WINNER model. These measurements were conducted in radio environments providing the best possible match with defined reference scenarios. For that purpose, the position and movement of communication terminals were chosen according to the typical usage pattern. The resulting scenario-specific model parameters sometimes also include results found in the literature, in order to come up with the most typical representatives for a targeted scenario.

#### 3. Structure of the WINNER Channel Model

WIM is a double-directional [11] geometry-based stochastic channel model, in which a time-variable channel impulse response is constructed as a finite sum of Multi-Path Components (MPCs). The MPCs are conveniently grouped into clusters, whose positions in multidimensional space are determined by Large-Scale Parameter (LSP) realizations. LSPs are controlling the distribution of the power (spreading) over the individual dimensions of the channel, as indicated in Table 2.

Within the modeling context LSPs are exploited to govern the evolution of the synthesized channel. The entire process of WIM parameter synthesis can be done in three hierarchy levels [9]. (1)On top level, large-scale parameters listed in Table 2 are drawn randomly from tabulated log-normal Probability Density Functions (PDFs). With the exception of XPR, all other LSPs are generated as *correlated* random variables. (2)On the second level, cluster parameters are determined. For the sake of the simplicity this level of freedom is reduced in SCM/WIM, since all clusters share the same (scenario dependent) intra-Cluster Angular Spread (CAS). (3)In order to further simplify cluster characterization, SCM/WIM does not deal with random placement of MPCs in delay or angular domains. Instead, the same, simple internal structure of the cluster is used and MPC parameters are calculated in a deterministic manner. Following the given WIM approximations the LSPs become the most important for the particular scenario characterization. We would therefore ignore the lower hierarchy levels when computing a scenario divergence in Section 5.

##### 3.1. Transformed LSP Domain

The WINNER model investigates LSP distributions and their correlations in transformed domain [6] where normal distributions for all transformed LSPs are assumed. For log-normally distributed LSPs (delay and angular spreads) the mapping is applied (Figure 2). The remaining LSPs (SF, XPR, and -factor) have Gaussian distributions when expressed in (dB). The WINNER tables specify marginal (per-dimension) TLSP distributions with mean and variance parameters and matrix of correlation coefficients . The entries of this matrix express pairwise correlations of the LSPs and in the form of the correlation coefficient: where . Figure 3 shows marginal distributions of delay spread in transformed domain for all WIM reference scenarios having NLoS propagation and additionally includes the maximum likelihood estimate of normal distributions for Ilmenau and Dresden measurements, which will be discussed in the following sections.

###### 3.1.1. Multivariate Normal Distribution of TLSPs

The multivariate normal probability density function of a -dimensional random vector can be expressed as [12] where is -dimensional mean vector and covariance matrix is Although the structuring of the WIM TLSP distribution parameters is slightly different, they basically represent maximum-likelihood estimates of a multivariate normal (MVN) distribution parameters (3) and (4). (The random variables showing very specific dependence are not jointly normally distributed even if their marginal distributions are normal. Since only dependence between continuous WIM LSPs is expressed by the correlation coefficient (1), we assume without explicit proof that the vector of LSPs will have jointly normal distribution.) It is therefore possible to reconstruct the full covariance matrix of MVN by the following scaling: where represents the diagonal covariance matrix of the uncorrelated LSPs. Accordingly, every WINNER scenario can be abstracted with (up to) 8-dimensional normally distributed random process where relevant dimensions describe different large-scale parameters listed in Table 2. In a given case, the multivariate normal process offers a straightforward approximation of WINNER scenario since the majority of them have identical cluster structure.

##### 3.2. Terminal Separation

Although WIM describes the propagation environment implicitly within LSP parametric space, it is still using distance to govern transmission loss and spatial variations of LSP realizations.

Local stationarity region represents a larger area where multipath structure of physical propagation channel does not change significantly (“local region of stationarity” [13], “drop” [6], and “channel segment” [7]), and it is therefore characterized by a single realization from multidimensional LSP distribution. In WINNER model it is conveniently assumed that the extent of local stationarity region can be represented by the scenario-dependent constant (*decorrelation distance*) that is independent from the LSP cross-correlations.

Both transmission loss and decorrelation distance are deterministic features in WIM. They do not impact MVN distribution of TLSPs and could be analyzed independently. Therefore, we investigate MVN process as joint model for WINNER LSP marginal distributions and cross-correlation coefficients. This representation of multidimensional channel, on the scenario scale, can be considered as a generalization of the 1D small-scale fading channel approach, where stochastic properties of instantaneous envelope are characterized by PDF.

#### 4. Representation of Measurements in WINNER Parametric Space

The multidimensional sounding enables the investigation of the complete spatiotemporal structure of a radio channel that, additionally to the temporal delay of incoming waves, includes their angular directions at transmission and at reception as well as their polarizations. This can be achieved by the specialized estimation algorithms as RIMAX [14] when calibration data of double-polarized measurement antenna array is available [15, 16]. The multidimensional sounding is performed by dedicated RF equipment that sequentially transmits and measures channel responses between multiple antennas on transmitter and receiver sides. This approach requires high reliability of the time referencing of measurement data on both sides of radiolink, which is typically achieved by highly stable rubidium or cesium clocks.

##### 4.1. Measurement Campaigns

In this paper data from two measurement experiments will be exploited: the first one is performed in Ilmenau, Germany, in 2008 and the second in Dresden, Germany in 2009. In the rest of the paper they will be conveniently labeled as IL and DR, respectively. Both measurements are performed with Medav RUSK channel sounder [17] at 2.53 GHz using frequency bandwidth of 100 MHz. To allow high-resolution parameter estimations of multipath structure, dedicated antenna arrays at transmit and receive sides are used, providing a total of 928 MIMO subchannels (Figure 4). The time necessary to record the responses of all wideband MIMO sub-channels, , was 12.1 ms for Ilmenau and 24.2 ms for Dresden measurement. This limits the maximally allowed speed between the measurement vehicle and other interacting objects, , to 17.6 km/h and 8.8 km/h, respectively ( cm denotes the wavelength of the carrier). Under these conditions channels are properly sampled in space-time.

In both campaigns, three well-separated base station locations within city centers are used (Figure 5), and mobile terminal is positioned on the rooftop of the car. The same macrocell measurement setup, including the configuration of the measurement equipment, provides the proper base for comparison of Ilmenau and Dresden propagation environments. Ilmenau is a small city compared to Dresden: whereas Ilmenau is characterized by 3-4 floor buildings, Dresden has buildings with 6–8 floors. Subsequently the base station height at Dresden (*~*50 m) was almost doubled compared to Ilmenau (*~*25 m).

##### 4.2. Estimated MVN Distribution Parameters

Data from both measurement locations is used to estimate the parameters of WINNER model. The parameters of marginal LSPs and corresponding correlation coefficients estimated from Ilmenau and Dresden measurements are given in Tables 8 and 9 together with parameters describing WINNER reference propagation scenarios. Additional details regarding Ilmenau and Dresden measurements and analysis can be found in [18, 19], respectively.

According to the description of WINNER reference propagation scenarios, both measurements should be assigned to the typical urban macrocell scenario, C2. Since both measurements are conducted after the publication of WIM-C2 parameters, they provide a proper test set for the validation of reported WIM parameters.

###### 4.2.1. Normality of Estimated LSPs in Transformed Domain

Figure 6 shows maximum likelihood fit of the empirical delay spread CDF from Dresden measurement with WINNER-like log-normal distribution. Although this approximation seems reasonable, the null hypothesis, that collected and transformed DS samples coming from normal distribution, is rejected by both Lilliefors [20] and Jarque-Bera [21] normality tests. In contrast to the one-sample Kolmogorov-Smirnov test, these tests are suitable when parameters of null distribution are unknown and must be estimated. The normality hypothesis is also rejected for other large-scale parameters estimated from Ilmenau and Dresden data, for both LoS and NLoS conditions. Therefore the representation of measurements in WINNER parametric space can be considered as maximum-likelihood approximation of empirical multivariate distribution with MVN process.

#### 5. Quantification of Scenario Distance

As showed in Section 3.1.1, the parameter set of WINNER model is equivalent to the parameters of the multivariate normal distribution. Therefore, all measurement projections and reference scenarios share the same parameter set and could be treated as multivariate probability distributions. This enables the introduction of a metric to quantify the divergence (distance) between different projections of measurements on the parameter set of the model, including the representatives of different reference scenarios. In order to illustrate the (dis)similarity between B3 and C2 WINNER scenarios under LoS and NLoS propagation, joint 2D PDFs of delay spread and shadow fading are presented in Figure 7. We can observe that these scenarios have differences, but some kind of similarity measure will be useful. Having in mind that we want to quantify the distance between two distributions and , it is possible to apply some form of relative entropy, for example, Kullback-Leibler (KL) divergence [22]: where and denote the densities of and . The computation of the KL divergence according to (7) would require multidimensional mapping of the subset into two PDFs: and . This approach may become impractical for a large number of dimensions: in the case of WINNER it is necessary to consider up to 8 dimensions (although the XPR is not correlated with other LSPs). The marginal PDFs of -factor are given only for scenarios with LoS propagation what reduces the dimensionality of MVN distribution for NLoS propagation for 1.

In a special case when considering divergence between two MVN distributions it is possible to construct an analytical expression that depends solely on distribution parameters. The Kullback-Leibler divergence from to , for nonsingular matrices and , is [23]

This metric enables a simple comparison of reference WINNER scenarios. Therefore we found that the original form of KL metric (7) is more suitable for this particular problem than its symmetrized form, the Jansen-Shannon divergence [24]. Since KL divergence is not symmetric, it is necessary to define some other symmetrized extension to obtain proper distance metric. We propose to use mean KL divergence:

##### 5.1. Negative Definite Covariance Matrices

In some cases negative or complex values are obtained for KL divergence, indicating that the matrix of correlation coefficients () is not positive semidefinite, that is, . The problem is manifesting only for scenarios with resolved elevation angles (Table 3) where the dimensionality of the MVN distribution is increased from 6 (LoS)/5 (NLoS) to 7/6 or 8/7. The problem is, however, not related to the number of dimensions or elevation parameters themselves since simple removal of elevation dimension(s) does not resolve it. This means that correlation coefficients between WINNER LSPs analyzed jointly do not form a proper correlation matrix (CM)—not even without elevations.

It is observed that the number of decimal places used for representation of CM elements cannot be arbitrarily reduced since the resulting matrix may become negative definite. Since individual coefficients in WINNER parameter tables are expressed using only one decimal place, it is possible that this lack of precision causes negative definite CM for scenarios with an increased number of dimensions.

In order to enable the comparison of problematic scenarios their correlation coefficients have to be slightly modified to form positive definite CM. The “real” correlation matrix is computed using alternate projections method (APM) [25, 26]. For a given symmetric matrix this method finds the nearest correlation matrix , that is, (semi) definite, and has units along the main diagonal. The solution is found in the intersection of the following sets of symmetric matrices and . The iterative procedure in th step applies updated Dykstra’s correction and subsequently projects intermediate result to both matrix sets, using projections and : The projection replaces all negative eigenvalues of the matrix with a small positive constant , and forces ones along the main diagonal. The procedure stops when Frobenius distance between projections from two consecutive iterations drops below a predefined tolerance . Note, however, that the small tolerance parameter does not insure that Frobenius distance (FD) from the original matrix is equally small.

The positive definite approximation obtained by APM will depend on the selected parameters and [27]: for , the effect of eigenvalue on Frobenius distance is illustrated in Table 3. The results show that FD decreases when a smaller value of is used to substitute originally negative eigenvalues. However, the selection of small will proportionally increase the eigenvalues and coefficients of the inverse correlation matrix, . This will consequently increase the KL divergence (7) to all other scenarios. As a compromise, the new WINNER correlation coefficients corresponding to positive definite matrix are recomputed for and and given in Table 9. The maximal absolute modification of original correlation coefficients per scenario, , where , is given in Table 3. The highest absolute correction is applied to C2-NLoS scenario.

The minimum number of decimal places required to keep positive definite is determined for different values of and listed in Table 3. The results show that smaller Frobenius distance requires higher precision for saving coefficients. For two decimal places are sufficient to express correlation coefficients for all scenarios (Table 9).

##### 5.2. Mean KL Divergence

In order to enable comparisons between LoS and NLoS scenarios where -factor is missing under NLoS, as well as other scenarios where certain parameters (dimensions) are missing, the reduction of dimensionality was necessary: only those dimensions existing in both scenarios are used to calculate the mean KL divergence. This means that scenarios with lower number of resolved dimensions could exhibit more similarity as a consequence of incomplete representation. A fair comparison would be possible only if all scenarios have the same number of dimensions. The respective mean KL divergences between all WINNER scenarios, including Ilmenau and Dresden measurements, are given in Table 4 (for WINNER scenarios that give two sets of LoS parameters, mean KL distances are computed for LoS parameters before breakpoint distance of transmission loss).

In order to simplify the analysis of obtained results, for each (scenario, propagation) combination the closest match is found and listed in Table 5. Divergences within the same WINNER scenario group, or having same propagation conditions, are not minimum as may have been expected. Table 5 shows that only 5 among 16 WINNER scenarios have the closest match within the same WINNER group (A, B, C, and D). This comes as consequence of subjective classification of similar environments, without previously introduced metric. The minimum distances from Table 5, , confirm some expectations: B4-NLoS (outdoor-to-indoor) is closest to A2-NLoS (indoor-to-outdoor) because these are reciprocal scenarios. Also, microcell and macrocell versions of outdoor-to-indoor (B4 and C4) are the closest although not belonging to the same group. Mean KL divergences from Table 5 suggest that there is a better way to group available scenarios.

The average distances between all scenarios from one WINNER group to all scenarios in the other groups are given in Table 6. If all groups gather the most similar scenarios, an average distance between any two groups will be higher than the average distance within a single group. From Table 6 we can see that this applies to groups A and D, which are, according to average intergroup distance, closest to themselves. This indicates that subjective WINNER grouping can be partially supported by mean KL distance. However, the deviations are observable for groups B and C where other groups appear to be closer (at a lower average distance) than other scenarios from the same group. This situation is possibly caused by inappropriate assignment of distant scenarios, outdoor-to-indoor microcell scenario B4-NLoS and suburban C1-NLoS, to the corresponding groups.

An inspection of Table 3 reveals that the majority of scenarios with increased dimensionality (resolved elevations) come from groups B and C. This means that larger intragroup distance may appear due to increased dimensionality. Additionally, their correlation matrices have been modified by APM method, so that the parameter impacts the absolute value of the mean KL distance. The joint effect of these phenomena is illustrated by mean KL distances along the main diagonal of Table 6: they are proportional to the number of modified group members that are listed in Table 3: , , , .

For 6 out of 16 (scenario, propagation) pairs the best match has the opposite propagation condition (LoS, instead of NLoS, and vice versa), indicating that WINNER LoS and NLoS parameters do not form disjunctive sets (Table 5). Calculation of the mean distance between all LoS and NLoS scenarios in Table 7 shows that lower average distance can be expected between scenarios having LoS propagation condition (they are more similar than different scenarios with NLoS propagation).

#### 6. Classification of Measurements

The same criterion, KL divergence, can be applied to classify measurements as well. For this purpose even empirical distributions of LSPs can be used since KL metric (7) supports that. However, the extraction of the corresponding WINNER parameters simplifies the comparison since analytical expression (8) can be applied. Therefore we use the latter approach to compare Ilmenau and Dresden measurements with other WINNER scenarios.

Since both measurements have been performed in urban environments with macrocell setup (antennas were elevated above rooftops), it is expected that the closest scenario will be WINNER C2, which represents typical urban macro-cells. These expectations are met for Ilmenau measurements, where WINNER C2-NLoS is the closest scenario for both LoS and NLoS conditions, with minimal distances and (Table 4). In the case of Dresden measurements minimal mean KL divergences (8.6 and 11.6) indicate that the closest WIM scenario is C1-LoS, for both LoS and NLoS propagation conditions. This resemblance of Dresden measurements to suburban propagation (WINNER C1) may come from dominant height of BS positions with respect to environment.

Figure 8 shows the 2D PDFs of the reference WINNER C2-NLoS scenario together with joint LSP realizations from Ilmenau and Dresden measurements. For the NLoS propagation condition Ilmenau and Dresden measurements are quite close to C2: C2-NLoS is the best match for Ilmenau-NLoS () and the second best match for Dresden-NLoS data (). Additionally, among all results presented in Table 4 the closest match of WIM C2-NLoS is just Dresden-NLoS (row showed in red).

For LoS conditions, distances from WINNER C2 and Ilmenau and Dresden measurements are larger (35.7 and 23.9) which classifies Ilmenau-LoS to C2-NLoS () and Dresden-LoS into C1-LoS (). Table 5 shows the increased similarity between LoS and NLoS propagation conditions in Ilmenau and Dresden measurements. This occurs also for WINNER B3, while other WINNER scenarios do not show this property. One possible interpretation comes from the data segmentation into LoS and NLoS classes: the actual propagation conditions for the LoS or NLoS-labeled data may actually correspond to, for example, obstructed line of sight (OLoS). The previous analysis demonstrates that mean KL divergence, additionally to the comparison of different measurements, enables the quantification of complex relations between different data segments of the same measurement, as long as they use the same LSP space representation.

The mean KL distances between Ilmenau and Dresden measurements (38.3-LoS and 22.8-NLoS) are higher than corresponding distances from these measurements to the reference WINNER-C2 scenario. This confirms that WINNER C2 parameters provide appropriate representation for a wide class of urban macro-cell environments.

#### 7. Conclusions

The paper presents the scenario concept of WINNER and proposes its abstraction to a multivariate normal distribution of large-scale parameters. Disregarding transmission loss and decorrelation distance removes the spatial extent from scenario definitions.

The generic property of the model is exploited to compare the large-scale parameters that describe different scenarios. For this purpose, a symmetrized extension of the Kullback-Leibler divergence is proposed. This enables the comparison of parameters between reference scenarios and measurements, as well as a direct comparison of empirical LSP distributions (measured or synthesized by channel model). The given approach can be also applied to other generic stochastic models if appropriate metrics are chosen that reflect models’ specifics.

The presented results indicate that, according to the mean Kullback-Leibler divergence, WINNER scenario groups or propagation classes do not ensure the minimum separation within the group/class. It appears that other criteria, for example, coverage range, were more significant for the WINNER taxonomy. Judged from the mean Kullback-Leibler divergence large similarity exists between the indoor-to-outdoor and outdoor-to-indoor scenarios (A2, B4, and C4), between macro-cell configurations for suburban, urban, and rural scenarios (C1, C2, and D1), and between the indoor/hotspot/microcellular scenarios (A1, B3, and B1).

It is demonstrated that the results of measurements could be associated with the closest WINNER scenario. As expected, typical urban macro-cell scenario C2 was the one closest to the Ilmenau measurements. For the Dresden measurements, suburban LSP distributions of C1 were closest. This measurement, however, appears to be at minimum distance from C2-NLoS, indicating a validity of assumed macro-cell measurement setup. The proper choice of WINNER C2 parameters for representation of the whole class of urban macrocells is confirmed by the Ilmenau and Dresden measurements: these measurements are closer to the C2 reference than to each other.

For those scenarios/measurements where the correlation coefficients form a negative definite symmetric matrix, the alternating projection method is exploited to determine the closest correlation matrix. Therefore, the paper also introduces the modified WINNER scenario parameters that enable a quantification of scenario divergence.

#### Appendix

#### Parameters of WINNER Channel Model Describing MVN Distributions

In order to ensure the traceability of the presented divergences, the relevant subset of WINNER parameters is given in Tables 8 and 9. They also include the MVN distribution parameters estimated from Ilmenau and Dresden measurements. Additionally, Table 9 contains the modified correlation coefficients that form positive definite correlation matrices. They are used for scenario representation instead of original coefficients.

#### Acknowledgment

The authors would like to acknowledge the contributions of their colleagues from EASY-C project during the conduction and analysis of measurements, especially to Carsten Jandura, Actix GmbH.

#### References

- A. F. Molisch, H. Asplund, R. Heddergott, M. Steinbauer, and T. Zwick, “The COST259 directional channel model-part I: overview and methodology,”
*IEEE Transactions on Wireless Communications*, vol. 5, no. 12, pp. 3421–3433, 2006. View at Publisher · View at Google Scholar · View at Scopus - “Digital land mobile radio communications,” Tech. Rep. COST207, 1989.
- L. Correia, Ed.,
*Wireless Flexible Personalised Communications: COST 259, European Co-Operation in Mobile Radio Research*, Wiley, 2001. - L. Correia, “Mobile broadband multimedia networks—techniques, models and tools for 4G,”
*Final Report*COST 273, Elsevier, Oxford, UK, 2006. View at Google Scholar - R. Verdone and A. Zanella,
*Pervasive Mobile and Ambient Wireless Communications*, Signals and Communication Technology, Springer, London, UK, 2012. - 3GPP TR25.996 V9.0.0, Spatial channel model for MIMO simulations,Third Generation Partnership Project (3GPP), Tech. Rep., December 2009, http://www.3gpp.org/.
- P. Kyösti, J. Meinilä, L. Hentilä et al., “IST-4-027756 WINNER II Deliverable 1.1.2. v.1.2, WINNER II Channel Models,” Tech. Rep. IST-WINNERII, 2007. View at Google Scholar
- Guidelines for evaluation of radio interface technologies for IMTAdvanced, Report ITU-R M.2135-1, International Telecommunication Union, Radiocommunication Sector (ITU-R), 2009, http://www.itu.int/pub/R-REP-M.2135-1-.
- M. Narandžić, C. Schneider, and R. S. Thomä, “WINNER wideband MIMO system-level channel model, comparison with other reference models,” in
*Proceedings of the Internationales Wissenschartliches Kolloquium (IWK '09)*, vol. 54, Ilmenau, Germany, 2009. - WINNER (Wireless world INitiative NEw Radio), http://projects.celtic-initiative.org/winner+/.
- M. Steinbauer, A. F. Molisch, and E. Bonek, “The double-directional radio channel,”
*IEEE Antennas and Propagation Magazine*, vol. 43, no. 4, pp. 51–63, 2001. View at Publisher · View at Google Scholar · View at Scopus - “Multivariate normal distribution,” http://en.wikipedia.org/wiki/Multivariate_normal_distribution.
- A. Gehring, M. Steinbauer, I. Gaspard, and M. Grigat, “Empirical channel stationarity in urban environments,” in
*Proceedings of the European Personal Mobile Communications Conference (EPMCC '01)*, Vienna, Austria, February 2001. - R. Thomä, M. Landmann, and A. Richter, “RIMAX—a maximum likelihood framework for parameter estimation in multidimensional channel sounding,” in
*Proceedings of the International Symposium on Antennas and Propagation*, Sendai, Japan, August 2004. - M. Landman, M. Narandžić, and R. S. Thomä, “Limits of experimental channel characterization related to antenna calibration,” in
*Proceedings of the 29th General Assembly of the International Union of Radio Science (URSI '08)*, Chicago, Ill, USA, August 2008. - M. Landmann, M. Käske, and R. S. Thomä, “Impact of incomplete and inaccurate data models on high resolution parameter estimation in multidimensional channel sounding,”
*IEEE Transactions on Antennas and Propagation*, vol. 60, no. 2, pp. 557–573, 2012. View at Google Scholar - Rusk channel sounder, http://www.channelsounder.de/.
- C. Schneider, M. Narandžić, M. Käske, G. Sommerkorn, and R. S. Thomä, “Large Scale parameter for the WINNER II channel model at 2. 53 GHz in urban macro cell,” in
*Proceedings of the IEEE 71st Vehicular Technology Conference (VTC '10)*, Taipei, Taiwan, May 2010. - M. Narandžić, C. Schneider, M. Käske, S. Jäckel, G. Sommerkorn, and R. S. Thomä, “Large-scale parameters of wideband MIMO channel in urban multi-cell scenario,” in
*Proceedings of the 5th European Conference on Antennas and Propagation (EUCAP '11)*, pp. 3759–3763, Rome, Italy, April 2011. View at Scopus - H. W. Lilliefors, “On the Kolmogorov-Smirnov test for normality with mean and variance unknown,”
*Journal of the American Statistical Association*, vol. 62, no. 318, pp. 399–402, 1967. View at Google Scholar - C. M. Jarque and A. K. Bera, “A test for normality of observations and regression residuals,”
*International Statistical Review*, vol. 55, no. 2, pp. 163–172, 1987. View at Google Scholar - S. Kullback and R. A. Leibler, “On information and sufficiency,”
*The Annals of Mathematical Statistics*, vol. 22, no. 1, pp. 79–86, 1951. View at Google Scholar - W. D. Penny and S. J. Roberts, “Variational bayes for generalised autoregressive models,” Tech. Rep. PARG-00-12, Department of Engineering Science, Oxford University, 2000. View at Google Scholar
- B. Fuglede and F. Topsøe, “Jensen-Shannon Divergence and Hilbert space embedding,” http://www.math.ku.dk/~topsoe/ISIT2004JSD.pdf.
- J. P. Boyle and R. L. Dykstra, “A method for finding projections onto the intersection of convex sets in Hilbert spaces,”
*Advances in Order Restricted Inference*, vol. 37, pp. 28–47, 1986. View at Google Scholar - N. J. Higham, “Computing the nearest correlation matrix—a problem from finance,”
*IMA Journal of Numerical Analysis*, vol. 22, no. 3, pp. 329–343, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. P. García Ariza, J. F. Monserrat, and L. Rubio, “Alternating projection method applied to indefinite correlation matrices for generation of synthetic MIMO channels,”
*International Journal of Electronics and Communications*, vol. 64, no. 1, pp. 1–7, 2010. View at Publisher · View at Google Scholar · View at Scopus