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

Aus Wikiversity

Wechseln zu: Navigation, Suche

Inhaltsverzeichnis

[Bearbeiten] Voraussetzung (b):

Seien die Zahl \xi \in \mathbb{R} und der Vektor \eta = (\eta_1, \ldots, \eta_m) \in \mathbb{R}^m mit m \in \mathbb{N} vorgegeben. Zu den festen positiven Konstanten a, b_1, \ldots, b_m \in (0, + \infty) betrachten wir wiederum das Rechteck

R := \{(x, y) = (x, y_1, \ldots, y_m) \in \mathbb{R}^{1 + m}: |x - \xi| \le a, |y_i - \eta_i| \le b_i, i = 1, \ldots, m\}.

Hierauf sind die beschränkten, stetigen Funktionen

f_i = f_i(x, y_1, \ldots, y_m): R \to \mathbb{R} \in C^0(R) mit |f_i| \le M in R

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

(1) \begin{matrix} |f_i(x, y_1, \ldots, y_m) - f_i(x, \overline{y}_1, \ldots, \overline{y}_m)| \le L \cdot \sum\limits^m_{k = 1} |y_k - \overline{y}_k| \\ \mathrm{\ f\ddot ur\ alle\ } (x, y_1, \ldots, y_m), (x, \overline{y}_1, \ldots, \overline{y}_m) \in R \text{ und } i = 1, \ldots, m \end{matrix}

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

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

(2) \begin{matrix} y_i = y_i(x): [\xi - h, \xi + h] \to [\eta_i - b_i, \eta_i + b_i], \\ y_i'(x) = f_i(x, y_1(x), \ldots, y_m(x)), \quad x \in [\xi - h, \xi + h], \\ y_i(\xi) = \eta_i\ \mathrm{f\ddot ur}\ i = 1, \ldots, m \end{matrix}

sowie

(3) \begin{matrix} \overline{y}_i = \overline{y}_i(x): [\xi - h, \xi + h] \to [\eta_i - b_i, \eta_i + b_i], \\ \overline{y}_i'(x) = f_i(x, \overline{y}_1(x), \ldots, \overline{y}_m(x)), \quad x \in [\xi - h, \xi + h], \\ \overline{y}_i(\xi) = \overline{\eta}_i\ \mathrm{f\ddot ur}\ i = 1, \ldots, m \end{matrix}

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

(4) y_i(x) = \eta_i + \int^x_\xi f_i(t, y_1(t), \ldots, y_m(t))\, dt, \quad x \in [\xi - h, \xi + h], \quad i = 1, \ldots, m

bzw.

(5) \overline{y}_i(x) = \overline{\eta}_i + \int^x_\xi f_i(t, \overline{y}_1(t), \ldots, \overline{y}_m(t))\, dt, \quad x \in [\xi - h, \xi + h], \quad i = 1, \ldots, m.

Wir ziehen nun diese beiden Gleichungen voneinander ab und erhalten

(6) \begin{matrix} [y_i(x) - \eta_i] - [\overline{y}_i(x) - \overline{\eta}_i] = \int^x_\xi \{f_i(t, y_1(t), \ldots, y_m(t)) - f_i(t, \overline{y}_1(t), \ldots, \overline{y}_m(t))\}\, dt, \\ x \in [\xi - h, \xi + h], \quad i = 1, \ldots, m. \end{matrix}

Die Lipschitz-Bedingung liefert

(7) \begin{matrix} |[y_i(x) - \eta_i] - [\overline{y}_i(x) - \overline{\eta}_i]| \le \int^x_\xi |f_i(t, y_1(t), \ldots, y_m(t)) - f_i(t, \overline{y}_1(t), \ldots, \overline{y}_m(t))|\, |dt| \\ \le \int^x_\xi L \cdot \sum^m_{k = 1} |y_k(t) - \overline{y}_k(t)|\, |dt|\ \mathrm{f\ddot ur\ alle}\ x \in [\xi - h, \xi + h], \quad i = 1, \ldots, m. \end{matrix}

