%《數字信號處理教程——MATLAB釋義與實現》第九章例9.3.1程序hc931
% 系數量化誤差對FIR濾波器特性的影響
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
b=[0.069735,0.388726,0.360530, 0.388726, 0.069735];                 % 原始的系數向量
[sos,G]=tf2sos(b,1);                                                % 原始二階級聯結構的系數向量
bq=brqtize([0.069735,0.388726,0.360530, 0.388726, 0.069735],4)      % 量化后的系數向量
sosq=brqtize(sos,4) ,Gq=brqtize(G,4)                                % 量化后的二階級聯結構的系數向量
subplot(1,2,1),
[H1,w1]=freqz(b,1);plot(w1,abs(H1)),hold on                         % 原始系統的頻率特性
[H2,w2]=freqz(bq,1);plot(w2,abs(H2),'-.'),                          % 量化后的系統的頻率特性
[H3,w3]=freqz(conv(sosq(1,1:3),sosq(2,1:3)),1);plot(w3,Gq*abs(H3),':'),% 量化后的二階級聯結構的頻率特性
subplot(1,2,2),                              
[H1,w1]=freqz(b,1);plot(w1,abs(H1)),hold on                         % 圖同上，只改變比例尺
[H2,w2]=freqz(bq,1);plot(w2,abs(H2),'-.'),
[H3,w3]=freqz(conv(sosq(1,1:3),sosq(2,1:3)),1);plot(w3,Gq*abs(H3),':'),
axis([1.8,2.2,0,0.2])
legend('理想系數','直接量化','級聯量化')
set(gcf,'color','w')                                                % 置圖形背景色為白