Kurs:Analysis II/Kapitel VI: Gewöhnliche Differentialgleichungen/Eindeutigkeit und sukzessive Approximation (§5)

Voraussetzung (b):

Seien die Zahl $\xi \in ℝ$ und der Vektor $\eta =({\eta }_1,\dots ,{\eta }_m)\in {ℝ}^m$ mit $m\in \mathbb{N}$ vorgegeben. Zu den festen positiven Konstanten $a,b_1,\dots ,b_m\in (0,+\mathrm{\infty })$ betrachten wir wiederum das Rechteck

$R:=\{(x,y)=(x,y_1,\dots ,y_m)\in {ℝ}^{1+m}:|x-\xi |\le a,|y_i-{\eta }_i|\le b_i,i=1,\dots ,m\}$.

Hierauf sind die beschränkten, stetigen Funktionen

$f_i=f_i(x,y_1,\dots ,y_m):R\to ℝ\in C^0(R)$ mit $|f_i|\le M$ in $R$

für $i=1,\dots ,m$ mit der Schranke $M>0$ gegeben. Weiter gebe es eine Lipschitz-Konstante $L>0$, so dass die Ungleichung

(1) $\begin{array}{c}|f_i(x,y_1,\dots ,y_m)-f_i(x,{\bar{y}}_1,\dots ,{\bar{y}}_m)|\le L\cdot \sum _{k=1}^m|y_k-{\bar{y}}_k|\\ \mathrm{f}\ddot{\mathrm{u}}\mathrm{r}\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{e}(x,y_1,\dots ,y_m),(x,{\bar{y}}_1,\dots ,{\bar{y}}_m)\in R\text{ und }i=1,\dots ,m\end{array}$

erfüllt ist. Wir fordern also, dass die Funktionen $f_i$ in den Variablen $y_1,\dots ,y_m$ einer Lipschitz-Bedingung genügen.

Auf dem Existenzintervall $[\xi -h,\xi +h]$ mit hinreichend kleinem $0<h\le a$ betrachten wir die Lösungen der Anfangswertprobleme

(2) $\begin{array}{c}{y}_i=y_i(x):[\xi -h,\xi +h]\to [{\eta }_i-b_i,{\eta }_i+b_i],\\ {y_i}^{\prime }(x)=f_i(x,y_1(x),\dots ,y_m(x)),x\in [\xi -h,\xi +h],\\ y_i(\xi )={\eta }_i\mathrm{f}\ddot{\mathrm{u}}\mathrm{r}i=1,\dots ,m\end{array}$

sowie

(3) $\begin{array}{c}{\bar{y}}_i={\bar{y}}_i(x):[\xi -h,\xi +h]\to [{\eta }_i-b_i,{\eta }_i+b_i],\\ {{\bar{y}}_i}^{\prime }(x)=f_i(x,{\bar{y}}_1(x),\dots ,{\bar{y}}_m(x)),x\in [\xi -h,\xi +h],\\ {\bar{y}}_i(\xi )={\bar{\eta }}_i\mathrm{f}\ddot{\mathrm{u}}\mathrm{r}i=1,\dots ,m\end{array}$

zu den Anfangswerten $\eta$ bzw. $\bar{\eta }$. Um diese beiden Lösungen $y(x)=(y_1(x),\dots ,y_m(x))$ und $\bar{y}(x)=({\bar{y}}_1(x),\dots ,{\bar{y}}_m(x))$ mit einander zu vergleichen, betrachten wir die äquivalenten Integralgleichungssysteme

(4) $y_i(x)={\eta }_i+\underset{\xi }{\overset{x}{\int }}{f}_i(t,y_1(t),\dots ,y_m(t))dt,x\in [\xi -h,\xi +h],i=1,\dots ,m$

bzw.

(5) ${\bar{y}}_i(x)={\bar{\eta }}_i+\underset{\xi }{\overset{x}{\int }}{f}_i(t,{\bar{y}}_1(t),\dots ,{\bar{y}}_m(t))dt,x\in [\xi -h,\xi +h],i=1,\dots ,m$.

