Um ein System von linearen Gleichungen zu lösen, ist die Verwendung der reduzierten Zeilen-Echelon-Form (RREF) ein bewährtes Verfahren. Dies wird häufig als Gauss-Jordan-Elimination bezeichnet. Dabei wird ein lineares Gleichungssystem durch eine Reihe elementarer Zeilenoperationen in eine vereinfachte Matrixform überführt, die die Lösung des Systems direkt ablesen lässt.

Betrachten wir zunächst ein allgemeines lineares Gleichungssystem:

a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2am1x1+am2x2++amnxn=bm\begin{aligned}
a_{11}x_1 + a_{12}x_2 + \dots + a_{1n}x_n &= b_1 \\ a_{21}x_1 + a_{22}x_2 + \dots + a_{2n}x_n &= b_2 \\ \vdots \\ a_{m1}x_1 + a_{m2}x_2 + \dots + a_{mn}x_n &= b_m \end{aligned}

Dieses System lässt sich in Matrixform als Ax=bA \cdot x = b darstellen, wobei AA die Koeffizientenmatrix des Systems ist, xx der Vektor der Unbekannten und bb der Vektor der konstanten Terme.

Erstellen der erweiterten Matrix

Um das System zu lösen, erweitern wir die Matrix AA mit dem Vektor bb zu einer sogenannten erweiterten Matrix:

[Ab][A | b]

Dabei wird der Vektor bb an die rechte Seite der Matrix angehängt. Das Ziel ist es nun, diese erweiterte Matrix durch elementare Zeilenoperationen in die reduzierte Zeilen-Echelon-Form zu überführen. Diese Form ermöglicht es uns, die Lösungen des Systems direkt abzulesen.

Elementare Zeilenoperationen

Es gibt drei grundlegende Arten von elementaren Zeilenoperationen:

  1. Vertauschen von Zeilen: Diese Operation tauscht zwei Zeilen der Matrix aus.

  2. Multiplizieren einer Zeile mit einer Konstante: Eine Zeile wird mit einer konstanten Zahl multipliziert, die nicht null ist.

  3. Addieren eines Vielfachen einer Zeile zu einer anderen Zeile: Diese Operation ermöglicht es, eine Zeile durch Addition eines Vielfachen einer anderen Zeile zu verändern.

Durch Anwendung dieser Operationen können wir die Matrix so umformen, dass sie die reduzierte Zeilen-Echelon-Form erreicht.

Reduzierte Zeilen-Echelon-Form

Eine Matrix befindet sich in reduzierter Zeilen-Echelon-Form, wenn sie folgende Bedingungen erfüllt:

  • Jede führende Eins in einer Zeile steht in einer neuen Spalte, die von allen vorherigen Zeilen führende Einsen enthält.

  • Alle Zeilen, die nur Nullen enthalten, stehen am Ende der Matrix.

  • Alle Einträge unter und über den führenden Einsen sind null.

Durch die Anwendung der Gauss-Jordan-Elimination können wir die Matrix in diese Form bringen. Sobald dies erreicht ist, können wir die Lösung des Gleichungssystems direkt ablesen. Wenn eine Zeile der erweiterten Matrix eine Gleichung wie 0=c0 = c (mit c0c \neq 0) ergibt, ist das System inkonsistent, und es existiert keine Lösung. Falls keine solche Gleichung auftritt, ist das System konsistent, und wir können die Werte für x1,x2,,xnx_1, x_2, \dots, x_n berechnen.

Konsistenz und Lösungen des Systems

Ein System von linearen Gleichungen kann zwei mögliche Ausgänge haben:

  1. Konsistentes System: Das System hat mindestens eine Lösung. In diesem Fall können wir die Lösung(en) in der reduzierten Zeilen-Echelon-Form direkt ablesen.

  2. Inkonsistentes System: Das System hat keine Lösung. Dies tritt auf, wenn während des Umformungsprozesses eine Zeile der Form 0=c0 = c entsteht, wobei c0c \neq 0 ist. In diesem Fall ist das Gleichungssystem widersprüchlich und daher unlösbar.

Beispiel: Lösen eines Systems von linearen Gleichungen

Nehmen wir an, wir haben das folgende lineare Gleichungssystem:

2y+4z=12x+4y+2z=13x+3y+z=1\begin{aligned}
2y + 4z &= 1 \\ 2x + 4y + 2z &= 1 \\ 3x + 3y + z &= 1 \end{aligned}

