李洁《数字信号处理》20058 ATLAB Ⅳ)编程思想及程序流程 function x= Myfft(x) [Xold, M, N]=Pinv(x)i 输入端第1级开始,逐级进 for I=l: M 行,共M级,所以L=1:M for J=0: B-1 第L级中每个蝶形 相距B=21个点 Np=exp(-1*2*pi*P/N); for kJJ: 2.L: N-1 x(k+1)=Xold(k+1)+Xold(k+B+l).*wNp; X(k+B+1)=Xo1d(k+1)-Xo1d(k+B+1) end k=1,81, 第L级共有21个不 Xold=X 同的旋转因子 end J=[0:B-们则旋转 因子指数分别为 个旋 P=2M-L*J 转因子的蝶 形相距2个 MATLAB function X= Myfft(x) TXold, M,N=Pinv(x); for l=l: M 八 for J=0: B-1 P=(2(ML)J; WNp=exp(-]*2*pi"P/N); for k=J: 2.L: N-1 X(k+l)=Xold(k+1)+Xold(k+B+l). * WNp; X(k+B+l)=Xold(k+1)-Xold(k+B+1). WNp; end XoldX d Digital Signal processing Jie Li 2005 7李洁《数字信号处理》2005® Digital Signal Processing__Jie Li 2005® 7 MATLAB Ⅳ）编程思想及程序流程 从输入端第1级开始，逐级进 行，共M级，所以L=1:M 第L级中每个蝶形 相距B=2L-1个点 第L级共有2L-1个不 同的旋转因子， J=[0:B-1]则旋转 因子指数分别为 P=2M-L*J 使用同一个旋 转因子的蝶 形相距2 L 个 点 function X = Myfft(x) [Xold,M,N]=Pinv(x); for L=1:M B=2.^(L-1); for J=0:B-1 P=(2.^(M-L)).*J; WNp=exp(-j*2*pi*P/N); for k=J:2.^L:N-1 X(k+1)=Xold(k+1)+Xold(k+B+1).*WNp; X(k+B+1)=Xold(k+1)-Xold(k+B+1).*WNp; end end Xold=X; end MATLAB function X = function X = Myfft(x Myfft(x) [Xold,M,N]=Pinv(x); for L=1:M B=2.^(L-1); for J=0:B-1 P=(2.^(M-L)).*J; L)).*J; WNp=exp(-j*2*pi*P/N); j*2*pi*P/N); for k=J:2.^L:N-1 X(k+1)=Xold(k+1)+Xold(k+B+1).* X(k+1)=Xold(k+1)+Xold(k+B+1).*WNp; X(k+B+1)=Xold(k+1) X(k+B+1)=Xold(k+1)-Xold(k+B+1).* Xold(k+B+1).*WNp; end end Xold=X; end