# data[].attributes.body
--> P0= [s^7+s^6+2*s^5+2*s^4-s^3-s^2-2*s-2]
--> roots(P0)

# data[].attributes.body
--> u=ones(t);
--> t=0:0.1:10;
--> plot(t,u)

# data[].attributes.body
--> ma=  1* (%s^2-%s+2)/(%s^2+2*%s+1)
--> h=1;
--> mf=ma/.h;
--> evans(mf)

# data[].attributes.body
Analisando a tabela, levando em consideração que originalmente linha  s5 é uma linha nula, obtém-se o polinômio:

P0= s^7+s^6+2*s^5+2*s^4-s^3-s^2-2*s-2

Na tabela há uma mudança de sinal, tendo assim uma raiz no SPD. Por simetria tem que ter a outra raiz no SPE. Como existem 6 raízes as outras 4 raízes são sobre o eixo
imaginário.

# data[].attributes.body
--> ma= 3.47/(%s^2 + 1.33*%s + 3.47)
-->  h=1
--> mf=ma/.h
--> t=0:0.1:10;
--> u=ones(t);
--> y= csim(u,t,mf);
--> plot(t,y)

# data[].attributes.body
--> g1= 5 / (%s + 5)
--> t=0:0.1:10;
--> u=ones(t);
--> y= csim(u,t,g1);
--> plot(t,y)

# data[].attributes.body
--> g2= 20 / (%s + 20)
--> t=0:0.1:10;
--> u=ones(t);
--> y= csim(u,t,g2);
--> plot(t,y)

# data[].attributes.body
sistema estável

# data[].attributes.body
constante de tempo = 0,20
tempo de subida = 0,44
tempo de assentamento =  0,79

# data[].attributes.body
constante de tempo = 0,073
tempo de subida = 0,12
tempo de assentamento =  0,19

este sistema é mais rápido devido a constante de tempo (A), onde quanto maior o valor, mais rápido a resposta e mais rápido o tempo de assentamento, uma vez que o tempo de assentamento é medido por 4/A.

# data[].attributes.body
--> t= [valores t da tabela];
--> s= [valores s da tabela];
--> plot(t,s);

# data[].attributes.body
--> g= 1654 / (%s^2 + 25.78*%s + 1654)
--> u=ones(t);
--> t=0:0.01:10;
--> y= csim(u,t,g);
--> plot(t,y)

# data[].attributes.body
--> ma=  2* (%s^2-%s+2)/(%s^2+2*%s+1)
--> h=1;
--> mf=ma/.h;
--> evans(mf)

# data[].attributes.body
--> ma= 5* (%s^2-%s+2)/(%s^2+2*%s+1)
--> h=1;
--> mf=ma/.h;
--> evans(mf)

# data[].attributes.body
sistema marginalmente estável

# data[].attributes.body
sistema instável

# data[].attributes.body
--> t0= [valores t da tabela];
--> s= [valores s da tabela];
--> g= 1654 / (%s^2 + 25.78*%s + 1654)
--> u=ones(t);
--> t=0:0.05:10;
--> y1=csim(u,t,g);
--> plot(t,y1,'r')
--> plot(t,u,'g')
--> plot(t0,s,'b')