n=0:1023;%截短
x=(1/2).^n;
w=-pi+n/1024*2*pi;
y=fftshift(fft(x)); 
y2=1./(1-0.5*exp(-i*w));
subplot(3,1,1);
plot(w,abs(y));
xlabel('頻譜分析');
legend('近似值');
xlabel('頻率（-pi——pi）');
subplot(3,1,2);
plot(w,abs(y2),'r');
legend('理論值');
xlabel('頻率（-pi——pi）');

subplot(3,1,3);
error=abs(y)-abs(y2);
plot(w,error);
xlabel('誤差');
xlabel('頻率（-pi——pi）');