Unit 4 - Digital Signal Processing - WWW - Rgpvnotes.in
Unit 4 - Digital Signal Processing - WWW - Rgpvnotes.in
Unit 4 - Digital Signal Processing - WWW - Rgpvnotes.in
UNIT IV
FFT algo ith s, deci atio i ti e algo ith , deci atio i f e ue c algo ith , deco positio fo N
composite number.
EFFICIENT COMPUTATION OF THE DFT
N-1
Xk = ∑ � −2� �/� , 0 ≤ k ≤ N-1
n=0
Let WN be the complex valued phase factor, which is an Nth root of unity expressed by,
WN = e j2 ∏ kn / N
Hence X (k) becomes
N-1
Xk = ∑ WNnk, 0 ≤ k ≤ N-1
n=0
RADIX-2 FFT
1. DECIMATION IN TIME (DITFFT)
N point sequence x (n) be decimated (broken) into two N/2 point data sequences f1(n) and f2(n). f1(n)
contains even numbered samples of x(n) and f2(n) contains odd numbered samples of x(n). This decimated
ope atio is called deci atio . “i ce it is do e o ti e do ai se ue ce it is called Decimation in Time .
Thus
N point DFT is given as
Since the sequence x (n) is decimated into even numbered and odd numbered samples, thus
X(k) = F1(k) + WNk F2(k)
N-1
Xk = ∑ WNnk, 0 ≤ k ≤ N-1
n=0
N/2-1 N/2-1
X(k) = ∑ x (2m) WN2mk + ∑ x (2m+1) WNk(2m+1)
m=0 m=0
x(0) X 0)
x(1) X 1)
x(2) X (2)
x(3) 8 Point DFT
X (3)
x(7) X (7)
EXAMPLE:
Given x (n) {0, 1, 2, 3}, find X (k) using DIT FFT algorithm.
SOLUTION:
Given N = 4
WNk = e-j(2π/N k
WN0 = 1 and WN1 = e-j(π/2 = -j
Using DIT FFT algorithm, X(k) from the given sequence x(n) is shown in figure below Therefore, X(k) =
{ 6, -2+2j, -2, -2, -2-j2}
x (0) = 0 2 X (0) = 6
W40=1 -1
x (1) = 1 4 X (1) = -2
W40=1 -1
x (3) = 3 -2 X (3) = -2-J2
W40=1 -1 W41=-j -1
N-1
X k = ∑ x (n) WNkn
n=0
N/2-1 N/2-1
Xk = ∑ WNkn + ∑ + N/2 WNk(n+N/2)
m=0 m=0
N/2-1 N/2-1
Xk = ∑ WNkn + WNkN/2 + ∑ + N/2 WNkn
m=0 m=0
N/2-1 N/2-1
Xk = ∑ WNkn + (-1) + ∑
k
+ N/2 WNkn
m=0 m=0
N/2-1
Xk = ∑ x (n) + (-1)k x(n + N/2) WNkn
N/2-1 m=0
Xk = ∑ x (n) + (-1)k x(n + N/2) WNkn
m=0
N/2-1
X(2k) = ∑ x (n) + (-1)2k x(n + N/2) WN2kn
m=0
N/2-1
X(2k) = ∑ x (n) + (-1)2k x(n + N/2) WN2kn
m=0
N/2-1
X(2k+1) = ∑ x (n)+(-1)(2k+1) x(n + N/2) WN(2k+1)n
m=0
a A= a + b
b B= (a –b)WNr
WNr
Fig 1. BUTTERFLY COMPUTATION
Fig 2 shows signal flow graph and stages for computation of radix-2 DIF FFT algorithm of N=4
x(0)
X(0)
A
x(1)
X(1)
w4 0
B
x(2)
X(2)
w4 0 C
x(3)
X(3)
w4 1 D w4 0
A composite or mixed radix FFT is used when N is a composite number which has more then one prime factor:
for example N = 6 or 12. For these cases also, efficient DIT and DIF algorithms can be developed. Here we give
a DIT FFT decomposition and hence N is written as a product of factors such as N = m 1, m2, ….. r.
If N =m1N1 where N = m1, m2, ….. r, the input sequence x(n) can be separated into m1 subsequences of N1
elements each. Then the DFT can be written as
DITFFT algorithms are based upon DIFFFT algorithms are based upon decomposition
decomposition of the input sequence into of the output sequence into smaller and smaller
smaller and smaller sub sequences. sub sequences.
In this input sequence x(n) is decimated into In this output sequence X(k) is considered to be
even and odd numbered samples decimated into even and odd numbered samples
Splitting operation is done on time domain Splitting operation is done on frequency domain
sequence. sequence.
In DIT FFT input sequence is in bit reversed In DIFFFT, input sequence is in natural order. And
order while the output sequence is in natural DFT should be read in bit reversed order.
order.