Introducción
En la entrada anterior definimos la exponencial de una matriz $\textbf{A}$ de coeficientes constantes, denotada por $\textbf{e}^{\textbf{A}}$, demostramos sus principales propiedades, y estudiamos la relación que guarda con el sistema lineal de ecuaciones $\dot{\textbf{X}}=\textbf{A}\textbf{X}$ y su matriz fundamental de soluciones. Con esta herramienta a nuestra disposición, podremos enunciar y demostrar el teorema de existencia y unicidad para sistemas lineales de primer orden con coeficientes constantes.
Como mencionamos en la entrada anterior, nuestra meta es tratar de generalizar la fórmula para soluciones a ecuaciones lineales de primer orden con condición inicial, la cual es de la forma $$y(t)=e^{-\int p(t) dt} \left[\int e^{\int p(t) dt}q(t)+k_{0}\right]$$ para cierta constante $k_{0}$, y encontrar una solución al problema de condición inicial $$\dot{\textbf{X}}=\textbf{A}\textbf{X}+\textbf{Q} \, \, \, \, \, \, \, \, \, \, ; \, \, \, \, \, \, \, \, \, \, \textbf{X}(0)=\textbf{C}$$ que se vea de la forma $$\textbf{X}(t)=\textbf{e}^{-\int \textbf{A}(t) dt} \left[\int \textbf{e}^{\int \textbf{A}(t) dt}\textbf{Q}(t)+\textbf{B}\right].$$
El teorema de existencia y unicidad para sistemas lineales de primer orden nos garantiza la existencia de tal solución. Además, una vez que definimos la exponencial de una matriz, ya no nos sorprenderá la notación de la fórmula anterior. Dividiremos el teorema y su demostración en dos casos: para sistemas homogéneos y para sistemas no homogéneos.
Teorema de existencia y unicidad para sistemas lineales homogéneos de primer orden con coeficientes constantes
En el primer video demostramos el teorema de existencia y unicidad para sistemas lineales homogéneos de primer orden con coeficientes constantes.
Teorema de existencia y unicidad para sistemas lineales no homogéneos de primer orden con coeficientes constantes
En el segundo video demostramos el mismo teorema pero ahora para sistemas lineales no homogéneos de primer orden con coeficientes constantes.
Tarea moral
Los siguientes ejercicios no forman parte de la evaluación del curso, pero te servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.
- Encuentra la solución al problema de condición inicial: $$\dot{\textbf{X}}=\begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix}\textbf{X} \, \, \, \, \, \, \, \, \, \, ; \, \, \, \, \, \, \, \, \, \, \textbf{X}(0)=\begin{pmatrix} 1 \\ 0 \end{pmatrix}.$$
- Encuentra la solución al problema de condición inicial: $$\dot{\textbf{X}}=\begin{pmatrix} 3 & 0 \\ 0 & 3 \end{pmatrix}\textbf{X} \, \, \, \, \, \, \, \, \, \, ; \, \, \, \, \, \, \, \, \, \, \textbf{X}(1)=\begin{pmatrix} 2 \\ 2 \end{pmatrix}.$$
- Encuentra la solución al problema de condición inicial: $$\dot{\textbf{X}}=\begin{pmatrix} 1 & 0 \\ 0 & 5 \end{pmatrix}\textbf{X}+ \begin{pmatrix} t \\ 1 \end{pmatrix} \, \, \, \, \, \, \, \, \, \, ; \, \, \, \, \, \, \, \, \, \, \textbf{X}(0)=\begin{pmatrix} 1 \\ 0 \end{pmatrix}.$$
- Encuentra la solución al problema de condición inicial: $$\dot{\textbf{X}}=\begin{pmatrix} 7 & 0 \\ 0 & 3 \end{pmatrix}\textbf{X}+ \begin{pmatrix} t \\ t \end{pmatrix} \, \, \, \, \, \, \, \, \, \, ; \, \, \, \, \, \, \, \, \, \, \textbf{X}(2)=\begin{pmatrix} 0 \\ 2 \end{pmatrix}.$$
Más adelante
Una vez que hemos encontrado formas explícitas para las soluciones a sistemas lineales con coeficientes constantes $\dot{\textbf{X}}=\textbf{A}\textbf{X}+\textbf{Q}$, debemos encontrar algún método para calcular eficientemente $\textbf{e}^{t\textbf{A}}$, sin pasar por el complicado camino de calcular cada serie que conforma a la exponencial de $t\textbf{A}$. El método que desarrollaremos es una aplicación de los eigenvalores y eigenvectores (o valores y vectores propios) que quizá hayas visto en cursos de álgebra lineal.
Es por eso que, aunque no estamos en un curso de álgebra lineal, haremos un alto en el camino y revisaremos de manera muy breve estos conceptos y demás herramientas que utilizaremos muy pronto. Iremos relacionando los conceptos con los temas que nos interesan, que son los de hallar una matriz fundamental de soluciones, la exponencial de una matriz, y por supuesto resolver sistemas lineales de primer orden.
Entradas relacionadas
- Ir a Ecuaciones Diferenciales I
- Entrada anterior del curso: La exponencial de una matriz
- Siguiente entrada del curso: Valores y vectores propios, polinomio característico y la exponencial de una matriz
- Notas escritas relacionadas con el tema: Teorema de existencia y unicidad para sistemas lineales
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»