hc433r.m

某本matlab书籍中的代码

%《数字信号处理教程——MATLAB释义与实现》第四章例4.3.3重印后修改的程序hc433
% 截取有限长来计算无限长序列频谱的误差
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=32; dw=2*pi/N;        % 给定N，确定频率分辨率
n=0:N-1;x=0.7.^n;       % 列出时间序列x
X=fft(x); 				% 求出32点数据的DFT数值解
Xt=1.0./(1-0.7.*exp(-j*n*dw));	% 求出解析法得出的32点频谱
e=max(abs(abs(X)-abs(Xt))); % 求数值解和解析解在32个频点上的误差
Xm=max(abs(X));
pe=e/Xm*100				% 求出最大相对误差（百分数）