Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic signal g(t) with period T is constructed by repeating the signal f (t) shown in Figure. (a) Find the coefficients C n of the complex Fourier series describing g(t). (b) Use the coefficients to deduce the Fourier transform G(ω) of g(t) directly from the standard formula for the FT of a complex Fourier series. (c) Now rewrite the periodic function g(t) as a convolution between f (t) and a train of δ-functions. (d) Determine G(ω) using the time-convolution property of Fourier transforms, and confirm that it agrees with that derived in part (b). f(t) T/ 2 a/ 2 0 a/ 2 T/ 2 t Figure : 2. The signal whose double-sided spectrum is shown in Figure 2 is sampled using a train of δ-functions before being passed through an ideal reconstruction filter. (a) What is the Nyquist limit? For sampling frequencies above and below the Nyquist limit, sketch and explain the appearance of the frequency spectra (i) after sampling and (ii) after reconstruction. (b) When is an anti-aliasing filter useful? Where is it placed in the processing chain? (c) Music is recorded on CDs at a sampling rate of 44. khz per channel, and players usually have reconstruction filters with a nearly flat frequency response over 0-20 khz, and rejection of greater than 70 db above 22 khz. Suggest why these various frequencies are chosen. ω m p(t) Figure 2: Reconstruction Filter
2 Time-Frequency Analysis II (HT 20) 2AH 3. A signal x(t) = cos(ω 0 t) is sampled at ω s with p(t) = δ(t kt s ) and then reconstructed from the sampled signal using an ideal low-pass filter with cut-off frequency of ω s /2. (a) Show that for a signal with spectrum X(ω) the sampled signal has spectrum Y (ω) = ω s + k= X(ω kω s ). and for this particular signal is Y (ω) = ω s 2 + k= [δ(ω kω s ω 0 ) + δ(ω kω s + ω 0 )]. (b) What is the reconstructed signal x r (t) when (i) ω 0 = ω s /6 (ii) ω 0 = 2ω s /3 (iii) ω 0 = ω s Illustrate your answers with a sketch showing the frequency spectra before and after filtering. (c) In each of these three cases, what would be the reconstructed signal were an anti-aliasing LPF with cut-off frequency ω s /2 used before sampling? 4. (a) Explain what is meant by a finite energy signal and contrast it with a finite power signal. (b) For finite energy signals Parseval s theorem states that the total energy in the time domain is equal to the total energy in the frequency domain. Prove the theorem. (c) By using a Fourier transform pair from HLT, verify the theorem by calculating, in both the time and frequency domains, the energy dissipated in a 47Ω resistor subjected to a time-varying voltage v(t) = 2u(t) exp( 2t) where u(t) is the unit step function. (See Useful Integral # at the end.)
Time-Frequency Analysis II (HT 20) 2AH 3 5. (a) By finding the energy in the time domain, confirm that f (t) = u(t)ae at, where a is a positive constant, is a finite energy signal, (b) For any finite energy signal f (t) derive the Wiener-Khinchin theorem that the Fourier transform of the autocorrelation R f f (τ) is equal to the energy spectral density E(ω). (c) By integration in the frequency domain, show that the total energy in f (t) is E = A 2 /2a and identical with that found in part (a). (d) Let f (t) be the input to a filter with impulse response h(t) = Be bt, b > 0, and let y(t) be the output of the filter. Compute the total energy of y(t). (See useful integral #2.) 6. (a) Explain the terms time-average, ensemble-average, stationary random function and ergodic random function. (b) A zero-mean stationary random process x(t) has autocorrelation R xx (τ) = Ae τ. Determine its average power and its power spectral density S xx (ω). (c) Determine the integral S xx(ω)dω. Is the value what you expected? (d) Let y(t) be the output response of a dynamic system transfer function H(s) = (/ + T s) when the system is driven by x(t). Determine the expected value and the variance of y(t). Show that the variance is also the average power of y(t).
4 Time-Frequency Analysis II (HT 20) 2AH 7. (a) Ideal white noise has a power spectral density S xx (ω) =. Determine its autocorrelation R xx (τ). (b) Ideal white noise x is band-limited to ω < ω 0 by an ideal low pass filter H(ω) = for ω < ω 0, and H(ω) = 0 elsewhere. Determine the autocorrelation of the filter s output. (c) Ideal white noise is now limited by the st order low pass filter shown in Figure 3 which has a half power (break) point at angular frequency ω 0. Derive expressions for the power spectral density S yy and autocorrelation R yy of the output. (d) Sketch the autocorrelations of parts (a), (b) and (c), and comment on your results. x(t) R C y(t) Figure 3: 8. Each pulse in a continuous binary train has a fixed duration, T, but takes the value, 0, or - with probabilities /4, /2 and /4 respectively, independently of the previous pulse. (An example is sketched in Figure 4.) (a) If x i is the random variable representing the height of pulse i, calculate E[x i ], E[x 2 i ], and E[x i x j ] for i j. (b) Show that the auto-correlation of the pulse train is R xx (τ) = (/2) ( τ /T ) for τ T, and zero elsewhere. (c) Hence find the power spectral density S xx (ω) of the pulse train. (d) Use the autocorrelation to find the mean power, and confirm that this agrees with the result from part (a). 0 t Figure 4:
Time-Frequency Analysis II (HT 20) 2AH 5 9. Be prepared to discuss the following with your tutor: (a) Explain how a sampled signal from can be reconstructed in the time domain, rather than using an analogue reconstruction filter. (b) What is over-sampling? Describe its effect on the frequency spectrum of a sampled signal. What benefit might it have for reconstruction in the time domain? (c) Describe a practical use for under-sampling.. ( ω 2 + a 2 ) dω = 2a 2. Two tan substitutions... ( ω 2 + a 2 ) ( Useful Integrals... a tan substitution. ) ω 2 + b 2 dω = = = b 2 a 2 b 2 a 2 2ab(a + b) [ [ 2a 2b ω 2 + a 2 ω 2 + b 2 ] ] dω Some Answers and Hints See the course page at www.robots.ox.ac.uk/ dwm/courses/2tf