ms13r3.m

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

x1=[8,7,6,5,4,3];     %建立x(n)序列
N1=length(x1);N=10;
x1=[x1,zeros(1,N-N1)]
n=0:N-1;k=0:N-1;
y1=x1(mod(-n,N)+1);   %建立x(-n),N=10序列
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);      %画x(n)
title('x(n)');
subplot(2,2,2),stem(n,y1);      %画x(-n)
title('x(-n)');
subplot(2,2,3),stem(n,abs(Xk)); %画|X(k)|
title('|X(k)|');
subplot(2,2,4),stem(n,abs(Yk)); %画|Y(k)|
title('|Y(k)|');