Sistemas de Ecuaciones Diferenciales

Homogéneos
E. Castilla y P.J. Chocano (Departamento de Matemática Aplicada)

La idea de esta entrega es aprender cómo nociones básicas de álgebra lineal tienen aplica-
ciones directas en otras áreas, en este caso, al estudio de ecuaciones diferenciales.

Introducción
Las aplicaciones de las ecuaciones diferenciales a la hora de modelizar ciertos fenómenos son
amplias: fı́sica, biologı́a, quı́mica, ingenierı́a, etc. Para el lector interesado en alguna de ellas
se recomiendan las referencias [2,3].
En esta entrega asumiremos conocimientos básicos de ecuaciones diferenciales ordinarias.
Nuestro objeto de estudio para esta entrega son los sistemas de ecuaciones diferenciales, es
decir, sistemas de la siguiente forma:


dy1

 dx
= a11 (x)y1 + a12 (x)y2 + ... + a1n (x)yn + b1 (x)

 dy2

= a21 (x)y1 + a22 (x)y2 + ... + a2n (x)yn + b2 (x)
dx
..


 .
 dyn

= an1 (x)y1 + an2 (x)y2 + ... + ann (x)yn + bn (x)
dx

donde las funciones aij y bi con i, j = 1, ..., n son funciones continuas en un cierto intervalo I.
A dicho sistema le llamaremos sistema lineal de orden n. Si b1 = b2 = ...bn = 0, diremos que
dicho sistema es homogéneo, en caso contrario, no homogéneo. Podemos reescribir dicho
sistema de manera matricial:
      
y10 a11 (x) a12 (x) . . . a1n (x) y1 b1 (x)
 y 0   a21 (x) a22 (x) . . . a2n (x)   y2   b2 (x) 
 2 
 ..  =  .. ..   ..  +  .. 
   
.. . .
.  . . . .  .   . 
0
yn an1 (x) an2 (x) . . . ann (x) yn bn (x)

Identificando
     
y1 a11 (x) a12 (x) . . . a1n (x) b1 (x)
 y2   a21 (x) a22 (x) . . . a2n (x)   b2 (x) 
y(x) =  ..  , A(x) =  .. ..  y b(x) =  .. 
     
.. . .
.  . . . .   . 
yn an1 (x) an2 (x) . . . ann (x) bn (x)

tenemos que el sistema de ecuaciones diferenciales queda:

y0 = A(x)y + b(x). (1)

Veamos un ejemplo concreto y sencillo.

1
Ejemplo 1
(
y10 = 4y1 − 3y2
y20 = 6y1 − 7y2

En este caso concreto, de manera matricial tenemos:

 0   
y1 4 −3 y1
= .
y20 6 −7 y2

Teorema 1 (Teorema de existencia y unicidad). Sean aij y bi funciones continuas en un

intervalo I donde i, j = 1, ..., n y sea x0 ∈ I e y10 , ..., yn0 ∈ R. Entonces existe una única
solución (y1 , ..., yn ) a el sistema (1) en el intervalo I que satisface las condiciones iniciales
y1 (x0 ) = y10 , ..., yn (x0 ) = yn0 .

La demostración de este resultado al igual que futuros resultados queda lejos del alcance para
los objetivos de esta entrega. El lector interesado puede consultar dicha demostración (y pos-
teriores) en la siguiente referencia [1].

Observación 1 Nótese que si el vector de coeficientes independientes b es nulo en el sistema

(1), entonces siempre tenemos que y = (0, ..., 0) es una solución a nuestro sistema.

Método de los autovalores para sistemas lineales homogéneos con

coeficientes constantes
El objetivo de esta sección es obtener soluciones a problemas de la siguiente forma:

y0 = Ay. (2)

La idea es extender el método que vimos para la resolución de ecuaciones diferenciales

lineales homogéneas, es decir, vamos a buscar soluciones de la forma:

y(x) = λeλx v

donde λ ∈ R y v es un vector de Rn . Si sustituimos nuestro candidato a solución en la ecuación

(2), obtenemos
λeλx v = A(eλx v)
dividiendo ambos miembros por la exponencial, llegamos a

Av = λv.

Por tanto, para que y(x) = λeλx v sea solución ha de verificarse la condición previa. Además,
también tenemos que dicha condición es suficiente. Ası́ pues, hemos convertido un problema
de hallar soluciones de un sistema de ecuaciones diferenciales en un problema algebraico de
encontrar autovalores y autovectores para una matriz A.
Ahora pasaremos a ver las posibles opciones que aparecen dependiendo de cómo sean los
autovalores.

Teorema 2 Si la matriz A tiene n autovalores, λ1 , ..., λn , reales distintos dos a dos, y v1 , ..., vn
son autovectores correspondientes a los autovalores, entonces la solución general del sistema (2)
viene dada por
y = c1 eλ1 x v1 + c2 eλ2 x v2 + ... + cn eλn x vn
donde c1 , ..., cn ∈ R son constantes cualesquiera.

2
Ejemplo 2 Recuperemos nuestro anterior ejemplo:
 0   
y1 4 −3 y1
0 = .
y2 6 −7 y2

Estudiamos los autovalores y autovectores de

 
4 −3
A= .
6 −7

Es sencillo comprobar que su polinomio caracterı́stico es λ2 − 3λ − 10. Por tanto, tenemos

dos autovalores: λ1 = 5 y λ2 = −2. Asociado a cada autovalor podemos obtener los siguientes
autovectores: vt1 = (2, 1) y vt2 = (1, −3). Por tanto, la solución general a dicho sistema viene
dada por:      
y1 (x) 2 5x 1
= c1 e + c2 e−2x .
y2 (x) 1 −3
Escrito de manera escalar:
(
y1 (x) = 2c1 e5x + c2 e−2x
y2 (x) = c1 e5x + −3c2 e−2x

Teorema 3 Si la matriz A tiene n autovalores distintos, λ1 , ..., λn−2m reales y α1 ±β1 i, ..., αm ±
βm i complejos, y v1 , ..., vn−2m son autovectores (reales) correspondientes a los autovalores reales
y a1 ± b1 i, ..., am ± bm i son autovectores correspondientes a los autovalores complejos, entonces
la solución general al sistema (2) viene dada por:
n−2m
X n
X n
X
λi x αj x
y(x) = ci e vi + e (aj cos(βk x)−bj sin(βj x))+ eαj x (aj cos(βk x)+bj sin(βj x))
i=1 j=n−2m j=n−2m

Ejemplo 3 Consideremos el sistema (2) donde

 
1 0 0
A = 0 2 −3
1 3 2

Los autovalores de A son 1, 2 ± 3i y sus autovectores asociados

   
10 0
−3 y  i 
−1 1

Por tanto, la solución general viene dada escrita de manera escalara por

y1 (x)
 = 10c1 ex
y2 (x) = −3c1 et + e2x (−c2 sin(3x) + c3 cos(3x))

y3 (x) = −c1 ex + e2x (c2 cos(3x) + c3 sin(3x))


El mismo resultado tenemos también en caso de que tengamos autovalores repetidos si

coindicen las multiplicidades geométricas y algebraicas.

3
Teorema 4 Si la matriz A tiene n autovalores λ1 , ..., λn−2m reales y α1 ±β1 i, ..., αm ±βm i com-
plejos de manera que las multiplicidades algebraicas y geométricas coinciden, y v1 , ..., vn−2m son
autovectores (reales) correspndientes a los autovalores reales y a1 ± b1 i, ..., am ± bm i son au-
tovectores correspondientes a los autovalores complejos, entonces la solución general al sistema
(2) viene dada por:
n−2m
X n
X n
X
y(x) = ci eλi x vi + eαj x (aj cos(βk x)−bj sin(βj x))+ eαj x (aj cos(βk x)+bj sin(βj x))
i=1 j=n−2m j=n−2m

Ejemplo 4 Consideremos el sistema (1) donde

 
9 4 0
A = −6 −1 0 .
6 4 3

Su polinomio caracterı́stico es (λ − 3)2 (λ − 5). Por tanto, tenemos dos autovalores teniendo
uno de ellos multiplicidad algebraica 2. La multiplicidad geométrica de dicho autovalor también
es 2 y tenemos que dos autovectores asociados son (2, −3, 0) y (0, 0, 1). Además, un autovector
asociado a λ = 5 es (1, −1, 1). Por tanto, la solución general a el sistema viene dada por

3x 5x
y1 (x) = 2c1 e + c3 e

y2 (x) = −3c1 e3x − c3 e5x

y3 (x) = c2 e3x + c3 e5x


Observación 2 Nótese que hasta ahora solamente hemos analizado los casos en los que la
matriz A es diagonalizable.

Para los casos en los que la multiplicidad algebraica no coincida con la multiplicidad
geométrica, basta con aplicar el proceso de obtención de forma canónica de Jordan. La base
que se emplea para llegar a la forma canónica de Jordan nos da la clave para obtener la solución
general del sistema. Para el lector interesado en el proceso completo se recomienda [1,3].

Caso General: idea

Nos falta analizar el caso más general dentro de los sistemas de ecuaciones con coeficientes
constantes, que es cuando el vector b es no nulo. Para dicho caso, tenemos “adaptaciones” de
lo que sucede para ecuaciones diferenciales lineales: método de variación de las constantes y
método de los coeficientes indeterminados. Ambos siguen la misma filosofı́a que los vistos en
clase. Aunque a nivel técnico la dificultad aumenta. No obstante la idea es la misma, resolver
el sistema homogéneo asociado y a continuación encontrar una solución particular.

Entrega
Supongamos que tenemos tres tanques de salmuera A, B y C cuyos volúmenes son 20, 40 y
60, respectivamente. Supongamos que se encuentran conectados en cascada mediante un tubo
siendo A el que está más alto, B el que se encuentra en medio y C el que está e una posición más
baja. Supongamos que un lı́quido entra en el tanque A, se desplaza al tanque B, llega al tanque
C y finalmente sale. Para simplificar el problema vamos a asumir que la proporción de lı́quido

4
en los cambios que se producen es constante. Además, asumimos que hay una concentración
de sal constante en cada tanque. Consideremos las siguiente variables:

y1 (x) = Cantidad de sal en tanque A a tiempo x

y2 (x) = Cantidad de sal en tanque B a tiempo x
y3 (x) = Cantidad de sal en tanque C a tiempo x
x= Tiempo en minutos

Supongamos que introducimos agua sin sal en el primer tanque (A). Por tanto, la sal en el
resto de tanques eventualmente se perderá. La cascada de tanques viene modelizada por la ley
quı́mica del balance. El problema queda modelizado por las siguientes ecuaciones diferenciales:

0 1
y1 = − 2 y1

y20 = 21 y1 − 14 y2

 0
y3 = 41 y2 − 16 y3 .

Ejemplo 5 Realiza las siguientes tareas:

1. Halla la solución general.

2. Suponiendo que y(x) = (y1 (x), y2 (x), y3 (x)) es la solución obtenida en el anterior apartado.
Estudia limx→∞ y(x). ¿Qué podemos decir de la concentración de sal en los tanques según
pasa el tiempo?

Referencias
1. Arnol’d, V. Ordinary Differential Equations Springer-Verlag, 1984.

2. Gustafson, G. B. Differential Equations and Linear Algebra. Undergraduate Mathematics

Science and Engineering, 2022.

3. Zill, D. y Cullen, M.R. Ecuaciones diferenciales con problemas con valores en la frontera.
Thomsom.