16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
5 x3 ~5 l7 r9 }, F
x=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
: `  I9 Y9 t! d" p* [X1=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来用一维数组来记录
for n0=0:3 %具体的求解思路在课本P170
, B' L  O8 N/ E   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
% F3 U0 V4 x" n  d6 V/ |   X1(n0+5)=W(2,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
   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
. g4 N) _& ^. ?& D
X2=zeros(16,1);
, `  C+ q) [! N$ p5 F; O; N$ Tfor k0=0:3  %X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
   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)]';
1 c4 ~" z# G, z5 m6 d, h   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)]';
9 \, P% h/ |- g8 m6 g   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)]';
- ?: A& I" N$ b   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)]';
6 K7 u, e  I% wend
, \" W1 J& \9 ^8 l2 w2 [% 倒位序
X3=zeros(16,1);
+ |' P/ u  v2 ~0 T. h! K  
/ o( b) @! @! W. C4 i9 Q; h- Cfor k=1:16
0 a9 b  G& k- E6 ?: i6 {   k0=fix((k-1)/4);
1 |$ W: X" V: ?) s) n+ J5 t- N* F7 n   k1=k-1-k0*4;
- {  Y. E) H3 X, u3 U   X3(k1*4+k0+1)=X2(k);
end