0
1.5kviews
Develop a radix 2 DIF -FFT algorithm for computing 8-point DFT.
1 Answer
0
39views

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]=N1n=0x[n]WknNx[k]=N/21n=0x[n]wknN+N1n=N/2x[n]wknN.

 Consider Second Sum: N1n=N/2x[n]WknN substitute n=m+N/2n=N2;m=N2N2=0n=N1;m=N1N2=N21

 By change of variable Nn=0x[n+N2]Wk(n+π/2)kx[k]=N1n=0x[n]wknH+N1n=0x[n+N2]wk(n+N/2H=n21n=0x[n]wkns+vk+kN/l21n=0N][n+N/2]wknNVKN/2H=ej2πMN2k

=[ej1π]k=(1)kx[k]=N21n=0x[n]wknn+(1)kN21n=0x[n+N/2]wknN

Now let us split (decimate) x[k] into the even 2 odd-numbered samples.

x[2k]=N21n=0[x[n]+(1)2kx[n+N2]]w2knN(1)2k=1 \& W2knH=WknN/2

x[2k+1]=N21n=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]=Nn=0[[x[n]x[n+N/2]]WnN]WN/2k=0,1N21

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,1N21

then

G1[k]=x[2k]=N21n=0g1[n]wknN/2G2[k]=x[2k+1]=N21n=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.

enter image description here

Second Stage Decimation:-

G1[k]=x[2k]=N/1n=0N/2g1[n]wknN/2=N/4n=0g1[n]wknN/2+kn=N/4g1[n]wknN/2g11[n]=g1[n]+g1[n+N/4]g12[n]=[g1[n]g1[n+N/4]]w2nNG11[k]=G1[2k]=N/4n=0g11[n]wknN/4G12[k]=G1[2k+1]=N/21N=1g12[n]wknN/4G2[k]=x[2k+1]=N/21n=0g2[n]wknN/2

g21[n]=g2[n]+g2[n+N/4]g22[n]=[g2[n]g2[n+N/4]]W2nNG21[k]=G2[2k]=N/4n=0g21[n]wknN/4G22[k]=G2[2k+1]=N/4n=1g122[n]wknN/4

Third Stage Decimation:

enter image description here

G11[k]=N1kn=0g11[n]wknN/h

G11[k]=1n=0g11[0]wkn2G11[k]=g11[0]+wk2g11[1]G11[0]=g11[0]+w02g11[1]G1f[1]=g11[0]+w2gq1[1]G11[1]=g11[0]g11[1]

Please log in to add an answer.