Transcript Document
The Discrete Fourier Transform
Quote of the Day
Mathematics may be defined as the subject in which
we never know what we are talking about, nor
whether what we are saying is true.
Bertrand Russell
Content and Figures are from Discrete-Time Signal Processing, 2e by Oppenheim, Shafer, and Buck, ©1999-2000 Prentice Hall
Inc.
Sampling the Fourier Transform
• Consider an aperiodic sequence with a Fourier transform
DTFT
x[n]
X ej
• Assume that a sequence is obtained by sampling the DTFT
~
Xk X e j
X e j2 / Nk
2 / Nk
• Since the DTFT is periodic resulting sequence is also periodic
• We can also write it in terms of the z-transform
~
Xk Xz z e2 / N k X e j2 / Nk
•
•
•
The sampling points are shown in figure
~
X k could be the DFS of a sequence
Write the corresponding sequence
1 N1 ~
~
x[n] Xk e j2 / Nkn
N k 0
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
2
Sampling the Fourier Transform Cont’d
• The only assumption made on the sequence is that DTFT exist
xme
Xe
j
~
Xk X e j2 / Nk
jm
m
1 N1 ~
~
x[n] Xk e j2 / Nkn
N k 0
• Combine equation to get
1 N 1
~
x[n] xme j2 / Nkm e j2 / Nkn
N k 0 m
1 N 1 j2 / Nk n m
~
xm e
x
m
p n m
m
N k 0
m
• Term in the parenthesis is
• So we get
1 N1 j2 / Nk nm
~
p n m e
n m rN
N k 0
r
~
x[n] xn
Copyright (C) 2005 Güner Arslan
r
r
n rN xn rN
351M Digital Signal Processing
3
Sampling the Fourier Transform Cont’d
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
4
Sampling the Fourier Transform Cont’d
• Samples of the DTFT of an aperiodic sequence
– can be thought of as DFS coefficients
– of a periodic sequence
– obtained through summing periodic replicas of original sequence
• If the original sequence
– is of finite length
– and we take sufficient number of samples of its DTFT
– the original sequence can be recovered by
~
x n 0 n N 1
xn
else
0
• It is not necessary to know the DTFT at all frequencies
– To recover the discrete-time sequence in time domain
• Discrete Fourier Transform
– Representing a finite length sequence by samples of DTFT
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
5
The Discrete Fourier Transform
• Consider a finite length sequence x[n] of length N
xn 0 outside of 0 n N 1
• For given length-N sequence associate a periodic sequence
~
x n
xn rN
r
• The DFS coefficients of the periodic sequence are samples of
the DTFT of x[n]
• Since x[n] is of length N there is no overlap between terms of
x[n-rN] and we can write the periodic sequence as
~
xn xn mod N xnN
• To maintain duality between time and frequency
~
– We choose one period of X k as the Fourier transform of x[n]
~
~
Xk 0 k N 1
Xk
Xk Xk mod N Xk N
0
else
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
6
The Discrete Fourier Transform Cont’d
• The DFS pair
~
Xk
N 1
~
x[n]e j2 / Nkn
n0
1 N1 ~
~
x[n] Xk e j2 / Nkn
N k 0
• The equations involve only on period so we can write
N1 ~
x[n]e j2 / Nkn 0 k N 1
Xk n 0
0
else
1 N1 ~
Xk e j2 / Nkn 0 k N 1
x[n] N k 0
0
else
• The Discrete Fourier Transform
Xk
N1
x[n]e
n0
j2 / Nkn
1 N1
x[n] Xk e j2 / Nkn
N k 0
• The DFT pair can also be written as
DFT
Xk
x[n]
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
7
Example
• The DFT of a rectangular pulse
• x[n] is of length 5
• We can consider x[n] of any
length greater than 5
• Let’s pick N=5
• Calculate the DFS of the periodic
form of x[n]
~
Xk
4
j2 k / 5 n
e
n0
1 e j2 k
1 e j2 k / 5
5 k 0,5,10,...
else
0
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
8
Example Cont’d
• If we consider x[n] of length 10
• We get a different set of DFT
coefficients
• Still samples of the DTFT but in
different places
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
9
Properties of DFT
• Linearity
x1 n
DFT
X1 k
x2 n
DFT
X2 k
ax1 n bx2 n DFT
aX1 k bX2 k
• Duality
DFT
xn
Xk
DFT
Xn
Nx k N
• Circular Shift of a Sequence
xn
DFT
Xk
DFT
xn mN 0 n N - 1
Xk e j2k / Nm
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
10
Example: Duality
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
11
Symmetry Properties
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
12
Circular Convolution
• Circular convolution of of two finite
length sequences
x3 n
x3 n
N 1
x mx n m
m 0
1
2
N
N 1
x mx n m
m 0
2
Copyright (C) 2005 Güner Arslan
1
N
351M Digital Signal Processing
13
Example
• Circular convolution of two rectangular pulses L=N=6
1 0 n L 1
x1 n x2 n
else
0
• DFT of each sequence
X1 k X2 k
N1
e
j
2
kn
N
n0
N k 0
0 else
• Multiplication of DFTs
N2 k 0
X3 k X1 k X2 k
0 else
• And the inverse DFT
N 0 n N 1
x3 n
else
0
Copyright (C) 2005 Güner Arslan
351M Digital Signal Processing
14
Example
• We can augment zeros to
each sequence L=2N=12
• The DFT of each sequence
X1 k X2 k
1e
j
1e
2 Lk
N
j
2 k
N
• Multiplication of DFTs
2 Lk
j
1 e N
X3 k
2 k
j
1e N
Copyright (C) 2005 Güner Arslan
2
351M Digital Signal Processing
15