Wir führen nun die Hilfsfunktion

(8) \Phi(x) := \sum^m_{i = 1} |[y_i(x) - \eta_i] - [\overline{y}_i(x) - \overline{\eta}_i]|, \quad x \in [\xi - h, \xi + h]

sowie die Hilfsgröße

(9) \varepsilon(\eta, \overline{\eta}) := \sum^m_{k = 1} |\eta_k - \overline{\eta}_k|

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

(10) \begin{matrix} |[y_i(x) - \eta_i] - [\overline{y}_i(x) - \overline{\eta}_i]| \le \int^x_\xi L \cdot \{\Phi(t) + \varepsilon(\eta, \overline{\eta})\}\, |dt| \\ \mathrm{f\ddot ur\ alle}\ x \in [\xi - h, \xi + h], \quad i = 1, \ldots, m. \end{matrix}

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

(11) \Phi(x) \le (m \cdot L) \int^x_\xi \{\Phi(t) + \varepsilon(\eta, \overline{\eta})\}\, |dt| für alle x \in [\xi - h, \xi + h].

[Bearbeiten] Satz 1 (Gronwallsches Lemma)

Die stetige Funktion f: [\xi - h, \xi + h] \to [0, + \infty) genüge der Integralungleichung
f(x) \le A \int^x_\xi (f(t) + \varepsilon)\, |dt| für alle x \in [\xi - h, \xi + h]
mit Konstanten A > 0 und \varepsilon \ge 0. Dann gilt für alle x \in [\xi - h, \xi + h] die Abschätzung
0 \le f(x) \le \varepsilon (e^{A |x - \xi|} - 1) = \varepsilon \sum^\infty_{k = 1} \frac{A^k}{k!} |x - \xi|^k.

[Bearbeiten] Beweis

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

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

Aus der Integralungleichung erhalten wir nämlich

f(x) \le MA |x - \xi| + \varepsilon 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 \varepsilon A |x - \xi| + A \int^x_\xi f(t)\, |dt|
\le \varepsilon A |x - \xi| + A \int^x_\xi \left\{ \varepsilon \sum^n_{k = 1} \frac{A^k}{k!} |x - \xi|^k + M \frac{A^n}{n!} |x - \xi|^n \right\}\, |dt|
= \varepsilon A |x - \xi| + \varepsilon \sum^n_{k = 1} \frac{A^{k + 1}}{(k + 1)!} |x - \xi|^{k + 1} + M \frac{A^{n + 1}}{(n + 1)!} |x - \xi|^{n + 1}
= \varepsilon \sum^{n + 1}_{k = 1} \frac{A^k}{k!} |x - \xi|^k + M \frac{A^{n + 1}}{(n + 1)!} |x - \xi|^{n + 1}.

Da nun

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

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

f(x) \le \varepsilon \sum^\infty_{k = 1} \frac{A^k}{k!} |x - \xi|^k = \varepsilon (e^{A |x - \xi|} - 1).

q.e.d.

[Bearbeiten] Satz 2 (Eindeutigkeit und Stabilität)

