idft.m

自己编写的自适应notch fileter算法

function[Xk]=idft(xn,N)
%计算离散傅立叶ni变换
%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/N