Informe 4 Laureano Mena Jhon
Informe 4 Laureano Mena Jhon
Informe 4 Laureano Mena Jhon
TEMA
CURSO
LABORATORIO DE CONTROL 1
DOCENTE
ALUMNO
CUI
20192206
2
Experiencia 4
Análisis en el Dominio del Tiempo
4.1. Objetivos. -
La ecuación anterior plantea que la salida c(t) es inicialmente cero y al final se vuelve
unitaria. Una característica importante de tal curva de respuesta exponencial c(t) es que,
para t=T, el valor de c(t) es 0.632, o que la respuesta c(t) alcanzo 63.2% de su cambio
total. Esto se aprecia con facilidad sustituyendo t=T en c(t). Es decir,
Si 0< <1, los polos en lazo cerrado son complejos conjugados y se encuentran en el
semiplano izquierdo del plano s. El sistema, entonces se denomina subamortiguado y la
respuesta transitoria es oscilatoria.
Ahora se obtendrá la respuesta transitoria para tres casos diferentes a una entrada
escalón unitaria.
a). Caso subamortiguado (0 < <1): en este caso, C(s)/ R(s) se escribe como
b). Caso críticamente amortiguado (=1): si los dos polos son casi iguales, el sistema
reaproxima mediante uno críticamente amortiguado. Los polos se encuentran ubicados
en:
c). Caso sobreamortiguado (>1): en este caso, los dos polos son reales, negativos y
diferentes.
1. Tiempo de retardo, td
2. Tiempo de levantamiento crecimiento, tr
3. Tiempo pico tp
4. Sobrepaso máximo, Mp
5. Tiempo de asentamiento, ts
2 1
2 s +1 s+ 3
9
4.4.1. Utilice Matlab para obtener la salida del sistema ante una entrada paso y una
entrada rampa del ejercicio 4.3.1, compare el resultado con el obtenido
analíticamente.
Entrada Rampa
11
4.4.2. A partir del ejercicio 4.3.1 genere las familias de curvas para los siguientes casos,
compare el tiempo de establecimiento, y el error en estado estable:
a) R1=R2=2kΩ C= 7μF
b) R1=1KΩ R2=2kΩ C=7μF
c) R1=2KΩ R2=1kΩ C=7μF
d) R1=1KΩ R2=2kΩ C=70μF
R1=8*10^3
R2=4*10^3
C=4*10^-6
num=[-R2];
den=[C*R2*R1 R1];
G1=tf(num,den)
t=0:0.00001:0.007;
hold on
y1=step(G1,t);
plot(t,y1,'b')
R11=1*10^3
R22=2*10^3
C11=7*10^-6
num=[-R2];
den=[C11*R22*R11 R11];
G2=tf(num,den)
12
%t=0:0.00001:0.007;
%hold on
y2=step(G2,t);
plot(t,y2,'r')
R111=1*10^3
R222=2*10^3
C111=7*10^-6
num=[-R2];
den=[C111*R222*R111 R111];
G3=tf(num,den)
%t=0:0.00001:0.007;
%hold on
y3=step(G3,t);
plot(t,y3,'g')
R1111=1*10^3
R2222=2*10^3
C1111=7*10^-6
num=[-R2222];
den=[C1111*R2222*R1111 R1111];
G4=tf(num,den)
%t=0:0.00001:0.006;
%hold on
y4=step(G4,t);
plot(t,y4,'m')
hold off
13
14
15
4.4.4. Mediante Matlab, obtenga los siguientes parámetros del ejercicio 4.3.3
a) Tiempo de retardo, td
b) Tiempo de levantamiento crecimiento, tr
c) Tiempo pico, tp
d) Sobrepaso máximo, Mp
e) Tiempo de asentamiento, ts
f) Error en estado estable.
>> Lab_4_5_2
cual es el numerador =2
a=
b=
cual es el coeficiente de s =5
c=
d=
18
sys =
2
---------------
2 s^2 + 5 s + 5
sys =
wn =
1.5811
E=
0.7906
wd =
0.9682
tp =
3.2446
Mp =
0.0173
ts =
3.2000
19
4.5. Informe. -
4.5.1. Utilizando simulink, dibuje la respuesta a una entrada impulso, escalón unitario,
rampa y una entrada tipo ruido del sistema del literal 4.4.4
20
a) Tiempo de retardo, td
b) Tiempo de levantamiento crecimiento, tr
c) Tiempo pico, tp
d) Sobrepaso máximo, Mp
e) Tiempo de asentamiento, ts
f) Error en estado estable.
21
>> Lab_4_5_2
cual es el numerador =4
a=
b=
cual es el coeficiente de s =7
c=
d=
sys =
4
---------------
5 s^2 + 7 s + 8
22
sys =
wn =
1.2649
E=
0.5534
wd =
1.0536
tp =
2.9819
Mp =
0.1240
ts =
5.7143
>> Lab_4_5_3
cual es el numerador =1
a=
b=
cual es el coeficiente de s =3
c=
d=
sys =
24
1
---------------
2 s^2 + 3 s + 4
sys =
wn =
1.4142
E=
0.5303
wd =
1.1990
tp =
2.6203
Mp =
0.1401
ts =
5.3333
25
4.6. Bibliografía. -
Manual de Matlab.