Chapter 4
Mobile Radio Propagation:
 Small-Scale Fading and
       Multipath
    4.1 Small-Scale Multipath Propagation
•   The three most important effects
     – Rapid changes in signal strength over a small travel distance or time
       interval
     – Random frequency modulation due to varying Doppler shifts on different
       multipath signals
     – Time dispersion caused by multipath propagation delays
•   Factors influencing small-scale fading
     –   Multipath propagation: reflection objects and scatters
     –   Speed of the mobile: Doppler shifts
     –   Speed of surrounding objects
     –   Transmission bandwidth of the signal
          • The received signal will be distorted if the transmission bandwidth is greater
            than the bandwidth of the multipath channel.
          • Coherent bandwidth: bandwidth of the multipath channel.
•   Doppler Shift
     – A mobile moves at a constant velocity v, along a path segment having
       length d between points X and Y.
     – Path length difference
             ∆l = d cos θ = v∆t cos θ
     – Phase change
                 2π∆l       2πv∆t
         ∆φ =           =           cosθ
                  λ          λ
     – Doppler shift
                  1 ∆φ v
          fd =      ⋅ = cosθ
                 2π ∆t λ
          4.2 Impulse Response Model of a
                 Multipath Channel
•   A mobile radio channel may be modeled as a linear filter with a time
    varying impulse response
     – time variation is due to receiver motion in space
     – filtering is due to multipath
•   The channel impulse response can be expressed as h(d,t). Let x(t)
    represent the transmitted signal, then the received signal y(d,t) at
    position d can be expressed as
                                                     ∞
              y ( d , t ) = x (t ) ⊗ h( d , t ) =∫  −∞
                                                       x (τ )h(d , t − τ )dτ
•   For a causal system                     t
                             y (d , t ) = ∫−∞
                                               x (τ )h(d , t − τ )dτ
•   The position of the receiver can be expressed as
                                                 d = vt
•   We have
                                         t
                   y (vt , t ) =     ∫−∞
                                             x(τ )h(vt , t − τ )dτ
•   Since v is a constant, y (vt , t ) is just a function of t.
                                t
                       y (t ) =
                                −∞    ∫
                                   x (τ )h(vt , t − τ )dτ
•   In general, the channel impulse response can be expressed h(t ,τ )
     –    t : time variation due to motion
     –    τ : channel multipath delay for a fixed value of t.
•   With the channel impulse response h(t ,τ ) , we may have the output
                              t
               y (t ) =   ∫
                          −∞
                                  x (τ )h(t ,τ )dτ = x(t ) ⊗ h(t ,τ )
•   For bandlimited bandpass channel, then h(t ,τ ) may be equivalently
    described by a complex baseband impulse response hb (t ,τ )
     – The equivalent baseband output
         1         1       1                                           1
           r (t ) = c(t ) ⊗ hb (t ,τ )                 or      r (t ) = c(t ) ⊗ hb (t ,τ )
         2         2       2                                           2
                                           1
x(t ) = Re{c(t ) exp( jωc t )}     r (t ) = c(t ) ⊗ hb (t ,τ )
                                           2
y (t ) = Re{r (t ) exp( jωc t )}
• Discretize the multipath delay axis τ into equal time delay segments
  called excess delay bins.
• The baseband response of a multipath channel can be expressed as
                      N −1
         hb (t ,τ ) = ∑ ai (t ,τ ) exp( j 2πf cτ i (t ) + jφ (t ,τ ) )δ (τ − τ i (t ))
                      i =0
          ai (t ,τ ) : amplitude of the ith multipath component
             τ i (t ) : excess delay of ith multipath component
• Define θ i (t ,τ ) = 2πf cτ i (t ) + φ (t ,τ )
•   If the channel impulse response is assumed to be time invariant, the
    channel impulse response may be simplified as
                            N −1
                   hb (τ ) = ∑ ai exp( jθ i )δ (τ − τ i )
                             i =0
•   The impulse response may be measured by using a probing pulse p(t )
    which approximates a delta function.
                           p(t ) ≈ δ (t − τ )
    4.2.1 Relationship Between Bandwidth
             and Received Power
