Lab3 Control
Lab3 Control
Lab3 Control
➢ Introducción
En este informe de laboratorio, abordaremos cuatro problemas con el
objetivo de utilizar la Transformada de Laplace para resolver sistemas de
ecuaciones diferenciales. Además, nos enfocaremos en conocer los
elementos que conforman el mapa generalizado de un sistema de lazo
abierto y las funciones de cada elemento. También exploraremos cómo
estructurar ambos tipos de sistemas utilizando la representación de
diagramas de bloques.
Además, utilizaremos la herramienta Scilab para graficar el
comportamiento de la solución de algunos problemas. Scilab es una
poderosa herramienta de software libre que nos permitirá visualizar y
analizar los resultados obtenidos a partir de la Transformada de Laplace
y la estructuración de los sistemas de control.
A lo largo de este informe, resolveremos problemas específicos que nos
permitirán aplicar los conceptos teóricos aprendidos y obtener una
comprensión práctica de la Transformada de Laplace y su aplicación en
sistemas de control. Cada problema se abordará paso a paso, desde la
formulación de las ecuaciones diferenciales hasta la obtención de la
solución y la representación gráfica utilizando Scilab.
➢ Desarrollo
Primera Parte
1. Obtener la transforma de Laplace de las siguientes funciones:
función A
1, 0≤𝑡<1
𝑔(𝑡) = {
−1, 2 >𝑡≥1
𝑇 1
∫0 𝑔(𝑡)𝑒 −𝑠𝑡 𝑑𝑡 − 𝑠 𝑒 −2𝑠
𝑙{𝑔(𝑡)} = =
1 − 𝑒 −𝑇𝑠 1 − 𝑒 −2𝑠
Grafica de la función con Scilab.
x = [0:0.001:5];
y = trozos(modulo(x,2));
title("Grafica de la funcion periodica g(t)")
xlabel("Tiempo t");
ylabel("Respuesta g(t)")
plot(x,y);
Función B
2𝑡, 0 ≤ 𝑡 < 0.5
𝑔(𝑡) = {
2𝑡 − 2, 0.5 > 𝑡 ≥ 1
𝑇 1 −0.5𝑠 2 −0.5𝑠 2 2 −𝑠
∫0 𝑔(𝑡)𝑒 −𝑠𝑡 𝑑𝑡 𝑒 + 𝑒 + 2 − 2𝑒
𝑠 𝑠 𝑠 𝑠
𝑙{𝑔(𝑡)} = =
1 − 𝑒 −𝑇𝑠 1 − 𝑒 −𝑠
Grafica de la función con Scilab.
endfunction
x = [0:0.001:5];
y = trozos(x);
title("Grafica de la funcion periodica g(t)")
xlabel("Tiempo t");
ylabel("Respuesta g(t)")
plot(x,y);
Segunda Parte
. Para las siguientes funciones, determinar la respuesta temporal para una
entrada tipo escalón unitario:
s+4
A. 𝐺(𝑠) = (s+2)(s2 +3s+1)
s+8
B. 𝐺(𝑠) = 2 2
𝑠 (s +5s+1)
s
C. 𝐺(𝑠) = (s+7)(2s2 +5s+2)
Respuesta temporal B
s+8 𝐴 𝐵 𝐶 𝐷 𝐹
A. 𝐺(𝑠) = = 𝑠 + 𝑠2 + 𝑠3 + 𝑠+0.2 + 𝑠+4.79
𝑠 2 (s2 +5s+1)
Función de transferencia B
k = lim 𝐺(𝑆) = 0
𝑠→0
Función de transferencia C
k = lim 𝐺(𝑆) = 0
𝑠→0
2.5 Graficar para una entrada tipo rampa e impulso, utilizando la respuesta en
el dominio de la frecuencia.
Graficas para entrada impulso
función de transferencia A
s = %s;
function u=rampa(t)
u = t;
endfunction
num = s+4;
den = (s+2)*(s^2 +3*s +1);
TF = syslin("c", num, den);
t = linspace(0,30,500);
yi =csim("impulse", t, TF);
yr = csim(rampa, t, TF);
𝑦̇ + 𝑢𝑦 + 𝑦 = 𝑢2 + 2
𝑦̇ + 𝑢𝑦 + 𝑦 − 𝑢2 − 2 = 0
13
Al linealizar se obtiene 𝑦̇ + 4𝑦 = 𝑢(𝑠)
4
13 13
Aplicando Laplace 𝑠𝑌(𝑠) + 4𝑌(𝑠) = 𝑢(𝑠) ; (𝑠 + 4)𝑌(𝑠) = 𝑢(𝑠)
4 4
13
0.2
Escalón 𝑢 = 𝑠
t=0:0.1:15;
y=ode(y0,0,t,f);
plot(t,y,'g')
function dy=f(t, y);
dy=-(4*y)+((13/4)*u);
endfunction
y0=1.5;
u0=2;
u=3;
y=15;
t=0:0.1:15;
y=ode(y0,0,t,f);
plot(t,y,'b')
Ilustración 10.Grafica con el comportamiento del modelo realizado. Cuando el valor de entrada es 3 (línea azul), y
3.2(línea verde
Cuarta Parte
𝐶(𝑆) 𝐺1 𝐺2
=
𝑅(𝑆) 1 + 𝐻1 𝐺2 + 𝐻1 𝐻2 𝐺1 𝐺2 + 𝐺1 𝐺2
Para encontrar la función de transferencia total del diagrama de bloques b,
primero se obtienen las siguientes relaciones:
A=R(S)-Y
C(S)=A*G1*G2
X=C(S)-Y
Y=X*H1
𝑑𝑤(𝑡) 𝑑𝑥(𝑡)
➢ Para + 10𝑤(𝑡) = + 8𝑥(𝑡)
𝑑𝑦 𝑑𝑡
∆𝑚 − 0.25∆𝑙 − 1.92∆𝑛 = 0
∆𝑚 = 0.25∆𝑙 + 1.92∆𝑛
En el dominio de la función se obtiene:
8
b) 0 = 𝑤0 − 5𝑦0 𝑦0 = 50 𝑥0 con el valor de 𝑤0 = 2.4 entonces
𝑦0 = 0.48
𝑚0
c) 𝑚0 = 𝑙0 𝑛0 𝑙0 = con los valores de 𝑚0 = 0.48 y 𝑛0 =
𝑛0
0.25 entonces 𝑙0 = 1.92
𝑧
d) 20𝑛0 = 𝑧0 𝑛0 = 200 con el valor de 𝑧0 = 5 entonces 𝑛0 =
0.25
e) 𝑚0 = 𝑦0
Por lo tanto, el diagrama de bloques del sistema sería el siguiente: