fft_1.m

分四种情况实现fft 分别是时间抽取递归 、时间抽取循环、频率抽取递归和频率抽取循环 以及附带dft的算法描述

function  b = fft_1(a) %逆序输出,正序输出 时间抽取的递归算法
N=length(a);
c=zeros(1,N);
x1=zeros(1,N);
c(1)=0;c(2)=N/2; %c(i)为逆序的序号
for i=1:(log2(N)-1)   
     c((2^i+1):2^(i+1))=c(1:2^i)+2^(log2(N)-1-i);
end

for i=1:N  %实现逆序输入
    x1(i)=a(c(i)+1);
end
tic;
b=fft_1d(x1);%递归调用fft_1-d()函数