dft.m

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

function b = dft(a)
N=length(a);
b=zeros(1,N);
for k=1:N
    for n=1:N
        b(k)=b(k)+a(n)*exp(-2*pi*j*(k-1)*(n-1)/N);
    end
end