METIS D1.4 v1.0
METIS D1.4 v1.0
METIS D1.4 v1.0
Project Name:
Mobile and wireless communications Enablers for the Twenty-twenty Information
Society (METIS)
Deliverable D1.4
Dissemination Level: PU
Contractual Date of Delivery: 28/02/2015
Security: Public
Status: Final
Version: 1.0
File Name: METIS_D1.4_v1.0.docx
Abstract
In this report, the end user scenarios, test cases and requirements envisioned by the METIS
project are mapped to propagation scenarios. Furthermore they are analysed for deriving a
new set of requirements relevant for radio channel and propagation modelling. Since none of
the existing channel models in the literature satisfies all these requirements, we develop new
channel models based on extensive measurement campaigns and analysis complemented by
computer simulations. The developed METIS channel models consist of a map-based model,
a stochastic model, and a hybrid model derived from both, to provide a flexible and scalable
channel modelling framework. In addition, guidelines including practical examples are
provided for utilizing the models in simulations.
Keywords
5G, channel model, stochastic, map-based, propagation, propagation scenarios, test cases,
channel measurements, D2D, M2M, V2V, backhaul, BS-UE, millimetre wave.
METIS Public ii
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Executive summary
METIS envisions a future where information access and sharing is available anywhere and
anytime to anyone and anything. The future information society will be provided with a wide
variety of applications and services including completely new industrial and professional
applications with diverse requirements. To meet these requirements, METIS sets the overall
technical goal to provide a system concept that, relative to today, supports a 1000 times
higher mobile data volume per area, a 10 to 100 times higher number of connected devices,
10 to 100 times higher user data rates, and 5 times reduced end-to-end latencies at a similar
cost and energy dissipation as today. In the early stage of the METIS project, 5 new scenarios
and 12 test cases have been identified which set new challenges for the envisioned 5G
technology components. A specific challenge is to provide the realistic and high quality radio
propagation models that a successful development and optimization of these technology
components require. The main objective of this deliverable is therefore to provide the METIS
channel models that are relevant in the future, particularly, for the 5G scenarios and test cases
envisioned by the METIS project.
To achieve this main objective, we start with the analysis of the test cases and requirements
from an end user perspective to derive a set of requirements relevant for radio channel and
propagation modelling. The identified requirements include an extremely wide frequency
range (up to 86 GHz and beyond), very high bandwidths (hundreds of MHz), fully three
dimensional and accurate polarization modelling, spherical wave modelling and high spatial
resolution, support of extremely large array antennas, dual-mobility for device-to-device (D2D),
machine-to-machine (M2M), vehicular-to-vehicular (V2V) communications, and spatial
consistency between link types (e.g. micro-/macro cellular, D2D) and between users/ devices
in a dense deployment. Through the literature study we conclude that none of the existing
channel models such as WINNER/ IMT-Advanced, COST 2100, and IEEE 802.11 fully satisfy
all these requirements and that consequently new channel models are needed.
To develop new channel models we have conducted extensive measurement campaigns. The
subsequent measurement data analyses have been complemented by computer simulations.
A comprehensive list of channel model parameters has been derived for diverse propagation
scenarios such as dense-urban macro, micro, indoor, shopping mall, D2D, and V2V links, with
a wide range of frequencies from 2 to 60 GHz based on the measurements.
The METIS channel models consist of a map-based model, a stochastic model, and a hybrid
model as a combination of both. This provides a flexible and scalable channel modelling
framework to meet diverse needs for simulations in terms of accuracy and computational
complexity. The map-based model is based on ray tracing using a simplified three dimensional
geometric description of the propagation environment and thus inherently accounts for
significant propagation mechanisms such as diffraction, specular reflection, diffuse scattering,
blocking, and so on. Therefore, the model provides accurate and realistic spatial channel
properties and is suitable for evaluating massive MIMO/ advanced beamforming, and also for
realistic pathloss modelling in case of D2D and V2V. Channel realizations are generated with
an implementation of the map-based model and are compared to the measurement results in
some selected scenarios by analysing propagation parameter statistics. The stochastic model
extends the geometry based stochastic channel model (GSCM), which has been further
developed from WINNER/3GPP, in order to provide multi-dimensional shadowing maps with
low complexity, millimetre-wave parameters, direct sampling of the power angular spectrum,
and frequency dependent pathloss models. The hybrid model provides a flexible and scalable
channel modelling framework. This is useful in balancing between the simulation complexity
and realism. For example, shadowing attenuation may be based on a map while small-scale
fading is stochastic. This deliverable also provides user guidelines for utilizing the channel
models in simulations. The guidelines include instructions on how to choose a channel model
depending on the propagation scenario of interest and the scope of the simulations by means
of some practical examples.
Contents
1 Introduction ...................................................................................................................... 1
1.1 Objective of the document ......................................................................................... 1
1.2 Structure of the document ......................................................................................... 3
2 5G Channel Model Requirements .................................................................................... 4
2.1 Wide range of propagation scenarios and network topologies ................................... 4
2.2 Spectrum related requirements ................................................................................. 4
2.3 Antenna related requirements ................................................................................... 4
2.3.1 Very large array antennas and Massive-MIMO .................................................. 4
2.4 Spatial consistency and mobility................................................................................ 5
2.5 Propagation related requirements ............................................................................. 5
2.5.1 Diffuse versus specular scattering ...................................................................... 5
2.6 General requirements ............................................................................................... 6
2.6.1 Complexity vs. accuracy..................................................................................... 6
2.7 D2D connections ....................................................................................................... 6
3 Propagation Scenarios ..................................................................................................... 7
3.1 Specified propagation scenarios ............................................................................... 7
4 Literature Review ............................................................................................................. 9
4.1 Available channel models.......................................................................................... 9
4.1.1 WINNER / IMT-Advanced .................................................................................. 9
4.1.2 COST 2100 ...................................................................................................... 10
4.1.3 IEEE 802.11 for 60 GHz ................................................................................... 10
4.1.4 Ray tracing ....................................................................................................... 10
4.2 Deficiencies in existing models................................................................................ 10
4.2.1 2D models ........................................................................................................ 10
4.2.2 Incomplete frequency coverage ....................................................................... 11
4.2.3 Bandwidth ........................................................................................................ 11
4.2.4 Limited number of scenarios ............................................................................ 11
4.2.5 Correlation inconsistency in D2D ..................................................................... 11
4.2.6 Inconsistent small-scale parameters ................................................................ 11
4.2.7 Dual mobility .................................................................................................... 11
4.2.8 Missing large array support .............................................................................. 12
4.2.9 Performance of the METIS models as regards the deficiencies........................ 12
4.3 Comparison of the models ...................................................................................... 12
4.4 Frequency dependent propagation effects .............................................................. 13
4.4.1 Material properties ........................................................................................... 13
4.4.2 Reflection from a wall and penetration though a wall ........................................ 13
4.4.3 Effect of surface roughness on reflection ......................................................... 13
4.4.4 Diffraction from the edge/ wedge...................................................................... 13
4.4.5 Diffuse scattering ............................................................................................. 14
4.4.6 Attenuation caused by oxygen and water vapour ............................................. 14
4.4.7 Rain attenuation ............................................................................................... 14
4.4.8 Vegetation loss ................................................................................................ 14
5 Findings from METIS Investigations, Measurements and Simulations ............................ 15
5.1 ITU-R UMi pathloss model fits our LOS measurements beyond 6 GHz but not
necessarily with our NLOS measurements ........................................................................ 15
5.2 A modified ITU-R UMi pathloss model covers the frequency range between 0.8 to
60 GHz .............................................................................................................................. 17
5.3 Scattered field patterns from a periodic surface show a critical frequency and
observation distance where the scattering effects start to be visible .................................. 20
5.4 Spatial channel model parameters available for dense UMi and D2D scenarios ..... 21
5.5 Supplemental parameters are derived for the 3GPP spatial channel model in UMa
scenarios at 2.3 GHz ......................................................................................................... 22
METIS Public iv
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
5.6 Elevation spread at the base station becomes smaller when elevating the base
station in UMi environments ............................................................................................... 24
5.7 3GPP spatial channel model parameters for inter-vehicular channels available at 2.3
and 5.25 GHz .................................................................................................................... 25
5.8 Simultaneous measurements of 10 and 60 GHz outdoor channels reveal consistent
delays of dominant propagation paths, while less multipaths were observed at 60 GHz .... 25
5.9 Diffraction can be a dominant propagation phenomenon in NLOS indoor 60 GHz
propagation........................................................................................................................ 27
5.10 Both specular and diffuse components are important in 60 GHz indoor channel ..... 27
5.11 The WINNER-type sub-path distribution model in a cluster is not suitable for 5G
channel models ................................................................................................................. 28
5.12 WINNER II channel model parameters are available for short-range 60-GHz links in
various environments ......................................................................................................... 29
5.13 Electrically very large array antenna makes the multi-user channels as ideal as the
i.i.d. condition..................................................................................................................... 31
5.14 Multiple-stream spatial transmission is possible at millimetre wave due to available
spatial degrees-of-freedom ................................................................................................ 33
6 Map-based Model .......................................................................................................... 35
6.1 Introduction ............................................................................................................. 35
6.2 Step-by-step instructions ......................................................................................... 35
6.3 Simplifications ......................................................................................................... 50
6.4 Future work ............................................................................................................. 51
7 Stochastic Model ............................................................................................................ 52
7.1 Parameter selection and procedure......................................................................... 52
7.1.1 Generate large-scale parameters ..................................................................... 54
7.1.2 Parameterisation (based on measurements / literature) ................................... 55
7.2 Antenna modelling .................................................................................................. 57
7.2.1 Spherical coordinate system ............................................................................ 57
7.2.2 Vector field rotation / mechanical tilting ............................................................ 58
7.2.3 Polarization transfer matrix ............................................................................... 59
7.2.4 LOS depolarization........................................................................................... 59
7.2.5 NLOS depolarization ........................................................................................ 60
7.3 Modelling V2V ......................................................................................................... 61
7.4 Sum-of-sinusoids approach for LSP generation ...................................................... 63
7.5 Ideas on dynamic modelling and spherical waves ................................................... 64
8 Guidelines for METIS Model Usage ............................................................................... 66
8.1 Model usage ........................................................................................................... 66
8.2 Map-based model ................................................................................................... 68
8.2.1 Massive deployment of sensors/actuators ........................................................ 69
8.2.2 Link-level simulation with precise route ............................................................ 69
8.2.3 Drop concept and signal level evaluations ....................................................... 69
8.2.4 Parameter ranges and default values ............................................................... 69
8.2.5 Simulation modes............................................................................................. 70
8.2.6 Simulator output ............................................................................................... 70
8.3 Stochastic model ..................................................................................................... 70
8.3.1 General settings ............................................................................................... 70
8.3.2 BS and UE locations ........................................................................................ 71
8.3.3 Antenna properties ........................................................................................... 71
8.3.4 System level description................................................................................... 71
8.4 Hybrid model ........................................................................................................... 74
9 Simulation Results ......................................................................................................... 75
9.1 Dense urban / Madrid .............................................................................................. 75
9.1.1 D2D.................................................................................................................. 75
9.1.2 Micro cell .......................................................................................................... 76
9.2 Specific LOS/NLOS simulations in urban microcell ................................................. 79
METIS Public v
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
List of Figures
Figure 5-1: Measurement environment and routes (left) and picture of LOS route from TX
antenna (right). ...................................................................................................................... 16
Figure 5-2: Comparison of measurement and M.2135 results. .............................................. 16
Figure 5-3: Frequency dependency of break point (BP) scaling factor (left) and pathloss offset
(right)..................................................................................................................................... 18
Figure 5-4: Example fits of the modified M.2135 UMi LOS model (5-4) and (5-5) at 26.36 GHz
(left) and 60 GHz (right); red curves are model, black dots are measurements. .................... 19
Figure 5-5: Example fits of the modified M.2135 UMi NLOS model (5-6) and (5-7) at
26.36 GHz (left) and 37.07 GHz (right). They were measured in two different streets with 𝒅𝟏
being 248 m (left) and 59.5 m (right); red curves are model, black dots are measurements. . 19
Figure 5-6: Frequency dependency of the pathloss offset (left) and shadow fading (right)..... 19
Figure 5-7: Analysis model (top left), frequency dependency with h = 0.1 m (bottom left), and
scattering effect on rough surface (right). .............................................................................. 21
Figure 5-8: Measurement environment during day- (left) and night-time (right)...................... 21
Figure 5-9: Received power in UMa O2I at 2.3 GHz; a prediction by a ray-tracing model is
overlaid. ................................................................................................................................ 24
Figure 5-10: ESD UMi O2O at 2.3 GHz (left), UMi O2I corridor at 2.3 GHz and 5.25 GHz
(right)..................................................................................................................................... 25
Figure 5-11: Location map (left), APDP in dBm for LOS measurement from P1 to P2 (right). 26
Figure 5-12: APDPs in dBm for LOS measurement from P3 to P1 (left) and NLOS
measurement from P1 to P4 (right). ...................................................................................... 26
Figure 5-13: Measurement scenario (left) and signal strength (right) relative to free space at
1 m distance for isotropic antennas measured and modelled at 2.4 GHz (blue dots) and
60 GHz (red dots) in a corridor of an indoor office scenario. .................................................. 27
Figure 5-14: Directional spread and angular spectra shown for delays indicated with numbers
1-5......................................................................................................................................... 28
Figure 5-15: Real measured [MAB+12] path directional distribution compared with WINNER
type (left) of distribution. The photograph (right) shows the measured paths (circles) as seen
from the base station location. ............................................................................................... 29
Figure 5-16: Power ordered MIMO channel eigenvalues for two different array antenna sizes
based on the WINNER model and the measured data. ......................................................... 29
Figure 5-17: Channel sounding with massive array antenna at the base station: array
antennas at the base station (top left), user terminal antenna (bottom left), and measurement
campaigns (right). ................................................................................................................. 31
Figure 5-18: Correlation properties of closely spaced devices under LOS condition; no user
proximity (top) and sum of normalized singular values with respect to the # of users (bottom
left). User and antenna numbering (bottom right). ................................................................. 32
Figure 5-19: CDF of the condition number for single-user MIMO channels (top) and channel
power variation over base station antenna elements and users (bottom). ............................. 33
Figure 5-20: SDoF of millimetre wave indoor channels: variation of the SDoF over different
antenna aperture size (left) and variation over different communication distances (right). ..... 34
Figure 5-21: Capacity for the single- and multi-stream spatial transmission for LOS channels
with 1λ2 (left) and 9λ2 (right) TX and RX antenna apertures. .................................................. 34
METIS Public ix
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 9-10: Reference pathloss plots between the access point and RX positions with no
obstructing/shadowing objects for AP antenna heights of 5 m (blue) and 10 m (red). ............ 81
Figure 9-11: Pathloss between access points and RX positions in the basic pedestrian
scenario for AP antenna heights of 5 m (blue) and 10 m (red). .............................................. 82
Figure 9-12: Mean pathloss between access points and RX positions in the basic pedestrian
scenario for AP antenna heights of 5 m (blue) and 10 m (red). .............................................. 82
Figure 9-13: Pathloss between the reference point (AP position in Figure 9-9) and RX
positions in the D2D user scenario. ....................................................................................... 83
Figure 9-14: Mean Pathloss between the reference point (AP position in Figure 9-9) and RX
positions in the D2D user scenario. ....................................................................................... 84
Figure A-1: PropSound system architecture. ......................................................................... 92
Figure A-2: TX ULA antenna for ME3, ME4 and ME5 scenarios at 2.3 GHz (left) and 5.25 GHz
(right)..................................................................................................................................... 93
Figure A-3: RX antenna for 2.3 GHz. ..................................................................................... 94
Figure A-4: RX antenna for 5.25 GHz. ................................................................................... 94
Figure A-5: Measurement setup for ME1. .............................................................................. 98
Figure A-6: TX leading and RX follows, vehicles move in the same directions. ..................... 98
Figure A-7: RX leading and TX follows, vehicles move to the same directions. ..................... 98
Figure A-8: TX (blue line) and RX (red line) routes moves to the opposite directions. ........... 99
Figure A-9: TX/RX (blue/red dot) stationary position and RX/TX (red/blue line) moves on the
cross street............................................................................................................................ 99
Figure A-10: TX antennas: 5.25 GHz dipole (SIMO) and 2.3 GHz ULA (MIMO). ................... 99
Figure A-11: RX measurement routes 1 and 2. ................................................................... 100
Figure A-12: RX measurement routes 3 and 4. ................................................................... 100
Figure A-13: RX measurement routes 5 and 6. ................................................................... 100
Figure A-14: RX (red line) measurement route 7. ................................................................ 100
Figure A-15: TX antenna on the top of building. .................................................................. 101
Figure A-16: Measurement spots in the hotel room for the ME3 scenario. ........................... 101
Figure A-17: The floor plan of the hotel with the measurement spots in the corridors on floors 3
to 5. ..................................................................................................................................... 102
Figure A-18: Magnified floor plan of the hotel corridors on floors 3 to 5. .............................. 102
Figure A-19: Measurement example, TX was located on the roof of neighbouring building. 103
Figure A-20: Views from the TX site towards the target building: from position 2 (left) and from
position 1 (right). ................................................................................................................. 103
Figure A-21: The hotel layout and measurement spots on the corridor of 2nd floor (left) and
6th floor (right). .................................................................................................................... 103
Figure A-22: Measurement spots in the hotel room and in the end of corridor at 2.3 GHz and
5.25 GHz. ............................................................................................................................ 104
Figure A-23: TX antenna distances at different heights in room measurements. ................. 104
Figure A-24: The TX mounted on an articulated crane at a height of 5 m (left) and 15 m (right).
............................................................................................................................................ 105
METIS Public x
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure A-25: TX antenna distances at different heights in corridor measurements. ............. 105
Figure A-26: Measurement spots on the hotel corridor for ME4 at 2.3 GHz and 5.25 GHz. . 105
Figure A-27: Measurement routes for ME5 at 2.3 GHz and 5.25 GHz. ................................ 106
Figure A-28: Pathloss, vehicles are moving to the same directions. .................................... 106
Figure A-29: Pathloss, vehicles are moving to the opposite directions................................. 106
Figure A-30: Received power, RX in the hotel room. ........................................................... 108
Figure A-31: Angular spreads for OLOS and NLOS cases. ................................................. 108
Figure A-32: UMi O2O pathloss at 2.3 GHz. ........................................................................ 111
Figure A-33: ESD vs. link distance in LOS scenario, TX height 5 m (left), TX height 10 m
(right)................................................................................................................................... 112
Figure A-34: ESD vs link distance in NLOS scenario, TX height 5 m (left), TX height 10 m
(right)................................................................................................................................... 113
Figure A-35: Measurement equipment: transmitter (left) and receiver (right). ...................... 114
Figure A-36: Measurement environment in front of Shibuya station. .................................... 115
Figure A-37: Measurement results. ..................................................................................... 116
Figure A-38: Pathloss measurement system at 26.35 GHz. ................................................ 118
Figure A-39: Map of Nihonbashi in Tokyo (typical Manhattan grid layout environment). ...... 119
Figure A-40: Photographs of the measurement environment. .............................................. 119
Figure A-41: Comparison of measurement and M.2135 results for multiple frequencies (same
as Figure 5-2). ..................................................................................................................... 120
Figure A-42: TX antenna height influence on the pathloss; M.2135 and measurement results
are compared at 37 GHz for different antenna heights. ....................................................... 121
Figure A-43: TX antenna height influence on the pathloss; M.2135 and measurement results
are compared at 4.7 GHz for different antenna heights. ...................................................... 121
Figure A-44: Measurement antenna patterns. ..................................................................... 122
Figure A-45: Measurement setup. ....................................................................................... 123
Figure A-46: Human body shadowing: Measurements vs. model at 60 GHz. ...................... 123
Figure A-47: Measurement locations for medium range measurement. ............................... 124
Figure A-48: Measurement setup for medium range measurement. .................................... 124
Figure A-49: Measured loss through door and window. ....................................................... 125
Figure A-50: Measurement locations for long range corridor measurement. ........................ 125
Figure A-51: Relative RX power measured and modelled at 2.4 GHz and 60 GHz for the long
range corridor measurement. .............................................................................................. 126
Figure A-52: Directional spread and angular spectra shown for delays indicated with numbers
1-5....................................................................................................................................... 127
Figure A-53: Measurement system and sounder configuration (left) and photograph of TX and
RX antennas (right). ............................................................................................................ 129
Figure A-54: Photographs of the measurement sites of the 60 GHz channel measurements at
the first and third floors of the Sello shopping mall............................................................... 129
METIS Public xi
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure A-55: Floor plan of the 1st floor of the Sello shopping mall with TX and RX locations.
............................................................................................................................................ 130
Figure A-56: Floor plan of the 3rd floor of the Sello shopping mall with TX and RX locations.
............................................................................................................................................ 130
Figure A-57: Number of measurements as a function of the TX-RX distance. ..................... 130
Figure A-58: Indoor cafeteria. .............................................................................................. 131
Figure A-59: Floor plan for cafeteria measurements. ........................................................... 131
Figure A-60: Open square measurement site. ..................................................................... 132
Figure A-61: Floor plan for open square measurements. ..................................................... 132
Figure A-62: Point cloud for open square in Kamppi, Helsinki. ............................................ 133
Figure A-63: Comparison of measured and predicted PDPs. .............................................. 133
Figure A-64: Examples of measured LOS and OLOS PDP’s with the detected peaks in the
shopping mall. ..................................................................................................................... 134
Figure A-65: Example of angular distribution of the detected peaks. ................................... 135
Figure A-66: Measured pathloss and pathloss model in shopping mall as a function of the link
distance............................................................................................................................... 135
Figure A-67: Simulated pathloss and pathloss model in cafeteria (left) and square (right) as a
function of the link distance. ................................................................................................ 136
Figure A-68: A CDF comparison to WINNER; measured delay spread in shopping mall LOS
(left) and OLOS (right). ........................................................................................................ 136
Figure A-69: A CDF comparison to WINNER; simulated delay spread in cafeteria LOS. ..... 136
Figure A-70: A CDF comparison to WINNER; simulated delay spread in square LOS (left) and
OLOS (right). ....................................................................................................................... 137
Figure A-71: A CDF comparison to WINNER; measured K-factor in shopping mall (left) and
cafeteria (right). ................................................................................................................... 137
Figure A-72: A CDF comparison to WINNER; simulated K-factor in square LOS. ............... 137
Figure A-73: A CDF comparison to WINNER; measured ASD in shopping mall LOS (left) and
OLOS (right). ....................................................................................................................... 138
Figure A-74: A CDF comparison to WINNER; simulated ASD in cafeteria LOS. .................. 138
Figure A-75: A CDF comparison to WINNER; simulated ASD in square LOS (left) and OLOS
(right)................................................................................................................................... 138
Figure A-76: A CDF comparison to WINNER; measured ASA in shopping mall LOS (left) and
OLOS (right). ....................................................................................................................... 139
Figure A-77: A CDF comparison to WINNER; simulated ASA in cafeteria LOS. .................. 139
Figure A-78: A CDF comparison to WINNER; simulated ASA in square LOS (left) and OLOS
(right)................................................................................................................................... 139
Figure A-79: Simulated ESD (left) and ASA (right) in cafeteria LOS and square LOS and
OLOS. ................................................................................................................................. 140
Figure A-80: Measured XPR, cross-polarisation measurement noise limit, and mean (black
solid line) and the 95 % tolerance interval limits (black dashed lines) fitted normal distribution.
............................................................................................................................................ 140
Figure A-81: a) Example of a measured power and b) cluster power over about 12 wavelength
long measurement path. ...................................................................................................... 141
Figure A-82: Average standard deviation and maximum-to-minimum difference of the cluster
power within 20 dB dynamic range as a function of cluster AS over a) four wavelengths and b)
12 wavelengths. .................................................................................................................. 141
Figure A-83: CDF of delay scaling parameter and per-cluster fading in shopping mall. ....... 142
Figure A-84: CDF of delay scaling parameter and per-cluster fading in cafeteria. ............... 142
Figure A-85: CDF of delay scaling parameter and per-cluster fading in square. .................. 142
Figure A-86: A CDF comparison to WINNER; measured maximum excess delay in shopping
mall LOS (left) and OLOS (right). ........................................................................................ 143
Figure A-87: A CDF comparison to WINNER; simulated maximum excess delay in cafeteria
LOS. .................................................................................................................................... 143
Figure A-88: A CDF comparison to WINNER; simulated maximum excess delay in square
LOS (left) and OLOS (right). ................................................................................................ 143
Figure A-89: Simulated correlation distances in cafeteria LOS. ........................................... 144
Figure A-90: Simulated correlation distances in square LOS (left) and OLOS (right). .......... 144
Figure A-91: A CDF comparison to WINNER; measured number of propagation paths in
shopping mall LOS (left) and OLOS (right). ......................................................................... 144
Figure A-92: A CDF comparison to WINNER; simulated number of propagation paths in
cafeteria LOS. ..................................................................................................................... 145
Figure A-93: A CDF comparison to WINNER; simulated number of propagation paths in
square LOS (left) and OLOS (right). .................................................................................... 145
Figure A-94: Channel sounder overview. ............................................................................. 146
Figure A-95: Horizontal (top) and vertical cuts (bottom) of the normalized antenna power
patterns in dB at 10 (left) and 60.4 GHz (right). ................................................................... 146
Figure A-96: Measurement environment including TX locations and RX routes. .................. 147
Figure A-97: On-site measurement setup. ........................................................................... 148
Figure A-98: Pathloss results at 10 and 60.4 GHz for LOS (left) and NLOS (right) condition.
............................................................................................................................................ 148
Figure A-99: Exemplary normalized APDP for 10 (left) and 60 GHz (right). ......................... 149
Figure A-100: Delay spread CCDF for LOS (left) and NLOS (right) condition. ..................... 149
Figure A-101: Planned measurement locations. .................................................................. 151
Figure B-1: Indoor office. [MET13-D61] ............................................................................... 153
Figure B-2: Madrid grid for Urban Micro and Urban Macro. [MET13-D61] ........................... 153
Figure B-3: Shopping mall. [MET13-D61] ............................................................................ 154
Figure B-4: Stadium [MET13-D61]. ..................................................................................... 154
Figure B-5: Highway. ........................................................................................................... 155
Figure C-1: Example of a street corner acting as a node (left), Manhattan map (middle),
topological example with four nodes (right). ........................................................................ 157
Figure C-2: Determination of interaction nodes. .................................................................. 158
Figure C-3: Determined pathways for the Berg recursive diffraction model. Here up to four
interactions of specular reflection and diffraction have been accounted for. ........................ 159
Figure C-4: Diffracted paths between TX and RX (left) and relative power over the RX route
(right)................................................................................................................................... 160
Figure C-5: Received power over the RX route for isotropic antennas and 0 dBm transmit
power. ................................................................................................................................. 160
Figure C-6: Path angles and propagation distances at RX along the route for direct paths only.
The power scale is relative to the strongest path for each RX location. ............................... 160
Figure C-7: Shadowing screen model.................................................................................. 161
Figure C-8: Diffracted paths between TX and RX (left), and, relative power over the RX route
(right). Obstructing objects are shown with black dots. ........................................................ 162
Figure C-9: Body blocking loss for a LOS link (4 m distance) at 60 GHz. ............................ 162
Figure C-10: Received power at 2 GHz over the RX route for isotropic antennas and 0 dBm
transmit power for 1.5 m RX height. The upper curve (blue) corresponds to the case with no
obstructing objects and the middle curve (red) to the case with obstructing objects when the
TX is at 10 m height and the lower curve (green) to the case when TX is at 1.5 m height. .. 163
Figure C-11: Schematic drawing of the scattering model. .................................................... 164
Figure C-12: Paths between TX and scatterers around one RX location (upper left) and paths
between one RX location and scatterers around TX (upper right), and, relative power over the
RX route due to scatterers around RX (lower left) and scatterers around TX (lower right). .. 164
Figure C-13: Distributions of path angles (left) and propagation distances (right) at RX for
paths between TX and scatterers around RX (lower) and paths between RX and scatterers
around TX (upper). The power is relative to the strongest path (LOS or diffracted) shown in
Figure C-6. .......................................................................................................................... 165
Figure C-14: Model for canonical problem of diffraction ....................................................... 166
Figure C-15: Example of determination of outdoor to indoor paths. In a) the paths are
identified with the building removed. In b) and c) the building is reintroduced in order to
determine penetration loss due to walls and floors. ............................................................. 168
Figure C-16: Combination of an outdoor map with an indoor layout. a) Example of
determination of outdoor to indoor paths. b) Virtual office layout (red drawing) inside a building
block of Madrid map. ........................................................................................................... 169
Figure C-17: Example of determination of over roof-top paths. In a) the paths are identified in
the horizontal plane assuming the building where the access point is located is removed. The
3D paths are blocked by the building and corresponding diffraction points are found at the
roof-top edge as shown on in b). ......................................................................................... 169
Figure D-1: Channel coefficient generation procedure [3GPP14-36873]. ............................ 172
Figure D-2: Cellular grid layout [ITUR09-2135]. ................................................................... 173
Figure D-3: 2D and 3D distances for outdoor (left) and indoor (right) UEs. .......................... 175
Figure D-4: Implementation of Laplacian shape. ................................................................. 178
Figure E-1: Reflection from a wall and penetration through a wall. ...................................... 187
Figure E-2: Reflection coefficient for concrete (45 deg incident angle). ............................... 188
Figure E-3: Penetration loss for a concrete wall of 5, 10 or 15 cm as a function of frequency for
three different angles of incidence. ...................................................................................... 189
Figure E-4: Attenuation due to knife-edge diffraction as function of 𝝂. ................................. 191
METIS Public xv
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
List of Tables
Table 3-1: Propagation scenarios. ........................................................................................... 8
Table 4-1: Comparison of different models. ........................................................................... 12
Table 5-1: Summary of major findings from METIS measurements and simulations. ............ 15
Table 5-2: Measurement results in the Tokyo downtown crowded area for UMi and D2D
scenarios............................................................................................................................... 22
Table 5-3: Parameters for UMa O2O scenario. ..................................................................... 23
Table 5-4: Spatial parameters for V2V scenario. ................................................................... 25
Table 5-5: WINNER parameterisation for 60 GHz range; parameters affecting validity range of
the derived parameters. ........................................................................................................ 30
Table 5-6: WINNER parameterisation for 60 GHz range; channel model parameters............ 30
Table 6-1: Parameter table for the map based model. ........................................................... 35
Table 7-1: Recommended PL and fading models for each propagation scenario. ................. 53
Table 7-2: Parameter ranges for UMi and UMa. .................................................................... 55
Table 7-3: Parameter ranges for D2D&V2V and Indoor. ........................................................ 56
Table 7-4: V2V channel model parameters for highway and rural scenario. .......................... 62
Table 8-1: Applicability of METIS channel models. ................................................................ 67
Table 8-2: Comparison of METIS models. ............................................................................. 68
Table 9-1: Cross-correlations of LS parameters (UMi LOS). .................................................. 77
Table 9-2: Cross-correlations of LS parameters (UMi NLOS). ............................................... 78
Table 9-3: Parameter table for the specific LOS/NLOS simulations with the map-based model.
.............................................................................................................................................. 80
Table 9-4: Excess pathloss due to obstructing/shadowing objects in Madrid grid scenario. ... 85
Table 9-5: Pathloss comparison between simulations and measurements in urban microcell.85
Table 9-6: Model parameter used in simulations. .................................................................. 86
Table 9-7: Simulation results for Madrid D2D and Madrid Micro scenario. ............................. 87
Table A-1: Overview of measurement campaigns within METIS............................................ 90
Table A-2: The summary of the measurement environments at 2.3 and 5.25 GHz in Oulu
downtown. ............................................................................................................................. 91
Table A-3: TX antenna properties.......................................................................................... 93
Table A-4: RX antenna properties. ........................................................................................ 93
Table A-5: Settings for measurements at 2.3 GHz................................................................. 94
Table A-6: Settings for measurements at 5.25 GHz............................................................... 96
Table A-7: V2V measurement results (ME1). ...................................................................... 107
Table A-8: Parameter for UMa O2I corridor scenario at 2.3 GHz (ME3). ............................. 108
Table A-9: Parameters for UMi O2I room scenario. ............................................................. 109
Table A-10: Parameters for UMi O2I corridor scenario. ....................................................... 110
Table A-11: Parameters UMi O2O scenario at 2.3 GHz. ..................................................... 111
List of Abbreviations
3GPP 3rd generation partnership project METIS Mobile and wireless communications
APDP Average power delay profile enablers for the twenty-twenty
AOA Azimuth angle of arrival information society
AOD Azimuth angle of departure MIMO Multiple input multiple output
AP Access point MMC Massive machine communications
ASA Azimuth Spread of Arrival MoM Moment method
ASD Azimuth Spread of Departure MPC Multipath component
BS Base station MU-MIMO Multi-user MIMO
BH Backhaul NLOS Non line of sight
CDF Cumulative density function O2I Outdoor to indoor
D2D Device-to-device O2O Outdoor to outdoor
DS Delay spread OLOS Obstructed LOS
EOA Elevation angle of arrival PDP Power delay profile
EOD Elevation angle of departure PL Pathloss
ESA Elevation Spread of Arrival PS Propagation scenario
ESD Elevation Spread of Departure RX Receiver
FBC First bounce cluster SAGE Space alternating generalized
GSCM Geometry based stochastic channel expectation-maximization
model SBC Single bounce cluster
HT Horizontal topic SCM Stochastic channel model
H-pol Horizontal polarization SDoF Spatial degrees of freedom
I2I Indoor to indoor SF Shadow Fading
IEEE Institute of electrical and electronics SS Small-scale
engineers TC Test case
i.i.d. Independent and identically distributed TX Transmitter
IMT- International mobile UE User equipment
Advanced telecommunications-advanced UMi Urban micro
ISD Inter-site distance UMa Urban macro
ITU International telecommunication union UTD Uniform geometrical theory of
ITU-R ITU radiocommunication sector diffraction
KF K-Factor V2V Vehicle-to-vehicle
LBC Last bounce cluster VNA Vector network analyser
LOS Line of sight VPL Vertical plane launch
LS Large-scale V-pol Vertical polarization
LSP Large-scale parameter WINNER Wireless world initiative new radio
M2M Machine-to-machine WLAN Wireless local area network
MED Maximum Excess Delay WP Work package
XPD Cross polarization discrimination
XPR Cross polarisation ratio
1 Introduction
1.1 Objective of the document
The main objective of the METIS (Mobile and wireless communications Enablers for the
Twenty-twenty Information Society) project is to lay the foundation of 5G, the next generation
mobile and wireless communications system for year 2020 and beyond. The project provides
technical enablers needed to address the requirements foreseen for this time frame.
METIS’ vision is a future where access to information and sharing of data is available
anywhere and anytime to anyone and anything. The future information society of private and
professional users will be provided with a wide variety of applications and services, ranging
from infotainment services, through increased safety and efficient usage of transportation, to
completely new industrial and professional applications [MET13-D11]. Realizing this vision
calls for solutions to challenges such as the provisioning of very high data rates, and the
handling of very dense user crowds, with higher requirements on the end-to-end performance
and user-experience. Other challenges that arise from new application areas are the
requirement on very low latency, very low energy consumption, hence cost, as well as the
support of a massive number of devices. As a consequence, the METIS overall technical goal
provides a system concept that, relative to today, supports 1000 times higher mobile data
volume per area (10 to 100 times higher number of connected devices and 10 to 100 times
higher user data rate), 10 times longer battery life for low power massive machine
communication (MMC), and 5 times reduced end-to-end latency, all of them at a similar cost
and energy dissipation as today.
To meet these demands, today’s spectrum can be used more efficiently, e.g. through higher
spectral efficiency in bits/s/Hz and increased spectrum utilization through dynamic spectrum
access. New frequency bands on cm-waves and millimetre waves (above 3 GHz up to 86 GHz
[MET13-D51]) provide much more spectrum than available today for 4G. Furthermore, system
capacity and coverage can be improved via new network topologies and technologies such as
moving networks, multi-hop communications, self-configuration networks, and direct device-to-
device (D2D) communications.
All the above-mentioned aspects set new requirements for radio channel and propagation
modelling as follows (more detailed requirements are described in Section 2).
Extremely wide frequency range from sub-1 GHz to 86 GHz and beyond
Very high bandwidth (> 500 MHz)
Full 3-Dimensional and accurate polarization modelling
Massive-MIMO: spherical waves instead of plane wave assumption and very high
spatial resolution
Extremely large arrays even beyond stationarity interval
Direct D2D, M2M and V2V Communication
Wide range of propagation scenarios and network topologies (from stationary to very
high speed, outdoor-to-indoor, from single antenna to massive arrays, from single link
to mesh network etc.)
Spatial Consistency between topologies and between users
o geographical locations of first and last bounce scatterers
o birth-death process and/or visibility regions for clusters
Importance of diffuse vs. specular scattering, especially in the high frequency domain
Currently recognized and widely used channel models, e.g. 3GPP/3GPP2 Spatial Channel
Model (SCM) [3GPP03-25996], WINNER [WIN208-D112; WIN+10-D53], and ITU-R IMT-
Advanced [ITUR09-2135], 3GPP 3D-UMi and 3D-UMa [3GPP14-36873], and IEEE 802.11ad
METIS Public 1
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
[MEP10] were found to be inadequate for 5G requirements [MBH+14]. The main deficiencies
of the existing models are in spatial consistency, frequency range, massive-MIMO support,
dual-mobility in moving environment (e.g. V2V communications).
While common channel models such as SCM, WINNER, and IMT-Advanced were designed
for frequencies of up to 6 GHz, there are other models available such as IEEE 802.11ad that
focus on the 60 GHz band. Whereas those models are only applicable for a specific frequency
range, the final METIS channel model shall cover the full frequency range from cellular bands
of below 6 GHz up to 86 GHz, which sets additional challenges to the propagation modelling.
The channel model investigation in the METIS project comprises the analysis of propagation
measurements conducted by METIS partners, extensive literature reviews, and ray tracing
simulations. The objective of the propagation research in METIS project is to ensure the
availability of relevant propagation models, especially for the new scenarios defined in METIS.
The radio channel and pathloss models should have the appropriate level of detail based on
the new scenarios defined in [MET13-D11] and meet the purpose of performance evaluation
of the technology components of respective horizontal topics (HTs) and system concepts in
the other work packages (WPs). Some of the models are derived from propagation
measurements conducted for relevant scenarios.
The key results of the METIS channel model investigation are as follows.
Identified 5G Requirements
Performed channel measurements at various bands between 2 GHz and 60 GHz
Provided channel model considering the requirements
o Map-based model
o Stochastic model for different frequency bands, including frequency-agile
pathloss model
Scalable (hybrid) channel model
o From accurate geometry specific to generic statistical model
o From high complexity to medium complexity
Key technical findings are discussed more in details in this deliverable, especially in Section 4.
The results of the measurement campaigns from the METIS project are provided based on the
propagation scenarios derived from the METIS scenarios and Test Cases in [MET13-D11].
The channel model consists of a stochastic model and a map-based model. The former model
is geometry-based stochastic channel model (GSCM) further developed from WINNER and
IMT-Advanced. The latter is ray-tracing model based on a simplified map. Scalability of the
model means a possibility in which users of the model may combine elements from map-
based model and partly utilise the stochastic model. Channel models are parameterised
based on propagation measurements, literature, and ray tracing simulations.
This deliverable is an updated version from the Initial channel models described in D1.2
[MET14-D12]. The updates are corrections and additions based on the feedback received
from other work-packages during the progress of the METIS project, analysis of measurement
data and literature studies. More specifically, the map-based model is extended to cover
indoor office, open air festival and highway scenarios. The map-based model is also simplified
to reduce the complexity, and partially validated via measurements and simulations. The
stochastic model is extended to provide a low complexity multi-dimensional shadowing map,
millimetre wave parameters (60 GHz), direct sampling of power angular spectrum and
frequency dependent pathloss models. A hybrid model providing scalability between the map-
based model and the stochastic model is also introduced. This report includes additional
measurement results to D1.2.
METIS Public 2
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 3
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
very large arrays and Massive-MIMO the following parameters have to be modelled
accurately: azimuth and elevation angles of paths, angle distributions, distance of the
first/large bounce scatterer for non-planar waves, correlation distance of large-scale
parameters, and polarization.
METIS Public 5
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 6
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
3 Propagation Scenarios
Propagation scenarios (PSs) define types of physical environments and relevant radio
propagation mechanisms that we need to consider. In order to classify the PSs, a number of
basic propagation environments, link types and topologies are defined as follows:
Basic propagation environments
o dense urban
o urban
o rural
o indoor (office and shopping mall)
o highway
Link types
o cellular base station to user equipment (BS-EU)
o backhaul (BH) BS-BS
o device-to-device (D2D)
Link topologies
o outdoor to outdoor (O2O)
o outdoor to indoor (O2I)
o indoor to indoor (I2I)
In addition to the basic propagation environments there are also more exceptional
environments such as stadiums, outdoor festivals and disaster environments, e.g., due to
earthquakes. For each propagation scenario, the following parameters defined:
Cell type, e.g., micro-cell, macro-cell, indoor, outdoor-to-indoor, and so on,
BS location(s) and antenna height(s) in relation to the local roof-top heights (outdoor),
and corresponding indoor parameters in relation to the building dimensions,
UE antenna height and location range(s) within the environment and UE velocity
range(s).
The requirements for the PSs are based on those for the corresponding METIS test cases
(TC1 – TC12) stated in [MET13-D11], which define the baseline propagation scenarios for
simulations. The detailed mapping between PSs and TCs are presented in Appendix B. The
usage of the test cases and PSs has been expressed in detail in [MET13-D61]. The text in this
section has been aligned with these documents.
METIS Public 7
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 8
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
4 Literature Review
This section, firstly, discusses available channel models in the literature with respect to the 5G
channel model requirements as described in Section 2. Available models like WINNER / IMT-
Advanced, COST 2100, and IEEE 802.11 are first discussed. Then the shortcomings of the
available models are emphasized.
In the METIS channel model, higher frequency bands up to millimetre wave should be
considered. Therefore, a general overview on the frequency dependency of propagation
parameters is given in this section.
The modelling of the vehicles and human beings as non-stationary (moving) objects are
needed for obtaining accurate channel modelling, especially in open air festival, stadium and
dense urban scenarios, which are summarized in Appendix C.2.
METIS Public 9
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 10
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 11
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
3GPP D2D
WINNER II
Advanced
3GPP 3D
802.11ad
IEEE
IMT-
Feature
stochastic map-based
Frequency
1–3 1–6 .45 – 6 1–4 1–4 60 – 66 up to 70 GHz up to 100 GHz
Range (GHz)
Bandwidth 100 MHz < 6 GHz, 10 % of the
5 100 100 100 100 2000
(MHz) 1 GHz @ 60 GHz centre frequency
Support
no limited no no limited yes no yes
massive-MIMO
Support
no no no no no no no yes
spherical waves
Support
extremely large
arrays beyond no no no no no no no yes
stationarity
interval
Support dual
no no no limited no no limited yes
mobility
Support Mesh
no no no no no no no yes
networks
Support 3D
no yes no no yes yes yes yes
(elevation)
Support mmW no no no no no yes partly yes
Dynamic very
no no no no limited no yes
modelling limited
METIS Public 12
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Spatial
no no no no no no SF only yes
consistency
The loss through a wall has been discussed in [ITUR13-2040] and consists of the
transmission through the air – wall surface, the attenuation due to the finite conductivity of the
wall material and the transmission through the wall – air surface on the opposite side. In
addition it is assumed that the effect of surface roughness in both of the surfaces affects
similarly as in the reflection.
4.4.3 Effect of surface roughness on reflection
Surface roughness is inversely proportional to the wavelength. This means that the effective
surface roughness is proportional to the carrier frequency. There is a critical height ℎ𝑐 that
divides the surface to smooth (ℎ < ℎ𝑐 ) and rough (ℎ > ℎ𝑐 ) when ℎ is the deviation of the
highest and lowest values of the surface height. In general, so-called Rayleigh criterion is
considered as the value of the critical height [LFR96]. In addition, when the height deviations
are assumed to follow the normal distribution (Gaussian rough surface), surface roughness
attenuation factor that is used to multiply the result for smooth surface reflection has been
clarified in [LFR96]. In [LFR93] and [LFR96] reflection from three different walls consisting of
different materials, brick wall, limestone and metal coated glass, were measured and
compared with theoretical values.
4.4.4 Diffraction from the edge/ wedge
The diffraction phenomenon can be understood as a shielding of the Fresnel zone, thus it is
frequency dependent. In [MH13] there is an approximation for the attenuation caused by one
knife-edge. It is used in the METIS map-based channel model and the equation can be found
in Section 6.2 Equation (6-6). Calculation of diffraction from the wedge is compared with
measurement at 5 GHz in ITU-R [ITUR13-2040] where the calculations of uniform geometrical
theory of diffraction (UTD) wedge diffraction have good matching with measurement results.
Material characteristics can be included in UTD calculation. Map-based model utilizes the
capability of UTD calculation with frequency dependent material properties. The details of
implementation are described in Section 6 and Appendix C.
METIS Public 13
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 14
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 5-1: Summary of major findings from METIS measurements and simulations.
Findings Measurements Relevance
5.1 ITU-R Urban Microcellular pathloss model fits our LOS
Explained in the
measurements beyond 6 GHz but not necessarily with our NLOS Stochastic model,
section
measurements. Section 7, Table
5.2 A modified ITU-R UMi LOS pathloss model covers the Appendices A.1, 7-1
frequency range between 0.8 to 60 GHz. A.3
5.3 Scattered field patterns from a periodic surface show a Diffuse scattering
critical frequency and observation distance where the scattering N/A in the map-based
effects start to be visible. model, Section 6.
5.4 Spatial channel model parameters available for dense-
Appendix A.2
urban microcellular and device-to-device scenarios.
5.5 Supplemental parameters are derived for the 3GPP spatial
channel model in urban macro cellular scenarios at 2.3 GHz. Stochastic model
5.6 Elevation spread at the base station becomes smaller Section 7,
when elevating the base station in urban microcellular Appendix A.1 Table 7-2.
environments.
5.7 3GPP spatial channel model parameters for inter-vehicular
channels available at 2.3 and 5.25 GHz.
5.8 Simultaneous measurements of 10 and 60 GHz outdoor
Stochastic model
channels reveal consistent delays of dominant propagation paths, Appendix A.6 Section 7
while less multipaths were observed at 60 GHz.
5.9 Diffraction can be a dominant propagation phenomenon in
NLOS indoor 60 GHz propagation. Map-based
5.10 Both specular and diffuse components are important in model, Section 6.
Appendix A.4
60 GHz indoor channel.
5.11 The WINNER-type sub-path distribution model in a cluster Stochastic model
is not suitable for 5G channel models. Section 7.
Stochastic model
5.12 WINNER II channel model parameters are available for
Appendix A.5 Section 7, Table
short-range millimetre wave links in various environments.
7-1 and Table 7-2.
5.13 Electrically very large array antenna makes the multi-user Massive MIMO
[MCN14]
channels as ideal as the i.i.d. condition. link design
5.14 Multiple-stream spatial transmission is possible at Millimetre-wave
[HGW13]
millimetre wave due to available spatial degrees-of-freedom. MIMO link design
5.1 ITU-R UMi pathloss model fits our LOS measurements beyond 6 GHz but
not necessarily with our NLOS measurements
The first finding concerns the applicability of the ITU-R pathloss model to wider range of radio
frequency beyond the present definition, i.e., 6 GHz. A multi-frequency pathloss measurement
was conducted in Hatchobori, Tokyo. Ten storey (approximately 40 m high) buildings surround
METIS Public 15
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
this area along the streets. A transmitter (TX) antenna was installed on a car roof with an
elastic pole to adjust the height while a receiver (RX) antenna was fixed on the roof of another
car. The received instantaneous power level was measured while driving the car with the RX
antenna along the routes of LOS, NLOS1, 2 and 3 as shown in Figure 5-1 (left). Frequencies
of 0.8, 2.2, 4.7, 8.45, 27 and 36 GHz and TX/ RX antenna heights of 10 m and 2.5 m with
omni-directional pattern in horizontal plane were used for the measurement. Figure 5-1 (right)
shows pictures of the measurement site where tall buildings surround the street.
32 m 15 m 32 m
Hatchobori- 10 m 5m 10 m
ni chome
NLOS 1
NLOS 2
NLOS 3
intersection Kayabacho
intersection
LOS Tx 32 m
12 m 10 m
242 m
57 m 169 m
Figure 5-1: Measurement environment and routes (left) and picture of LOS route from TX
antenna (right).
The measured data was post-processed in the following manner: (1) Reference points were
established every meter. (2) The data at the reference point and other data recorded within
5 m before/after the reference point were considered as a data set for the reference point. (3)
The median of the data set was determined as the measurement results at the reference point.
60 60 60
80 80 80
Path loss [dB]
Path loss [dB]
METIS Public 16
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
In Figure 5-2, the measurement and M.2135 [ITUR09-2135] results are shown for each
frequency. M.2135 shows decent matching with measurement in LOS route, even though the
model was originally developed for below 6 GHz frequencies, while large deviations are
observed in NLOS routes at all frequencies. As an example, the predicted pathloss is too high,
whereas for NLOS 2 route the prediction is too low for NLOS 1 and 3 routes. The result
supports the use of the M.2135 UMi LOS model even for the higher radio frequencies than
6 GHz (see Table 7-1 for a list of recommended pathloss model in METIS).
5.2 A modified ITU-R UMi pathloss model covers the frequency range between
0.8 to 60 GHz
Inspired by the applicability of the ITU-R M.2135 UMi pathloss model [ITUR09-2135] to a
diverse range of radio frequencies summarized in Section 5.1, the model is further tested with
additional channel measurement sets detailed in Appendices A.1, A.3 and A.6 that cover a
radio frequency range from 0.8 to 60 GHz obtained from three cities. The nominal gain of the
TX and RX antennas was subtracted from the measured pathloss. For the LOS scenario,
comparison of the model with the extended channel data sets suggests slight modifications of
the pathloss model with respect to the following aspects.
Break point distance: the breakpoint distances were found to be much smaller than that
specified in the M.2135 model in many cases, and therefore, a breakpoint scaling factor
𝛼BP is introduced; 𝛼BP is a function of the radio frequency as illustrated in Figure 5-3 (left)
and given by
𝑓𝑐
log 10 ( )
1 GHz
𝛼BP = 0.37 exp (− 1.13
). (5-1)
It must be noted, however, that (5-1) is only valid for elevated base stations above 5 m.
Figure 5-3 (left) shows that the V2V measurements did not follow this trend and hence
was not taken into account in deriving (5-1). In case of the V2V scenario, the BP scaling
factor was 7.5 and 1.3 at 2.3 and 5.25 GHz, respectively.
Pathloss offset: adding an offset, 𝑃𝐿1 , to the model improved the overall agreement with
measurements. Since the original M.2135 model is close to the free space pathloss, the
initial pathloss reflects the effect of surrounding scattering environments. The offset is
given by
𝑐 𝑓
𝑃𝐿1|dB = −1.38 log10 (1 GHz) + 3.34. (5-3)
as shown in Figure 5-3 (right). The offset tends to decrease as the frequency increases,
showing that the higher frequency tends to follow the free space loss better.
Other parameters of the model such as power decay constant before and after the breakpoint,
𝑛1 = 2.2 and 𝑛2 = 4.0, and the shadow fading 𝜎S = 3.1 are found almost frequency
independent for the tested sets of the measurements. The modified pathloss formula is
defined as:
𝑑1 𝑐 𝑓
𝑃𝐿(𝑑1 )|dB = 10𝑛1 log10 (1 m ) + 28.0 + 20 log10 (1 GHz) + 𝑃𝐿1|dB (5-4)
METIS Public 17
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
𝑑1 ℎ ℎ 𝑓
𝑃𝐿(𝑑1 )|dB = 10 𝑛2 log10 (1 m ) + 7.8 − 18 log10 (1 BS
m1m
UT 𝑐
) + 2 log10 (1 GHz ) + 𝑃𝐿(𝑑BP )|dB (5-5)
for 𝑑𝐵𝑃 < 𝑑1 < 500 m where (5-4) and (5-5) represent before and after the breakpoint,
respectively. It is important to note the valid TX-RX distance range since the METIS
measurements do not necessarily cover the entire distance range of the original M.2135
model, leading to the valid range only up to 500 m. The last term 𝑃𝐿(𝑑BP ) in (5-5) is derived
from (5-4) by substituting 𝑑BP into 𝑑 to ensure consistency of the pathloss at the breakpoint.
Exemplary comparison of the modified pathloss model with measurements is shown in Figure
5-4. The fitting at 26.36 GHz shows significant underestimation of the pathloss for the TX-RX
distance shorter than 40 m. This is because of the model assuming close to free space while
the measurements having a large difference in TX and RX antenna heights, leading to a
mismatch in antenna directivities.
Figure 5-3: Frequency dependency of break point (BP) scaling factor (left) and pathloss offset
(right).
We based our NLOS pathloss model on the Manhattan grid layout given by
𝑑2 𝑓c
𝑃𝐿|dB = 𝑃𝐿LOS (𝑑1 )|dB + 17.9 − 12.5𝑛𝑗 + 10𝑛𝑗 log10 ( )+ 3log10 ( )+ 𝑃𝐿2|dB (5-6)
1m 1 GHz
where 𝑑1 and 𝑑2 are the distances from the BS to a cross section of streets and from the UE
to the cross section. The BS and UE are located in the two crossing streets; 𝑃𝐿LOS (𝑑1 ) is the
pathloss between the BS and the cross section derived from (5-4) and (5-5), and the last term
𝑃𝐿2 is the pathloss offset. The model (5-6) is simplified relative to the original model; while the
original model calculates the pathloss from BS to UE and from UE to BS, and takes the
smaller of them [ITUR09-2135], the present model calculates the pathloss from BS to UE only.
Though the pathloss between BS and UE must be the same regardless of the direction of
signal transmission due to reciprocity, the pathloss models do not necessarily hold the
reciprocity. It was found that the simplified model works as good as the original model
according our available data sets. Figure 5-5 shows exemplary fits of the measured pathloss
with the simplified NLOS model. The right figure shows significant overestimation of the
measured pathloss, suggesting a room for improvement of the power decay constant.
METIS Public 18
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 5-4: Example fits of the modified M.2135 UMi LOS model (5-4) and (5-5) at 26.36 GHz (left)
and 60 GHz (right); red curves are model, black dots are measurements.
Figure 5-5: Example fits of the modified M.2135 UMi NLOS model (5-6) and (5-7) at 26.36 GHz
(left) and 37.07 GHz (right). They were measured in two different streets with 𝒅𝟏 being 248 m
(left) and 59.5 m (right); red curves are model, black dots are measurements.
Figure 5-6: Frequency dependency of the pathloss offset (left) and shadow fading (right).
Finally, Figure 5-6 shows the estimated pathloss offset 𝑃𝐿2 and shadow fading 𝜎S . Different
marks represent estimates from different streets. The figures reveal that the pathloss offset
and shadow fading does not change over frequencies, while shows notable dependence on
streets, reflecting the fact that the coupling from the main to perpendicular streets depends
highly on the artefact and vegetation at the street crossing. It is therefore possible to model
these parameters as random variables representing different streets. The mean and standard
deviation of the pathloss offset are -9.1 and 6.1 dB, while those of the shadow fading are 3.0
and 1.3 dB, respectively. A normal distribution may work as a distribution to reproduce the
parameter values, but needs more datasets to ensure statistical validity.
METIS Public 19
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
5.3 Scattered field patterns from a periodic surface show a critical frequency
and observation distance where the scattering effects start to be visible
The literature review in Section 3 shows inadequacy in the knowledge of frequency
dependence in scattering. To fill the knowledge gap, frequency dependency of scattering
effect from rough surface is investigated. Since the electrical size of scatterer becomes large
in high frequency bands, the effect of surface roughness may impact scattering phenomena
which may be diffusing rather than specular reflection. Figure 5-7 (top left) shows the analysis
model. The overall size of a plate is 10 m by 10 m, surface roughness is considered as
periodical structure where 1 m by 1 m tiles are arranged with (grey tile) and without (white tile)
a height offset Δℎ, which appears periodically.
Figure 5-7 (right) shows the calculation results when normal plane wave incidence of 𝜃0 =
0 deg and scattering 𝜃𝑠 = 0 deg. The horizontal axis represents frequency. The vertical axis
represents distance D from the origin in Figure 5-7 (bottom left). Note the distance D is a
variable identical to the distance along z axis from the origin. The colour represents scattered
field strength from rough surface with 𝛥ℎ = 0.1 m. Here, the scattered field is normalized by
the scattered field strength from smooth surface aligned in xy plane (equivalently all tiles with
Δℎ = 0 𝑚) and is illustrated in relative level. The physical optics is used for the calculation. 𝑅1
and 𝑅2 are radii of the first and second Fresnel zones, respectively. The radius of the 𝑛-th
Fresnel zone, 𝑅𝑛 , is defined by 𝑅𝑛 = √𝑛𝜆𝐷 (here 𝜆 is wave length). That is, the Fresnel radius
becomes large as distance D increases. Figure 5-7 (right) illustrates that relative level is very
small when 𝑅2 < 0.5 m. This means that the size (or diameter) of the second Fresnel zone is
smaller than the size of a roughness element (a 1 m by 1 m tile in our model). The scattering
behaves similarly to specular reflection.
Furthermore, frequency 𝑓0 is illustrated at the relative level of 3 dB when 𝐷 = 0.1 m. The
scattering can be considered as specular reflection for lower frequencies than 𝑓0 . In other
words, the frequency 𝑓0 indicates a boundary that the surface with 𝛥ℎ is considered as rough
surface or not.
In general, when 𝜃0 = 0 deg and 𝜃𝑠 = 0 deg, the surface roughness is described as [MEE+14]:
𝜆 𝐶
𝛥ℎ < 𝛼 (= 𝛼𝑓) (5-8)
METIS Public 20
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
r1=0.5
r1=5
r2=05
r2=5
y 100
x R1=5m
Distance, D (m)
10m
R2=5m
1m Top view 1
y
h
1m x R1=0.5m
D R2=0.5m
q0 qs 0.1
Incident wave 0.1 1 10
Observation point f0 Frequency (GHz)
z
Figure 5-7: Analysis model (top left), frequency dependency with h = 0.1 m (bottom left), and
scattering effect on rough surface (right).
5.4 Spatial channel model parameters available for dense UMi and D2D
scenarios
Spatial channel models need to be complemented as new test cases and propagation
scenarios emerge. To this end, wideband radio channel measurements were performed in one
of the most crowded areas in the vicinities of Shibuya railway station in Japan, for METIS’
dense urban test case (TC2). Here, in order to analyse influence of shadowing of pedestrians,
the measurements were performed in daytime and midnight-time as shown in Figure 5-8. In
the investigated area, there were few pedestrians in midnight-time. The measurements were
conducted by DOCOMO channel sounder [KSO+09] at 2.225 GHz centre frequency for urban
microcell (UMi) and D2D scenarios. The user equipment (UE) antenna height was set to 1.5 m
and the base station (BS) antenna heights for UMi and D2D scenarios were set to 3 and 1.5 m,
respectively. A sleeve antenna and a slotted cylinder antenna were used to transmit vertically
and horizontally polarized wave, respectively. Here, these antennas were manually switched.
The signal was received by cylindrical array antenna, which has 96 dual-polarized patch
antenna elements (192 feeds) that allows directional analysis. The distribution both at UE and
BS side were measured by placing the cylindrical receive array antenna on the intended side.
In data processing, space alternating generalized expectation-maximization (SAGE) algorithm
was applied to extraction of paths on the angular-delay-polarization domain. Details of the
measurement campaign are described in Appendix A.2.2.
In both UMi and D2D scenarios, the average received power in daytime is about 5 dB lower
than that in midnight-time (here, TX and RX polarizations are vertical). In UMi scenario in
daytime, the median cross polarization ratio (XPR) from V-pol. to H-pol. was 8 dB, and the
median XPR from H-pol. to V-pol. was 5 dB. The other obtained results are shown in Table
5-2, which are usable directly in a spatial stochastic channel model detailed in Section 7. The
angular spread values are shown in a log-scale in the table because of the compatibility with
the 3GPP and WINNER models. The values in a linear-scale are summarized in Appendix
A.2.3.
Figure 5-8: Measurement environment during day- (left) and night-time (right).
METIS Public 21
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 5-2: Measurement results in the Tokyo downtown crowded area for UMi and D2D
scenarios.
5.5 Supplemental parameters are derived for the 3GPP spatial channel model
in UMa scenarios at 2.3 GHz
Spatial channel model parameters for three-dimensional radio propagation are complemented
through extensive outdoor macro cellular measurements. The measurements were conducted
for outdoor to outdoor (O2O) and outdoor to indoor (O2I) scenario by EB PropSound channel
sounder [Ele04] at 2.3 GHz in downtown Oulu, Finland. The height of TX antenna was 18 m.
The RX antenna was installed on the rooftop of car at the height of 2.5 m in O2O
measurements. In O2I measurements, the RX antenna was located 1.6 m above floor level on
the different floors of Hotel Scandic Oulu. The antenna configurations were 30 (TX) x 16 (RX)
and 30 x 56 for O2O and O2I measurements, respectively. Details of the measurement
campaigns are summarized in Appendix A.1. Table 5-3 summarizes new parameters obtained
from UMa O2O measurements for the 3GPP channel model framework. These parameters
are seen as a supplement to the already available parameters in [3GPP14-36873] intended
for the spatial channel model detailed in Section 7.
METIS Public 22
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Partly insufficient SNR, and therefore the loss of multipath components, causes smaller
angular statistic in the NLOS case than the LOS case.
Figure 5-9 presents the received power for the measurement over the building rooftop in O2I
scenario. The TX antenna had two different positions. The results are overlaid by ray-based
prediction that considers the knife-edge diffraction at the building rooftop and the transmission
loss caused by a window.
METIS Public 23
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 5-9: Received power in UMa O2I at 2.3 GHz; a prediction by a ray-tracing model is
overlaid.
5.6 Elevation spread at the base station becomes smaller when elevating the
base station in UMi environments
Three-dimensional spatial channel model also requires insights into the dependence of
channel model parameters on the height of base station antennas. To clarify the height
dependency, especially with respect to the elevation angular spread at the base station, urban
microcell (UMi) measurements were performed for O2O and O2I scenarios at 2.3 GHz and
5.25 GHz at downtown Oulu, Finland. The antenna configurations were the same as in UMa
scenarios at 2.3 GHz as described in the previous section. For 5.25 GHz, the antenna
configurations were 30 (TX) x 18 (RX) and 30 x 50 for O2O and O2I measurements,
respectively. Details of the measurements are available in Appendix A.1.
Figure 5-10 (left) presents the cumulative density function (CDF) for the ESD for O2O
measurements at 2.3 GHz. The difference between the LOS and the NLOS scenario is clearly
seen from Figure 5-10 (left). The contribution of direct path in the LOS propagation condition is
significantly higher that of multipath components, which leads to smaller ESD in the LOS
scenario. In the NLOS scenario, the ESD differs remarkably between the two TX antenna
heights. One reason for this is the fact that the stronger MPCs are detected from the ground
reflection with the lower TX antenna height. Figure 5-10 (right) presents CDF for the ESD for
UMi O2I corridor measurements at 2.3 GHz and 5.25 GHz. The heights of TX antenna were
5 m, 10 m and 15 m. The RX was located on the different floors of Hotel Scandic Oulu and
several measurement spots were recorded in each floor. It can be seen that ESD decreases,
when the TX antenna height increases. Partly insufficient dynamic range in the measurement,
and therefore a loss of multipath components, causes smaller ESD at 5.25 GHz than at
2.3 GHz.
METIS Public 24
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 5-10: ESD UMi O2O at 2.3 GHz (left), UMi O2I corridor at 2.3 GHz and 5.25 GHz (right).
METIS Public 25
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The measurement campaign was conducted in Kreuzberg, Berlin, Germany. This is a typical
residential and commercial area. The streets are limited by 5 to 6 story buildings to both sides,
thus forming a street canyon. The transmitter was placed on the sidewalk at a height of 5 m
above ground. The receiver was mounted on a mobile cart at a height of 1.5 m. Figure 5-11
(left) shows a map of the scenario and the positions P1 to P4 where the transmitter and
receiver were placed during the campaign.
Measurements were performed with obstruction-free line-of-sight (LOS) between the
transmitter and the receiver as well as with blocked LOS (NLOS). For the LOS measurements
the transmitter was placed at P1 and the receiver at P2 at a distance of ca. 28 m. Another
LOS measurement was performed with the transmitter at P3 and the receiver at P1 at a
distance of ca. 142 m. A NLOS measurement was performed around the corner of a house
with the transmitter at P1 and the receiver at P4 at a distance of ca. 28 m. Details of the
measurement campaign is available in Appendix A.6.
Comparing the average power delay profile (APDP) magnitudes between the two frequencies
can be misleading due to the difference in the effective antenna aperture. When considering
the free space pathloss (according to the Friis transmission equation), one obtains a
difference in pathloss of 15.6 dB for a given reference distance.
Figure 5-11: Location map (left), APDP in dBm for LOS measurement from P1 to P2 (right).
Figure 5-12: APDPs in dBm for LOS measurement from P3 to P1 (left) and NLOS measurement
from P1 to P4 (right).
Comparing the LOS measurements reveals some similarities between the measurements.
The peaks of the LOS multipath component (MPC) in Figure 5-11 (right) at 122 ns differ by
15.5 dB between 60 GHz and 10 GHz. In Figure 5-12 (left) the difference is 15.7 dB. This is in
line with the predicted difference. The MPCs with a greater delay are being caused by
reflections from buildings, cars and other objects. Some of these components are clearly
METIS Public 26
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
visible at both frequencies. In both measurements, however, overall there appear to be less or
weaker multipath components at 60 GHz than at 10 GHz, partly due to smaller dynamic range
of the 60 GHz measurements.
The NLOS scenario in Figure 5-12 (right) exhibits a larger difference between the two
frequencies. While some MPCs are clearly visible at 10 GHz and 60 GHz, some others are
only present at one frequency. In this case as well there seem to be less or weaker MPCs at
60 GHz than at 10 GHz.
Figure 5-13: Measurement scenario (left) and signal strength (right) relative to free space at 1 m
distance for isotropic antennas measured and modelled at 2.4 GHz (blue dots) and 60 GHz (red
dots) in a corridor of an indoor office scenario.
5.10 Both specular and diffuse components are important in 60 GHz indoor
channel
In Figure 5-14 are shown results from a directional channel measurement and modelling using
an extremely large array antenna. The array antenna was realized using a virtual array
antenna principle, where a single antenna element was placed at 25 different positions on the
x-, y-, and z-axes to form a cubic array collectively consisting of 253 = 15625 antenna
positions. The measurement was performed in LOS condition in a large office room with 1.5 m
distance between the transmitter and receiver. The results are based on the same
measurement set-up as described in Section 5.9 and Appendix A.4.2. Full space angle
directional spectrum is provided by Fourier transformation from space to direction domain.
Interesting findings are that there are some distinct spikes in the power delay profile (PDP) for
which the directional spectrum also consists of single spikes. This is interpreted as specular
reflections. Between the spikes of the PDP the directional spectrum is rich and substantially
spread out suggesting that the scattering is due to a multitude of smaller objects or rough
surfaces. Moreover, the directional spread (according to the definition by Fleury as described
in Appendix A.4.2) with respect to the vertical z-axis and one horizontal y-axis decays faster
METIS Public 27
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
than with respect to the horizontal x-axis. The x-dimension is the largest and the z-dimension
the smallest of the room which may explain the different decay times. More details are
provided in Appendix A.4.2.
Figure 5-14: Directional spread and angular spectra shown for delays indicated with numbers
1-5.
METIS Public 28
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 5-15: Real measured [MAB+12] path directional distribution compared with WINNER type
(left) of distribution. The photograph (right) shows the measured paths (circles) as seen from
the base station location.
Further analysis of corresponding MIMO channel has been performed as shown in Figure
5-16. Here the super resolved path estimates from [MAB+12] and the WINNER channel
model, using the same angular spread and a single cluster, is used for determination of the
corresponding MIMO channel singular value distributions. Two different array sizes, 30x30 cm
and 4x4 m, with half wavelength element spacing have been simulated at 2 GHz carrier
frequency. The resulting MIMO channel eigenvalue distributions are similar for the small array
antenna size. However, for the larger array a drastically better performing MIMO channel is
provided by the WINNER channel model. This is an effect of that each single sub-path is
resolved by the large array antenna which has a dramatic impact on the MIMO channel
characteristics. Furthermore, this result clearly demonstrates that the WINNER channel model
is not realistic for 5G in the case of large array antennas.
Figure 5-16: Power ordered MIMO channel eigenvalues for two different array antenna sizes
based on the WINNER model and the measured data.
5.12 WINNER II channel model parameters are available for short-range 60-GHz
links in various environments
Present spatial channel model lacks parameters for higher frequencies than 6 GHz. To fill the
gap, millimetre wave spatio-temporal radio channel measurements have been performed in
the 60 GHz band in an indoor shopping mall, an indoor cafeteria, and on an outdoor open
square, that correspond to TC3 "Indoor shopping mall", TC1 "Indoor office", and TC2 "Dense
METIS Public 29
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
urban". A full WINNER II [WIN208-D112] parameter table is derived for shopping mall LOS
and obstructed-LOS (OLOS) due to pillars, cafeteria LOS, and square LOS and OLOS due to
lampposts and human bodies. Additionally, the elevation parameters for the WINNER+
[WIN+10-D53] are provided for the cafeteria and for the square. The parameter tables are
given in Table 5-5 and Table 5-6. The parameterisation is based on channel measurements
with a 4 GHz bandwidth which gives good delay resolution and allows us to detect the
propagation paths directly from the measurement results. Furthermore, it is assumed in the
parameterisation that the each propagation path is considered as a single cluster because we
did not find clustering effects of propagation paths. For the cafeteria and the square, the point
cloud field prediction method [JH14], calibrated with the measurement results, is used to
generate the channel data used in the parameterisation. Detailed procedure of the
measurements and parameterisation is given in Appendix A.5. The parameters are fully
compatible and hence usable with the METIS stochastic model detailed in Section 7 to
reproduce 60 GHz double-directional wideband channels.
Table 5-5: WINNER parameterisation for 60 GHz range;
parameters affecting validity range of the derived parameters.
Shopping mall Cafeteria Square
LOS OLOS LOS LOS OLOS
min 1.5 4.0 1.0 6.4 6.4
BS-MS distance [m]
max 13.4 16.1 13.1 36.3 36.3
BS 2 2 2 6 6
Antenna height [m]
MS 2 2 1 1 1
Bandwidth [GHz] max 4 4 4
Centre frequency [GHz] 63 63 63
Dynamic Range [dB] 20 20 20
Table 5-6: WINNER parameterisation for 60 GHz range; channel model parameters.
Shopping mall Cafeteria Square
Symbol LOS OLOS LOS LOS OLOS
2D 3D 3D
A 18.4 3.59 15.4 20.3 26.2
PL =A log10(d/1 m) + B
B 68.8 94.3 67.1 67.5 70.5
𝜇lgDS -8.28 -7.78 -8.24 -8.82 -7.72
lgDS = log10(DS/1 ns)
𝜎lgDS 0.32 0.10 0.18 0.37 0.32
𝜇lgASD 1.09 1.61 1.63 1.10 1.49
lgASD = log10(ASD/1 deg)
𝜎lgASD 0.43 0.11 0.25 0.75 0.35
𝜇lgASA 1.19 1.62 1.56 0.24 1.31
lgASA = log10(ASA/1 deg)
𝜎lgASA 0.47 0.14 0.19 0.54 0.44
𝜇lgESD N/A N/A 1.31 0.43 0.74
lgESD = log10(ESD/1 deg)
𝜎lgESD N/A N/A 0.16 0.29 0.32
𝜇lgESA N/A N/A 1.28 0.77 0.95
lgESA = log10(ESA/1 deg)
𝜎lgESA N/A N/A 0.23 0.96 1.19
SF in dB 𝜎SF 1.2 2.1 0.9 0.3 3.5
𝜇KF 7.9 -2.5 8.4
KF in dB
𝜎KF 5.8 2.4 2.2
ASD[°] vs DS[s] 0.4 0.5 0 -0.1 0.4
ASA[°] vs DS[s] 0.2 0.3 0.6 0.3 0.3
ASA[°] vs SF[dB] 0 0.1 0.2 -0.2 0.1
ASD[°] vs SF[dB] 0 -0.1 -0.3 -0.6 0
DS[s] vs SF[dB] 0.2 -0.4 0.5 0.1 0.5
Cross-Correlations
METIS Public 30
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
5.13 Electrically very large array antenna makes the multi-user channels as
ideal as the i.i.d. condition
Finally, we present two link design guidelines that we learnt from METIS channel
measurements. The first guideline described in this section concerns user separability in
massive antenna channels [MCN14; SCP14]. A channel sounding campaign at 5.8 GHz was
performed to investigate the impact of the size of a massive array in a large indoor
environment having high ceiling. The massive array comprises 64 antennas arranged in three
different shapes of different size: 1) a linear array with very large aperture (6 m), 2) a linear
array with large aperture (2 m), and 3) 2D compact array (25cm x 28cm). Eight user terminals
equipped with two antennas are considered. Both LOS and NLOS channels were covered,
and furthermore, various separation distance and orientations of the user terminals were
considered as shown in Figure 5-17.
Figure 5-17: Channel sounding with massive array antenna at the base station: array antennas
at the base station (top left), user terminal antenna (bottom left), and measurement campaigns
(right).
Figure 5-18 (top) shows channel correlation properties for inter- and intra-device settings. The
inter-device correlation (numbered 1 to 8) was much higher for the compact array and for
METIS Public 31
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
aligned devices in the elevation domain. The intra-device correlation (indexed as a and b)
shows that the very large array resolved the two antennas in the same user terminal, while
that was harder in the large and compact arrays. Figure 5-18 (bottom left) depicts the sum of
normalized singular values to discuss inter-link orthogonality. The impact of antenna aperture
size is mainly visible when the users are closely grouped. Even if the number of users
increases, the very large array is able to hold performance as close as the i.i.d. channel.
User
User
b b b
5a 5a 5a 0.5
b b b
6 ab 6a 6a
b b
7a 7a 7a
b b b
8a 8a 8a
b b b
ab ab a b a b a b a b a b ab a b a b ab a b a b a b a b a b a b a b a b ab a b a b a b a b 0
1 2 3 4 5 6 7 8 1 2 3 4 5 6 7 8 1 2 3 4 5 6 7 8
User User User
Figure 5-18: Correlation properties of closely spaced devices under LOS condition; no user
proximity (top) and sum of normalized singular values with respect to the # of users (bottom
left). User and antenna numbering (bottom right).
Figure 5-19 (top) illustrates the condition number of various single-user MIMO channels to see
potential of capacity increase through spatial eigenmode transmission. The condition number
remains low for the very large array thanks to its excellent spatial channel resolvability.
Furthermore, our measurements found that the condition numbers are largely affected by the
user hand grip. Finally, Figure 5-19 (bottom) shows channel power variations across the array.
The largest variations are seen in the very large array, suggesting that different parts of the
array antenna see a different propagation environment, notably different shadowing. In
contrast, even for the compact array, we can see power variations by greater than 10 dB.
METIS Public 32
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
22 22
29 29
-65
36 36
43 43
50 50 -70
57 57
64 a b a b a b a b a b a b a b a b 64 a b a b a b a b a b a b a b a b
1 2 3 4 5 6 7 8 1 2 3 4 5 6 7 8
User User
Figure 5-19: CDF of the condition number for single-user MIMO channels (top) and channel
power variation over base station antenna elements and users (bottom).
channel matrix that exceed a noise level [HGW13]. The noise level is determined to be 5, 10,
15, or 20 dB below the strongest singular value. Figure 5-20 (left) shows the SDoF of a single-
polarized 60 GHz radio channel measured in a conference room environment. Identical
antenna aperture size is assumed for the TX and RX sides. The SDoF is plotted for four
different SNR levels ranging from 5 to 20 dB with 5 dB steps. The result shows that increasing
the antenna aperture size leads to larger SDoF of the channel, more apparently in NLOS
scenarios.
Figure 5-20 (right) illustrates the SDoF of LOS channels with varying communication distance,
showing that more than two SDoFs are available if the distance is beyond 2 m and the SNR is
higher than 15 dB. The results demonstrate that the mm-wave radio channel offers multiple
SDoFs both in LOS and NLOS scenarios such that the spatial multi-stream transmission can
improve the channel capacity, provided that the SNR is sufficiently high to utilize them.
With -10 dBm of the transmit power and 2 GHz of signal bandwidth, Figure 5-21 plots the
capacity for the single-stream gain focusing and multi-stream Eigenmode transmission for
measured LOS channels labelled “BF” and “MUX”, respectively. The figure indicates that the
multi-stream eigenmode transmission outperforms the conventional gain focusing when high
receive SNR is guaranteed, i.e., the antenna aperture size is as large as 9λ2. When the link
budget is limited due to electrically small antennas and/or long TX-RX distances, the
conventional single-stream gain focusing approximates the capacity of the multi-stream
transmission as shown in Figure 5-21 (right).
20
5dB Solid lines: LOS
10dB Dashed lines: NLOS
15 15dB
20dB
SDoF
10
0
0 2 4 6 8
2
Antenna aperture size [l ]
Figure 5-20: SDoF of millimetre wave indoor channels: variation of the SDoF over different
antenna aperture size (left) and variation over different communication distances (right).
Figure 5-21: Capacity for the single- and multi-stream spatial transmission for LOS channels
2 2
with 1λ (left) and 9λ (right) TX and RX antenna apertures.
METIS Public 34
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
6 Map-based Model
6.1 Introduction
The map-based model is intended for cases where accurate and realistic spatial channel
properties are required, for example when studying massive MIMO and advanced beam-
forming techniques. It is also suitable for realistic modelling of pathloss in the case of D2D and
V2V. The model is based on ray tracing using a simplified 3D geometric description of the
propagation environment. The significant propagation mechanisms i.e. diffraction, specular
reflection, diffuse scattering, blocking are accounted for. Building walls are modelled as
rectangular surfaces with specific electromagnetic material properties. The complexity is
scalable as different components, like specular paths and diffuse scattering, may be turned on
or off. A simplified diffraction modelling is also provided in order to further reduce the
complexity. It should be noted that the map based model is significantly less complex than
corresponding stochastic modelling in terms of inter-parameter correlations (i.e. correlation
between path-loss, angle distributions, delay distributions etc.).
The model parameters are given in Table 6-1 below and are explained in the following step by
step description.
Table 6-1: Parameter table for the map based model.
Symbol Propagation scenario
Parameter
[unit] PS1 PS2 PS4 PS5 PS6 PS7 PS8
2 0.1 (traffic 0.1 (traffic
Object density D [1/m ] 0.1 0.05 0.05 4 4
jam) jam)
Object height h [m] 1.5/4 1.5/4 1.5 1.5 1.5/4 1.5 1.5
Object width w [m] 0.5/3 0.5/3 0.5 0.5 3 0.5 0.5
Scatterer absorption
𝛼 0 0 0 0 0 0 0
coefficient
Specular/ diffuse power ratio 𝛽 0.5 0.5 0.5 0.5 0.5 0.5 0.5
Angle dependency factor 𝑞𝜆 [m−1 ] 0.031 0.031 0.031 0.031 0.031 - 0.031
Angle dependency exponent v 3.5 3.5 3.5 3.5 3.5 - 3.5
Angle dependency factor
𝛾 0.95 0.95 0.95 0.95 0.95 - 0.95
(HH)
METIS Public 35
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
1. 2. 3. 4. 5. 6.
Define map Draw random Define point Define Tx and Rx Determine pathways: Calculate path
objects source distribution locations interaction types and lengths and arrival&
for diffuse coordinate points departure
scattering directions
xyz coordinates of interaction points, interaction type per path segment, arrival & departure wave vectors, path lengths
11.
scattering
where 𝑋𝑛 ~𝑈(0, 𝑋𝑚𝑎𝑥 ), 𝑌𝑛 ~𝑈(0, 𝑌𝑚𝑎𝑥 ), 𝑍𝑛 = 1, 𝑛 = 1, … , (𝑋𝑚𝑎𝑥 𝑌𝑚𝑎𝑥 𝐷), 𝑋𝑚𝑎𝑥 and 𝑌𝑚𝑎𝑥 are the
edges of the map on x and y axis, respectively. Different distribution densities are used for
streets (Madrid grid), indoor office, shopping mall, stadium and outdoor festival. Prevent
objects being too close, i.e. closer than half object width, to the TX, the RX, walls or to each
other. The shadowing screen sizes (object height & width) correspond to different objects, i.e.
humans, vehicles, trees, lamp posts and so on (the values provided in Table 6-1 correspond
to people and vehicles).
If a moving environment is considered, define the motion of each object. For the case of true
motion, specify the trajectory and the speed for each object. In case of virtual motion, allocate
only a velocity vector for each object in order to model Doppler shift component for a possible
scattering interaction.
Step 3: Define point source distributions for diffuse scattering over planar surfaces i.e. exterior
and indoor walls, floors, ground etc. These distributions should be as dense as required by the
angular resolution (aperture) of the antenna used in the simulations. For low resolution a tile
size of 10 m x 10 m is used and for high resolution a tile size of 5 m x 5 m is used. For each
surface specify xyz coordinates of point sources. They can be either drawn randomly or based
on a regular pattern (tiles) as shown in Figure 6-10.
METIS Public 36
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The centre points of the tiles, for a rectangular wall, can be determined, e.g., as follows. The
number of tile centres in vertical (z-axis) direction and that in horizontal direction are
respectively:
wall height
𝑁𝑧 = ⌈ ⌉
√∆𝑆
{ wall width
(6-2)
𝑁𝑥𝑦 = ⌈ ⌉
√∆𝑆
where ∆𝑆 is the upper limit of tile area. The z-coordinates of tile centres are determined by
dividing the range of z coordinates of the wall to 𝑁𝑧 equal segments and taking the centre of
each segment as a z coordinate of a tile. The same procedure is followed for the x- and y-
coordinates using 𝑁𝑥𝑦 segments. Finally the wall contains 𝑁𝑧 𝑁𝑥𝑦 tiles and tile centres, new tile
wall height∙wall width
area is calculated as ∆𝑆 = 𝑁𝑧 𝑁𝑥𝑦
.
Steps 1-3 are performed only once. After this step the procedure is fully deterministic.
TX and RX locations:
Step 4: Define a single location or a trajectory, in xyz co-ordinates, for each transmitter and
receiver antenna element. The roles of TX and RX are interchangeable, but for simplicity we
use terms TX and RX throughout this description. Multiple radio links can be modelled
consistently to the environment specified by Steps 1-3 by repeating Steps from 4 to 12. To
simplify the notation the following steps are described for a single radio link only.
Define position vectors for receiver antenna elements 𝑢 and transmitter antenna elements 𝑠
𝐫𝑢RX = [𝑥𝑢 𝑦𝑢 𝑧𝑢 ]𝑇
{ (6-3)
𝐫𝑠TX = [𝑥𝑠 𝑦𝑠 𝑧𝑠 ]𝑇
METIS Public 37
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 6-2: Example of possible pathways with diffraction based on Berg’s model and diffraction
based on UTD (right).
The output of Step 5 is a set of parameter vectors 𝛙𝑘 = {𝜓𝑘𝑖 } = {𝑥𝑘𝑖 , 𝑦𝑘𝑖 , 𝑧𝑘𝑖 , 𝑇𝑘𝑖 }, 𝑘 = 1, … , 𝐾,
𝑖 = 1, … , 𝐼𝑘 , where 𝐾 is the number of pathways, 𝐼𝑘 (max 5 recommended for reduced
complexity) is the number of path segments, 𝑥𝑘𝑖 , 𝑦𝑘𝑖 , 𝑧𝑘𝑖 are x, y and z coordinates of 𝑖 th
interaction point of 𝑘th pathway, 𝑇𝑘𝑖 is the interaction type {direct, reflection, diffraction, object
scattering, diffuse scattering}. The last path segment ends always to a terminating node (RX /
TX) and the zeroth node is TX (or RX).
Determination of different interaction types is described here with more details:
1) Determine specular paths according ray optics principles. These paths are equivalent
to the “direct paths” using specular images of the TX and RX nodes. Use single order
reflection of TX and RX. All surfaces seen by either RX or TX will result in respective
mirror images. Combining single order of both TX and RX provides second order
specular paths.
2) Identify building corners for diffraction at vertical edges based on a 2D layout as
demonstrated in Figure 6-2. Corners either in LOS or accessible via single specular
reflection, to both TX and RX node, are taken into account. When a diffracting corner is
identified the z-coordinate of the actual diffraction point can be determined with basic
trigonometry only after the whole pathway from TX to RX is identified. Notice that
different pathways are allowed whether diffraction model of Berg’s recursive model or
uniform theory of diffraction UTD is utilized. In the former only diffraction to the shadow
region is allowed. In the latter diffraction also to the lit region is enabled which typically
results a higher number of pathways. The difference in possible pathways between the
two methods is illustrated in Figure 6-2, where cyan lines represent pathways between
example TX and RX locations in the Madrid layout.
3) Identify building roof top edges for diffraction at horizontal edges if either TX or RX is
above surrounding buildings. Follow here the principles of the vertical-plane-launch
(VPL) method described in [LB98]. The principle is simple if only diffraction and
reflection in vertical direction are accounted for. This is done by launching a vertical
plane via TX and RX locations. Then diffraction point coordinates are obtained by
identifying crossing points of the plane and horizontal edges of roof tops (See Figure
C-17). Reflection points from opposite buildings or the ground can be found utilizing
ray optics principles assuming diffraction points as a source radiating to all directions
within the plane. VPL supports also multiple successive diffraction and reflection both
METIS Public 38
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
in vertical planes and in 2D (horizontal) dimensions. For simplicity this option is left out
from the proposed model.
4) Identify pathways: 1) via scatterers (both on diffusing wall surfaces and other objects)
which are in LOS and near to either TX or RX node, and 2) via scatterers, which are in
LOS to two nodes, which are in NLOS should be accounted for. For case 1) discard
pathways via weak scatterers i.e. those for which
𝑅 𝑑direct
20 log10 ( 2𝑑1 𝑑2
) < −30 dB (6-4)
where 𝑅 is the radius of the scatterer, 𝑑direct is length for the segment of line between
the node before (node 1) the scatterer the node after (node 2) the scatterer, and 𝑑1 and
𝑑2 are the distances for the two path segments from node 1 to the scatterer and from
the scatterer to node 2, respectively.
Step 5:b) Repeat Step 5 a) for each identified secondary node treating them as TX nodes.
Step 5:c) Repeat Step 5 a) for each identified third node treating them as TX nodes, but find
pathways only to the RX location. Pathways not terminating to the RX location are discarded.
Principles of determining pathways for outdoor-to-indoor case are described in Appendix B.
Step 6: Determine arrival and departure directions for each path 𝑘 in form of wave vectors
from the geometry. Wave vector 𝐤 TX 𝑘,𝑢,𝑠 is pointing from TX to the first interaction point and
𝐤 RX
𝑘,𝑢,𝑠 from the last interaction point to RX. Determine the length 𝑑𝑘,𝑖,𝑢,𝑠 of each segment 𝑖 of
each path 𝑘 by calculating the Euclidean distance from the previous interaction point to the 𝑖th
𝐼
interaction point. Determine the total path length 𝑑𝑘,𝑢,𝑠 = ∑𝑖𝑘 𝑑𝑘,𝑖,𝑢,𝑠 and propagation delay
𝜏𝑘,𝑢,𝑠 = 𝑑𝑘,𝑢,𝑠 /𝑐 for each path 𝑘, where 𝑐 is the speed of light.
Determination of propagation channel matrices for path segments:
Notice that when the Berg’s recursive model is utilized Steps 8, 9 and 11 are not performed for
any path 𝑘 containing diffraction (see details in Step 10a).
Step 7: Determine shadowing due to objects for path segments. Each blocking object is
approximated by a rectangular screen as illustrated in Figure 6-3. The screen is vertical and
perpendicularly oriented with respect to the line connecting the two nodes of the link in the
projection from above. This means that as either node is moving the screen turns around a
vertical axis the centred at the screen so that it is always perpendicular to the line connecting
node 1 (before the screen) and node 2 (behind the screen). The corresponding shadowing
loss is modelled using a simple knife edge diffraction model for the four edges of the screen
as
where 𝐹ℎ1 ,𝐹ℎ2 and 𝐹𝑤1 ,𝐹𝑤2 account for knife edge diffraction at the four edges corresponding
to the height, ℎ, and width, 𝑤, of the screen (see Figure 6-3). The shadowing for a single edge
is given by
𝜋 𝜋
atan(± √ (𝐷1 +𝐷2 −𝑟))
2 𝜆
𝐹= 𝜋
(6-6)
where 𝜆 is the wave length, 𝐷1 and 𝐷2 are the projected distances (according to the
projections from side and from above in Figure 6-3) between the nodes and the edges of the
screen and 𝑟 is the projected distance between the nodes. The plus sign refers to the shadow
zone for each projection (i.e. it is possible that one projection is in LOS and the other in
NLOS). When the link is in NLOS the plus sign apply to both edges. For LOS conditions the
METIS Public 39
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
edge farthest from the link is in the shadow zone (plus sign) and the other in the LOS zone
(minus sign) as shown in Figure 6-3.
𝑟𝑤
𝑟ℎ
For a sparse distribution of shadowing objects, e.g. cars in a street, the pathloss due to
multiple screens is simply given by the sum of the individual losses of all screens in dB units.
In case of a dense distribution of screens, e.g. in the open air festival scenario, the
corresponding sum will result in unrealistic excessive loss. For this case the blocking model is
complemented with the Walfisch-Bertoni model [WB88]. First the dominating shadowing
screen in the vicinity of the RX is identified as shown in Figure 6-4. For this screen the
shadowing 𝐿sh1 is determined as described above. Then the additional loss due to multiple
screens is determined. For this purpose all the screens, which are blocking the LOS between
the RX and TX are identified. The additional diffraction loss 𝐿md due to multiple screens is
given by
𝜃 𝑑
𝑔𝑝 = √ (6-8)
𝜃0 𝜆
with 𝜃0 = 1 rad and 𝑑 is the average distance between the screens and 𝜃 ∈ [0, 𝜋/2] is the
elevation angle from the upper edge of the main blocking screen and the TX. The total
resulting shadowing loss is then given by
METIS Public 40
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
TX
RX
In the case where both TX and RX are low in height, as illustrated in Figure 6-5, a modified
model is used. In this case the angle has a lower limit
max(ℎoff ; ℎTX −ℎsc )
𝜃= 𝐷2𝑅𝑋
𝜃0 (6-10)
where
𝑑 −0.9
ℎoff = (2.35 𝜆 ) 𝑑0 (6-11)
and 𝑑0 is the distance between the TX and the closest blocking screen. Again, the loss due to
the screens closest to RX and TX are added to the multiscreen loss, resulting in
Tx
Rx
Figure 6-5: Shadowing of multiple screens when both TX and RX are at low height.
Step 8: For LOS path segments 𝑘, 𝑖 the 2x2 polarization transfer matrix is given by
1 0
𝐡𝑘,𝑖,𝑢,𝑠 = [ ] (6-13)
0 −1
METIS Public 41
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Step 9: Accounting for all path segments 𝑘, 𝑖 with specular reflection the 2x2 propagation
matrix at node 𝑖 is given by
𝐡𝑘,𝑖,𝑢,𝑠 = 𝛽𝐀ref
𝑘,𝑖,𝑢,𝑠 (6-14)
where 𝛽 is the ratio of reflected and scattered power. The 𝛽 parameter can be either taken
from Table 6-1 or alternatively calculated based on the surface roughness effect as defined in
Appendix E.2.4. The divergence factor for reflection is 𝐹 ref = 𝑠𝑟 ⁄(𝑠in + 𝑠𝑟 ) , where 𝑠in is the
distance from the TX to the reflection point, and 𝑠𝑟 is the distance from the reflection point to
the RX node. The 𝐀 matrix is the polarimetric reflection coefficient defined as follows:
In Figure 6-6 point A is a transmitter, B is a reflection point on a surface (plane W), and C is a
receiver. Vector ⃗⃗⃗⃗⃗
𝐴𝐵 is the direction of the incidence ray and vector ⃗⃗⃗⃗⃗
𝐵𝐶 is the direction of the
reflected ray.
eq i eq r
e i e r
e ||
e C
A (rx)
(tx)
⃗⃗⃗⃗⃗ ‖ + ‖𝐵𝐶
ray. The total path length is ‖𝐴𝐵 ⃗⃗⃗⃗⃗ ‖. Angles of arrival and departure, composed of the
azimuth and elevation directions, are given by the vectors ⃗⃗⃗⃗⃗ 𝐴𝐵 and ⃗⃗⃗⃗⃗
𝐶𝐵. The polarimetric
reflection coefficient matrix
𝛼𝜃𝜃 𝛼𝜃𝜙
𝐀 = [𝛼 𝛼𝜙𝜙 ] (6-15)
𝜙𝜃
defines the relation between the incidence and reflected electric field components as follows
𝐸𝜃𝑟 𝐸𝜃𝑖
[𝐸 ] = 𝐀 [𝐸 ] (6-16)
𝜙𝑟 𝜙𝑖
where 𝜃 and 𝜙 denote the polarization components in 𝒆𝜽 -direction and 𝒆𝝓 -direction with
respect to the propagation direction 𝒆𝒓 of the path. Entries of matrix 𝐀 are calculated as
where the unit vectors and Fresnel reflection coefficient are defined as
METIS Public 42
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
𝑥𝑖
𝒆𝑟𝑖 = [𝑦𝑖 ] is the unit vector in the direction of the incidence ray (𝒆𝑟𝑖 ∥ ̅̅̅̅
𝐴𝐵)
𝑧𝑖
𝑥𝑟
𝒆𝑟𝑟 = [𝑦𝑟 ] is the unit vector in the direction of the reflected ray (𝒆𝑟𝑟 ∥ ̅̅̅̅
𝐵𝐶 )
𝑧𝑟
𝑥𝑒𝑛𝑑 − 𝑥𝑏𝑒𝑔𝑖𝑛 0 𝑥𝑒𝑛𝑑 − 𝑥𝑏𝑒𝑔𝑖𝑛 0
𝒏 = [𝑦𝑒𝑛𝑑 − 𝑦𝑏𝑒𝑔𝑖𝑛 ] × [0]⁄‖[𝑦𝑒𝑛𝑑 − 𝑦𝑏𝑒𝑔𝑖𝑛 ] × [0]‖ is the unit normal vector of the
0 1 0 1
surface to direction of the incidence plane, × is the vector cross product
−𝑦𝑖
√𝑥𝑖 2 +𝑦𝑖 2
𝒆𝜙𝑖 = 𝑥𝑖
and 𝒆𝜃𝑖 = 𝒆𝑟𝑖 × 𝒆𝜙𝑖 are unit vectors representing the directions of
√𝑥𝑖 +𝑦𝑖 2
2
[ 0 ]
incidence for polarized electric fields 𝐸𝜙𝑖 and 𝐸𝜃𝑖 at the reflection point B
𝒆 ×(𝒏×𝒆 )
𝒆⊥𝑖 = 𝒆𝑟𝑖 × 𝒆∥𝑖 and 𝒆∥𝑖 = ‖𝒆𝑟𝑖×(𝒏×𝒆𝑟𝑖)‖ are unit vectors of incident rays that have electric
𝑟𝑖 𝑟𝑖
fields perpendicular and parallel to the plane ABC
𝒆 ×(𝒏×𝒆𝑟𝑟 )
𝒆∥𝑟 = ‖𝒆𝑟𝑟 and 𝒆⊥𝑟 = 𝒆𝑟𝑟 × 𝒆∥𝑟 are unit vectors of reflected rays that have
𝑟𝑟 ×(𝒏×𝒆𝑟𝑟 )‖
electric fields parallel and perpendicular to the plane ABC
−𝑦𝑟
√𝑥𝑟 2 +𝑦𝑟 2
𝒆𝜙𝑟 = 𝑥𝑟
and 𝒆𝜃𝑟 = 𝒆𝑟𝑟 × 𝒆𝜑𝑟 are unit vectors representing the directions of
√𝑥𝑟 2 +𝑦𝑟 2
[ 0 ]
reflection for polarized electric fields 𝐸𝜙𝑟 and 𝐸𝜃𝑟 in the reflection point B
𝜖𝑟 sin 𝜓−√𝜖𝑟 −cos2 𝜓
𝑅∥ = denotes the Fresnel reflection coefficients for electric fields
𝜖𝑟 sin 𝜓+√𝜖𝑟 −cos2 𝜓
parallel to the plane ABC, where 𝜓 is the grazing angle, i.e. the angle between the
plane W and the incidence ray and 𝜖𝑟 is the relative permittivity of the plane (See
Table E-1 and Table E-2) [VA03]
Step 10a (Berg’s recursive option): For path ways 𝑘 with diffraction calculate 2x2 propagation
matrices with simplified method called the Berg’s recursive model
𝜃𝜃 𝜃𝜙
𝜃𝜙
𝐿𝜃𝜃
𝑘 e
𝑗Φ𝑘,𝑖
𝐿𝑘 𝑒 𝑗Φ𝑘,𝑖
𝐡dif
𝑘,𝑖,𝑢,𝑠 = [ 𝜙𝜃 𝜙𝜙 ], (6-18)
𝜙𝜃 𝜙𝜙
𝐿𝑘 e𝑗Φ𝑘,𝑖 𝐿𝑘 e𝑗Φ𝑘,𝑖
where terms 𝐿𝑘 are defined below and i.i.d. random phases Φ follow uniform distribution in
range [0,2𝜋]. Notice, if a pathway 𝑘 contains more than one diffraction, the effect of all
diffractions is calculated jointly and set to the first diffraction segment. Furthermore, specular
reflection points are neglected in this step, i.e. the path distance between diffraction nodes is
based on all (if any) intermediate specular reflections. For the successive diffraction segments
Equation (6-18) is substituted by unit matrices (all one entries for diagonal elements,
otherwise zero). Notice, as mentioned in the Step 5a) this simplified method can be utilized
only for diffraction producing paths to the shadow region.
METIS Public 43
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The LOS and diffracted pathways are described by the Berg recursive model [Ber95]. It is
based on the assumption that a street corner appears as a secondary source when a
propagating radio wave turns around it. The corners of buildings and the antennas represent
nodes (See Figure 6-7 left).
Along a propagation path each node contributes a loss which depends on the change in
direction 𝜃. The total loss at a specific node 𝑗 is given by the well-known expression for free
space loss between isotropic antennas where a fictitious distance 𝑑𝑗 is used, i.e.
4𝜋𝑑𝑗
𝐿𝑗|dB = 20 log10 ( 𝜆
) (6-19)
where is the wave length. It should be noted that the fictitious distance corresponds to the
real distance but multiplied by a factor at each diffraction node. The result is that the fictitious
distance 𝑑𝑗 becomes longer than the real distance meaning that it accounts for diffraction loss
when used in the free space loss, cf. Equation (6-19). An example with four nodes is shown in
Figure 6-7 (right). At each node 𝑗, the fictitious distance is given by the following recursive
expression
with
where 𝑠𝑗 is the real distance between node 𝑗 and its following node (𝑗 + 1), 𝑞𝑗 is a function of
𝜃𝑗 (See Equation (6-22) and Figure 6-7 middle). The initial values are 𝑑0 = 0 and 𝑘0 = 1.
TX
y [m]
x [m]
Figure 6-7: Example of a street corner acting as a node (left). Manhattan map (middle).
Topological example with four nodes (right).
The angle dependence for the fictitious propagation distance extension is given by the
following expression
𝜃𝑗 𝑣
𝑞𝑗 = 𝑞90 (90 deg) , (6-22)
where 𝜃𝑗 ∈ [0, 180] deg , 𝑞90 and 𝜈 are parameters determined by fitting the model to
measurement data. The parameter 𝑞90 accounts for the amount off diffraction loss caused by
each node. A larger value results in larger diffraction loss. The corresponding frequency
dependency is given by
𝑞
𝑞90 = √ 𝜆𝜆 (6-23)
METIS Public 44
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
where 𝑞𝜆 is a model parameter given in Table 6-1. The parameter 𝜈 accounts for how fast the
loss changes in the transition zone between LOS and NLOS. Further, there is a corresponding
depolarization matrix:
1 ∞
𝐐90 = 𝑞90 [ ], (6-24)
∞ 𝛾
where 𝑞90 is the corresponding co-polarized (𝜃𝜃) depolarization coefficient and 𝑞90 𝛾 is the co-
polarized (𝜙𝜙) coefficient, 𝑞90 and 𝛾 are both given in Table 6-1 (for rooftop diffraction the
diagonal elements should be swapped).
The Berg recursive model ignores all types of interactions except diffraction, i.e. if a pathway
contains diffraction, no reflection or scattering is calculated in Steps 9 and 11. This is
illustrated by Figure 6-8 where the corresponding distances 𝑠𝑗 for the segments 𝑠𝑔𝑖 , are
𝑠0 = 𝑠𝑔1 , 𝑠1 = 𝑠𝑔2 + 𝑠𝑔3, 𝑠2 = 𝑠𝑔4 and 𝑠3 = 2𝑠𝑔5 𝑠𝑔6⁄𝑅, where 𝑅 is the scatterer radius used
in Equation (6-44). Then the resulting set of segments 𝑠𝑗 , is used in Equation (6-20) for
providing 𝐡dif
𝑘,𝑖,𝑢,𝑠 in Equation (6-18).
Figure 6-8: Example of a path with diffraction, specular reflection and object scattering.
Step 10b (UTD option): Determine a diffractions coefficient matrix 𝐀 with the uniform theory of
diffraction (UTD) formulas and calculate 2x2 propagation matrix as
𝐡dif dif
𝑘,𝑖,𝑢,𝑠 = 𝐀 𝑘,𝑖,𝑢,𝑠 . (6-25)
The divergence factor for diffraction is 𝐹 dif = √𝑠in⁄𝑠D (𝑠in + 𝑠D ), where 𝑠in is the distance from
the TX to the diffraction point, and 𝑠D is the distance from the diffraction point to the RX node.
Calculation of the diffractions coefficient matrix 𝐀 is described in the following. Figure 6-9
illustrates canonical model for wedge diffraction where the spread angle of the wedge is
(2 − 𝑛)𝜋 with the incident ray impinges at skew angle. A unit vector 𝒍 parallel to the direction
of the wedge is also defined. Unit vectors 𝒔in and 𝒔D indicate direction of propagation for the
incident ray and the diffracted ray with respect to diffraction point D, respectively. These unit
vectors are related based on the Fermat principle:
(𝒔in − 𝒔D ) ∙ 𝒍 = 0. (6-26)
METIS Public 45
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Equation (6-26) indicates the extended Snell`s law of 𝛽D = 𝛽in . The diffracted ray should be
defined within a cone of the spread angle 𝛽D with an apex at D: the diffracted ray propagates
with arbitrary angles on the xy-plane (Figure 6-9 right). Hence the unit vector of direction of
propagation for diffracted ray can be determined by
lˆ 2 n x
in “0” 面
y
sin “n”面
in z
in 2 n
n0
“0” 面 D D D
“n” 面
x
z
y sD D
nn
D
Figure 6-9: Model for canonical problem of diffraction.
𝐷 𝐷𝑏
𝐀=[ 𝑎 ] (6-28)
𝐷𝑐 𝐷𝑑
defines relation between the incidence and diffracted electric field components as follows
𝐸𝛽D 𝐸𝛽in
] = 𝐀[ ]. [ (6-29)
𝐸𝜉D 𝐸𝜉in
The unit vectors of electric filed are defined as
𝒍−𝒔in cos 𝛽in
𝒆𝛽in = sin 𝛽in
is unit vector representing the directions of incidence for polarized
electric fields 𝐸𝛽in at diffraction point D
𝒍−𝒔in cos 𝛽in 𝒔 ×𝒍
𝒆𝛽in = sin 𝛽in
and 𝒆𝜉in = sinin𝛽 are unit vectors of polarized electric fields 𝐸𝛽in and
in
𝐸𝜉in of incidence ray at D, respectively,
𝒍−𝒔D cos 𝛽D 𝒔 ×𝒍
𝒆𝛽D = sin 𝛽D
and 𝒆𝜉D = sinD 𝛽 are unit vectors of polarized electric fields 𝐸𝛽D and 𝐸𝜉D
D
of diffracted ray at D, respectively.
There are various expressions of the diffraction coefficient. Especially UTD is well-known
method formulated by R. G. Kouyoumjian and P. H. Pathak [KP74] where the wedge is
considered to be conductor. In this section, one of the practical methods by Luebbers [Lue84]
is described. His expression of diffraction coefficient accounts for finite conductive wedge
which extended the capability of UTD for more practical use. When the wedge is long enough
compared to wavelength, the diffraction coefficient in Figure 6-9 of the canonical problem is
defined as
𝐷𝑎 𝐷𝑏
[ ] = 𝑰2×2 (𝐷+ (𝜉D − 𝜉in ) + 𝐷 − (𝜉D − 𝜉in )) + 𝐑 0 𝐷 − (𝜉D + 𝜉in ) + 𝐑 𝑛 𝐷+ (𝜉D + 𝜉in ) (6-30)
𝐷𝑐 𝐷𝑑
with
METIS Public 46
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
𝜋
− exp(−𝑗 ) 𝜋+Δ𝜉
𝐷 ± (Δ𝜉) = 4
cot ( ) 𝐹 (𝑘𝐿𝑎± (Δ𝜉)), (6-31)
2𝑛√2𝜋𝑘 sin 𝛽in 2𝑛
Here 𝑰2×2 is the 2x2 identity matrix, and 𝐑 0 and 𝐑 𝑛 are matrices of reflection coefficient for
“0”face and “n” face defined in Figure 6-9 respectively as expressed by
𝒆𝛽in ∙ 𝒆⊥𝑟 𝒆𝛽in ∙ 𝒆∥𝑟 𝑅⊥ 0 𝒆⊥𝑖 ∙ 𝒆𝛽in 𝒆⊥𝑖 ∙ 𝒆𝜉in
𝐑 0 = [𝒆 ∙ 𝒆 𝒆𝜉in ∙ 𝒆∥𝑟 ] [ 0 𝑅∥ ] [ 𝒆∥𝑖 ∙ 𝒆𝛽in 𝒆∥𝑖 ∙ 𝒆𝜉in ], (6-32)
𝜉in ⊥𝑟
where definitions of unit vectors 𝒆⊥𝑖 , 𝒆∥𝑖 , 𝒆⊥𝑟 , 𝒆∥𝑟 and reflection coefficients 𝑅⊥ , 𝑅∥ are same as
in Step 9. 𝐹(𝑥) in Equation (6-31) and denotes the Fresnel integral
+∞ 2
𝐹(𝑥) = 2𝑗√𝑥𝑒 𝑗𝑥 ∫√𝑥 𝑒 −𝑗𝜏 𝑑𝜏 (6-34)
2𝑛𝜋𝑁 ± −Δ𝜉
𝑎± (Δ𝜉) = 2 cos2 ( ). (6-36)
2
Here the closest integer value satisfying the equation is chosen. Note, a method for numerical
approximation of the Fresnel integral is described in [ACL12].
In a special case when TX and RX are under Shadow Boundary (i.e. 𝜉D = 𝜉in + 𝜋 ) or
Reflection Boundary (i.e. 𝜉D = 𝜋 − 𝜉in ) condition, the cotangent function in Equation (6-31)
becomes singular. In this case, the term containing cotangent function should be replaced by
𝜋±𝜉 𝑗𝜋 𝑗𝜋
cot ( 2𝑛 ) 𝐹 (𝑘 𝐿 𝑎± (𝜉)) ≈ 𝑛 [√2𝜋 𝑘 𝐿 sgn(𝜀) − 2 𝑘 𝐿 𝜀 exp ( 4 )] exp ( 4 ) (6-38)
with 𝜀 defined by
𝜀 = 𝜋 ± (𝜉 − 2𝑛 𝜋 𝑁 ± ). (6-39)
For further explanations, the terms in (6-31) – (6-37) and their formulations can be also found
in [Bal89].
Step 11: For path 𝑘 with scattering from objects or diffuse scattering point sources calculate
2x2 propagation matrices as follows
𝜃𝜃 𝜃𝜙
exp(𝑗Φ𝑘,𝑖 ) exp (𝑗Φ𝑘,𝑖 )
𝑠𝑐
𝐡𝑘,𝑖,𝑢,𝑠 = √𝐺𝑘,𝑖,𝑢,𝑠 [ 𝜙𝜃 𝜙𝜙
] (6-40)
exp (𝑗Φ𝑘,𝑖 ) exp (𝑗Φ𝑘,𝑖 )
METIS Public 47
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
where random phases Φ follow uniform distribution in range [0, 2𝜋]. The divergence factor for
scattering is 𝐹 sc = 1⁄𝑠𝑠 , where 𝑠𝑠 is distance from the scattering point to the receiver node.
𝑠𝑐_𝑜𝑏𝑗
The scattering gain 𝐺𝑘,𝑖,𝑢,𝑠 of the scattered wave from an object is modelled based on the
scattering cross section for a perfectly conducting sphere as
𝑅𝐶𝑆 = 𝜋𝑅 2 . (6-41)
The received power when substituting the effective aperture 𝐴𝑒 = 𝐺𝑟 𝜆2⁄4𝜋 of the RX antenna
is given by
𝑃 𝐺𝑡 𝐺𝑟 𝜆2 𝜆𝑅 2 1 2
𝑃𝑟 = 𝑆𝑟 𝐴𝑒 = 𝜋𝑅 2 (4𝜋𝑅𝑡 2 = 𝑃𝑡 𝐺𝑡 𝐺𝑟 (8𝜋) (𝑅 ) . (6-44)
1 𝑅2 ) 4𝜋 1 𝑅2
1 2 2
(𝑅 ) = (𝐹 LOS 𝐹 sc ) , (6-45)
1 𝑅2
which are taken into account in Step 12. Now, we can define the scattering gain from Equation
(6-44) as
𝜆𝑅 2
𝐺 𝑠𝑐 = ( ) . (6-46)
8𝜋
Further, it is assumed that the scattered power is also shadowed, i.e. it is reduced according
to the shadowing model (6-5). The corresponding scattering gain is given by
2
𝐺 𝑠𝑐_𝑜𝑏𝑗 = 𝐺 𝑠𝑐 (1 − 𝛼)(1 − (𝐹ℎ1 + 𝐹ℎ2 )(𝐹𝑤1 + 𝐹𝑤2 )) . (6-47)
where 𝛼 is the absorption coefficient of the scatterer given in Table 6-1. In order to match the
size of the screen with the cross section of the sphere the radius is set to
𝑤ℎ
𝑅=√𝜋 . (6-48)
As the expression is reciprocal, the RX and TX nodes may be switched depending on which
node is in the vicinity of the scatterer.
METIS Public 48
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
∆𝑆
R2
R1
Figure 6-10: Schematic drawing of the scattering model (left). Random point sources
approximating a rough surface (right).
Diffuse scattering occurs also from surfaces accounting for the surface roughness. According
to [DGM+04] the received power scattered by a surface may be expressed in a similar way as
for the spherical scatterer, i.e. as
𝜆 2
𝑃𝑟 = 𝑆𝑠𝑐 𝐺𝑟 ∆𝑆(1 − 𝛽) ( ) cos(𝜃𝑖 ) cos(𝜃𝑠 ), (6-49)
4𝜋𝑅2
where ∆𝑆 the size of a fraction of the surface area corresponding to the point source, 𝛽 is the
relative amount of specular power and 𝜃𝑖 and 𝜃𝑠 are the angles of the incoming and scattered
paths relative to the normal of the surface (see Figure 6-10 right). Reformulating (6-49) yields
𝜆2 1 2
𝑃𝑟 = 𝑃𝑡 𝐺𝑡 𝐺𝑟 ∆𝑆(1 − 𝛽) cos(𝜃𝑖 ) cos(𝜃𝑠 ) ( ) (6-50)
64𝜋3 𝑅1 𝑅2
𝜆2
𝐺 𝑠𝑐_𝑜𝑏𝑗 = ∆𝑆(1 − 𝛽) cos(𝜃𝑖 ) cos(𝜃𝑠 ), (6-51)
64𝜋3
The parameter 𝛽 can be either taken from Table 6-1 or alternatively calculated based on the
surface roughness effect as defined in Appendix E.2.4.
If virtually moving scattering objects, e.g. cars, are modelled each object has a path specific
𝑜𝑏𝑗
Doppler shift parameter 𝜔𝑘 in rad/s. The Doppler shift may be either drawn randomly or
calculated based on a velocity and an incidence angle. A path with virtually moving scattering
objects has time dependent channel matrix
𝑠𝑐 𝑜𝑏𝑗
𝐡𝑘,𝑖,𝑢,𝑠 (𝑡) = √𝐺𝑘,𝑖,𝑢,𝑠 exp(𝑗𝑡𝜔𝑘 ). (6-52)
METIS Public 49
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
where 𝐠 RX TX
𝑢 and 𝐠 𝑠 are the complex polarimetric antenna pattern vectors, of RX element 𝑢 and
𝑇𝑘𝑖
TX element 𝑠, for the direction and frequency represented by the wave vectors, and 𝐹𝑘,𝑖,𝑢,𝑠 is
the divergence factor (see [DGM+04]) defined in Steps 8-11 for the corresponding path
segment. When calculating a divergence factor the path length to the interaction point is
determined cumulatively starting from the initial node. With Berg recursive option (Step 10a)
𝑇𝑘𝑖
the divergence factor 𝐹𝑘,𝑖,𝑢,𝑠 ≡ 1 for pathways 𝑘 containing diffraction. In this step the product
operator of 2x2 polarization matrices 𝐡 is defined as the element-wise matrix product. In this
step the product operator of 2x2 polarization matrices 𝐡 is defined as the element-wise matrix
product. Divergence factor (see [DGM+04]) 𝐹 𝑇𝑘𝑖 ≡ 1 with Berg option (Step 10a) for pathways
𝑘 containing diffraction. The time parameter 𝑡 is interchangeable to a parameter indicating TX
and RX locations through their moving speed (except if moving scattering objects are
introduced having time dependent channel coefficients with temporal variation independent on
TX and RX locations).
In case of virtual motion (assuming the UE has a velocity that causes small-scale effects like
Doppler, but is fixed in the large-scale geometry) there is no time evolution of large-scale
parameters and the complex impulse response is given by
RX 𝑇 𝑗2𝜋𝑑𝑘,𝑢,𝑠 𝐼𝑘 𝑇𝑘𝑖
𝐾 𝐠 RX
𝑢 (−𝐤 𝑘,𝑢,𝑠 ) 𝑒 𝜆 (∏𝑖=1 𝐡𝑘,𝑖,𝑢,𝑠 (𝑡)𝐹𝑘,𝑖,𝑢,𝑠 )
𝐻𝑢,𝑠 (𝑡, 𝜏) = ∑𝑘=1 𝐷
, (6-54)
TX
𝐠 TX
𝑠 (𝐤 𝑘,𝑢,𝑠 )𝑒 𝑗𝑡𝜔𝑘
𝛿(𝜏 − 𝜏 𝑘,𝑢,𝑠 )
where 𝜔𝑘𝐷 is the Doppler frequency in rad/s of path 𝑘. If both array antennas are small in
physical size, i.e. plane wave assumption is valid, the equation can be further simplified by
approximation
𝑇 𝑗2𝜋𝑑𝑘 𝐼 𝐷
𝐻𝑢,𝑠 (𝑡, 𝜏) = ∑𝐾 RX RX
𝑘=1 𝐠 𝑢 (−𝐤 𝑘 ) 𝑒 𝜆 (∏𝑖=1
𝑘
𝐡𝑘,𝑖 (𝑡)𝐹 𝑇𝑘𝑖 ) 𝐠 TX TX 𝑗𝑡𝜔𝑘
𝑠 (𝐤 𝑘 )𝑒 𝛿(𝜏 − 𝜏𝑘 ), (6-55)
where propagation matrices are determined in reference to centre points of TX and RX arrays
and therefore antenna indices 𝑢 and 𝑠 can be dropped off from polarization matrices 𝐡.
Notice 1, pay attention on path segment lengths and their use in computing attenuation of
matrices 𝐡 in Steps 8-11 such that the pathloss is never calculated twice per any segment.
Notice 2, in the case of Outdoor-to-indoor each path 𝑘 is attenuated by penetration loss as
described in Appendix B “Modelling O2I”.
6.3 Simplifications
The complexity may be reduced by means of accounting only for diffracted paths, which are
within a lower order of Fresnel zones with respect to the corresponding specular or direct
path. If both the specular and the direct paths are shadowed then the diffracted path should
always be accounted for. This is illustrated in Figure 6-11 for diffraction around a corner of a
building at 2 GHz. In Figure 6-11 a) the diffraction corner should be accounted for, while in
Figure 6-11 b) it should not be accounted for. Figure 6-11 c) shows the effect of accounting for
diffracted paths within two Fresnel zones for UTD and horizontal polarization. It is clear that
removing diffraction paths outside the two most inner Fresnel zones has very minor impact on
the corresponding signal strength and channel response. For details more see [LB98].
METIS Public 50
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
a) b)
TX
TX
RX RX
c)
Relative power [dB]
angle in a)
angle in b)
Planned as well as already performed measurement campaigns are being, and planned to be,
compared with model simulations in order to further validate the model.
METIS Public 51
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
7 Stochastic Model
This section describes the stochastic model, i.e. a geometry-based stochastic channel model
(GSCM) further developed from WINNER/3GPP. Within METIS a multitude of channel models
per propagation scenario was proposed. A coarse classification of the proposed channel
models can be done by their requirement for explicit building structures or scene models
containing objects like furniture or vegetation. In case explicit building/scene models are
available it is possible to calculate propagation paths based on this deterministic environment,
cf. Section 6. Without these building/scene models propagation paths need to be generated
fully stochastically. Besides these two extreme cases, it might be possible to utilize explicit
building/scene models within the stochastic modelling approach in order to get more site-
specific results. This third approach is a hybrid of the map-based and the stochastic model.
Within this section we will refer to it as hybrid model, wherever it is applicable.
For all models it shall be possible to apply various antenna models. The channel models are
therefore independent of the antennas. Section 7.1 gives some basic information about
antenna modelling, especially on how to correctly deal with 3D antenna patterns.
Appendix D gives the steps for the WINNER based channel model approach including new
findings from 3GPP’s study items on 3D channel modelling and D2D, as well as further
investigations in METIS.
The output of the channel model shall be complex MIMO channel impulse responses per user
link depending on the user location. In some evaluations it may be adequate to investigate
only the pathloss and shadowing effect. In these cases, the channel model output is the sum
effect of pathloss and shadowing, which can be determined by PL & SF models of WINNER
based approach or by a total path attenuation of the map-based model, without considering
polarizations or antenna models. For some cases the time evolution of the channel is of
importance, therefore it has to be clarified whether the users are mobile or stationary or if the
drop-based simulation is used, when the stabile users are in virtual motion. In case of mobile
users, user routes need to be defined in terms of time and space.
METIS Public 52
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 7-1: Recommended PL and fading models for each propagation scenario.
Freq.
In/ Link
# PS band Pathloss Model Fading Channel Model
Out type
[GHz]
0.45 - 6 PL, steps and parameters: 3GPP 3D-UMi LOS/NLOS [3GPP14-
BS-UE
36873], except in LOS condition PL model is Section 5.2
0.45 - 6 Parameters:1) For small to medium
sized city: Table 5-2,
LOS: Section 5.2
For large city: Table 5-4
O2O NLOS: WINNER+ B1
Urban These measured parameters need to
D2D [WIN+10-D53] Manhattan
Micro be complemented with parameters
1 /V2V layout with 10 dB extra
loss from [3GPP14-36843] Annex A.
Steps: [3GPP14-36873], except
Doppler as in [3GPP14-36843] Annex
A
O2I BS-UE 0.45 - 6 PL, steps and parameters: 3GPP 3D-UMi O2I [3GPP14-36873]
Squar 50 - 70 Parameters: Table 5-6,
O2O BS-UE Table 5-6
e Steps: [3GPP14-36873]*
0.45 - 6 PL, steps and parameters: 3GPP 3D-UMa LOS/NLOS [3GPP14-
BS-UE 36873]
O2O Supplemental parameters from Table 5-3 can be used**
Urban
2 0.45 - 6 PL, steps and parameters: IMT-Advanced UMa [ITUR09-2135]
Macro
BH
(ℎUE = 10 m, ℎBS = 25 m)
O2I BS-UE 0.45 - 6 PL, steps and parameters: 3GPP 3D-UMa O2I [3GPP14-36873]
BS-UE 0.45 - 6 PL, steps and parameters: IMT-Advanced RMa [ITUR09-2135]
0.45 - 6 WINNER+ C2 [WIN+10- Parameters & steps: WINNER+ C2
D2D D53] [WIN+10-D53], except Doppler and
O2O /V2V (ℎUE = 1.5 m, ℎBS = AOA/ASA as in [3GPP14-36843]
1.5 m) Annex A
Rural
3 0.45 - 6 PL, steps and parameters: IMT-Advanced RMa [ITUR09-2135]
Macro BH
(ℎUE = 10 m, ℎBS = 25 m)
BS-UE 0.45 - 6 PL, parameters & steps: 3GPP 3D-UMa O2I [3GPP14-36873]
O2I 0.45 - 6 Parameters & steps: IMT-Advanced-
D2D WINNER+ C4 [WIN+10-
UMa [ITUR09-2135], except AOD/ASD
/V2V D53], page 75, Table 4-1
have same distributions with AOA/ASD
Indoor 0.45 - 6 PL, steps and parameters: WINNER II A1-rr (room-to-room)
BS-UE
Office [WIN208-D112]
4 I2I
Cafete 50 - 70 Parameters: Table 5-6
BS-UE Table 5-6
ria Steps: [3GPP14-36873]
Indoor BS-UE 0.45 - 6 PL, steps and parameters: WINNER II A1 [WIN208-D112]
Shoppi 50 - 70
5 I2I Parameters: Table 5-6
ng BS-UE Table 5-6
mall Steps: [WIN208-D112]
Highw BS-UE 0.45 - 6 PL, steps and parameters: IMT-Advanced RMa [ITUR09-2135]
6 O2O
ay V2V 0.45 - 6 PL, steps and parameters: Karedal [KTC+09]
BS-UE 0.45 - 6 PL, steps and parameters: IMT-Advance RMa [ITUR09-2135]
Open 0.45 - 6 WINNER+ C2 [WIN+10-
D53] Parameters & steps: WINNER+ C2
Air
7 O2O D2D [WIN+10-D53], except AOD/ASD have
Festiv (ℎUE = 1.5 m, ℎBS = same distributions with AOA/ASD
al 1.5 m)
BH 0.45 - 6 PL, steps and parameters: IMT-Advance RMa [ITUR09-2135]
METIS Public 53
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
For the hybrid model approach there is also another possible option to generate site-specific
LSPs for example by extracting the deterministic shadow fading and generating the remaining
LSPs taking the inter-parameter correlations into account. Figure 7-1 shows an example
utilizing the Madrid map. The left figure shows the path gain in dB that was calculated using a
ray tracing simulation within METIS. The figure in the middle shows the shadow fading in dB,
i.e. the path gain minus the log-distance pathloss. Finally the figure on the right shows the
delay spread in dB, which was generated by using the inter-parameter correlation between
shadow fading and delay spread on a random map following a Gaussian distribution.
METIS Public 54
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
ASD vs ASA -0.30 0.56 -0.20 0.59 -0.50 0.20 -0.30 0.30 0.00 0.40 0.00 0.20
ASD vs KF -0.50 0.30 -0.10 0.10
ASA vs KF -0.30 0.01 -0.50 -0.20
DS vs KF -0.70 0.20 -0.40 -0.40
SF vs KF -0.10 0.90 0.00 0.50
ESD vs SF -0.50 0.20 -0.10 0.10 -0.40 0.30 -0.30 0.00 -0.20 0.00 -0.70 0.20
ESA vs SF 0.00 0.30 -0.70 0.00 -0.80 -0.80 -0.80 -0.40 -0.30 0.00
ESD vs KF -0.50 0.00 -0.20 0.00
ESA vs KF 0.00 0.20
ESD vs DS -0.50 0.00 -0.50 0.10 -0.60 0.50 -0.50 0.30 -0.50 0.00 -0.60 0.20
ESA vs DS -0.20 0.20 -0.20 0.00 -0.40 0.50 -0.20 0.50
METIS Public 55
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
ESD vs ASD 0.40 0.70 0.10 0.60 -0.50 0.50 0.40 0.50 0.20 0.50 -0.20 0.50
ESA vs ASD -0.40 0.50 -0.10 0.50 -0.60 0.20 -0.40 -0.10 0.00 0.10
ESD vs ASA -0.20 0.10 -0.30 0.00 -0.60 0.40 -0.30 0.00 -0.20 0.00
ESA vs ASA 0.00 0.20 -0.20 0.20 0.10 0.70 0.40 0.40 0.50 0.80
ESD vs ESA 0.00 0.10 0.00 0.20 -0.40 0.60 -0.30 0.50
Delay scaling
𝑟𝜏 1.90 3.20 1.70 4.80 1.40 2.20 2.10 2.50 1.90 2.30 1.70 2.20
parameter
𝜇XPR 9.00 29.0 8.00 29.0 6.00 10.0 8.00 12.0 7.00 10.0 9.00 12.0
XPR [dB]
𝜎XPR 3.00 9.00 3.00 10.0 9.00 12.0 4.00 9.00 3.00 10.0 10.0 11.0
Number of clusters 𝑁 4.00 12.0 16.0 25.0 12.0 12.0 8.00 12.0 20.0 20.0 12.0 12.0
Cluster ASD in
𝐶𝐴𝑆𝐷 0.50 5.99 0.50 10.0 5.00 8.00 5.00 6.00 2.00 2.00 5.00 8.00
[deg]
Cluster ASA in
𝐶𝐴𝑆𝐴 0.50 18.0 0.50 22.0 5.00 8.00 11.0 12.0 15.0 15.0 5.00 8.00
[deg]
Cluster ESA in
𝐶𝐸𝑆𝐴 0.50 7.00 0.50 7.00 3.00 3.00 7.00 7.00 7.00 7.00 3.00 3.00
[deg]
Cluster ESD in
𝐶𝐸𝑆𝐷 0.50 3.64 0.50 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
[deg]
Per cluster
𝜁 3.00 7.00 3.00 6.00 3.00 6.00 3.00 6.00 3.00 6.00 4.00 5.00
shadowing std [dB]
DS 1.30 9.00 1.40 10.0 10.0 10.0 30.0 40.0 40.0 40.0 10.0 10.0
ASD 8.00 23.0 3.20 10.0 11.0 11.0 15.0 18.0 50.0 50.0 11.0 15.0
Decorrelation ASA 1.00 12.0 8.10 9.00 17.0 17.0 15.0 15.0 50.0 50.0 17.0 35.0
distance in the
SF 7.80 14.0 3.60 13.0 7.00 7.00 37.0 45.0 50.0 50.0 7.00 7.00
horizontal plane
[m] K 7.90 15.0 12.0 12.0
ESA 12.0 15.0 1.80 10.0 25.0 25.0 15.0 15.0 50.0 50.0 25.0 25.0
ESD 12.0 17.0 0.80 10.0 25.0 25.0 15.0 15.0 50.0 50.0 25.0 25.0
frequency range 0.45-63 GHz 0.45-6 GHz
UE height range 1.5 - 22.5 m 1.5 - 22.5 m
BS height range 3 m - 10 m 5 m - 10 m 5 m - 15 m 18 m - 25 m
METIS Public 56
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 57
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
𝒆𝝓
𝒆𝒓
𝜃
𝑟 𝒆𝜽
𝜙 𝑦
0 𝑇
𝜃 = arccos {(0) 𝒆𝒓 (𝜃, 𝜙)} (7-5)
1
and
1 𝑇
𝜙 = arg {( 𝑗 ) 𝒆𝒓 (𝜃, 𝜙)}. (7-6)
0
7.2.2 Vector field rotation / mechanical tilting
When an antenna is rotated around its principal axes the radiation field is rotated likewise.
Since this vector field rotation is not too trivial we will discuss it in the following. We define a
global coordinate system (GCS) with coordinates (𝜃, 𝜙) and a local coordinate system (LCS)
with coordinates (𝜃′, 𝜙′) . With 𝑹 ∈ ℝ3x3 being an arbitrary rotation between those two
coordinate systems with
0 𝑇
θ′ = arccos {(0) 𝑹𝑻 𝒆𝒓,𝑮𝑪𝑺 (𝜃, 𝜙)}. (7-9)
1
From (7-8) and (7-6) it follows that
1 𝑇
ϕ′ = arg {( 𝑗 ) 𝑹𝑻 𝒆𝒓,𝑮𝑪𝑺 (𝜃, 𝜙)} (7-10)
0
METIS Public 58
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Since the basis vectors 𝒆𝜽 (𝜃, 𝜙) and 𝒆𝝓 (𝜃, 𝜙) are different in both coordinate systems, there
has to be a transformation from LCS to GCS:
1 0
𝑴𝐿𝑂𝑆 = ( ) (7-16)
0 −1
The factor −1 is due to the opposing directions of 𝒆𝝓,𝐭𝐱 and 𝒆𝝓,𝐫𝐱 as can be seen in Figure 7-3.
METIS Public 59
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Due to LOS time delay there is a deterministic phase ΦLOS which needs to be applied,
resulting in the following expression for the LOS component:
𝑇
𝐹𝜃,𝐺𝐶𝑆,rx (𝜃 𝑎 , 𝜙 𝑎 ) +exp(𝑗ΦLOS ) 0 𝐹𝜃,𝐺𝐶𝑆,tx (𝜃 𝑑 , 𝜙 𝑑 )
( ) ( )( ). (7-17)
𝐹𝜙,𝐺𝐶𝑆, rx (𝜃 𝑎 , 𝜙 𝑎 ) 0 −exp(𝑗ΦLOS ) 𝐹𝜙,𝐺𝐶𝑆,tx (𝜃 𝑑 , 𝜙 𝑑 )
Note that this is different to WINNER II, in which a random phase offset between vertical-to-
vertical (VV) and horizontal-to-horizontal (HH) polarization is introduced to the LOS path.
7.2.5 NLOS depolarization
In case of NLOS 𝑴 is the product of various polarization filters and coordinate transformations
according to the incident planes of passed interactions.
𝑇
𝐹𝜃,𝐺𝐶𝑆,rx (𝜃 𝑎 , 𝜙 𝑎 )
𝑚𝜃𝜃 𝑚𝜃𝜙 𝐹𝜃,𝐺𝐶𝑆,tx (𝜃 𝑑 , 𝜙 𝑑 )
( 𝑎 𝑎
) (𝑚 𝑚𝜙𝜙 ) (𝐹 𝑑 𝑑
), (7-18)
𝐹𝜙,𝐺𝐶𝑆,rx (𝜃 , 𝜙 ) 𝜙𝜃
𝜙,𝐺𝐶𝑆,tx (𝜃 , 𝜙 )
with 𝑚𝜃𝜃 , 𝑚𝜃𝜙 , 𝑚𝜙𝜃 , 𝑚𝜙𝜙 being the NLOS depolarization coefficients. These coefficients can
only be determined if there is additional data available about the interactions. Such data
includes electromagnetic properties of the scatterer material, their locations as well as their
plane normals. With the GSCM approach this is not the case, therefore a stochastic approach
is used to randomly generate the depolarization coefficients. With the map based model the
depolarization coefficients can be determined by calculating the reflection and diffraction
coefficients as described in Section 6.2.
7.2.6 Modelling array antennas
METIS Public 60
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The complex impulse response for V2V channel is presented as double-directional and time-
variant channel as the superposition of 𝑁𝑃 paths, i.e., the contributions from the different types
scatterers [Mol04]
𝑁𝑃
ℎ(𝑡, 𝜏) = ∑𝑖=1 𝑎𝑖 𝑒 𝑗𝑘𝑑𝑖(𝑡) 𝛿(𝜏 − 𝜏𝑖 ) × 𝛿(ϕ𝑅 − ϕ𝑅,𝑖 )𝛿(ϕ 𝑇 − ϕ 𝑇,𝑖 )𝑔𝑅 (ϕ𝑅 ) 𝑔𝑇 (ϕ 𝑇 ), (7-19)
where 𝑎𝑖 is the complex amplitude associated with path 𝑖, 𝜏𝑖 , ϕ𝑅,𝑖 and ϕ 𝑇,𝑖 are the excess
delay, azimuth angle of arrival (AOA), and azimuth angle of departure (AOD) of path 𝑖, 𝑔𝑇 (ϕ 𝑇 )
and 𝑔𝑅 (ϕ𝑅 ) are the TX and the RX antenna patterns, 𝑒 𝑗𝑘𝑑𝑖(𝑡) is the corresponding distance-
induced phase shift and 𝑘 = 2𝜋/𝜆 is the wave number.
The distinctions between discrete scatters are made by examining their Doppler shifts. The
MD causes large Doppler shift whereas SD causes small Doppler shift. The impulse
responses of diffuse scatterers are derived from the measurement data by subtracting the
LOS component and the discrete scatters.
The impulse response of (7-19) can be divided into four parts: LOS, discrete components
originating from interaction with mobile objects, discrete components from interaction with
static objects, and diffuse scattering. Thus, the complex impulse response for V2V channel
can presented as (omitting the AOA and AOD notation for convenience)
where 𝑃 is the number of mobile discrete scatterers, 𝑄 is the number of static discrete
scatterers and 𝑅 is the number of diffuse scatterers. The complex path amplitude for each
discrete scatterer 𝑝 (MD or SD) is modelled as
METIS Public 61
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
𝑛𝑝
𝑗𝜙𝑝 𝑑ref 2
𝑎𝑝 (𝑑𝑝 ) = 𝑔𝑠,𝑝 𝑒 √𝐺0,𝑝 ( 𝑑 ) , (7-21)
𝑝
where 𝑑𝑝 is the sum of distances from the transmitting vehicle to scatterer and from the
receiving vehicle to scatterer, 𝐺0,𝑝 is the received power at the reference distance 𝑑𝑟𝑒𝑓 , 𝑛𝑝 is
the pathloss exponent and 𝑔𝑠,𝑝 is the real-valued, slowly varying, stochastic amplitude gain of
the scatterer. For each discrete scatterer is assigned its own values for 𝑛𝑝 , 𝐺0,𝑝 . For the LOS
path, the model for complex amplitude is the same as for discrete scatterers with the
exception that subindex 𝑝 is replaced by LOS. The complex path amplitude of diffuse scatter 𝑟
is modelled as
𝑛DI
𝑑ref 2
𝑎𝑟 = √𝐺0,DI 𝑐𝑟 ( 𝑑 ) , (7-22)
𝑟
The densities of scatterers are given by 𝜒MD , 𝜒SD and 𝜒DI designating the number of
scatterers per meter. Assuming the moving direction on the x-axis, the y-coordinate of mobile
discrete scatterers are drawn from a uniform discrete probability density function (PDF) where
the possible number of outcomes equals the number of road lanes, 𝑁𝑙𝑎𝑛𝑒𝑠 . The initial x–
coordinates are modelled by a (continuous) uniform distribution over the length of the road
strip, i.e., 𝑥𝑝,0 ~𝑈[𝑥min , 𝑥max ]. Each mobile scatterer is assigned a constant velocity along the
x–axis given by a truncated Gaussian distribution.
The x–coordinates of static discrete scatterers and diffuse scatterers are also modelled
through 𝑥𝑞 ~𝑈[𝑥min , 𝑥max ] and 𝑥𝑟 ~𝑈[𝑥min , 𝑥max ]. The parameters for V2V channel model are
presented in [KTC+09].
Table 7-4: V2V channel model parameters for highway and rural scenario.
Scenario Parameter LOS MD SD DI
G0 [dB] -5 −89 + 24n 104
n 1.8 U [0, 3.5] 5.4
σ
μ 6.8 9.4 6.3
μc [m] 7.2 5.4 4.9
𝑑𝑐𝑚𝑖𝑛 [m] 4.4 1.1 1.0
Highway 𝜒 [1/m] 0.005 0.005 1
y1 [m] -13.5 -13.5
y2 [m] 13.5 13.5
WDI [m] 5
WROAD [m] 18
NLANES 4
G0 [dB] -9 −89 + 24n 23
n 1.6 U [0, 3.5] 3.0
σ
μ 11.7 15.1 14.8
μc [m] 8 8.3 2.5
𝑑𝑐𝑚𝑖𝑛 [m] 5.4 2.5 1.4
Rural 𝜒 [1/m] 0.001 0.05
y1 [m] -9.5 -9.5
y2 [m] 9.5 9.5
WDI [m] 5
WROAD [m] 8
NLANES 2
METIS Public 62
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
2
2𝜎𝑆𝐹
𝑆𝐹 = √ ∑𝐾 ̅ ̅
𝑘=1 sin(𝐷 ∙ 𝛽𝑘 + 𝜃𝑘 ), (7-23)
𝐾
where K is the number of waves, 𝜎𝑆𝐹 is the targeted standard deviation of SF in decibels,
̅ = [𝑥𝑢 , 𝑦𝑢 , 𝑧𝑢 , 𝑥𝑣 , 𝑦𝑣 , 𝑧𝑣 ] is the location vector of the TX/RX pair in 6D space, 𝑢 and 𝑣 indices
𝐷
refer to the TX and RX number, respectively, 𝛽𝑘̅ is the wave vector of the kth wave, 𝜃𝑘 is a
random initial phase in the range [0,2𝜋]. Directions of wave vectors 𝛽𝑘̅ are drawn randomly
from uniform distribution. In practice this is done by drawing randomly each of the six
elements of vector 𝛽𝑘̅ from a uniform distribution 𝒰(−1, 1). The norm of all K wave vectors is
scaled to
1
2𝜋 √2 arg(𝐽0 (𝑥)=𝑒)
‖𝛽̅ ‖ = 𝑑
, (7-24)
𝑐𝑜𝑟
where 𝑑𝑐𝑜𝑟 is the target correlation distance in meters and 𝐽0 (𝑥) is the zeroth order Bessel
function of the first kind. For arg(𝐽0 (𝑥) = 1𝑒), where 𝐽0 (𝑥) is the 1D auto-correlation function, the
value 𝑥 can be found numerically.
The method is efficient in computational complexity and especially in memory consumption.
Only 𝐾(6 + 1) real numbers have to be stored in the memory. An example output in 3D space
(instead of 6D for visualization purposes only) is illustrated in Figure 7-4. In principle all the
other spatial correlated large-scale parameters could be generated with the method as well.
METIS Public 63
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 7-4: Example of 3D shadowing map, colour indicates shadowing value in dB.
METIS Public 64
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
After fixing the physical locations, drifting of LS and SS parameters is enabled for a short
distance movement as illustrated in Figure 7-6. Implementation of the drifting is
straightforward and is fully based on the geometry (for each impulse response, phase, delay,
and angle of arrival is recalculated).
movement
Figure 7-6: Propagation parameter drifting due to small movement of the UE.
METIS Public 65
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 66
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The used stochastic models are discussed in more detail in section 7 and Table 7-1. The map-
based model listed above is discussed in section 6. The hybrid model can be applied in the
same propagation scenarios and frequency ranges as the corresponding stochastic model.
If the target PS can be modelled with more than one option the model user should make the
selection with his/her own judgement. The model choice depends also on the frequency range
on scope of simulations, and the desired accuracy versus the simulation time. The simulation
time depends mostly on the scope of the simulated system like the number of BSs and UEs.
The scope of the simulations determines whether, e.g. spherical waves are needed (with
massive MIMO) or spatially consistent angular properties between users (MU-MIMO) are
required, etc. These matters will be discussed below.
One aspect needs to be mentioned: in case of accurate link level simulations with large array
antennas in a specific radio environment the map-based model can be recommended. The
propagation paths, their directions and powers can be calculated efficiently taking into account
wall materials, windows etc., which is important especially when using high carrier
frequencies. Fast fading properties can also be obtained by coherent combining of
propagation paths at each user antenna position. This is very important when evaluating array
antenna receive/transmit algorithms, for example in 3D beam forming case.
Main differences in map-based and stochastic models can be listed as follows:
- The map-based model relies on the geometrical description of the environment
including e.g. wall materials with roughness and ray tracing calculations. Stochastic
model is based on parameter distributions extracted from the measurements. Map-
based model gives a possibility to model cases that cannot easily be covered by
measurements, e.g. very low signal strengths or signals in places hard to reach like
towers. On the other hand, the measurements can extract phenomena which the map-
based method cannot capture, like actually almost all real environments consisting of
plenty of details practically impossible to be modelled accurately.
- Frequency range of the available parameterization for the stochastic model is typically
below 6 GHz due to limitations of measurement equipment. However, the
METIS Public 67
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
parameterization for 60 GHz is available for the indoor environments and an outdoor
square.
- Frequency range of the map-based model can be extended up to 100 GHz, but care
should be taken when applying the material parameters, especially the surface
roughness.
- Complexity is the main limitation of the map-based model. Some methods to alleviate
this are:
o Calculate only the path-loss for the users, if adequate.
o Remove weak paths in the calculations (e.g. take 20 dB dynamic range instead
of the default 30 dB).
o Remove the surrounding cars and pedestrians from the simulations.
o Disregard the movement, as in the Stadium Test Case.
o Consider pre-calculating and storing the channel coefficients in a memory and
reading these in the simulator when needed.
To ease the selection of the model there is a comparison between the METIS models in Table
8-2.
METIS Public 68
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
it is recommended to carefully understand the requirements set by the simulations where the
model is to be utilized and to choose a model option (stochastic/map-based) with an adequate
level of accuracy among the choices defined in Table 8-1. In the following we discuss whether
simplifications could be used in the following example cases:
Massive deployment of sensors/actuators
link-level simulation with a precise route
drop concept and signal level evaluations
In the following we discuss these aspects in more detail.
8.2.1 Massive deployment of sensors/actuators
Antenna locations can be specified independently for each element in a coordinate system as
described in Step 4 of Section 6.2. As a consequence the model determines propagation
paths and path parameters, such as path delays, directions and shadowing, independently for
each TX and RX element pair (s,u). This capability is beneficial in the case of very large array
antennas, because antenna elements are not necessarily within a stationary interval of the
electric field and wide sense stationary (WSS) assumption of the fading coefficients may not
hold between elements. Thus when simulating very large arrays it is recommended to specify
an explicit position vector for each receiver antenna element, 𝑢, and transmitter antenna
element 𝑠. The drawback is the complexity of determining propagation paths. It will increase
𝑢 ∙ 𝑠 times higher compared to doing it only for the array phase centres (measurement
centres). Obviously antenna radiation patterns have to be defined (obtained from simulations
or measurements) per element phase centres (measurement centres).
8.2.2 Link-level simulation with precise route
Channel realizations for a trajectory of TX, RX or both are constructed following the principle
of true motion as described above. Step 12 of Section 6.2 is utilized for this purpose.
Trajectories have to be sampled dense enough such that the Nyquist sampling criterion of the
fading process is fulfilled. If the simulation system requires higher sampling rate an
interpolation in time domain is performed. In order to avoid tracking of propagation paths in
time domain it may be beneficial to represent channel coefficient in frequency domain, instead
of delay domain, for the interpolation.
8.2.3 Drop concept and signal level evaluations
For system level simulations with a high number of radio links it is typically not necessary to
model precisely any motion or trajectories. In this case the option of virtual motion described in
Section 6.2 is followed. TX and RX locations are drawn randomly or they are specified
according to a deployment scenario. To reduce complexity Berg’s recursive model is used for
diffraction instead of UTD, as described in Step 10 in Section 6.2.
If only attenuation (i.e. pathloss and shadowing without considering antenna effects) for each
radio link is needed, it can be calculated simply. Follow Steps 1 to 11. In Step 12 choose
Equation (6-55) and substitute it with
𝐼 2
𝐻𝑢,𝑠 = ∑𝐾
𝑘=1|∏𝑖=1 𝐡𝑘,𝑖,𝑢,𝑠 |
𝑘
(8-1)
If antenna effects and fast fading are needed then Equation (6-54) should be utilized.
8.2.4 Parameter ranges and default values
Parameter values used in generating example channel model outputs in Section 9 are
specified in Table 9-6. Outputs are illustrated in Sections 9.1 and statistics of simulated
scenarios marked with * are listed in Table 9-7.
METIS Public 69
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 70
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Pairing matrix A is in the example case below a 6x3 matrix with values n,m {0,1}. Value 0
stands for link BS cell n to UE m is not modelled and value 1 for link is modelled.
χc1,ms1 ⋯ χc1,ms3
𝑨=[ ⋮ ⋱ ⋮ ] (8-2)
χc6,ms1 ⋯ χc6,ms3
The pairing matrix can be applied to select which radio links will be generated and which will
not.
METIS Public 71
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
ms1M
c1
ms21
ms12
ms21
c2
ms11
ms2N
Here, M locations of mobile station 1, and N locations of mobile station 2 are defined yielding
a total of M+N points or labels. The resulting procedure is as follows.
1. Set BS c1 and c2 locations and array orientations according to layout.
2. Set MS locations ms11 to ms2N and array orientations according to layout.
3. Set the links to be modelled to 1 in the pairing matrix.
4. Generate all the radio links at once to obtain correct correlation properties. It is
possible to generate more channel realizations, i.e. time samples, for each channel
segment afterwards. This can be done by applying the same values of small-scale
parameters and restoring final phases of the rays.
5. Simulate channel segments in parallel or consecutively according to the desired
motion of the mobiles.
Typically, the links between the MSs and the links between the BSs are not of interest.
Cellular systems are traditionally networks where all traffic goes through one or more BS. The
BS themselves again only talk to a BS hub and not between them.
Multihop and relaying connections break with this limitation as well as D2D connections. In
multihop networks, the data can take a route over one or more successive MS. Relaying
networks, on the other hand, employ another level of network stations, the relays, which
depending on the specific network, might offer more or less functionality to distribute traffic
intelligently. The WINNER channel model can be used to obtain the channels for multihop or
relaying scenarios, as described below.
METIS Public 72
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
MS3 BS3
(BS4)
MS2
MS1
BS2
BS1
(MS4)
X
Figure 8-2: Multihop and relaying scenarios.
In the example figure above the signal from MS2 to BS3 is transmitted via MS3 and BS2 acts
as a repeater for BS1. These scenarios can be generated by introducing a BS-MS pair into
position of a single BS serving as a relay or into position of a single MS serving as a multihop
repeater. The resulting procedure is as follows:
1. Set base station BS1 to BS3 locations and array orientations according to layout.
2. Set mobile locations MS1 to MS3 and array orientations according to layout.
3. Add extra base station BS4 to position of MS3 and extra mobile MS4 to position of
BS2 with same array orientations and array characteristics as MS3 and BS2
respectively.
4. Set the BS x MS pairing matrix to
0 0 0 1
1 0 0 0
A
0 0 1 0
0 1 0 0
5. Generate all the radio links at once.
6. Simulate the channel segments in parallel.
METIS Public 73
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 74
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
9 Simulation Results
Outputs of METIS map-based model are illustrated in this section with example plots and
propagation parameter statistics. The purpose is to give insight to the model output statistics,
like pathloss exponents, delay spreads etc., in a set of example cases. Model outputs are also
compared to measurement results in selected cases. Layouts of building and street maps are
not modified according to the measured environments. Only parameters like antenna heights
and frequencies are aligned with measurement settings. Thus a perfect match is not looked
for in the comparison. The channel realizations are generated with an implementation of the
METIS map-based model using the Berg recursive model for diffraction, except in the
transition scenario of Figure 9-6 where UTD model for diffraction was used.
Figure 9-1: Modelled pathloss in Madrid D2D scenario at 5 GHz (left). Layout with TX locations
denoted with blue circles and RX locations denoted with green dots (right).
METIS Public 75
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 9-2: Comparison of modelled and measured V2V LOS pathloss data measured by Oulu
University at 5.25 GHz (left) and at 2.3 GHz (right).
Figure 9-3: Layout of dense urban micro cell with TX locations denoted with blue circles and RX
locations denoted with green dots (right).
METIS Public 76
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 9-4: Modelled pathloss in Madrid micro scenario with a TX height of 10 m at 5 GHz (left)
and at 26.45 GHz (right).
Figure 9-5: Measured and modelled pathloss at 26.45 GHz in NLOS condition (left) and LOS
condition (right).
Cross correlations of large-scale parameters of LOS links are shown in Table 9-1:
Table 9-1: Cross-correlations of LS parameters (UMi LOS).
DS ASA ASD SF K
DS 1 0.735 0.5351 0.1809 -0.1017
ASA 0.735 1 0.5603 0.231 0.0127
ASD 0.5351 0.5603 1 0.3264 -0.002
SF 0.1809 0.231 0.3264 1 0.0057
K -0.1017 0.0127 -0.002 0.0057 1
Cross correlations of large-scale parameters of NLOS links are shown in Table 9-2:
METIS Public 77
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
An example of a dynamic scenario is shown in Figure 9-6 (left). It shows the TX (red triangle),
RX (green dots) and scattering/shadowing objects (blue dots). The centre frequency is 2 GHz,
the TX height is 15 m and the RX height 1.6 m. The RX is moving along a uniformly sampled
linear route starting from the south and travelling towards north. The time evolution of path
length and angular parameters is depicted in Figure 9-6 (right), Figure 9-7 and Figure 9-8. The
number of paths increases in and around the line-of-sight area, mostly because of diffuse
scattering from the surrounding walls. Arrival angles evolve smoothly along the route, while
departure angles are concentrated to the two street opening directions of the TX location.
Simulation parameters are defined in Table 9-6 in the column “Madrid micro-cell (transition)”.
Figure 9-6: Transition scenario layout (left) and resulting path length with the colour indicating
the path gain in dB.
METIS Public 78
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure 9-7: Resulting AOA (left) and AOA (right) for the transition scenario with the colour
indicating the path gain in dB.
Figure 9-8: Resulting EOA (left) and EOD (right) for the transition scenario with the colour
indicating the path gain in dB.
METIS Public 79
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 9-3: Parameter table for the specific LOS/NLOS simulations with the map-based model.
Parameter Symbol PS1 PS1
[unit] Pedestrian D2D
Map model ‘Madrid’ ‘Madrid’
AP antenna height [m] 5/10 1.5
User antenna
[m] 1.5 1.5
height
Carrier frequency [GHz] 2.6/6/30/60 2.6/6/30/60
2
Object density D [1/m ] 0.05 0.05
Object height h [m] 1.5 1.5
Object width w [m] 0.5 0.5
Scatterer absoption
𝛼 0 0
coefficient
Specular/ diffuse
𝛽 0.5 0.5
power ratio
Angle dependency
𝑞90 0.5 0.5
factor
Angle dependency
v 3.5 3.5
exponent
Angle dependency
𝛾 0.95 0.95
factor (HH)
Diffraction model Berg Berg
Figure 9-9: Basic simulation test scenario for typical pedestrian user scenario with AP
connections and for D2D connections in an urban microcell (Madrid grid).
The simulation results are shown in Figure 9-10, which represents a reference case with
empty streets, and in Figure 9-11 and Figure 9-12 with shadowing pedestrians with a density
of 5 persons per 10x10 m2 area. Blue curves, referring to antenna height of 5 m, and red
METIS Public 80
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
curves referring to antenna height of 10 m, are practically the same. The corner effect is
sharper the higher the carrier frequency is. In Figure 9-11 and Figure 9-12 the PL values vary
roughly 5-15 dB between different distributions of obstructing/scattering objects. The average
pathloss with AP height of 10 m is approximately 1- 6 dB smaller than with AP height of 5 m.
60 GHz 30 GHz
Figure 9-10: Reference pathloss plots between the access point and RX positions with no
obstructing/shadowing objects for AP antenna heights of 5 m (blue) and 10 m (red).
METIS Public 81
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
60 GHz 30 GHz
Figure 9-11: Pathloss between access points and RX positions in the basic pedestrian scenario
for AP antenna heights of 5 m (blue) and 10 m (red).
60 GHz 30 GHz
Figure 9-12: Mean pathloss between access points and RX positions in the basic pedestrian
scenario for AP antenna heights of 5 m (blue) and 10 m (red).
METIS Public 82
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
60 GHz 30 GHz
Figure 9-13: Pathloss between the reference point (AP position in Figure 9-9) and RX positions
in the D2D user scenario.
METIS Public 83
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
60 GHz 30 GHz
Figure 9-14: Mean Pathloss between the reference point (AP position in Figure 9-9) and RX
positions in the D2D user scenario.
It is noted that due to the limited statistics (only 10 random shadowing object configurations
are covered) the simulation results have some degree of variation and are therefore mainly
indicative examples.
METIS Public 84
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 9-4: Excess pathloss due to obstructing/shadowing objects in Madrid grid scenario.
Microcell Madrid grid
Excess pathloss vs. reference case [dB] Frequency
Distance 80 m 200 m 202 m 224 m [GHz]
1 4 1 1 60
Pedestrian 2 4 1 3 30
AP height 10 m 1 3 2 4 6
1 4 3 2 2.6
1 6 3 3 60
Pedestrian 1 5 4 4 30
AP height 5 m 4 7 6 7 6
3 5 6 6 2.6
8 17 8 6 60
Pedestrian 11 18 12 7 30
D2D 11 21 23 17 6
10 19 19 18 2.6
Table 9-5: Pathloss comparison between simulations and measurements in urban microcell.
Microcell Madrid grid
Pathloss in LOS and non-LOS cases [dB]
LOS NLOS Frequency
Distance 80 m 150 m 200 m 202 m 224 m [GHz]
Simulation 76 84 86 105 120 2.6
Measurement 75 83 92 100 110 2.2
Simulation 83 93 93 108 130 6
Measurement 83 90 95 100 110 4.7
Simulation 97 106 108 125 140 30
Measurement 98 104 107 120 125 26
9.3 Summary
The parameter values used for the simulations are marked with * in Table 9-6. The results of
simulated scenarios are listed in Table 9-7. All possible combinations of environments (maps),
frequencies and antenna heights would result in a huge number of result tables. Thus only a
few selected example cases are reported here.
Comparing the statistics of Madrid micro to median values of WINNER Urban micro large-
scale parameters in Table 4-8 of [WIN208-D112] we observe, e.g., the following. The delay
spread is higher in WINNER, in LOS 36 ns vs. 25 ns, and in NLOS 76 ns vs. 44 ns. The
angular spread at the base station side (ASD) is clearly lower in WINNER LOS, 3 vs. 14,
while it is at the same level in NLOS, 15 vs. 15. When comparing the pathloss exponent of
WINNER urban micro LOS condition in Table 4-4 of [WIN208-D112] to Madrid micro LOS we
find a good match, 2.27 in WINNER vs. 2.3 in the table below.
Measured and simulated pathloss exponents in urban environment can be compared based
on figures in sub-sections of 9.1. In urban micro cell type of deployment at 26.4 GHz
METIS Public 85
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
frequency the measured / simulated PL exponents were 3.2 / 4.2 in NLOS and 1.6 / 2.0 in
LOS condition. In vehicle to vehicle LOS case the corresponding values are 2.4 /. 3.0 at
5.25 GHz and 2.1 / 2.9 at 2.3 GHz frequency. Further on, simulated pathloss were compared
to the microcell measurement results conducted in Hatchobori. Simulation and measurement
results for the LOS case match very well throughout the frequency range. The NLOS results
near the street corner match also quite well. However, further from the street corner the
simulation gives significantly larger pathloss values. In reality, there are probably more
reflectors/ scatterers, like traffic lights and street signs, in the street crossings which carry
signal power across the corner.
A very limited set of simulations are shown here and compared with a few measurement
results and to parameters of an existing channel model. Reasonable behaviour of simulated
output parameters can be observed. In principle the validation against measurements is
difficult, because identical environments were not modelled. The same applies to model
parameters of the existing channel model. Anyhow similarities and matching trends can be
found. As a future work outputs of the METIS map-based model must be generated
systematically in different environments, frequency bands, and deployment scenarios. A more
comprehensive validation against measured channel characteristics, both from METIS
measurement campaigns and the literature, must be performed. It should be pointed out that
many of the METIS map-based model parameters are not well calibrated, particularly at
millimetre wave frequencies. Therefore available measurement data should also be used for
improvements by calibration of the model parameters.
Table 9-6: Model parameter used in simulations.
Parameter Madrid Madrid V2V Madrid Madrid micro-cell Madrid micro-cell
D2D* (UOULU) micro-cell* (Docomo) (transition)
Frequency [GHz] 5.25 5.25 5.0 26.45 2.0
Map layout Madrid Madrid Madrid Madrid Madrid
Diffraction model Berg Berg Berg Berg UTD
# TX locations 60 32 32 10 1
# RX locations 501 501 501 501 501
# radio links 30060 16032 16032 5010 2029
Object density [1/m2] 0.01 0.01 0.05 0.05 0.001
Object width [m] 3 3 0.5 0.5 0.5
Object height [m] 1.5 1.5 1.5 1.5 1.5
TX height [m] 1 1.6 10 10 15
RX height [m] 1.5 2.5 1.5 1.5 1.6
Scatterer absorption 0 0 0 0 0
coefficient 𝛼
Specular/ diffuse 0.5 0.5 0.5 0.5 0.5
power ratio
Angle dependency 0.5 0.5 0.5 0.5 N/A
factor 𝑞90
Angle dependency 3.5 3.5 3.5 3.5 N/A
exponent v
Angle dependency 0.95 0.95 0.95 0.95 N/A
factor (HH) 𝛾
METIS Public 86
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Table 9-7: Simulation results for Madrid D2D and Madrid Micro scenario.
Madrid D2D Madrid Micro
Parameter
LOS NLOS LOS NLOS
μlgDS -7.56 -7.39 -7.71 -7.59
lgDS = log10(DS/1 ns)
𝜎lgDS 0.32 0.43 0.31 0.45
μlgASD 1.29 1.27 1.08 1.00
lgASD = log10(ASD/1 deg)
𝜎lgASD 0.27 0.37 0.25 0.39
𝜇lgESD 0.32 0.13 0.27 -0.69
lgESD = log10(ESD/1 deg)
𝜎lgESD 0.33 0.73 0.37 0.83
𝜇lgASA 1.33 0.97 1.37 0.93
lgASA = log10(ASA/1 deg)
𝜎lgASA 0.27 0.48 0.24 0.46
𝜇lgESA 0.32 0.13 0.69 -0.75
lgESA = log10(ESA/1 deg)
𝜎lgESA 0.33 0.73 0.34 0.85
SF in dB 𝜎SF 1.7 17.1 2.1 20.6
𝜇KF -1.8 - 0 -
KF in dB
𝜎KF 2.4 - 3.1 -
ASD vs DS - - 0.5351 0.2723
ASA vs DS - - 0.735 0.278
ASA vs SF - - 0.231 0.3673
ASD vs SF - - 0.3264 0.4247
DS vs SF - - 0.1809 0.1731
Cross-Correlations
ASD vs ASA - - 0.5603 0.5932
ASD vs KF - - -0.002 -
ASA vs KF - - 0.0127 -
DS vs KF - - -0.1017 -
SF vs KF - - 0.0057 -
Pathloss exponent n 2.2 5.6 2.3 5.4
𝜇N 56.9 10.2 47.1 9.0
Number of paths (>-30 dB)
𝜎N 49.2 15.3 28.6 12.4
+
Note! LOS path may be obstructed. Power ratio of direct path vs. all other paths.
METIS Public 87
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
10 Conclusion
This deliverable introduced 5G channel model requirements, propagation scenarios,
measurement results, and the METIS channel models. More specifically, the scenarios and
test cases that have been identified from an end user perspective in an early stage of the
METIS project were analysed and mapped to the METIS propagation scenarios. A new set of
requirements relevant for the radio channel and propagation modelling was derived. Those
propagation scenarios and requirements have been investigated by studying the literature,
conducting extensive measurement campaigns, and exploring several new modelling
approaches. The final METIS channel models consist of a map-based model, a stochastic
model, and their hybrid model to provide a flexible and scalable channel modelling framework.
A comprehensive list of channel model parameters has been derived for diverse propagation
scenarios such as dense-urban macro-cell, micro-cell, indoor, shopping mall, D2D, and V2V
links, with a wide range of frequencies from 2 to 60 GHz. User guidelines, including some
practical examples, are provided for utilizing the models in simulations.
The developed channel models are intended to account for all the requirements and radio
channel characteristics that were identified as important for 5G mobile communications.
However, due to the substantially wider scope of requirements, compared to previous models,
it is very challenging and time consuming to fully validate the developed models by
simulations and measurements. Nevertheless, validation efforts have been made in selected
scenarios by generating channel realizations and comparing the results of these simulations
with the conducted measurements. The channel model outputs showed reasonable and
expected behaviour in comparison to reality, in terms of measured channels and model
parameters of conventional channel models.
As future work, the METIS map-based and stochastic model should be tested systematically
on different environments, frequencies, and deployment scenarios. Based on the literature
survey done during the METIS project, it was recognized that propagation measurement
results between 6 and 60, and above 70 GHz are very limited. Additionally, most of the
measurements have been done at a single frequency band, which means that an
understanding of the frequency dependency of certain propagation parameters is still
emerging. More comprehensive parameterization and validation against measured
propagation characteristics, both from METIS and the literature, is therefore still crucial, in
order to improve the applicable range and accuracy of the models further. Another important
issue is the trade-off between model accuracy and complexity that significantly influences the
scalability; all the METIS models are designed with utmost conscious on this trade-off so that
they are usable in as many concrete test cases as possible with best reliability. Being able to
tackle this trade-off in a very general and holistic test case still remains a challenge.
METIS Public 88
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Appendices
The following content is intended to be used as supplementary reference material. The
content is largely repeated, but substantially more detailed, information from the main
document. The material is intentionally repeated to provide a self-contained material. The
reader who is not interested in supplementary details is referred to the main part of this
document.
METIS Public 89
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 90
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
A.1 Channel measurements at 2.3 GHz and 5.25 GHz in Oulu downtown
The measurements have been executed in the Oulu downtown area. Five measurement
environments (MEs) were covered for METIS Dense Urban Test Case 2 (TC 2):
ME1: Urban Vehicle to Vehicle (V2V), SIMO
ME2: Urban Macrocell Outdoor (O2O), SIMO and MIMO
ME3: Urban Macrocell Outdoor to Indoor (O2I), MIMO
ME4: Urban Microcell Outdoor to Indoor (O2I), MIMO
ME5: Urban Microcell Outdoor (O2O), MIMO
The measurement environments are summarised in Table A-2.
Table A-2: The summary of the measurement environments at 2.3 and 5.25 GHz in Oulu
downtown.
Measurement
ME1 ME2 ME3 ME4 ME5
environment
Link topology UE-UE BS-UE BS-UE BS-UE BS-UE
Propagation
Urban V2V UMa O2O UMa O2I UMi O2I UMi O2O
Scenario
Centre 2.3 GHz and 2.3 GHz and 2.3 GHz and 2.3 GHz and
2.3 GHz
Frequencies 5.25 GHz 5.25 GHz 5.25 GHz 5.25 GHz
TX: Vertical
(dipole) /dual
TX Vertical, RX Dual pol. Dual pol.
Polarisations pol. ±45° (array Dual pol. (±45°)
dual pol. (45°) (±45°) (±45°)
antenna), RX
dual pol. (±45°)
2 Stationary
Oulu One
TX location / One stationary locations, on 2 Stationary
downtown, 0- Stationary
velocity location the roof of locations
30 km/h location
building
TX height 5 m, 10 m, 15 m
above ground 1.6 m 18 m 17.45 m (only for corridor 5 m, 10 m
level measurements)
Stationary, a) Oulu
Hotel room :5 downtown,
Hotel corridor: 7
Oulu Oulu spots/TX Stationary
RX location / spots/floor, Hotel
downtown, 0- downtown, 0- location //floor (vehicle stops
velocity room: 5
30 km/h 30 km/h b) Hotel for
spots/floor
corridor: 17 measurement
spots/floor )
RX height 2.5 m, RX on
1.6 m+ floor 1.6 m + floor
above ground 2.5 m 2.5 m the rooftop of
height height
level a car
TX-RX
20–250 m 20–150 m 40–100 m 15–50 m 30–200 m
distance
RX az. ±180°, RX az.
RX az. ±180°,
TX az. -70 – ±180°, TX az.
Antenna RX az. ±180°, RX az. ±180°, TX az. -70 –
+70°, -70 – +70°,
beamwidth, TX Omnidir., TX Omnidir., +70°,
RX vertical - RX vertical -
Azimuth / RX el. -70° – RX el. -70° – RX vertical -70 –
70 – +90°, TX 70 – +90°,
Elevation +90°, TX ±80° +90°, TX ±80° +90°, TX vertical
vertical -70° – TX vertical -
-70° – +90°
+90° 70° – +90°
METIS Public 91
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Several long
measurements
Several long covering
measurements streets and
covering street corners
Number of
streets and 2.3 GHz:
measurement ~ 30 / floor / > 50 /
street corners 583.567 Hz 27 / floor
s, Channel frequency frequency
2.3 GHz: (SIMO),
sampling rate
583.567 Hz 192.99 Hz
5.25 GHz: (MIMO)
978.847 Hz 5.25 GHz:
978.847 Hz
(SIMO),
2.3 GHz
30 x 16
measurements
(2.3 GHz),
for SIMO (1 x
30 x 18
56) (antennas UE in several
UE in several (5.25 GHz),
TX and RX (TX/RX)) and spots in
spots in different Limited
Remarks move MIMO (30 x different
floors, varying MIMO
simultaneously 16) setup, floors, varying
link length configuration
5.25 GHz link length
due to
measurements
moving
only for SIMO
scatterers
(1 x 50) setup
Campaign
2 days 2 days 2 days 2 days 2 days
duration
Amount of
89.9 GB 43.78 GB 43.67 GB 94.4 GB 4.41 GB
data stored
A.1.1 Measurement equipment and antennas
The measurements were conducted by EB PropSound channel sounder [Ele04]. The
measurement device consists of two separate units: transmitter (TX) and receiver (RX). Both
the units use the same intermediate frequency (IF) of 1.45 GHz and have replaceable RF
units for different RF frequency bands. The architecture of PropSound system is presented in
Figure A-1.
TRANSMITTER RECEIVER
Display & Control
Control Notebook Notebook
Control Control Slave Control Slave Control Slave Control Slave Control
Master Master
AGC
RTCU GPS ASU ASU GPS MMI RTCU
Demodulation
Timing Timing
Code Generator Up Down I/Q
TSIU Modulation TIU Conversion TRU Conversion RRU RIU RSIU
I/Q
TRIG
Digital
Data
Signal
Acquisition
Processing
Storag
e
RBPU
PPS
PropSound uses direct sequence spread spectrum technique and BPSK modulation for
channel sounding. Impulse responses of channel samples are obtained by correlating the
received signal with the same spreading code as was used in the transmission. Sounding in
the spatial domain is employed by switching through multiple antennas in the time domain.
METIS Public 92
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The antenna elements are switched almost instantaneously so that the channel response
remains practically constant within antenna switching period.
A.1.1.1 Antennas
Dipole antennas were used as the TX antenna in ME1 and ME2 scenarios. Their properties
are presented in Table A-3.
Table A-3: TX antenna properties.
Antenna designation Dipole_2G45 Dipole_5G25 (Modified 5G version of
AV1433 WLAN antenna)
Frequency / Bandwidth 2.25 -2.55 GHz 5.2—5.3 GHz
Uniformly spaced linear array antennas (ULAs) illustrated in Figure A-2 were used as TX
antennas in ME3, ME4 and ME5 scenarios.
Figure A-2: TX ULA antenna for ME3, ME4 and ME5 scenarios at 2.3 GHz (left) and 5.25 GHz
(right).
METIS Public 93
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
The RX antennas for 2.3 GHz and 5.25 GHz are presented in Figure A-3 and Figure A-4,
respectively.
METIS Public 94
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 95
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 96
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
METIS Public 97
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
RX antenna
ARA 118/A for
frequency spectrum
scanning
TX antenna
TX RX
Figure A-6: TX leading and RX follows, vehicles move in the same directions.
Figure A-7: RX leading and TX follows, vehicles move to the same directions.
METIS Public 98
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure A-8: TX (blue line) and RX (red line) routes moves to the opposite directions.
Figure A-9: TX/RX (blue/red dot) stationary position and RX/TX (red/blue line) moves on the
cross street.
TX antenna
Figure A-10: TX antennas: 5.25 GHz dipole (SIMO) and 2.3 GHz ULA (MIMO).
A.1.2.2.1 Measurement routes
TX antenna location (65° 0.647' LAT, 25° 28.987' LON, blue circle) and RX routes (red line)
are presented in Figure A-11 to Figure A-14.
METIS Public 99
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure A-16: Measurement spots in the hotel room for the ME3 scenario.
17 spots were measured on the hotel corridor on floors 2 to 5. In addition, 9 spots were
measured on the 6th floor. The layout of the hotel corridor is similar for the 3rd, 4th and 5th
floor. Figure A-17 presents the measured spots in the corridor of these floors.
Figure A-17: The floor plan of the hotel with the measurement spots in the corridors on floors 3
to 5.
Figure A-18 shows the detailed distance information of the corridor and Figure A-19 shows the
example of the measurements on the 8th measurement spot (cf. Figure A-17) on the 5th floor.
Figure A-20 presents the view from the TX location. Figure A-21 shows the layout and
collected distance information on the 2nd floor and 6th floor.
Figure A-19: Measurement example, TX was located on the roof of neighbouring building.
Figure A-20: Views from the TX site towards the target building: from position 2 (left) and from
position 1 (right).
Figure A-21: The hotel layout and measurement spots on the corridor of 2nd floor (left) and 6th
floor (right).
A.1.2.4 ME4: Urban microcell outdoor-to-indoor (MIMO)
The measurements were performed within a multi-storey building in the same manner as in
ME3. The RX location was either the hotel room or the corridor of hotel. The TX unit and
antenna were located on an articulated crane. The TX and the RX were stationary during the
recordings. Emphasis was on collecting data for the TX elevation angle spread analysis.
A.1.2.4.1 Measurement spots
Five measurement spots were selected in the hotel room and one measurement spot in the
end of the corridor facing to Saaristonkatu. The measurements were repeated in the same
spots on floors 2 to 5 in the corresponding rooms (room numbers 258, 358, 458 and 558,
Figure A-22). The height of the TX antenna was 5 m and 10 m and the distance in relation to
the reference corners of the target building is depicted in Figure A-23.
Figure A-22: Measurement spots in the hotel room and in the end of corridor at 2.3 GHz and
5.25 GHz.
In the second part of the ME4 measurements, 6 spots were measured in the hotel corridor on
floors 2 to 5. Furthermore, one measurement was executed in the elevator, having the TX
directed towards the elevator. The measurements were performed with TX heights of 5 m,
10 m and 15 m (Figure A-24). It should be noted that the 6th measurement spot was skipped
in 5.25 GHz measurements due to lack of signal. Figure A-26 presents the distance
information collected from these measurements. Notice that the image orientation of Figure
A-26 is rotated 90 degrees CCW with respect to the previous indoor sketches.
METIS Public 104
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
Figure A-24: The TX mounted on an articulated crane at a height of 5 m (left) and 15 m (right).
Figure A-26: Measurement spots on the hotel corridor for ME4 at 2.3 GHz and 5.25 GHz.
A.1.2.5 ME5: Urban microcell outdoor (MIMO)
In ME5, the measurements were performed using a single RX vehicle in the same manner as
in ME2. The TX unit and antenna were located on an articulated crane at 5 and 10 m above
the ground. The RX setup is similar as in ME1. The RX is moving close to the TX location
gathering pathloss information along the street and around the street corners. The
measurement parameters are selected balancing between moderate scatter velocities and
maximum MIMO array size. Emphasis was on collecting data for TX elevation angle spread
analysis. The measurement routes are presented in Figure A-27. The TX position is marked
as a blue dot. Eight to twenty measurement spots, i.e., vehicle stops for data recording, were
recorded on the three different measurement routes.
Figure A-27: Measurement routes for ME5 at 2.3 GHz and 5.25 GHz.
A.1.3 Measurement results
A.1.3.1 Urban V2V (ME1)
Figure A-28 and Figure A-29 present the pathloss for the case measurements where the
vehicles are moving to the same directions and the opposite directions.
Similar statistic for the decorrelation distance of large-scale parameters has not been found for
V2V channel in the existing literature. Therefore, we compare the obtained results of
decorrelation distances to the UMi BS-MS LOS scenario in [ITUR09-2135], which is the
closest to V2V scenario by means of antenna heights. There are no significant differences in
the correlation distance for delay spread and shadow fading in comparison with parameters in
[ITUR09-2135]. However, the obtained decorrelation distance for K-factor is much smaller in
comparison to the results of [ITUR09-2135]. Since TX and RX are moving simultaneously, the
propagation channel is more dynamic and the propagation characteristics in terms of small-
scale fading are different in comparison with a traditional BS-MS scenario. Moving scatters,
i.e., other cars, cause strong reflections to the measured V2V channel, which has likely
shortened the decorrelation distance for the K-factor [RJM+14].
A.1.3.2 UMa O2O (ME2)
The results for UMa O2O MIMO measurements are presented in Section 5.5.
A.1.3.3 UMa O2I (ME3)
RX in the hotel room
The signal-to-noise ratio was around 15 dB in the measurements where the RX was located
on the 2nd floor and 3rd floor. Therefore, the spatial analysis cannot be carried out for the
recorded data of lowest floors.
Figure A-30 presents the received power for the measurement over the building rooftop in O2I
scenario. The TX antenna had two different positions on the roof (cf. Figure A-16). The knife-
edge diffraction and transmission loss caused by a window (cf. Section 3.3.2.1) was taken into
account in the model.
Table A-8: Parameter for UMa O2I corridor scenario at 2.3 GHz (ME3).
All Floors [3GPP14-
Parameter Symbol
36843]
OLOS NLOS All
𝜇lgDS -7.76 -7.52 -7.65 -6.62
lgDS = log10(DS/1 ns)
𝜎lgDS 0.2 0.19 0.22 0.32
𝜇lgASD 0.28 0.4 0.34 1.25
lgASD = log10(ASD/1 deg)
𝜎lgASD 0.23 0.21 0.23 0.42
METIS Public 108
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
ASA vs K
DS vs K
SF vs K
ESD vs SF 0 -0.7 -0.7 0
ESA vs SF -0.1 -0.3 -0.3 0
ESD vs K
ESA vs K
ESD vs DS -0.1 -0.3 -0.5 -0.6
ESA vs DS 0.5 0 0 -0.2
ESD vs ASD 0.3 0.3 0.2 -0.2
ESA vs ASD 0.1 0 0 0
ESD vs ASA -0.2 -0.1 -0.2 0
ESA vs ASA 0.8 0.5 0.6 0.5
ESD vs ESA -0.3 -0.1 0 0.5
Delay Distribution Exp.
AoD and AoA Distribution Wrapped Gaussian
EoD and EoA Distribution Laplacian
𝑟𝜏
Delay Scaling parameter 2.2 1.7 1.9 2.2
𝜇XPR 12 10 11 9
XPR in dB
𝜎XPR 11 10 11 11
Per cluster shadowing in dB 𝜁 5 5 5 4
EoD in deg 𝜇EOD -11.5 -10.8 -11.2
** Parameters taken from [WIN+10-D53]
𝜎XPR 10 10 11 10 10 11 10 10
Per cluster
𝜁 6 3 6 3 6 3 6 3
shadowing in dB
EoD in deg 𝜇EOD 4.9 7.7 -3.3 -1.4 -16 -12 -4.5 -1.3
For the TX antenna height of 5 m, the SNR was insufficient for receiving signal beyond the
mound due to diffraction loss.
Table A-11 summarizes the results obtained for UMi O2O measurements at 2.3 GHz.
The ESD with respect to the distance between the TX and the RX is tested against linear
model and negative exponential model. Figure A-33 shows ESD with respect to the distance
between the TX and the RX for the LOS scenario.
Figure A-33: ESD vs. link distance in LOS scenario, TX height 5 m (left), TX height 10 m (right).
Figure A-34 shows ESD with respect to the distance between the TX and the RX for the NLOS
scenario.
Figure A-34: ESD vs link distance in NLOS scenario, TX height 5 m (left), TX height 10 m (right).
The distance dependency of ESD was observed to follow negative exponential model and
linear model in the LOS and the NLOS case [RHM+14]. The parameters for ESD distance
dependency model are given in Table A-12.
𝛾 is the standard deviation between the model and the measured ESDs
parameter spreads, we classified the paths into 10 clusters by the K-Power-Means based
clustering algorithm [SKI+12]. In the clustering algorithm, the normalization parameters of the
Multipath-Component-Distance calculation were ToA 1.0s and AoA ( AoD ) 360 deg .
Sleeve antenna
2.9m
1.45m
UE was mounted on a carriage and the channel was measured continuously every 1 s while
UE moved along the courses at about 1 m/s that is the typical walking speed.
P1
N
pedestrian
crossing BS(UMi-Scenario)
BS(D2D-Scenario)
A5
P2
(B)The view of station plaza
ticket
A4 A3 in the daytime (10:00-23:00)
gate
A2
ticket
A1 gate
measurement course
UMi scenario
D2D scenario
D2D scenario
UMi scenario
D2D scenario
UMi scenario
UMi scenario
UMi scenario
The spreads of the propagation parameters for all measurement conditions are summarized in
Table A-14. For reference, the large-scale parameters of urban micro cell scenario of ITU-R
M.2135 channel model are also shown. Although the spreads of propagation parameters
almost corresponded to the parameters of ITU-R M.2135 channel model, the delay spreads
are smaller and the ASDs are larger than the spreads of the channel model. In UMi scenario
in the daytime, there was no significant difference between VV-polarized channel and HH-
polarized channel in regard to the received power, ASA, ESA, ASD, and ESD. The median
XPR from V-pol. to H-pol. was 8 dB, and the median XPR from H-pol. to V-pol. was 5 dB.
In order to investigate the frequency dependence of the pathloss, the measurement system at
4.7 GHz was also utilized for a measurement at 2.2 GHz.
Figure A-39: Map of Nihonbashi in Tokyo (typical Manhattan grid layout environment).
Some photographs of the measurement environment are shown in Figure A-40. The
measurement sites are surrounded by tall buildings with heights of approximately 40 m, which
is considered to be a Manhattan-like grid layout.
60 60 60
80 80 80
Figure A-41: Comparison of measurement and M.2135 results for multiple frequencies (same as
Figure 5-2).
From Figure A-41, it is found that, the M.2135 model shows a good match with the
measurement results for all frequencies, while they do not match well for the NLOS routes.
Especially, when the distance between TX and the intersection of the route is large, the
difference between M.2135 and the measurement results becomes relatively large (e.g.,
NLOS routes 1 and 3). Figure A-42 and Figure A-43 show the TX height influence on the
pathloss at 37 and 4.7 GHz, respectively.
60 60 60
80 80 80
Figure A-42: TX antenna height influence on the pathloss; M.2135 and measurement results are
compared at 37 GHz for different antenna heights.
60 60 60
80 80 Path loss [dB] 80
Path loss [dB]
Path loss [dB]
Figure A-43: TX antenna height influence on the pathloss; M.2135 and measurement results are
compared at 4.7 GHz for different antenna heights.
Figure A-42 illustrates the TX antenna height influence on the pathloss at 37 GHz. It is verified
that the differences between M.2135 calculation and the measured pathloss values are
reasonably matched for NLOS route 1. On the other hand, they do not match well for NLOS
routes (especially, NLOS routes 1 and 3).
Figure A-43 illustrates the TX height influence on the pathloss at 4.7 GHz. The difference
between M.2135 calculation and measurement results for a TX height of 1.5 m is large
compared to the other two different TX heights of 6 and 10 m. However, the trend of the
difference between M.2135 calculation and measurement results may be similar to the
observation for 37 GHz of Figure A-42.
All in all, M.2135 shows good estimation for the LOS while there are relatively large
differences in case of NLOS. No clear frequency and TX height dependencies are observed.
In order to evaluate the pathloss in high SHF/EHF bands, a different calculation method such
as ray-tracing similar to the map-based model approach may provide a better estimation.
-5 dB -5 dB
-10 dB -10 dB
-15 dB -15 dB
The measured wall, window and door losses were 7.5, 1.0 and 11.5 dB respectively as
indicated by Figure A-49. The loss in excess of free space loss for the case corridor-to-room
(TX location 1) and room-to-room (TX location 2) was 29 dB and 48 dB which is double the
loss at 2.4 GHz in dB units. For TX location 1 the antenna was directed along the corridor and
for TX location 2 the antenna was directed towards the corridor.
60 GHz
2.4 GHz
Figure A-51: Relative RX power measured and modelled at 2.4 GHz and 60 GHz for the long
range corridor measurement.
̂ )𝑑𝛺 = 1.
∫ P(𝒖 (A-1)
𝝁𝒖̂ = ∫ 𝒖 ̂ )𝑑Ω ,
̂ P(𝒖 (A-2)
It is also possible to calculate the angle spread for a certain direction, e.g. parallel to the x, y,
or z-axis as:
180
[𝜎𝑑𝑖𝑟 ]𝑛 =
𝜋
√∫[𝒖
̂ − 𝝁𝒖̂ ]2𝑛 P(𝒖
̂ )𝑑Ω , 𝑛 = 𝑥, 𝑦, 𝑧 . (A-4)
Interesting findings are that there are some distinct spikes in the power delay profile (PDP) for
which the directional spectrum also consists of single spikes. This is interpreted as specular
reflections. Between the spikes of the PDP the directional spectrum is rich and substantially
spread out suggesting that the scattering is caused by a multitude of small objects or rough
surfaces. Moreover, the directional spread with respect to the vertical z-axis and one
horizontal y-axis decay faster than with respect to the horizontal x-axis. The x-dimension is the
largest and the z-dimension the smallest of the room which may explain the different decay
times.
Figure A-52: Directional spread and angular spectra shown for delays indicated with numbers
1-5.
The RF part of the used measurement system is illustrated in Figure A-53 (left). The system is
based on using a vector network analyser (VNA) with an intermediate frequency (IF) of 5 –
9 GHz. The RF frequencies 61 – 65 GHz are generated with up and down converters and a
LO operating at 14 GHz. Both TX and RX sides are connected by cables to the VNA.
The up converter and the TX antenna are on a rotator. A 20 dBi standard gain horn is used as
the TX antenna. An omnidirectional biconical horn antenna is used as the RX antenna. A
photograph of the antennas is presented in Figure A-53 (right). The TX antenna is rotated in
the azimuth direction from 0° to 360° with 3° steps. Three samples of the amplitude and phase
are measured at each direction with 2001 frequency points. The 4 GHz IF band width leads to
a 0.25 ns delay resolution and the maximum delay is 500 ns. A direct connection back-to-back
calibration is used to compensate the transfer function of the measurement system. Both, TX
and RX antennas have relatively narrow elevation plane radiation patterns which limits the
measurements to the azimuth plane. Both antennas are vertically polarised. Cross-polarisation
is measured by rotating the horn antenna by 90°.
Figure A-53: Measurement system and sounder configuration (left) and photograph of TX and
RX antennas (right).
Figure A-54: Photographs of the measurement sites of the 60 GHz channel measurements at the
first and third floors of the Sello shopping mall.
Escalator
Stage
Figure A-55: Floor plan of the 1st floor of the Sello shopping mall with TX and RX locations.
Figure A-56: Floor plan of the 3rd floor of the Sello shopping mall with TX and RX locations.
Tx11
18
Tx7
13.5
Rx1 Tx1
-4.5
The scattering model contains two parameters, a scattering coefficient S and a scattering lobe
width αR, which relate to the material properties of the local surface. Since it is not feasible to
obtain material parameter values for all local surfaces, we calibrate S and αR such that the
predicted and measured rms delay spread agree as well as possible using the same
parameter values for all points in the prediction of a single channel.
In the present work, point cloud-based field prediction has been used in two scenarios, an
indoor cafeteria and an outdoor square, in which measurement campaigns were conducted in
order to calibrate the scattering model parameters for the field prediction tool. For the
cafeteria, the LOS scenario was calibrated with 3 measured links resulting in the parameters
𝑆 = 0.9, 𝛼𝑅 = 20. A comparison of the measured and predicted PDPs is shown in Figure
A-63. In the open square, the scattering model parameters were found to be 𝑆 = 1.0 ,
𝛼𝑅 = 50.
-70 Measured
Predicted
-80
Amplitude [dB]
-90
-100
-110
-120
-130
0 20 40 60 80 100 120
Delay [ns]
Figure A-63: Comparison of measured and predicted PDPs.
A.5.4 Post-processing and detection of propagation paths
Post-processing is same for both the measurement results from the shopping mall and for the
point cloud results from the cafeteria and for the square. The radio propagation channel is
characterised as a sum of multipaths. In large spaces specular mechanisms represented by
propagation paths dominate over diffuse scattering [HJK+14], and therefore it is sufficient to
characterise only the resolvable multipaths from results with large bandwidth.
Multipath components are identified from PDP. Multipath components are local maxima peaks
in the PDP. Point on the PDP curve 𝑃(𝜏) is considered as a peak if
1 𝜏+𝛥/2
𝑃(𝜏) > 𝛥 ∫𝜏−𝛥/2 𝑃(𝑡) 𝑑𝑡, (A-5)
and
τ ≥ 𝑑/𝑐0 , (A-9)
where 𝛥 = 2.0 ns is the length of sliding window over delays, 𝛥𝜏 = 0.25 ns is the delay
resolution, 𝑑 is the distance from TX to RX, and 𝑐0 is the speed of light. The length of the
sliding window, compared to the delay resolution, needs to be selected carefully in order to
detect all the peaks. An example of PDPs with the detected peaks is presented in Figure
A-64.
The parameterisation is done based on 20 dB dynamic range, i.e., only the propagation paths
within 20 dB of the strongest path are used.
Figure A-64: Examples of measured LOS and OLOS PDP’s with the detected peaks in the
shopping mall.
Angles for the detected propagation paths are found from the PADP after the peak detection
from the PDP. The path amplitudes are calculated from the path PADP amplitudes taking in to
account the effects of antenna gains and the windowing loss. The detected path amplitudes,
delays, and angles can be illustrated in a discrete form of the PADP, as in Figure A-65.
Only the arrival angles are available from the channel measurements in the shopping mall.
The departure angles are can be calculated directly from the angels and amplitudes of the
detected paths. Knowing the arrival angles, excess delays, and BS and MS locations, the
scattering points can be identified and assuming single-bounce condition the angles can be
calculated from the BS point of view.
N
Ptot pi . (A-10)
i 1
The sign of the shadow fading term is defined so that increasing values of SF correspond to
increasing received power [WIN208-D112]. The measured and simulated pathloss and the
fitted pathloss models are illustrated in Figure A-66 and Figure A-67.
Figure A-66: Measured pathloss and pathloss model in shopping mall as a function of the link
distance.
Figure A-67: Simulated pathloss and pathloss model in cafeteria (left) and square (right) as a
function of the link distance.
A.5.5.2 Delay spread
Cumulative distribution functions of measured and simulated delay spreads are compared to
delay spreads generated by WINNER implementation in Figure A-68, Figure A-69 and Figure
A-70.
Figure A-68: A CDF comparison to WINNER; measured delay spread in shopping mall LOS (left)
and OLOS (right).
Figure A-69: A CDF comparison to WINNER; simulated delay spread in cafeteria LOS.
Figure A-70: A CDF comparison to WINNER; simulated delay spread in square LOS (left) and
OLOS (right).
A.5.5.3 K-factor
For the LOS-cases K-factor is calculated as the power ratio between the LOS component and
the total power is all the other detected propagation paths within the 20 dB dynamic range.
Cumulative distribution functions of measured and simulated K-factors are compared to the
normal distributions generated by WINNER implementation in Figure A-71, and Figure A-72.
Figure A-71: A CDF comparison to WINNER; measured K-factor in shopping mall (left) and
cafeteria (right).
The constant CAS that is a scaling factor related to the total number of clusters N and is given
in [WIN208-D112] is given only up to 𝑁 = 20. For 𝑁 < 20 simple linear interpolation is used.
For 𝑁 > 20 extrapolation is used
Based on the CDF comparisons in Figure A-73 - Figure A-78 the azimuth angular spreads are
quite similar in the measurements and simulation results as the ones calculated from the
clusters generated by WINNER II.
Figure A-73: A CDF comparison to WINNER; measured ASD in shopping mall LOS (left) and
OLOS (right).
Figure A-75: A CDF comparison to WINNER; simulated ASD in square LOS (left) and OLOS
(right).
Figure A-76: A CDF comparison to WINNER; measured ASA in shopping mall LOS (left) and
OLOS (right).
Figure A-78: A CDF comparison to WINNER; simulated ASA in square LOS (left) and OLOS
(right).
Figure A-79: Simulated ESD (left) and ASA (right) in cafeteria LOS and square LOS and OLOS.
Figure A-80: Measured XPR, cross-polarisation measurement noise limit, and mean (black solid
line) and the 95 % tolerance interval limits (black dashed lines) fitted normal distribution.
(in dBs) of the cluster power as a function of the WINNER-type cluster AS is presented in
Figure A-82.
The measured average standard deviation of the cluster power is 0.65 dB and 1.1 dB over
four and 12 wavelengths, respectively, corresponding to cluster AS of 0.63° and 0.38°,
respectively. The measured average maximum-to-minimum difference of the cluster power is
2.3 dB and 4.6 dB over four and 12 wavelengths, respectively, corresponding to cluster AS of
0.67° and 0.48°, respectively.
Based on these results we propose to use cluster AS of 0.5° at the 60-GHz range. This cluster
AS generates realistic variation of channel coefficients for each cluster and each receiver and
transmitter element over an array antenna that is approximately four to 12 wavelengths
across.
Figure A-81: a) Example of a measured power and b) cluster power over about 12 wavelength
long measurement path.
Figure A-82: Average standard deviation and maximum-to-minimum difference of the cluster
power within 20 dB dynamic range as a function of cluster AS over a) four wavelengths and b)
12 wavelengths.
where 𝑟𝜏 is the delay scaling parameter, 𝑡𝑚𝑎𝑥 is the maximum excess delay, and 𝜎𝜏 is the
delay spread. In LOS-case, the maximum excess delays 𝑡𝑚𝑎𝑥 in need to be replaced by
𝐿𝑂𝑆 𝐿𝑂𝑆
𝑡𝑚𝑎𝑥 𝐷, where 𝑡𝑚𝑎𝑥 are the measured or simulated maximum excess delay and D is a scaling
factor that depends on the K-factor (see Equation (4.3) in [WIN208-D112]).
The resulting single slope exponential power delay profile is used to determine the per-cluster
fading by optimising the associated distribution with log-likelihood function as in [GBT14]. The
dynamic range limit taken into account in a similar manner as the noise level with
measurement results. The average delay scaling parameter and per-cluster fading is then
used for each scenario. CDFs of delay scaling parameter and per-cluster fading are given in
Figure A-83, Figure A-84 and Figure A-85. CDFs of maximum excess delay are given in
Figure A-86, Figure A-87, and Figure A-88.
Figure A-83: CDF of delay scaling parameter and per-cluster fading in shopping mall.
Figure A-84: CDF of delay scaling parameter and per-cluster fading in cafeteria.
Figure A-85: CDF of delay scaling parameter and per-cluster fading in square.
Figure A-86: A CDF comparison to WINNER; measured maximum excess delay in shopping mall
LOS (left) and OLOS (right).
Figure A-87: A CDF comparison to WINNER; simulated maximum excess delay in cafeteria LOS.
Figure A-88: A CDF comparison to WINNER; simulated maximum excess delay in square LOS
(left) and OLOS (right).
Figure A-90: Simulated correlation distances in square LOS (left) and OLOS (right).
Figure A-92: A CDF comparison to WINNER; simulated number of propagation paths in cafeteria
LOS.
Figure A-93: A CDF comparison to WINNER; simulated number of propagation paths in square
LOS (left) and OLOS (right).
A.5.5.11 Parameters
The WINNER parameters are given in Section 5.12. In the WINNER II implementation those
parameters go to a file called ScenParTables.m. In addition to those parameters it is important
to use 'IntraClusterDsUsed' = 'no' and 'DelaySamplingInterval' = 0.125 ns (= 1/(2·band width))
in a file called wimparset.m. Naturally, also the centre frequency needs to be changed to
63 GHz.
The presented WINNER II channel model parameters have been validated by comparison of
the cumulative distributions of delay spreads and azimuth angular spreads.
the two transmit antennas. At the receiver side, two fully parallel RF and baseband chains are
used.
Figure A-95 shows the antenna power patterns of the utilized antennas at 10 and 60.4 GHz,
respectively. Both antenna types are omni-directional, vertically polarized and show a similar
characteristic. The same type of antenna is used at the TX and RX ends. The commercially
available Flann Model MD249-AA antenna [Fla15] was used for the 60 GHz band, while a
proprietary quarter-wavelength dipole antenna was used for the 10 GHz band.
Figure A-95: Horizontal (top) and vertical cuts (bottom) of the normalized antenna power
patterns in dB at 10 (left) and 60.4 GHz (right).
Other important parameters for the channel sounder are listed in Table A-18.
Figure A-98: Pathloss results at 10 and 60.4 GHz for LOS (left) and NLOS (right) condition.
Figure A-98 shows the scatter plot of all pathloss values based on the APDP results, divided
by frequency and propagation condition. The solid and dashed lines show the linear least
squares fit to the data samples. For the pathloss fitting the following formula was used:
𝑑
𝑃𝐿(𝑑)|dB = 𝑃𝐿(𝑑0 )|dB + 10𝑛 log10 (𝑑 ) , (A-13)
0
where 𝑃𝐿(𝑑0 )|dB denotes the pathloss at reference distance 𝑑0 , and 𝑛 denotes the pathloss
exponent. In the LOS case the intercept was fixed to be equal to the free space pathloss,
while in the NLOS case both, i.e. intercept and exponent value, were obtained using
parameter fitting. As indicated in the right figure, not all data samples were used for the fitting,
since the influence of the noise floor was too big for distances greater than 100 m.
Table A-19 shows the intercept and exponent values for the different frequencies and
propagation conditions. The main difference in the PL intercept for 10 and 60 GHz is due to
the different effective antenna apertures.
The delay spread was calculated as the second moment of the APDPs. Figure A-99 shows
exemplary APDPs for 10 and 60 GHz. As indicated in the figures, a threshold of 25 dB was
used to remove the noise floor.
Figure A-99: Exemplary normalized APDP for 10 (left) and 60 GHz (right).
Figure A-100: Delay spread CCDF for LOS (left) and NLOS (right) condition.
The complementary cumulative density function (CCDF) of the delay spreads for the LOS and
the NLOS cases are shown in Figure A-100. It can be seen, that the delay spread for 60 GHz
is in general smaller compared to 10 GHz. Since both antenna pairs were co-located, the
difference has to be explained by the difference in frequency and therefore relative roughness
of the surrounding object surfaces. Another explanation might be the additional attenuation at
60 GHz due to the Oxygen absorption.
However, because of the relatively narrow bandwidth the two centremost bands (between
18.1 GHz and 19.3 GHz) are not used to decrease the total measurement time. The set of
measurements is limited to the indoor measurements.
A.7.3 Measurement locations
The map below shows the preliminary plan for the measurement locations. There are quite
many measurement cases and probably some of them will be dropped because of time
constraints.
The idea below is to perform three different sets of measurements. The legends with the red
symbols and letters R indicate the receiver at the height of 1.5 m. The blue legends (T) mean
the transmitters at the height of 1.5 m (user is moving). The transmitters are in two different
environments; the open office area T1 to T7 and the open café area T9 to T10. The green
legends (S) indicate the case that the receivers and transmitters are on the table (height 1 m).
T2 R4
T3 = Receiver @ 1.5 m
R3 T1 S4 = Transmitter with 0° dir @ 1.5 m
S2 S3 = Transmitter / receiver @1 m
R2 R5
S1 S5 T5
T6
R1 T4
T7
T9 R7
T8
R6
T10
During 2014, WP6 of METIS project has refined the assumptions for different propagation
environments in order to streamline the system simulations for the twelve TCs [MET15-D65].
Table B-1 summarizes the assumptions. Efficient deployment of the Madrid grid is feasible in
most of the TCs.
Table B-1: Refined test cases.
TC # TC name Carrier freq. Propagation model Notes
[GHz]
TC1 Indoor office 60 & 3.5 METIS indoor office High data rates
grid
TC2 Dense urban 3.5 & 2 Madrid grid High data rates, D2D important
TC3 Shopping mall 60 METIS shopping mall Large array antennas (64x16
grid MIMO), large data rates
TC4 Stadium 60 & 2.6 NA D2D important, high data rates
TC5 Teleprotection 3.5 & 2 Madrid grid
TC6 Traffic jam 2.6 & 0.8 Madrid grid D2D important
TC7 Blind spot 2 Madrid grid Macro cell
TC8 Real-time 3.5 & 2 Madrid grid
computing
TC9 Open air festival 60 & 2-6 METIS model D2D important
TC10 Emergency 3.5 & 2 & 0.8 Madrid grid Like TC2 but ruined environment,
communications I2O important
TC11 Massive use of 2 Madrid grid D2D important, low data rates,
sensors/actuator wide coverage
s
TC12 Traffic efficiency 5.9 Madrid grid D2D important
and safety
Since one TC can cover several PSs, mapping between PSs and TCs is needed as presented
in Table B-2. Lower case ‘x’ means that the PS is for below 6 GHz frequencies and capital ‘X’
means that the PS is for both below and above 6 GHz. For brevity the backhaul link type is not
shown. Drawings of some selected PSs are presented below in Figure B-1, Figure B-2, Figure
B-3, and Figure B-4 [MET13-D61].
Test Case
Propagation Scenario
1 2 3 4 5 6 7 8 9 10 11 12
Urban Micro O2O X X x x x x X X x
Urban Micro O2I x x x
Urban Macro O2O x X x x x x X x x x
Urban Macro O2I x x x x
BS-UE
Rural O2O x x x x X x
Indoor Office X X x
Indoor Shopping mall X X
Open Air Festival x
Stadium x
Urban O2O (also V2V) x x x x x x x x
Urban O2I x x x
D2D
Rural O2O x x x
Indoor Office X X
Indoor Shopping mall X x
Highway V2V x x
Open Air Festival x
Stadium x
Figure B-2: Madrid grid for Urban Micro and Urban Macro. [MET13-D61]
200 m
300 m
Figure B-3: Shopping mall. [MET13-D61]
The simulation region in the stadium is shown above as the brown 30x40 m area. The slope
and dimensions of the spectator stand are shown in the rightmost figure. Other data is given
below:
The user relative height is 1 m above tribune level.
The horizontal distance between rows is 1 m and the seats are situated in 0.5 m
intervals within a row.
For D2D transmission additional propagation loss, e.g. 3 dB/m, is added to account for
human body loss attenuation.
Although no mobility of users is assumed for Stadium, a velocity of 3 km/h may be
used to account for small-scale effects.
A generic layout of the highway scenario is shown in Figure B-5. Parameters for the highway
scenario (two cases) are shown in Table B-3.
C.1 Model
C.1.1 LOS and diffracted pathways
The LOS and diffracted pathways are described by the Berg recursive model [Ber95]. This is a
semi-empirical model designed for signal strength prediction along streets in an urban
environment. It is semi-empirical in the sense that it reflects physical propagation mechanisms
without being strictly based on electromagnetics theory. It is based on the assumption that a
street corner appears like a source of its own when a propagating radio wave turns around it.
The corners of buildings and the antennas represent nodes (See Figure C-1 left).
Any two subsequent nodes must be in LOS with respect to each other. Moreover, for any
three subsequent nodes the middle node is blocking the LOS between the first and the third
node. Along a propagation path each node contributes a loss which depends on the change in
direction 𝜃. The loss at a specific node 𝑗 is a result of all previous nodes and is given by the
well-known expression for free space pathloss between isotropic antennas where a fictitious
distance 𝑑𝑗 is used, i.e.
4𝜋𝑑𝑗
𝐿𝑗|dB = 20 log10 ( 𝜆 ) (C-1)
where 𝜆 is the wave length. It should be noted that the fictitious distance corresponds to the
real distance multiplied by a factor at each diffraction node. The result is that the fictitious
distance 𝑑𝑗 becomes longer than the real distance meaning that it accounts for diffraction loss
when used in the free space loss expression (C-1). An example with four nodes is shown in
Figure C-1 (right). At each node, the fictitious distance is given by the following recursive
expression
where 𝑠𝑗 is the real distance between node 𝑗 and its following node (𝑗 + 1), 𝑞𝑗 is a function of
𝜃𝑗 (See Figure C-1 middle). The initial values are 𝑑0 = 0 and 𝑘0 = 1.
TX
y [m]
x [m]
Figure C-1: Example of a street corner acting as a node (left), Manhattan map (middle),
topological example with four nodes (right).
𝜃𝑗 𝑣
𝑞𝑗 = 𝑞90 (90 deg) , (C-4)
where 𝜃𝑗 ∈ [0, 180] deg , 𝑞90 and 𝜈 are parameters determined by fitting the model to
measurement data. The parameter 𝑞90 accounts for the amount of diffraction loss caused by
each node. A larger value results in larger diffraction loss. The corresponding frequency
dependency is given by
𝑞
𝑞90 = √ 𝜆𝜆 , (C-5)
where 𝑞𝜆 is a model parameter given in Table 6-1. The parameter 𝜈 accounts for how fast the
loss changes in the transition zone between LOS and NLOS. These two parameters are
known for frequencies below 6 GHz but need calibration by measurements for higher
frequencies. Further, there is a corresponding polarimetric matrix for these parameters
expressed in terms of horizontal and vertical polarizations as described in Appendix C.1.8.2.
In order to provide the received power for each path the antenna gains of each end of the link
are applied
TX TX RX RX 𝐽𝑛
𝑃RX TX
𝑛|dB = 𝑃|dB + 𝐺|dB (𝐤𝑛 ) + 𝐺|dB (𝐤𝑛 ) − 𝐿𝐽𝑛 |dB ({𝜃𝑗 } ) (C-6)
𝑗=0
TX RX
where 𝐺TX RX
|dB (𝐤𝑛 ) and 𝐺|dB (𝐤𝑛 ) are the TX and RX antenna gains for the corresponding wave
vectors 𝐤TX RX
𝑛 and 𝐤𝑛 , while 𝐽𝑛 denotes the number of interactions of path 𝑛.
circles. The corresponding reflection surfaces are marked with green lines. This procedure
may be repeated to achieve any number of diffraction and specular reflection interactions.
When repeated, the nodes of previous steps act as TX / RX node in the first step.
TX
y [m]
x [m]
TX
y [m]
RX
x [m]
Figure C-3: Determined pathways for the Berg recursive diffraction model. Here up to four
interactions of specular reflection and diffraction have been accounted for.
y [m]
RX
x [m] x [m]
Figure C-4: Diffracted paths between TX and RX (left) and relative power over the RX route
(right).
Power [dBm]
Distance [m]
Figure C-5: Received power over the RX route for isotropic antennas and 0 dBm transmit power.
Propagation distance [m]
Azimuth angle [deg]
Figure C-6: Path angles and propagation distances at RX along the route for direct paths only.
The power scale is relative to the strongest path for each RX location.
perpendicularly oriented with respect to the line connecting the two nodes of the link in the
projection from above. This means that as either node is moving the screen turns around a
vertical line through the centre of the screen such that it is always perpendicular to the line
connecting TX and RX. The corresponding shadowing loss is modelled using a simple knife
edge diffraction model for the four edges of the screen as
where 𝐹ℎ1 ,𝐹ℎ2 and 𝐹𝑤1 ,𝐹𝑤2 account for knife edge diffraction at the four edges corresponding
to the height, ℎ, and width, 𝑤, of the screen (See Figure 6-3). The shadowing for a single edge
is given by
𝜋 𝜋
atan(± √ (𝐷1 +𝐷2 −𝑟))
2 𝜆
𝐹= 𝜋
, (C-8)
where 𝜆 is the wave length, 𝐷1 and 𝐷2 are the projected distances (according to the
projections from side and from above in Figure C-7) between the nodes and the edges of the
screen and 𝑟 is the projected distance between the nodes.
The plus sign refers to the shadow zone for each projection (i.e. it is possible that one
projection is in LOS and the other in NLOS). When the link is in NLOS the plus sign applies to
both edges. For LOS conditions the edge farthest from the link is in the shadow zone (plus
sign) and the other in the LOS zone (minus sign) as shown in Figure C-7. Figure C-9 shows a
comparison of the model with measurements at 60 GHz. For multiple blocking screens see
Section 6.2 Step 7.
𝑟𝑤
𝑟ℎ
A random distribution of shadowing objects has been introduced along the streets as shown in
Figure C-8. In this case all blocking objects have a height of 2.5 m and a width of 3 m. The
shadowing model is applied between all pairs of connected nodes and for all screens. Figure
C-10 shows the corresponding result on the received power for two different TX heights of
1.5 m and 10 m. For the 10 m TX height the impact of shadowing objects is relatively small.
For the lower TX height of 1.5 m however the impact of shadowing is significant (about 7 dB
loss). Figure C-9 shows body blocking loss for a LOS link (4 m distance) at 60 GHz. The
signal strength was measured as one person was walking back and forth crossing the LOS
link. A model screen width of 30 cm is used to approximate the width of the person. It should
be noted that the severe shadowing loss shown in the figure is due to the high frequency and
that the blocking object is close to the antenna. The results shown in Figure C-10 are for
2 GHz and in average considerable longer distances to the blocking objects.
TX TX
y [m]
y [m]
RX
Figure C-8: Diffracted paths between TX and RX (left), and, relative power over the RX route
(right). Obstructing objects are shown with black dots.
Shadowing [dB]
Time [s]
Figure C-9: Body blocking loss for a LOS link (4 m distance) at 60 GHz.
Distance [m]
Figure C-10: Received power at 2 GHz over the RX route for isotropic antennas and 0 dBm
transmit power for 1.5 m RX height. The upper curve (blue) corresponds to the case with no
obstructing objects and the middle curve (red) to the case with obstructing objects when the TX
is at 10 m height and the lower curve (green) to the case when TX is at 1.5 m height.
𝑅𝐶𝑆 = 𝜋𝑅 2 . (C-9)
The power density at the scatterer, according to Figure C-11, is given by
𝑃TX 𝐺 TX
𝑆 𝑠𝑐 = 4𝜋(𝑅 2 (C-10)
1)
The power at the receiver is then given by
𝜆 2
𝑃RX = 𝑆 𝑠𝑐 𝑅𝐶𝑆 𝐺 RX (4𝜋𝑅 ) . (C-11)
2
Inserting the expression for RCS gives
𝜆𝑅 2
𝑃RX = 𝑃TX 𝐺 TX 𝐺 RX ( ) , (C-12)
8𝜋𝑅1 𝑅2
𝜆 2 2
𝑃RX = 𝑆 𝑠𝑐 𝑅𝐶𝑆 𝐺 RX (4𝜋𝑅 ) (1 − 𝛼)(1 − (𝐹ℎ1 + 𝐹ℎ2 )(𝐹𝑤1 + 𝐹𝑤2 )) . (C-13)
2
where 𝛼 is the absorption coefficient of the scatterer. In order to match the size of the screen
with the cross section of the sphere the radius is set to
𝑤ℎ
𝑅=√ . (C-14)
𝜋
As the expression is reciprocal, receive and transmit nodes may be switched depending on
which node is in the vicinity of the scatterer. Moreover, if the scatterer is located between two
intermediate nodes, e.g. diffraction corners, the shadowing loss and the scattered paths are
determined according to (C-14)) except that the receive antenna gain is used only at the RX
node i.e.
1 2 2
𝑃𝑛𝑜𝑑 = 𝑆 𝑠𝑐 𝑅𝐶𝑆 (𝑅 ) (1 − 𝛼)(1 − (𝐹ℎ1 + 𝐹ℎ2 )(𝐹𝑤1 + 𝐹𝑤2 )) . (C-15)
2
Figure C-12 shows the pathways between the TX and significant scatterers around the RX,
and, RX and significant scatterers around TX. The corresponding distributions of path angles
and path propagation distances are shown in Figure C-13.
TX TX
y [m]
y [m]
RX RX
x [m] x [m]
y [m]
y [m]
Figure C-12: Paths between TX and scatterers around one RX location (upper left) and paths
between one RX location and scatterers around TX (upper right), and, relative power over the RX
route due to scatterers around RX (lower left) and scatterers around TX (lower right).
Figure C-13: Distributions of path angles (left) and propagation distances (right) at RX for paths
between TX and scatterers around RX (lower) and paths between RX and scatterers around TX
(upper). The power is relative to the strongest path (LOS or diffracted) shown in Figure C-6.
Only scatterers which are in LOS to either TX and/or RX or in LOS to the first and the last of
three consecutive nodes are selected. Moreover, the power of the scattered wave should be
more than -40 dB relative to the strongest path.
C.1.6 Specular paths
Specular reflection is described in Section 6.2 Step 9
C.1.7 Diffuse paths
In addition to the specular paths a set of random point sources can be introduced over the
surfaces of the exterior walls accounting for the surface roughness. A uniform distribution of
point sources may be used. Moreover the distribution should be as sparse as possible in order
reduce the model complexity. Typically the density of the point sources should not be higher
than what can be resolved by the antennas used in the simulation. Details can found in
Section 6.2 Step 11.
C.1.8 Diffraction
C.1.8.1 Calculation of electric field with ray-method
When the ray interacts N times (such as reflection, transmission and diffraction) between
source and observation points, let ̅𝑾 ̅̅𝑛 (𝑛 = 1~𝑁) be a dyadic coefficient which accommodates
the ray interactions, in general, the electric field of the ith ray can be derived as
launched ray from the source to the first diffraction point. The antenna pattern of transmitter 𝑭𝒕
is introduced with gain pattern 𝐺𝜃 and 𝐺𝜙 as
Here 𝒖̂ 𝒕𝜽 , 𝒖
̂ 𝒕𝝓 are unit vectors with the origin at the source point. In this section, calculation of
diffraction is of interest and the concept of UTD (Uniform Geometrical Theory of Diffraction) is
focused.
C.1.8.2 Propagation of diffraction
Figure C-14 illustrates canonical model for wedge diffraction where the spread angle of the
wedge is (2 − 𝑛)𝜋 with the incident ray impinges at skew angle. Unit vector 𝒍̂ parallel to the
direction of the wedge is also defined as shown. Unit vectors 𝒓̂in and 𝒓̂𝐷 indicate direction of
propagation for the incident ray and the diffracted ray with respect to diffraction point D,
respectively.
These unit vectors are related based on the Fermat principle;
̅̅̅𝑛 = 𝑫
𝑾 ̅ 𝐴(𝑟𝑖𝑛 , 𝑟𝐷 ), (C-20)
where 𝐴(𝑟𝑖𝑛 , 𝑟𝐷 ) is so-called the spreading factor. For example, when the incident ray is a
spherical wave, the spreading factor is defined as
𝑟𝑖𝑛
𝐴(𝑟𝑖𝑛 , 𝑟𝐷 ) = √(𝑟 +𝑟
. (C-21)
𝑖𝑛 𝐷 ) 𝑟𝐷
Here the path length between a source point (or a former diffraction point when multiple
diffractions are considered) and a diffraction point is 𝑟𝑖𝑛 . The path length between the
diffraction point to an observation point (or a latter diffraction point when multiple diffractions
are considered) is 𝑟𝐷 . The dyadic diffraction coefficient 𝑫̅ is expressed with scalar diffraction
coefficient
𝐷𝑎 𝐷𝑏 𝒖 ̂D
𝜷
̅ = (𝒖
𝑫 ̂ in
𝜷 ̂ in
𝒖 𝝓)( ) ( D) (C-22)
𝐷𝑐 𝐷𝑑 𝒖 ̂𝝓
̂ in
Here, 𝒖 𝜷,𝒖̂ in ̂D
𝝓, 𝒖 ̂D
𝜷 , and 𝒖𝝓 are defined as
𝒓̂ ×𝒍̂
̂ in
𝒖 in
𝝓 = sin 𝛽 (C-24)
in
𝒍̂−𝒓̂D cos 𝛽D
̂D
𝒖𝜷 = (C-25)
sin 𝛽D
𝒓̂ ×𝒍̂
̂D
𝒖 𝐷
𝝓 = sin 𝛽 (C-26)
𝐷
When multiple ray interactions (e.g. reflection, transmission and diffraction) occurred N times,
the dyadic coefficients ̅𝑾
̅̅1~𝑁 should be defined respectively. After all, the coefficients are
multiplied by the inner product as in (C-16) in order to derive electric field at an observation
point.
Simple polarization model
In order to provide low complexity modelling of polarization the following approach is
proposed. The parameters of the Berg recursive model [Ber95] are generalized to account for
polarization by providing a polarization matrix for q90 i.e.
qVV VH
q90
Q90 90
HV HH
. (C-27)
q90 q90
Each element accounts for the polarization coupling between horizontal (H) and vertical (V)
polarized TX and RX antennas. If the propagation path is elevated the V-polarization is in the
direction perpendicular to the path and the horizontal direction (H). The parameter is
assumed to be independent of polarization. Literature indicates cross polarization
discrimination (XPD) values in the range 10-20 dB. For diffraction it is assumed that the XPD
is large. Lower XPD values are most likely caused by diffuse scattering and specular reflection
resulting in paths having significant elevation angles.
where 𝐿𝑤𝑒 is a constant for the exterior wall loss, 𝛼 [dB/m] is a constant for the interior wall
and floor penetration loss, and, 𝑑𝑖𝑛 is the path distance inside the building.
2) The penetration loss for path 𝑖 is determined by
𝑓 𝑁 𝑁
𝐿𝑖 = ∑𝑛=1 𝐿𝑛 + ∑𝑚=1
𝑤
𝐿𝑚 (C-29)
where 𝑁𝑓 and 𝑁𝑤 are the number of penetrated floors each having the loss 𝐿𝑛 and 𝐿𝑚
respectively. Either constant values (frequency dependent) are used for each type of wall and
floor or the more accurate definition of Appendix E.2.3, where the loss depends on thickness,
material, frequency and angle of incidence.
In the case (i) of no indoor layout, the propagation paths are determined from TX to RX as
illustrated in Figure C-15 based only on the outdoor map. Then for each pathway the
penetration loss is calculated by (C-28) based on the distance 𝑑𝑖𝑛 from a penetration point on
an exterior wall to the indoor link end. Finally per path penetration losses are multiplied to path
gains and channel matrix coefficients.
In the second case (ii) it is assumed that the indoor layout is confined and located inside a
building block. Now, different to the previous case, there may be interactions with indoor
structures also. Again propagation pathways from outdoor to indoor are determined. Exterior
walls are transparent similarly to the previous case. The walls, floor and ceiling confining the
indoor layout are transparent from outdoor-to-indoor, but not transparent from indoor-to-
outdoor. This is illustrated by Figure C-16 a where the light blue pathways would be present
also in the previous case, and the red pathways result from interaction with the outdoor space
and indoor space. When pathways are determined the penetration loss is calculated by
Equation (C-29) or alternatively Equation (C-28) and multiplied to path gains and channel
matrix coefficients.
a) b)
c)
Figure C-15: Example of determination of outdoor to indoor paths. In a) the paths are identified
with the building removed. In b) and c) the building is reintroduced in order to determine
penetration loss due to walls and floors.
b)
a)
Figure C-16: Combination of an outdoor map with an indoor layout. a) Example of determination
of outdoor to indoor paths. b) Virtual office layout (red drawing) inside a building block of
Madrid map.
a)
b)
Figure C-17: Example of determination of over roof-top paths. In a) the paths are identified in the
horizontal plane assuming the building where the access point is located is removed. The 3D
paths are blocked by the building and corresponding diffraction points are found at the roof-top
edge as shown on in b).
only for outdoor scenarios. However, the indoor scenario should be possible to model
essentially without additional effort as compared with the outdoor scenario.
Coefficient generation:
In D2D simulations, we need to generate locations for pairs of UEs for each link instead of BS
locations. This can be done by increasing the number of UEs and generating the locations as
described in the following subsection.
2𝐷 2
𝑑f = , (D-1)
𝜆
where 𝐷 denotes the largest dimension of the radiator and 𝜆 the wavelength of the transmitted
radio signal. In the past this was ensured by placing the UEs at minimum 2D distances of
35 m in case of UMa and 10 m in case of UMi. The radiator dimensions are usually
proportional to the radio frequency, but they also depend on the desired antenna
characteristics, e.g. antenna beam widths. Generally, smaller desired beam widths require
more antenna elements, which increase the antenna size. The LOS direction of departure
𝑑 𝑑 𝑎 𝑎
(𝜃LOS , 𝜙LOS ) and arrival (𝜃LOS , 𝜙LOS ) result directly from the generated BS and UE locations.
1 𝑑2D ≤ 18 m
Prob(LOS|outdoor UE) = { 18 m 18 m 𝑑2D
+ (1 − ) exp (− ) 𝑒𝑙𝑠𝑒
𝑑2D 𝑑2D 36 m
1 𝑑2D−out ≤ 18 m
Prob(LOS|indoor UE) = { 18 m 18 m 𝑑2D−out
+ (1 − ) exp (− ) 𝑒𝑙𝑠𝑒
𝑑2D−out 𝑑2D−out 36 m
Figure D-3: 2D and 3D distances for outdoor (left) and indoor (right) UEs.
1 𝑑2D ≤ 18 m
Prob(LOS|outdoor UE) = {18 m 18 m 𝑑2D
+ (1 − ) exp (− ) 𝑒𝑙𝑠𝑒
𝑑2D 𝑑2D 63 m
Assuming that ℎUE < 13 m.
Note that 𝑑2D ≥ 𝑑2D,min with 𝑑2D,min = 35 m
Prob(LOS|indoor UE)
1 𝑑2D−out ≤ 18 m
={ 18 m 18 m 𝑑2D−out
( + (1 − ) exp (− )) (1 + 𝐶(𝑑2D−out , ℎUE )) 𝑒𝑙𝑠𝑒
𝑑2D−out 𝑑2D−out 63 m
where
0 hUE < 13 𝑚
𝐶(𝑑2D−out , ℎUE ) = { hUE − 13 m 1.5
( ) 𝑔(𝑑2D−out ) 13 m ≤ hUE ≤ 23 m
10 m
and
𝑑2D−out 3 𝑑2D−out
𝑔(𝑑2D−out ) = 1.25 ∙ 10−6 ( ) exp (− )
1m 150 m
A linear approximation considering the BS antenna height yields the following general
expression for the LOS probability of outdoor UEs
1 𝑑2D ≤ 18 m
Prob(LOS|outdoor UE) = {18 m 18 m 5𝑑2D . (D-3)
+ (1 − ) exp (− ) 𝑒𝑙𝑠𝑒
𝑑2D 𝑑2D 9ℎBS +90 m
where 𝐷𝑆 is the RMS delay spread, 𝑟𝜏 the delay distribution proportionality factor, 𝑋𝑛 ~𝒰(0,1),
and cluster index 𝑛 ∈ [1, 𝑁].
Given a uniform delay distribution, the delay values 𝜏𝑛′ are drawn from the corresponding
range.
Finally, the delays are being normalised and sorted in ascending order, s.t. 𝜏1 = 0 and
𝜏𝑛−1 ≤ 𝜏𝑛 for 𝑛 ∈ [2, 𝑁].
𝜏𝑛 = 𝑠𝑜𝑟𝑡(𝜏𝑛′ − min({𝜏𝑖′ }𝑁
𝑖=1 )) (D-5)
In the case of LOS condition, additional scaling of delays is required to compensate for the
effect of LOS peak addition to the delay spread.
𝜏
𝜏𝑛𝐿𝑂𝑆 = 𝐶 𝑛 . (D-6)
𝐷𝑆
The heuristically determined scaling constant 𝐶𝐷𝑆 depends on the Ricean K-factor and is
defined as
𝐾𝐹 𝐾𝐹 2 𝐾𝐹 3
𝐶𝐷𝑆 = 0.7705 − 0.0433 (dB) + 0.0002 (dB) + 0.000017 (dB) , (D-7)
with 𝐾𝐹 being the Ricean K-factor in dB scale defined in Table 7-2. The scaled delays are not
to be used in cluster power generation.
Note that the resulting delay values are relative to the LOS delay. In case absolute delay
values are of interest, the propagation time of the LOS path needs to be added.
2 )
where 𝐶𝑆𝐹𝑛 ~𝒩(0, 𝜎𝐶𝑆𝐹 is the per cluster shadow fading (CSF) term in [dB]. Average the
power, so that the sum power of all cluster powers is equal to one, i.e.,
𝑃′
𝑃𝑛 = ∑𝑁 𝑛 ′ . (D-9)
𝑖=1 𝑃𝑖
In the case of LOS condition an additional specular component is added to the first cluster.
The power of the LOS component is given by
𝐾𝐹
𝑃1,𝐿𝑂𝑆 = 𝐾𝐹+1, (D-10)
1 𝑃′
𝑃𝑛 = 𝐾+1 ∑𝑁 𝑛𝑃′ + 𝛿(𝑛 − 1)𝑃1,𝐿𝑂𝑆 , (D-11)
𝑖=1 𝑖
where 𝛿(𝑛) is Dirac’s delta function and 𝐾𝐹 is the Ricean K-factor given in Table 7-2 and
converted to linear scale. These power values are used only in (D-10) and (D-10), but not in
(D-25).
Assign the power of each subpath 𝑚 within a cluster 𝑛 is given by
𝑃𝑛
𝑃𝑛,𝑚 = 𝑀
, (D-12)
𝐴𝑆𝐴 𝑃𝑛
𝜙̂𝑛𝑎 = √− ln ( ), (D-13)
0.7 𝐶𝐴𝑆 max({𝑃𝑖 }𝑁
𝑖=1 )
where 𝐴𝑆𝐴 denotes the azimuth spread of arrival as generated in Section 7.1.1.
In (D-13) the constant 𝐶𝐴𝑆 is a scaling factor related to the total number of clusters 𝑁 and is
given in Table D-3.
Table D-3: Scaling factors for AOA, AOD generation.
𝑁 4 5 8 10 11 12 14 15 16 19 20
𝐶𝐴𝑆 0.779 0.860 1.018 1.090 1.123 1.146 1.190 1.211 1.226 1.273 1.289
In the LOS case, constant 𝐶𝐴𝑆 is also dependent on the Ricean K-factor. Constant 𝐶𝐴𝑆 in
𝐿𝑂𝑆
Equation (D-13) is substituted by 𝐶𝐴𝑆 . An additional scaling of the angles is required in order
to compensate for the effect of the additional LOS. The heuristically determined Ricean K-
factor dependent scaling constant is given by
𝐿𝑂𝑆 𝐾𝐹 𝐾𝐹 2 𝐾𝐹 3
𝐶𝐴𝑆 = 𝐶𝐴𝑆 (1.1035 − 0.028 (dB) − 0.002 (dB) + 0.0001 (dB) ), (D-14)
where 𝐾𝐹 is the Ricean K-factor in dB scale as given in Table 7-2. More details on the
𝐿𝑂𝑆
derivation of 𝐶𝐴𝑆 and 𝐶𝐴𝑆 can be found in [3GPP13].
Finally the azimuth angles of arrival are given by
𝑋𝑛 𝜙̂𝑛𝑎 + 𝑌𝑛 + 𝜙𝐿𝑂𝑆
𝑎
for NLOS links
𝜙𝑛𝑎 = { , (D-15)
𝑋𝑛 𝜙̂𝑛𝑎 + 𝑌𝑛 − 𝑋1 𝜙̂1𝑎 − 𝑌1 + 𝜙𝐿𝑂𝑆
𝑎
for LOS links
𝐴𝑆𝐴 2 𝑎
with 𝑋𝑛 ~ 𝒰{−1, +1} being a uniformly distributed random variable and 𝑌𝑛 ~𝒩 (0, ( 7
) ). 𝜙𝐿𝑂𝑆
is the LOS direction defined in the network layout description.
Finally add offset angles 𝛼𝑚 from (D-4) to the cluster angles
𝑎
𝜙𝑛,𝑚 = 𝜙𝑛𝑎 + 𝐶𝐴𝑆𝐴 𝛼𝑚 , (D-16)
where 𝐶𝐴𝑆𝐴 is the cluster-wise RMS azimuth spread of arrival angles (cluster ASA) as given
in Table 7-2.
Table D-4: Ray offset angles within a cluster, given for 1 RMS angle spread [WIN208-D112].
Ray number 𝒎 Offset angles 𝜶𝒎
1,2 ± 0.0447
3,4 ± 0.1413
5,6 ± 0.2492
7,8 ± 0.3715
9,10 ± 0.5129
11,12 ± 0.6797
13,14 ± 0.8844
15,16 ± 1.1481
17,18 ± 1.5195
19,20 ± 2.1551
The above table assumes relatively coarse spatial resolution at transmitter/receiver. It is not
suitable for massive MIMO. More accurate way to generate the ray offsets is to use direct
sampling of the Laplacian shape as described in [FJN+15].
𝐸𝑆𝐴 𝑃
𝜃̂𝑛𝑎 = − 𝐶 ln (𝑚𝑎𝑥({𝑃𝑛 }𝑁 ), (D-17)
𝐸𝑆 𝑖 𝑖=1 )
where 𝐸𝑆𝐴 denotes the elevation spread of arrival as generated in Section 7.1.1.
In (D-17) the constant 𝐶𝐸𝑆 is a scaling factor related to the total number of clusters 𝑁 and is
𝐿𝑂𝑆
given in . More details on the derivation of 𝐶𝐸𝑆 and 𝐶𝐸𝑆 can be found in [3GPP13].
𝐿𝑂𝑆
In the LOS case, constant 𝐶𝐸𝑆 in (D-17) is substituted by 𝐶𝐸𝑆 given by:
𝐿𝑂𝑆 𝐾𝐹 𝐾𝐹 2 𝐾𝐹 3
𝐶𝐸𝑆 = 𝐶𝐸𝑆 (1.3086 + 0.0339 (𝑑𝐵) − 0.0077 (𝑑𝐵) + 0.0002 (𝑑𝐵) ) (D-18)
𝐸𝑆𝐴 2
with 𝑋𝑛 ~ 𝒰{−1, +1} being a uniformly distributed random variable and 𝑌𝑛 ~𝒩 (0, ( 7
) ).
Furthermore
where 𝐶𝐸𝑆𝐴 is the cluster-wise RMS elevation angular spread of arrival (cluster ESA) in Table
𝑎 𝑎 𝑎
7-2. Assuming that 𝜃𝑛,𝑚 is wrapped within [0°, 360°], if 𝜃𝑛,𝑚 ∈ [180°, 360°], then 𝜃𝑛,𝑚 is set to
𝑎
(360° − 𝜃𝑛,𝑚 ).
D.14.4 Generate elevation angles of departure
𝑑
The generation of elevation angles of departure 𝜃𝑛,𝑚 follows a similar procedure as the one
presented in Section D.14.3 except for (D-19) which should be replaced by (D-22)
𝐸𝑆𝐷 2
with 𝑋𝑛 ~ 𝒰{−1, +1} being a uniformly distributed random variable and 𝑌𝑛 ~𝒩 (0, ( 7
) ),
where 𝐸𝑆𝐷 denotes the elevation spread of departure as generated in Section 7.1.1.
Furthermore 𝜇offsetEOD is a function of distance and receiver height and is given in Table D-6,
and Table D-7. In addition (D-21) should be replaced by (D-23)
3
𝑑
𝜃𝑛,𝑚 = 𝜃𝑛𝑑 + 8 𝛼𝑚 10𝜇lgESD (D-23)
Where 𝜇lgESD is the mean of the ESD log-normal distribution and also a function of distance
and receiver height, cf. Table D-6, Table D-7.
𝐸𝑆𝐷 2
lg ~𝒩(𝜇lgESD , 𝜖lgESD )
1deg
with
𝑑2D ℎUE
max (−0.5, −2.1 ( ) − 0.01 ( ) + 0.765) for LOS links
𝜇lgESD ={ 1 km 1m
𝑑2D ℎUE
max (−0.5, −2.1 ( ) − 0.01 ( ) + 0.915) for NLOS links
1 km 1m
𝐸𝑆𝐷 2
lg ~𝒩(𝜇lgESD , 𝜖lgESD )
1deg
with
𝑑2D ℎUE − ℎ𝐵𝑆
max (−0.5, −2.1 ( ) + 0.01 | | + 0.75) for LOS links
𝜇lgESD ={ 1 km 1m
𝑑2D ℎUE − ℎ𝐵𝑆
max (−0.5, −2.1 ( ) + 0.01 𝑚ax ( , 0) + 0.9) for NLOS links
1 km 1m
𝑑
c. Couple randomly the azimuth angles of departure 𝜙𝑛,𝑚 with the elevation angles of
𝑑
departure 𝜃𝑛,𝑚 within a cluster 𝑛 or within a sub-cluster in the case of two strongest
clusters (Section D.18 and Table D-8).
𝜃𝜙
exp(𝑗Φ𝑛,𝑚 )
𝑇 𝜃𝜃
𝑎
𝐹𝜃,𝐺𝐶𝑆,𝑟𝑥,𝑢 (𝜃𝑛,𝑚 , 𝑎
𝜙𝑛,𝑚 ) exp(𝑗Φ𝑛,𝑚 )
𝑁𝐿𝑂𝑆 (𝑡) 𝑃𝑛 𝑀 √𝜅𝑛,𝑚
𝐻𝑢,𝑠,𝑛 = √ ∑𝑚=1 ( 𝑎 𝑎
) 𝜙𝜃
𝑀 𝐹𝜙,𝐺𝐶𝑆,𝑟𝑥,𝑢 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) exp(𝑗Φ𝑛,𝑚 ) 𝜙𝜙
exp (𝑗Φ𝑛,𝑚 )
√𝜅𝑛,𝑚
( )
𝑑 𝑑
𝐹𝜃,𝐺𝐶𝑆,tx,s (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 2𝜋 𝑎 𝑎 𝑇
∙( 𝑑 𝑑
) exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝒅𝒓𝒙,𝒖 )) , (D-25)
𝐹𝜙,𝐺𝐶𝑆,tx,s (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝑐
2𝜋 𝑑 𝑑 𝑇 𝑎 2𝜋
𝑎 𝑇
∙ exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝒅𝒕𝒙,𝒔 )) exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝒗𝒓𝒙 𝑡))
𝑐 𝑐
2𝜋 𝑑 𝑑 𝑇
∙ exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝒗𝒕𝒙 𝑡))
𝑐
where 𝐹𝜃,𝐺𝐶𝑆,𝑟𝑥,𝑢 and 𝐹𝜙,𝐺𝐶𝑆,𝑟𝑥,𝑢 are the radiation field patterns in the direction of the spherical
basis vectors, 𝒆𝜽 and 𝒆𝝓 respectively of receive antenna element 𝑢 , while 𝐹𝜃,𝐺𝐶𝑆,tx,s and
𝐹𝜙,𝐺𝐶𝑆,tx,s are the radiation field patterns in the direction of the spherical basis vectors, 𝒆𝜽 and
𝒆𝝓 respectively of transmit antenna element s. 𝒆𝒓 is the spherical unit vector defined in Section
7.2.1. 𝒅𝒓𝒙,𝒖 and 𝒅𝒕𝒙,𝒔 are the position vectors of the receive antenna element 𝑢 and transmit
antenna element s, given in GCS. 𝜆𝑐 is the wavelength of the carrier frequency. The Doppler
component in (D-25) is calculated by using the velocity vector 𝒗𝒓𝒙 , which describes the RX’s
direction of movement in GCS as well as its magnitude in velocity.
For the two strongest clusters, say n = 1 and 2, rays are spread in delay to three sub-clusters
(per cluster), with fixed delay offset {0,5,10 ns} (see Table D-5). The delays of the sub-clusters
are
𝜏𝑛,1 = 𝜏𝑛 + 0 𝑛s
𝜏𝑛,1 = 𝜏𝑛 + 5 𝑛s (D-26)
𝜏𝑛,1 = 𝜏𝑛 + 10 𝑛s
Twenty rays of a cluster are mapped to sub-clusters as presented in Table D-8 below. The
corresponding offset angles are taken from Table D-4 with mapping of Table D-8.
Table D-8: Sub-cluster information for intra cluster delay spread clusters.
Sub-cluster # Mapping to rays Power Delay offset
1 1,2,3,4,5,6,7,8,19,20 10/20 0 ns
2 9,10,11,12,17,18 6/20 5 ns
3 13,14,15,16 4/20 10 ns
In the LOS case, determine the channel coefficients by adding a single line-of-sight ray and
scaling down the other channel coefficients generated in Equation (C-21). The channel
coefficients are given by
𝐿𝑂𝑆 (𝑡) 1 𝐾
𝐻𝑢,𝑠,𝑛 = √ 𝐻 𝑁𝐿𝑂𝑆 (𝑡) +
𝐾+1 𝑢,𝑠,𝑛
𝛿(𝑛 − 1)√𝐾+1
𝑎 𝑎 𝑇
𝐹𝜃,𝐺𝐶𝑆,𝑟𝑥,𝑢 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) exp(𝑗ΦLOS ) 0
∙( ) ( )
𝑎 𝑎
𝐹𝜙,𝐺𝐶𝑆,𝑟𝑥,𝑢 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 0 −exp(𝑗ΦLOS )
𝑑 𝑑
𝐹𝜃,𝐺𝐶𝑆,tx,s (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 2𝜋 𝑎 𝑎 𝑇
∙ ( 𝑑 𝑑
) exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 𝒅𝒓𝒙,𝒖 )) , (D-27)
𝐹𝜙,𝐺𝐶𝑆,tx,s (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 𝑐
2𝜋 𝑑 𝑑 𝑇 2𝜋
𝑎 𝑎 𝑇
∙ exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 𝒅𝒕𝒙,𝒔 )) exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) 𝒗𝒓𝒙 𝑡))
𝑐 𝑐
2𝜋 𝑑 𝑑 𝑇
∙ exp (𝑗 𝜆 (𝒆𝒓 (𝜃𝑛,𝑚 , 𝜙𝑛,𝑚 ) 𝒗𝒕𝒙 𝑡))
𝑐
where 𝛿(𝑛) is the Dirac’s delta function and 𝐾 is the Ricean K-factor defined in Table 7-2 in
linear scale.
𝑑 𝑑 𝑎 𝑎
Note that 𝒆𝒓 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) = −𝒆𝒓 (𝜃𝐿𝑂𝑆 , 𝜙𝐿𝑂𝑆 ) .
Applicability
𝑑2D−in range
𝑃𝐿O2I |𝑑𝐵 = 𝑃𝐿LOS/NLOS | (𝑑3D−out + 𝑑3D−in ) + 20 + 0.5 ( )
𝑑𝐵 1m 10 m
< 𝑑2D−out
𝑑2D−in is assumed to be uniformly distributed between 0 and 25 m. + 𝑑2D−in < 1 km
Shadow fading std in [dB] is assumed to be 𝜎𝑆𝐹 = 7 according to [3GPP10-36814] 0 m < 𝑑2D−in
< 25 m
default antenna
heights
ℎ𝐵𝑆 = 10 m
1.5 m ≤ ℎUE
≤ 22.5 m
′ 𝑓
with𝑑𝐵𝑃 = 4(ℎ𝐵𝑆 − ℎ𝐸 )(ℎUE − ℎ𝐸 ) 𝑐, where 𝑐 denotes the speed of light.
𝑐
In the event that the link is determined to be LOS, ℎ𝐸 = 1 m with a
probability equal to 1 ⁄(1 + 𝐶(𝑑2𝐷 , ℎUE )) and chosen from a discrete
uniform distribution 𝒰{12, 15, … , (ℎ𝑈𝐸 − 1.5)} otherwise. The function
𝐶(𝑑2𝐷 , ℎUE )is defined in Table D-2.
Shadow fading std in [dB] is assumed to be 𝜎𝑆𝐹 = 4 according to [3GPP10-
36814]
Applicability range
𝑃𝐿NLOS |𝑑𝐵 = max(𝑃𝐿LOS |𝑑𝐵 , 𝑃𝐿A |𝑑𝐵 ) 10 m < 𝑑2𝐷 < 5 km
10 m ≤ ℎ𝐵𝑆 ≤ 150 m
𝑤street ℎ̅building 1.5 m ≤ ℎUE ≤ 22.5 m
𝑃𝐿A |𝑑𝐵 = 161.94 − 7.1 log10 ( ) + 7.5 log10 ( ) 5 m ≤ 𝑤street ≤ 50 m
1m 1m
2 5 m ≤ ℎ̅building ≤ 50 m
ℎ̅building ℎ𝐵𝑆
− (24.37 − 3.7 ( ) ) log10 ( ) default values
ℎBS 1m 𝑤street = 20 m
ℎ𝐵𝑆 𝑑3𝐷 ℎ̅building = 20 m
+ (43.42 − 3.1 log10 ( )) (log10 ( ) − 3)
1m 1m ℎ𝐵𝑆 = 25 m
𝑓𝑐 ℎUE
+ 20 log10 ( ) − 0.6 ( )
1 GHz 1m
with 𝑤street being the street width, and ℎ̅building being the average building
height. Shadow fading std in [dB] is assumed to be 𝜎𝑆𝐹 = 6 according to
[3GPP10-36814]
Applicability range
𝑑2D−in 10 m < 𝑑2D−out +
𝑃𝐿O2I |𝑑𝐵 = 𝑃𝐿LOS/NLOS | (𝑑3D−out + 𝑑3D−in ) + 20 + 0.5 ( )
𝑑𝐵 1m 𝑑2D−in < 1 km
0 m < 𝑑2D−in < 25 m
𝑑2D−in is assumed to be uniformly distributed between 0 and 25 m. default antenna heights
Shadow fading std in [dB] is assumed to be 𝜎𝑆𝐹 = 7 according to [3GPP10- ℎ𝐵𝑆 = 25 m
36814]
The formulas for the frequency dependence of the relative permittivity 𝜀𝑟′ and the conductivity
𝜎 are in [ITUR13-2040] given as
𝜀𝑟′ = 𝑎 𝑓 𝑏 (E-1)
and
𝜎 = 𝑐 𝑓𝑑 (E-2)
In addition values for the forest and concrete blocks (rubble) may be needed. Because they
actually are not mere materials, the appropriate information can be found in the literature and
should be discussed in connection of the TC10 (Emergency Communications.)
In [Saa06] there are measurement results from real asphalt roads. The author gives the
asphalt relative permittivity values 4 – 8 and even values 8 – 15, if steel slag is present, that is
rather common in many countries. Imaginary parts of the permittivity were not measured.
Metal coated glass is nowadays a common wall and window material. Such a glass is coated
by metal oxide and/or metal layers of the total depth about 100 nm. In [KKO+07] there are
theoretical results for the material showing that the coated glass is almost a perfect reflector
while the transmission coefficient is near -30 dB. In [KKO+08] the theoretical results are
validated with transmission measurements. In Table E-2 a surface conductivity is given.
Because of the multilayer construction the simple formulas should not be used. For reflection
the reflection coefficient of perfect conductor is recommended, and, for penetration the
attenuation from [KKO+07] is recommended.
According to [MPH06] the conductivity and permittivity of the human body are extremely
frequency dependent. For example, the conductivity σ increases from 1 S/m to 10 S/m and the
relative permittivity 𝜀𝑟 goes from 80 to 30 when the frequency increases from 100 MHz to
10 GHz. In this frequency range there were only few measurements so it should be utilized
with caution. However, the measurements are well in line with the other measurements in the
lower frequencies.
σ
where 𝜀𝑟 is the complex relative permittivity, 𝜀𝑟 = 𝜀𝑟′ − j𝜀𝑟′′ = 𝜀𝑟′ − j ω𝜀′ , σ is the conductivity of
0
the material where the wave is arriving from the air, 𝜔 is the angular frequency, 𝜀0 is the
permittivity of the vacuum and 𝛾 is the angle between the normal of the surface between air
and the other material and the direction of the wave.
where 𝑍in is the input impedance and 𝑍 is the characteristic impedance (of air in this case).
The input impendence in the case of a single layer dielectric material when the wall material is
lossy is calculated by
where 𝑍𝑤 is the characteristic impedance of the medium, 𝑑 is the width of the dielectric slab
and 𝛿 = 𝑘0 √𝜀𝑟 − sin2 𝛾, 𝜀𝑟 is the complex permittivity (𝜀𝑟′ − j𝜀𝑟′′ ), 𝛾 is the incident angle and
2π
𝑘0 = 𝜆 is the wave number (here, λ is the wave length in the vacuum).
Wall
Reflected Transmitted
Normal
Incident
d
Figure E-1: Reflection from a wall and penetration through a wall.
With positive real and negative imaginary parts, the wave impedances for transverse electric
TE and transverse magnetic TM will also have complex values such that
TE 𝜂
𝑍𝑤 = , (E-7)
√𝜀𝑟′ −j𝜀𝑟′′ −sin2 𝜃
TM 𝜂
𝑍𝑤 = 𝜀′ −j𝜀′′ √𝜀𝑟′ − j𝜀𝑟′′ − sin2 𝜃 , (E-8)
𝑟 𝑟
𝜂
𝑍 TE = , (E-9)
√1−sin2 𝜃
𝑍 TM = √1 − sin2 𝜃 , (E-10)
𝑍 −𝑍TE
𝛤TE = 𝑍in +𝑍TE , (E-12)
in
𝑍 −𝑍TM
𝛤TM = 𝑍in +𝑍TM . (E-13)
in
The specific values of complex permittivity of particular materials are shown in the Table E-1
and E-2. Note that the imaginary part of the complex permittivity is calculated by
𝜎
𝜀𝑟′′ = 17.98 𝑓GHz
, (E-14)
where 𝜎 is the conductivity.
1
TE polarisation
TM polarisation
0.8
0.6
||
0.4
0.2
0
0 1 2
10 10 10
Frequency [GHz]
Figure E-2: Reflection coefficient for concrete (45 deg incident angle).
The loss through a wall has been discussed in [ITUR13-2040] and consists of the
transmission through the air – wall surface, the attenuation due to the finite conductivity of the
wall material and the transmission through the wall – air surface on the opposite side. In
addition it is assumed that the effect of surface roughness in both of the surfaces affects
similarly as in the reflection.
If we assume that the building material is a homogeneous dielectric single-layer plate with a
smooth surface, with the geometry shown in Figure E-1, we can express the transmission
coefficient, 𝑇, of the building material as given in [ITUR13-2040]:
(1−𝛤 2 )exp(−𝑗(𝛿−𝑘0 𝑑)
𝑇= 1−𝛤 2 exp(−𝑗2𝛿)
, (E-15)
2π
where 𝛿 = 𝑘0 √𝜀𝑟 − sin2 𝛾 , 𝑘0 = 𝜆 𝑑 is the thickness of the wall, λ is the wave length in the
vacuum, 𝜀𝑟 is the complex permittivity ( 𝜀𝑟′ − j𝜀𝑟′′ ) and 𝛾 is the angle between the normal of the
wall surface and the direction of the propagation in the air of the incoming radio wave.
The equation is applicable for both the TE and the TM waves when 𝛤 means the Fresnel
reflection coefficient for the TE and TM waves (which are defined as (E-3) and (E-4)),
respectively. Figure E-3 shows the penetration loss for a concrete wall of 5, 10 or 15 cm as a
function of frequency for three different angles of incidence. The material parameter values
are taken from the Table E 1.
Polarization TE, Thickness 5 (top), 10 and 15 cm (bottom) Polarization TM, Thickness 5 (top), 10 and 15 cm (bottom)
0 0
Angle 0 Angle 0
-10 30 -10 30
60 deg 60 deg
-20 -20
-30 -30
-40 -40
Loss [dB]
Loss [dB]
-50 -50
-60 -60
-70 -70
-80 -80
-90 -90
-100 -100
0 10 20 30 40 50 60 70 80 90 100 0 10 20 30 40 50 60 70 80 90 100
Frequency [GHz] Frequency [GHz]
Figure E-3: Penetration loss for a concrete wall of 5, 10 or 15 cm as a function of frequency for
three different angles of incidence.
Here 𝜆 is the wave length in free space and 𝜃𝑖 is the angle of incidence for the incoming wave.
In the Gaussian rough surface scattering model the height deviations are assumed to follow
the normal distribution. We can calculate the coefficient 𝜌𝑆 due to surface roughness that is
used to multiply the result for smooth surface reflection to take the surface roughness into
account [LFR96].
𝜋𝜎ℎ 𝑐𝑜𝑠𝜃𝑖 2
𝜌𝑆 = exp [−8 ( 𝜆
) ], (E-17)
where 𝜎ℎ is the standard deviation of the zero mean height of the surface and other notations
are like in the previous formula.
In [LFR93] and [LFR96] reflection from three different walls consisting of different materials,
brick wall, limestone and metal coated glass, are measured and compared with theoretical
values.
Two wall materials, brick and metal coated glass were smooth. Measurements were
conducted at 4 GHz for parallel and perpendicular polarisations. Agreement of the results with
Fresnel reflection formulas was good for both materials.
The third wall consisted of limestone blocks and had a very rough surface. It was measured at
two frequencies 1.9 and 4 GHz and for parallel and perpendicular polarisations. The
measured reflection coefficients were smaller than for the Fresnel reflection. However the
coefficients were greater than predicted by Gaussian rough surface scattering model that
takes into account the surface roughness. The authors propose in [LFR96] their own formula
that fits better in their measurement results. The method is simply to take an average of the
Fresnel and Gaussian rough surface results.
It seems that for the case of the rough surface of the limestone the reflection coefficient is
rather randomly scattered around the expected value at each incident angle. In such a case
the proposed method in [LFR96] would certainly give better results. In the channel model this
could be utilised by introducing random fluctuation on the reflection coefficients, when rough
surfaces are modelled.
E.3 Diffraction
E.3.1 Diffraction from knife edge
The power that the radio waves carry from a transmitter to a receiver can be understood to
move through an ellipsoid, named Fresnel ellipsoid. The size of the ellipsoid is defined by the
wave-length of the radio signal and the distances from the transmitter and the receiver to the
observation point. The formula is expressed below for the radius of the first Fresnel zone
(cross section of the first Fresnel ellipsoid). There is an infinite number of Fresnel zones of
higher order, but their importance is marginal [Ber00].
Radius of the first Fresnel zone 𝑟𝑓 is:
𝜆𝑑1 𝑑2
𝑟𝑓 = √𝑑 , (E-18)
1 +𝑑2
where 𝑑1 is the distance from the TX to the diffraction point 𝑑2 is the distance from the
diffraction point to the RX and 𝜆 is the wavelength.
If the first Fresnel zone is free of obstacles, the radio wave is received without extra loss. If the
first Fresnel zone is obstructed by a knife-edge, the signal is attenuated as indicated by the
factor 𝐹(𝜈) as given in [Par00].
2
1+𝑗 ∞ −𝑗𝜋𝑡
𝐹(𝜈) = ∫ 𝑒 2 𝑑𝑡, (E-19)
2 𝜈
2 1 1
𝜈 = ℎ √𝜆 (𝑑 + 𝑑 ). (E-20)
1 2
There is a good approximation for the 𝐹(𝜈) defined in decibels, cf. [ITUR13-52613]:
Both functions 𝐽(𝜈) and 𝐹(𝜈) are illustrated in dB in Figure E-4. Parameters h, d1 and d2 are
illustrated below.
In [MH13] there is an approximation for the attenuation caused by one knife-edge. It is used in
the METIS map-based channel mode and the equation can be found in Section 6.2 Equation
(6-6).
It is also possible to calculate the knife-edge diffraction using UTD. The general formula is
given in Section 6.2 Equation (6-30). Letting the wedge angle go to zero we get the knife-edge
situation. This equation can model the two different polarisations. However, the knife-edge
equations can be calculated faster and this may justify their usage.
| [dB]
0
freespace
|E / Efreespace
Green: 20 GHz Green: 20 GHz
0.75 m Magenta: 40 GHz
Magenta: 40 GHz -20
-20
1.636 m
Tx -30 -30
10
20log10
-50 -50
-60 -60
0 2 4 6 8 10 0 2 4 6 8 10
• The heights of Tx and Rx are 0.75 m. d [m]
d [m]
Figure E-6 illustrates frequency dependency of wedge diffraction. For example, the relative
power level at d = 10 m shows about 16 dB differences between 1 and 40 GHz for both
polarisations in this specific scenario. Map-based model utilizes the capability of UTD
calculation with frequency dependent material properties. The details of implementation are
described in Section 6 and Appendix C.
Table E-3: Attenuation caused by oxygen and water vapour (rounded maximum values visually
estimated from [ITUR13-67610].
Frequency attenuation (dB/km) influence
22 GHz 0.2 water vapour
60 GHz 15 oxygen
𝛾𝑅 = 𝑘𝑅 𝛼 , (E-22)
in which
2
log10 𝑓−𝑏𝑗
log10 𝑘 = ∑4𝑗=1 𝑎𝑗 exp [− ( 𝑐𝑗
) ] + 𝑚𝑘 log10 𝑓 + 𝑐𝑘 , (E-23)
and
2
log10 𝑓−𝑏𝑗
𝛼 = ∑5𝑗 𝑎𝑗 exp [− ( 𝑐𝑗
) ] + 𝑚𝛼 log10 𝑓 + 𝑐𝛼 , (E-24)
where
𝑓 : frequency (GHz)
𝑘 : either 𝑘𝐻 or 𝑘𝑉
𝛼 : either 𝛼𝐻 or 𝛼𝑉
The parameter values are available from [ITUR05-8383].
Based on the above equations and the parameter values from the ITU-R recommendation, we
can calculate frequency dependency of rain attenuation (see Figure E-7). These curves show
the extra attenuation caused by rain and should not be mixed up with pathloss model.
V-pol, 1 mm/h
H-pol, 1 mm/h
30 V-pol, 3 mm/h
H-pol, 3 mm/h
V-pol, 10 mm/h
H-pol, 10 mm/h
25 V-pol, 30 mm/h
H-pol, 30 mm/h
V-pol, 100 mm/h
attenuation [dB/km]
15
10
0
0 10 20 30 40 50 60 70 80 90 100
frequency [GHz]
𝑃𝐿𝑂𝑆
𝐾𝐹 = ∑𝑖≠𝐿𝑂𝑆 𝑃𝑖
(F-1)
This is the simplest way to calculate K-factor. This method was used in the data analysis of
60 GHz measurements (Section 5.12 and Appendix A.5.5) and the analysis of simulation
results (Section 9.3)
(𝐸[𝑃𝑁𝐵 ])2
𝑚= (F-2)
𝐸[(𝑃𝑁𝐵 −𝐸[𝑃𝑁𝐵 ])2 ]
where PNB presents the power of narrowband signal. The corresponding Ricean K-factor can
be approximated as 𝑚 = 0.4998528 𝐾 + 0.7622159, when 𝐾 ≥ 2 [AS96], which leads to K-
factor approximation
𝑚−0.7622159
𝐾≈ (F-3)
0.4998528
This method is valid only for narrowband K-factor calculations, i.e. calculated values are
appropriate for frequency segments narrower than the correlation bandwidth Bc of the
channel, when averaging period is longer then correlation time of the short-term fluctuations
[GME99]. For selective fading channels, separate analysis should be performed per frequency
bins narrower than Bc. The Moment-method was used in the data analysis in the Section 5.5
and Appendix A.1.3. It should be noted that the K-factor is used only for LOS channel.
However, in theory it is possible to obtain K-factor also for NLOS channel by using the MoM.
References
[3GPP03-25996] 3GPP TR 25.996, “Spatial channel model for multiple input multiple
output (MIMO) simulations”, V6.1.0, September 2003.
[3GPP10-36814] 3GPP TR 36.814, “Further advancements for E-UTRA physical layer
aspects”, V9.0.0, March 2010.
[3GPP13] 3GPP R1-135895, “Remaining issues of EoD/EoA modeling”,
Qualcomm Incorporated, RAN1#75, San Francisco, USA, November
2013.
[3GPP14] 3GPP R1-143352, “Spatial Correlation Accuracy with WINNER II Large-
scale Parameter Generation”, 3GPP TSG RAN1#78, 2014.
[3GPP14-36843] 3GPP TR 36.843, “Study on LTE Device to Device Proximity Services”,
V12.0.1, March 2014.
[3GPP14-36873] 3GPP TR 36.873, “Study on 3D channel model for LTE”, V12.0.0,
September 2014.
[ACL12] M. Alazah, S. N. Chandler-Wilde and S. L. Porte. "Computing Fresnel
integrals via modified trapezium rules", in eprint arXiv: 1209.3451
Numerische Mathematik, pp. 1-27, 2012.
[ALS+14] M. R. Akdeniz, Y. Liu, M. K. Samimi, S. Sun, S. Rangan, T. Rappaport
and E. Erkip, “Millimeter wave channel modeling and cellular capacity
evaluation”, IEEE J. Sel. A. Commun., no.6, pp.1164-1179, Sept. 2014.
[AP09] P. Agrawal and N. Patwari, “Correlated link shadow fading in multihop
wireless networks”, IEEE Trans. Wireless Commun., vol. 8, no. 8, pp.
4024–4036, Aug. 2009.
[AS96] S. A. Abbas and A. U. Sheikh, “A Geometric Theory of Nakagami
Fading Multipath Mobile Radio Channel with Physical Interpretations”,
IEEE 46th Vehicular Technology Conference, Atlanta, USA, Vol. 2, pp.
637–641, April 1996.
[Bal89] C. A. Balanis, “Advanced engineering electromagnetics”, Vol. 20, New
York, Wiley, 1989.
[BE10] H. Buddendick and T. F. Eibert, "Parallelized physical optics
computations for scattering center models in radio channel simulations",
in Proc. IEEE Vehicular Networking Conference (VNC), 2010.
[Ber00] H. L. Bertoni, “Radio Propagation for Modern Wireless Systems”,
Prentice Hall PTR, Upper Saddle River, HJ 07458, 2000.
[Ber95] J.-E. Berg, "A recursive method for street microcell pathloss
calculations", Sixth IEEE International Symposium on Personal, Indoor
and Mobile Radio Communications, PIMRC'95, Wireless: Merging onto
the Information Superhighway. Vol. 1, 1995.
[BV11] M. Boban and T. T. V. Vinhoza, “Modeling and Simulation of Vehicular
Networks: towards Realistic and Efficient Models”, Mobile Ad-Hoc
Networks: Applications, Prof. Xin Wang (Ed.), January 2011.
[BVF+11] M. Boban, T. T. V. Vinhoza, M. Ferreira, J. Barros, and O. K. Tonguz,
"Impact of vehicles as obstacles in vehicular ad hoc networks", IEEE
Journal on Selected Areas in Communications, Vol 29, no 1, pp 15-28,
2011.
[ITUR05-8383] ITU-R P.838-3, “Specific attenuation model for rain for use in prediction
methods”, International Telecommunication Union,
Radiocommunication Sector, 2005.
[ITUR09-2135] ITU-R M.2135-1, “Guidelines for evaluation of radio interface
technologies for IMT-Advanced”, International Telecommunication
Union (ITU), Geneva, Switzerland, Technical Report, December 2009.
[ITUR13-52613] ITU-R P.526-13, “Propagation by Diffraction”, International
Telecommunication Union, Radiocommunication Sector, 2013.
[ITUR13-2040] ITU-R P.2040, “Effects of building materials and structures on
radiowave propagation above about 100 MHz”, International
Telecommunication Union Radiocommunication Sector ITU-R, 2013.
[ITUR13-67610] ITU-R P.676-10, “Attenuation by atmospheric gases”, International
Telecommunication Union Radiocommunication Sector ITU-R, 2013.
[Jal10] N. Jalden, “Analysis and Modelling of Joint Channel Properties from
Multi-site, Multi-Antenna Radio Measurements”, PhD thesis, Royal
Institute of Technology KTH Sweden, 2010.
[JH14] J. Järveläinen and K. Haneda, “Sixty gigahertz indoor radio wave
propagation prediction method based on full scattering model”, Radio
Science, vol. 49, no. 4, pp. 293-305, 2014.
[JK15] T. Jämsä and P. Kyösti, “Device-to-device extension to geometry-based
stochastic channel models”, submitted to EuCAP, Lisbon, Portugal,
2015.
[JKK+14] J. Järveläinen, M. Kurkela, A. Karttunen, K. Haneda and J. Putkonen,
“70 GHz Radio Wave Propagation Prediction in a Large Office”, in Proc.
10th Loughborough Antenna & Propagation Conference (LAPC2014),
Loughborough, UK, 2014.
[JPM+11] M. Jacob, S. Priebe, A. Maltsev, V. Erceg and T. Kürner, “A ray tracing
based stochastic human blockage model for the IEEE 802.11 ad
60 GHz channel model”, in Proc. 5th European Conference on
Antennas and Propagation (EUCAP), IEEE, 2011.
[JRB+14a] S. Jaeckel, L. Raschkowski, K. Börner, L. Thiele, F. Burkhardt and E.
Eberlein, "QuaDRiGa- Quasi Deterministic Radio Channel Generator,
User Manual and Documentation", Fraunhofer Heinrich Hertz Institute,
Tech. Rep. v1.2.0-283, 2014.
[JRB+14b] Jaeckel, S.; Raschkowski, L.; Börner, K.; Thiele, L., “QuaDRiGa: A 3-D
Multicell Channel Model with Time Evolution for Enabling Virtual Field
Trials,” IEEE Transactions on Antennas and Propagation, 2014.
[KKO+07] G. I. Kiani, A. Karlsson, L. G. Olsson and K. P. Esselle, “Glass
characterization for designing frequency selective surfaces to improve
transmission through energy saving glass windows,” in Proc. Asia-
Pacific Microwave Conference, (APMC), 2007.
[KKO+08] G. I. Kiani, A. Karlsson, L. G. Olsson and K. P. Esselle, “Transmission
Analysis of Energy Saving Glass Windows for the Purpose of Providing
FSS Solutions at Microwave Frequencies”, in Proc. Antennas and
Propagation Society International Symposium, 2008.
[KKP+13] W. Keusgen, A. Kortke, M. Peter and R. Weiler, “A highly flexible digital
radio testbed and 60 GHz application examples”, in European
Microwave Conference (EuMC), pp. 740–743, 2013.
METIS Public 198
Document: FP7-ICT-317669-METIS/D1.4
Date: 24/02/2015 Security: Public
Status: Final Version: 1.0
[MCN14] A.O. Martinez, E. De Carvalho and J.O. Nielsen, “Towards very large
aperture massive MIMO: a measurement based study,” in Proc. Global
Commun. Conf. (GLOBECOM 2014), Austin, TX, December 2014.
[MEE+14] L. Minghini, R. D’Errico, V. D. Esposti and E. M. Vitucci,
"Electromagnetic simulation and measurement of diffuse scattering from
building walls", in Proc. 8th European Conference on Antennas and
Propagation (EuCAP), 2014.
[MEP10] A. Maltsev, V. Erceg and E. Perahia, “Channel models for 60 GHz
WLAN systems”, Document IEEE 802.11-09/0334r8, 2010.
[MET13-D11] METIS D1.1, “Future radio access scenarios, requirements and KPIs”,
V1.0, ICT-317669, METIS project, 2013.
[MET13-D51] METIS D5.1, “Intermediate description of the spectrum needs and
usage principles”, V1.0, ICT-317669, METIS project, August 2013.
[MET13-D61] METIS D6.1, “Simulation guidelines,” V1.0, ICT-317669, METIS project,
2013.
[MET14-D12] METIS D1.2, “Initial channel models based on measurements”, V1.0,
ICT-317669, METIS project, 2014.
[MET15-D65] METIS D6.5, “Report on simulation results and evaluations”, V1.0, ICT-
317669, METIS project, 2015.
[MH13] J. Medbo and F. Harrysson, "Channel modeling for the stationary UE
scenario", in Proc. 7th European Conference on Antennas and
Propagation (EuCAP), pp. 2811-2815, 2013.
[MKA02] H. Masui, T. Kobayashi and M. Akaike, “Microwave path-loss modeling
in urban line-of-sight environments”, IEEE J. Sel. Areas Commun., vol.
20, no. 6, pp. 1151-1155, August 2002.
[Mol04] A. F. Molisch, “A generic channel model for MIMO wireless propagation
channels in macro- and microcells”, IEEE Trans. Signal Processing, vol.
52, no. 1, pp. 61–71, January 2004.
[MPH06] D. Miklavčič, N. Pavšelj and F. X. Hart, "Electric properties of tissues",
Wiley encyclopedia of biomedical engineering, DOI:
10.1002/9780471740360.ebs0403, Accessed 04.04.2006.
[Par00] J. D. Parsons, “The Mobile Radio Propagation Channel”, John Wiley &
Sons, 2000.
[PB04] P. Pongsilamanee and H. L. Bertoni, "Specular and nonspecular
scattering from building facades", IEEE Transactions on Antennas and
Propagation, Vol 52, no 7, pp. 1879-1889, 2004.
[PSH+11] J. Poutanen, J. Salmi, K. Haneda, V. Kolmonen and P. Vainikainen,
“Angular and shadowing characteristics of dense multipath components
in indoor radio channels”, IEEE Transactions on Antennas and
Propagation, Vol 59, no 1, 245-253, 2011.
[Qin02] Q. Qineti, “A Generic Model of 1-60 GHz Radio Propagation through
Vegetation”, Final Report, 2002.
[RBM+12] T. Rappaport, E. Ben-Dor, J. Murdock and Q. Yijun, “38 GHz and
60 GHz angle-dependent propagation for cellular & peer-to-peer
wireless communications”, Proc. IEEE Int. Conf. Commun. (ICC 2012),
pp. 4568-4573, June 2012.