fft64.m

基于蝶形运算的64点快速傅立叶变换

function y=fft64(x);
% x=zeros(1,64);
% x(1)=1;
x1=reshape(x,8,8);
y=zeros(8);
y1=zeros(8);
y2=zeros(8);
y3=zeros(8);
for i=0:7
    y1(:,i+1)=fft64_8([x(i*4+1:i*4+4);x(32+i*4+1:32+i*4+4)]);
end

t1=y1(1:4,:);
t2=y1(5:8,:);
z1=[t1(:);t2(:)];

for i=0:7
    y2(:,i+1)=fft64_4([z1(i*2+1:i*2+2);z1(16+i*2+1:16+i*2+2);z1(32+i*2+1:32+i*2+2);z1(48+i*2+1:48+i*2+2)]);
end
t1=y2(1:2,:);
t2=y2(3:4,:);
t3=y2(5:6,:);
t4=y2(7:8,:);
z2=[t1(:);t2(:);t3(:);t4(:)];


for i=0:7
    y3(:,i+1)=fft64_2([z2(i+1:8:i+57)]);
end
y3=y3';
z3=y3(:);

wn=exp(-j*2*pi/64);
% ww=wn*ones(8,1);
ww=wn.^([0:7]'*[0:7]);
www=ww(:,[1 5 3 7 2 6 4 8]);
% www=ww*[0:7];
www=www(:);

z3=z3.*www;

for i=0:7;
    y(:,i+1)=fft2_4_8(z3(8*i+1:8*i+8));
end
a=[ 0 32 16 48 8 40 24 56]+1;
a=a';
a=[a a+4 a+2 a+6 a+1 a+5 a+3 a+7];
% y=y';
% y=y(:,[0 4 2 6 1 5 3 7]+1);
% y=y';
y=y(a(:));