ms14_3.m

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

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