Laboratório de controle 1 - semana 2 by Marcela Couto Barcellos

marcelabarcellos — 2021-05-10 22:21:28 UTC

f=(%s^4 + 5*%s^3 + 6*%s^2 + 9*%s + 30)/(%s^4 + 6*%s^3 + 21*%s^2 + 46*%s+30)
pfss(f)

marcelabarcellos — 2021-05-10 22:39:02 UTC

B=poly([-1 -2],"s","roots");
A=poly([0 -4 -6], "s","roots");
FT=5*B/A

marcelabarcellos — 2021-05-15 16:11:56 UTC

K1=1600
//K2=1000
//K3=1386
MA1=(K1)/(%s*(%s+7)*(%s+11))
//MA2=(K2)/(%s*(%s+7)*(%s+11))
//MA3=(K3)/(%s*(%s+7)*(%s+11))
H=1
MF1=MA1/.H
//MF2=MA2/.H
//MF3=MA3/.H
t=0:0.1:100;
u=ones(t);
y= csim(u,t,MF1);
plot(t,y)
d1=MF1('den')
//d2=MF2('den')
//d3=MF3('den')
roots(d1)
//roots(d2)
//roots(d3)

marcelabarcellos — 2021-05-15 16:16:56 UTC

marcelabarcellos — 2021-05-15 16:23:50 UTC

//K1=1600
K2=1000
//K3=1386
//MA1=(K1)/(%s*(%s+7)*(%s+11))
MA2=(K2)/(%s*(%s+7)*(%s+11))
//MA3=(K3)/(%s*(%s+7)*(%s+11))
H=1
//MF1=MA1/.H
MF2=MA2/.H
//MF3=MA3/.H
t=0:0.1:100;
u=ones(t);
y= csim(u,t,MF2);
plot(t,y)
//d1=MF1('den')
d2=MF2('den')
//d3=MF3('den')
//roots(d1)
roots(d2)
//roots(d3)

marcelabarcellos — 2021-05-15 16:27:25 UTC

//K1=1600
//K2=1000
K3=1386
//MA1=(K1)/(%s*(%s+7)*(%s+11))
//MA2=(K2)/(%s*(%s+7)*(%s+11))
MA3=(K3)/(%s*(%s+7)*(%s+11))
H=1
//MF1=MA1/.H
//MF2=MA2/.H
MF3=MA3/.H
t=0:0.1:100;
u=ones(t);
y= csim(u,t,MF3);
plot(t,y)
//d1=MF1('den')
//d2=MF2('den')
d3=MF3('den')
//roots(d1)
//roots(d2)
roots(d3)

marcelabarcellos — 2021-05-15 16:28:50 UTC

marcelabarcellos — 2021-05-15 16:29:45 UTC

marcelabarcellos — 2021-05-15 20:57:48 UTC

F1=21/%s
F2=1/(%s+2)
MA=F1*F2
H=1
MF=MA/.H

marcelabarcellos — 2021-05-15 20:59:23 UTC

marcelabarcellos — 2021-05-15 21:59:34 UTC

F1=21/%s
F2=1/(%s+2)
MA=F1*F2
H=1
MF=MA/.H
N=MF('num')
OmegaN=sqrt(21)
zeta=2/(2*OmegaN)
pi=%pi
Mp=exp(-pi*zeta/(sqrt(1-(zeta*zeta))))

marcelabarcellos — 2021-05-15 22:00:45 UTC

marcelabarcellos — 2021-05-15 22:03:40 UTC

F1=21/%s
F2=1/(%s+2)
MA=F1*F2
H=1
MF=MA/.H
N=MF('num')
OmegaN=sqrt(21)
zeta=2/(2*OmegaN)
pi=%pi
//Mp=exp(-pi*zeta/(sqrt(1-(zeta*zeta))))
tpico=pi/(OmegaN*sqrt(1-zeta*zeta))

marcelabarcellos — 2021-05-15 22:04:36 UTC

marcelabarcellos — 2021-05-15 22:26:46 UTC

F1=21/%s
F2=1/(%s+2)
MA=F1*F2
H=1
MF=MA/.H
N=MF('num')
OmegaN=sqrt(21)
zeta=2/(2*OmegaN)
pi=%pi
Mp=exp(-pi*zeta/(sqrt(1-(zeta*zeta))))
tpico=pi/(OmegaN*sqrt(1-zeta*zeta))
t=0:0.1:100
u=ones(t)
y= csim(u,t,MF)
xlabel('tempo(s)')
ylabel('y')
plot(t,y)
drawaxis(x=(tpico), y=(1+Mp),dir='r',val= 'Mp',fontsize=1, textcolor=5, sub_int=11) 
drawaxis(x=(tpico), y=(1+Mp),dir='u',val= 'tpico',fontsize=1, textcolor=5, sub_int=11)

