Modelagem e Controle de Sistemas

Universidade do Estado do Rio de Janeiro
Instituto Politécnico

Curso de Graduação em Engenharia Mecânica

Modelagem e controle de sistemas

Diego Pablo Vieira de Juarez

Nova Friburgo/RJ,
Dezembro de 2014

ATIVIDADE 1:

Sendo

P1=4
P2=11
P3=3 e As entradas X e X1:

Codigo no Matlab:

clc clear P1 = 4
P2 = 11
P3 = 3 syms s t
G1 = P1/(P2*s+1)
X = P3/s
Y = G1*X y = ilaplace(Y,s,t) t=[0:P2*5] ezplot(y,t) xlabel('TEMPO (t)') ylabel('SAIDA') title('GRAFICOS') datacursormode on

clc clear P1 = 4
P2 = 11
P3 = 3 syms s t
G1 = P1/(P2*s+1)
X = P3/s
Y = G1*X y = ilaplace(Y,s,t) t=[0:P2*5] ezplot(y,t) xlabel('TEMPO (t)') ylabel('SAIDA') title('GRAFICOS') datacursormode on

clc clear P1 = 4
P2 = 11
P3 = 3 syms s t
G1 = P1/(P2*s+1)
X = P3/s
Y = G1*X y = ilaplace(Y,s,t) X1 = 0.5*P3/s
Y1 = G1*X1 y1 = ilaplace(Y1,s,t) t=[0:.25:P2*5] ezplot(y1,t) set(ezplot(y1,t),'Color','green') hold on ezplot(y,t) legend('X1(s)','X(s)') xlabel('TEMPO') ylabel('SAIDA') title('GRAFICOS') datacursormode on

Atividade 2:
Sendo

P4=5
P5=7
P6=6
P7=8
e

Logo:

Codigo no Matlab:

G2 = tf([1 12 35],[1 14 48]) rlocus(G2) E sendo

P8=8
P9=-9
e

Logo:

Codigo no Matlab:

G3 = tf([1 8],[ 1 -9 0]) rlocus(G3) Atividade 3:
Sendo
P10:25
P11:5
P12:0
P13:1
e

Logo

Codigo no matlab:

G4 = tf([25 5],[ 1 1 0]) bode(G4) Conclusão:

Conclui-se após os experimentos realizados a grande eficiência do programa Matlab em representar graficamente as menores perturbações em um sistema de equações. Assim como na primeira atividade onde a comparação é inevitável entre as diferentes entradas X e X1, as quais diferem inicialmente em pequena escala, mas através dos gráficos, nota-se a grande discrepância entre as duas entradas.
Já quanto as atividades 2 e 3 a precisão