Wireless Channel: Dr. Manjunatha. P
Wireless Channel: Dr. Manjunatha. P
Dr. Manjunatha. P
[email protected]
Professor
Dept. of ECE
March 5, 2015
Wireless Channel Syllabus
Wireless Channel
Physical modeling for wireless channels
1 Free space, fixed transmit and receive antennas
2 Free space, moving antenna
3 Reflecting wall, fixed antenna
4 Reflecting wall, moving antenna
5 Reflection from a ground plane
6 Power decay with distance and shadowing
7 Moving antenna, multiple reflectors
Input /output model of the wireless channel
1 The wireless channel as a linear time-varying system
2 Baseband equivalent model
3 A discrete-time baseband model
Time and frequency coherence
1 Doppler spread and coherence time
2 Delay spread and coherence bandwidth
Statistical channel models
1 Modeling philosophy
2 Rayleigh and Rician fading
3 Tap gain auto-correlation function
Introduction
The transmission path between transmitter and receiver is varying due to obstruction by
buildings, mountains and foliage (trees).
Wired channels are stationary and predictable, but radio channels are extremely random
and difficult to analzye.
Modeling the radio channels has been one of the most difficult task of mobile radio
system design, and is typically done in a statical method based on measurements [2].
The electromagnetic wave propagation is attributed to reflection, diffraction and
scattering.
The high rise buildings cause severe diffraction loss.
The received signal strength decreases as the distance between the transmitter and
receiver increase.
Due to multiple reflectors from objects, the electromagnetic wave travel along different
paths of varying lengths and causes multipath fading.
The propagation models that predict the mean signal strength as a function of distance
are called large-scale fading.
The propagation models that characterize the rapid fluctuations of the received signal
strength as a function of short distance or short time durations are called sall-scale fading.
Diffraction of signals: Occurs when the path between transmitter and receiver is
obstructed by a sharp edges of an objects.
Reflection of signals: Occurs when EM wave impinging upon a object which has very large
in dimension when compared to the wavelength.
Scattering of signals: Occurs when the medium through which the wave travels consists of
objects which are small in dimension in compared to wavelength and very large number of
obstacles per unit volume is large.
Physical modeling for wireless channels Free space, fixed transmit and receive antennas
Consider a fixed antenna radiating its electromagnetic field into the free space.
The electric and magnetic field at any given location are perpendicular to each other and
are to the direction of propagating antenna.
Consider an electromagnetic sinusoid wave with frequency f i.e, cos2πft is transmitted.
The electric far field received at a point u with a distance r is
1
The electric field is proportional to the distance r and given as E ∝ r
and the power per
square meter in the free space wave decreases as r −2 .
If we considered electric field as concentric spheres of increasing radius r around the
antenna, the total power radiated through the sphere remains constant, but the surface
area increases as r 2 .
For a fixed receive antenna at the location u = (r , θ, ψ), the received electric field is
α(θ, ψ, f )e −j2πfr /c
H(f ) =
r
Er (f , t, u) = [H(f )e j2πft ]
The sinusoid at frequency f has been converted to a sinusoid of frequency −fv /c there
has been a Doppler shift of −fv /c due to the motion of the observation point.
The electric field at a point u(t), is
Transmit Reflecting
antenna Receive Wall
antenna
When the phase difference is an integer multiple of 2π, the two waves add constructively,
and the received signal is strong.
When the phase difference is an odd integer multiple of π, the two waves add
destructively, and the received signal is weak.
As a function of r, this translates into a spatial pattern of constructive and destructive
interference of the waves.
The distance from a peak to a valley is called the coherence distance:
λ
Δxc =
4
where λ = c/f is the wavelength of the signal. At distances much smaller than Δxc , the
received signal at a particular time does not change appreciably.
The distance in space over which a fading channel appears to be unchanged.
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 11 / 49
Physical modeling for wireless channels Reflecting wall, fixed antenna
The constructive and destructive interference pattern also depends on the frequency f: for
a fixed r, if f changes by
1 2d − r r −1
−
2 c c
2d − r r
Td = −
c c
The constructive and destructive interference pattern does not change appreciably if the
frequency changes by an amount much smaller than 1/Td (coherence bandwidth).
Coherence bandwidth: The range of frequencies over which the channel can be
considered flat. GSM is of 200 KHz bandwidth, if coherence bandwidth is less than 200
KHz then interference will occur.
Transmit Reflecting
antenna Receive Wall
antenna
r(t)
Speed v
d
The first term, the direct wave, is a sinusoid at frequency f (1 − v /c), experiencing a
Doppler shift D1 := −fv /c. The second is a sinusoid at frequency f (1 + v /c),
experiencing a Doppler shift D2 := +fv /c. The Dopper shift is positive, if the mobile is
moving toward the direction of arrival of the wave and is negative, if the mobile is moving
away from the direction of arrival of the wave. The Doppler spread is defined as
Ds = D2 − D1
Example, Consider a mobile is moving at 60 km/h and the carrier frequency of the signal
is f = 900 MHz, then the Doppler spread is
D1 = fv /c = [(900x106 )x16.66]/(3x108 ) = 49.8 50Hz
Ds = D2 − D1 = 50 + 50 = 100Hz
Consider a mobile is moving at 60 miles/h and the carrier frequency is f = 1850 MHz,
then the Doppler spread is
D1 = fv /c = [(1850x106 )x26.82]/(3x108 ) = 49.8 16Hz
Ds = D2 − D1 = 16 + 16 = 32Hz
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 14 / 49
When the attenuations are roughly same for both paths, then approximate the
denominator of the second term by r = r0 + vt.
Then, combining the two sinusoids,
This is the product of two sinusoids, one at the input frequency f, which is typically of the
order of GHz, and the other one at fv /c = Ds /2, which might be of the order of 50 Hz.
Thus, the response to a sinusoid at f is another sinusoid at f with a time-varying envelope,
with peaks going to zeros around every 5 ms.
The envelope is at its widest when the mobile is at a peak of the interference pattern and
at its narrowest when the mobile is at a valley.
Thus, the Doppler spread determines the rate of traversal across the interference pattern
and is inversely proportional to the coherence time of the channel.
Coherence time: The time duration in which channel impulse response is invariant.
1
Tc =
fm
Where fm is maximum doppler shift. If the reciprocal of the baseband of the signal is
greater than the coherence time then the channel will change during the transmission of
the baseband message thus causes the distortion at the receiver.
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 15 / 49
Physical modeling for wireless channels Reflection from a ground plane
r1
Direct path Receive
hs
r2 antenna
Reflected
path
Ground plane hr
E0 d0 cos[ωc (t − d /c)] E0 d0 cos[ωc (t − d /c)]
Er (f , t) = + (−1)
d d
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 16 / 49
Physical modeling for wireless channels Power decay with distance and shadowing
Shadowing
The density of obstacles between the transmit and receive antennas depends very much
on the physical environment.
For example, outdoor plains have very little by way of obstacles while indoor environments
pose many obstacles.
This randomness in the environment is captured by modeling the density of obstacles and
their absorption behavior as random numbers; the overall phenomenon is called shadowing.
The effect of shadow fading differs from multipath fading in an important way.
The duration of a shadow fade lasts for multiple seconds or minutes, and hence occurs at
a much slower time-scale compared to multipath fading.
Physical modeling for wireless channels Power decay with distance and shadowing
Dealing with multiple reflectors, using the technique of ray tracing, is in principle simply a
matter of modeling the received waveform as the sum of the responses from the different
paths rather than just two paths.
The reflected field model is valid only at distances from the wall that are small relative to
the dimensions of the wall.
At very large distances, the total power reflected from the wall is proportional to both
d −2 and to the area of the cross section of the wall.
The power reaching the receiver is proportional to [d − r (t)]−2 .
Thus, the power attenuation from transmitter to receiver (for the large distance case) is
proportional to [d − r (t)]−2 rather than to [2d − r (t)]−2 .
Scattering
The type of reflection which occur in the atmosphere or in reflections from very rough
objects is known as scattering.
Here there are a very large number of individual paths, and the received waveform is better
modeled as an integral over paths with infinitesimally small differences in their lengths.
y (t) = ai (f , t)φ[t − τi (f , t)]
i
Input/output model of the wireless channel The wireless channel as a linear time-varying system
When the transmitter, receiver and the environment are all stationary, ai (t) and τi (t) do
not depend on time t, then linear time-invariant channel with an impulse response
h(τ ) = ai (t)δ[t − τi ]
i
For the time-varying impulse response h(τ, t) the time varying frequency response
∞
H(f , t) = h(τ, t)e −j2πf τ dτ = ai (t)e −j2πf τi (t)
−∞ i
1
[ sb (t )]
[ sb (t )]
W W
2 2
s (t )
1
[ sb (t )]
[ sb (t )]
W W
2 2
2 sin 2S f ct 2 sin 2S f ct
Let xb (t) and yb (t) be the complex baseband equivalents of the transmitted signal x(t)
and the received signal y (t), respectively.
Figure 10 shows the system diagram from xb (t) to yb (t).
This implementation of a passband communication system is known as quadrature
amplitude modulation (QAM).
The signal [xb (t)] is sometimes called the in-phase component I and [xb (t)] the
quadrature component Q (rotated by π/2).
2 cos2S fct 2 cos2S fct
1
[ xb (t)]
[ yb (t )]
W W
2 2
x(t ) y (t )
h(W , t )
1
[ yb (t )]
[ xb (t )]
W W
2 2
√ √
x(t) = 2[xb (t)e j2πfc t ] and y (t) = 2[yb (t)e j2πfc t ]
[yb (t)e j2πfc t ] = ai (t)[xb (t − τi (t))e j2πfc (t−τi (t)) ]
i
= ai (t)xb (t − τi (t))e −j2πfc τi (t) e j2πfc t
i
Similarly,
[yb (t)e j2πfc t ] = ai (t)xb (t − τi (t))e −j2πfc τi (t) e j2πfc t
i
where
aib (t) = ai (t)e −j2πfc τi (t)
The baseband output is the sum, over each path, of the delayed replicas of the baseband
input.
The magnitude of the ith such term is the magnitude of the response on the given path.
The phase is changed by π/2 when the delay on the path changes by 1/(4fc ), or
equivalently, when the path length changes by a quarter wavelength, i.e., by c/(4fc ).
If the path length is changing at velocity v, the time required for such a phase change is
c/(4fc v ).
The Doppler shift D at frequency f is fv /c, and noting that f ≈ fc for narrowband
communication, the time required for a π/2 phase change is 1/(4D).
For the single reflecting wall example, this is about 5 ms (assuming fc = 900MHz and v =
60km/h). The phases of both paths are rotating at this rate but in opposite directions.
S 3S
S 2S
2 2
T
1
T
fc
S
Phase shifted by
2
1
4 fc
sin(πt)
sinc(t) =
πt
The sampled outputs at multiples of 1/W, y [m] = yb (m/W ), are then given by
y [m] = x[n] aib (m/W )sinc[m − n − τi (m/W )W ]
n i
Let = m − n.
y [m] = x[m − ] aib (m/W )sinc[ − τi (m/W )W ]
i
= h [m]x[m − ]
n
When the gains aib (t) and the delays τi (t) of the paths
are time-invariant, then the channel is linear
time-invariant.
h [m] = aib sinc[ − τi W ]
i Figure 11: Plot of h[m]
The th tap can be interpreted as the sample (/W )th of the low-pass filtered baseband channel
response hb (τ ) convolved with sinc(W τ ).
[ xb (t )] 1 [ yb (t )] [ y[m]]
[ x[m]]
sin c(Wt n)
W W
2 2
x(t ) y (t )
h(W , t )
[ xb (t )] 1 [ yb (t )] [ y[ m]]
[ x[m]]
sin c(Wt n)
W W
2 2
y (t) = ai (t)x(t − τi (t)) + w (t)
i
y [m] = hl [m]x[m − l] + w [m]
l
[ xb (t )] 1 [ yb (t )] [ y[m]]
[ x[m]] w(t )
sin c(Wt n)
W W
y (t ) 2 2
x(t )
h(W , t )
[ xb (t )] 1 [ yb (t )] [ y[m]]
[ x[m]]
sin c(Wt n)
W W
2 2
where w[m] is the low-pass filtered noise at the sampling instant m/W.
The white noise w(t) is down-converted, filtered at the baseband and ideally sampled.
∞
(w [m]) = w (t)ψm,1 (t)dt
−∞
∞
(w [m]) = w (t)ψm,2 (t)dt
−∞
where
√
ψm,1 (t) = 2W cos(2πfc t)sinc(Wt − m)
√
ψm,2 (t) = − 2W sin(2πfc t)sinc(Wt − m)
From this expression the significant changes in ai occur over periods of seconds.
Significant changes in the phase of the ith path occur at intervals of 1/(4Di ), where
Di = fc τi (t) is the Doppler shift for that path.
Due to the different Doppler shifts from different paths will have significant changes in
the magnitude of h [m].
The Doppler spread Ds is defined as the largest difference between the Doppler shifts:
Ds = max fc |τi (t) − τj (t)|
Coherence Time
The coherence time Tc of a wireless channel is defined as the interval over which h [m] changes
significantly as a function of m.
1
Tc =
4Ds
If the symbol period of a baseband signal is greater than the coherence time, then the
signal will distort, since channel will change during the transmission of the signal.
Coherence time definition implies that two signals arriving with a time separation greater
than Tc are affected differently by the channel.
If the coherence time is defined as the time over which the time correlation function is
above 0.5, then it is approximated as
9 0.179
Tc = =
16fm fm
v
where fm = maximum Doppler shift and is given by fm = fd, max = λ
= v cf
The channels are categorized as fast fading and slow fading.
The channel is fast fading if the coherence time Tc is much shorter than the delay
requirement of the application, and slow fading if Tc is longer.
A channel is fast or slow fading depends not only on the environment but also on the
application.
For example in voice communication, typically has a short delay requirement of less than
100 ms, while some types of data applications can have a laxer delay requirement
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 35 / 49
Time and frequency coherence Delay spread and coherence bandwidth
The multipath delay spread, Td , is defined as the difference in propagation time between
the longest and shortest path, counting only the paths with significant energy.
Td = max|τi (t) − τj (t)|
Modeling philosophy
All analytical work is done with simplified models, like white Gaussian noise is often
assumed in communication models but the model is valid only for small frequency bands.
Doppler spread, multipath spread, etc. are defined for wireless channel with probabilistic
models, but these channels are very different from each other and cannot be characterized
by probabilistic models.
Consider a continuous time multipath fading channel defined as
y (t) = ai (t)x(t − τi (t)) + w (t)
i
where h [m]
h [m] = ai (m/W )e −j2πfc τ (m/W ) sinc[ − τ (m/W )W ]
i
−5
−10
−15
−20
−25
−30
−35
−40
0 200 400 600 800 1000 1200 1400 1600 1800 2000
The phase of the ith path is 2πfcτi modulo 2π. It is assumed that the phase for each
path is uniformly distributed between 0 and 2π and that the phases of different paths are
independent.
fc τi = di /λ
where di is the distance traveled by the ith path and λ is the carrier wavelength.
Since the reflectors and scatterers are far away relative to the carrier wavelength, i.e.,
di λ,
It is assumed that the phase for each path is uniformly distributed between 0 and 2π and
that the phases of different paths are independent.
It follows that (h [m]) is the sum of many small independent real random variables, and
so by the Central Limit Theorem, it can reasonably be modeled as a zero-mean Gaussian
random variable.
Similarly, because of the uniform phase,(h [m]e jφ ) is Gaussian with the same variance
for any fixed φ.
The magnitude h [m] of the th tap is a Rayleigh random variable with probability density
function (pdf) is given by
x x2
p(x) = 2 exp − 2 x ≥0
σ 2σ
Ricean Distribution
When there is a dominant stationary (non-fading) signal component present, such as a
line-of-sight propagation path, the small-scale fading envelope distribution is Ricean
The pdf of the Ricean distribution is given by
x x 2 + A2 Ax
p(x) = 2 exp − I0 A ≥ 0, x ≥ 0
σ 2σ 2 σ2
k 1
h [m] = σl e jθ + CN(0, σl 2 )
k +1 k +1
The first term corresponding to specular path and second term corresponding to the
aggregation of the large number of reflected and scattered paths.
The parameter k is the ratio of the energy in the specular path to the energy in the
scattered paths.
Dr. Manjunatha. P (JNNCE) Wireless Channel[1, 2] March 5, 2015 43 / 49
Statistical channel models Tap gain auto-correlation function
Modeling each h [m] as a complex random variable provides the statistical description.
A statistical quantity that models the rate of channel variation is known as the tap gain
auto-correlation function, R [n]. It is defined as
For each tap , R [n] gives the auto-correlation function of the sequence of random
variables assuming that this is not a function of time m.
h [m] is independent of h [m ] for all = and all m, m . This final assumption is
intuitively plausible since paths in different ranges of delay contribute to h [m] for
different values of .
Clarkes model
This is a popular statistical model for flat fading.
The transmitter is fixed, the mobile receiver is moving at speed v, and the transmitted
signal is scattered by stationary objects around the mobile.
There are K paths, the ith path arriving at an angle θi : 2πi/K , i = 0 . . . K − 1, with
respect to the direction of motion.
The scattered path arriving at the mobile at the angle θ has a delay of τθ (t) and a time
invariant gain aθ (t), and the input/output relationship is given by
K −1
y (t) = aθi x(t − τθi (t))
i=0
The received power distribution p(θ) and the antenna gain pattern α(θ) are functions of
the angle θ
It is assumed as uniform
√ power distribution and isotropic antenna gain pattern, i.e., the
amplitudes aθ = a/ K for all angles θ.
The details of the description is as shown in Figure 17 known as the one ring model.
√
The amplitude of each path scaled by K so that the total received energy along all
paths is a2 .
When the communication bandwidth W is much smaller than the reciprocal of the delay
spread, then the complex baseband channel can be represented by a single tap at each
time:
y [m] = h0 [m]x[m] + w [m]
The phase of the signal arriving at time 0 from an angle θ is 2πfc τθ (0) mod 2π, where fc
is the carrier frequency.
Making the assumption that this phase is uniformly distributed in [0, 2π] and
independently distributed across all angles θ, the tap gain process hl [m] is a sum of many
small independent contributions, one from each angle.
By the Central Limit Theorem, it is modeled the process as Gaussian with an
autocorrelation function R0 [n] as shown in Figue 18given by:
R0 [n] = 2a2 πJ0 (nπDs /W )
where J0 (.) is the zeroth-order Bessel function of the first kind
π
Rx
J0 (x) = e jxcosθ dθ
0
Figure 18: Auto-correlation function
Figure 19: Doppler spectrum
Sinwave
1
Amplitude
0.5
−0.5
−1
0 1 2 3 4 5 6
Time, [s]
Autocorrelation function of the sinewave
1000
Autocorrelation
500
−500
0 200 400 600 800 1000 1200
Lags
An aircraft is heading towards a control tower with 500 kmph, at an elevation of 20◦ .
Communication between aircraft and control tower occurs at 900 MHz. Find out the expected
Doppler shift.
Solution
v = 500 kmph
The horizontal component of the velocity is
v = vcosθ = 500 × cos20. = 130m/s
900×106
Hence, it can be written that λ = 3×108
= 13m
130
fd = 1/3
= 390Hz
If the plane banks suddenly and heads for other direction, the Doppler shift change will be
390Hz to -390Hz.
References