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