Wir ziehen nun diese beiden Gleichungen voneinander ab und erhalten

(6) $\begin{array}{c}[y_i(x)-{\eta }_i]-[{\bar{y}}_i(x)-{\bar{\eta }}_i]=\underset{\xi }{\overset{x}{\int }}\{f_i(t,y_1(t),\dots ,y_m(t))-f_i(t,{\bar{y}}_1(t),\dots ,{\bar{y}}_m(t))\}dt,\\ x\in [\xi -h,\xi +h],i=1,\dots ,m.\end{array}$

Die Lipschitz-Bedingung liefert

(7) $\begin{array}{c}|[y_i(x)-{\eta }_i]-[{\bar{y}}_i(x)-{\bar{\eta }}_i]|\le \underset{\xi }{\overset{x}{\int }}|f_i(t,y_1(t),\dots ,y_m(t))-f_i(t,{\bar{y}}_1(t),\dots ,{\bar{y}}_m(t))||dt|\\ \le \underset{\xi }{\overset{x}{\int }}L\cdot \sum _{k=1}^m|y_k(t)-{\bar{y}}_k(t)||dt|\mathrm{f}\ddot{\mathrm{u}}\mathrm{r}\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{e}x\in [\xi -h,\xi +h],i=1,\dots ,m.\end{array}$

Wir führen nun die Hilfsfunktion

(8) $\mathrm{\Phi }(x):=\sum _{i=1}^m|[y_i(x)-{\eta }_i]-[{\bar{y}}_i(x)-{\bar{\eta }}_i]|,x\in [\xi -h,\xi +h]$

sowie die Hilfsgröße

(9) $\epsilon (\eta ,\bar{\eta }):=\sum _{k=1}^m|{\eta }_k-{\bar{\eta }}_k|$

ein. Wir entnehmen dann (7) die Abschätzung

(10) $\begin{array}{c}|[y_i(x)-{\eta }_i]-[{\bar{y}}_i(x)-{\bar{\eta }}_i]|\le \underset{\xi }{\overset{x}{\int }}L\cdot \{\mathrm{\Phi }(t)+\epsilon (\eta ,\bar{\eta })\}|dt|\\ \mathrm{f}\ddot{\mathrm{u}}\mathrm{r}\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{e}x\in [\xi -h,\xi +h],i=1,\dots ,m.\end{array}$

Summation von $1,\dots ,m$ liefert schließlich die Ungleichung

(11) $\mathrm{\Phi }(x)\le (m\cdot L)\underset{\xi }{\overset{x}{\int }}\{\mathrm{\Phi }(t)+\epsilon (\eta ,\bar{\eta })\}|dt|$ für alle $x\in [\xi -h,\xi +h]$.

Satz 1 (Gronwallsches Lemma)

Die stetige Funktion $f:[\xi -h,\xi +h]\to [0,+\mathrm{\infty })$ genüge der Integralungleichung

$f(x)\le A\underset{\xi }{\overset{x}{\int }}(f(t)+\epsilon )|dt|$ für alle $x\in [\xi -h,\xi +h]$

mit Konstanten $A>0$ und $\epsilon \ge 0$. Dann gilt für alle $x\in [\xi -h,\xi +h]$ die Abschätzung

$0\le f(x)\le \epsilon (e^{A|x-\xi |}-1)=\epsilon \sum _{k=1}^{\mathrm{\infty }}\frac{A^k}{k!}|x-\xi {|}^k$.

Beweis

Wir setzen $M:=\max \{f(x):\xi -h\le x\le \xi +h\}$ und zeigen durch vollständige Induktion

$f(x)\le \epsilon \sum _{k=1}^n\frac{A^k}{k!}|x-\xi {|}^k+M\frac{A^n}{n!}|x-\xi {|}^n,x\in [\xi -h,\xi +h]$.

Aus der Integralungleichung erhalten wir nämlich

$f(x)\le MA|x-\xi |+\epsilon A|x-\xi |$ für alle $x\in [\xi -h,\xi +h]$,

so dass der Fall $n=1$ gesichert ist. Gilt nun obige Abschätzung für ein $n\in \mathbb{N}$, so finden wir