•   Consider a pulsed, transmitted signal of the form
                          x(t ) = Re{p (t ) exp( j 2πf c t )}
                                            p(t)
                            real response
                                                   Tbb                   t
                   imaginary response
                                                         TREP
•   The signal p(t) is a repetitive baseband pulse train with very narrow
    pulse width Tbb and repetition period TREP , with TREP >> τ max .
•   Now, let
                      p (t ) = 2 τ max / Tbb                    0 ≤ t ≤ Tbb
     •       The channel output r(t) closely approximates the impulse response and
             is given by
                                            N −1
                                            ∑ a exp(− jθ ) ⋅ p(t −τ )
                                        1
                             r (t ) =                 i              i             i
                                        2   i =0
                                        N −1
                                                                         τ max
                                   =    ∑
                                        i =0
                                               ai exp(− jθ i ) ⋅
                                                                          Tbb
                                                                                    T
                                                                                       2
                                                                                                 
                                                                              rect t − bb − τ i 
                                                                                                 
     •       Instantaneous multipath power delay profile
                     τ max
               1
                     ∫
     2
r (t0 ) =                r (t )r * (t )dt
             τ max   0
                                  N −1            N −1                                                                    
                                        ∑∑
                     τ max
               1             1 
         =
             τ max   ∫
                     0
                              Re 
                             4  j =0
                                                  i =0
                                                          a j (t0 )ai (t0 ) p (t − τ j ) p (t − τ i ) exp(− j (θ j − θ i ))dt
                                                                                                                            
•   If all the multipath components are resolved by the probe p(t), then
                           τ j − τ i > Tbb                   ∀ j≠i
•   Then we have
                                           1                   
                                                  N −1
                                                  ∑
                                   τ max
                          1
                                 ∫          a (t0 ) p (t − τ k ) dt
                 2
           r (t0 ) =                                     2      2
                     τ max 0               4                   
                                                         k
                                     k =0                        
                           N −1            τ max 
                                                   τ                             
                                 ∑                                           
                       1                                            T
                   =                             
                                                    ∫    rect t −      − τ k  dt
                                 2                   max             bb
                                ak (t0 )         
                     τ max k =0           0       Tbb              2        
                         N −1
                     =   ∑
                         k =0
                                ak2 (t0 )
•   The total receiving power is related to the sum of the powers in the
    individual multipath components.
•   Assuming that the received power from the multipath components
    forms a random process where each component has a random
    amplitude and phase at any time t, the average small-scale received
    power is
                                    N −1                2
                                                                N −1
             Ea ,θ [ PWB ] = Ea ,θ ∑
                                    i =0
                                           ai exp( jθ i )  =
                                                          
                                                                ∑
                                                                i =0
                                                                       ai2
•   Now, consider a CW signal which is transmitted into the exact same
    channel, and let the complex envelope be given by c(t)=2. Then the
    received signal can be expressed as
                                    N −1
                         r (t ) =   ∑ a exp( jθ (t,τ ))
                                    i =0
                                             i       i
•   The instantaneous power is given by
                                      N −1                      2
                              2
                         r (t ) =    ∑ a exp( jθ (t,τ ))
                                      i =0
                                                 i       i
•   In a local area, ai varies little, but θ i will vary greatly due to changes
    in propagation distance over space, resulting in large fluctuations of
    r(t).
•   The average received power over a local area is given by
                                    N −1                         
                                                                  2
            Ea ,θ [PCW ] = Ea ,θ 
                                    ∑      ai exp( jθ i (t ,τ )) 
                                                                  
                                    i =0
                                                                  
                            N −1             N −1     N
                        ≈   ∑
                            i =0
                                   ai2 + 2  ∑∑ r cos(θ −θ )
                                             i =0 i , j ≠i
                                                             ij   i   j
    where
                                        rij = Ea [ai a j ]
•   The received power for CW wave has large fluctuations than that for
    WB signal.
