Introducción

Vamos a concluir la tercera unidad del curso revisando el teorema de existencia y unicidad para sistemas de ecuaciones de primer orden, en su forma general, es decir, para sistemas lineales y no lineales que satisfagan las hipótesis del teorema. Hasta el momento únicamente demostramos el teorema de existencia y unicidad para sistemas lineales con coeficientes constantes, pero es importante demostrar la versión general al igual que hicimos para las ecuaciones de primer orden.

Lo primero que veremos es que un sistema de ecuaciones de la forma $$\begin{alignedat}{4} \dot{x}_{1} &= F_{1}(t,x_{1},x_{2},…,x_{n}) \\ \dot{x}_{2} &= F_{2}(t,x_{1},x_{2},…,x_{n}) \\ & \; \; \vdots \notag \\ \dot{x}_{n} &= F_{n}(t,x_{1},x_{2},…,x_{n}) \end{alignedat}$$ se puede escribir en forma abreviada como sigue: $$\dot{\textbf{X}}(t)=\textbf{F}(t,\textbf{X}(t))$$ donde $\textbf{F}$ es el vector conformado por las funciones $F_{i}$ del sistema, con $i \in \{1,…,n\}$. Si además agregamos la condición inicial $\textbf{X}(t_{0})=\textbf{Y}$, entonces podemos ver que el sistema se reduce a una expresión muy similar al problema de condición inicial $$\frac{dy}{dt}=f(t,y(t)) \,\,\,\,\, ; \,\,\,\,\, y(t_{0})=y_{0}$$ salvo que ahora $\textbf{X}$ es una función que toma valores en $\mathbb{R}^{n}$, y $\textbf{F}$ es una función de $\mathbb{R}^{n+1}$ a $\mathbb{R}^{n}$.

Afortunadamente la mayoría de los lemas y teoremas que usamos para demostrar el teorema de existencia y unicidad para ecuaciones de primer orden se pueden extender a funciones de varias variables, por lo que la demostración será muy similar a la demostración de este último teorema.

Antes de iniciar te dejo la entrada correspondiente al teorema de existencia y unicidad de Picard, para que te familiarices con él y te sea más fácil ver los videos de esta entrada.

El teorema de existencia y unicidad para sistemas de ecuaciones de primer orden. Ecuación integral asociada

Enunciamos el teorema de existencia y unicidad para sistemas de ecuaciones de primer orden, analizamos las similitudes que existen con el teorema de existencia y unicidad de Picard, y vemos que resolver el problema de condición inicial es equivalente a resolver la ecuación integral $$\textbf{X}(t)=\textbf{Y}+\int_{t_{0}}^{t} \textbf{F}(s, \textbf{X}(s)) \, ds.$$

Demostración de la existencia de la solución al problema de condición inicial

Demostramos la existencia de una solución al problema de condición inicial estudiando bajo qué circunstancias converge uniformemente la sucesión de iteraciones de Picard del problema. En dado caso que esto último suceda, la función a la cual convergen las iteraciones será solución a la ecuación integral del video anterior.

Demostración de la unicidad de la solución al problema de condición inicial

Concluimos la demostración del teorema probando la unicidad de la solución al problema de condición inicial.

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.

• Sea $\textbf{F}(t,\textbf{X}(t))$ continua en un dominio $E \subseteq \mathbb{R}^{n+1}$ que contenga a $(t_{0},\textbf{Y})$. Demuestra que $\textbf{X}(t)$ es solución al problema de condición inicial $$\dot{\textbf{X}}(t)=\textbf{F}(t,\textbf{X}(t)) \,\,\,\,\, ; \,\,\,\,\, \textbf{X}(t_{0})=\textbf{Y}$$ si y sólo si es solución a la ecuación integral $$\textbf{X}(t)=\textbf{Y}+\int_{t_{0}}^{t} \textbf{F}(s,\textbf{X}(s)) \, ds.$$

• Considera el problema de condición inicial $$\dot{\textbf{X}}=\begin{pmatrix} 0 & 1\\ 1 & 0 \end{pmatrix} \textbf{X} + \begin{pmatrix} t \\ t \end{pmatrix} \,\,\,\,\, ; \,\,\,\,\, \textbf{X}(0)=\begin{pmatrix} 1 \\ -1 \end{pmatrix}.$$ Calcula las iteraciones de Picard correspondientes al problema. ¿Convergen a alguna función? En caso afirmativo, muestra que dicha función es solución al problema de condición inicial.

• Supongamos que $\textbf{F}(t,\textbf{X}(t))$ es continua en $$R:=\{(t,x_{1},…,x_{n}) \in \mathbb{R}^{n+1} : |t-t_{0}| \leq a, \lVert \textbf{X}(t) – \textbf{Y} \rVert \leq b, \, \, a, b \in \mathbb{R}\}.$$ Demuestra que existe $M > 0$ y $h \in \mathbb{R}$ tal que $$\lVert \textbf{X}^{n}(t)-Y \rVert \leq M |t-t_{0}|, \forall n \in\mathbb{N}, \forall t \in I_{h} \subseteq \mathbb{R}.$$ Recuerda que $\textbf{X}^{n}(t)$ es la $n$-ésima iteración de Picard correspondientes al problema de condición inicial que estudiamos a lo largo de la entrada. (Hint: La prueba es similar al lema análogo que probamos en este video para el teorema de existencia y unicidad de Picard).

• Consideremos el problema de condición inicial $$a\frac{d^{2}y}{dt^{2}}+b\frac{dy}{dt}+cy=0 \,\,\,\,\, ; \,\,\,\,\, y(t_{0})=y_{0} \,\,\,\,\, ; \,\,\,\,\, \frac{dy}{dt}(t_{0})=y_{1}$$ con $a,b,c$ constantes. ¿Si el sistema de ecuaciones asociado satisface el teorema de existencia y unicidad, entonces el problema de condición inicial original tiene una única solución?

Más adelante

Con este teorema finalizamos la tercera unidad del curso. En la cuarta unidad comenzaremos con la teoría cualitativa de los sistemas de ecuaciones de primer orden.

Veremos que los sistemas tienen puntos de equilibrio, los clasificaremos según su estabilidad. En virtud de esto vamos a analizar el comportamiento de las soluciones cerca de puntos de equilibrio y dibujaremos el plano fase de un sistema.

Abordaremos sistemas no lineales, y aunque no los resolveremos explícitamente, veremos el comportamiento de sus soluciones cerca de sus puntos de equilibrio.

Finalmente, veremos algunos sistemas que satisfacen propiedades interesantes, como los sistemas Hamiltonianos, los disipativos, entre otros.

Entradas relacionadas

• Ir a Ecuaciones Diferenciales I
• Entrada anterior del curso: Sistemas de ecuaciones lineales no homogéneas. Solución por variación de parámetros
• Siguiente entrada del curso: Geometría de soluciones a sistemas de dos ecuaciones de primer orden. Plano fase y campo vectorial asociado

• Notas escritas relacionadas con el tema: Teorema de existencia y unicidad para sistemas de ecuaciones diferenciales de primer orden

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»