Channel Model Tutorial
Channel Model Tutorial
V1.0 February 21, 2007 Please send comments/corrections/feedback to Raj Jain, [email protected] Please send comments to [email protected]
V1 Created on 2/21/2007
TABLE CONTENTS
CHANNEL MODELS: A TUTORIAL..................................................................................................................................3 CHANNEL MODELS: A TUTORIAL..................................................................................................................................3 A.1 BASIC CONCEPTS .............................................................................................................................................................3 A.1.1 Channel ....................................................................................................................................................................3 A.1.2 Path Loss..................................................................................................................................................................4 A.1.3 Shadowing ................................................................................................................................................................4 A.1.4 Multipath ..................................................................................................................................................................5 A.1.5 Tapped Delay Line Model ........................................................................................................................................7 A.1.6 Doppler Spread ........................................................................................................................................................7 A.2 EMPIRICAL PATH LOSS MODELS ......................................................................................................................................8 A.2.1 Hata Model...............................................................................................................................................................8 A.2.2 COST 231 Extension to Hata Model ........................................................................................................................9 A.2.3 COST 231-Walfish-Ikegami Model ..........................................................................................................................9 A.2.4 Erceg Model ...........................................................................................................................................................12 A.2.5 Stanford University Interim (SUI) Channel Models...............................................................................................14 A.2.6 ITU Path Loss Models............................................................................................................................................18 REFERENCES ......................................................................................................................................................................21
Table of Figures
FIGURE A.1.1: CHANNEL FIGURE A.1.2: PATH LOSS, SHADOWING, AND MULTIPATH [GOLDSMITH2005] FIGURE A.1.3: SHADOWING FIGURE A.1.4: MULTIPATH FIGURE A.1.5: MULTIPATH POWER DELAY PROFILE FIGURE A.1.6: TAPPED DELAY LINE MODEL FIGURE A.2.1: PARAMETERS OF THE COST-231 W-I MODEL [MOLISCH2005] FIGURE A.2.2: STREET ORIENTATION ANGLE [CICHON] FIGURE A.2.5.1: GENERIC STRUCTURE OF SUI CHANNEL MODELS 3 6 5 5 6 7 10 11 15
List of Tables
TABLE A.1.1: TYPICAL DOPPLER SPREADS AND COHERENCE TIMES FOR WIMAX [ANDREWS2007] 8 TABLE A.2.1: PARAMETERS OF THE ERCEG MODEL 13 TABLE A.2.5.1: TERRAIN TYPE AND DOPPLER SPREAD FOR SUI CHANNEL MODELS 14 TABLE A.2.5.1: SCENARIO FOR SUI CHANNEL MODELS 15 TABLE A.2.5.1: SUI 1 CHANNEL MODEL 16 TABLE A.2.5.2: SUI 2 CHANNEL MODEL 16 TABLE A.2.5.3: SUI 3 CHANNEL MODEL 16 TABLE A.2.5.4: SUI 4 CHANNEL MODEL 17 TABLE A.2.5.5: SUI 5 CHANNEL MODEL 17 TABLE A.2.5.6: SUI 6 CHANNEL MODEL 18 TABLE A.2.6.1: ITU CHANNEL MODEL FOR INDOOR OFFICE 18 TABLE A.2.6.2: ITU CHANNEL MODEL FOR OUTDOOR TO INDOOR AND PEDESTRIAN TEST ENVIRONMENT 19 TABLE A.2.6.3: ITU CHANNEL MODEL FOR VEHICULAR TEST ENVIRONMENT 19 TABLE A.2.6.4: PERCENTAGE OCCURRENCE AND ASSOCIATED RMS DELAY SPREAD FOR ITU CHANNEL MODELS 19
Page 2 of 21
V1 Created on 2/21/2007
Channel
Subscriber Station
The characteristics of wireless signal changes as it travels from the transmitter antenna to the receiver antenna. These characteristics depend upon the distance between the two antennas, the path(s) taken by the signal, and the environment (buildings and other objects) around the path. The profile of received signal can be obtained from that of the transmitted signal if we have a model of the medium between the two. This model of the medium is called channel model. In general, the power profile of the received signal can be obtained by convolving the power profile of the transmitted signal with the impulse response of the channel. Convolution in time domain is equivalent to multiplication in the frequency domain. Therefore, the transmitted signal x, after propagation through the channel H becomes y: y(f)=H(f)x(f)+n(f) Here H(f) is channel response, and n(f) is the noise. signal frequency f. Note that x, y, H, and n are all functions of the
The three key components of the channel response are path loss, shadowing, and multipath as explained below.
Page 3 of 21
V1 Created on 2/21/2007
A.1.2 Path Loss The simplest channel is the free space line of sight channel with no objects between the receiver and the transmitter or around the path between them. In this simple case, the transmitted signal attenuates since the energy is spread spherically around the transmitting antenna. For this line of sight (LOS) channel, the received power is given by: Gl2 Pr=Pt 4d Here, Pt is the transmitted power, Gl is the product of the transmit and receive antenna field radiation patterns, is the wavelength, and d is the distance. Theoretically, the power falls off in proportion to the square of the distance. In practice, the power falls off more quickly, typically 3rd or 4th power of distance. The presence of ground causes some of the waves to reflect and reach the transmitter. These reflected waves may sometime have a phase shift of 180 and so may reduce the net received power. A simple two-ray approximation for path loss can be shown to be: Gt Gr ht2 hr2 d4 Here, ht and hr are the antenna heights of the transmitter and receiver, respectively. Note that there are three major differences from the previous formula. First, the antenna heights have effect. Second, the wavelength is absent and third the exponent on the distance is 4. In general, a common empirical formula for path loss is: Pr = Pt
d Pr = Pt P0 0 d
Where P0 is the power at a distance d 0 and is the path loss exponent. The path loss is given by:
d PL(d ) dB = PL(d 0 ) + 10 log d0 Here PL(d 0 ) is the mean path loss in dB at distance d 0 . The thick dotted line in Figure A.1.2 shows the received power as a function of the distance from the transmitter.
A.1.3 Shadowing
If there are any objects (such buildings or trees) along the path of the signal, some part of the transmitted signal is lost through absorption, reflection, scattering, and diffraction. This effect is called shadowing. As shown in Figure A.1.3, if the base antenna were a light source, the middle building would cast a shadow on the subscriber antenna. Hence, the name shadowing.
Page 4 of 21
V1 Created on 2/21/2007
Here is a normally (Gaussian) distributed random variable (in dB) with standard deviation . represents the effect of shadowing. As a result of shadowing, power received at the points that are at the same distance d from the transmitter may be different and have a lognormal distribution. This phenomenon is referred to as lognormal shadowing.
A.1.4 Multipath
The objects located around the path of the wireless signal reflect the signal. Some of these reflected waves are also received at the receiver. Since each of these reflected signals takes a different path, it has a different amplitude and phase.
Page 5 of 21
V1 Created on 2/21/2007
Depending upon the phase, these multiple signals may result in increased or decreased received power at the receiver. Even a slight change in position may result in a significant difference in phases of the signals and so in the total received power. The three components of the channel response are shown clearly in Figure A.1.4. The thick dashed line represents the path loss. The lognormal shadowing changes the total loss to that shown by the thin dashed line. The multipath finally results in variations shown by the solid thick line. Note that signal strength variations due to multipath change at distances in the range of the signal wavelength.
Path loss Shadow + Path loss Multipath + Shadowing + Path loss Log(Pr/Pt)
Log(d)
Figure A.1.4: Path loss, shadowing, and Multipath [Goldsmith2005]
Since different paths are of different lengths, a single impulse sent from the transmitter will result in multiple copies being received at different times as shown in Figure A.1.5
Power Transmitted
Power Received
Delay
Delay
Page 6 of 21
V1 Created on 2/21/2007
The maximum delay after which the received signal becomes negligible is called maximum delay spread max. A large max indicates a highly dispersive channel. Often root-mean-square (rms) value of the delay-spread rms is used instead of the maximum.
h(t , ) = ci (t ) ( i )
i =1
Note that the impulse response h varies with time t. The coefficients ci(t) vary with time. There are N coefficients in the above model. The selection of the N and delay values i depends upon what is considered a significant level. This model represents the channel by a delay line with N taps. For example, the channel shown in Figure A.1.5 can be represented by a 4-tap model as shown in Figure A.1.6. c1 Delay Line 1 2 c2 c3 c4
If the transmitter, receiver, or even the other objects in the channel move, the channel characteristics change. The time for which the channel characteristics can be assumed to be constant is called coherence time. This is a simplistic definition in the sense that exact measurement of coherence time requires using the autocorrelation function. For every phenomenon in the time domain, there is a corresponding phenomenon in the frequency domain. If we look at the Fourier transform of the power delay profile, we can obtain the frequency dependence of the channel characteristics. The frequency bandwidth for which the channel characteristics remain similar is called coherence bandwidth. Again, a more strict definition requires determining the autocorrelation of the channel characteristics. The coherence bandwidth is inversely related to the delay spread. The larger the delay spread, less is the coherence bandwidth and the channel is said to become more frequency selective.
Page 7 of 21
V1 Created on 2/21/2007
the channel. The Doppler power spectrum is nonzero for (f-fD, f+fD), where fD is the maximum Doppler spread or Doppler spread. The coherence time and Doppler spread are inversely related: CoherenceTime 1 Doppler Spread
Thus, if the transmitter, receiver, or the intermediate objects move very fast, the Doppler spread is large and the coherence time is small, i.e., the channel changes fast. Table A.1.1 lists typical values for the Doppler spread and associated channel coherence time for two WiMAX frequency bands. Note that at high mobility, the channel changes 500 times per second, requiring good channel estimation algorithms
Table A.1.1: Typical Doppler Spreads and Coherence Times for WiMAX [Andrews2007] Carrier Freq Speed Max Doppler Spread Coherence Time 2.5 GHz 2 km/hr 4.6 Hz 200 ms 2.5 GHz 45 km/hr 104.2 Hz 10 ms 2.5 GHz 100 km/hr 231.5 Hz 4 ms 5.8 GHz 2 km/hr 10.7 Hz 93 ms 5.8 GHz 45 km/hr 241.7 Hz 4 ms 5.8 GHz 100 km/hr 537 Hz 2 ms
PL ,urban (d )dB = 69.55 + 26.16log10( f c ) 13.82log10(ht ) a(hr ) + (44.9 6.55log10( ht ))log10( d ) Here, f c is the carrier frequency, ht is the height of the transmitting (base station) antenna, hr is the height of the receiving (mobile) antenna, and a (hr ) is a correction factor for the mobile antenna height based on the size of the coverage area.
Page 8 of 21
V1 Created on 2/21/2007
Hata model well approximates the Okumura model for distances greater than 1km. This model is intended for large cells with BS being placed higher than the surrounding rooftops. Both models are designed for 150-1500 MHz and are applicable to the first generation cellular systems. They may not work well for WiMAX systems with smaller cell sizes and higher frequencies.
COST 321-Hata model is designed for large and small macro-cells, i.e., base station antenna heights above rooftop levels adjacent to base station.
Page 9 of 21
10
V1 Created on 2/21/2007
This model distinguishes between the line-of-sight (LOS) and non-line-of-sight (NLOS) cases. For LOS, the total path loss is: PL dB = 42.6 + 26 log(d ) + 20 log( f c ) for d 0.02 km. Here d is in units of kilometers, and f c is the carrier frequency in MHz. For the non-LOS case, path loss consists of three terms: the free space path loss L0 , the multi-screen loss Lmsd along the propagation path, and attenuation from the last roof edge to the MS, Lrts (rooftopto-street diffraction and scatter loss):
P + Lrts + Lmsd PL dB = L 0 PL 0
The free space path loss is: PL 0 dB = 32.4 + 20 log d + 20 log f c Ikegami derived the diffraction loss Lrts as: Lrts = 16.9 10 log w + 10 log f c + 20 log hm + Lori Here w is the width of the street in meters, and hm is the difference between the building height hRoof and the height of the MS hm hm = hRoof hm Orientation of the street is taken into account by an empirical correction factor Lori :
Page 10 of 21
11
V1 Created on 2/21/2007
Here is the angle between the street orientation and the direction of incidence in degrees as shown in Figure A.2.2. While Lrts expression is from Ikegami, the expressions for Lori are different from those proposed by Ikegami.
The multi-screen loss Lmsd is obtained by modeling building edges as screens. The multi-screen loss is [Walfish1988]: Lmsd = Lbsh + ka + kd log d + k f log f c 9 log b Here b is the distance between two buildings (in meters), and:
18log(1 + hb ) Lbsh = 0
for for
hb = hb hRoof and hb is the height of the BS. The path loss depends on frequency and distance as given via the
Page 11 of 21
12
V1 Created on 2/21/2007
parameters kd and k f :
18 kd = 18 15hb /hRoof
fc 0.7 925 1 for medium size cities kf = suburban areas with average vegetation density f 1.5 c 1 for metropolitan areas 925 The validity range for this model is: Carrier frequency f c Height of BS antenna hb Height of MS antenna hm Distance d 8002,000 MHz 450m 13m 0.025km
The model assumes a Manhattan street grid (streets intersecting at right angles), constant building height, and flat terrain. The model does not include the effect of wave guiding through street canyons, which can lead to an underestimation of the received field strength. The COST 231-WI model has been accepted by ITU-R and is included in Report 567-4. The estimation of the path loss agrees well with measurements for base station antenna heights above rooftop level. The prediction error becomes large for hBase = hRoof compared to situations where hBase ! hRoof . Also, the performance of the model is poor for hBase " hRoof . The prediction error for micro-cells may be quite large. The reliability of the estimation decreases also if terrain is not flat or the land cover is inhomogeneous.
Page 12 of 21
13
V1 Created on 2/21/2007
terrain with light tree density or flat terrain with moderate-to-heavy tree density. Category B has an intermediate path loss. For all three categories, the median path loss at distance d > d 0 is given by: PL dB = 20log10(4 d 0 / ) + 10 log10(d/d 0 ) + s for d > d 0 Here, is the wavelength in meters, is the path-loss exponent with:
= a bhb + d/hb
hb is the height of the base station in meters (between 10 m and 80 m), d 0 = 100 m, and a, b, c are constants dependent on the terrain category. These parameters are listed in the table below.
Table A.2.1: Parameters of the Erceg Model Model Parameter Terrain Type A Terrain Type B Terrain Type C
a b c
4 0.0065 17.1
3.6 0.005 20
s represents the shadowing effect and follows a lognormal distribution with a typical standard deviation of 8.2 to 10.6 dB.
The above model is valid for frequencies close to 2 GHz and for receive antenna heights close to 2 m. For other frequencies and antenna heights (between 2 m and 10 m), the following correction terms are recommended [Molisch2005]: PLmodified = PL + PL f + PLh Here, PL, is the path loss given earlier, PL f is the frequency term, and PLh is the receive antenna height correction terms given as follows: PL f = 6log10( f/2000)
Page 13 of 21
14
V1 Created on 2/21/2007
C C B B A A
The terrain type A, B,C are same as those defined earlier for Erceg model. The multipath fading is modeled as a tapped delay line with 3 taps with non-uniform delays. The gain associated with each tap is characterized by a Rician Distribution and the maximum Doppler frequency. In a multipath environment, the received power r has a Rician distribution, whose pdf is given by:
pdf (r ) = r
e 2
[ fracr 2 + A2 2 2 ]I
rA 2
0r
Here, I 0 ( x) is the modified Bessel function of the first kind, zero order. A is zero if there is no LOS component and the pdf of the received power becomes:
r 2 2 r pdf (r ) = 2 e 2 0r
This is the Raleigh distribution. The ratio K = A2 /(2 2 ) in the Rician case represents the ratio of LOS component to NLOS component and is called the "K-Factor" or "Rician Factor." For NLOS case, K-factor is zero and the Rician distribution reduces to Raleigh Distribution. The general structure for the SUI channel model is as shown below in Figure A.2.5.1. This structure is for Multiple Input Multiple Output (MIMO) channels and includes other configurations like Single Input Single Output (SISO) and Single Input Multiple Output (SIMO) as subsets. The SUI channel structure is the same for the primary and interfering signals.
Page 14 of 21
15
V1 Created on 2/21/2007
Receiver
This part models correlation between input signals if multiple transmitting antennas are used.
A.2.5.2 Tapped Delay Line Matrix
This part models the multipath fading of the channel. The multipath fading is modeled as a tappeddelay line with 3 taps with non-uniform delays. The gain associated with each tap is characterized by a distribution (Rician with a K-factor > 0, or Raleigh with K-factor = 0) and the maximum Doppler frequency.
A.2.5.3 Output Mixing Matrix
This part models the correlation between output signals if multiple receiving antennas are used. Using the above general structure of the SUI Channel and assuming the following scenario, six SUI channels are constructed which are representative of the real channels.
A.2.5.4 Scenario for modified SUI channels Table A.2.5.1: Scenario for SUI Channel Models
Cell size BTS Antenna Height Receive Antenna Height BTS Antenna beam Width Receive Antenna Beam Width Vertical Polarization Only 90% cell coverage with 99.9% reliability at each location covered.
Page 15 of 21
16
V1 Created on 2/21/2007
In the following modles, the total channel gain is not normalized. Before using a SUI model, the specified normalization factors have to be added to each tap to arrive at 0dB total mean power. The specified Doppler is the maximum frequency parameter. The Gain Reduction Factor (GRF) is the total mean power reduction for a 30 antenna compared to an omni antenna. If 30 antennas are used the specified GRF should be added to the path loss. Note that this implies that all 3 taps are affected equally due to effects of local scattering. K-factors have linear values, not dB values. K-factors for the 90% and 75% cell coverage are shown in the tables, i.e., 90% and 75% of the cell locations have Kfactors greater or equal to the K-factor value specified, respectively. For the SUI channels 5 and 6, 50% K-factor values are also shown.
Table A.2.5.1: SUI 1 Channel Model
Tap 1 Tap 2 Tap 3 Units 0 0.4 0.9 s Delay 0 -15 -20 dB Power (omni ant.) 4 0 0 90% K-factor (omni) 20 0 0 75% K-factor (omni) 0 -21 -32 dB Power (30 ant.) 16 0 0 90% K-factor (30) 72 0 0 75% K-factor (30) 0.4 0.3 0.5 Hz Doppler Terrain Type: C Antenna Correlation: ENV = 0.7 Omni antenna: RMS = 0.111 s, Gain Reduction Factor: GRF = 0 dB Normalization Factor: Fomni = -0.1771 dB, overall K: K = 3.3 (90%); K = 10.4 (75%) F30 = -0.0371 dB 30 antenna: RMS = 0.042 s, overall K: K = 14.0 (90%); K = 44.2 (75%)
Page 16 of 21
17
V1 Created on 2/21/2007
0 0.4 Delay 0 -5 Power (omni ant.) 1 0 90% K-factor (omni) 7 0 75% K-factor (omni) 0 3 19 -11 Power (30 ant.) 0 90% K-factor (30) 0 75% K-factor (30) 0.4 0.3 Doppler Antenna Correlation: ENV = 0.4 Gain Reduction Factor: GRF = 3 dB Normalization Factor: Fomni = -1.5113 dB, F30 = -0.3573 dB
s dB dB
0.5 Hz Terrain Type: B Omni antenna: RMS = 0.264 s, overall K: K = 0.5 (90%); K = 1.6 (75%) 30 antenna: RMS = 0.123 s, overall K: K = 2.2 (90%); K = 7.0 (75%)
Page 17 of 21
18
V1 Created on 2/21/2007
The MATLAB code for implementing these six SUI models is available in [Erceg2001].
Tap
Channel A Relative delay (ns) 0 50 110 170 290 310 Average power (dB) 0 3.0 10.0 18.0 26.0 32.0
Channel B Relative delay (ns) 0 100 200 300 500 700 Average power (dB) 0 3.6 7.2 10.8 18.0 25.2
Doppler spectrum
1 2 3 4 5 6
Page 18 of 21
19
V1 Created on 2/21/2007
Table A.2.6.2: ITU Channel Model for Outdoor to Indoor and Pedestrian Test Environment
Tap
Channel A Relative delay (ns) 0 110 190 410 Average power (dB) 0 9.7 19.2 22.8
Channel B Relative delay (ns) 0 200 800 1 200 2 300 3 700 Average power (dB) 0 0.9 4.9 8.0 7.8 23.9
Doppler spectrum
1 2 3 4 5 6
Tap
Channel A Relative delay (ns) 0 310 710 1 090 1 730 2 510 Average power (dB) 0.0 1.0 9.0 10.0 15.0 20.0
Channel B Relative delay (ns) 0 300 8.900 12 900 17 100 20 000 Average power (dB) 2.5 0 12.8 10.0 25.2 16.0
Doppler spectrum
1 2 3 4 5 6
Each of these 6 cases is expected to occur some percentage of time. The following table gives the expected percentage of occurance and associated average RMS delay spread as specified in the recommendation:
Table A.2.6.4: Percentage Occurrence and Associated RMS Delay Spread for ITU Channel Models
Channel A Test environment Indoor office Outdoor to indoor and pedestrian r.m.s. (ns) 35 45
Channel B
P (%)
45 55
Page 19 of 21
20
V1 Created on 2/21/2007
370
40
4 000
55
WiMAX Forum recommends using just two of these six models. These are Ped-A and Veh-B models.
Page 20 of 21
21
V1 Created on 2/21/2007
References
[M.1225] ITU-R Recommendation M.1225, "Guidelines for evaluation of radio transmission technologies for IMT-2000," 1997 [Okumura1968] T. Okumura, E. Ohmore, and K. Fukuda, "Field strength and its variability in VHF and UHF land mobile service," Rev. Elec. Commun. Lab., pp. 825-73, September-October 1968. [Hata1980] M. Hata, "Empirical formula for propagation kiss ub kabd nibuke radui services," IEEE Trans. Veh. Tech., pp. 317-25, August 1980. [COST231] European Cooperative in the Field of Science and Technical Research EURO-COST 231, "Urban transmission loss models for mobile radio in the 900- and 1,800 MHz bands (Revision 2)," COST 231 TD(90)119 Rev. 2, The Hague, The Netherlands, September 1991, available at http://www.lx.it.pt/cost231/final_report.htm [Ikegami1984] F. Ikegami, S. Yoshida, T. Takeuchi, and M. Umehira, Propagation factors controlling mean field strength on urban streets, IEEE Trans. on Antennas Propagation, vol. 32, no. 8, pp. 822 829, Aug. 1984. [Walfisch1988] J. Walfisch, H. L. Bertoni, "A theoretical model of UHF propagation in urban environments," IEEE Trans. on Antennas and Propagation, vol. 36, no. 12, pp. 1788-1796, Dec. 1988 [Erceg1999] V. Erceg et. al, An empirically based path loss model for wireless channels in suburban environments, IEEE JSAC, vol. 17, no. 7, July 1999, pp. 1205-1211. [Hari2000] K.V. S. Hari, K.P. Sheikh, and C. Bushue, Interim channel models for G2 MMDS fixed wireless applications, IEEE 802.16.3c-00/49r2, November 15, 2000, available as www.ieee802.org/16/tg3/contrib/802163c-00_49r2.pdf [Erceg2001] V. Erceg, et al, "Channel Models for Fixed Wireless Applications," IEEE 802.16.3c01/29r4, July 2001, available as www.ieee802.org/16/tg3/contrib/802163c-01_29r4.pdf [Molisch2005] Andreas F. Molisch, "Wireless Communications," Wiley 2005, 622 pp. and associated web page http://www.wiley.com/legacy/wileychi/molisch/supp/appendices/Chapter_7_Appendices.pdf [Goldsmith2005] Andrea Goldsmith, "Wireless Communications," Cambridge University Press, 2005, 644 pp. [Andrews2007] Jeffrey G. Andrews, Arunabha Ghosh, and Rias Muhamed, "Fundamentals of WiMAX: Understanding Broadband Wireless Networking,", Prentice-Hall, February 2007, 496 pp. [Cichon] Dieter J. Cichon, Thomas Kurner, "Propagation Prediction Models," Chapter 4, Available at www.it.lut.fi/kurssit/04-05/010651000/Luennot/Chapter4.pdf
Page 21 of 21