program_9_2.m

这是国外一本经典教材附带的程序

% Program 9_2
% Coefficient Quantization Effects on the
% Frequency Response of a Cascade Form IIR Filter
%
clf;
[z,p,k] =3D ellip(5,0.4,50,0.4);
[b,a] =3D zp2tf(z,p,k);
[h,w] =3D freqz(b,a,512); g =3D 20*log10(abs(h));
sos =3D zp2sos(z,p,k);
sosq =3D a2dT(sos,5);
R1 =3D sosq(1,:); R2 =3D sosq(2,:); R3 =3D sosq(3,:);
b1 =3D conv(R1(1:3),R2(1:3)); bq =3D conv(R3(1:3),b1);
a1 =3D conv(R1(4:6),R2(4:6)); aq =3D conv(R3(4:6),a1);
[hq,w] =3D freqz(bq,aq,512); gq =3D 20*log10(abs(hq));
plot(w/pi,g,'b',w/pi,gq,'r:');grid
axis([0 1 -70 5]);
xlabel('\omega/\pi');ylabel('Gain, dB');
title('original - solid line, quantized - dashed = line');