dft.m

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

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