16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
" y5 f7 ~. C: W7 P! d6 L" J
3 \! w. c0 J2 z! S4 Q" b2 \" bx=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
) p" }4 e# n- E) S1 s- ^( w& a# jX1=zeros(16,1);
5 \! r7 n; M6 y8 V. E) i4 W9 PW=[1 1 1 1;1 -i -1 i;1 -1 1 -1;1 i -1 -i];
& N( p0 U# \: h2 t%X1(k0,n0)利用4*k0+n0+1来用一维数组来记录
for n0=0:3 %具体的求解思路在课本P170
2 d2 p# z6 \- n3 s5 R9 h   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
; `4 }0 L8 N2 }' _& U! B   X1(n0+5)=W(2,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
9 {$ L! L; `: D/ f; n3 A   X1(n0+9)=W(3,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
- h# D! P" G( R   X1(n0+13)=W(4,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
end
: Z+ k7 f" F/ }( d! W
X2=zeros(16,1);
" p; j7 N6 q) i% A9 x4 a, [for k0=0:3  %X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
8 ?% }+ Y% |; E, c( S! Z8 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)]';
: v" [4 S* N3 `: }   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)]';
* P+ {" \" t7 I; a/ M   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)]';
% o4 F1 X# i8 }) |! \+ U. ~2 P   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)]';
; }" u; C# z- u+ Zend
% 倒位序
! x6 G( u* m5 j% a, K$ Y+ C- b; nX3=zeros(16,1);
* Q5 f# b7 i# |8 S& k* W; B  \  
6 m7 h4 m8 k% W" Y1 i8 Lfor k=1:16
   k0=fix((k-1)/4);
   k1=k-1-k0*4;
7 j* R- g3 ^- z3 ]   X3(k1*4+k0+1)=X2(k);
end