4.3 Small-Scale Multipath Measurement
•   Multipath channel measurement techniques
     – Direct pulse measurements
     – Spread spectrum sliding correlator measurements
     – Swept frequency measurements
           4.3.1 Direct RF Pulse System
•   Direct RF pulse system
     – This system transmits a repetitive pulse of width τ bb , and uses a
        receiver with a wideband filter with bandwidth BW = 2 / τ bb
     – Envelope detector to detect the amplitude response.
•   Minimum resolvable delay τ bb
•   No phase information can be measured.
    4.3.2 Spread Spectrum Sliding Correlator
                         Channel Sounding
•   System description
     – A carrier is spread over a large bandwidth by using a pseudo-noise
       sequence having chip duration Tc and a chip rate Rc .
     – Despread using a PN sequence identical to that used at the transmitter.
•   The power spectrum envelope of the transmitted
•   The probing signal is wideband.
•   Use a narrowband receiver preceded by a wideband mixer.
•   The transmitter chip clock is run at a slightly faster rate than the
    receiver chip clock – sliding correlator.
•   The time resolution of multipath components using a spread spectrum
    system with sliding correlation is
                                            1
                                 ∆τ = 2Tc =
                                            Rc
•   The time between maximum correlation can be calculated
                                              rl
                              ∆T = Tc rl =
                                              Rc
              Tc : chip period             r : sliding factor
              Rc : chip rate               l : sequence length
•   The sliding factor can be expressed as
                                       α
                                 r=
                                      α −β
    α : transmitter chip clock rate             β : receiver chip clock rate
•   The incoming signal is mixed with a PN sequence that is slower than
    the transmitter sequence. The signal is down converted to a low-
    frequency narrow band signal.
•   The observed time scale on the oscilloscope using a sliding correlator
    is related to the actual propagation time scale by
                                             Observed Time
             Actual Propagation Time =
                                                   r
                  actual channel response
                                     τ
                                         expansion by a factor of r
                                                                             t
                 display from oscilloscope
    4.3.3 Frequency Domain Channel Sounding
•   Dual relationship between time domain and frequency domain.
•   It is possible to measure the channel impulse response in the frequency
    domain.
•   Measure the frequency domain response and then converted to the time
    domain using inverse discrete Fourier transform (IDFT).
      4.4 Prameters of Mobile Multipath
                  Channels
•   Power delay profiles for different types of channels are different
             Outdoor                                 Indoor
       4.4.1 Time Dispersion Parameters
•   Time dispersion parameters
     – mean excess delay
     – RMS delay spread
     – excess delay spread
•   Mean excess delay
                               ∑a τ              2
                                                 k k     ∑ P(τ )τ               k       k
                         τ   =       k
                                                       =        k
                               ∑a        k
                                                  2
                                                  k      ∑ P(τ )    k
                                                                                    k
•   RMS delay spread
                     σ τ = τ 2 − (τ 2 )
    where
                           ∑a τ          2 2
                                         k k            ∑ P(τ )τ        k
                                                                                    2
                                                                                    k
                 τ   2
                         =   k
                                                      = k
                           ∑a    k
                                             2
                                             k          ∑ P(τ )
                                                            k
                                                                            k
•   Depends only on the relative amplitude of the multipath components.
•   Typical RMS delay spreads
     – Outdoor: on the order of microseconds
     – Indoor: on the order of nanoseconds
•   Maximum excess delay (X dB) is defined to be the time delay during
    which multipath energy falls to X dB below the maximum.
                      excess delay = τ X − τ 0
    τ X : maximum delay at which a multipath component is within X dB
    τ 0 : delay for the first arriving signal
•   Example of an indoor power delay profile; rms delay spread, mean
    excess delay, maximum excess delay (10dB), and the threshold level
    are shown
              4.4.2 Coherent Bandwidth
•   Coherent bandwidth, Bc , is a statistic measure of the range of
    frequencies over which the channel can be considered to be “flat”.
•   Two sinusoids with frequency separation greater than Bc are affected
    quite differently by the channel.
