3、转因子,每求出一个旋转因子,就计算完它对应的所有2M-L个蝶形。这样,我们可用三重循环程序实现DIT-FFT运算,程序框图如下图所示。其中,倒序部分的程序框图如下图所示:4.实验内容及步骤按照实验原理中的算法,用MATLAB编写实现算法的程序。5.实验程序及结果:程序代码:xn=[0,1,2,3,4,5,6,7]; N=length(xn); A=xn; %DIT_FFT NI=N/2; for I=1:N-1 if I
4、 end T=N/2; while NI>=T NI= NI-T; T=T/2; end NI= NI+T; end disp('逆序x[n]:'),disp(A); %butterfly WN=exp(-i*2*pi/N); v=floor(log2(N)); for m=1:v for k=0:2^m:N-1 for K=0:2^(m-1)-1 p=k+K; q=p+2^(m-1); r=2^(v-m)*mod(p,2^
5、m); B(p+1)=A(p+1)+A(q+1)*WN^r; B(q+1)=A(p+1)-A(q+1)*WN^r; end end A=B; end disp('FFT_X[k]:'),disp(A);运行结果: