TALLER N.

º 4
SOLUCIÓN DE UN SISTEMA DE ECUACIONES DIFERENCIALES.
Matemáticas Aplicadas a la ingeniería química
Nombre: Ronni Bermúdez Espinosa
Profesor Francisco Muñoz Paba M. Sc
Fecha de entrega: octubre 25 del 2021
SOLUCION
Punto (a)
Para desarrollar las ecuaciones que modelan ambos tanques, hay que realizar un balance sobre el
sistema. En este caso cada tanque representa 1 sistema. La forma general del balance es:
Acumulación = Entrada −Salida
Para el tanque 1 se obtiene:
𝑑𝑉1
= 𝐹0 − 𝐹1 (1)
𝑑𝑡
Para el tanque 2 se obtiene:
𝑑𝑉2
= 𝐹1 − 𝐹2 (2)
𝑑𝑡
Ahora, teniendo en cuenta que en forma general el volumen es igual al producto del área de la base
con la altura (V = A∗h), se puede reemplazar para cada uno de los tanques:
𝑑(𝐴1 ℎ1 )
= 𝐹0 − 𝐹1 (3)
𝑑𝑡
𝑑(𝐴2 ℎ2 )
= 𝐹1 − 𝐹2 (4)
𝑑𝑡
Como se especifica el área constante, esta puede salir del diferencial de la ecuación, tanto para el
tanque 1 como el tanque 2. Como el flujo de salida en cada tanque tiene una relación lineal con la
altura del mismo, significa que tiene la forma:
𝐹𝑖= 𝑎𝑖 − ℎ𝑖 (5)
Donde α representa una constante de proporcionalidad. Con esto, se puede obtener la ecuación para
cada uno de los tanques:

𝑑ℎ2
𝐴2 = 𝑎ℎ1 − 𝛽ℎ2 (6)
𝑑𝑡
Reorganizando los términos de las ecuaciones se llega a:
𝑑ℎ1 𝐹0 − 𝑎ℎ1
= (7)
𝑑𝑡 𝐴1
𝑑ℎ2 𝑎ℎ1 − 𝛽ℎ2
= (8)
𝑑𝑡 𝐴2

Dado el pasado sistema de ecuaciones diferenciales, de la ecuación 8 despejamos h1:

𝑎ℎ1 − 𝛽ℎ2
ℎ2′ = (9)
𝐴2
 𝑎ℎ1 𝛽ℎ2
ℎ2′ = − (10)
𝐴2 𝐴2
𝑎ℎ1 𝛽ℎ2
= ℎ2′ + (11)
𝐴2 𝐴2
𝐴2 ′ 𝛽ℎ2
ℎ1 = ℎ2 + (12)
𝑎 𝑎
Derivamos h1
𝐴2 ′′ 𝛽 ′
ℎ1′ = ℎ2 + ℎ2 (13)
𝑎 𝑎
Reemplazamos h1 y h1’ en la ecuación (7).
𝐴2 ′ 𝛽ℎ2
𝐴2 ′′ 𝛽 ′ 𝐹0 − 𝑎( 𝑎 ℎ2 + 𝑎 )
ℎ2 + ℎ2 = (14)
𝑎 𝑎 𝐴1
𝐴2 ′′ 𝛽 ′ 𝐹𝑜 − 𝐴2ℎ2′ − 𝛽ℎ2
ℎ2 + ℎ2 = (15)
𝑎 𝑎 𝐴1
Obtenemos la ecuación de segundo orden en función de h2:
𝐴2 ′′ 𝛽 ′ 𝐴2ℎ2′ + 𝛽ℎ2 𝐹𝑜
ℎ2 + ℎ2 + = (16)
𝑎 𝑎 𝐴1 𝐴1

Punto (b)
Para resolver el sistema de ecuaciones desarrollado anteriormente para los tanques 1 y 2, se utiliza
la transformada de Laplace. Entonces, tomando la ecuación para el tanque 1:
𝑑ℎ1 𝐹0 − 𝑎ℎ1
= (17)
𝑑𝑡 𝐴1
Reemplazando los valores correspondientes se llega a:

𝑑ℎ1
= 0.2 − 0.1 ∗ ℎ1 (18)
𝑑𝑡
Si se aplica transformada de Laplace a ambos lados de la ecuación se obtiene:
𝑑ℎ1
𝐿( ) = 𝐿(0.2) − 0.1𝐿 (ℎ1) (19)
𝑑𝑡
0.2
𝑆𝐻1(𝑆) − ℎ1(0) = − 0.1𝐻1(𝑆) (20)
𝑆
Reemplazado el valor inicial de h1 (h1(0) = 1), se obtiene:

0.2
𝑆𝐻1(𝑆) − 1 = − 0.1𝐻1(𝑆) (21)
𝑆
Como en esta función se encuentra solamente H1(S), se puede resolver de la siguiente manera:
0.2
𝑆𝐻1(𝑆) + 0.1𝐻1(𝑆) = + 1 (22)
𝑆
0.2 + 𝑆
𝐻1(𝑆)(𝑆 + 0.1) = (23)
𝑆
0.2 + 𝑆 1
𝐻1(𝑆) = ( )( ) (24)
𝑆 𝑆 + 0.1
0.2 + 𝑆
𝐻1(𝑆) = ( ) (25)
𝑆(𝑆 + 0.1)
Esta última ecuación obtenida, corresponde a la transformada de Laplace de la ecuación que
describe la altura del tanque 1 con respecto al tiempo. Esta ecuación se puede separar en fracciones
parciales, para esto se aplica el método de los límites.
𝐴 𝐵
𝐻1(𝑆) = + (26)
𝑆 𝑆 + 0.1
0.2 + 𝑆
𝐴 = lim =2 (27)
𝑆→0 𝑆 + 0.1

