dft.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 8 行

function[Xk]=dft(xn,N)
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;