dft.m

自制的快速傅立叶变换函数

function sw=dft(st)
% DFT function
M=length(st);
N=M;
WN=exp(2*pi*j/N);

%main loop
for n=0:N-1
    temp=0;
    for m=0:M-1
        temp=temp+st(m+1)*WN^(-n*m);
    end
    sw(n+1)=temp;
end