0.2 + 𝑆
𝐵 = lim = −1 (28)
𝑆→−0.1 𝑆
Ahora se reemplazan los valores de A y B:
2 1
𝐻1(𝑆) = − (29)
𝑆 𝑆 + 0.1
A esta ecuación se le puede aplicar la transformada inversa de Laplace para llevarla nuevamente al
dominio del tiempo, se obtiene:
2 1
𝐿−1 (𝐻1(𝑆)) = 𝐿−1 ( ) − 𝐿−1 ( ) (30)
𝑆 𝑆 + 0.1
ℎ1 (𝑡) = 2 − 𝑒 −0.1𝑡 (31)
Ahora, para obtener la ecuación que describe la altura del tanque 2 con respecto al tiempo, se sigue
un procedimiento similar con la ecuación diferencial del tanque 2:
𝑑ℎ2 𝑎ℎ1 𝛽ℎ2
= − (31)
𝑑𝑡 𝐴2 𝐴2
Reemplazando los valores correspondientes:

𝑑ℎ2
= 0.05ℎ1 − 0.05ℎ2 (32)
𝑑𝑡
Ahora se aplica transformada de Laplace, con la condición inicial h2(0) = 0.
𝑑ℎ2
𝐿( ) = 𝐿 (0.05ℎ1) − 𝐿(0.05ℎ2) (33)
𝑑𝑡
 𝑆𝐻2(𝑆) − ℎ2(0) = 0.05𝐻1(𝑆) − 0.05𝐻2(𝑆) (34)
𝑆𝐻2(𝑆) = 0.05𝐻1(𝑆) − 0.05𝐻2(𝑆) (35)
0.05𝐻1 (𝑆)
𝐻2(𝑆) = (36)
𝑆 + 0.05
Ahora se reemplaza H1(S) por su correspondiente función:
0.05 0.2 + 𝑆
𝐻2(𝑆) = ( )( ) (37)
𝑆 + 0.05 𝑆(𝑆 + 0.1)
0.01 + 0.05(𝑆)
𝐻2(𝑆) = ( ) (38)
𝑆(𝑆 + 0.05)(𝑆 + 0.1)
Esta ecuación se puede separar en fracciones parciales:

𝐴 𝐵 𝐶
𝐻2(𝑆) = + + (39)
𝑆 𝑆 + 0.05 𝑆 + 0.1
Se hallan los valores de A, B y C por medio de límites:

0.01 + 0.05𝑆
𝐴 = lim =2 (40)
𝑆→0 (𝑆 + 0.05)(𝑆 + 0.1)

0.01 + 0.05𝑆
𝐵= lim = −3 (41)
𝑆→−0.05 𝑆(𝑆 + 0.1)

0.01 + 0.05𝑆
𝐶 = lim =1 (42)
𝑆→−0.1 𝑆(𝑆 + 0.05)

Reemplazando los valores para la ecuación:

2 3 1
𝐻2(𝑆) = − + (43)
𝑆 𝑆 + 0.05 𝑆 + 0.1
Se aplica transformada inversa y se obtiene:
2 3 1
𝐿−1 (𝐻2(𝑆)) = 𝐿−1 ( ) − 𝐿−1 ( ) + 𝐿−1 ( ) (44)
𝑆 𝑆 + 0.05 𝑆 + 0.1
ℎ2 (𝑡) = 2 − 3𝑒 −0.05𝑡 + 𝑒 −0.1𝑡 (45)
Procedemos a realizar la gráfica en Matlab:
 En el editor de Matlab:
t= [0 5];
h1=(2 - exp(-0.1.*t));
h2=(2-3.*exp(-0.05.*t)+ exp(-0.1.*t));
plot(t,h1,'b:+',t,h2,'b:*')
xlabel('t')
ylabel('h')
legend('h1','h2')
disp([t;h1;h2])

 Command Window:
>>Punto2
0 5.0000
1.0000 1.3935
0 0.2701

El gráfico obtenido en Matlab es el siguiente:

Punto (c)
Ahora resolvemos el sistema de ecuaciones diferenciales:
Primeramente, realizamos un cambio de variables, donde decimos que z1=h1 y z2=h2.
𝑑ℎ1 𝐹0 − 𝑎𝑧1 2 − 1𝑧1
= = (46)
𝑑𝑡 𝐴1 10
𝑑ℎ2 𝑎𝑧1 − 𝛽𝑧2 1𝑧1 − 1𝑧2
= = (47)
𝑑𝑡 𝐴2 20
Aplicamos la función teniendo en cuenta el cambio de variable en Matlab y la guardamos.

 En el editor de Matlab:
function xyd=pd(t,z)
xyd=[0.2-0.1.*z(1); 0.05*z(1)-0.05.*z(2)];
Ahora en otro editor realizamos la operación con ode45.

 En el editor de Matlab:
[t,z]=ode45('pd', [0:0.5:5], [1 0]);
plot(t,z(:,1), '-b',t,z (:,2), 'b')
disp([t,z])

 Command Window:
>>Punto3C

0 1.0000 0
0.5000 1.0488 0.0253
1.0000 1.0952 0.0511
1.5000 1.1393 0.0775
2.0000 1.1813 0.1042
2.5000 1.2212 0.1313
3.0000 1.2592 0.1587
3.5000 1.2953 0.1863
4.0000 1.3297 0.2141
4.5000 1.3624 0.2421
5.0000 1.3935 0.2701
El gráfico obtenido en Matlab es el siguiente:

Como nos damos cuenta las gráficas son similares, pero en los resultados se presenta una diferencia
en los decimales de las alturas de los tanques, esto se debe a que en ocasiones no se tienen en cuenta
todas las cifras significativas necesarias.