16点基4-FFT的程序
按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙x=;
X1=zeros(16,1);
W=;
%X1(k0,n0)利用4*k0+n0+1来用一维数组来记录
for n0=0:3 %具体的求解思路在课本P170
X1(n0+1)=W(1,:)*';
X1(n0+5)=W(2,:)*';
X1(n0+9)=W(3,:)*';
X1(n0+13)=W(4,:)*';
end
X2=zeros(16,1);
for k0=0:3%X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
X2(4*k0+1)=W(1,:)*';
X2(4*k0+2)=W(2,:)*';
X2(4*k0+3)=W(3,:)*';
X2(4*k0+4)=W(4,:)*';
end
% 倒位序
X3=zeros(16,1);
for k=1:16
k0=fix((k-1)/4);
k1=k-1-k0*4;
X3(k1*4+k0+1)=X2(k);
end
页:
[1]
