%《數字信號處理教程——MATLAB釋義與實現》第九章例9.3.2程序hc932
% 系數量化誤差對IIR濾波器極點分布的影響
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
clear,close all
b=[1,0.4,-0.03,0.232];              % 分子系數向量
a=[1,-2.4760,2.9541,-1.9631,0.4900];% 分母系數向量
aq=brqtize(a,7), bq=brqtize(b,7);   % 量化后的分母、分子系數向量 
ra=(aq-a)./abs(a),                  % 量化后的分母系數的相對誤差
p=roots(a),pq=roots(aq),	        % 分母系數量化前、后的極點
z=roots(b);zq=roots(bq);	        % 分子系數量化前、后的零點
rp=abs(pq-p)./abs(p)                % 量化前、后極點的相對誤差
zplane(z,p), hold on,               % 畫出量化前零、極點位置(o和x)
plot(pq,'*'),plot(zq,'+'),          % 畫出量化后零、極點位置(+和*)
set(gcf,'color','w')                % 置圖形背景色為白
pause,
[sos,G]=tf2sos(b,a);                % 化為二階級聯形式
a1=sos(1,[4:6]),a2=sos(2,[4:6])     % 取出sos中的二階分母系數向量       
a1q=bqtize(a1,7);a2q=bqtize(a2,7);  % 量化后的二階分母系數向量 
p1q=roots(a1q),p2q=roots(a2q),      % 二階分母系數量化后的極點
abs(p1q),abs(p2q)                   % 極點的模
[H,w]=freqz(b,a);                   % 量化前的頻率特性
[Hq,w]=freqz(bq,aq);                % 量化后的頻率特性
figure(2),
subplot(2,1,1),plot(w,abs(H))
subplot(2,1,2),plot(w,abs(Hq))
set(gcf,'color','w')                 % 置圖形背景色為白