Dr. David A. Clifton Group Leader Computational Health Informatics (CHI) Lab Lecturer in Engineering Science, Balliol College

Transcription

1 Dr. David A. Clifton Group Leader Computational Health Informatics (CHI) Lab Lecturer in Engineering Science, Balliol College

2 1. Introduction to Fourier analysis, the Fourier series 2. Sampling and Aliasing 3. Discrete Fourier methods, and Applications

3 1. Amplitude modulation 2. Sampling 3. Recovering a signal after sampling 4. Aliasing 5. Anti-aliasing

4 Amplitude modulation means to modulate the amplitude of some carrier wave according to some interesting signal: E.g., a high-frequency sinusoidal carrier wave injected through a patient s chest might be modulated by their respiration The envelope of the result then corresponds to the interesting signal (here, respiration)

5 Amplitude modulation is merely the multiplication of one signal g(t) (the carrier) by some other signal f(t) (the information wave) We have previously seen that if we multiply in the time domain, this corresponds to convolution in the frequency domain: f(t) x g(t) F(ω) * G(ω)...where we generally try to avoid convolution, primarily because it is terrifying:

6 And yet we can avoid all of the convolution difficulties. Recall that our information signal f(t) is being multiplied by a sinusoidal carrier wave g(t) f(t) x g(t) F(ω) * G(ω) We saw previously that the FT of a sinusoid was just a delta function:...so multiplying f(t) by a sinusoidal wave is the same as convolving with a delta function

7 Happily, convolving with a delta function is the one special case of convolution that we can easily handle: f(t) x g(t) F(ω) * G(ω) f(t) x cos(ω 0 t) F(ω) * π[δ(ω±ω 0 )]...because convolving with a delta function δ(ω±ω 0 ) just does the following to F(ω):

8 We could recover our original information from the modulated signal using demodulation, which just selects the relevant frequency information (using a filter)...and we can then recover our original time-signal using FT -1 (Hence the term modem: modulation-demodulation)

9 The concept of amplitude modulation is very useful, and employed whenever we wish to turn continuous signals (such as those we measure from patients) into discrete quantities that can be manipulated by computers: digitisation, a.k.a. sampling Whenever we sample, we multiply our original signal by a sampling function, p(t):

10 What is the sampling function, p(t)?...a.k.a. a train of pulses / Dirac deltas That s my pulse train

11 The FT of p(t) is: We know what happens to F(ω) when we convolve it with δ(ω±ω s )... This time, we re convolving F(ω) with δ(ω), δ(ω±ω s ), δ(ω±2ω s ), etc. I was the Lucasian Professor of Mathematics at a certain minor university, a position latterly held by Stephen Hawkin

12 Convolving F(ω) with δ(ω), δ(ω±ω s ), δ(ω±2ω s ), has the following effect:...where the frequency spectrum of our interesting signal f(t) is repeated at multiples (harmonics!) of the sampling frequency

13 Just as we previously recovered our original signal after amplitude modulation, we can pick out our original signal after sampling using a filter: (This is a reconstruction filter)...allowing us to recover our time-series by performing FT -1, as before

14 What could possibly go wrong with this idyllic situation? Sufficiently high ω s : Insufficiently high ω s : If the sampling frequency ω s is too low, the images of the frequency spectrum begin to overlap, termed aliasing

15 Use of our reconstruction filter would include components from higher- and lower- frequency images, hence we would not obtain the original signal: Aliasing can be avoided if ω s > 2 ω m (the Nyquist-Shannon theorem) Information Theory was mine, btw. Hej, I m actually Swedish [NIK-VIST]

16 The corresponding behaviour in the time domain is obvious is we consider a sinusoid of frequency f m : f s < 2 f m f s < 2 f m f s = f m Hence we have to sample at least twice every period in order to disambiguate, and so f s > 2 f m, or equivalently ω s > 2 ω m

17 Our previous signal was very well-behaved, in that its frequency spectrum F(ω) was limited to a certain band of frequencies this made it easy to separate from its images that appear after sampling. It was band-limited. However, most signals are not so well-behaved, and contain a wide spread of frequencies. We need to decide the highest frequency of interest in F(ω), and then suppress all higher frequencies using an anti-aliasing filter

18 We then sample the anti-aliased signal, which produces non-overlapped (unaliased) images in the frequency domain:

19 Finally, we can apply our reconstruction filter to pick out the original signal, after sampling: but this leads to potential problem #3...

20 Sadly, real filters are not ideal they are not brick wall filters, and have some roll-off slope If the filter is a low order, some components from the images will once again appear in the window further aliasing! Thus, the less ideal the filter, the greater the sampling frequency must be in order to adequately separate the images

21 Amplitude modulation results in: Sampling, which is a form of amplitude modulation, results in: Aliasing can be avoided by: A sampling rate > than the Nyquist-Shannon frequency Using an anti-aliasing filter to band-limit our signal Understanding that our reconstruction filter is non-ideal

22 Sampling allows us to represent our continuous time-domain signal using discrete pulses......but the corresponding frequency spectrum remains continuous! For analysis in a computer, we need to be able to store the frequency spectrum in a discrete form. We need the Discrete Fourier Transform.

23 Dr. David A. Clifton Group Leader Computational Health Informatics (CHI) Lab Lecturer in Engineering Science, Balliol College