Ftyjtk

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

y′(t)=f(t,y(t)){displaystyle y'(t)=f(t,y(t))}

und einer zusätzlichen Anfangsbedingung

y(t0)=y0{displaystyle y(t_{0})=y_{0}}

besteht, mit

• dem Anfangswert y0{displaystyle y_{0}} und


• einem Zeitpunkt t0∈R{displaystyle t_{0}in mathbb {R} }.

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

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

y(t)=y0+∫t0tf(s,y(s))ds{displaystyle y(t)=y_{0}+int _{t_{0}}^{t}f(s,y(s)),ds}

für alle t{displaystyle t} im Definitionsintervall gilt.^[1]

Anfangswertproblem k-ter Ordnung |

Gegeben seien k∈N{displaystyle kin mathbb {N} } und eine Funktion f:D→Rn{displaystyle fcolon Drightarrow mathbb {R} ^{n}}. Ihr Definitionsbereich D{displaystyle D} sei hierbei eine Teilmenge von I×Rn×k{displaystyle Itimes mathbb {R} ^{ntimes k}}, worin I⊂R{displaystyle Isubset mathbb {R} } ein Intervall bezeichnet, welches t0{displaystyle t_{0}} umfasst. Dann heißt

{y(k)=f(t,y(t),y′(t),…,y(k−1)(t))y(i)(t0)=yimiti=0,…,k−1{displaystyle {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-1end{cases}}}

ein Anfangswertproblem k{displaystyle k}-ter Ordnung. Jedes Anfangswertproblem k{displaystyle k}-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⋅sgn(y(t))⋅|y(t)| , y(0)=0 ,{displaystyle y'(t)=2cdot {rm {sgn}}(y(t))cdot {sqrt {|y(t)|}} , y(0)=0 ,}

welches zu

f(t,x):=2⋅sgn(x)⋅|x|{displaystyle f(t,x):=2cdot {rm {sgn}}(x)cdot {sqrt {|x|}}}

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

y(t)≡0{displaystyle y(t)equiv 0}

auch noch für jedes c≥0{displaystyle cgeq 0} die Lösungen

y(t)={0 ,falls t<c ,(t−c)2 ,falls t≥c ,{displaystyle y(t)=left{{begin{array}{ll}0 ,&{textrm {falls}} t<c ,\(t-c)^{2} ,&{textrm {falls}} tgeq c ,\end{array}}right.}

sowie

y(t)={0 ,falls t<c ,−(t−c)2 ,falls t≥c .{displaystyle y(t)=left{{begin{array}{ll}0 ,&{textrm {falls}} t<c ,\-(t-c)^{2} ,&{textrm {falls}} tgeq c .\end{array}}right.}

Damit Anfangswertprobleme eindeutige Lösungen besitzen, sind Zusatzeigenschaften (an f{displaystyle f}) 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 n{displaystyle n} Veränderliche x1,…,xn{displaystyle x_{1},dotsc ,x_{n}}, so erhält man partielle Differentialgleichungen. Im Folgenden stehe α∈N0n{displaystyle alpha in mathbb {N} _{0}^{n}} für einen Multiindex der Länge n{displaystyle n}. Beachte, dass es genau (n+k−1k){displaystyle {tbinom {n+k-1}{k}}} Multiindizes mit |α|:=α1+⋯+αn≤k{displaystyle |alpha |:=alpha _{1}+dotsb +alpha _{n}leq k} gibt. Es sei weiter eine Funktion F{displaystyle F} in n+(n+k−1k){displaystyle n+{tbinom {n+k-1}{k}}} Variablen gegeben. Beim allgemeinen Cauchy-Problem sucht man nach Funktionen u{displaystyle u}, die von n{displaystyle n} Variablen x1,…,xn{displaystyle x_{1},dotsc ,x_{n}} abhängen und die Gleichung

(1) F(x,(∂αu(x))|α|≤k))=0{displaystyle F(x,(partial ^{alpha }u(x))_{|alpha |leq k)}),=,0}

erfüllen. Beachte, dass die Stelligkeit von F{displaystyle F} gerade so gewählt wurde, dass man x=(x1,…,xn){displaystyle x=(x_{1},dotsc ,x_{n})} und alle partiellen Ableitungen ∂αu(x){displaystyle 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 S{displaystyle S} eine Hyperfläche der Klasse C^k mit Normalenfeld ν{displaystyle nu }. Mit ∂νj{displaystyle partial _{nu }^{j}} seien die Normalenableitungen bezeichnet. Sind dann φ0,…,φk−1{displaystyle varphi _{0},dotsc ,varphi _{k-1}} vorgegebene auf S{displaystyle S} definierte Funktionen, so fordert man beim allgemeinen Cauchy-Problem, dass die Funktionen u{displaystyle u} zusätzlich die Bedingungen

(2) u=φ0,∂νu=φ1,…,∂νk−1u=φk−1{displaystyle u=varphi _{0},,partial _{nu }u=varphi _{1},,dotsc ,,partial _{nu }^{k-1}u=varphi _{k-1}} auf S{displaystyle S}

erfüllen. Die Funktionen φj{displaystyle varphi _{j}} heißen die Cauchy-Daten des Problems, jede Funktion u{displaystyle u}, 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=(x1,…,xn);xn=0}{displaystyle 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 t{displaystyle t} (lateinisch tempus = Zeit) um, die Anfangsbedingungen beschreiben dann die Verhältnisse zum Zeitpunkt t=0{displaystyle t=0}. Die Variablen sind also x1,…,xn−1,t{displaystyle x_{1},dotsc ,x_{n-1},t}. Da die betrachtete Hyperebene durch die Bedingung t=0{displaystyle t=0} gegeben ist, wird die Normalenableitung einfach zur Ableitung nach t{displaystyle t}. Schreibt man abkürzend x=(x1,…,xn−1){displaystyle x=(x_{1},dotsc ,x_{n-1})} und α=(α1,…,αn−1){displaystyle alpha =(alpha _{1},dotsc ,alpha _{n-1})}, so lautet das Cauchy-Problem nun

(1') F(x,t,(∂xα∂tju(x,t))|α|+j≤k))=0{displaystyle F(x,t,(partial _{x}^{alpha }partial _{t}^{j}u(x,t))_{|alpha |+jleq k)}),=,0}

(2') u(x,0)=φ0(x),∂tu(x,0)=φ1(x),…,∂tk−1u(x,0)=φk−1(x){displaystyle 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

∂t2u−c2⋅Δu=f{displaystyle partial _{t}^{2}u-c^{2}cdot Delta u=f}

u(x,0)=φ0(x),∂tu(x,0)=φ1(x){displaystyle u(x,0)=varphi _{0}(x),,partial _{t}u(x,0)=varphi _{1}(x)},

wobei c{displaystyle c} eine Konstante, f{displaystyle f} eine vorgegebene Funktion und Δ=∂x12+∂x22+∂x32{displaystyle Delta =partial _{x_{1}}^{2}+partial _{x_{2}}^{2}+partial _{x_{3}}^{2}} der Laplace-Operator seien.

Ist u{displaystyle u} eine Lösung, was gleichzeitig ausreichende Differenzierbarkeit implizieren soll, so sind alle Ableitungen ∂xα∂tju(x,0){displaystyle partial _{x}^{alpha }partial _{t}^{j}u(x,0)} mit |α|+j≤k,j<k{displaystyle |alpha |+jleq k,j<k} bereits durch die Cauchy-Daten vorgegeben, denn es ist ∂xα∂tju(x,0)=∂xαφj{displaystyle partial _{x}^{alpha }partial _{t}^{j}u(x,0)=partial _{x}^{alpha }varphi _{j}}. Lediglich die Ableitung ∂tku{displaystyle 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 ∂tku{displaystyle partial _{t}^{k}u} auflösen kann. Das Cauchy-Problem hat dann die Form

(1") ∂tku(x,t)=G(x,t,(∂xα∂tju(x,t))|α|+j≤k,j<k){displaystyle partial _{t}^{k}u(x,t)=G(x,t,(partial _{x}^{alpha }partial _{t}^{j}u(x,t))_{|alpha |+jleq k,j<k})}

(2") u(x,0)=φ0(x),∂tu(x,0)=φ1(x),…,∂tk−1u(x,0)=φk−1(x){displaystyle 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 G{displaystyle G} eine geeignete Funktion der Stelligkeit n−1+(n+k−1k){displaystyle n-1+{tbinom {n+k-1}{k}}} sei.
In der zuletzt gegebenen Formulierung haben alle auftretenden Ableitungen eine Ordnung ≤k{displaystyle leq k}, und die k{displaystyle k}-te Ableitung nach t{displaystyle t} 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 k{displaystyle k} daher auch die Ordnung des Cauchy-Problems.
Das obige Beispiel der dreidimensionalen Wellengleichung ist offenbar leicht in diese Form zu bringen,

∂t2u=f+c2⋅Δu{displaystyle partial _{t}^{2}u=f+c^{2}cdot Delta u}

u(x,0)=φ0(x),∂tu(x,0)=φ1(x){displaystyle u(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.

Abstraktes Cauchy-Problem |

Seien X{displaystyle X} ein Banachraum und A:D(A)⊂X→X{displaystyle Acolon D(A)subset Xrightarrow X} ein linearer oder nichtlinearer Operator. Die Fragestellung, ob bei gegebenem T>0{displaystyle T>0}, u0∈X{displaystyle u_{0}in X} und f:(0,T)→X{displaystyle fcolon (0,T)rightarrow X} eine differenzierbare Funktion u:[0,T)→X{displaystyle ucolon [0,T)rightarrow X} mit u(t)∈D(A){displaystyle u(t)in D(A)} für alle T>t>0{displaystyle T>t>0} existiert, die das Anfangswertproblem

u′(t)+A(u(t))=f(t),T>t>0u(0)=u0{displaystyle {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.


• Isao Miyadera, Choong Yun Cho: Nonlinear Semigroups. American Mathemat. Soc., Providence, RI 1992, ISBN 0-8218-4565-9.


• Amnon Pazy: Semigroups of Linear Operators and Applications to Partial Differential Equations. Springer-Verlag, New York 1983, ISBN 0-387-90845-5.


• Gerald B. Folland: Introduction to Partial Differential Equations. Princeton University Press, 1976, ISBN 0-691-08177-8. (insbesondere Kapitel 1.C. für das allgemeine Cauchy-Problem)

Einzelnachweise |

1. 
   ↑ Rannacher, Rolf: Numerik 1. Numerik gewöhnlicher Differentialgleichungen. Heidelberg 2017, S. 13. 
   

Weblinks |

• Gert Lube: Anfangswertaufgaben. (Skript, Universität Göttingen)


• Clemens Brand: Illustrationen zu einem einfachen Anfangswertproblem. (Skript, Uni Leoben)


• 
  taramath Online-Tool zur Lösung von Anfangswertproblemen.