fft_1.m

来自「分四种情况实现fft 分别是时间抽取递归 、时间抽取循环、频率抽取递归和频率抽」· M 代码 · 共 14 行

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()函数