dft.m

来自「自己编写的自适应notch fileter算法」· M 代码 · 共 11 行

function[Xk]=dft(xn,N)
%计算离散傅立叶变换
%Xk 0<=n<=N-1,xn为N点有限持续时间序列

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