•   If the coherent bandwidth is defined as the bandwidth over which the
    frequency correlation function is above 0.9, then the coherent
    bandwidth is approximately
                                    1
                           Bc ≈
                                  50σ τ
•   If the frequency correlation function is above 0.5
                                   1
                           Bc ≈
                                  5σ τ
4.4.3 Doppler Spread and Coherent Time
•   Doppler spread and coherent time are parameters which discribe the
    time varying nature of the channel in a small-scale region.
•   When a pure sinusoidal tone of f c is transmitted, the received signal
    spectrum, called the Doppler spectrum, will have components in the
    range f c − f d and f c + f d , where f d is the Doppler shift.
                                  Channel
           fc                                          fc − fd   fc    fc + fd
•    f d is a function of the relative velocity of the mobile, and the angle
    between the direction of motion of the mobile and direction of arrival
    of the scattered waves
•   Coherent time TC is the time domain dual of Doppler spread.
•   Coherent time is used to characterize the time varying nature of the
    frequency dispersiveness of the channel in the time domain.
                                   1
                              TC ≈
                                   fm
    f m : maximum Doppler shift given by f m = v / λ
      v : speed of the mobile    λ : speed of the light
•   Two signals arriving with a time separation greater than TC are
    affected differently by the channel
•   A statistic measure of the time duration over which the channel
    impulse response is essentially invariant.
•   If the coherent time is defined as the time over which the time
    corrleation function is above 0.5, then
                                      9
                               TC ≈
                                    16πf m
        4.4 Types of Small-Scale Fading
•   Multipath delay spread leads to time dispersion and frequency selective
    fading.
•   Doppler spread leads to frequency dispersion and time selective fading.
•   Multipath delay spread and Doppler spread are independent of one
    another.
                      4.5.1 Flat Fading
•   If the channel has a constant gain and linear phase response over a
    bandwidth which is greater than the bandwidth of the transmitted
    signal, the received signal will undergo flat fading.
•   The received signal strength changes with time due to fluctuations in
    the gain fo the channel caused by multipath.
•   The received signal varies in gain but the spectrum of the transmission
    is preserved.
•   Flat fading channel is also called amplitude varying channel.
•   Also called narrow band channel: bandwidth of the applied signal is
    narrow as compared to the channel bandwidth.
•   Time varying statistics: Rayleigh flat fading.
•   A signal undergoes flat fading if
                            BS << BC
    and
                            TS >> σ τ
     TS : reciprocal bandwidth (symbol period)
    BS : bandwidth of the transmitted signal
    BC : coherent bandwidth
    σ τ : rms delay spread
       4.5.1 Frequency Selective Fading
•   If the channel possesses a constant-gain and linear phase response over
    a bandwidth that is smaller than the bandwidth of transmitted signal,
    then the channel creates frequency selective fading.
                           signal spectrum S ( f )
                        channel response
                                                          f
                                  BC
                     received signal spectrum
                                                          f
•   Frequency selective fading is due to time dispersion of the transmitted
    symbols within the channel.
     – Induces intersymbol interference
•   Frequency selective fading channels are much more difficult to model
    than flat fading channels.
•   Statistic impulse response model
     – 2-ray Rayleigh fading model
     – computer generated
     – measured impulse response
•   For frequency selective fading
                              BS > BC
    and
                              TS > σ τ
•   Frequency selective fading channel characteristic
     4.5.2 Fading Effects Due to Doppler
                   Spread
•   Fast Fading: The channel impulse response changes rapidly within the
    symbol duration.
     – The coherent time of the channel is smaller then the symbol period of the
       transmitted signal.
     – Cause frequency dispersion due to Doppler spreading.
•   A signal undergoes fast fading if
                             TS > TC
    and
                             BS < BD
•   Slow Fading: The channel impulse response changes at a rate much
    slower than the transmitted baseband signal s(t).
     – The Doppler spread of the channel is much less then the bandwidth of the
       baseband signal.
•   A signal undergoes slow fading if
                           TS << TC
    and
                           BS >> BD
    4.6 Rayleigh and Ricean Distributions
