dft.m

实现快速傅里叶变换

function y=DFT(x)
N=length(x)
k=0:N-1
W=exp((-2*pi*j/N)*(k'*k))
y=x*W
plot(k,abs(y),'*')
return