ms16_1.m

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

Ts=1;N0=[3,5,10];
for r=1:3;
   N=N0(r);
	D=2*pi/(Ts*N);
	kn=floor(-(N-1)/2:-1/2);
	kp=floor(0:(N-1)/2);
	w=[kp,kn]*D;
   X=3+2*exp(-j*w)+exp(-j*2*w)+2*exp(-j*3*w)+3*exp(-j*4*w);
   n=0:N-1;
   x=abs(ifft(X,N))
   subplot(1,3,r);stem(n*Ts,x,'filled','k');
end