16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
2 o2 x9 V. X( ]8 _; w) b, u1 w* ?
) v7 g- w2 Z* v- r( ^x=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
4 B; w( A$ P3 j% n' N, UX1=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
! m& I  h, E# N- [2 X8 Y   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
: Y; \8 ]1 l! w, U* {9 N6 j   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)]';
# h0 S# }% \- o* v( C   X1(n0+13)=W(4,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
; w, _  d# K6 N- A2 n% send

' o# O- ]' L, E' y0 ?$ P3 d8 kX2=zeros(16,1);
for k0=0:3  %X2(k0,k1)利用4*k0+k1+1来用一维数组来记录，具体的计算矩阵在课本P171
+ o" Z9 o5 \# d7 G+ }   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" P' n9 u" L3 c  ^   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)]';
, K7 M$ _3 X* \% }& Z6 o, W0 k9 C! j   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)]';
, P0 I  N1 z/ s+ 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)]';
# B3 h/ I4 L2 W5 kend
% 倒位序
3 I' _; L- [$ b. {  mX3=zeros(16,1);
  
for k=1:16
; B" ^: ]/ g; I1 l, L   k0=fix((k-1)/4);
9 Z1 Q- H" E2 R$ t   k1=k-1-k0*4;
   X3(k1*4+k0+1)=X2(k);
  I/ L/ \; T9 u5 cend