# data[].attributes.body
clc;clear all
s = %s;
valores = [1600,1000,1386];
for i = 1:3
k = valores(i);
MA = syslin ('c', (k)/(s*(s+7)*(s+11)));
H=1;
MF = MA/.H;
d = MF.den ;
polos = roots(d);
mprintf("\n Polos do sistema para K = %d:",valores(i))
disp(polos);
end

# data[].attributes.headline
FT Do sistema em malha fechada e Critério de Routh Hurwitz

# data[].attributes.subject
FT Do sistema em malha fechada e Critério de Routh Hurwitz

# data[].attributes.body
Código similar foi utilizado durante todo o desenolvimento, alterando apenas os valores de w e z.

# data[].attributes.body
Sistema Instável
Polos: -18.509986 + 0.i
0.2549929 + 9.2938048i
0.2549929 - 9.2938048i

# data[].attributes.body
Sistema Estável
Polos: -16.939395 + 0.i
-0.5303027 + 7.6650353i
-0.5303027 - 7.6650353i

# data[].attributes.body
Sistema Marginalmente Instável
Polos: -18.       + 0.i
-5.551D-17 + 8.7749644i
-5.551D-17 - 8.7749644i

# data[].attributes.body
Neste caso temos uma resposta oscilatório, podemos inferir que os polos deste sistema estão sobre o eixo imaginário.

# data[].attributes.body
Temos neste caso uma resposta subamortecido, em razão dos polos estarem no semiplano esquerdo.

# data[].attributes.body
Novamente temos uma resposta subamortecida, em razão da localização dos polos no semiplano esquerdo.

# data[].attributes.body
Portanto, o sistema em malha fechada será estável para K = 1, devido a posição dos polos no semiplano esquerdo.

Para K = 2 Temos um sistema marginalmente estável;
Para K = 5 Temos um sistema instável.