16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
5 J# I/ n! R1 T
x=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
# w0 W8 V. x( J" b/ ~& ^8 rX1=zeros(16,1);
0 ]4 [2 F& k5 ~& H  n7 IW=[1 1 1 1;1 -i -1 i;1 -1 1 -1;1 i -1 -i];
%X1(k0,n0)利用4*k0+n0+1来用一维数组来记录
for n0=0:3 %具体的求解思路在课本P170
   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
7 s9 _( N6 A: M0 ?+ S2 b   X1(n0+5)=W(2,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
- ]; n. o( |3 Y7 G   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)]';
- c% _8 \; U+ g/ N: k$ \- }end

8 z2 F  J- L0 X' f2 pX2=zeros(16,1);
for k0=0:3  %X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
4 h: k' N1 h: e+ Q8 X3 b   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)]';
# n' q; {. p; {( a$ g) d* c5 ?   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)]';
# Q3 k& W: ?$ K+ `   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)]';
   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
% i8 y  I' X& N' Y- w% 倒位序
/ O7 m) {% \3 r& D- q. y! ?X3=zeros(16,1);
  
  C$ k/ ~. o' z4 z2 M9 z  zfor k=1:16
   k0=fix((k-1)/4);
   k1=k-1-k0*4;
   X3(k1*4+k0+1)=X2(k);
end