bilt — bilinear or biquadratic transform SISO system given by a zero/poles representation
[npl,nzr,ngn] = bilt(pl,zr,gn,num,den)
a vector, the poles of the given system.
a vector, the zeros of the given system.
a polynomial with degree equal to the degree of
den, the numerator of the
transform.
a polynomial with degree 1 or 2, the denominator of the transform.
a vector, the poles of the transformed system.
a vector, the zeros of the transformed system.
a scalar, the gain of the transformed system.
function for calculating the gain poles and zeros which result from a bilinear transform or from a biquadratic transform. Used by the functions iir and trans.
Hlp=iir(3,'lp','ellip',[0.1 0],[.08 .03]);
pl=roots(Hlp.den);
zr=roots(Hlp.num);
gn=coeff(Hlp.num,degree(Hlp.num))/coeff(Hlp.den,degree(Hlp.den));
z=poly(0,'z');
a=0.3;
num=z-a;
den=1-a*z;
[npl,nzr,ngn] = bilt(pl,zr,gn,num,den)
Hlpt=ngn*poly(nzr,'z','r')/poly(npl,'z','r')
//comparison with horner
horner(Hlp,num/den)