hc931.m

来自「这个是西安电子科技大学的陈怀琛教授编写的《数字信号处理教程——MATLAB释义与」· M 代码 · 共 25 行

%《数字信号处理教程——MATLAB释义与实现》(第2版)第九章例9.3.1程序hc931
% 系数量化误差对FIR滤波器特性的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
b=[0.069735,0.388726,0.360530, 0.388726, 0.069735];
[sos,G]=tf2sos(b,1);
V1=2^ceil(log2(max(abs(b)))),  			 % 求最大量化值
bq=bqtize(b,7,V1)
V2=2^ceil(log2(max(max(max(abs(sos))),G))), 	 % 求最大量化值
sosq=bqtize(sos,7,V2) ,Gq=bqtize(G,7,V2)
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),':'),
legend('理想系数','直接量化','级联量化')
grid on
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),':'),
legend('理想系数','直接量化','级联量化')
axis([1.8,2,0,0.2])
grid on
set(gcf,'color','w')                                                % 置图形背景色为白