Digital Processing of Continuous-Time Signals

Digital Processing of Continuous-Time Signals

Program that uses the function zp2sos

% Program

% Conversion of a rational transfer function

% to its factored form

num = input(’Numerator coefficient vector = ’);

den = input(’Denominator coefficient vector = ’);

[A, B] = eqtflength(num, den);

[z,p,k] = tf2zp(A, B);

sos = zp2sos(z,p,k)

Program to develop a cascade realization of the following causal IIR transfer

Function

% Program

% Parallel Form Realizations of an IIR Transfer Function

num = input(’Numerator coefficient vector = ’);

den = input(’Denominator coefficient vector = ’);

[r1,p1,k1] = residuez(num,den);

[r2,p2,k2] = residue(num,den);

disp(’Parallel Form I’)

disp(’Residues are’);

disp(r1);

disp(’Poles are at’);

disp(p1);

disp(’Constant value’);

disp(k1);

disp(’Parallel Form II’)

disp(’Residues are’);

disp(r2);

disp(’Poles are at’);

disp(p2);

disp(’Constant value’); disp(k2);