%《數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)》(第2版)第九章例9.3.2程序hc932
% 系數(shù)量化誤差對(duì)IIR濾波器極點(diǎn)分布的影響
% 電子工業(yè)出版社出版  陳懷琛編著 2008年10月
%
b=[1,0.4,-0.03,0.232];              				% 分子系數(shù)向量
a=[1, -2.4760, 2.9541,-1.9631,0.4900];		% 分母系數(shù)向量
V=2^ceil(log2(max(abs([a,b]))))				% 求兩系數(shù)向量最大量化值
aq=bqtize(a,7,V), bq=bqtize(b,7,V);   			% 量化后的分母、分子系數(shù)向量 
ra=(aq-a)./abs(a),                  	% 量化后的分母系數(shù)的相對(duì)誤差
p=roots(a),pq=roots(aq),	        		% 分母系數(shù)量化前、后的極點(diǎn)
z=roots(b);zq=roots(bq);	        		% 分子系數(shù)量化前、后的零點(diǎn)
pause
[sos,G]=tf2sos(b,a);                				% 化為二階級(jí)聯(lián)形式
V2=2^ceil(log2(max(max(max(abs(sos))),G))) 	% 求系數(shù)矩陣的最大量化值
sosq=bqtize(sos,7,V2), Gq=bqtize(G,7,V2)		% 對(duì)各系數(shù)進(jìn)行量化
a1q=sosq(1,[4:6]),a2q=sosq(2,[4:6])     		% 取出sos中的二階分母系數(shù)向量 
p1q=roots(a1q), p2q=roots(a2q)				% 求系數(shù)向量的量化后的根