ms11r3.m

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

xn1=[1,1,0,0,0,0,0,0,0,0];
N=length(xn1);
xn=[xn1,xn1,xn1,xn1];
n=0:4*N-1;
k=0:4*N-1;
Xk=xn*exp(-j*2*pi/N).^(n'*k);       %离散傅里叶级数变换
subplot(2,2,1),stem(xn1,'k');
title('x(n)');axis([0,N,1.1*min(xn),1.1*max(xn)]);
subplot(2,2,2),stem(n,xn,'k');%显示逆变换结果
title('x(n+4r)');
axis([-1,4*N,1.1*min(xn),1.1*max(xn)]);
subplot(2,2,3),stem(k,abs(Xk),'k');%显示序列的幅度谱
title('|X(k)|');
axis([-1,4*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]);
subplot(2,2,4),stem(k,angle(Xk),'k');%显示序列的相位谱
title('arg|X(k)|');
axis([-1,4*N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]);