Pomana

function [] = lab4()
    A = [0 0 1; 0 -1 6; 0 1 0]
    B = [0; 1; 0]
    C = [-2 0 1]
    [n n] = size(A);
    syms s k z h t l;

    Hs = simplify( C * inv(s * eye(n) - A) * B)

    yft = ilaplace(Hs / s)

    tmp = ilaplace(Hs / s, s, k * h);
    Hdz = simplify(((z - 1) / z) * ztrans(tmp, k, z))

    [N, D] = numden(Hdz);
    beta = coeffs(N, z);
    alfa = coeffs(D, z);

    alfa = alfa / alfa(1);
    AC = [ 0 *eye(2) 1*eye(2); -alfa(3) * eye(2) -alfa(2) * eye(2)]
    BC = [0*eye(2) ; 1*eye(2)]
    CC = [ beta(2) * eye(2) beta(1) * eye(2)]
end