function ya=neichahanshu(xa,Na,wa)
ka=0:Na-1;
XK1a=fft(xa,Na);
for W = 1:length(wa)
num1=1-exp(-j*wa(W)*Na);
num2=Na*(1-exp(-j*(wa(W)-2*pi*ka/Na)));
num3=num1./num2;
y(W)=sum(XK1a.*num3);
ya=abs(y);
end
end
定义dft重构dtft内插函数
最新推荐文章于 2024-07-07 01:39:48 发布