Anfangswertproblem

Als Anfangswertproblem (abgekürzt AWP), manchmal auch Anfangswertaufgabe (abgekürzt AWA) oder Cauchy-Problem genannt, bezeichnet man in der Analysis eine wichtige Klasse von Differentialgleichungen. Die Lösung eines Anfangswertproblems ist die Lösung der Differentialgleichung unter zusätzlicher Berücksichtigung eines vorgegebenen Anfangswertes.

In diesem Artikel wird das Anfangswertproblem zunächst für gewöhnliche Differentialgleichungen und später auch für partielle Differentialgleichungen erklärt.

Gewöhnliche Differentialgleichungen

Anfangswertproblem 1. Ordnung

Ein Anfangswertproblem erster Ordnung ist ein Gleichungssystem, das aus einer gewöhnlichen Differentialgleichung erster Ordnung

und einer zusätzlichen Anfangsbedingung

$y(t_{0})=y_{0}$

besteht, mit

dem Anfangswert $y_{0}$ und
einem Zeitpunkt $t_{0}\in {\mathbb {R}}$ .

Eine konkrete Funktion ist eine Lösung des Anfangswertproblems, wenn sie beide Gleichungen erfüllt.

Gesucht ist also eine Funktion , die die Bedingungen der Differentialgleichung und des Anfangswertes erfüllt. Ist die Funktion stetig, so ist dies nach dem Hauptsatz der Integralrechnung genau dann der Fall, wenn

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

für alle im Definitionsintervall gilt.

Anfangswertproblem k-ter Ordnung

Gegeben seien $k\in\N$ und eine Funktion $f\colon D\rightarrow \mathbb{R} ^{n}$ . Ihr Definitionsbereich sei hierbei eine Teilmenge von $I\times \mathbb{R} ^{{n\times k}}$ , worin $I\subset {\mathbb {R}}$ ein Intervall bezeichnet, welches $t_{0}$ umfasst. Dann heißt

