Lektion
Die Theorie — Gaußsches Eliminationsverfahren262 Wörter
Die Elimination überführt ein System in ein äquivalentes Dreieckssystem. Daraus kannst du die Lösung dann von unten nach oben ablesen. Ihre Bedeutung liegt nicht darin, dass sie funktioniert. Das tun andere Verfahren auch. Entscheidend ist, was sie kostet. Denn das bestimmt, welche Methode eine Maschine tatsächlich ausführt.
Was die einzelnen Symbole bedeuten
n- die Anzahl der Unbekannten und die einzige Größe, von der die Kosten abhängen.
pivot- der Koeffizient, durch den du eine Zeile teilst. Die vom Algorithmus festgelegte Reihenfolge macht das Verfahren entscheidungsfrei. Welche Zahl dort landet, entscheidet darüber, ob das Verfahren die Gleitkommaarithmetik übersteht.
- Setzt voraus
- Ein vollbesetztes System ohne ausnutzbare Struktur. Dünnbesetzte, Band- oder symmetrische Systeme sind günstiger. Die Wahl des passenden Spezialverfahrens macht den Großteil der numerischen linearen Algebra aus.
- Versagt, wenn
- Der Aufwand ist kubisch, und kubisch ist eine Wand. Die Elimination bei n Unbekannten erfordert etwa
n³/3Multiplikationen und ebenso viele Additionen. Zählst du bei einer echten Elimination direkt mit, kommst du bein = 100auf 333.300 Multiplikationen. Die Schätzung liefert 333.333. Dieser Exponent bedeutet eine feste Strafe für jedes Wachstum. Verdoppelst du die Unbekannten, vervielfacht sich die Arbeit um genau 8. Für 100 Unbekannte brauchst du 333.300 Multiplikationen. Bei 200 sind es schon 2.666.600. Eine Maschine, die tausend Unbekannte in einer Sekunde abarbeitet, braucht für zweitausend folglich acht Sekunden. Bei zwanzigtausend vergehen etwa zweieinviertel Stunden. Das ist die Gestalt dieser Wand. Deshalb lautet die spannende Frage in der numerischen linearen Algebra fast nie, wie du ein vollbesetztes System löst. Die Frage ist vielmehr, wie du vermeidest, überhaupt eines vor dir zu haben.
Aufgaben vollständig gelöst
-
Lösen von 2x + y − z = 8 und der Genauigkeitsverlust 7 Schritte
Löse 2x + y − z = 8, −3x − y + 2z = −11, −2x + y + 2z = −3 und berechne anschließend, was dich die Elimination an Genauigkeit gekostet hat.
-
Schreibe das System als erweiterte Koeffizientenmatrix auf. Die Buchstaben leisten hier keine echte Arbeit. Lass sie also weg und behalte nur die Spalten.
-
Nutze den größten Eintrag der ersten Spalte als Pivot. Das ist die −3 in der zweiten Zeile. Tausche diese Zeile nach oben. Da jeder Zeilentausch das Vorzeichen der Determinante umdreht, solltest du hier mitzählen.
-
Räume den Rest der ersten Spalte leer, indem du Vielfache der Pivotzeile abziehst.
-
Wiederhole das für die zweite Spalte. Lies dann die unterste Zeile ab: Sie enthält genau eine Unbekannte und ihren Wert.
-
Substituiere rückwärts nach oben. Die letzte Zeile liefert z, die mittlere liefert y, und aus der obersten erhältst du x.
-
Multipliziere die Pivots und wende das Vorzeichen für zwei Vertauschungen an, um die Determinante zu erhalten. Prüfe dein Ergebnis, indem du es in die ursprünglichen Gleichungen einsetzt und nicht in die reduzierten.
-
Das Residuum liegt bei etwa 10⁻¹⁶. Das ist ungleich null. Es entspricht exakt der Größenordnung jener Rundungsfehler, die ein Computer schlichtweg nicht vermeiden kann. Und genau das sollte man ehrlicherweise berichten. Eine Lösung, die von sich behauptet, absolut exakt zu sein, stammt entweder aus der Ganzzahlarithmetik oder schaut einfach nicht genau hin. Ändere nun bei ausgeschalteter Pivotisierung den führenden Koeffizienten auf 10⁻¹⁷. Das Residuum springt auf 1. Das ist dieselbe Warnung, nur brüllt sie dich diesmal an.
Antwort
x = 2, y = 3, z = −1, Determinante −1. Ein Residuum von etwa 10⁻¹⁶ zeigt, dass die Arithmetik genau so exakt arbeitet, wie es mit Fließkommazahlen eben möglich ist. Wenn du die Pivotisierung bei einem schlecht skalierten System abschaltest, ist genau diese Zahl das Signal dafür, dass deine Lösung falsch ist.
-
-
Zweimal eine Determinante von 0, mit zwei verschiedenen Bedeutungen 7 Schritte
Diese Voreinstellung hat keine Lösung und die nächste unendlich viele, und beide geben eine Determinante von 0 aus. Rechne heraus, was die beiden tatsächlich trennt, und entscheide, was ein Ändern der rechten Seite jemals bringen könnte.
-
Die Spaltenpivotisierung sucht in der ersten Spalte den größten Koeffizienten. Dort steht die 3 in der untersten Zeile, der erste Zug ist also ein Tausch.
-
Räum die Spalte leer: Zieh ⅔ der neuen ersten Zeile von der zweiten ab und ⅓ davon von der dritten.
-
In der zweiten Spalte steht 4/3 bereits über 2/3, es muss also nichts getauscht werden. Zieh die Hälfte von Zeile 2 von Zeile 3 ab.
-
Zeile 3 ist links nun leer und trägt rechts −½, was sich als 0 = −½ liest. Es wurden nur zwei Pivots gelegt, 3 und 4/3, wo drei Unbekannte drei brauchen, und ein fehlendes Pivot macht das Produkt null.
-
Der Rang ist das, was die Determinante verbirgt. Die Koeffizientenmatrix hat zwei unabhängige Zeilen. Die erweiterte Matrix hat drei, weil jenes übrig gebliebene −½ eine Zeile ist, die nichts wegkürzen kann. Rang 2 gegen Rang 3 ist genau das, was „keine Lösung“ bedeutet.
-
Die Voreinstellung Eine Gleichung dreimal rechnet dieselbe Arithmetik und sagt etwas anderes. Jede Zeile ist ein Vielfaches der ersten, also sind beide Ränge 1 und die Determinante wieder 0. Die Lösungsmenge hat 3 − 1 = 2 freie Richtungen: die ganze Ebene x + y + z = 3. Das Werkzeug sagt „unendlich viele“, ohne zu sagen, in wie vielen Richtungen.
-
Ändere eine Zahl und sieh, welche Größe sich bewegt. Setz eine 6 dorthin, wo rechts in Zeile 2 die 7 steht. Die Koeffizientenmatrix bleibt unberührt, die Determinante rührt sich also nicht, aber der Rang der erweiterten Matrix fällt auf 2, und aus keiner Lösung wird eine ganze Gerade von Lösungen.
Antwort
Die Determinante meldet nur, ob die Koeffizientenmatrix einen vollständigen Satz Pivots hat. Die rechte Seite kommt darin nie vor, und deshalb kann sie die beiden Voreinstellungen nicht trennen, während die beiden Ränge es können: keine Lösung, wenn sie sich widersprechen, und eine Lösungsmenge der Dimension 3 − Rang, wenn sie übereinstimmen. Ist die Determinante erst einmal 0, kauft dir keine rechte Seite mehr eine einzelne Antwort; du bekommst keine oder eine ganze Familie, und die Wahl zwischen beiden ist alles, worüber b entscheidet. Die Residuumskarte bleibt aus demselben Grund dunkel, aus dem es sie gibt: Ein Residuum braucht eine Lösung, die man wieder in die Gleichungen einsetzen kann.
-
Lernpfad
Nach x auflösen, vom Minuszeichen aufwärts
Quellen (1)
- Two centuries of the method being reinvented and misattributed, including how Gauss's name came to be on it: J. F. Grcar, "How ordinary elimination became Gaussian elimination." Historia Mathematica 38:2 (2011), 163–218.