Polytech Projets Ge » Polytech Ge Sous Traitance

h1. Scilab Bilineaire

h1. Bilinéaire inverse

<pre>

function Sysc=dls2cls(Sysd,T)

    s=poly(0,'s');

    z=poly(0,'z');

    numdeg = degree(Sysd.num);

    dendeg = degree(Sysd.den);

    if  numdeg>dendeg

        error('The discrete-time system is improper');

    end

    f = -(s*T+2);

    g = (s*T-2);

    bn = coeff(Sysd.num);

    an = coeff(Sysd.den);

    npad = dendeg - numdeg;

    if npad>0

        for i=0:npad

            bn = [bn 0];

        end

    end

    nums = 0;

    dens = 0;

    for i=0:dendeg

        nums = nums + bn(i+1)*(f^i)*(g^(dendeg-i));

        dens = dens + an(i+1)*(f^i)*(g^(dendeg-i));

    end

    Sysc = syslin('c',nums,dens);

endfunction

</pre>