Grundsätzliche Matrix-Operationen

für lineare Gleichungssysteme A  =
[size=100]Im CAS lege ich eine Matrix A an und verwende Einträge an denen sich leicht die Zeilen und Spalten-Positionen ablesen lassen.
Ich verwende Zahleneinträge, damit sich die Matrixelemente einfach abgreifen lassen: 
A(n) steht dann für die Matrixzeile n und A(z,s) steht für das Element in Zeile z, Spalte s![/size]
Im CAS lege ich eine Matrix A an und verwende Einträge an denen sich leicht die Zeilen und Spalten-Positionen ablesen lassen. Ich verwende Zahleneinträge, damit sich die Matrixelemente einfach abgreifen lassen: A(n) steht dann für die Matrixzeile n und A(z,s) steht für das Element in Zeile z, Spalte s!
[size=100]Tauschmatrix zum Vertauschen von Zeilen oder Spalten. Sie entstehen aus der Einheitsmatrix in dem die entsprechenden Zeilen vertauscht wurden: T12, T13, T14, T23, T24, T34 wäre ein kompletter Satz an Matrizen um alle möglichen Zeilen- oder Spalten-Vertauschungen einer 4x4 Matrix zu beschreiben.[/size]
Tauschmatrix zum Vertauschen von Zeilen oder Spalten. Sie entstehen aus der Einheitsmatrix in dem die entsprechenden Zeilen vertauscht wurden: T12, T13, T14, T23, T24, T34 wäre ein kompletter Satz an Matrizen um alle möglichen Zeilen- oder Spalten-Vertauschungen einer 4x4 Matrix zu beschreiben.
[size=100]Anwenden der Tauschmatrix T12 zum Zeilentausch[/size]
Anwenden der Tauschmatrix T12 zum Zeilentausch
[size=100]Anwenden der Tauschmatrix T12 zum Spaltentausch[/size]
Anwenden der Tauschmatrix T12 zum Spaltentausch
[size=100]Inverse: Wende ich eine Tauschmatrix auf sich selber an, dann tausche ich einen Satz vertauschter Zeilen wieder zurück und ich erhalte die Einheitsmatrix.[/size]

[size=150]Umformen von linearen Gleichungssystemen[/size]
Inverse: Wende ich eine Tauschmatrix auf sich selber an, dann tausche ich einen Satz vertauschter Zeilen wieder zurück und ich erhalte die Einheitsmatrix. Umformen von linearen Gleichungssystemen
[size=100]Das GeoGebra-Ergebnis A x = b etwas schöner aufgeschrieben:[/size]
[math]\begin{matrix}11r+12s+13u+14v\\21r+22s+23u+24v\\31r+32s+33u+34v\\41r+42s+43u+44v\end{matrix}=\begin{matrix}a\\b\\c\\d\end{matrix}[/math]
Das GeoGebra-Ergebnis A x = b etwas schöner aufgeschrieben:
[size=100]Das GeoGebra-Ergebnis T12 A x = T12 b etwas schöner aufgeschrieben:[/size]
[math]\begin{matrix}
21r+22s+23u+24v\\11r+12s+13u+14v
\\31r+32s+33u+34v
\\41r+42s+43u+44v\end{matrix}=\begin{matrix}b\\a\\c\\d\end{matrix}[/math]
Das GeoGebra-Ergebnis T12 A x = T12 b etwas schöner aufgeschrieben:
A x = b Eine Matrix-Gleichung kann einfach von Links mit einer Zeilentauschmatrix multipliziert werden um eine Zeilenoperation in der Matrix A durchzuführen. T12 A x = T12 b Um an der Matrix A eine Spaltenoperation durchzuführen kann die Matrix-Gleichung nicht einfach von rechts multipliziert werden, da die Tauschmatrix dann bei x einen Spaltentausch durchführem will wo keine Spalten sind. A x T12 = b T12 Bringen wir zwischen A x die Einheitsmatrix in Form von zwei Tauschmatrizen T12 T12 = E ein, so tauschen wir in A die Spalten und in x die Zeilen, was dann das GLS erhält - wir haben ja nur mit der Einheitsmatrix E multipliziert. Ergebnis siehe oben!
Teilbetrachtung Zeilentausch und Spaltentausch
 T12 A T12                          T12 [math]\vec{x}[/math]
