Title: The Discrete Fourier Transform
1The Discrete Fourier Transform
2Content
- Introduction
- Representation of Periodic Sequences
- DFS (Discrete Fourier Series)
- Properties of DFS
- The Fourier Transform of Periodic Signals
- Sampling of Fourier Transform
- Representation of Finite-Duration Sequences
- DFT (Discrete Fourier Transform)
- Properties of the DFT
- Linear Convolution Using the DFT
3The Discrete Fourier Transform
4Signal Processing Methods
Fourier Series
Continuous-Time Fourier Transform
DFS
Discrete-Time Fourier Transform and z-Transform
5Frequency-Domain Properties
Continuous Aperiodic
Discrete Aperiodic
Continuous Periodic (2?)
?
6Frequency-Domain Properties
Relation?
Relation?
Relation?
Relation?
7Frequency-Domain Properties
Discrete Periodic
8The Discrete Fourier Transform
- Representation of Periodic Sequences --- DFS
9Periodic Sequences
- Notation a sequence with period N
where r is any integer.
10Harmonics
Facts
Each has Periodic N.
N distinct harmonics e0(n), e1(n),, eN?1(n).
11Synthesis and Analysis
Notation
Both have Period N
Synthesis
Analysis
12Example
A periodic impulse train with period N.
13Example
14Example
15Example
16DFS vs. FT
17Example
18Example
19Example
20The Discrete Fourier Transform
21Linearity
22Shift of A Sequence
Change Phase (delay)
23Shift of Fourier Coefficient
Modulation
24Duality
25Periodic Convolution
Both have Period N
26Periodic Convolution
Both have Period N
27The Discrete Fourier Transform
- The Fourier Transform of Periodic Signals
28Fourier Transforms of Periodic Signals
Sampling
Sampling
29Fourier Transforms of Periodic Signals
30The Discrete Fourier Transform
- Sampling of
- Fourier Transform
31Equal Space Sampling of Fourier Transform
32Equal Space Sampling of Fourier Transform
33Equal Space Sampling of Fourier Transform
34Example
35Example
Time-Domain Aliasing
36Time-Domain Aliasing vs. Frequency-Domain Aliasing
- To avoid frequency-domain aliasing
- Signal is bandlimited
- Sampling rate in time-domain is high enough
- To avoid time-domain aliasing
- Sequence is finite
- Sampling interval (2?/N) in frequency-domain is
small enough
37DFT vs. DFS
- Use DFS to represent a finite-length sequence is
call the DFT (Discrete Fourier Transform). - So, we represent the finite-duration sequence by
a periodic sequence. One period of which is the
finite-duration sequence that we wish to
represent.
38The Discrete Fourier Transform
- Representation of
- Finite-Duration Sequences --- DFT
39Definition of DFT
Synthesis
Analysis
40Example
41Example
42The Discrete Fourier Transform
43Linearity
44Circular Shift of a Sequence
45Circular Shift of a Sequence
46Duality
47Example
Choose N10
Rex1(n) ReX(n)
ReX(k)
Imx1(n) ImX(n)
ImX(k)
X1(k) 10x((?k))10
48Linear Convolution (Review)
49Circular Convolution
both of length N
50Example
n02, N5
51Example
n02, N7
52Example
LN6
53Example
54The Discrete Fourier Transform
- Linear Convolution Using the DFT
55Why Using DFT for Linear Convolution?
- FFT (Fast Fourier Transform) exists.
- But., we have to ensure that circular convolving
nature of DFT gives the linear convolving result.
56The Procedure
- Let x1(n) and x2(n) be two sequences of length L
and P, respectively. - Compute N-point (N ?) DFTs X1(k) and X2(k).
- Let X3(k) X1(k) X2(k), 0 ? k ? N?1.
- Let x3(n) DFT?1X3(k) x1(n) ? x2(n).
x1(n) x2(n) x1(n) ? x2(n)?
57Linear Convolution ofTwo Finite-Length Sequences
x3(?1) 0
x3(n) ? 0
n 0, 1, ?, LP?2
x3(LP?1) 0
58Linear Convolution ofTwo Finite-Length Sequences
Length of x3(n) x1(n)x2(n) LP?1
59Circular Convolution as Linear Convolution with
Time Aliasing
x1(n) ?? Length L x2(n) ?? Length P
x(n) x1(n)x2(n) ?? Length LP?1
60Circular Convolution as Linear Convolution with
Time Aliasing
61N LP?1
62For Finite Sequences
x p(n) x1(n)x2(n) x1(n)?x2(n), 0 ?
n ? LP?2
63N L
64N L
Corrupted (P?1) points
Uncorrupted (L?P1) points
65FIR Filter for Indefinite-Length Signals
h (n)
x (n)
- Overlap-Add Method
- Overlap-Save Method
Block Convolution
66Overlay-Add Method
67Overlay-Add Method
68Overlay-Add Method
Set N LP?1 for each block convolution
69Overlay-Save Method
- Each block is of length L.
- P?1 samples are overlaid btw. adjacent blocks.
- Set N LP?1 for each block convolution.
- Save the last L?P1 values for each block
convolution.