ms16r2.m

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

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