clear

clc

num=[2 6 8]

den=conv([0 1 2],[1 0 4])

[r,p,k]=residue(num,den)

M=abs(r(1))

fi=angle(r(1))*180/pi

t=[0:0.01:10];

w=2

y=r(1)*exp(p(1)*t)+2*M*cos(w*t+fi*pi/180)

plot(t,y)

grid on

clear

clc

A=[-2 -5 -6;1 -1 1;0 1 0]

B=[1;-1;1]

C=[1 3 6]

D=0

[num,den]=ss2tf(A,B,C,D)

y=tf(num,den)