Local Multipath Model Parameters for Generating 5G Millimeter-Wave 3GPP-like Channel Impulse Response

Local Multipath Model Parameters for Generating 5G Millimeter-Wave 3GPP-like Channel Impulse Response

Abstract

This paper presents 28 GHz and 73 GHz empirically-derived large-scale and small-scale channel model parameters that characterize average temporal and angular properties of multipaths. Omnidirectional azimuth scans at both the transmitter and receiver used high gain directional antennas, from which global 3GPP modeling parameters for the mean global azimuth and zenith spreads of arrival were found to be 22 and 6.2 at 28 GHz, and 37.1 and 3.8 at 73 GHz, respectively, in non-line of sight (NLOS). Small-scale spatial measurements at 28 GHz reveal a mean cross-polar ratio for individual multipath components of 29.7 dB and 16.7 dB in line of sight and NLOS, respectively. Small-scale parameters extracted using the KPowerMeans algorithm yielded on average 5.3 and 4.6 clusters at 28 GHz and 73 GHz, respectively, in NLOS. The time cluster - spatial lobe (TCSL) modeling approach uses an alternative physically-based binning procedure and recreates 3GPP model parameters to generate channel impulse responses, as well as new parameters like the RMS lobe angular spreads useful in quantifying millimeter-wave directionality. The TCSL algorithm faithfully reproduces first- and second-order statistics of measured millimeter-wave channels.

\newtoggle

conference \togglefalseconference

M. K. Samimi, T. S. Rappaport, “Local Multipath Model Parameters for Generating 5G Millimeter-Wave

3GPP-like Channel Impulse Response,” in the 10th European Conference on Antennas and Propagation

(EuCAP’2016), April 2016.

{IEEEkeywords}

28 GHz; 73 GHz; millimeter-Wave; multipath; angular spread; RMS delay spread; shadow fading; cluster; cross-correlation, cross-polar ratio; 3GPP; WINNER II; channel impulse response; 5G; time clusters; spatial lobes; TCSL.

\IEEEpeerreviewmaketitle

1 Introduction

The channel impulse response of a radio-propagation channel is composed of multipath components, whose local average temporal and angular characteristics can be obtained from large-scale and small-scale parameters. Large-scale parameters usually denote the azimuth spread (AS), the root-mean-square (RMS) delay spread, the shadow fading (SF), and the Rician -factor [1, 2], but must also be extended to include zenith (i.e., elevation) spreads, as directionality in azimuth and elevation is expected to drive future mmWave systems through multi-antenna arrays [3]. While large-scale parameters for channel impulse responses describe typical (average) local properties of a radio-channel (over a local area of tens of wavelengths), it is not to be confused with large-scale path loss which accounts for signal level fluctuations over several thousands of wavelengths resulting from large obstructions in the wireless channel (such as buildings). The 3GPP and WINNER II spatial channel models (SCMs) model the power delay profile (PDP) by utilizing the large-scale parameters as first- and second-order inputs to statistical distributions that govern the statistics of small-scale parameters of multipaths. Small-scale parameters refer to the properties of a single multipath component, such as path delay, and angles of arrival and departure. By way of contrast, the COST 2100 model does not explicitly generate large-scale parameters, but instead statistically generates an environment that is independent of mobile station (MS) location, and synthesizes the large-scale parameters based on cluster scattering [4]. 5G millimeter-wave (mmWave) channel models do not yet exist in the 3GPP and ITU standards, driving the need for channel model standards for future mmWave technologies, with large bandwidths (bit rates) and directional antenna gains (multi-element antenna arrays) [3].

