%《數字信號處理教程——MATLAB釋義與實現》第四章例4.3.3程序hc433
% 截取有限長來計算無限長序列頻譜的誤差
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
T=0.5; N=32; dw=2*pi/(N*T); % 給定N，確定頻率分辨率
n=0:N-1;x=0.7.^n;       % 列出時間序列x
X=fft(x); 				% 求出32點數據的DFT數值解
Xt=1.0./(1-0.7.*exp(-j*0.5.*n*dw));	% 求出解析法得出的32點頻譜
e=max(abs(abs(X)-abs(Xt))); % 求數值解和解析解在32個頻點上的誤差
Xm=max(abs(X));
pe=e/Xm*100				% 求出最大相對誤差（百分數）