[math]\left(\begin{array}{rrrr}22&21&23&24\\12&11&13&14\\32&31&33&34\\42&41&43&44\\\end{array}\right)[/math] [math]\left( \begin{matrix}s\\r\\u\\v\end{matrix} \right) [/math]
Teilbetrachtung Zeilentausch und Spaltentausch T12 A T12                T12

Zeilenumformung zur Lösung eines linearen GLS

Beim Lösen eines GLS mittels Gauss-Algorithmus oder bei der Erzeugung der LR-Zerlegung werden die Zeilen der Matrix A durch Zeilenumformungen auf die Einheitsmatrix oder die Dreiecksmatrizen L R zurückgeführt. Das Standardverfahren addiert ein Faktor q-faches der ersten Zeile zur allen folgenden Zeilen, so dass in der ersten Spalte Nullen entstehen: Die Zeilenumformungsmatrix T2plus1 soll ein q faches der ersten Zeile zur zweiten und ein p-faches der zweiten Zeile zur dritten addieren.
[size=100]Allgemein werden in Zeilentransformationsmatrizen die Faktoren für die erste Zeile unterhalb der Diagonalen in der ersten Spalte und die Faktoren für die zweite Zeile in der zweiten Spalte hinterlegt. Weitere Zeilenfaktoren analog eintragen ...[/size]
Allgemein werden in Zeilentransformationsmatrizen die Faktoren für die erste Zeile unterhalb der Diagonalen in der ersten Spalte und die Faktoren für die zweite Zeile in der zweiten Spalte hinterlegt. Weitere Zeilenfaktoren analog eintragen ...
[size=100]Im Ergebnis erhalten wir Zeile2=q*Zeile1+Zeile2 und Zeile3=p*Zeile2+Zeile3
Die Zeilenumformung Zeile1+Zeile2 soll im Lösungsschritt Gauss/LR-Zerlegung in der ersten Spalte 0 ergeben:
11q + 21 = 0
q = -21/11 oder allgemein q = -A(2,1)/A(1,1)[/size]
Je nach Versionszustand muss die Matrixindizierung auch via Befehl ELEMENT erfolgen: A(2,1) = ELEMENT(A,2,1)!
Im Ergebnis erhalten wir Zeile2=q*Zeile1+Zeile2 und Zeile3=p*Zeile2+Zeile3 Die Zeilenumformung Zeile1+Zeile2 soll im Lösungsschritt Gauss/LR-Zerlegung in der ersten Spalte 0 ergeben: 11q + 21 = 0 q = -21/11 oder allgemein q = -A(2,1)/A(1,1) Je nach Versionszustand muss die Matrixindizierung auch via Befehl ELEMENT erfolgen: A(2,1) = ELEMENT(A,2,1)!
[size=100]Den Vorzeichenwechsel der Matrix-Elemente unter der Diagonalen lässt sich rechnerisch darstellen: die Einheitsmatrix mit 2 multiplizieren und dann die Matrix L1 abziehen: 
L1^-1 = 2 E - L1 
[/size]
Den Vorzeichenwechsel der Matrix-Elemente unter der Diagonalen lässt sich rechnerisch darstellen: die Einheitsmatrix mit 2 multiplizieren und dann die Matrix L1 abziehen: L1^-1 = 2 E - L1

Zusammenfassung

In den diskutierten Beispielen habe ich einige Matrix-Operationen und Eigenschaften der verwendeten Matrizen zusammen gestellt, die auf dem Weg zur Lösung von linearen GLS zur Anwendung kommen. Ich beginne eine neues Arbeitsblatt um ein konkretes GLS mittels Gauss-Algorithmus bzw. mittels LR-Zerlegung zu lösen. Siehe KonvertGLS2Matrix.ggb, um aus einem in x1,x2,x3... gegebenen GLS die Matrixgleichung A x = b heraus zu ziehen.

Rang, Kern in-,sur-,bijektiv

Rang der Matrix = Anzahl der Spalten ⇔ Abbildungs-Matrix ist injektiv. Rang der Matrix = Anzahl der Zeilen ⇔ Abbildungs-Matrix ist surjektiv. Kern der Matrix = 0 ⇔ Abbildungs-Matrix ist injektiv. Bei quadratischen Abbildungs-Matrizen folgt aus der Surjketivität die Injektivitiät und umgekehrt aus der Injektivität die Surjektivität. Wenn die Determinante einer quadratischen Matrix ≠0 ist, kann man die Matrix invertieren, also gibt es eine Umkehrabbildung. Die Abbildungs-Matrix ist dann bijektiv. Das gilt auch umgekehrt.