dft_2.m

应用傅里叶变换DFT

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）');