fft_2.m

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

function  b = fft_2(a) %正序输出,逆序输出 频率抽取的递归算法
N=length(a);
c=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
a1=fft_2d(a);
for i=1:N %将输出的逆序转成正序
    b(i)=a1(c(i)+1);
end