Un sistema de ecuaciones diferenciales es un conjunto de varias ecuaciones diferenciales con varias funciones incógnitas y un conjunto de condiciones de contorno. Una solución del mismo es un conjunto de funciones diferenciables que satisfacen todas y cada una de las ecuaciones del sistema. Según el tipo de ecuaciones diferenciales puede tenerse un sistema de ecuaciones diferenciales ordinarias o un sistema de ecuaciones en derivadas parciales.
Sistemas de ecuaciones diferenciales ordinarias
En un sistema de ecuaciones diferenciales ordinarias de cualquier orden, puede ser reducido a un sistema equivalente de primer orden, si se introducen nuevas variables y ecuaciones. Por esa razón en este artículo sólo se consideran sistemas de ecuaciones de primer orden. Un sistema de ecuaciones diferenciales ordinarias de primer orden escrito en forma explícita es un sistema de ecuaciones de la forma:
Reducción a un sistema de primer orden
Dado un sistema de ecuaciones diferenciales de orden n con m ecuaciones:
Existe un sistema equivalente de primer orden con a lo sumo (n+1) x m ecuaciones. Para ver esto consideremos un sistema en el que intervienen m funciones incógnitas xi y sus n derivadas, e introduzcamos un nuevo conjunto de variables yi,k definidos de la siguiente manera:
El sistema de primer orden equivalente en las variables yi,k resulta ser:
Como ejemplo de reducción de un sistema de ecuaciones diferenciales podemos considerar las ecuaciones de movimiento de la mecánica newtoniana de una partícula que es un sistema de segundo orden con tres ecuaciones:
Si se introducen tres funciones incógnita nuevas que representan la velocidad, el sistema anterior se puede reducir a un sistema de primer orden y seis ecuaciones:
Existencia y unicidad de la solución
Un sistema de ecuaciones diferenciales general tiene la forma
(*)
,
donde:
es una función vectorial.
es una función matricial.
Aplicando el Teorema de Picard-Lindelöf a los sistemas de ecuaciones diferenciales ordinarias de la forma (*), se establece la existencia y unicidad de la solución de dicho sistema en las que tanto la matriz
como la función
sean continuas en un intervalo compacto
.
La demostración se basa en el hecho de que la función
es Lipschitz respecto a
, para todo
si
. A partir del resultado del teorema se puede concluir que la sucesión de Picard
converge a la solución.
Sistemas lineales de ecuaciones diferenciales ordinarias
Un sistema lineal de ecuaciones diferenciales con coeficientes constantes es un sistema de la forma:
Que de forma abreviada denotaremos como:
Cuando el sistema anterior verifica que
decimos que el sistema anterior es homogéneo, en caso contrario diremos que es no homogéneo.
Llamamos sistema fundamental de ecuaciones al conjunto de
soluciones independientes del sistema homogéneo. (Obsérvese que dichas soluciones son vectores).
La matriz que cuyas columnas representan cada una de las soluciones anteriores recibe el nombre de matriz fundamental
y verifica que
.
Si dada una matriz fundamental existe un
tal que la matriz resultante es la identidad diremos que dicha matriz es la matriz principal en
.
La resolución de estos sistemas es muy diversa. Sin embargo, la resolución de sistemas de ecuaciones diferenciales se simplifica mucho si las ecuaciones son de coeficientes constantes, o lo que es lo mismo, un sistema autónomo.
Resolución de sistemas de n ecuaciones lineales de primer orden con coeficientes constantes
El problema homogéneo
En el caso de los sistemas lineales, si la matriz del sistema es de coeficientes constantes podemos resolver el sistema usando el método de los autovalores.
Sea el sistema lineal homogéneo:
, con
Comenzamos encontrando los puntos críticos del sistema, es decir, aquellos valores
donde:
Si
entonces solo hay un punto crítico. Pero si
, entonces existe
tal que
,
es una recta por lo que tenemos rectas de puntos críticos. Estas son las soluciones fijas, no cambian con el paso del tiempo.
Por otra parte, si
es el autovalor de
correspondiente al vector
tenemos
Cada uno de los
autovalores de la matriz nos dará información para construir una solución
para el sistema.
La solución anterior depende de la naturaleza del autovalor:
Raíz real única: En este caso la solución para un valor
concreto viene dada directamente por
.
Raíz real de multiplicidad
: En este caso no podemos expresar la solución como antes, ya que si lo hacemos tendríamos información redundante.
En este caso las soluciones para un valor de
se construyen multiplicando la función exponencial por
como sigue:
La solución es la combinación lineal de las anteriores:
Raíces complejas: Hay que tener en cuenta que siempre vienen en parejas de números conjugados.
En este caso usamos distinguimos la parte real y la imaginaria de la raíz y posteriormente usamos la fórmula de Euler y añadiendo constantes.
La solución general al sistema homogéneo será, finalmente una combinación lineal de todas ellas multiplicadas cada una de ellas por su autovector correspondiente:
, donde
son constantes.
La solución anterior puede reescribirse de forma matricial como:
Si se trata de un problema de valor inicial, los valores de
se determinan a partir del dato inicial.
Ejemplo del problema homogéneo
Consideremos el siguiente problema de valor inicial:
Los valores propios de la matriz son
y
con los autovectores asociados
;
La matriz fundamental asociada a dicho sistema homogéneo sería:
Y la matriz fundamental principal quedaría:
Otro método para resolver este tipo de sistemas es utilizar la exponencial matricial:
Sea
:
Entonces la solución viene dada por la exponencial de la matriz A; para calcularla, utilizando las iteraciones de Picard, distinguiremos los siguientes casos:
La solución general será
.
- A no es diagonal, pero diagonaliza:
Si la matriz A diagonaliza puede expresarse de la forma:
donde D es una matriz diagonal y así tendríamos
.
En este caso la solución general a nuestro problema será
Si la matriz A no diagonaliza en general se puede escribir de la forma
donde J es la forma de Jordan de la matriz A (Forma canónica de Jordan) y ahora sólo se necesita calcular la exponencial de un bloque de Jordan para terminar calculando la exponencial de la matriz J como la matriz que tiene por bloques la exponencial de cada uno de los bloques de Jordan.
Primero descompondremos cada bloque de Jordan en dos matrices más simples:
donde
es una matriz diagonal
que contiene a
en sus elementos diagonales y
una matriz nilpotente
de la forma:
Utilizando lo anterior, podemos definir:
Por tanto,
Quedando así:
. Por tanto, la solución general al problema será:
El problema no homogéneo
Sea el problema no homogéneo:
, con
Entonces la solución general es de la forma
, donde
es la solución general del problema homogéneo y
es una solución particular del no homogéneo.
Una vez obtenida una solución general para el problema homogéneo
, (recordemos que:
, siendo
un vector de constantes), para calcular
utilizamos el método de variación de las constantes, donde buscamos convertir los coeficientes
del vector
en funciones.
Queremos soluciones de la forma
, por lo que imponemos que sea solución de nuestro problema.
Como
es solución fundamental,
.
Si ahora calculamos dicha derivada por la regla de la cadena y utilizando la propiedad anterior de
, tenemos que:
.
Luego volviendo a nuestra ecuación y comparando, nos queda que
. Despejamos:
. Integramos:
Por lo que hemos obtenido que la solución particular del problema no homogéneo es
.
Finalmente:
(Recordemos que los valores de
se determinan a partir del dato inicial del PVI).
Ejemplo del problema no homogéneo
Por último, estudiaremos un problema no homogéneo con autovalores imaginarios para ilustrar otro ejemplo.
Primero hallamos la solución del caso homogéneo:
Los valores propios de la matriz son
. La matriz fundamental principal asociada a este problema es:
Ahora resolvemos la parte no homogénea:
Usando las fórmulas vistas, nuestro vector
resulta:
Por lo visto anteriormente nuestra solución general del problema no homogéneo tendrá la siguiente forma sustituyendo los datos antes calculados:
Referencias
Bibliografía
- Zill, Dennis G. (2006). Ecuaciones diferenciales con aplicaciones (2ª edición). Grupo Editorial Iberoamérica.
- Aranda Iriarte, José Ignacio (2008). Apuntes de ecuaciones diferenciales I. Universidad Complutense de Madrid. Archivado desde el original el 15 de febrero de 2017. Consultado el 2 de mayo de 2023.
- Birkhoff, Garrett; Rota, Gian-Carlo (1978). Ordinary Differential Equations. John Wiley and Sons, Inc. ISBN 0-471-07411-X.
- Gershenfeld, Neil (1999). The Nature of Mathematical Modeling. Cambridge University Press. ISBN 978-0-521-57095-4.
- Robinson, James C. (2004). An Introduction to Ordinary Differential Equations. Cambridge University Press. ISBN 0-521-82650-0.
- Simmons, G.F. (2016). Differential equations with applications and historical notes. CRC Press.
- Teschl, G. (2012). Ordinary equations and dynamical systems. American Mathematical Society.
- Hubbard, J.H.; West, B.H. (2013). Differential equations: A dynamical systems approach. Springer.
- Blanchard, P.; Devaney, R.L. (2012). Differential equations. Cengage Learning.
- Gil, O. (1999). Curso introductorio a las ecuaciones diferenciales. Facultad de ingeniería de la Universidad de la República de Uruguay.
- Gutermuth, Denise. Picard’s existence and uniqueness theorem.
Véase también
Datos: Q6787870