ms12r41.m

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

xn=[7,6,5,4,3,2];
N=length(xn);
n=0:N-1;
w=linspace(-2*pi,2*pi,500);%将[-2π,2π]区间分割为500份
X=xn*exp(-j*n'*w);            %离散时间傅里叶变换
subplot(3,1,1),stem(n,xn);
ylabel('x(n)');
subplot(3,1,2),plot(w,abs(X));  %显示序列的幅度谱
axis([-2*pi,2*pi,1.1*min(abs(X)),1.1*max(abs(X))]);
ylabel('幅度谱');
subplot(3,1,3),plot(w,angle(X));%显示序列的相位谱
axis([-2*pi,2*pi,1.1*min(angle(X)),1.1*max(angle(X))]);
ylabel('相位谱');