Lokaler Diskretisierungsfehler

Ein lokaler Diskretisierungsfehler ( ${\displaystyle \delta }$) oder Abschneidefehler ist in der Numerik gewöhnlicher Differentialgleichungen ein bestimmter Fehler, der innerhalb eines Verfahrensschritts auftritt.

Dieser Artikel wurde auf der Qualitätssicherungsseite des Portals Mathematik eingetragen. Dies geschieht, um die Qualität der Artikel aus dem Themengebiet Mathematik auf ein akzeptables Niveau zu bringen. Bitte hilf mit, die Mängel dieses Artikels zu beseitigen, und beteilige dich bitte an der Diskussion! (Artikel eintragen)

Dieser Beitrag beschränkt sich auf Anfangswertprobleme gewöhnlicher Differentialgleichungen (AWP) sowie auf numerische Einschrittverfahren (ESV), die von einer Diskretisierungsschrittweite (${\displaystyle h}$) abhängen. Damit das ESV einen Bezug zum gegebenen AWP besitzt, muss gezeigt werden, dass dieses numerische Ersatzproblem für ${\displaystyle h\rightarrow 0}$ in das gegebene AWP übergeht. Hierbei spielt der lokale Diskretisierungsfehler ${\displaystyle \delta (\cdot )}$ eine wichtige Rolle. Allgemein gesprochen ist der lokale Diskretisierungsfehler das durch die Schrittweite ${\displaystyle h}$ dividierte Residuum, welches sich ergibt, wenn man die exakte Lösung des AWPs in die Formel des ESVs einsetzt. Offensichtlich geht es hier nicht darum, dass die aus dem ESV berechnete Lösung gegen die exakte Lösung des AWPs konvergiert, sondern es wird erst einmal gefordert, dass das Ersatzproblem in einer gewissen Beziehung zum gegebenen AWP steht. Des Weiteren lassen sich auf der Basis des lokalen Diskretisierungsfehlers Algorithmen zur automatischen Steuerung der Schrittweite im jeweiligen ESV konstruieren.

Der oben genannte Sachverhalt soll nun etwas genauer dargestellt werden. Zur Lösung des folgenden Anfangswertproblems (AWPs) einer gewöhnlichen Differentialgleichung

${\displaystyle {\dot {x}}(t):={\frac {d}{dt}}x(t)=f(t,x(t)),x(t_{0})=x_{0}}$,

mit ${\displaystyle t\in J\subset \mathbb {R} }$ und ${\displaystyle x(t)\in \Omega \subset \mathbb {R} ^{n}}$,

sei auf dem Gitter ${\displaystyle J_{h}:=\{t_{0},t_{1},...,t_{N}\}}$ das implizite Einschrittverfahren (ESV)

${\displaystyle x_{i+1}^{h}=x_{i}^{h}+h\Phi (t_{i},x_{i}^{h},x_{i+1}^{h};h)}$

gegeben. Es handelt sich dabei um eine Differenzenapproximation des AWPs.

Der Differenzenoperator

