Splitting-Verfahren

In der numerischen Mathematik sind Splitting-Verfahren iterative Verfahren zum Lösen linearer Gleichungssysteme Ax=bAx=b mit einer Matrix A∈Cn×nA\in\mathbb{C}^{n\times n} und rechter Seite b∈Cn b\in\mathbb{C}^n\, Im Unterschied zu direkten Verfahren nähert man sich dabei ausgehend von einer Startnäherung schrittweise der gesuchten Lösung an und bricht ab, falls die Genauigkeit hoch genug ist.

Beschreibung

Das Verfahren ergibt sich über ein Splitting der Systemmatrix A=B+(A−B)A = B + (A-B) mit einer invertierbaren Matrix B∈Cn×nB\in\mathbb{C}^{n\times n}.
Ax=b⇔B−1(B+(A−B))x=B−1b Ax=b \Leftrightarrow B^{-1}(B + (A-B))x=B^{-1}b\,
Daraus erhält man die Fixpunktgleichung
x=B−1(B−A)x+B−1b x = B^{-1}(B-A)x + B^{-1}b.
Mit M=B−1(B−A)M = B^{-1}(B-A) ergibt sich die Fixpunktiteration
  1. Wähle einen Startvektor x0∈Cnx_0\in\mathbb{C}^n
  2. Setze xk=Mxk−1+B−1b=(I−B−1A)xk−1+B−1bx_k=Mx_{k-1}+B^{-1}b = (I - B^{-1}A)x_{k-1} + B^{-1}b
Man kann die Iteration abbrechen, falls die Norm des Residuums rk=b−Axkr_k=b-Ax_k eine vorgegebene Fehlerschranke unterschreitet. Das Verfahren konvergiert genau dann, wenn der Spektralradius der Matrix MM kleiner 11 ist. Mit Hilfe des Banachschen Fixpunktsatzes folgt ferner die lineare Konvergenzgeschwindigkeit der gesamten Verfahrensklasse. Je kleiner der Spektralradius ist, umso schneller konvergiert das Verfahren. Falls sich BB und AA nur wenig unterscheiden kann man mit dem Störungslemma zeigen, dass auch der Spektralradius von MM klein ist. Damit ergibt sich ein Gegensatz von schneller Konvergenz (BB approximiert AA sehr gut) zu geringen Kosten pro Iteration (BB ist einfach invertierbar). Insgesamt sind diese Verfahren für viele praktische Probleme zu langsam. Allerdings stellen sie aufgrund ihrer einfachen Anwendbarkeit gute Vorkonditionierer dar. Darüberhinaus sind viele Splitting-Verfahren als Glätter in einem Mehrgitterverfahren geeignet.

Beispiele

Modifikationen

Man unterscheidet stationären Verfahren - in dem Fall ist die Iterationsmatrix konstant und instationäre Verfahren - hier dürfen die Matrizen MM vom Index ii abhängen. τ∈R\tau\in\mathbb{R} ist ein Regularisierungsparameter der vom konkreten Problem und vom Verfahren abhängt.

Literatur

 
 

An Archimedes wird man sich erinnern, wenn Aischylos vergessen ist - weil zwar die Sprachen sterben, nicht aber die mathematischen Ideen.

Godfrey Harold Hardy

Copyright- und Lizenzinformationen: Diese Seite basiert dem Artikel Splitting-Verfahren aus der frеiеn Enzyklοpädιe Wιkιpеdιa und stеht unter der Dοppellizеnz GNU-Lιzenz für freie Dokumentation und Crеative Commons CC-BY-SA 3.0 Unportеd (Kurzfassung). In der Wιkιpеdιa ist eine Listе dеr Autorеn des Originalartikels verfügbar. Da der Artikel geändert wurde, reicht die Angabe dieser Liste für eine lizenzkonforme Weiternutzung nicht aus!
Anbieterkеnnzeichnung: Mathеpеdιa von Тhοmas Stеιnfеld  • Dοrfplatz 25  •  17237 Blankеnsее  • Tel.: 01734332309 (Vodafone/D2)  •  Email: cο@maτhepedιa.dе
Datenschutzerklärung

 

Numerische Verfahren