fft_1d.m

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

function  b = fft_1d(a) %逆序输出,正序输出 时间抽取的递归算法
N=length(a);

if( N ~=2 )
    c(1:N/2)=fft_1d(a(1:N/2));
    c(N/2+1:N)=fft_1d(a(N/2+1:N));
    for k=1:N/2
        b(k)=c(k)+c(k+N/2)*exp(-j*2*pi/N*(k-1));
        b(k+N/2)=c(k)-c(k+N/2)*exp(-j*2*pi/N*(k-1));
    end
else
    b(1)=a(1)+a(2);
    b(2)=a(1)-a(2);
end