The cross-correlation of large-scale parameters enhances the spatial consistency and accuracy of system-level simulations in which many users within close proximity must experience realistic correlated channels as observed in [5, 1, 2]. In the 3GPP model, the shadow fading values on the omnidirectional received powers for two users separated in space at the same time instant are generated from two Gaussian random variables (in dB) with a correlation coefficient of 0.5 [1], based on observations in [5]. The WINNER II urban microcellular (UMi) models generate the path delays from an exponential distribution whose mean is a function of the RMS delay spread (see Eq. (4.1) in [2]), and generate path angles using a wrapped Gaussian function whose standard deviation is a function of second-order measured angular spreads (see Eq. (4.8) in [2]). Both large-scale and small-scale path parameters are thus critical in generating wideband channel impulse responses that recreate the statistics of a large ensemble of collected measurements. The large-scale and small-scale parameters of 28 GHz and 73 GHz mmWave channels are provided, obtained from measurements carried out in 2012, 2013 [6, 3], and 2015 [7] in New York City.

2 Measurement Descriptions

Two outdoor propagation measurement campaigns were performed at 74 RX and 36 RX locations at 28 GHz and 73 GHz [6, 3], respectively, using a 400 megachips-per-second broadband sliding correlator channel sounder, and a pair of 24.5 dBi and 27 dBi directional antennas that provided over 12,000 directional PDPs to study AOD and AOA statistics of the mmWave channel [3, 8]. Additional 28 GHz small-scale fading measurements investigated the statistics of cross-polarization ratio (XPR) of individual multipath component amplitudes over an outdoor local area at one TX and four RX locations in line of sight (LOS), LOS-to-NLOS, and NLOS environments for distances ranging from 8 m to 12.9 m, using a pair of 15 dBi (28.8 and 30 half-power beamwidths in azimuth and elevation, respectively) gain horn antennas at the TX and RX, for vertical-to-vertical and vertical-to-horizontal polarization configurations. The measurements emulated a single-input multiple-output (SIMO) realistic base-to-mobile communication link, with the TX and RX placed 4 m and 1.4 m above ground, respectively. The PDPs were spatially sampled every 5.35 mm on a linear track by placing the RX antenna over two axes of a cross to emulate a virtual array (shown as two orthogonal arrows in Fig. 1), and capturing a PDP at static linear track position. Fig. 1 shows a map of the environment where the small-scale track measurements were collected [7].

Figure 1: Map of the environment in which the 28 GHz small-scale fading measurements were collected. One TX on Bridge (BRI) Street and four RX locations were selected in LOS (SS 1 and SS 2), LOS-to-NLOS (SS 3), and NLOS (SS 4) scenarios. The two orthogonal arrows represent the two axes of a cross, over which a PDP was sampled every 5.35 mm to emulate a virtual array of directional antenna elements.

3 Large-Scale Parameters