•   Rayleigh Fading Distribution
     – The sum of two quadrature Gaussian noise signals
•   Consider a carrier signal at frequency ω0 and with an amplitude a
                            s(t ) = a exp( jω0t )
•   The received signal is the sum of n waves
                    n
     sr (t ) = ∑ ai exp( jω0t + θi ) = r exp[ j (ω0t + θ )] = r exp( jθ ) exp( jω0t )
                   i =1
                                          n
    where
                           r exp( jθ ) = ∑ ai exp(θi )
                                         i =1
    define                       n                     n
                   r exp( jθ ) = ∑ ai cos(θi ) + j ∑ ai sin(θi ) = x + jy
                                i =1                  i =1
    We have
             n                                                    n
     x ≡ ∑ ai cos(θi ) = r cos(θ )              and          y ≡ ∑ ai sin(θi ) = r sin(θ )
            i =1                                                 i =1
•   It can be assumed that x and y are Gaussian random variables with
    mean equal to zero due to the following reasons
      – n is usually very large.
      – The individual amplitude ai are random.
      – The phases θ have a uniform distribution.
                      i
•   Because x and y are independent random variables, the joint distribution p(x,y)
    is
                                             1       x2 + y2 
            p ( x, y ) = p ( x ) p ( y ) =       exp −         
                                           2πσ 2         2σ 
                                                             2
                                                    
• The distribution p ( r , θ ) can be written as a function of p( x, y )
                           p( r , θ ) = J p( x, y )
                     ∂x / ∂r ∂x / ∂θ         cosθ     − r sin θ
               J =                       =                        =r
                     ∂y / ∂r ∂y / ∂θ         sin θ    r cosθ
•   We have
                                  r    r2 
                 p( r, θ ) =       exp − 2 
                             2πσ 2     2σ 
                                           
•   The Rayleigh distribution has a pdf given by
                                          r     r2 
                      2π                  2 exp − 2    r≥0
         p( r ) = ∫        p( r, θ )dθ = σ      2σ 
                  0                                  
                                               0       otherwise
                                         
•   pdf of Rayleigh distribution
                          r    r2 
                          exp − 2  0 ≤ r ≤ ∞
                p( r ) = σ 2   2σ 
                              0          r<0
                         
    σ: rms value of the received signal before envelop detection
    σ 2: time - average power of the received signal before envelop detection
•   Cumulative distribution function (CDF)
                                    R                       R2 
         P( R ) = Pr( r ≤ R ) = ∫       p( r )dr = 1 − exp − 2 
                                   0
                                                            2σ 
•   The mean value of the Rayleigh distribution is given by
                               ∞                   π
              rmean = E[r ] = ∫ rp( r )dr = σ          = 1.2533σ
                              0                    2
•   The variance of the Rayleigh distribution is given by
                                           ∞             σ 2π
             σ = E[r ] − E [r ] = ∫ r p( r )dr −
               2
               r
                        2     2                2
                                           0               2
                             π
                   = σ 2  2 −  = 0.4292σ 2
                             2
•   Ricean Fading 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.
                      Scattered waves              Direct wave
               sr (t ) = r exp[ j (ω0t + θ )] + A exp( jω0t )
                     = [( x + A) + jy ] exp( jω0t )
                  r 2 = ( x + A)2 + y 2
               x + A = r cosθ
                   y = r sin θ
•   By following similar steps described in Rayleigh distribution, we
    obtain
                    r    r 2 + A2   Ar 
                    exp −        I
                                 2  0  2 
                                                       A ≥ 0, r ≥ 0
          p( r ) = σ 2      2σ   σ 
                              0                            r<0
                   
     where
                     Ar  1        2π       Ar cosθ 
                 I0  2  =
                     σ  2π
                                ∫0
                                         exp
                                             σ
                                                  2   dθ
                                                      
     is the modified Bessel function of the first kind and zero-order.
•   The Ricean distribution is often described in terms of a parameter K
    which is defined as the ratio between the deterministic signal power
    and the variance of the multipath. It is given by K = A2 /(2σ 2 ) or in
    terms of dB                         2
                                           A
                    K (dB) = 10 log               dB
                                          2σ 2
