# data[].attributes.body
clc;
clear all

s=%s;

k=1;

MA = syslin('c', k*(s^2-s+2)/(s^2+2*s+1));

H = 1

MF = MA/.H

disp("Eq em malha fechada (K=1)")
disp(MF)
disp("Polos da função")
disp(roots(MF.den))

k=2;

MA =  k*(s^2-s+2)/(s^2+2*s+1);

H = 1

MF = MA/.H

disp("Eq em malha fechada (K=2)")
disp(MF)
disp("Polos da função")
disp(roots(MF.den))

k = 5;

MA =  k*(s^2-s+2)/(s^2+2*s+1);

H = 1

MF = MA/.H

disp("Eq em malha fechada (K=5)")
disp(MF)
disp("Polos da função")
disp(roots(MF.den))

# data[].attributes.body
clc
clear all

s=%s;//Declara a variável complexa s

C=syslin('c',(1)/(s+2));
G=syslin('c',(1)/((s^2)+(2*s)+4));
H=1;
MF=syslin('c',(C*G)/.H)
n=MF.num
ztf=roots(n);
d=MF.den
ptf=roots(d)
disp("Zeros:")
disp(ztf)
disp("Polos:")
disp(ptf)
disp("Sistema em malha fechada")
disp(MF)
k=0;
r = routh_t(MF,k)
disp("Criterio de Routh-Hurwitz:")
disp(r)
disp("")
disp("Desta forma podemos entender como um sistema estável, tendo em vista que temos 3 polos no spe e nao ha variacao do sinal na primeira coluna")

# data[].attributes.body
clc;clear all
s=%s;

C = 21/(s)

P = 1/(s+2)

MA = C*P

H = 1

MF = MA/.H

disp(MF)

t = 0.1:0.1:20

u = ones(t)

y = csim(u,t,MF)

plot(t,y)

# data[].attributes.body
clc
clear all
s=%s;

wn=2;
zeta=0;
MA = syslin('c',(wn^2 )/(s^2+2*zeta*wn*s +wn^2))
H=1;

MF=syslin('c',MA/.H)

d = MF.den
polosT = roots ( d );
n = MF.num
zerosT = roots ( n );

t=0.1:0.1:100;
y=csim('impuls',t,MF);
plot(t,y)

# data[].attributes.body
Resposta oscilatório, o que nos permite concluir que os polos estão sobre o eixo imaginário.

# data[].attributes.body
Com base no que foi feito anteriormente, é possível concluir que o sistema é estável apenas para K = 1. Pois somente os polos de K= 1 estão do lado esquerdo do plano

# data[].attributes.body
clc
clear all
s=%s;

wn=2;
zeta=0.1;
MA = syslin('c',(wn^2 )/(s^2+2*zeta*wn*s +wn^2))
H=1;

MF=syslin('c',MA/.H)

d = MF.den
polosT = roots ( d );
n = MF.num
zerosT = roots ( n );

t=0.1:0.1:100;
y=csim('impuls',t,MF);
plot(t,y)

# data[].attributes.body
Para as demais questões, alteramos apenas o numerador

# 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
Resposta Subamortecida

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

# data[].attributes.body
clc
clear all
s=%s;

wn=0.5;
zeta=0.9;
MA = syslin('c',(wn^2 )/(s^2+2*zeta*wn*s +wn^2))
H=1;

MF=syslin('c',MA/.H)

d = MF.den
polosT = roots ( d );
n = MF.num
zerosT = roots ( n );

t=0.1:0.1:100;
y=csim('impuls',t,MF);
plot(t,y)

# data[].attributes.body
Resposta subamortecida

# data[].attributes.body
Por configuração do Scilab os gráficos foram sobrepostos, tentei alterar, porém não consegui