Fehler beim Parsen (SVG (MathML kann über ein Browser-Plugin aktiviert werden): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „http://localhost:6011/de.wikipedia.org/v1/“:): {\displaystyle \triangle_h \, u(t_{i+1}) := \frac{1}{h}\,[u(t_{i+1}) - u(t_i) - h\, \Phi(t_i,u(t_i),u(t_{i+1});h)]}

werde für ${\displaystyle i=0,1,\ldots ,N}$ auf eine Gitterfunktion ${\displaystyle \{u(t_{i})\}_{i=0}^{N},t_{i}\in J_{h},}$ angewendet, für die ${\displaystyle u(t_{0})}$ fixiert ist. Des Weiteren sei ${\displaystyle x^{h}}$ eine approximierende Gitterfunktion, die den Wert ${\displaystyle x_{i}^{h}}$ an der Stelle ${\displaystyle t_{i},i=0,1,\ldots ,N}$, annimmt. Offensichtlich erfüllt das obige ESV die Beziehung

${\displaystyle \triangle _{h}\,x^{h}(t_{i+1})=0}$.

Die numerische Analyse gewöhnlicher Differentialgleichungen beschäftigt sich schwerpunktmäßig mit den Fehlern, die bei jedem Integrationsschritt durch die verwendete Differenzenapproximation hervorgerufen werden und wie diese sich während der Rechnung akkumulieren. Ein Maß für den Fehler, der in einem Schritt auf dem Gitter zu verzeichnen ist, stellt der lokale Diskretisierungsfehler ${\displaystyle \delta (\cdot )}$ dar. Er ist das Residuum des Differenzenoperators ${\displaystyle \triangle _{h}}$, wenn dieser auf die exakte Lösung des AWPs angewendet wird, d. h.

${\displaystyle \delta (t_{i+1},x(t_{i+1});h):=\triangle _{h}\,x(t_{i+1})={\frac {1}{h}}\,[x(t_{i+1})-x(t_{i})-h\,\Phi (t_{i},x(t_{i}),x(t_{i+1});h)].}$

Somit gibt der lokale Diskretisierungsfehler an, wie gut der Differenzenoperator den Differentialoperator approximiert. Es erweist sich, dass eine grundlegende Voraussetzung, die ein ESV erfüllen muss, dessen Konsistenz ist. Diese wichtige Eigenschaft der ESVn wird auf der Basis des lokalen Diskretisierungsfehlers definiert. So ist das ESV konsistent (bezüglich des gegebenen AWPs) mit der Konsistenzordnung ${\displaystyle p}$, wenn ${\displaystyle p}$ die größte ganze Zahl bezeichnet, für die der lokale Diskretisierungsfehler die Beziehung

${\displaystyle \delta (t_{i+1},x(t_{i+1});h)=O(h^{p}),\quad \;h\rightarrow 0,}$

erfüllt. Dabei impliziert ${\displaystyle O(h^{p})}$ die Existenz endlicher Konstanten ${\displaystyle C}$ und ${\displaystyle h_{0}>0}$, so dass

${\displaystyle \delta (t_{i+1},x(t_{i+1});h)\leq C\,h^{p}}$

ist, für ${\displaystyle h\leq h_{0}.}$ Konsistenz bedeutet im Allgemeinen, dass für die Ordnung des Verfahrens ${\displaystyle p\geq 1}$ gilt.

Das folgende Beispiel soll den dargestellten Sachverhalt illustrieren. Betrachtet werde das Euler(vorwärts)-Verfahren

${\displaystyle x_{i+1}^{h}=x_{i}^{h}+hf(x_{i}^{h})}$

zur Lösung des AWPs einer skalaren und autonomen Differentialgleichung erster Ordnung

${\displaystyle {\dot {x}}(t)=f(x(t)),\quad x(t_{0})=x_{0}.}$

Unter hinreichenden Glattheitsvoraussetzungen an die Funktion ${\displaystyle f}$ lässt sich die zugehörige Lösung unter Verwendung einer Taylor-Entwicklung in der Form

${\displaystyle x(t_{i+1})=x(t_{i}+h)=x(t_{i})+h{\dot {x}}(t_{i})+O(h^{2})}$

darstellen. Verwendet man den über die Differentialgleichung bestehenden Zusammenhang zwischen der ersten Ableitung von ${\displaystyle x}$ und ${\displaystyle f(x)}$, so resultiert

${\displaystyle x(t_{i+1})=x(t_{i})+hf(x(t_{i}))+O(h^{2}).}$

Damit ergibt sich für den lokalen Diskretisierungsfehler des Euler(vorwärts)-Verfahrens

${\displaystyle {\begin{aligned}\delta (t_{i+1},x(t_{i+1});h)&={\frac {1}{h}}[x(t_{i+1)}-x(t_{i})-hf(x(t_{i}))]\\&={\frac {1}{h}}[x(t_{i})+hf(x(t_{i}))+O(h^{2})-x(t_{i})-hf(x(t_{i}))]\\&=O(h),\end{aligned}}}$(ESV) im Grenzwert in das AWP übergeht.

d. h. das Euler(vorwärts)-Verfahren besitzt die Konsistenzordnung ${\displaystyle p=1}$. Insbesondere geht der lokale Diskretisierungsfehler linear in ${\displaystyle h}$ für ${\displaystyle h\rightarrow 0}$ gegen Null, was zumindest garantiert, dass das numerische Ersatzproblem (d. h. das ESV) im Grenzwert in das AWP übergeht.