hc932.m

书籍中包含的第三部分实例的matlab代码

%《数字信号处理教程——MATLAB释义与实现》(第2版)第九章例9.3.2程序hc932
% 系数量化误差对IIR滤波器极点分布的影响
% 电子工业出版社出版  陈怀琛编著 2008年10月
%
b=[1,0.4,-0.03,0.232];              				% 分子系数向量
a=[1, -2.4760, 2.9541,-1.9631,0.4900];		% 分母系数向量
V=2^ceil(log2(max(abs([a,b]))))				% 求两系数向量最大量化值
aq=bqtize(a,7,V), bq=bqtize(b,7,V);   			% 量化后的分母、分子系数向量 
ra=(aq-a)./abs(a),                  	% 量化后的分母系数的相对误差
p=roots(a),pq=roots(aq),	        		% 分母系数量化前、后的极点
z=roots(b);zq=roots(bq);	        		% 分子系数量化前、后的零点
pause
[sos,G]=tf2sos(b,a);                				% 化为二阶级联形式
V2=2^ceil(log2(max(max(max(abs(sos))),G))) 	% 求系数矩阵的最大量化值
sosq=bqtize(sos,7,V2), Gq=bqtize(G,7,V2)		% 对各系数进行量化
a1q=sosq(1,[4:6]),a2q=sosq(2,[4:6])     		% 取出sos中的二阶分母系数向量 
p1q=roots(a1q), p2q=roots(a2q)				% 求系数向量的量化后的根