$f(x)\le \epsilon A|x-\xi |+A\underset{\xi }{\overset{x}{\int }}f(t)|dt|$

$\le \epsilon A|x-\xi |+A\underset{\xi }{\overset{x}{\int }}\{\epsilon \sum _{k=1}^n\frac{A^k}{k!}|x-\xi {|}^k+M\frac{A^n}{n!}|x-\xi {|}^n\}|dt|$

$=\epsilon A|x-\xi |+\epsilon \sum _{k=1}^n\frac{A^{k+1}}{(k+1)!}|x-\xi {|}^{k+1}+M\frac{A^{n+1}}{(n+1)!}|x-\xi {|}^{n+1}$

$=\epsilon \sum _{k=1}^{n+1}\frac{A^k}{k!}|x-\xi {|}^k+M\frac{A^{n+1}}{(n+1)!}|x-\xi {|}^{n+1}$.

Da nun

$\underset{n\to \mathrm{\infty }}{\lim }\frac{(A|x-\xi |{)}^{n+1}}{(n+1)!}=0$

richtig ist, folgt durch Grenzübergang in obiger Abschätzung

$f(x)\le \epsilon \sum _{k=1}^{\mathrm{\infty }}\frac{A^k}{k!}|x-\xi {|}^k=\epsilon (e^{A|x-\xi |}-1)$.

q.e.d.

Satz 2 (Eindeutigkeit und Stabilität)

Unter der Voraussetzung (b) lösen die Funktionen $y_i(x)$ und ${\bar{y}}_i(x)$ für $i=1,\dots ,m$ die Anfangswertprobleme (2) bzw. (3). Dann gilt für die assoziierten Funktionen aus (8) und (9) die Ungleichung

$0\le \mathrm{\Phi }(x)\le \epsilon (\eta ,\bar{\eta })\cdot \{\exp (mL|x-\xi |)-1\}$ für alle $x\in [\xi -h,\xi +h]$.

Somit stimmen die Lösungen bei gleichen Anfangswerten überein und sie hängen überdies stetig von diesen Anfangswerten ab.

Unter der Voraussetzung (b) erklären wir die Größe

$h:=\min \{a,\frac{b_1}{M},\dots ,\frac{b_m}{M}\}$

und erhalten nach dem Peanoschen Existenzsatz genau eine Lösung des Anfangswertproblems (2). Unter Verwendung der Lipschitzbedingung werden wir diese mit dem Verfahren der sukzessiven Approximation von Picard und Lindelöf auf ganz anderem Wege konstruieren. Wir setzen

(12) $f(x,y):=\left(\begin{array}{c}{f}_1(x,y_1,\dots ,y_m)\\ ⋮\\ f_m(x,y_1,\dots ,y_m)\end{array}\right)$ für $(x,y)\in R$

und konstruieren die Funktionenfolge

(13) $\stackrel{k}{y}(x):=\left(\begin{array}{c}{\stackrel{k}{y}}_1(x)\\ ⋮\\ {\stackrel{k}{y}}_m(x)\end{array}\right),x\in [\xi -h,\xi +h]=:I,k=0,1,2,\dots$

wie folgt:

(14) $\stackrel{0}{y}:=\eta =\left(\begin{array}{c}{\eta }_1\\ ⋮\\ {\eta }_m\end{array}\right)$ und $\stackrel{k}{y}:=\eta +\underset{\xi }{\overset{x}{\int }}f(t,\stackrel{k-1}{y}(t))dt$ für $k=1,2,\dots$.

Hilfssatz 1

Es gilt $(x,\stackrel{k}{y}(x))\in R$ für alle $x\in I$ und $k=0,1,2,\dots$.

Beweis

$𝐤\mathbf{=}\mathrm{𝟎}\mathbf{:}$ $\stackrel{0}{y}=\eta .$
$𝐤\mathbf{-}\mathrm{𝟏}\mathbf{\to }𝐤\mathbf{:}$

