N=100;
xn=[0,1,2,3,4,5,6,7,zeros(1,N-8)];
n=0:4N-1;k=0:4N-1;
xnl=xn(mod(n,N)+1);
Xk=xnlexp(-j2*pi/N).^(n'*k);
subplot(2,2,1),stem(xn);
title('原主值信号x(n)');
subplot(2,2,2),stem(n,xnl);
title('周期序列信号');
subplot(2,2,3),stem(k,abs(Xk));
title('|X(k)|');
subplot(2,2,4),stem(k,angle(Xk));
title('arg|X(k)|');