white noise
Turbulence velocity
Linear Filter (‘coloured noise')
(Dryden)
Figure 4.4: Modelling atmospheric turbulence as filtered white noise
50 Chapter 4. External atmospheric disturbances
To apply these relations, the spatial spectral density functions of the turbulence velocities
must be transformed to functions of the circular frequency w. This can be
done, because Taylor’s hypothesis implies that the circular frequency w [rad s−1] encountered
at the aircraft’s center of gravity is related directly to the spatial frequency
W [rad m−1]:
w = WV (4.17)
The resulting transformation is given by:
S(w) = 1
V
S (W) (4.18)
The extra term 1/V in the spectral density function is due to the fact that we now
integrate over the time instead of over the distance in the Fourier transform.
The Dryden spectra were developed to approximate the von Kármán turbulence
spectra by means of rational functions. This makes it possible to apply relation (4.18)
for the generation of turbulence velocity components from white noise inputs. From
the definitions of the Dryden spectra in equations (4.11) to (4.13) and relation (4.18)
the following expressions are found:
|Hugw1 (w)|2 = 2su
2 Lu
V
1
1 +
Lu
w
V
2 (4.19)
|Hvgw2 (w)|2 = sv
2 Lv
V
1 + 3
Lv
w
V
2
1 +
Lv
w
V
2
2 (4.20)
|Hwgw3 (w)|2 = sw
2 Lw
V
1 + 3
Lw
w
V
2
1 +
Lw
w
V
2
2 (4.21)
Solving equations (4.19) to (4.21) yields the following candidate functions for the
frequency responses of the forming filters:
Hugw1 (w) = su
r
2Lu
V
1
1 ± Lu
V jw
(4.22)
Hvgw2 (w) = sv
r
Lv
V
1 ±
p
3 Lv
V jw
1 ± Lv
V jw
2 (4.23)
Hwgw3 (w) = sw
r
Lw
V
1 ±
p
3 Lw
V jw
1 ± Lw
V jw
2 (4.24)
In these equations w1, w2, and w3 are independent white noise signals. Choosing the
minus sign in the denominators would lead to unstable filters and hence should be
rejected for physical reasons. And choosing the minus sign in the numerators leads
to non-minimum phase systems, so we shall use positive signs in both the numerator
and denominator [27].
The Laplace transforms H(s) are obtained by replacing the imaginary variable iw
by the more general complex variable s, and the resulting transfer functions can subsequently
be converted to state-space systems using the technique from section 6.4.2.
4.2. Stochastic disturbances 51
中国航空网 www.aero.cn
航空翻译 www.aviation.cn
本文链接地址:FDC 1.4 – A SIMULINK Toolbox for Flight Dynamics and Contro(29)