用Matlab编写fft

FFT算法编写1、程序function[A]=myfft(A,M)N=2^M;LH=N/2;J=LH;N1=N-2;forI=1:1:N1ifIJT=A(I+1);A(I+1)=A(J+1);A(J+1)=T;endK=LH;whileJ=KJ=J-K;K=K/2;endJ=J+K;endforL=1:1:MB=2^(L-1);forJ=0:B-1p=J*2^(M-L);fork=J:2^L:N-1T=A(k+1)+A(k+B+1)*exp(-i*2*pi*p/N);A(k+B+1)=A(k+1)-A(k+B+1)*exp(-i*2*pi*p/N);A(k+1)=T;endendendend2、程序检验：在Matlab中输入如下代码，将自己编写的FFT算法与Matlab自带算法的计算结果进行比较。A=[01234567];myfft(A,3)fft(A,8)计算结果如下图：由结果可以看出，编写的myFFT算法得到的结果与Matlab自带算法得到的结果相同，说明编写正确。