Lista 3 de exercicios- Thátilla by Thátilla Curty

scilab

thatillacurty — 2021-09-28 20:47:31 UTC

//exercicio 1
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)

observamos que é um sistema estável, uma vez que nao há variação do sinal na primeira coluna.

Resultado

thatillacurty — 2021-09-28 20:53:38 UTC

scilab

thatillacurty — 2021-09-28 20:56:54 UTC

//questao 2

clc;
clear all
s = %s;
k = 1600;  // modificar os valores      
MA = syslin ('c', (k)/(s*(s+7)*(s+11))); 
H=1;
MF = MA/.H;

t=0.1:0.1:100;
u = ones(t)

y=csim(u,t,MF);
plot(t,y)

letra A k=1600

thatillacurty — 2021-09-28 20:59:06 UTC

Sistema Instável
Polos: -18.509986 + 0.i
0.2549929 + 9.2938048i
0.2549929 - 9.2938048i

letra B k=1000

thatillacurty — 2021-09-28 21:03:28 UTC

thatillacurty — 2021-09-28 21:06:10 UTC

Sistema Estável
Polos: -16.939395 + 0.i
-0.5303027 + 7.6650353i
-0.5303027 - 7.6650353i

letra C k=1386

thatillacurty — 2021-09-28 21:06:36 UTC

thatillacurty — 2021-09-28 21:07:57 UTC

Sistema Marginalmente Instável
Polos: -18. + 0.i
-5.551D-17 + 8.7749644i
-5.551D-17 - 8.7749644i

a) scilab

thatillacurty — 2021-09-28 21:08:39 UTC

//questao numero 3 - lista 3

clc
clear all
s=%s; //declara a variavel complexa s 

w=2;
l=0;
MA = syslin('c',(w^2 )/(s^2+2*l*w*s +w^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)

disp(polosT)
disp(zerosT)

a) grafico

thatillacurty — 2021-09-28 21:09:06 UTC

0. + 2.82842712474619i
0. - 2.82842712474619i

b)scilab

thatillacurty — 2021-09-28 21:09:36 UTC

//questao numero 3 - lista 3

clc
clear all
s=%s; //declara a variavel complexa s 

w=2;
l=0,1;
MA = syslin('c',(w^2 )/(s^2+2*l*w*s +w^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)

disp(polosT)
disp(zerosT)

b) gráfico

thatillacurty — 2021-09-28 21:10:10 UTC

-0.2 + 2.8213472i
-0.2 - 2.8213472i

c)scilab

thatillacurty — 2021-09-28 21:10:32 UTC

//questao numero 3 - lista 3

clc
clear all
s=%s; //declara a variavel complexa s 

w=0,5;
l=0,9;
MA = syslin('c',(w^2)/(s^2+2*l*w*s +w^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)

disp(polosT)
disp(zerosT)

c)scilab

thatillacurty — 2021-09-28 21:11:19 UTC

-0.45 + 0.5454356i
-0.45 - 0.5454356i

thatillacurty — 2021-09-28 21:12:15 UTC

Na letra A, por termos polos no eixo puramente imaginário, obtemos um gráfico semelhante ao destacado de vermelho.
E na letra B e C, por termos parte real e imaginária, o gráfico foi oscilando ate obter estabilidade num certo valor, como o gráfico destacado de azul.

scilab

thatillacurty — 2021-09-28 21:14:15 UTC

// questao 4

C = 21/(%s)
P = 1/(%s +2)
MA = C*P
H = 1

MF = MA/.H
t = 0.1:0.1:20
u = ones(t)
y = csim(u,t,MF)
plot(t,y)

letra a)

thatillacurty — 2021-09-28 21:14:49 UTC

21
----------
21 +2s +s²

letra b)

thatillacurty — 2021-09-28 21:15:35 UTC

letra c)

thatillacurty — 2021-09-28 21:23:59 UTC

scilab

thatillacurty — 2021-09-28 21:25:39 UTC

//questao 5

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("Equação em malha fechada -- K=1")
disp(MF)
disp("Polo 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("Equação em malha fechada--K=2")
disp(MF)
disp("Polo 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("Equação em malha fechada--K=5")
disp(MF)
disp("Polo da função")
disp(roots(MF.den))

resultado

thatillacurty — 2021-09-28 21:28:37 UTC

O sistema é estável para k=1, como podemos conferir pelas raízes da equação característica.
Apenas as raizes de k=1 estao do lado esquerdo do plano.