RREF

Als RREF wird die reduzierte Zeilenstufenform bezeichnet. Die Abkürzung RREF kommt vom Englischen “Reduced Row Echelon Form”.

Eine RREF kann die folgende Form haben:

Alle Kästchen, die in in weiss sind müssen zwingend mit einer gefüllt sein. Kästchen, welche grau hinterlegt sind können jeden beliebigen Wert annehmen, sofern kein Wert eingetragen ist.

In jeder Spalte steht ein Standardeinheitsvektor. Die Spalten sind dabei genau die Spalten, bei welchen der graue Bereich eine “Stufe” nach unten macht.

Beispiele

  • Der RREF von ist: RREF()
  • Der RREF von ist RREF()
Link to original

Gauss-Jordan-Elimination

Bei der Gauss-Jordan-Elimination wird wie schon bei der Gauss-Elimination das Gleichungssystem in ein anderes schöneres System umgewandelt, welches dann einfach zu lösen ist. Die folgende Umwandlung wird bei der Gauss-Jordan-Elimination durchgeführt:
wobei in Reduzierter Zeilenstufenform (RREF) geschrieben wird.

Direkte Lösung

Die direkte Lösung von , in RREF()
Die Resultate können direkt aus der Matrix In Zeilenstufenform abgelesen werden. So ist die Lösung einfach das Skalar, mit welchem der Standardeinheitsvektor an der benötigten Spalte multipliziert werden muss, damit man auf den Wert von gelangt. Alle anderen Einträge von sind .

Bsp.
Bei muss der Standardeinheitsvektor mit multipliziert werden, damit man auf den gewünschten Wert in kommt.

Elimination

Wenn noch nicht in RREF ist, muss die folgende Transformation stattfinden:
Die folgende Illustration zeigt, wie das Durchführen der Gauss-Jordan-Elimination und die Transformation der Matrix zur Matrix funktioniert. Alle Schritte müssen auch auf der rechten Seite durchgeführt werden.

Das Vorgehen ist nun gleich wie bei der Gauss-Elimination (unter dem Element der Diagonalen dürfen in dieser Spalte nur die stehen). Im Gegensatz zur Gauss-Elimination wird hier aber zuerst das Pivot-Element auf gebracht.

Bei der Gauss-Elimination würde dieser Fall nun zum Scheitern und Aufgeben führen, da es keine Zeilen mehr gibt, welche für Spalte sind. Bei der Gauss-Jordan-Elimination setzten wird dies aber nicht voraus und akzeptieren auch, dass es eine Spalte gibt, in welcher dort eine steht.

In Spalte 3 gibt es jetzt einen Eintrag , er steht aber noch in Zeile 3 und nicht 2 wo er sein müsste. Das ist aber kein Problem, die Zeilen dürfen vertauscht werden. Wie oben gezeigt muss jetzt das Pivot Element auf gebracht werden, was hier durch Division mit erzielt wird.

Bei der Gauss-Elimination musste man sich nur darum kümmern, dass die Einträge unterhalb des Pivot Elements sind. Hier ist es so, dass alle anderen Elemente in dieser Spalte sein müssen. Die Zeile muss also einmal von Zeile abgezogen werden.

Das wird nun auch noch für Spalte ausgeführt und es entsteht die folgende Matrix in RREF:

Dies ist nun die aus gewonnene Matrix , aus welcher wie in Direkte Lösung gezeigt die Lösung der Gleichung abgelesen werden kann.

Link to original

Der folgende Algorithmus zeigt auf, wie das ganze programmatisch funktioniert und ausserdem, wie mit der rechten Seite vorgegangen wird:

Gaus-Jordan-Elimination-Algorithmus

Link to original

Die Laufzeit des obigen Algorithmus ():

  • Es müssen Zeilenoperationen müssen ausgeführt werden
  • Pro Zeilenoperation gibt es Einzelschritte

Die Gesamtlaufzeit beträgt daher:

Im Falle läuft das auf hinaus. Somit hat das Gauss-Jordan Verfahren für den Fall, in welchem das Gauss Verfahren angewendet werden kann die gleiche Laufzeit.

ist die RREF-Standardform von und gibt uns auch die CR-Zerlegung nach

Theorem 3.18

Sei eine Matrix. Es gibt genau eine Matrix (die aus Theorem 3.17) mit den folgenden Eigenschaften:

  • (i) mit einer invertierbaren Matrix
  • (ii) ist in RREF
Link to original
.
Genauer: ist in RREF(), wobei die Indizes er unabhängigen Spalten von sind, und

mit der eindeutigen Matrix, so dass (Submatrix der unabhängigen Spalten) in Theorem 2.44.

Berechnung der Inversen von mit Gauss-Jordan() nach Theorem 3.19