ms16_1n.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 12 行

Ts=1;N0=[3,5,10];
for r=1:3;
   N=N0(r);
	D=2*pi/(Ts*N);
	k=floor(-(N-1)/2:(N-1)/2);
	w=k*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