Unidad 5 Ecuaciones lineales de segundo orden con coecientes variables

5.1 Ecuación de Bessel

Sistemas de ecuaciones lineales de primer orden

Ejemplo Dos grandes tanques, cada uno de los cuales contiene 24 litros de una solución salina, están
conectados entre sí mediante unos tubos, como se muestra en la gura

El tanque A recibe agua pura a razón de 6 litros/minuto y el líquido sale del tanque B con la misma
razón; además, se bombean 8 litros/minuto de líquido del tanque A al tanque B y 2 litros/minuto
del tanque B al tanque A. Los líquidos dentro de cada tanque se mantienen bien revueltos, de modo
que cada mezcla es homogénea. Si en un principio la solución salina en el tanque A contiene x0
kg de sal y la del tanque B contiene inicialmente y0 kg de sal, determinar la masa de sal en cada
tanque en el instante t > 0.
Solución Observe que el volumen de líquido en cada tanque es constante e igual a 24 litros, debido al
equilibrio entre las razones de entrada y salida. Por lo tanto, tenemos dos funciones incógnitas de
t: la masa de sal x(t) en el tanque A y la masa de sal y(t) en el tanque B. Si centramos nuestra
atención en un tanque a la vez, podemos deducir dos ecuaciones que relacionen estas incógnitas.
Como el sistema es alimentado con agua pura, esperamos que el contenido de sal de cada tanque
tienda a cero cuando t → ∞.
Para formular las ecuaciones de este sistema, igualamos la razón de cambio de sal en cada tanque
con larazón neta con la que se transere la sal a ese tanque. La concentración de sal en el tanque
 
x(t) 8x
A es [kg/litro], de modo que el tubo superior saca sal del tanque A a razón de
24   24
2y
[kg/minuto]; de manera similar, el tubo inferior lleva sal al tanque A a razón de [kg/minuto]
hyi 24
(la concentración de sal en el tanque B es [kg/litro)]. El ujo de agua pura, por supuesto no
24
transere sal (simplemente mantiene el volumen del tanque A en 24 litros). Por nuestra premisa,
dx
=razón de entrada-razón de salida,
dt
de modo que la razón de cambio de la masa de sal en el tanque A es
dx 2 8 1 1
= y− x= y− x
dt 24 24 12 3
La razón de cambio de sal en el tanque
  B se determina mediante los mismos tubos de conexión y
6y
por el tubo de drenado, que saca [kg /minuto]:
24

dy 8 2 6 1 1
= x− y− y = x− y
dt 24 24 24 3 3

Facultad de Ciencias UNAM Prof. Esteban Rubén Hurtado Cruz

Ecuaciones Diferenciales 1 1
Unidad 5 Ecuaciones lineales de segundo orden con coecientes variables
5.1 Ecuación de Bessel

Así, los tanques interconectados quedan descritos mediante un sistema de ecuaciones diferenciales:
x0 = − 31 x + 12
1
y
0 1 1 (1)
y = 3x − 3y

Aunque ambas incógnitas x(t) y y(t) aparecen en cada una de las ecuaciones (están acopladas), la
estructura es tan transparente que podemos obtener una ecuación sólo en términos de y, despejando
x en la segunda ecuación,
x = 3y 0 + y (2)
y sustituyendo (1) en la primera ecuación para eliminar x:
1 1
(3y 0 + y)0 = − (3y 0 + y) + y
3 12
00 0 0 1 1
3y + y = −y − y + y
3 12
o
1
3y 00 + 2y 0 + y = 0
4
Esta última ecuación, que es lineal con coecientes constantes, se resuelve fácilmente mediante la
ecuación auxiliar
1
3r2 + 2r + =0
4
1 1
que tiene raíces − , − y cuya solución genral es
2 6

(3)
t t
y(t) = c1 e− 2 + c2 e− 6

Una vez determinada y, la usamos en la ecuación (2) para deducir una fórmula para x:
  1 1
(4)
t t t t t t
x(t) = 3 −c1 e− 2 − c2 e− 6 + c1 e− 2 + c2 e− 6 = − c1 e− 2 + c2 e− 6
2 2
Las fórmulas contienen dos parámetros indeterminados, c1 y c2 , que podemos ajustar para cumplir
con las condiciones iniciales dadas:
1 1
x(0) = − c1 + c2 = x0 , y(0) = c1 + c2 = y0
2 2
o
y0 − 2x0 y0 + 2x0
c1 = , c2 =
2 2
Así, las masas de sal en los tanques A y B en el instante t son, respectivamente,
y0 −2x0 t y0 +2x0 t
e− 2 e− 6



