ms14r23.m

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

Ts=0.2;
C=[8,32,64];                 %输入不同的N值 
for r=0:2;
   N=C(r+1);
   n=0:N-1;
   xn=0.5.^n;   %建立x(n)
	D=2*pi/(N*Ts);
	k=floor(-(N-1)/2:(N-1)/2);
	X=fftshift(fft(xn,N));
	subplot(3,2,2*r+1);plot(k*D,abs(X),'k');
	subplot(3,2,2*r+2);stairs(k*D,angle(X),'k');
end