Dieses System kann als Matrix geschrieben werden:

A=(024121331)undb=(111)A = \begin{pmatrix}
0 & 2 & 4 \\ 1 & 2 & 1 \\ 3 & 3 & 1 \end{pmatrix} \quad \text{und} \quad b = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}

Die erweiterte Matrix ist daher:

(024112113311)\left( \begin{array}{ccc|c}
0 & 2 & 4 & 1 \\ 1 & 2 & 1 & 1 \\ 3 & 3 & 1 & 1 \end{array} \right)

Durch Anwenden der Gauss-Jordan-Elimination erhalten wir die reduzierte Zeilen-Echelon-Form:

(10014010000112)\left( \begin{array}{ccc|c}
1 & 0 & 0 & \frac{1}{4} \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & \frac{1}{2} \end{array} \right)

Die Lösung des Systems lautet also:

x=14,y=0,z=12x = \frac{1}{4}, \quad y = 0, \quad z = \frac{1}{2}

Inverse einer Matrix und Anwendungen

Neben der Lösung von Gleichungssystemen ist die Inverse einer Matrix ein weiteres zentrales Konzept, das oft mit der reduzierten Zeilen-Echelon-Form in Verbindung gebracht wird. Eine Matrix AA ist invertierbar, wenn es eine Matrix A1A^{ -1} gibt, die die Eigenschaft AA1=IA \cdot A^{ -1} = I erfüllt, wobei II die Einheitsmatrix ist.

Die Inverse einer Matrix lässt sich durch das gleiche Verfahren finden, das auch zur Lösung eines Gleichungssystems verwendet wird. Wir erweitern die Matrix AA mit der Einheitsmatrix II und führen elementare Zeilenoperationen durch, um AA in die Einheitsmatrix zu überführen. Die resultierende Matrix auf der rechten Seite der erweiterten Matrix ist dann die Inverse von AA.

Der Determinant eines Produkts von Matrizen

Ein weiteres wichtiges Ergebnis ist, dass der Determinant des Produkts zweier Matrizen gleich dem Produkt der Determinanten der einzelnen Matrizen ist. Formal ausgedrückt:

det(AB)=det(A)det(B)\det(AB) = \det(A) \cdot \det(B)

Dies gilt für quadratische Matrizen AA und BB derselben Größe. Wenn zum Beispiel eine der Matrizen nicht invertierbar ist (also det(A)=0\det(A) = 0), dann ist auch das Produkt ABAB nicht invertierbar und det(AB)=0\det(AB) = 0.

Diese Eigenschaft ist besonders nützlich, wenn man mit Matrizenoperationen arbeitet und die Invertierbarkeit von Matrizen und die Berechnung ihrer Determinanten berücksichtigt.

Wie hängen charakteristisches Polynom, Minoren und die Jordan-Normalform zusammen?

Sei f(λ)=λnb1λn1+b2λn2b3λn3++(1)nbnf(\lambda) = \lambda^n - b_1 \lambda^{n-1} + b_2 \lambda^{n-2} - b_3 \lambda^{n-3} + \cdots + (-1)^n b_n das charakteristische Polynom einer Matrix AMn(F)A \in M_n(F). Es lässt sich zeigen, dass die Koeffizienten bib_i genau die Summen der Haupt-ii-Minoren von AA sind. Diese Hauptminoren entstehen durch Auswahl von Teilmatrizen, die auf denselben Zeilen- und Spaltenindizes basieren: Genauer sind dies die i×ii \times i-Untermatrizen, die aus der Auswahl von Zeilen r1<r2<<rir_1 < r_2 < \cdots < r_i und denselben Spalten gebildet werden. Somit spiegelt das charakteristische Polynom nicht nur algebraische Eigenschaften von AA wider, sondern ist auch eng mit den Kombinatorik-Strukturen der Matrix verbunden.

Ein zentraler Schritt zum tieferen Verständnis linearer Abbildungen ist die Untersuchung der sogenannten Jordan-Normalform. Im Gegensatz zur rationalen Normalform, welche algebraisch sehr allgemein ist, bietet die Jordan-Form eine oft „simplere“ Darstellung, die eng mit den Eigenwerten und deren algebraischer Vielfachheit verbunden ist. Voraussetzung für die Existenz einer Jordan-Form ist, dass das charakteristische Polynom vollständig in lineare Faktoren über dem zugrundeliegenden Körper FF zerfällt. Dies wird durch den Begriff des algebraisch abgeschlossenen Körpers präzisiert: Ein Körper FF ist algebraisch abgeschlossen, wenn jedes Polynom über FF in lineare Faktoren zerfällt. Die komplexen Zahlen C\mathbb{C} erfüllen diese Eigenschaft, die reellen Zahlen R\mathbb{R} dagegen nicht.

