Kurs:Analysis (Osnabrück 2014-2016)/Teil II/Vorlesung 56

Differential- und Integralgleichungen

Mit dem Begriff des Integrals einer Kurve kann man Differentialgleichungen auch als Integralgleichungen schreiben.

Lemma

Es sei ${}V$ ein endlichdimensionaler reeller Vektorraum, ${}I\subseteq \mathbb {R}$ ein reelles Intervall, ${}U\subseteq V$ eine offene Menge und

$f\colon I\times U\longrightarrow V,\,(t,v)\longmapsto f(t,v),$

ein stetiges Vektorfeld auf ${}U$. Es sei ${}(t_{0},w)\in I\times U$ vorgegeben.

Dann ist eine stetige Abbildung

$v\colon J\longrightarrow U,\,t\longmapsto v(t),$

auf einem Intervall ${}J\subseteq I$ mit ${}t_{0}\in J$ genau dann eine Lösung des Anfangswertproblems  (insbesondere muss ${}v$ differenzierbar sein)

$v'=f(t,v){\text{ und }}v(t_{0})=w,$

wenn ${}v$ die Integralgleichung

${}v(t)=w+\int _{t_{0}}^{t}f(s,v(s))\,ds\,$

erfüllt.

Beweis

Es sei die Integralbedingung erfüllt. Dann ist

${}v(t_{0})=w\,$

und aufgrund des Hauptsatzes der Infinitesimalrechnung gilt ${}v'(t)=f(t,v(t))$. Insbesondere sichert die Integralbedingung, dass ${}v$ differenzierbar ist.
Wenn umgekehrt ${}v$ eine Lösung des Anfangswertproblems ist, so ist ${}v'(s)=f(s,v(s))$ und daher

${}w+\int _{t_{0}}^{t}f(s,v(s))\,ds=w+\int _{t_{0}}^{t}v'(s)\,ds=w+v(t)-v(t_{0})=v(t)\,.$

$\Box$

Der Satz von Picard-Lindelöf

Wir kommen nun zum wichtigsten Existenz- und Eindeutigkeitssatz für die Lösungen von gewöhnlichen Differentialgleichungen.

Satz

Es sei ${}V$ ein endlichdimensionaler reeller Vektorraum, ${}I\subseteq \mathbb {R}$ ein reelles Intervall, ${}U\subseteq V$ eine offene Menge und

$f\colon I\times U\longrightarrow V,\,(t,v)\longmapsto f(t,v),$

ein Vektorfeld auf ${}U$. Es sei vorausgesetzt, dass dieses Vektorfeld stetig sei und lokal einer Lipschitz-Bedingung genüge.

Dann gibt es zu jedem ${}(t_{0},w)\in I\times U$ ein offenes Intervall ${}J$ mit ${}t_{0}\in J\subseteq I$ derart, dass auf diesem Intervall eine eindeutige Lösung für das Anfangswertproblem

$v'=f(t,v){\text{ und }}v(t_{0})=w$

existiert.

Beweis

Nach Lemma 56.1 ist eine stetige Abbildung

$v\colon J\longrightarrow V$

genau dann eine Lösung des Anfangswertproblems, wenn ${}v$ die Integralgleichung

${}v(t)=w+\int _{t_{0}}^{t}f(s,v(s))\,ds\,$

erfüllt. Wir wollen die Existenz und Eindeutigkeit einer Lösung für diese Integralgleichung unter Verwendung des Banachschen Fixpunktsatzes dadurch erweisen, dass wir für die Abbildung (man spricht von einem Funktional)

$\psi \longmapsto (t\mapsto w+\int _{t_{0}}^{t}f(s,\psi (s))\,ds)$

einen Fixpunkt finden. Hierbei stehen links und rechts Abbildungen in ${}t$ (aus einem gewissen Teilintervall von ${}I$ mit Werten in ${}V$). Die Fixpunkteigenschaft ${}H(\psi )=\psi$ bedeutet gerade, dass ${}\psi (t)=w+\int _{t_{0}}^{t}f(s,\psi (s))ds$ ist. Um den Fixpunktsatz anwenden zu können müssen wir ein Definitionsintervall festlegen, und eine Metrik auf dem Abbildungsraum nach ${}V$ definieren, diesen metrischen Raum dann als vollständig und das Funktional als stark kontrahierend nachweisen. Aufgrund der Voraussetzung über die lokale Lipschitz-Bedingung gibt es eine offene Umgebung

