16点基4-FFT的程序

zhaosl

按照书上的流程编的，但是最后的输出顺序不对，值是对了，就是顺序不对，求各位大牛帮个忙
7 m# P% ~. @. Q
5 g3 _% H4 F& {. m7 |x=[1 exp(2/2) exp(4/2) exp(6/2) 0 0 0 0 0 0 0 0 0 0 0 0];
+ }! D5 ]- ?+ n' e5 E- Z1 Z% F6 LX1=zeros(16,1);
$ H2 m/ _: R! a# a4 D0 k  zW=[1 1 1 1;1 -i -1 i;1 -1 1 -1;1 i -1 -i];
%X1(k0,n0)利用4*k0+n0+1来用一维数组来记录
6 }9 _, B8 `4 \for n0=0:3 %具体的求解思路在课本P170
- m% o$ h) T) `# I+ Z   X1(n0+1)=W(1,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';   
7 }; ~/ ]( G# C2 u   X1(n0+5)=W(2,:)*[x(n0+1),x(5+n0),x(9+n0),x(13+n0)]';
3 I& ^# _$ r1 u+ f/ f4 S* b   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)]';
% D! B2 m; p4 l, \5 |/ ]4 |end

X2=zeros(16,1);
for 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)]';
   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)]';
   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)]';
# [0 F# K8 U8 W- M3 t5 t) ~end
; x" ?' W" _, ^( p$ s  Q/ I% 倒位序
X3=zeros(16,1);
  
for k=1:16
   k0=fix((k-1)/4);
2 n1 C& d4 Q( I' x   k1=k-1-k0*4;
8 M0 i* d) Y& |, I" q) {7 E+ c   X3(k1*4+k0+1)=X2(k);
8 U' Z" O. ~) X! eend