ms13_41.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 17 行

x1=[1,2,3,4,5,6,7,0,0,0];   %建立x(n),N=10序列
N=length(x1);
n=0:N-1;k=0:N-1;
y1=x1(mod(-n,N)+1);   %建立x(-n),N=7序列
Xk=x1*exp(-j*2*pi/N).^(n'*k) %求x(n)的DFT
Yk=y1*exp(-j*2*pi/N).^(n'*k) %求y(n)的DFT
subplot(2,2,1),stem(n,x1,'k');      %画x(n)
title('x(n)');
subplot(2,2,2),stem(n,y1,'k');      %画x(-n)
title('x(-n)');
subplot(2,2,3),stem(n,abs(Xk),'k'); %画x(n)的DFT
title('x(n)的DFT');
subplot(2,2,4),stem(n,abs(Yk),'k'); %画x(-n)的DFT
title('x(-n)的DFT');