$|\stackrel{k}{y}-{\eta }_i|=\left|\underset{\xi }{\overset{x}{\int }}{f}_i(t,\stackrel{k-1}{y}(t))dt\right|\le \underset{\xi }{\overset{x}{\int }}\left|f_i(t,\stackrel{k-1}{y}(t))\right||dt|\le M|x-\xi |\le Mh\le b_i$ für alle $x\in I$ und $i=1,\dots ,m$.

Somit folgt die Behauptung.

q.e.d.

Wir betrachten nun die Funktionenfolge

$\stackrel{k}{\psi }(x):=\sum _{i=1}^m|{\stackrel{k+1}{y}}_i(x)-{\stackrel{k}{y}}_i(x)|,x\in I$ für $k=0,1,2,\dots$.

Hilfssatz 2

Es gilt

$\stackrel{k}{\psi }(x)\le mM\frac{(mL{)}^k|x-\xi {|}^{k+1}}{(k+1)!}$

für alle $x\in I$ und $k=0,1,2,\dots$.

Beweis

$𝐤\mathbf{=}\mathrm{𝟎}\mathbf{:}$ Wir sehen $|{\stackrel{1}{y}}_i(x)-{\stackrel{0}{y}}_i(x)|=|\underset{\xi }{\overset{x}{\int }}{f}_i(t,\eta )dt|\le M|x-\xi |$ für $i=1,\dots ,m$ ein und erhalten $\stackrel{0}{\psi }(x)\le mM|x-\xi |$ für alle $x\in I$.
$𝐤\mathbf{\to }𝐤\mathbf{+}\mathrm{𝟏}\mathbf{:}$ Wir ermitteln

$|{\stackrel{k+2}{y}}_i(x)-{\stackrel{k+1}{y}}_i(x)|=\left|\underset{\xi }{\overset{x}{\int }}[f_i(t,\stackrel{k+1}{y}(t))-f_i(t,\stackrel{k+1}{y}(t))]dt\right|\le L\underset{\xi }{\overset{x}{\int }}\stackrel{k}{\psi }(t)|dt|$

$\le (mM)(mL{)}^{k}L\frac{1}{(k+1)!}\underset{\xi }{\overset{x}{\int }}|t-\xi {|}^{k+1}|dt|=(mM)(mL{)}^{k}L\frac{1}{(k+2)!}|x-\xi {|}^{k+2}$

für $i=1,\dots ,m$. Wir erhalten dann

$\stackrel{k+1}{\psi }(x)\le mM\frac{(mL{)}^{k+1}|x-\xi {|}^{k+2}}{(k+2)!}$ für alle $x\in I$.

q.e.d.

Die Funktionenreihe

$\sum _{k=0}^{\mathrm{\infty }}(\stackrel{k+1}{y}(x)-\stackrel{k}{y}(x)),x\in I$

hat somit die konvergente Majorante

$\frac{M}{L}\sum _{k=0}^{\mathrm{\infty }}\frac{(mL{)}^k|x-\xi {|}^k}{k!}=\frac{M}{L}(\exp (mL|x-\xi |)-1)$

Satz 3 (Sukzessive Approximation nach Picard und Lindelöf)

Unter der Voraussetzung (b) konvergiert die in (14) definierte Funktionenfolge

$\stackrel{k}{y}(x):I\to {ℝ}^m$ für $k=0,1,2,\dots$

gleichmäßig auf dem Intervall $I=[\xi -h,\xi +h]$ gegen eine Lösung

$y(x)=\left(\begin{array}{c}{y}_1(x)\\ ⋮\\ y_m(x)\end{array}\right)\in C^1(I,{ℝ}^m)$

des Anfangswertproblems (2).

Beweis

Es gilt

$y(x):=\underset{k\to \mathrm{\infty }}{\lim }\stackrel{k+1}{y}(x)=\eta +\underset{k\to \mathrm{\infty }}{\lim }\sum _{l=1}^k(\stackrel{l+1}{y}(x)-\stackrel{l}{y}(x)),x\in I$.

Wir vollziehen den Grenzübergang in der Integralgleichung (14) und erhalten

$y(x)=\eta +\underset{\xi }{\overset{x}{\int }}f(t,y(t))dt,x\in I$.

q.e.d.