Las matemáticas son el arte de dar el mismo nombre a diferentes cosas.
– Henri Poincare
Introducción
Ya conocemos las propiedades de los sistemas de ecuaciones diferenciales lineales de primer orden, así como las de sus soluciones. Mucho de lo que vimos en las dos entradas anteriores es bastante similar a lo que desarrollamos en las dos primeras entradas de la unidad 2, sin embargo en esta entrada aprenderemos algo completamente nuevo, definiremos una herramienta matemática que nos será de bastante utilidad. ¡Se trata de la exponencial de una matriz!.
En esta entrada definiremos lo que significa , donde es una matriz de con componentes reales constantes.
Así mismo, estudiaremos algunas de sus propiedades y su utilidad en la resolución de sistemas lineales.
Operaciones sobre componentes de matrices
Muchas de las operaciones que se pueden hacer hacía una matriz son aplicables sobre cada una de las componentes que conforman a dicha matriz. Para comprender este hecho es conveniente definir lo que significa la derivada e integral de una matriz, esto nos permitirá ganar intuición.
Consideremos por un momento una matriz de compuesta de funciones.
donde , son funciones sobre algún intervalo común . Comencemos por definir la derivada de una matriz.
Definición: Sea una matriz de con componentes , . Si cada componente es derivable en un intervalo , entonces se define la derivada de la matriz como la matriz con entradas conformadas por las derivadas de las funciones , es decir,
y se denota por
Algunas reglas de derivación se enuncian a continuación.
Teorema: Sean y matrices con componentes dependientes de y sea una matriz constante. Entonces se cumplen las siguientes operaciones.
De manera equivalente se puede definir la integral de una matriz.
Definición: Sea una matriz de con componentes , . Si cada componente es integrable en un intervalo , entonces se define la integral de la matriz como la matriz con entradas conformadas por las integrales de las funciones , es decir,
y se denota por
Ejemplo: Calcular la derivada de la matriz
Solución: Aplicamos la derivada sobre cada componente de la matriz.
De manera similar se puede hacer la integral de una matriz.
Definamos lo que es una serie de matrices. En este caso consideremos matrices constantes.
Definición: Se define una serie de matrices de tamaño con coeficientes constantes como la matriz de series
donde es la entrada de la matriz y se denota por
Si tenemos series de matrices es claro que algunas pueden o no converger. A continuación definimos la convergencia en matrices.
Definición: Si cada serie es convergente, diremos que la serie converge a la matriz conformada por los valores a los cuales convergen las series .
Con esto en mente puede resultar más comprensible la definición de la exponencial de una matriz.
Exponencial de una matriz
Recordemos que la función escalar de la exponencial se define como
Con las definiciones anteriores podemos extender la serie de la exponencial anterior a una serie de matrices.
Definición: Sea una matriz de con coeficientes constantes (). Se define la exponencial de como la serie convergente
Se puede demostrar que la serie () converge, sin embargo se requiere de un poco más de teoría que queda fuera de nuestro interés.
Definición: Si la exponencial depende de , se define la exponencial de la matriz como
Con la matriz identidad, también y recordemos que , etcétera.
Veamos un ejemplo en el que determinemos la exponencial de una matriz.
Ejemplo: Determinar la matriz , en donde
Solución: Para determinar la matriz usemos directamente la definición (). Sabemos que
Ahora bien,
Sustituimos en ().
¡Uff!. En la última igualdad lo que hicimos fue multiplicar por un uno en la forma , esto nos permitió hacer que la potencia de los sea la misma que en el factorial.
Escribamos la última matriz como series infinitas.
Nos gustaría hacer que las series comiencen en . Notemos que, de acuerdo a la forma de la serie, el termino daría como resultado un , considerando esto la expresión anterior la podemos escribir como
Para las componentes de la matriz en las que la serie aún no comienza en sumamos y restamos un , así el puede ser incluido en la serie.
Ahora todas las series comienzan en . Sabemos que la serie converge a un número.
Por definición de convergencia en matrices, se tiene
Por lo tanto, la matriz que buscamos es
Como se puede notar, calcular la exponencial de una matriz usando la definición puede ser una tarea bastante tediosa. Por su puesto existen métodos que nos permiten calcular este tipo de matrices de forma más sencilla, más adelante revisaremos uno de ellos.
Algunas propiedades de la exponencial de una matriz se enuncian a continuación.
Teorema: Sean y matrices constantes de y sean y constantes, entonces
Con el vector cero e la matriz identidad.
Con la matriz identidad.
Siempre que .
No demostraremos este teorema ya que nuestro principal interés está en conocer como estos conceptos y propiedades se pueden aplicar en nuestro estudio sobre sistemas lineales.
A continuación mostraremos un resultado importante e interesante y es que la función () ¡es solución del sistema lineal homogéneo !.
La exponencial de una matriz y los sistemas lineales
Teorema: Sea una matriz de con coeficientes constantes. La función es solución del sistema lineal homogéneo .
Demostración: Consideremos la función . Apliquemos la derivada término a término de la definición ().
Como las matrices y son constantes, entonces se tiene lo siguiente.
Con esto hemos mostrado que
Es decir,
La ecuación () no sólo prueba que es solución del sistema lineal, sino que además muestra cuál es la derivada de la matriz .
Veamos un problema de valores iniciales.
Teorema: Sea una matriz de con coeficientes constantes. La función
es solución del problema de valores iniciales
Donde es un vector de constantes.
Demostración: Consideremos la función
con un vector constante, si la derivamos obtenemos lo siguiente.
En donde se ha hecho uso del resultado (). Esto muestra que la función es solución del sistema .
Si tomamos y considerando que , se tiene
Es decir, se satisface la condición inicial .
Nota: Es claro que la matriz es una matriz de , mientras que el vector constante es un vector de , así que es importante el orden de las matrices, la función del teorema anterior es la cual esta bien definida y el resultado es una matriz de , mientras que la operación no está definida de acuerdo al algoritmo de multiplicación de matrices. Cuidado con ello.
Para concluir esta entrada retomaremos el concepto de matriz fundamental de soluciones visto en la entrada anterior.
Matriz fundamental de soluciones
Retomemos la definición de matriz fundamental de soluciones.
Definición: Sean , soluciones linealmente independientes del sistema lineal homogéneo . Se le denomina matriz fundamental de soluciones del sistema a la matriz conformada por los vectores solución.
En esta sección denotaremos por a una matriz fundamental de soluciones.
El objetivo de esta sección es mostrar que se puede determinar directamente la matriz a partir de cualquier matriz fundamental de soluciones. Antes de llegar a este resultado veamos unos resultados previos.
Teorema: Una matriz es una matriz fundamental de soluciones de si y sólo si satisface al sistema y el Wronskiano es distinto de cero.
Demostración:
) Por demostrar: satisface al sistema y se cumple que .
Supongamos que es una matriz fundamental de soluciones del sistema conformada por los vectores solución , es decir,
Si es la matriz de coeficientes (), entonces
La -ésima columna del resultado de multiplicar estas matrices es
Identificamos que esta matriz corresponde a la derivada ya que cada es solución del sistema, es decir,
Entonces tiene por columnas a los vectores . Por lo tanto
es decir, la matriz fundamental satisface al sistema.
Como es una matriz fundamental de soluciones, entonces los vectores que la componen son linealmente independientes, de manera que el Wronskiano es distinto de cero.
en particular se cumple para , es decir,
Demostremos el regreso.
) Por demostrar: es una matriz fundamental de soluciones.
Sea una matriz compuesta por los vectores y cuya derivada es
entonces las columnas satisfacen
así es solución del sistema
para .
Por otro lado, por hipótesis
Como es distinto de cero en un punto, entonces lo debe ser en todo el dominio, es decir, se cumple que
lo que significa que son linealmente independientes.
De ambos resultados concluimos que la matriz es una matriz fundamental de soluciones.
Veamos un resultado interesante. Ya vimos que la matriz es solución del sistema , pero no solo resulta ser solución, sino que además ¡es una matriz fundamental de soluciones!.
Teorema: La función es una matriz fundamental de soluciones del sistema .
Demostración: Anteriormente mostramos que
lo que prueba que es solución del sistema .
Supongamos que está compuesta por la matriz de vectores . Si , se tiene que
y además el determinante es distinto de cero.
o bien,
Por el teorema anterior concluimos que es una matriz fundamental de soluciones del sistema lineal.
Veamos un resultado más antes de llegar a nuestro objetivo.
Teorema: Si y son matrices fundamentales de soluciones del sistema , entonces existe una matriz constante , tal que
Demostración: Sean y matrices fundamentales del sistema . Supongamos que
y
Donde la -ésima columna de las matrices anteriores son
Como ambas matrices son matrices fundamentales de soluciones, entonces cada y cada son linealmente independientes. Esto nos permite escribir cualquier columna de como combinación lineal de las columnas de , es decir, existen constantes , tales que
Donde el índice de las constantes indica el número de columna de la matriz , es decir, si definimos el vector
podemos escribir la combinación lineal () como
Definamos la matriz constante como
En forma matricial la operación corresponde al producto de las siguientes matrices de .
Observemos con cuidado que el resultado () corresponde a la -ésima columna de hacer el producto . Por lo tanto, las ecuaciones () son equivalentes al sistema
Hemos llegado al resultado final. Dicho resultado involucra el concepto de matriz inversa, recordemos este concepto de álgebra lineal.
Definición: Sean y dos matrices de con componentes constantes. Suponiendo que
Con la matriz identidad. Entonces se llama la matriz inversa de y se denota por .
La ecuación anterior queda como
Demostremos el resultado que nos muestra cómo obtener la matriz a partir de una matriz fundamental de soluciones.
Teorema: Si es una matriz fundamental de soluciones del sistema , entonces
donde es la matriz inversa de evaluada en .
Demostración: Sabemos que y son matrices fundamentales de soluciones de , de acuerdo al teorema anterior ambas funciones se relacionan de la siguiente forma.
para alguna matriz constante .
Tomemos , por un lado
Por otro lado, de ()
De ambas ecuaciones tenemos que
Esta ecuación obliga que
Sustituyendo en () concluimos que
Finalicemos con un ejemplo.
Ejemplo: Determinar la matriz , donde
Solución: Imagina lo complicado que sería este problema si lo intentáramos resolver usando la definición ().
En la entrada anterior vimos que una matriz fundamental de soluciones del sistema lineal
es
Determinemos la matriz usando la expresión ().
Calcular la matriz inversa puede ser una tarea larga y debido a que no es el objetivo llevar a cabo estas operaciones se puede comprobar que la matriz inversa de es
Basta probar que . Para calcular la inversa se puede hacer uso de algún programa computacional. Si en la matriz inversa evaluamos en , se tiene
Por lo tanto,
Haciendo la multiplicación de matrices correspondiente obtenemos finalmente que
Existen otras formas de calcular la exponencial de una matriz, una de ellas es usando la transformada de Laplace y otra puede ser diagonalizando matrices. Si lo deseas puedes investigar sobre estos métodos en la literatura, por nuestra parte serán temas que no revisaremos ya que, más que obtener exponenciales de matrices, estamos interesados en obtener soluciones de sistemas de ecuaciones diferenciales, aunque como vimos, están estrechamente relacionados.
Tarea moral
Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.
- Sea
- Demostrar que la derivada de la matriz es
- Demostrar que la integral de a de la matriz es
- Demostrar que
Donde,
con una constante.
- Obtener la matriz para los siguientes casos:
- Sea una matriz fundamental de soluciones del sistema . Demostrar que
- Una matriz fundamental del sistema
es
- Demostrar que la matriz anterior en efecto es una matriz fundamental de soluciones del sistema.
- Demostrar que la matriz está dada por
Más adelante…
En estas tres primeras entradas de la unidad 3 establecimos la teoría básica que debemos conocer sobre los sistemas lineales de primer orden compuestos por ecuaciones diferenciales lineales de primer orden. En particular, esta entrada es de interés, pues más adelante la exponencial de una matriz volverá a aparecer cuando estudiemos métodos de resolución y cuando justifiquemos los teoremas de existencia y unicidad.
En las siguientes entradas comenzaremos a desarrollar los distintos métodos de resolución de estos sistemas lineales. En particular, en la siguiente entrada desarrollaremos el método de eliminación de variables, éste método en realidad es muy sencillo, útil y práctico en muchas ocasiones, aunque también es un método muy limitado.
Entradas relacionadas
Agradecimientos
Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»