x(t) = − +
4
y0 −2x0 t
4
y0 +2x0 t (5)
e− 2 e− 6



y(t) = − 2 + 2

Facultad de Ciencias UNAM Prof. Esteban Rubén Hurtado Cruz

Ecuaciones Diferenciales 1 2
Unidad 5 Ecuaciones lineales de segundo orden con coecientes variables
5.1 Ecuación de Bessel

Supóngase que una partícula de masa m se mueve en el plano XY debido a que actúa sobre ella una
fuerza tal que su dirección está en el plano XY, y que su magnitud depende de la posición instantenea
(x, y) de la partícula en l tiempo t, si la pertícula está iniciañmente en reposo en algún punto, digamos
el origen, es natural preguntar donde estará la particula en cualquier tiempo posterior. Si consideramos
la fuerza en dos componentes F1 y F2 en las direcciones positivas X y Y como se ve en la gura

entonces por la ley de Newton se tiene que:

d2 d2
m x = F1 (x, y, t), m y = F2 (x, y, t) (6)
dt dt
con las condiciones iniciales x(0) = y(0) = x0 (0) = y 0 (0) = 0. A la expresión (6) se le denomina Sistema
de ecuaciones diferenciales, siendo su forma más general la expresión
Fm (t, f1 , f 0 (t), ...f n , ..., fm , fm
0 n
, ..., fm ) (7)
El orden de la mayor derivada que aparece en el sistema de ecuaciones diferenciales dene el orden del
sistema, por ejemplo el sistema (6) es de segundo orden.

ejemplo Considérese el sistema resorte-masa de la gura. Las dos masas se mueven sobre un supercie
sin fricción bajo la inuencia de las fuerzas externas F1 (t) y F2 (t) y también están restringidas por
los tres resortes cuyas constantes son k1 , k2 y k3 , respectivamente.

Se encuentran las siguientes ecuaciones para las coordenadas x1 y x2 de las dos masas
d2
m1 x1 = k2 (x2 − x1 ) − k1 x1 + F1 (t)
dt2
= −(k1 + k2 )x1 + k2 x2 + F1 (t)

d2
m2 x2 = −k3 x2 − k2 (x2 − x1 ) − k1 x1 + F2 (t)
dt2
= k2 x1 − (k1 + k2 )x2 + F2 (t)

Facultad de Ciencias UNAM Prof. Esteban Rubén Hurtado Cruz

Ecuaciones Diferenciales 1 3
Unidad 5 Ecuaciones lineales de segundo orden con coecientes variables
5.1 Ecuación de Bessel

Ejemplo A continuación, considere el circuito LRC paralelo que se muestra en la Figura. Sea V la caida
de voltaje a través del capacitor e I la corriente que pasa por la inductancia.

Podemos mostrar que el voltaje y la corriente están descritas por el sistema de ecuaciones.
dI V
=
dt L
dV I V
=− −
dt C RC
en donde L es la inductancia, C la capacitancia y R la resistencia
Una razón por la que los sistemas de ecuaciones de primer orden son particularmente importantes es que
Las ecuaciones de orden superior siempre se pueden transformar en tales sistemas.
Ejemplo El movimiento de cierto sistema resorte-masa se describe por la ecuación diferencial de segundo
orden
u00 + 0,125u0 + u = 0 (8)
Escribir de nuevo esta ecuación como un sistema de ecuaciones de primer orden.
Solución Sean x1 = u y x2 = u0 ; entonces se deduce que x01 = x2 . además, u00 = x02 . Al sustituir u, u0 y
u de la ecuación (8) se obtiene
00

x01 = x2
x02 = −x1 − 0,125x2

ejemplo La ecuación general de movimiento de un sistema resorte-masa,

mu00 + γu0 + k u = 0 (9)
puede transformarse en un sistema de ecuaciones de primer orden de la misma manera. Si se hace
x1 = u y x2 = u0 y se procede como en el ejemplo anterior y se obtiene el sistema

x01 = x2
 
0 k γ F (t)
x2 = x1 − x2 +
m m m

Facultad de Ciencias UNAM Prof. Esteban Rubén Hurtado Cruz

Ecuaciones Diferenciales 1 4