16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
) f& ?- T9 G6 g1 D2 r
x=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
3 |+ F! i: v' y+ x7 [2 X9 Z" FX1=zeros(16,1);
W=[1 1 1 1;1 -i -1 i;1 -1 1 -1;1 i -1 -i];
%X1(k0,n0)利用4*k0+n0+1来用一维数组来记录
' r' T1 ~' A' W. u% yfor n0=0:3 %具体的求解思路在课本P170
0 W6 r/ G: z1 a- _- ?& S" g+ k! D   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
$ {4 o4 z% V, c1 b2 X8 L   X1(n0+5)=W(2,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
- t' ?: \9 S) `8 G2 @   X1(n0+9)=W(3,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
   X1(n0+13)=W(4,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
end

6 }4 J0 E0 p6 N9 {% r" _X2=zeros(16,1);
for k0=0:3  %X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
  V7 }0 ?7 F% E   X2(4*k0+1)=W(1,:)*[X1(4*k0+1),X1(4*k0+2)*exp(-i*pi/8*k0),X1(4*k0+3)*exp(-i*pi/4*k0),X1(4*k0+4)*exp(-i*pi/8*3*k0)]';
" E# ?: t& |! A* N   X2(4*k0+2)=W(2,:)*[X1(4*k0+1),X1(4*k0+2)*exp(-i*pi/8*k0),X1(4*k0+3)*exp(-i*pi/4*k0),X1(4*k0+4)*exp(-i*pi/8*3*k0)]';
   X2(4*k0+3)=W(3,:)*[X1(4*k0+1),X1(4*k0+2)*exp(-i*pi/8*k0),X1(4*k0+3)*exp(-i*pi/4*k0),X1(4*k0+4)*exp(-i*pi/8*3*k0)]';
- V; t. ]+ y- J, T: e$ ~( w   X2(4*k0+4)=W(4,:)*[X1(4*k0+1),X1(4*k0+2)*exp(-i*pi/8*k0),X1(4*k0+3)*exp(-i*pi/4*k0),X1(4*k0+4)*exp(-i*pi/8*3*k0)]';
end
% 倒位序
X3=zeros(16,1);
6 E/ X2 O6 {1 K1 x: q3 k" O) Y  
; j/ l1 s) x4 ?. D! r* Efor k=1:16
2 [+ w1 N2 X" ?4 V& N   k0=fix((k-1)/4);
   k1=k-1-k0*4;
   X3(k1*4+k0+1)=X2(k);
end