marcelabarcellos — 2021-05-15 22:27:09 UTC

marcelabarcellos — 2021-05-16 00:37:47 UTC

OmegaN1=2
Zeta1=0
//OmegaN2=2
//Zeta2=0,3
//OmegaN3=0,5
//Zeta3=1
FT1=(OmegaN1*OmegaN1)/((%s*%s)+(2*Zeta1*OmegaN1*%s)+(OmegaN1*OmegaN1))
//FT2=(OmegaN2*OmegaN2)/((%s*%s)+(2*Zeta2*OmegaN2*%s)+(OmegaN2*OmegaN2))
//FT3=(OmegaN3*OmegaN3)/((%s*%s)+(2*Zeta3*OmegaN3*%s)+(OmegaN3*OmegaN3))

marcelabarcellos — 2021-05-16 00:43:11 UTC

marcelabarcellos — 2021-05-16 00:45:27 UTC

//OmegaN1=2
//Zeta1=0
OmegaN2=2
Zeta2=0.3
//OmegaN3=0.5
//Zeta3=1
//FT1=(OmegaN1*OmegaN1)/((%s*%s)+(2*Zeta1*OmegaN1*%s)+(OmegaN1*OmegaN1))
FT2=(OmegaN2*OmegaN2)/((%s*%s)+(2*Zeta2*OmegaN2*%s)+(OmegaN2*OmegaN2))
//FT3=(OmegaN3*OmegaN3)/((%s*%s)+(2*Zeta3*OmegaN3*%s)+(OmegaN3*OmegaN3))

marcelabarcellos — 2021-05-16 00:46:58 UTC

marcelabarcellos — 2021-05-16 01:21:08 UTC

Para zeta <1 nota se o sistema subamortecido.

marcelabarcellos — 2021-05-16 01:22:19 UTC

//OmegaN1=2
//Zeta1=0
//OmegaN2=2
//Zeta2=0.3
OmegaN3=0.5
Zeta3=1
//FT1=(OmegaN1*OmegaN1)/((%s*%s)+(2*Zeta1*OmegaN1*%s)+(OmegaN1*OmegaN1))
//FT2=(OmegaN2*OmegaN2)/((%s*%s)+(2*Zeta2*OmegaN2*%s)+(OmegaN2*OmegaN2))
FT3=(OmegaN3*OmegaN3)/((%s*%s)+(2*Zeta3*OmegaN3*%s)+(OmegaN3*OmegaN3))

marcelabarcellos — 2021-05-16 01:23:53 UTC

marcelabarcellos — 2021-05-16 01:29:24 UTC

Para zeta =1 nota se o sistema criticamente amortecido.

marcelabarcellos — 2021-05-16 01:31:14 UTC

Para zeta = 0 nota se o sistema não amortecido.

Polos sobre o eixo imaginario.

marcelabarcellos — 2021-05-16 01:55:18 UTC

Polos no Semi Plano esquerdo

marcelabarcellos — 2021-05-16 02:00:48 UTC

Neste caso o grafico deveria tender a zero pois os polos estão sobre o eixo real negativo... alguma coisa deu errado...rsrsrs

marcelabarcellos — 2021-05-16 02:04:30 UTC

Está correto a questão é o sinal de entrada. O bloco é outro para o sinal degrau unitário.

obs

seleneandrade — 2021-05-19 19:39:29 UTC

exercício correto

obs

seleneandrade — 2021-05-19 19:49:26 UTC

exercício correto

obs

seleneandrade — 2021-05-19 19:57:23 UTC

exercício correto

obs

seleneandrade — 2021-05-19 19:58:01 UTC

exercício correto

obs

seleneandrade — 2021-05-19 19:58:39 UTC

exercício correto

obs

seleneandrade — 2021-05-19 20:09:48 UTC

não foi possível confirmar pelo gráfico apresentado, se os valores calculados foram encontrados graficamente também ( - 0,2 pontos)

obs

seleneandrade — 2021-05-19 20:11:11 UTC

exercício correto

obs

seleneandrade — 2021-05-19 20:11:19 UTC

exercício correto

obs

seleneandrade — 2021-05-19 20:11:51 UTC

exercício correto