written 2.4 years ago by | • modified 2.4 years ago |
Solution:
This algorithm is obtained by using the divide \& conquer approach.
To derive this algorithm we begin with 2 splitting the DFT formula into hoo summation g of which one involves the sum over the first N/2 data points and the second sum over the last N/2 data points.
First Stage Decimation:-
By the definition of DFT.
x[k]=N−1∑n=0x[n]WknNx[k]=N/2−1∑n=0x[n]wknN+N−1∑n=N/2x[n]wknN.
Consider Second Sum: N−1∑n=N/2x[n]WknN substitute n=m+N/2∴n=N2;m=N2−N2=0. n=N−1;m=N−1−N2=N2−1
By change of variable ∴N∑n=0x[n+N2]Wk(n+π/2)kx[k]=N−1∑n=0x[n]wknH+N−1∑n=0x[n+N2]wk(n+N/2H=n2−1∑n=0x[n]wkns+vk+kN/l2−1n=0N][n+N/2]wknNVKN/2H=e−j2πM⋅N2⋅k
=[e−j1π]k=(−1)kx[k]=N2−1∑n=0x[n]wknn+(−1)kN2−1∑n=0x[n+N/2]wknN
Now let us split (decimate) x[k] into the even 2 odd-numbered samples.
x[2k]=N2−1∑n=0[x[n]+(−1)2kx[n+N2]]w2knN(−1)2k=1 \& W2knH=WknN/2
x[2k+1]=N2−1∑n=0[x[n]+(−1)2k+1x[n+N2]]W(2k+1)nN(−1)2k+1=−1,W(2k+1)nN=w2knNWnN=WknN/2WnNx[2k+1]=N∑n=0[[x[n]−x[n+N/2]]WnN]WN/2k=0,1…N2−1
Let us define N2 point sequences g1[n]&g2[n] as
g1[n]=x[n]+x[n+N2]g2[n]=[x[n]−x[n+N2]]WNNn=0,1…N2−1
then
G1[k]=x[2k]=N2−1∑n=0g1[n]wknN/2G2[k]=x[2k+1]=N2−1∑n=0g2[n]wknN/2
At this stage, we have decimated the sequence of N point DFT into two N2 points BETs. given by eqn(7)&8
Let us consider N=8 &\ N=4
Here we get even and odd 0/p of x[k].
Thus decimation of sequence in frequency. Therefore it is called Decimation in Frequency Radix - 2 FFT algorithm.
Second Stage Decimation:-
G1[k]=x[2k]=N/1∑n=0N/2g1[n]wknN/2=N/4∑n=0g1[n]wknN/2+k∑n=N/4g1[n]wknN/2g11[n]=g1[n]+g1[n+N/4]g12[n]=[g1[n]−g1[n+N/4]]w2nN−G11[k]=G1[2k]=N/4∑n=0g11[n]wknN/4G12[k]=G1[2k+1]=N/2−1∑N=1g12[n]wknN/4G2[k]=x[2k+1]=N/2−1∑n=0g2[n]wknN/2
g21[n]=g2[n]+g2[n+N/4]g22[n]=[g2[n]−g2[n+N/4]]W2nN−G21[k]=G2[2k]=N/4∑n=0g21[n]wknN/4G22[k]=G2[2k+1]=N/4∑n=1g−122[n]wknN/4
Third Stage Decimation:
G11[k]=∑N−1kn=0g11[n]wknN/h
G11[k]=1∑n=0g11[0]wkn2G11[k]=g11[0]+wk2g11[1]G11[0]=g11[0]+w02g11[1]G1f[1]=g11[0]+w′2gq1[1]G11[1]=g11[0]−g11[1]