Unter der Voraussetzung (b) lösen die Funktionen yi(x) und \overline{y}_i(x) für i = 1, \ldots, m die Anfangswertprobleme (2) bzw. (3). Dann gilt für die assoziierten Funktionen aus (8) und (9) die Ungleichung
0 \le \Phi(x) \le \varepsilon(\eta, \overline{\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 \left\{ a, \frac{b_1}{M}, \ldots, \frac{b_m}{M} \right\}

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) := \begin{pmatrix} f_1(x, y_1, \ldots, y_m) \\ \vdots \\ f_m(x, y_1, \ldots, y_m) \end{pmatrix} für (x, y) \in R

und konstruieren die Funktionenfolge

(13) \stackrel{k}{y}(x) := \begin{pmatrix} \stackrel{k}{y}_1(x) \\ \vdots \\ \stackrel{k}{y}_m(x) \end{pmatrix}, \quad x \in [\xi - h, \xi + h] =: I, \quad k = 0, 1, 2, \ldots

wie folgt:

(14) \stackrel{0}{y} := \eta = \begin{pmatrix} \eta_1 \\ \vdots \\ \eta_m \end{pmatrix} und \stackrel{k}{y} := \eta + \int^x_\xi f \left(t, \stackrel{k - 1}{y}(t)\right)\, dt für k = 1, 2, \ldots.

[Bearbeiten] Hilfssatz 1

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

[Bearbeiten] Beweis

\mathbf{k = 0:} \stackrel{0}{y} = \eta.
\mathbf{k - 1 \to k:}

\left|\stackrel{k}{y} - \eta_i \right| = \left| \int^x_\xi f_i \left( t, \stackrel{k - 1}{y}(t) \right)\, dt \right| \le \int^x_\xi \left| f_i \left( t, \stackrel{k - 1}{y}(t) \right) \right|\, |dt| \le M |x - \xi| \le Mh \le b_i für alle x \in I und i = 1, \ldots, m.

Somit folgt die Behauptung.

q.e.d.

Wir betrachten nun die Funktionenfolge

\stackrel{k}{\psi}(x) := \sum^m_{i = 1} \left| \stackrel{k + 1}{y}_i(x) - \stackrel{k}{y}_i(x) \right|, \quad x \in I für k = 0, 1, 2, \ldots.

[Bearbeiten] 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, \ldots.

[Bearbeiten] Beweis

\mathbf{k = 0:} Wir sehen \left| \stackrel{1}{y}_i(x) - \stackrel{0}{y}_i(x) \right| = \left| \int^x_\xi f_i(t, \eta)\, dt \right| \le M |x - \xi| für i = 1, \ldots, m ein und erhalten \stackrel{0}{\psi}(x) \le mM |x - \xi| für alle x \in I.
\mathbf{k \to k + 1:} Wir ermitteln

\left| \stackrel{k + 2}{y}_i(x) - \stackrel{k + 1}{y}_i(x) \right| = \left| \int^x_\xi \left[ f_i \left( t, \stackrel{k + 1}{y}(t) \right) - f_i \left( t, \stackrel{k + 1}{y}(t) \right) \right]\, dt \right| \le L \int^x_\xi \stackrel{k}{\psi}(t)\, |dt|
\le (mM) (mL)^k L \frac{1}{(k + 1)!} \int^x_\xi |t - \xi|^{k + 1}\, |dt| = (mM) (mL)^k L \frac{1}{(k + 2)!} |x - \xi|^{k + 2}

für i = 1, \ldots, 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^\infty_{k = 0} \left( \stackrel{k + 1}{y}(x) - \stackrel{k}{y}(x) \right), \quad x \in I

hat somit die konvergente Majorante

\frac{M}{L} \sum^\infty_{k = 0} \frac{(mL)^k |x - \xi|^k}{k!} = \frac{M}{L} \Bigl(\exp(mL |x - \xi|) - 1 \Bigr)

[Bearbeiten] 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 \mathbb{R}^m für k = 0, 1, 2, \ldots
gleichmäßig auf dem Intervall I = [ξ − h,ξ + h] gegen eine Lösung
y(x) = \begin{pmatrix} y_1(x) \\ \vdots \\ y_m(x) \end{pmatrix} \in C^1(I, \mathbb{R}^m)
des Anfangswertproblems (2).

[Bearbeiten] Beweis

Es gilt

y(x) := \lim_{k \to \infty} \stackrel{k + 1}{y}(x) = \eta + \lim_{k \to \infty} \sum^k_{l = 1} \left(\stackrel{l + 1}{y}(x) - \stackrel{l}{y}(x)\right), \quad x \in I.

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

y(x) = \eta + \int^x_\xi f(t, y(t))\, dt, \quad x \in I.

q.e.d.

Von „http://de.wikiversity.org/wiki/Kurs:Analysis_II/Kapitel_VI:_Gew%C3%B6hnliche_Differentialgleichungen/Eindeutigkeit_und_sukzessive_Approximation_(%C2%A75)“
Persönliche Werkzeuge