ms14_3_1.m

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

Fs=10;N=1000;
xn=[1,2,3,2,1];Nx=length(xn);
xn=[1,2,3,2,1,zeros(1,N-Nx-1)];
D=2*pi*Fs/N;
k=floor(-(N-1)/2:(N-1)/2);
X=fftshift(fft(xn,N));
subplot(1,2,1);plot(k*D,abs(X),'k');
axis([-40,40,-1,10]);
title('幅度频谱');xlabel('rad/s');
subplot(1,2,2);plot(k*D,angle(X),'k'); 
title('相位频谱');xlabel('rad/s');