Trunkierungsfehler (numerische Integration) – Wikipedia
Fehler bei der numerischen Integration
Kürzungsfehler bei der numerischen Integration gibt es zwei Arten:
- lokale Kürzungsfehler – der durch eine Iteration verursachte Fehler und
- globale Kürzungsfehler – der kumulative Fehler, der durch viele Iterationen verursacht wird.
Definitionen[edit]
Angenommen, wir haben eine stetige Differentialgleichung
und wir wollen eine Näherung berechnen
der wahren Lösung
in diskreten Zeitschritten
. Nehmen Sie der Einfachheit halber an, dass die Zeitschritte gleichmäßig verteilt sind:
Angenommen, wir berechnen die Folge
mit einer einstufigen Methode der Form
Die Funktion
heißt der Inkrementfunktion, und kann als Schätzung der Steigung interpretiert werden
.
Lokaler Kürzungsfehler[edit]
Die lokaler Kürzungsfehler
ist der Fehler, den unsere Inkrement-Funktion,
, Ursachen während einer einzelnen Iteration, vorausgesetzt, dass die wahre Lösung bei der vorherigen Iteration vollkommen bekannt ist.
Formaler ist der lokale Trunkierungsfehler,
, bei Schritt
berechnet sich aus der Differenz zwischen linker und rechter Seite der Gleichung für das Inkrement
:
- [1][2]
Die numerische Methode ist konsistent wenn der lokale Kürzungsfehler ist
(das bedeutet, dass für jeden
so dass
für alle
; siehe Little-O-Notation). Wenn die Inkrement-Funktion
stetig ist, dann ist die Methode genau dann konsistent, wenn
.[3]
Weiterhin sagen wir, dass die numerische Methode Auftrag
falls für eine hinreichend glatte Lösung des Anfangswertproblems der lokale Abschneidefehler
(was bedeutet, dass es Konstanten gibt
und
so dass
für alle
).[4]
Globaler Kürzungsfehler[edit]
Die globaler Kürzungsfehler ist die Ansammlung von lokaler Kürzungsfehler über alle Iterationen unter der Annahme einer perfekten Kenntnis der wahren Lösung im anfänglichen Zeitschritt.[citation needed]
Formaler ist der globale Trunkierungsfehler,
, zum Zeitpunkt
ist definiert durch:
- [5]
Die numerische Methode ist konvergent wenn der globale Kürzungsfehler auf null geht, wenn die Schrittgröße auf null geht; mit anderen Worten, die numerische Lösung konvergiert gegen die exakte Lösung:
.[6]
Zusammenhang zwischen lokalen und globalen Kürzungsfehlern[edit]
Manchmal ist es möglich, eine Obergrenze für den globalen Trunkierungsfehler zu berechnen, wenn wir den lokalen Trunkierungsfehler bereits kennen. Dies erfordert, dass sich unsere Inkrementfunktion ausreichend gut verhält.
Der globale Trunkierungsfehler erfüllt die Wiederholungsbeziehung:
Dies folgt unmittelbar aus den Definitionen. Nehmen wir nun an, dass die Inkrementfunktion im zweiten Argument Lipschitz-stetig ist, d. h. es existiert eine Konstante
so dass für alle
und
und
, wir haben:
Dann erfüllt der globale Fehler die Schranke
- [7]
Aus der obigen Schranke für den globalen Fehler folgt, dass wenn die Funktion
in der Differentialgleichung ist im ersten Argument stetig und im zweiten Argument Lipschitz-stetig (die Bedingung aus dem Satz von Picard-Lindelöf), und die Inkrementfunktion
in allen Argumenten stetig und im zweiten Argument Lipschitz stetig ist, dann geht der globale Fehler als Schrittweite gegen Null
gegen Null geht (mit anderen Worten, die numerische Methode konvergiert gegen die exakte Lösung).[8]
Erweiterung auf lineare Mehrschrittverfahren[edit]
Betrachten Sie nun eine lineare Mehrschrittmethode, gegeben durch die Formel
Somit wird der nächste Wert für die numerische Lösung berechnet nach
Die nächste Iteration einer linearen Mehrschrittmethode hängt von der vorherigen ab S iteriert. Somit wird bei der Definition für den lokalen Trunkierungsfehler nun davon ausgegangen, dass der vorherige S iteriert alle entsprechen der exakten Lösung:
- [9]
Auch hier ist die Methode konsistent, wenn
und es hat ordnung P wenn
. Auch die Definition des globalen Trunkierungsfehlers bleibt unverändert.
Die Beziehung zwischen lokalen und globalen Trunkierungsfehlern unterscheidet sich geringfügig von der einfacheren Einstellung von Einschrittverfahren. Für lineare Mehrschrittverfahren wird ein zusätzliches Konzept namens Nullstabilität benötigt, um die Beziehung zwischen lokalen und globalen Trunkierungsfehlern zu erklären. Lineare Mehrschrittverfahren, die die Bedingung der Nullstabilität erfüllen, haben dieselbe Beziehung zwischen lokalen und globalen Fehlern wie Einschrittverfahren. Mit anderen Worten, wenn ein lineares Mehrschrittverfahren nullstabil und konsistent ist, dann konvergiert es. Und wenn eine lineare Mehrschrittmethode nullstabil ist und einen lokalen Fehler hat
, dann ist sein globaler Fehler erfüllt
.[10]
Siehe auch[edit]
- ^ Gupta, GK; Sacks-Davis, R.; Tischer, PE (März 1985). “Ein Überblick über die jüngsten Entwicklungen bei der Lösung von ODEs”. Computerumfragen. 17 (1): 5–47. CiteSeerX 10.1.1.85.783. mach:10.1145/4078.4079.
- ^ Süli & Mayers 2003, S. 317, Anrufe der Kürzungsfehler.
- ^ Süli & Mayers 2003, S. 321 & 322
- ^ Iserles 1996, p. 8; Süli & Mayers 2003, S. 323
- ^ Süli & Mayers 2003, S. 317
- ^ Iserles 1996, p. 5
- ^ Süli & Mayers 2003, p. 318
- ^ Süli & Mayers 2003, p. 322
- ^ Süli & Mayers 2003, S. 337, verwendet eine andere Definition und dividiert diese im Wesentlichen durch h
- ^ Süli & Mayers 2003, S. 340
Recent Comments