MATLAB实现的FFT函数

写于09年08月29日，函数功能，用MATLAB实现的FFT函数。

 

完整程序 function xn=myfft(xn,N,M) mybitrevorder(xn,N); for m=1:M B=2^(m-1); for r=0:B-1 P=2^(M-m)*r; for k=r:2^m :N-1 f= xn(k+1); xn(k+1)= xn(k+1)+ xn(k+B+1)*exp(-i*pi*P/N); xn(k+B+1)= f-xn (k+B+1)*exp(-i*pi*P/N); end end end 倒序 function xn=mybitrevorder(xn,N) LH=N/2; r=LH; N1=N-2; for t=1:N1 if t<r T=xn(t+1); xn(t+1)=xn(r+1); xn(r+1)=T; end K=LH; while r>=K r=r-K; K=K./2; end r=r+K; end  

综合起来是：

fft程序 function xn=myfft(xn,N,M) LH=N/2; r=LH; N1=N-2; for t=1:N1 if t<r T=xn(t+1); xn(t+1)=xn(r+1); xn(r+1)=T; end K=LH; while r>=K r=r-K; K=K./2; end r=r+K; end for m=1:M B=2^(m-1); for r=0:B-1 P=2^(M-m)*r; for k=r:2^m :N-1 f= xn(k+1); xn(k+1)= xn(k+1)+ xn(k+B+1)*exp(-i*pi*P/N); xn(k+B+1)= f-xn(k+B+1)*exp(-i*pi*P/N); end end end