•   The parameter K is known as the Ricean factor and completely
    specifies the Ricean distribution.
•   As A → 0 , we have K → −∞ dB. The dominant path decrease in
    amplitude, the Ricean distribution degenerates to a Rayleigh
    distribution.
4.7 Statistical Models for Multipath
          Fading Channels
       4.7.1 Clarke’s Models for Flat
                  Fading
•   Clark developed a model where the statistical characteristics of the
    electromagnetic fields of the received signal are deduced from
    scattering.
•   The model assumes a fixed transmitter with a vertically polarized
    antenna.
•   The received antenna is assumed to comprise of N azimuthal plane
    waves with arbitrary carrier phase., arbitrary angle of arrival, and each
    wave having equal average amplitude.
•   Equal amplitude assumption is based on the fact that in the absence of
    a direct line-of-sight path, the scattered components arriving at a
    receiver will experience similar attenuation over small-scale distance.
•   Doppler shift due to the motion of the receiver.
•   Assume no no excess delay due to multipath.
     – Flat fading assumption.
•   For the nth wave arriving at an angle α n to the x-axis, the Doppler
    shift is given by
                      ν
               fn =     cos α n
                      λ
•   The vertically polarized plane waves arriving at the mobile have E
    field components given by (assume a single tone is transmitted)
                                 N
                  E z (t ) = E0 ∑ Cn cos(2πf ct + θ n )
                                 n =1
    E0 : real amplitude of local average E - field (constant)
    Cn : real random variable representing the amplitude of nth arriving wave.
     f c : carrier frequency.
    θ n : random phase of the nth arriving wave.
•   The random arriving phase is given by
                                θ n = 2πf nt + φn
•   The amplitude of E-field is normalized such that
                                        N
                                        ∑ n =1
                                         C 2
                                        i =1
•   E z (t ) can be modeled as a Gaussian random process if N is sufficient
    large.
•   Since the Doppler shift is very small when compared to the carrier
    frequency, the three field components may be modeled as narrow band
    random process.
                E z (t ) = Tc (t ) cos(2πf ct ) + Ts (t ) sin(2πf ct )
    where                        N
                  Tc (t ) = E0 ∑ Cn cos(2πf nt + φn )
                                iN
                                 =1
                  Ts (t ) = E0 ∑ Cn sin(2πf nt + φn )
                                i =1
•   Tc (t ) and Ts (t ) are Gaussian random processes which are denoted as Tc
    and Ts , respectively.
• Tc (t ) and Ts (t ) are uncorrelated zero-mean Gaussian random variable
  with equal variance given by
                                         2
                      Tc2 = Tc2 = E z = E02 / 2
•   The envelope of the received E-field is given by
                 E z (t ) = Tc2 (t ) + Ts2 (t ) = r (t )
•   It can be shown that the random received signal envelope r has a
    Rayleight distribution given by
                         r    r2 
                         exp − 2  0 ≤ r ≤ ∞
              p ( r ) = σ 2   2σ 
                             0          r<0
                        
    where σ 2 = E02 / 2
•   Let p(α )dα denote the function of the total incoming power within dα
    of the angle α , and let A denote the average received power with
    respect to an isotropic antenna.
•   As N → ∞ , p(α )dα approached a continuous distribution.
•   If G (α ) is the azimuthal gain pattern of the mobile antenna as a
    function of the angle of arrival, the total received power can be
    expressed as
                     2π
               Pr = ∫ AG(α ) p(α )dα
                     0
•   The instantaneous frequency of the received signal arriving at an angle α
    is given by:
                           v
            f (α ) = f =       cos(α ) + f c = f m cos α + f c
                           λ
    where f m is the maximum Doppler shift.
•   If S(f) is the power spectrum of the received signal, the differential
    variation of received power with frequency is given by
            S ( f ) | df |= A[ p(α )G (α ) + p (−α )G (−α )] dα