${}(t_{0},w)\in J'\times U{\left(w,\epsilon \right)}\subseteq I\times U\,$

und ein ${}L\in \mathbb {R} _{\geq 0}$ mit

${}\Vert {f(t,v)-f(t,{\tilde {v}})}\Vert \leq L\Vert {v-{\tilde {v}}}\Vert \,$

für alle ${}t\in J'$ und ${}v,{\tilde {v}}\in U{\left(w,\epsilon \right)}$. Durch Verkleinern der Radien können wir annehmen, dass der Abschluss von ${}J'\times U{\left(w,\epsilon \right)}$, also das Produkt des abgeschlossenen Intervalls mit der abgeschlossenen Kugel, ebenfalls in ${}I\times U$ liegt. Aufgrund von Satz 36.12 gibt es ein ${}M\in \mathbb {R} _{+}$ mit

$\Vert {f(t,v)}\Vert \leq M{\text{ für alle }}(t,v)\in J'\times U{\left(w,\epsilon \right)}$

(da diese Beschränktheit auf dem Abschluss gilt). Wir ersetzen nun ${}J'$ durch ein kleineres Intervall

${}J=[t_{0}-\delta ,t_{0}+\delta ]\subseteq J'\,$

mit ${}\delta >0$, ${}\delta \leq \epsilon /M$ und ${}\delta \leq 1/(2L)$.
Wir betrachten nun die Menge der stetigen Abbildungen

${}{\begin{aligned}C&={\left\{\psi :J\rightarrow V\mid \psi {\text{ stetig}},\,\Vert {\psi (t)-w}\Vert \leq \epsilon {\text{ für alle }}t\in J\right\}}\\&={\left\{\psi :J\rightarrow V\mid \psi {\text{ stetig}},\,\Vert {\psi -w}\Vert \leq \epsilon \right\}}.\end{aligned}}$

Dabei wird also ${}C$ mit der Maximumsnorm auf ${}J$ versehen. Dieser Raum ist nach Satz 55.9 und nach Aufgabe 36.16 wieder ein vollständiger metrischer Raum.
Wir betrachten nun auf diesem konstruierten Intervall ${}J$ bzw. der zugehörigen Menge ${}C$ die Abbildung

$H\colon C\longrightarrow C,\,\psi \longmapsto H(\psi )=(t\mapsto w+\int _{t_{0}}^{t}f(s,\psi (s))\,ds).$

Dazu müssen wir zunächst zeigen, dass ${}H(\psi )$ wieder zu ${}C$ gehört. Für ${}t\in J$ ist aber nach Satz 39.1

${}{\begin{aligned}\Vert {H(\psi )(t)-w}\Vert &=\Vert {\int _{t_{0}}^{t}f(s,\psi (s))\,ds}\Vert \\&\leq {op:Betrag\int _{t_{0}}^{t}\Vert {f(s,\psi (s))}\Vert \,ds}\\&\leq \vert {t-t_{0}}\vert M\\&\leq {\frac {\epsilon }{M}}M\\&=\epsilon ,\end{aligned}}$

und ${}H(\psi )$ ist stetig, da es durch ein Integral definiert wird.
Zum Nachweis der Kontraktionseigenschaft seien ${}\psi _{1},\psi _{2}\in C$ gegeben. Für ein ${}t\in J$ ist

${}{\begin{aligned}\Vert {H(\psi _{1})(t)-H(\psi _{2})(t)}\Vert &=\Vert {\int _{t_{0}}^{t}f(s,\psi _{1}(s))\,ds-\int _{t_{0}}^{t}f(s,\psi _{2}(s))\,ds}\Vert \\&=\Vert {\int _{t_{0}}^{t}(f(s,\psi _{1}(s))-f(s,\psi _{2}(s)))\,ds}\Vert \\&\leq \vert {\int _{t_{0}}^{t}\Vert {f(s,\psi _{1}(s))-f(s,\psi _{2}(s))}\Vert \,ds}\vert \\&\leq \vert {\int _{t_{0}}^{t}L\Vert {\psi _{1}(s)-\psi _{2}(s)}\Vert \,ds}\vert \\&=L\cdot \vert {\int _{t_{0}}^{t}\Vert {\psi _{1}(s)-\psi _{2}(s)}\Vert \,ds}\vert \\&\leq L\cdot \vert {\int _{t_{0}}^{t}\Vert {\psi _{1}-\psi _{2}}\Vert \,ds}\vert \\&\leq L\vert {t-t_{0}}\vert \cdot \Vert {\psi _{1}-\psi _{2}}\Vert \\&\leq {\frac {1}{2}}\Vert {\psi _{1}-\psi _{2}}\Vert .\end{aligned}}$

