program_8_3.m

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

% Program 8_3
% Illustration of Cascade Realization of IIR Filters
% and Structure Verification
%
format long
x =3D [1 zeros(1,6)];
b1 =3D 0.0662272*[1 1];
a1 =3D [1 -0.2593284];
y1 =3D filter(b1,a1,x,0);
si =3D [0 0];
b2 =3D [1 2 1];
a2 =3D [1 -0.6762858 0.3917468];
y2 =3D filter(b2,a2,y1,si);
[p,d] =3D strucver(y2,3);
disp('Actual numerator coefficients are ');disp(p');
disp('Actual denominator coefficients are =');disp(d');