fft_2d.m

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

function  b = fft_2d(a) %正序输出,逆序输出 频率抽取的递归算法
N=length(a);
if( N ~=2 )
    for k=1:N/2
        c(k)=a(k)+a(k+N/2);
        c(k+N/2)=( a(k)-a(k+N/2) )*exp(-j*2*pi/N*(k-1));
    end
    b(1:N/2)=fft_2d(c(1:N/2));
    b(N/2+1:N)=fft_2d(c(N/2+1:N));
else
    b(1)=a(1)+a(2);
    b(2)=a(1)-a(2);
end