Da dies für jedes ${}t\in J$ gilt, folgt aus dieser Abschätzung direkt

${}\Vert {H(\psi _{1})-H(\psi _{2})}\Vert \leq {\frac {1}{2}}\Vert {\psi _{1}-\psi _{2}}\Vert \,,$

d.h. es liegt eine starke Kontraktion vor. Daher besitzt ${}H$ ein eindeutiges Fixelement ${}\psi \in C$, und diese Abbildung löst die Differentialgleichung. Dies gilt dann erst recht auf jedem offenen Teilintervall von ${}J$.
Damit haben wir insbesondere bewiesen, dass es in ${}C$ nur eine Lösung geben kann, wir wollen aber generell auf dem Intervall ${}J$ Eindeutigkeit erhalten. Für eine Lösung ${}v\colon J\rightarrow V$ gilt aber wegen der Integralbeziehung wieder

${}v(t)=w+\int _{t_{0}}^{t}f(s,v(s))ds\,$

und die gleichen Abschätzungen wie weiter oben zeigen, dass die Lösung zu ${}C$ gehören muss.

$\Box$

Die Picard-Lindelöf-Iteration

Der Beweis des Satzes von Picard-Lindelöf ist prinzipiell konstruktiv. Darauf beruht die Picard-Lindelöf-Iteration, mit der man Lösungen approximieren kann. Die Güte der Approximationen wird dabei durch geeignete Normen auf Funktionenräumen gemessen, was wir nicht ausführen.

Bemerkung  

Es sei ${}V$ ein endlichdimensionaler reeller Vektorraum, ${}I\subseteq \mathbb {R}$ ein reelles Intervall, ${}U\subseteq V$ eine offene Menge und

$F\colon I\times U\longrightarrow V,\,(t,v)\longmapsto F(t,v),$

ein Vektorfeld auf ${}U$. Es sei ${}(t_{0},w)\in I\times U$ eine Anfangsbedingung. Es sei vorausgesetzt, dass dieses Vektorfeld stetig sei und lokal einer Lipschitz-Bedingung genüge. In der Picard-Lindelöf-Iteration definiert man iterativ eine Folge von Funktionen

$\varphi _{n}\colon I\longrightarrow V$

durch ${}\varphi _{0}=w$ (dies ist also die konstante Funktion mit dem Wert ${}w$) und durch

${}\varphi _{n+1}(t)=w+\int _{t_{0}}^{t}F(s,\varphi _{n}(s))\,ds\,.$

Dann gibt es ein Teilintervall ${}{]a,b[}\subseteq I$ mit ${}t_{0}\in {]a,b[}$ derart, dass für ${}t\in {]a,b[}$ die Folge ${}\varphi _{n}(t)$ gegen einen Punkt ${}\varphi (t)$ konvergiert, wobei gleichmäßige Konvergenz vorliegt. Die Grenzfunktion ${}\varphi$ ist dann eine Lösung des Anfangswertproblems

$v'=f(t,v){\text{ und }}v(t_{0})=w.$

Bei einer linearen Differentialgleichung mit stetigen Koeffizientenfunktionen konvergiert dieses Verfahren auf ganz ${}I$.

Bemerkung  

Zu einem ortsunabhängigen Vektorfeld

$F\colon I\times \mathbb {R} ^{n}\longrightarrow \mathbb {R} ^{n},\,(t,v)\longmapsto F(t,v)=F(t),$

und der Anfangsbedingung ${}v(t_{0})=w$ führt die erste Picard-Lindelöf-Iteration auf

${}\varphi _{1}(t)=w+\int _{t_{0}}^{t}F(s,w)ds=w+\int _{t_{0}}^{t}F(s)ds=w+G(t)\,,$