${\begin{cases}y^{{(k)}}&=f(t,y(t),y'(t),\dotsc ,y^{{(k-1)}}(t))\\y^{{(i)}}(t_{0})&=y_{i}\qquad {\mathrm {mit}}\;i=0,\dotsc ,k-1\end{cases}}$

ein Anfangswertproblem -ter Ordnung. Jedes Anfangswertproblem -ter Ordnung lässt sich umschreiben in ein Anfangswertproblem 1. Ordnung.

Ein spezielles Anfangswertproblem ist das Riemann-Problem, bei dem die Anfangsdaten konstant sind bis auf eine Unstetigkeitsstelle.

Anfangswertprobleme treten z.B. in den Naturwissenschaften auf, wenn ein mathematisches Modell für natürliche Prozesse gesucht wird.

Lösbarkeit

Wichtige Sätze, die die Lösbarkeit von Anfangswertproblemen für gewöhnliche Differentialgleichungen betreffen, sind der (lokale) Existenzsatz von Peano und der Existenz- und Eindeutigkeitssatz von Picard-Lindelöf. Ein Hilfsmittel ist die grönwallsche Ungleichung.

Beispiel

Das Anfangswertproblem

$y'(t)=2\cdot \operatorname {sgn}(y(t))\cdot {\sqrt {|y(t)|}}\ ,\ y(0)=0\ ,$

welches zu

$f(t,x):=2\cdot \operatorname {sgn}(x)\cdot {\sqrt {|x|}}$

korrespondiert, hat unendlich viele Lösungen, nämlich neben der trivialen Lösung

$y(t)\equiv 0$

auch noch für jedes $c\geq 0$ die Lösungen

$y(t)=\left\{{\begin{array}{ll}0\ ,&{\text{falls}}\ t<c\ ,\\(t-c)^{2}\ ,&{\text{falls}}\ t\geq c\ ,\\\end{array}}\right.$

sowie

$y(t)=\left\{{\begin{array}{ll}0\ ,&{\text{falls}}\ t<c\ ,\\-(t-c)^{2}\ ,&{\text{falls}}\ t\geq c\ .\\\end{array}}\right.$

Damit Anfangswertprobleme eindeutige Lösungen besitzen, sind Zusatzeigenschaften (an ) nachzuweisen. Dies kann beispielsweise über den Satz von Picard-Lindelöf geschehen, dessen Voraussetzungen in diesem Beispiel jedoch nicht erfüllt werden.

Numerische Lösungsmethoden

Zur numerischen Lösung von Anfangswertproblemen werden Einschritt- oder Mehrschrittverfahren eingesetzt. Dabei wird die Differentialgleichung mittels einer Diskretisierung approximiert.

Partielle Differentialgleichungen

Verallgemeinert man das Cauchy-Problem auf mehrere Veränderliche, etwa Veränderliche $x_{1},\dotsc ,x_{n}$ , so erhält man partielle Differentialgleichungen. Im Folgenden stehe $\alpha \in \N_0^n$ für einen Multiindex der Länge . Beachte, dass es genau ${\tbinom {n+k-1}{k}}$ Multiindizes mit $|\alpha |:=\alpha _{1}+\dotsb +\alpha _{n}\leq k$ gibt. Es sei weiter eine Funktion in $n+{\tbinom {n+k-1}{k}}$ Variablen gegeben. Beim allgemeinen Cauchy-Problem sucht man nach Funktionen , die von Variablen $x_{1},\dotsc ,x_{n}$ abhängen und die Gleichung

(1) $F(x,(\partial ^{\alpha }u(x))_{{|\alpha |\leq k)}})\,=\,0$

erfüllen. Beachte, dass die Stelligkeit von gerade so gewählt wurde, dass man $x=(x_{1},\dotsc ,x_{n})$ und alle partiellen Ableitungen $\partial ^{\alpha }u(x)$ einsetzen kann. Darüber hinaus fordert man, dass die gesuchten Funktionen den im Folgenden beschriebenen sogenannten Anfangs- bzw. Randbedingungen genügen. Zu deren Formulierung sei eine Hyperfläche der Klasse C^k mit Normalenfeld $\nu$ . Mit $\partial _{\nu }^{j}$ seien die Normalenableitungen bezeichnet. Sind dann $\varphi _{0},\dotsc ,\varphi _{{k-1}}$ vorgegebene auf definierte Funktionen, so fordert man beim allgemeinen Cauchy-Problem, dass die Funktionen zusätzlich die Bedingungen

(2) $u=\varphi _{0},\,\partial _{\nu }u=\varphi _{1},\,\dotsc ,\,\partial _{\nu }^{{k-1}}u=\varphi _{{k-1}}$ auf

erfüllen. Die Funktionen $\varphi _{j}$ heißen die Cauchy-Daten des Problems, jede Funktion , die beide Bedingungen (1) und (2) erfüllt, heißt eine Lösung des Cauchy-Problems.

Durch eine geeignete Koordinatentransformation kann man sich auf den Fall $S=\{x=(x_{1},\dotsc ,x_{n});\,x_{n}=0\}$ zurückziehen. Dann spielt die letzte Variable eine Sonderrolle, denn die Anfangsbedingungen sind dort gegeben, wo diese Variable 0 ist. Da diese Variable in vielen Anwendungen als Zeit interpretiert wird, benennt man sie gern in (lateinisch tempus = Zeit) um, die Anfangsbedingungen beschreiben dann die Verhältnisse zum Zeitpunkt t=0 . Die Variablen sind also $x_{1},\dotsc ,x_{{n-1}},t$ . Da die betrachtete Hyperebene durch die Bedingung t=0 gegeben ist, wird die Normalenableitung einfach zur Ableitung nach . Schreibt man abkürzend $x=(x_{1},\dotsc ,x_{{n-1}})$ und $\alpha =(\alpha _{1},\dotsc ,\alpha _{{n-1}})$ , so lautet das Cauchy-Problem nun

(1') $F(x,t,(\partial _{x}^{\alpha }\partial _{t}^{j}u(x,t))_{{|\alpha |+j\leq k)}})\,=\,0$

(2') $u(x,0)=\varphi _{0}(x),\,\partial _{t}u(x,0)=\varphi _{1}(x),\,\dotsc ,\,\partial _{t}^{k-1}u(x,0)=\varphi _{k-1}(x)$ .

Ein typisches Beispiel ist etwa die dreidimensionale Wellengleichung

$\partial _{t}^{2}u-c^{2}\cdot \Delta u=f$

$u(x,0)=\varphi _{0}(x),\,\partial _{t}u(x,0)=\varphi _{1}(x)$ ,

wobei eine Konstante, eine vorgegebene Funktion und $\Delta =\partial _{{x_{1}}}^{2}+\partial _{{x_{2}}}^{2}+\partial _{{x_{3}}}^{2}$ der Laplace-Operator seien.

Ist eine Lösung, was gleichzeitig ausreichende Differenzierbarkeit implizieren soll, so sind alle Ableitungen $\partial _{x}^{\alpha }\partial _{t}^{j}u(x,0)$ mit $|\alpha |+j\leq k,j<k$ bereits durch die Cauchy-Daten vorgegeben, denn es ist $\partial _{x}^{\alpha }\partial _{t}^{j}u(x,0)=\partial _{x}^{\alpha }\varphi _{j}$ . Lediglich die Ableitung $\partial _{t}^{k}u$ ist nicht durch (2') festgelegt, hier kann also nur (1') eine Bedingung stellen. Damit (1') tatsächlich eine nicht-triviale Bedingung und damit das Cauchy-Problem nicht von vornherein schlecht gestellt ist, wird man fordern, dass man die Gleichung (1') nach $\partial _{t}^{k}u$ auflösen kann. Das Cauchy-Problem hat dann die Form

(1") $\partial _{t}^{k}u(x,t)=G(x,t,(\partial _{x}^{\alpha }\partial _{t}^{j}u(x,t))_{{|\alpha |+j\leq k,j<k}})$

(2") $u(x,0)=\varphi _{0}(x),\,\partial _{t}u(x,0)=\varphi _{1}(x),\,\dotsc ,\,\partial _{t}^{k-1}u(x,0)=\varphi _{k-1}(x)$ ,

wobei eine geeignete Funktion der Stelligkeit $n-1+{\tbinom {n+k-1}{k}}$ sei. In der zuletzt gegebenen Formulierung haben alle auftretenden Ableitungen eine Ordnung $\leq k$ , und die -te Ableitung nach tritt tatsächlich auf, denn dies ist gerade die linke Seite von (1") und sie kommt nicht auf der rechten Seite von (1") vor. Man nennt daher auch die Ordnung des Cauchy-Problems. Das obige Beispiel der dreidimensionalen Wellengleichung ist offenbar leicht in diese Form zu bringen,

$\partial _{t}^{2}u=f+c^{2}\cdot \Delta u$

$fu(x,0)=\varphi _{0}(x),\,\partial _{t}u(x,0)=\varphi _{1}(x)$

es liegt daher ein Cauchy-Problem der Ordnung 2 vor.

Sind alle Cauchy-Daten analytisch, so sichert der Satz von Cauchy-Kowalewskaja eindeutige Lösungen des Cauchy-Problems.

Bestimmung der Integrationskonstante

In der Schulmathematik wird die Bestimmung der Integrationskonstante eines unbestimmten Integrals für einen gegebenen Punkt als Anfangswertproblem bezeichnet.

Beispiel

Gesucht ist die Stammfunktion $F_{C}$ der gebrochenrationalen Funktion gegeben durch

$f(x)={\frac {1}{x^{2}-2x+1}}$ ,

die durch den Punkt $P(4|5)$ geht.

Zunächst faktorisieren wir den Nenner:

$f(x)={\frac {1}{(x-1)^{2}}}=(x-1)^{-2}$ .

Nun können wir substituieren:

$\int u^{-2}\cdot 1\,\mathrm {d} u=-u^{-1}=-{\frac {1}{x-1}}+C$ .

Als nächstes müssen wir die x-Koordinate des Punktes einsetzen und den Term mit dem y-Wert gleichsetzen

$-{\frac {1}{4-1}}+C=5\quad \Leftrightarrow \quad C={\frac {16}{3}}$ .

Die gesuchte Stammfunktion lautet demnach:

$F_{C}(x)=-{\frac {1}{x-1}}+{\frac {16}{3}}$ .

Abstraktes Cauchy-Problem

Seien ein Banachraum und $A\colon D(A)\subset X\rightarrow X$ ein linearer oder nichtlinearer Operator. Die Fragestellung, ob bei gegebenem T>0 , $u_{0}\in X$ und $f\colon (0,T)\rightarrow X$ eine differenzierbare Funktion $u\colon [0,T)\rightarrow X$ mit $u(t)\in D(A)$ für alle T>t>0 existiert, die das Anfangswertproblem

${\begin{matrix}u'(t)+A(u(t))&=&f(t),&\quad T>t>0\\u(0)&=&u_{0}&\end{matrix}}$

erfüllt, bezeichnet man als abstraktes Cauchy-Problem. Zu ihrer Lösbarkeit benötigt man die Theorie der stark stetigen Halbgruppen bzw. der analytischen Halbgruppen. Zu den verschiedenen Anfangsbedingungen und Operatoren gibt es verschiedene Arten des Lösungsbegriffes, im linearen distributionelle Lösungen, im nichtlinearen die integrale Lösung. Mit klassisch differenzierbaren, beziehungsweise fast überall differenzierbaren Lösungen beschäftigt sich die nachgelagerte Regularitätstheorie.

Literatur

Wolfgang Walter: Gewöhnliche Differentialgleichungen: Eine Einführung. 7. Auflage. Springer, 2000, ISBN 3-540-67642-2.
Martin Hermann: Numerik gewöhnlicher Differentialgleichungen. Band 1: Anfangswertprobleme und lineare Randwertprobleme. 2., überarbeitete und erweiterte Auflage. Walter de Gruyter Verlag, Berlin und Boston, 2017, ISBN 978-3-11-050036-3.

Basierend auf einem Artikel in:

Wikipedia.de