•   Differentiation f = f m cos α + f c
           df
                = f m − sin α     ⇒        df = dα − sin α f m
          dα
                                      f − fc 
                                          −1
  On the other hand, we have α = cos         
                                         f
                                      m 
• This implies
                                                2
                                   f − fc 
                      sinα = 1 −         
                                   fm 
•   Finally, we have
                         A[ p (α )G (α ) + p ( −α )G (−α )]
              S( f ) =
                                                      2
                                         f − fc 
                               f m 1 −         
                                         fm 
     where S ( f ) = 0,        f − fc > fm
•   The spectrum is centered on the carrier frequency and is zero outside
    the limits f c ± f m .
•   Each of the arriving waves has its own carrier frequency (due to its
    direction of arrival) which is slightly offset from the center frequency.
•   Vertical λ / 4 antenna ( G (α ) = 1.5 ).
•   Uniform distribution p (α ) = 1 /( 2π ) over 0 to 2π .
•   The output spectrum
                                        1 .5
                    S( f ) =
                                                      2
                                         f − fc 
                               πf m 1 −        
                                         fm 
               S( f )
    4.7.2 Simulation of Clarke Fading
                  Model
•   Produce a simulated signal with spectral and temporal characteristics
    very close to measured data.
•   Two independent Gaussian low pass noise are used to produce the in-
    phase and quadrature fading branches.
•   Use a spectral filter to sharp the random signal in the frequency
    domain by using fast Fourier transform (FFT).
•   Time domain waveforms of Doppler fading can be obtained by using
    an inverse fast Fourier transform (IFFT).
•    Smith simulator using N carriers to generate fading signal
    1. Specify the number of frequency domain points N used to
        represent S ( f ) and the maximum Doppler frequency
        shift f m .
    2. Compute the frequency spacing between adjacent spectral lines
        as ∆f = 2 f m /( N − 1) . This defines the time duration of a fading
        waveform, T = 1 / ∆f .
    3. Generate complex Gaussian random variables for each of the
        N/2 positive frequency components of the noise source.
    4. Construct the negative frequency components of the noise
        source by conjugating positive frequency and assigning these at
        negative frequency values.
    5. Multiply the in-phase and quadrature noise sources by the fading
        spectrum S ( f ) .
    6. Perform an IFFT on the resulting frequency domain signal from
        the in-phase and quadrature arms, and compute the sum of the
        squares of each signal.
    7. Take the square root of the sum.
•   Frequency selection fading model
     4.7.3 Level Crossing and Fading
                 Statistics
•   The level crossing rate (LCR) is defined as the expected rate at which
    the Rayleigh fading envelope crosses a specified level in a positive-
    going direction.
•   Useful for designing error control codes and diversity.
•   Relate the time rate of change of the received signal to the signal level
    and velocity of the mobile.
•   The number of level crossing per second to the level R is given by
                 ∞
         N R = ∫ r&p ( R, r&)dr& = 2π f m ρe   −ρ 2
                                                                     (A)
                0
    r& : time derivation of r(t) (slope)
    p ( R, r&) : joint density function of r and r& at r = R.
    f m : maximum Doppler frequency
    ρ = R / Rrms : value of the specified level R, normalized to the rms
    amplitude of the fading envelope.
•   Average fading duration is defined as the average period of time for
    which the received signal is below a specified level R.
•   For a Rayleigh Fading signal, this is given by
                                1
                        τ =        Pr[ r < R ]                      (B)
                                NR
    with                                      1
                              Pr[ r < R ] =     ∑
                                              T i
                                                  τi
     where τ i is the duration of the fade and T is the observation interval.
•   For Rayleigh distribution
                          R
           Pr[ r < R ] = ∫ p (r )dr = 1 − exp( − ρ 2 )               (C)
                         0
•   Average fading duration, (using (A), (B), (C))
                                      2
                                   eρ −1
                              τ =
                                  ρf m 2π
•   The average duration of a signal fading helps determine the most likely
    number of signaling bits that nay be lost during a fade.
•    Average fade duration primarily depends upon the speed of the mobile,
    and decreases as the maximum Doppler frequency f m becomes large.