wobei ${}G(t)$ eine Stammkurve zu ${}F(t)$ mit ${}G(t_{0})=0$ sei. Die erste Iteration liefert hier also direkt die Lösung. Die kontrahierende Abbildung im Beweis zu Satz 56.2 ist in dieser Situation konstant.

Wir wenden dieses approximative Verfahren auf eine Differentialgleichung mit getrennten Variablen an, für die wir die Lösung schon kennen (siehe Aufgabe 30.6).

Beispiel

Wir wenden die Picard-Lindelöf-Iteration auf die Differentialgleichung

${}y'=F(t,y)=ty\,$

mit der Anfangsbedingung

${}y(0)=1\,$

an (die Lösung ist $e^{{\frac {1}{2}}t^{2}}$). Daher ist ${}\varphi _{0}=1$. Die erste Iteration liefert

${}\varphi _{1}(t)=1+\int _{0}^{t}sds=1+{\frac {1}{2}}t^{2}\,.$

Die zweite Iteration liefert

${}{\begin{aligned}\varphi _{2}(t)&=1+\int _{0}^{t}F{\left(s,\varphi _{1}(s)\right)}ds\\&=1+\int _{0}^{t}F{\left(s,1+{\frac {1}{2}}s^{2}\right)}ds\\&=1+\int _{0}^{t}s+{\frac {1}{2}}s^{3}ds\\&=1+{\frac {1}{2}}t^{2}+{\frac {1}{8}}t^{4}.\end{aligned}}$

Die dritte Iteration liefert

${}{\begin{aligned}\varphi _{3}(t)&=1+\int _{0}^{t}F{\left(s,\varphi _{2}(s)\right)}ds\\&=1+\int _{0}^{t}F{\left(s,1+{\frac {1}{2}}s^{2}+{\frac {1}{8}}s^{4}\right)}ds\\&=1+\int _{0}^{t}s+{\frac {1}{2}}s^{3}+{\frac {1}{8}}s^{5}ds\\&=1+{\frac {1}{2}}t^{2}+{\frac {1}{8}}t^{4}+{\frac {1}{48}}t^{6}.\end{aligned}}$

Dabei stimmt die ${}i$-te Iteration mit der Taylor-Entwicklung der Ordnung ${}2i$ der Lösung überein.

Bemerkung  

Es sei

${}v'=Mv\,$

eine lineare Differentialgleichung mit konstanten Koeffizienten auf dem ${}\mathbb {R} ^{n}$ und es sei eine Anfangsbedingung ${}v(0)=w$ gegeben. Wir behaupten, dass die ${}n$-te Picard-Lindelöf-Iteration gleich

${}\varphi _{n}(t)=\sum _{k=0}^{n}{\frac {1}{k!}}M^{\circ k}(w)t^{k}\,$

ist, wobei ${}M^{\circ k}$ die ${}k$-fache Potenz der Matrix bezeichnet. Diese Aussage zeigen wir durch Induktion nach ${}n$. Für ${}n=0$ steht rechts einfach die konstante Kurve ${}w$. Es sei die Aussage nun für ${}n$ schon bewiesen. Dann ist

${}{\begin{aligned}\varphi _{n+1}(t)&=w+\int _{0}^{t}M(\varphi _{n}(s))ds\\&=w+\int _{0}^{t}M{\left(\sum _{k=0}^{n}{\frac {1}{k!}}M^{\circ k}(w)s^{k}\right)}ds\\&=w+\sum _{k=0}^{n}{\frac {1}{k!}}\int _{0}^{t}s^{k}M{\left(M^{\circ k}(w)\right)}ds\\&=w+\sum _{k=0}^{n}{\frac {1}{k!}}\int _{0}^{t}s^{k}M^{\circ k+1}(w)ds\\&=w+\sum _{k=0}^{n}{\frac {1}{k!}}{\frac {1}{(k+1)}}t^{k+1}M^{\circ k+1}(w)\\&=w+\sum _{k=0}^{n}{\frac {1}{(k+1)!}}t^{k+1}M^{\circ k+1}(w)\\&=w+\sum _{\ell =1}^{n+1}{\frac {1}{\ell !}}t^{\ell }M^{\circ \ell }(w)\\&=\sum _{\ell =0}^{n+1}{\frac {1}{\ell !}}t^{\ell }M^{\circ \ell }(w)\end{aligned}}$