The global azimuth spread quantifies angular dispersion in the omnidirectional azimuth plane, computed as in (1[1]:

(1)

where,

(2)
(3)

where is the power azimuth spectrum, is the azimuth AOA of path , and . Since azimuth spreads are circular from 0 to , taking the minimum value with respect to the dummy variable removes the ambiguity. The global zenith spread can be computed from (1), (2), and (3), by replacing the azimuth angles with the elevation angles .

The global (i.e., omnidirectional) RMS delay spread is a measure of channel temporal dispersion, and is defined as the second central moment of a PDP [3]. The Rican -factor specifies the ratio of the strongest multipath power to the sum of powers of the other weaker multipaths, as in (4):

(4)

where is the total received power from all multipath components (e.g. area under the PDP curve). The -factors were computed to determine whether strong multipath components exist in LOS and NLOS environments.

The global zenith spreads of arrival (ZSA) were observed to exhibit a dependence on transmitter-receiver (T-R) separation , motivating a local mean model. The local means of the were modelled following the 3GPP approach [9]:

(5)

where the coefficients , , and are shown in Table 1, obtained using the minimum mean square error (MMSE) method, by minimizing the error between (5) and the local mean of values. The ZSDs were not computed because of the limited TX elevation diversity in the measurements [3].

The global parameters in Table 1 and cross-correlation coefficients in Table 2, and small-scale parameters in Table 3 were computed from a subset of previously used data presented in [8, 3], where all RX locations selected in this analysis had both an available omnidirectional PDP and TX/RX power angular spectra. The computed LSPs were obtained from 3 and 5 locations in LOS, and 13 and 19 locations in NLOS [10], at 28 GHz and 73 GHz, respectively.

Scenario Urban Microcell (UMi)
LOS NLOS
Parameters 28-73 GHz Combined 28 GHz 73 GHz
DS med (ns) 17.5 29.9 44.8
(ns) 26.6 42.1 45.9
log(DS) -7.71 -7.64 -7.53
0.34 0.50 0.51
ASD med 18.5 30.9 26.0
32.3 33.7 28.96
log(ASD) 1.28 1.38 1.34
0.50 0.41 0.39
ASA med 50.9 22.0 37.1
56.9 22.0 37.1
log(ASA) 1.69 1.39 1.50
0.27 0.39 0.20
ZSA med 4.0 6.1 3.5
4.0 6.2 3.8
log(ZSA) 0.6 0.72 0.55
0.09 0.30 0.15
E[ZSA] a 0.05 -0.002 -0.06
b -4.08 2.3 3.58
c 0.81 0.66 0.54
-factor 2.4 -0.4 1.5
2.0 4.3 6.8
Delay scaling 3.9 2.7 3.2
parameter 2.1 3.6 4.3
Table 1: Table summarizing large-scale parameters obtained from 28 GHz and 73 GHz channel measurements.
Scenario Urban Microcell (UMi)
LOS NLOS
Parameters 28-73 GHz Combined 28 GHz 73 GHz
ASD vs DS 0.32 -0.051 0.021
ASA vs DS 0.49 0.153 0.264
ASA vs SF 0.54 -0.637 0.044
ASD vs SF -0.04 0.051 0.008
DS vs SF 0.35 -0.508 -0.187
ASD vs ASA 0.72 0.405 -0.257
ASD vs K -0.16 -0.217 0.162
ASA vs K 0.07 -0.069 -0.428
DS vs K -0.46 -0.133 -0.449
SF vs K -0.03 -0.278 0.029
ZSA vs SF 0.16 -0.480 -0.327
ZSA vs K -0.37 -0.077 -0.105
ZSA vs DS 0.44 0.347 0.144
ZSA vs ASD 0.95 0.042 -0.027
ZSA vs ASA 0.72 0.323 0.081
Table 2: Cross-correlation coefficients obtained from 28 GHz and 73 GHz ultrawideband channel measurements.

4 Small-Scale Parameters

The small-scale parameters specify the characteristics of individual propagation paths, and are used to generate path powers, path delays, AODs, and AOAs, and are commonly obtained using high-resolution joint delay-angle multipath parameter extraction algorithms, such as the SAGE [11] and KPowerMeans algorithms [12]. The KPowerMeans algorithm provides an optimum assignment of multipath components into joint delay-angle clusters given a desired number of clusters from multi-dimensional channel impulse response data, using the mathematically-based power-weighted multipath component distance (MCD) metric [13]. The optimum number of clusters is then determined from two optimal criteria, the Caliñski-Harabasz and the Davies-Bouldin indices [13]. The KPowerMeans algorithm was run 50 times to remove the effects of initialization of centroid starting positions. The parameters , , and were set to 2, 0.9, and 0.9, respectively, for the combineValidate and shapePruning steps. The average number of clusters and cluster subpaths were determined to be 5.0 and 12.4 for the combined 28-73 GHz LOS scenario, respectively, 5.3 and 12.8 for the 28 GHz NLOS environment, and 4.6 and 13.2 for the 73 GHz NLOS environments, as summarized in Table 3.

Scenario Urban Microcell (UMi)
LOS NLOS
Combined 28-73 GHz 28 GHz 73 GHz
Clustering KPowerMeans
# of clusters 5.0, 3.5 5.3, 2.4 4.6, 3.3
# of subpaths 12.4, 14.8 12.8, 21.4 13.2, 22.1
Cluster ASD () 1.5, 2.2 3.0, 4.5 2.1, 6.9
Cluster ASA () 6.7, 16.1 9.6, 20.1 5.2, 12.1
Cluster ZSD () 0.8, 1.0 - 0.8, 1.15
Cluster ZSA () 1.8, 2.0 1.6, 3.5 1.5, 1.9
Per-cluster 13.6 16.1 17.4
shadowing (dB)
Table 3: Table summarizing small-scale model parameters, extracted using the KPowerMeans algorithm from 28 GHz and 73 GHz ultrawideband channel measurements.

A flaw of the KPowerMeans algorithm is that it fails to converge to one final solution as a result of the random initialization procedure, that assigns mutipath components to random delay-angle clusters. Every time the algorithm is run, the first assignment of multipaths to clusters is arbitrary, that can lead to differences in the final cluster partition. In addition, the Caliñski-Harabasz index carries an undefined numerator for a value of 1 cluster, as shown below [13]:

(6)

where,

(7)
(8)

where is an integer representing the desired number of clusters, is the total number of paths considered, is the number of paths in the th cluster, is the multipath component distance between vectors and , is the global centroid of the data, and is the centroid of the th cluster. It is clear that (6) is undefined for , thus suggesting that cluster can never be optimum for an arbitrary dataset. The fine-tuning of parameters , , , and the number of total runs of the algorithm are not explicitly discussed and remain open issues in [13]. In addition, the shapePruning step in [13] discards outlier data that do not significantly affect the properties of final clusters.

The time cluster - spatial lobe (TCSL) algorithm [8] offers an alternative physically-based binning procedure, where a time cluster corresponds to a group of traveling multipaths with similar delays but with potentially varying AOAs, and a spatial lobe denotes a strong direction of arrival (or departure) where energy is received contiguously in the azimuth and/or elevation dimensions. Statistics of time clusters are obtained using a physically-based 25 ns minimum inter-cluster void interval [8] in the time domain (representing 8 m in propagation distance, typical minimum spatial voids between narrow streets or buildings in New York City) by partitioning omnidirectional PDPs based on time of arrivals, since multipaths tend to separate due to free space gaps in the environment between buildings and other larger reflectors. Separately, the characteristics of spatial lobes are extracted by defining a -10 dB or -20 dB lobe power threshold with respect to the maximum peak received power in the 3-D power angular spectrum, where all contiguous angular segments above the lobe power threshold constitute one 3-D spatial lobe. After performing the separate time and spatial domain clustering, by randomly allocating the time delay multipaths to particular spatial lobes, the first- and second-order temporal and spatial statistics are well reproduced [14, 8]. Table 4 presents key parameters describing time clusters and spatial lobes obtained from all available omnidirectional PDPs and TX/RX angular spectra, used to generate mmWave channel coefficients [8].

Scenario Urban Microcell (UMi)
LOS NLOS
Parameters 28 GHz 73 GHz Combined 28-73 GHz 28 GHz 73 GHz
Number of clusters 3.0, 2.0 1.8, 0.8 2.3, 1.4 2.1, 1.4 2.7, 1.4
Number of subpaths 7.4, 5.8 7.8, 5.0 7.6, 5.2 9.1, 10.9 5.7, 6.0
Cluster decay constant (ns) 38.6 17.5 25.9 49.4 56.0
Per-cluster shadowing (dB) 1 3
Subpath decay constant (ns) 25.2 13.0 16.9 16.9 15.3
Per-subpath shadowing (dB) 6
# of AOD spatial lobes 3.3, 0.6 1.0, 0 1.9, 1.2 1.6, 1.8 1.5, 0.7
# of AOA spatial lobes 2.3, 1.5 1.4, 0.5 1.8, 1.0 1.6, 0.7 2.5, 1.1
AOD/AOA azimuth angles Uniform(0,360)
AOD/AOA elevation angles Gaussian (AOD) / Laplacian (AOA)
RMS lobe ASD (deg) 6.0, 3.2 4.7, 1.0 5.6, 2.7 6.2, 3.3 4.9, 2.7
RMS lobe ESD (deg) - - - - 2.2, 0.8
RMS lobe ASA (deg) 10.1, 3.9 3.9, 0.8 7.0, 4.2 6.8, 4.8 3.7, 2.3
RMS lobe ESA (deg) 11.2, 2.7 2.9, 0.7 7.0, 4.7 6.7, 2.3 2.2, 1.7
Table 4: Key parameters describing time clusters and spatial lobes to generate mmWave channel coefficients [8, 14]. A ‘-’ indicates that the data is unavailable from the measurements, which considered one omnidirectional azimuth scan for one fixed TX elevation downtilt at 28 GHz (LOS and NLOS), and 73 GHz LOS [3].

5 Number of Multipath Components for Directional Pointing Angles

Table 5 shows the number of multipath components measured at unique TX-RX directional pointing angles at 28 GHz and 73 GHz, obtained using a peak detection algorithm. The mean number of multipath components per unique pointing angle was obtained by considering all directional PDPs with at least 5 dB signal-to-noise (SNR) ratio, and then averaging the number of resolvable multipaths in all profiles over all RX locations. The mean number of multipath components per unique pointing angle was 7.2 and 3.8 at 28 GHz and 73 GHz in LOS, respectively, and 5.2 and 3.3 at 28 GHz and 73 GHz in NLOS, respectively. The mean number of multipath components at 28 GHz is larger than at 73 GHz in LOS environments, and may result from the larger beamwidth used at 28 GHz (10.9 beamwidth), and from the difference in measurement procedures, where the 28 GHz measurements used rigid TX and RX pre-determined angle combinations, whereas the 73 GHz angles were determined based on strongest received power from exhaustive initial searching during the field measurements. Since the 28 GHz measurements did not always consider strongest elevation planes, diffuse scattering was found to be more prominent at the weaker measured angles at 28 GHz, carrying many more multipath components than for the 73 GHz measurements, where fewer, but stronger multipath components were detected. Note that in [3], a different methodology was used to compute the mean number of multipath components at unique pointing angles: the numbers of detected multipath components in all directional PDPs satisfying at least a 5 dB SNR were first averaged over all pointing angles at each RX location, and then averaged again over all RX locations, yielding means of 4.7 and 3.3 at 28 GHz and 73 GHz in NLOS, respectively.

Number of Multipaths (Directional)
Scenario LOS NLOS
28 GHz 7.2, 5.3 5.2, 4.6
73 GHz 3.8, 3.1 3.3, 2.7
Table 5: Mean and standard deviation of the number of multipath components for directional TX-RX arbitrary pointing angles [3].

6 Cross-Polar Ratio (XPR) Measurements

28 GHz XPR small-scale fading measurements investigated the effects of polarization on resolvable multipath components by sampling PDPs at half-wavelength spatial increments over a 33-wavelength long linear track. For each 2.5 ns time delay bin, the XPR was extracted by dividing the received power (mW) from vertically-polarized TX and RX antennas, by the received power (mW) from vertically-polarized TX and horizontally-polarized RX antennas. Fig. 2 shows the cumulative distribution functions (CDFs) for the measured XPR in LOS, LOS-to-NLOS, and NLOS environments, and simulated Gaussian CDFs with values in dB. The mean and standard deviation of the XPR are summarized in Table 6, and were found to be 28.7 dB and 6.0 dB in LOS, 29.2 dB and 5.5 dB in LOS-to-NLOS, and 16.7 dB and 8.8 dB in NLOS, respectively. The XPR curves closely follow a Gaussian distribution (in dB), with the following probability density function (PDF):

(9)

where and are set to the mean and standard deviation (in dB) of the measurement set. Note that in Fig 2, the simulated data is obtained by taking the maximum with respect to 0 dB. In NLOS, the mean XPR is 12 dB smaller than in LOS, indicating that the propagating radio-waves experience depolarization from reflections and scattering in the environment.

Figure 2: Cumulative distribution function of 28 GHz ultrawideband cross-polarization ratio (XPR) measurements for LOS, LOS-to-NLOS, and NLOS environments, and corresponding simulated Gaussian distributions.
Scenario Outdoor Street Canyon
LOS LOS-to-NLOS NLOS
XPR (dB) 28.7 29.2 16.7
6.0 5.5 8.8
Truncated Distribution
Table 6: Means () and standard deviations () of the XPR obtained from small-scale fading track measurements at 28 GHz.

7 Conclusion

This paper presented 28 GHz and 73 GHz large-scale and small-scale model parameters required to generate the temporal and angular characteristics of multipath components in a PDP. The RMS DS and SF were found to be negatively correlated at both 28 GHz and 73 GHz, with a value of -0.508 and -0.187 in NLOS, respectively, while the ASA and SF were found to be negatively correlated at 28 GHz in NLOS with a value of -0.637 and nearly uncorrelated at 73 GHz in NLOS with a value of 0.044. Joint delay-angle clustering approaches, such as the KPowerMeans algorithm, use a mathematically-based clustering scheme to extract cluster statistics. Instead, the TCSL algorithm employs a physically-based thresholding technique separately in delay and angle to extract multipath properties [8]. The statistics presented here can be used in realistic system-level simulations and air-interface design [15] of next generation mmWave communication systems.

References

  1. “Spatial Channel Model for Multiple Input Multiple Output (MIMO) Simulations,” Tech. Rep. 3GPP 25.996 V12.0.0, Sept. 2014.
  2. P. Kyosti et al., “WINNER II channel models,” European Commission, IST-WINNER, Tech. Rep. D1.1.2, Sept. 2007.
  3. T. S. Rappaport et al., “Wideband Millimeter-Wave Propagation Measurements and Channel Models for Future Wireless Communication System Design,” IEEE Transactions on Communications, vol. 63, no. 9, pp. 3029–3056, Sept. 2015.
  4. L. Liu et al., “The COST 2100 MIMO channel model,” IEEE Wireless Communications, vol. 19, no. 6, pp. 92–99, December 2012.
  5. A. Algans et al., “Experimental analysis of the joint statistical properties of azimuth spread, delay spread, and shadow fading,” IEEE Journal on Selected Areas in Communications, vol. 20, no. 3, pp. 523–531, Apr. 2002.
  6. T. S. Rappaport et al., “Millimeter Wave Mobile Communications for 5G Cellular: It Will Work!” IEEE Access, vol. 1, pp. 335–349, 2013.
  7. M. K. Samimi and T. S. Rappaport, “28 GHz Millimeter-Wave Ultrawideband Small-Scale Fading Models in Wireless Channels,” in 2016 IEEE Vehicular Technology Conference (VTC-2016 Spring), May 2016.
  8. ——, “Statistical Channel Model with Multi-Frequency and Arbitrary Antenna Beamwidth for Millimeter-Wave Outdoor Communications,” in 2015 IEEE Global Telecommunications Conference (GLOBECOM), Workshop, Dec. 2015.
  9. “Study on 3D channel model for LTE,” Tech. Rep. 3GPP 36.873 V12.2.0, June 2015.
  10. G. R. MacCartney Jr. et al., “Exploiting directionality for millimeter-wave wireless system improvement,” 2015 IEEE International Conference on Communications (ICC), June 2015.
  11. H. B. Fleury et al., “Channel parameter estimation in mobile radio environments using the SAGE algorithm,” IEEE Journal on Selected Areas in Communications, vol. 17, no. 3, pp. 434–450, Mar. 1999.
  12. N. Czink et al., “A Framework for Automatic Clustering of Parametric MIMO Channel Data Including Path Powers,” in 2006 IEEE 64th Vehicular Technology Conference (VTC-2006 Fall), Sept. 2006, pp. 1–5.
  13. ——, “Improving clustering performance using multipath component distance,” Electronics Letters, vol. 42, no. 1, Jan. 2006.
  14. M. K. Samimi and T. S. Rappaport, “3-D statistical channel model for millimeter-wave outdoor mobile broadband communications,” in 2015 IEEE International Conference on Communications (ICC), June 2015, pp. 2430–2436.
  15. V. Fung et al., “Bit error simulation for pi/4 DQPSK mobile radio communications using two-ray and measurement-based impulse response models,” IEEE Journal on Selected Areas in Communications, vol. 11, no. 3, pp. 393–405, Apr. 1993.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minumum 40 characters
Add comment
Cancel
Loading ...
105458
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description