ms16_2.m

《数字信号处理实验(MATLAB版)》这

wm=62.8;Ts=pi/wm;
N0=[8,20];
for r=1:2;
   N=N0(r);
	D=2*pi/(Ts*N);
   k=[0:N-1]+eps;
   Omg=k*D;
   X=sin(0.275*Omg)./sin(0.025*Omg);
   n=0:N-1;
   x=abs(ifft(X,N))
   subplot(1,2,r);stem(n*Ts,x,'filled','k');
end