"Objectives" of Lecture #DSP
The need for DSP Aliasing & windowing Introduction to FFT.
Why Digital Signal Processing?
Allows fast & complex calculations, usually without loss of information.
E.G.
MRIs Communication (TV,CD,MP3, etc.) Signal diagnostics (machine health, looking for hidden details)
2
A/D converters (analogue/digital)
This lecture, 3 concepts & 1 process.
Aliasing (sampling rate) A/D resolution Windowing (leakage)
Time domain FFT Frequency domain
3
Getting it into the computer! (A/D)
x(t)
time, t We now have a sequence of numbers, i.e., a vector.
4
Mention of one other issue: Resolution
voltage
For a specific voltage range there will only be a certain number of "divisions" available (see later).
Some definitions
x(t) Sampling rate fs=1/t t Window length or Time Record length, T.R.
t T.R.
6
Example of a A/D spec.
12bit A/D converter with a sampling rate of 15kHz over a range of -/+ 5 volts Voltage increments: Time increment:
voltage = 10/212 = 10/4096 = 2.441mV t = 1/15,000 sec.
7
Aliasing
How fast do we have to sample? What happens if we don't sample fast enough? What can we do about it?
Fast enough?
1Hz analogue signal 10Hz sampling frequency
9
Fast enough?
1Hz analogue signal 5Hz sampling frequency
10
1Hz analogue signal or a 4Hz analogue signal 5Hz sampling frequency
11
We must sample AT LEAST twice as fast as the highest frequency that is present. e.g. here, 4Hz signal therefore sample 8Hz.
12
To avoid aliasing, the analogue signal is first FILTERED (anti-aliasing filters) to ensure that all frequencies higher than 1/2 the sampling frequency have been removed. The maximum detectable frequency is sometimes called the Nyquist frequency.
13
Windowing
Finite record length, hence "stop" and "start"
T.R.
14
Windowing
15
Windowing function
16
Recap Fourier series.
Periodic signal of period ,
and where See eq(1.2.1-3)
17
Fourier series and FFT
Recall that for the following square wave: 1 x(t) t
18
Another way of viewing this:
x(t) t Time domain
FFT
Magnitude of Sine component
Frequency domain
frequency
19
Lab. Example
20
FFT summary
t Mag. 1/T.R. time
* * *
* *
frequency
0 T.R. FFT N data points
1/2t N/2 Mag. data points N/2 Phase. data points
21