Die Äquivalenz folgender Aussagen für einen linearen Endomorphismus T:VVT: V \to V auf einem endlichdimensionalen FF-Vektorraum VV ist grundlegend: Das charakteristische Polynom von TT zerfällt in lineare Faktoren über FF, genau dann wenn auch das Minimalpolynom sowie alle invarianten Faktoren in lineare Faktoren zerfallen, und alle elementaren Teiler von TT die Form (λr)e(\lambda - r)^e für rFr \in F und eZ+e \in \mathbb{Z}^+ haben.

Der Begriff des Eigenwerts rFr \in F wird durch Eigenvektoren v0v \neq 0 charakterisiert, für die Tv=rvTv = r v gilt. Die Menge aller Eigenvektoren zu einem Eigenwert rr bildet den sogenannten Eigenraum Er=ker(TrI)E_r = \ker(T - rI), der ein Untervektorraum von VV ist.

Ein wesentlicher Baustein der Jordan-Form sind sogenannte Jordan-Blöcke. In der Spezialfallbetrachtung eines zyklischen F[λ]F[\lambda]-Moduls mit Minimalpolynom (λr)e(\lambda - r)^e erhält man eine Basis, deren Darstellungsmatrix von TT die Form eines Jordan-Blocks annimmt: Ein Block mit rr auf der Hauptdiagonalen, Einsen auf der Nebendiagonalen und Nullen sonst. Diese Form veranschaulicht die verallgemeinerte Eigenstruktur, bei der neben Eigenvektoren auch verallgemeinerte Eigenvektoren eine Rolle spielen.

Die allgemeine Jordan-Normalform entsteht durch Zerlegung des Vektorraums in direkte Summen zyklischer F[λ]F[\lambda]-Untermoduln, zu denen jeweils ein Jordan-Block gehört. Die Jordan-Normalform ist dann eine Blockdiagonalmatrix aus Jordan-Blöcken, wobei die Vielfachheit und Größen der Blöcke durch die elementaren Teiler und damit durch das Minimal- und charakteristische Polynom bestimmt sind.

Für Körper, die algebraisch abgeschlossen sind, existiert für jede Matrix über diesem Körper eine Jordan-Form. Das bedeutet: Über C\mathbb{C} kann man jede lineare Abbildung durch eine Basiswahl so darstellen, dass die Matrix eine möglichst einfache, strukturierte Blockform besitzt. Im Gegensatz dazu ist über R\mathbb{R} dies nicht immer möglich, da Polynomzerlegungen nicht notwendigerweise linear sind (Beispiel: λ2+4\lambda^2 + 4 ist irreduzibel über R\mathbb{R}).

Die Jordan-Form ist für viele Anwendungen in der linearen Algebra, Differentialgleichungen und theoretischen Physik von zentraler Bedeutung, weil sie eine systematische Beschreibung der verallgemeinerten Eigenstruktur erlaubt. Die Suche nach der Jordan-Form und der zugehörigen Basis erfolgt über die Bestimmung der elementaren Teiler, des Minimalpolynoms und der Eigenräume.

Die Beispiele zu Jordan-Formen verdeutlichen auch, dass sich durch Veränderung des zugrundeliegenden Körpers die Darstellbarkeit einer Matrix erheblich ändert. Über C\mathbb{C} besitzt jede Matrix eine Jordan-Form, während über R\mathbb{R} manche Matrizen nur rationale Normalformen besitzen, aber keine Jordan-Form.

Wichtig ist zu verstehen, dass die Koeffizienten des charakteristischen Polynoms (die Summen der Hauptminoren) tief in der Struktur der Matrix verwurzelt sind und dass der Zerfall dieses Polynoms über dem Körper die Existenz der Jordan-Form bedingt. Somit verbindet die Theorie der Minoren, Polynomzerlegung und der Jordan-Form fundamentale algebraische Eigenschaften linearer Abbildungen. Dies impliziert auch, dass die Wahl des Grundkörpers eine entscheidende Rolle für die Darstellbarkeit und Struktur von Matrizen spielt.