und die Aussage ist auch für ${}n+1$ richtig. Diese Approximationen sind die Anfangsglieder in der „Exponentialreihe in dem Ausdruck“ ${}tM$. Man kann zeigen, dass diese Exponentialreihe auf ${}\mathbb {R}$ konvergiert und in der Tat die Lösung des Anfangswertproblems ist (der Satz von Picard-Lindelöf sichert nur die Konvergenz auf einer Intervallumgebung).

Beispiel

Wir wenden die Picard-Lindelöf-Iteration auf das Anfangswertproblem

${}v(0)=(0,0)\,$

zum Vektorfeld

$F\colon \mathbb {R} \times \mathbb {R} ^{2}\longrightarrow \mathbb {R} ^{2},\,(t,x,y)\longmapsto (t^{2}y-xy^{2},x^{3}-y+t)$

an. Es ist

${}\varphi _{0}(t)=(0,0)\,.$

Daher ist

${}\varphi _{1}(t)=\int _{0}^{t}(0,s)ds=\left(0,\,{\frac {1}{2}}t^{2}\right)\,.$

Es ist

${}F(s,\varphi _{1}(s))=\left(s^{2}{\frac {1}{2}}s^{2},\,-{\frac {1}{2}}s^{2}+s\right)=\left({\frac {1}{2}}s^{4},\,-{\frac {1}{2}}s^{2}+s\right)\,$

und daher

${}\varphi _{2}(t)=\int _{0}^{t}\left({\frac {1}{2}}s^{4},\,-{\frac {1}{2}}s^{2}+s\right)ds=\left({\frac {1}{10}}t^{5},\,-{\frac {1}{6}}t^{3}+{\frac {1}{2}}t^{2}\right)\,.$

Wegen

${}{\begin{aligned}F(s,\varphi _{2}(s))&=\left(s^{2}{\left(-{\frac {1}{6}}s^{3}+{\frac {1}{2}}s^{2}\right)}-{\frac {1}{10}}s^{5}{\left(-{\frac {1}{6}}s^{3}+{\frac {1}{2}}s^{2}\right)}^{2},\,{\left({\frac {1}{10}}t^{5}\right)}^{3}-{\left(-{\frac {1}{6}}s^{3}+{\frac {1}{2}}s^{2}\right)}+s\right)\\&=\left(-{\frac {1}{6}}s^{5}+{\frac {1}{2}}s^{4}-{\frac {1}{10}}s^{5}{\left({\frac {1}{36}}s^{6}-{\frac {1}{6}}s^{5}+{\frac {1}{4}}s^{4}\right)},\,{\frac {1}{1000}}s^{15}+{\frac {1}{6}}s^{3}-{\frac {1}{2}}s^{2}+s\right)\\&=\left(-{\frac {1}{360}}s^{11}+{\frac {1}{60}}s^{10}-{\frac {1}{40}}s^{9}-{\frac {1}{6}}s^{5}+{\frac {1}{2}}s^{4},\,{\frac {1}{1000}}s^{15}+{\frac {1}{6}}s^{3}-{\frac {1}{2}}s^{2}+s\right)\end{aligned}}$

und daher

${}{\begin{aligned}\varphi _{3}(t)&=\int _{0}^{t}\left(-{\frac {1}{360}}s^{11}+{\frac {1}{60}}s^{10}-{\frac {1}{40}}s^{9}-{\frac {1}{6}}s^{5}+{\frac {1}{2}}s^{4},\,{\frac {1}{1000}}s^{15}+{\frac {1}{6}}s^{3}-{\frac {1}{2}}s^{2}+s\right)ds\\&=\left(-{\frac {1}{4320}}t^{12}+{\frac {1}{660}}t^{11}-{\frac {1}{400}}t^{10}-{\frac {1}{36}}t^{6}+{\frac {1}{10}}t^{5},\,{\frac {1}{16000}}t^{16}+{\frac {1}{24}}t^{4}-{\frac {1}{6}}t^{3}+{\frac {1}{2}}t^{2}\right).\end{aligned}}$

<< | Kurs:Analysis (Osnabrück 2014-2016)/Teil II | >> PDF-Version dieser Vorlesung Arbeitsblatt zur Vorlesung (PDF)