hc752a.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.5.2程序hc752a
% （逐步逼近）用等波动法设计低通滤波器2
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=11;L=floor((N-1)/2);              % 滤波器及方程阶数
w0=[0, 0.65, 1, 1.5, 1.8, 2.4,  3.14]; % 初设频率向量
D=[1,1,1,0,0,0,0];                  % 预期幅频特性
P=cos(w0'*[0:L]);                   % 符幅特性谐波向量
Wt=[0.5,0.5,0.5,1,1,1,1];           % 权重函数
WT=(-1).^[0:6]./Wt;                 % 求系数中的WT项
Pa=[P,WT'];                          % 最小最大方程系数向量
d=Pa\ D'                            % 求解最小最大方程
h=[flipud(d(2:6)/2)',d(1),d(2:6)'/2]  % 按类型I的FIR滤波器由d求h
[H,w]=freqz(h,1);				% 由脉冲响应求频率特性
plot(w,abs(H)),hold on					% 画出实际的幅频特性
H0=freqz(h,1,w0);abs(H0)-D		% 求给定频点上的幅特性，并求点上的误差
plot(w0,abs(H0),'.')   		% 画出设定频点上的幅特性