ms14_4.m

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

Ts=0.5;C=[32,64,512];                   %输入不同的N值 
for r=0:2;
    N=C(r+1);
xn=[ones(1,6),zeros(1,N-11),ones(1,5)]; % 建立x(n)
D=2*pi/(N*Ts);
k=floor(-(N-1)/2:(N-1)/2);
X=fftshift(fft(xn,N));
subplot(3,2,2*r+1);plot(k*D,abs(X),'k');
subplot(3,2,2*r+2);stairs(k*D,angle(X),'k'); 
end