Majorización de ODEs: comparación de soluciones

Mayorización de Sistemas de EDO: Comparación de Soluciones mediante los Lados Derechos

En el análisis de sistemas dinámicos, las soluciones exactas de las EDO a menudo son inalcanzables. El teorema de mayorización nos permite estimar el comportamiento de un sistema "complejo" a través de la solución de uno "simple" con propiedades conocidas. Si el lado derecho y las condiciones iniciales de un sistema son mayores componente a componente que las de otro, entonces las trayectorias de las soluciones preservan esta desigualdad.

Para vectores $x=(x^1,\ldots,x^m)$, $y=(y^1,\ldots,y^m) \in \mathbb{R}^m$, denotamos $x \le y$ si $x^i \le y^i$ para todo $i=1,\dots,m$. El cubo $[x,y]:=\{z\in\mathbb{R}^m\mid x\le z\le y\}$.

Consideremos un dominio $D \subset \mathbb{R}^m$ tal que para $x,y\in D$, $x\le y$ implica $[x,y]\subset D$. Funciones $f_1,f_2\in C(I\times D,\mathbb{R}^m)$ con jacobianas continuas $\frac{\partial f_i}{\partial x}\in C(I\times D,\mathbb{R}^{m\times m})$, donde $I=[t_0,t_1]$.

Dos problemas de Cauchy:

$\dot x = f_i(t,x)$, $x(t_0)=\hat x_i$, $i=1,2$.

Formulación del Teorema Principal

Sea $f_1(t,x)\le f_2(t,x)$ para todo $(t,x)\in I\times D$ y $\hat x_1\le\hat x_2$. Existen matrices diagonales $U_i(t,x)$ con elementos no negativos y matrices $B_i(t,x)$ con elementos no negativos:

$\frac{\partial f_i}{\partial x}=U_i + B_i$, $i=1,2$.

Enunciado: Si $x_i(t)$ son soluciones de los problemas de Cauchy en $I$, entonces $x_1(t)\le x_2(t)$ para todo $t\in I$.

Esta propiedad nos permite mayorizar un sistema insoluble con uno simple cuyas propiedades son conocidas.

Estructura de la Demostración

Introducimos una homotopía $y=y(t,s)$ para $s\in[0,1]$:

$y_t=(1-s)f_1(t,y)+s f_2(t,y)$, $y(t_0,s)=(1-s)\hat x_1 + s\hat x_2$.

En $s=0$, obtenemos $x_1(t)$; en $s=1$, obtenemos $x_2(t)$. Basta mostrar que $s\mapsto y(t,s)$ es no decreciente para $t$ fijo, es decir, $y_s\ge 0$.

Diferenciación con Respecto al Parámetro

Diferenciamos con respecto a $s$:

$$\begin{aligned} y_{st} &= f_2 - f_1 + (u + b)y_s, \\ u(t,s) &= sU_2(t, y(t,s)) + (1-s)U_1(t, y(t,s)), \\ b(t,s) &= sB_2(t, y(t,s)) + (1-s)B_1(t, y(t,s)), \\ y_s|_{t=t_0} &= \hat x_2 - \hat x_1 \ge 0.\end{aligned}$$

Región de Monotonía

$M=A=\{t\in I\mid y_s(\xi,s)\ge0,\ \xi\in[t_0,t],\ s\in[0,1]\}$. El conjunto es no vacío ($t_0\in A$), cerrado. Sea $\tau=\sup A$.

Supongamos $\tau<t_1$. En forma integral:

$$y_s(t,s)=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)+\int_\tau^t e^{\int_\xi^t u\,d\lambda}(f_2-f_1+by_s)\,d\xi=:\Phi(y_s).$$

Método de Contracción

Por el principio de la aplicación contractiva, la ecuación $\Phi(y_s)=y_s$ se resuelve para $t-\tau>0$ pequeño. Aproximación cero $Y_0=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)\ge0$. Inductivamente $Y_{k+1}=\Phi(Y_k)\ge0$, ya que el operador $\Phi$ preserva la no negatividad. El límite $y_s\ge0$ en $[\tau,\tau+\varepsilon]\times[0,1]$, contradiciendo la definición de $\tau$.

Aplicación: Desigualdad de Gronwall-Bellman

Para $w(t)\ge0$ con $w(t)\le C+\int_{t_0}^t v(\xi)w(\xi)\,d\xi$, $v\ge0$.

Auxiliar $\phi(t)=C+\int_{t_0}^t v w\,d\xi$, $w\le\phi$, $\dot\phi=vw\le v\phi$.

Sistema 1: $\dot x=f_1(t,x):=vx-(v\phi-\dot\phi)$, $x(t_0)=C$. Los paréntesis son no negativos.

Mayorizando: $\dot x=f_2(t,x):=vx$, $x(t_0)=C$, solución $x(t)=Ce^{\int_{t_0}^tv\,d\xi}$.

Por el teorema $w\le\phi\le Ce^{\int_{t_0}^tv\,d\xi}$.

Puntos Clave:

• La mayorización funciona cuando $f_1\le f_2$ y $\hat x_1\le\hat x_2$.
• Las jacobianas se descomponen en diagonal $U$ (no negativa) y $B$ (no negativa).
• Demostración mediante homotopía y método de contracción.
• Conexión directa con la desigualdad de Gronwall para estimaciones integrales.
• Aplicación: estimación de soluciones de EDO no lineales con exponenciales simples.

Conclusión

El teorema proporciona una herramienta para el análisis cualitativo de la dinámica sin integración explícita. Las condiciones sobre las jacobianas aseguran la monotonía de la homotopía, lo cual es crítico para preservar la desigualdad.

— Editorial Team