Majorisation von ODEs: Vergleich der Lösungen

Majorisierung von ODE-Systemen: Lösungsvergleich über rechte Seiten

In der Analyse dynamischer Systeme sind exakte Lösungen von gewöhnlichen Differentialgleichungen oft nicht erreichbar. Der Majorisierungssatz erlaubt es uns, das Verhalten eines "komplexen" Systems durch die Lösung eines "einfachen" Systems mit bekannten Eigenschaften abzuschätzen. Wenn die rechte Seite und die Anfangsbedingungen eines Systems komponentenweise größer sind als die eines anderen, dann bewahren die Lösungstrajektorien diese Ungleichung.

Für Vektoren $x=(x^1,\ldots,x^m)$, $y=(y^1,\ldots,y^m) \in \mathbb{R}^m$ schreiben wir $x \le y$, wenn $x^i \le y^i$ für alle $i=1,\dots,m$ gilt. Der Würfel $[x,y]:=\{z\in\mathbb{R}^m\mid x\le z\le y\}$.

Betrachten Sie ein Gebiet $D \subset \mathbb{R}^m$, sodass für $x,y\in D$, $x\le y$ impliziert $[x,y]\subset D$. Funktionen $f_1,f_2\in C(I\times D,\mathbb{R}^m)$ mit stetigen Jacobi-Matrizen $\frac{\partial f_i}{\partial x}\in C(I\times D,\mathbb{R}^{m\times m})$, wobei $I=[t_0,t_1]$.

Zwei Cauchy-Probleme:

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

Formulierung des Hauptsatzes

Sei $f_1(t,x)\le f_2(t,x)$ für alle $(t,x)\in I\times D$ und $\hat x_1\le\hat x_2$. Es existieren Diagonalmatrizen $U_i(t,x)$ mit nicht-negativen Elementen und Matrizen $B_i(t,x)$ mit nicht-negativen Elementen:

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

Aussage: Wenn $x_i(t)$ Lösungen der Cauchy-Probleme auf $I$ sind, dann gilt $x_1(t)\le x_2(t)$ für alle $t\in I$.

Diese Eigenschaft erlaubt es uns, ein unlösbares System mit einem einfachen System, dessen Eigenschaften bekannt sind, zu majorisieren.

Beweisstruktur

Einführung einer Homotopie $y=y(t,s)$ für $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$.

Bei $s=0$ erhalten wir $x_1(t)$; bei $s=1$ erhalten wir $x_2(t)$. Es genügt zu zeigen, dass $s\mapsto y(t,s)$ für festes $t$ nicht-fallend ist, d.h. $y_s\ge 0$.

Differentiation nach Parameter

Differenzieren nach $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}$$

Monotoniebereich

$M=A=\{t\in I\mid y_s(\xi,s)\ge0,\ \xi\in[t_0,t],\ s\in[0,1]\}$. Die Menge ist nicht-leer ($t_0\in A$), abgeschlossen. Sei $\tau=\sup A$.

Angenommen $\tau<t_1$. In Integralform:

$$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).$$

Kontraktionsmethode

Nach dem Kontraktionsprinzip wird die Gleichung $\Phi(y_s)=y_s$ für kleine $t-\tau>0$ gelöst. Nullnäherung $Y_0=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)\ge0$. Induktiv $Y_{k+1}=\Phi(Y_k)\ge0$, da der Operator $\Phi$ Nicht-Negativität erhält. Der Grenzwert $y_s\ge0$ auf $[\tau,\tau+\varepsilon]\times[0,1]$, im Widerspruch zur Definition von $\tau$.

Anwendung: Gronwall-Bellman-Ungleichung

Für nicht-negative $w(t)\le C+\int_{t_0}^t v(\xi)w(\xi)\,d\xi$, $v\ge0$.

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

System 1: $\dot x=f_1(t,x):=vx-(v\phi-\dot\phi)$, $x(t_0)=C$. Die Klammer ist nicht-negativ.

Majorisierend: $\dot x=f_2(t,x):=vx$, $x(t_0)=C$, Lösung $x(t)=Ce^{\int_{t_0}^tv\,d\xi}$.

Nach dem Satz $w\le\phi\le Ce^{\int_{t_0}^tv\,d\xi}$.

Wichtige Punkte:

• Majorisierung funktioniert, wenn $f_1\le f_2$ und $\hat x_1\le\hat x_2$.
• Jacobi-Matrizen zerlegen in diagonale $U$ (nicht-negativ) und $B$ (nicht-negativ).
• Beweis über Homotopie und Kontraktionsmethode.
• Direkter Zusammenhang zur Gronwall-Ungleichung für Integralabschätzungen.
• Anwendung: Abschätzung von Lösungen nichtlinearer ODEs mit einfachen Exponentialfunktionen.

Fazit

Der Satz bietet ein Werkzeug für qualitative Analyse der Dynamik ohne explizite Integration. Bedingungen an die Jacobi-Matrizen sichern die Monotonie der Homotopie, was entscheidend für die Erhaltung der Ungleichung ist.

— Editorial Team