hc931.m

数字信号处理,对与初雪着很有用,注释非常详细,调试都通过了!

%《数字信号处理教程——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')                                                % 置图形背景色为白