Matrix (Mathematik)

In der linearen Algebra ist eine Matrix (pl.: Matrizen) eine 2-dimensionale Anordnung von Zahlenwerten (aber auch anderen Objekten wie Operatoren) in Tabellenform. Man spricht von den Spalten und Zeilen der Matrix, und bezeichnet selbige auch als Vektoren (d.h. Zeilenvektoren und Spaltenvektoren). Die Objekte, die in der Matrix angeordnet sind, nennt man Komponenten oder Elemente der Matrix.

Wenn die Matrix m Zeilen und n Spalten besitzt, spricht man von einer m × n-Matrix, und nennt m und n die Dimensionen der Matrix. Ist m = n spricht man von einer n × n oder quadratischen Matrix. Die Komponente, die in der i-ten Zeile an j-ter Stelle steht, hat die Indices i,j. Eine allgemeine 2 × 3 Matrix A sieht zum Beispiel so aus:

${\begin{pmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\end{pmatrix}}=(a_{ij})_{i=1\ldots 2,j=1\ldots 3}$

Die Menge aller m × n-Matrizen über einer Menge K bezeichnet man mit K^m×n oder K^m,n, selten auch mit ^mKⁿ.

Hat eine Matrix nur eine einzige Spalte oder Zeile, dann nennt man sie Vektor. Man unterscheidet dabei zwischen einem Zeilenvektor (mit nur einer Zeile) oder einem Spaltenvektor (mit nur einer Spalte). Oft benötigt man diese Unterscheidung nicht (wenn man sie nicht mit Matrizen multiplizieren muss) und bezeichnet die Menge aller n-stelligen Vektoren über einer Menge K mit Kⁿ statt K^1×n oder K^n×1.

Im Zusammenhang mit Matrizen oft auftretenden Begriffe sind der Rang und die Determinante einer Matrix.

"Kippt" man die Matrix A an der Hauptdiagonalen (die von links oben nach rechts unten geht, ihre Komponenten sind a_ii), dann erhält man die zu A transponierte Matrix A^T.

Zwei Matrizen A und B gleicher Dimension mit Komponenten in einer Zahlenmenge (z.B. den reellen Zahlen) kann man komponentenweise addieren. Stimmt dagegen die Zeilenanzahl von A mit der Spaltenanzahl von B überein, dann kann man das Matrixprodukt A*B berechnen. Siehe dazu die folgenden Beispiele.

Rechnen mit Matrizen

Addieren von Matrizen

Zwei Matrizen A und B werden addiert, indem man die in den Matrizen an entsprechender Stelle stehenden Komponenten addiert:

${\begin{pmatrix}1&3&2\\1&0&0\\1&2&2\end{pmatrix}}+{\begin{pmatrix}0&0&5\\7&5&0\\2&1&1\end{pmatrix}}={\begin{pmatrix}1+0&3+0&2+5\\1+7&0+5&0+0\\1+2&2+1&2+1\end{pmatrix}}={\begin{pmatrix}1&3&7\\8&5&0\\3&3&3\end{pmatrix}}$

Genauso werden auch Vektoren addiert.

${\begin{pmatrix}1\\2\\9\end{pmatrix}}+{\begin{pmatrix}3\\7\\2\end{pmatrix}}={\begin{pmatrix}1+3\\2+7\\9+2\end{pmatrix}}={\begin{pmatrix}4\\9\\11\end{pmatrix}}$

${\begin{pmatrix}1&2&9\end{pmatrix}}+{\begin{pmatrix}3&7&2\end{pmatrix}}={\begin{pmatrix}1+3&2+7&9+2\end{pmatrix}}={\begin{pmatrix}4&9&11\end{pmatrix}}$

Die Matritzen müssen die selben Dimensionen aufweisen.

Vervielfachen von Matrizen (Skalarmultiplikation)

Eine Matrix oder ein Vektor A wird mit einer Zahl r vervielfacht (multipliziert), indem man jede Komponente von A mit r multipliziert:

$2\cdot {\begin{pmatrix}1&8&-3\\4&-2&5\end{pmatrix}}={\begin{pmatrix}2\cdot 1&2\cdot 8&2\cdot (-3)\\2\cdot 4&2\cdot (-2)&2\cdot 5\end{pmatrix}}={\begin{pmatrix}2&16&-6\\8&-4&10\end{pmatrix}}$

Diese Rechenoperation nennt man Skalarmultiplikation, das Ergebnis ist ein skalares Produkt, es ist zu unterscheiden vom Skalarprodukt zweier Vektoren.

Multiplizieren von Matrizen

Zwei Matrizen A und B werden miteinander multipliziert, indem jeweils die Zeilenelemente der ersten Matrix mit den entsprechenden Spaltenelementen der zweiten Matrix multipliziert werden. Die Multiplikation von Matrizen ist nur dann möglich, wenn die Länge der Zeilen (= die Anzahl der Spalten) der ersten Matrix mit der Länge der Spalten (= Anzahl der Zeilen) der zweiten Matrix übereinstimmt. Ist A eine l×m-Matrix und B eine m×n-Matrix, so ist das Produkt eine l×n-Matrix. Aber Achtung: Bei Matrizen gilt Kommutativgesetz NICHT!

${\begin{pmatrix}1&2&3\\4&5&6\\\end{pmatrix}}\cdot {\begin{pmatrix}6&-1\\3&2\\0&-3\end{pmatrix}}={\begin{pmatrix}(1\cdot 6+2\cdot 3+3\cdot 0)&(1\cdot -1+2\cdot 2+3\cdot -3)\\(4\cdot 6+5\cdot 3+6\cdot 0)&(4\cdot -1+5\cdot 2+6\cdot -3)\\\end{pmatrix}}={\begin{pmatrix}12&-6\\39&-12\end{pmatrix}}$

Dabei können A und B auch Vektoren sein, solange die Formate passen (siehe dazu auch den Abschnitt Vektor-Vektor-Produkte).

Formal definiert ist die Matrixmultiplikation für Matrizen A = (a_ij) und B = (b_ij) der Formate l×m und m×n als die Matrix C = (c_ij) mit den Komponenten

$c_{ij}=\sum _{k=1}^{m}a_{ik}b_{kj}\quad i=1,\ldots l,j=1,\ldots n$

Die transponierte Matrix

Die Transponierte der Matrix A = (a_ij) vom Format m×n ist die Matrix A^T = (a_ji) vom Format n×m.

${\begin{pmatrix}1&8&-3\\4&-2&5\end{pmatrix}}^{T}={\begin{pmatrix}1&4\\8&-2\\-3&5\end{pmatrix}}$

Die inverse Matrix

Ist (R,+,*,0) ein Ring, dann bildet die Menge R^n×n der quadratischen Matrizen vom Format n×n ebenfalls einen Ring mit der oben definierten Matrixaddition und -multiplikation. Das Nullelement ist die Nullmatrix 0, deren Komponenten alle 0 sind. Hat R ein Einselement 1, dann ist die Einheitsmatrix E das Einselement des Matrixrings. Sie hat auf der Hauptdiagonalen 1 und sonst 0.

Ist K ein Körper, dann sind im Ring K^n×n genau diejenigen Matrizen invertierbar (regulär), deren Determinante ungleich 0 ist. Man kann die zur Matrix A inverse Matrix A^-1 zum Beispiel mit dem Gauss-Algorithmus bestimmen. Dazu löst man das lineare Gleichungssystem A*X = E. Die Matrix E ist die Einheitsmatrix, die Matrix X ist dann das Inverse von A.

${\begin{pmatrix}1&2&3\\1&1&4\\1&4&5\end{pmatrix}}\cdot {\begin{pmatrix}a&b&c\\d&e&f\\g&h&i\end{pmatrix}}={\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}}$

Es entsteht ein LGS mit 9 Unbekannten und 9 Gleichungen:

${\begin{matrix}1a+1b+1c+0d+0e+0f+0g+0h+0i=1\\2a+1b+4c+0d+0e+0f+0g+0h+0i=0\\3a+4b+5c+0d+0e+0f+0g+0h+0i=0\\0a+0b+0c+1d+1e+1f+0g+0h+0i=0\\0a+0b+0c+2d+1e+4f+0g+0h+0i=1\\0a+0b+0c+3d+4e+5f+0g+0h+0i=0\\0a+0b+0c+0d+0e+0f+1g+1h+1i=0\\0a+0b+0c+0d+0e+0f+2g+1h+4i=0\\0a+0b+0c+0d+0e+0f+3g+4h+5i=1\\\end{matrix}}$

Bei näherer Betrachtung stellt man aber sehr schnell fest, dass man dieses Gleichungssystem immer auch als drei getrennte Gleichungssysteme mit je drei Gleichungen und Unbekannte zerlegen kann:

${\begin{matrix}1a+1b+1c=1\\2a+1b+4c=0\\3a+4b+5c=0\\\end{matrix}}$

${\begin{matrix}1d+1e+1f=0\\2d+1e+4f=1\\3d+4e+5f=0\\\end{matrix}}$

${\begin{matrix}1g+1h+1i=0\\2g+1h+4i=0\\3g+4h+5i=1\\\end{matrix}}$

Es muss also dreimal dasselbe Gleichungssystem mit unterschiedlichen rechten Seiten (nämlich den drei Einheitsvektoren) gelöst werden. Dies gilt analog bei Matrizen höherer Dimension. Effizient geht dies über elementare Zeilenumformungen, also über das Gauß-Verfahren. Das Ergebnis sollte lauten:

${\begin{pmatrix}a&b&c\\d&e&f\\g&h&i\end{pmatrix}}={\begin{pmatrix}2,75&-0,5&-1,25\\0,25&-0,5&0,25\\-0,75&0,5&0,25\end{pmatrix}}$

Berechnung der Inversen mittels elementarer Zeilenumformungen:

Wir schreiben rechts neben die zu invertierende Matrix die Einheitsmatrix, z.B.

${\begin{pmatrix}1&2&3&|&1&0&0\\4&5&6&|&0&1&0\\7&8&7&|&0&0&1\end{pmatrix}}$

Jetzt formen wir die Matrix so lange mit Elementaren Zeilenumformungen um, bis auf der linken Seite die Einheitsmatrix steht. Das ist nicht besonders kompliziert, im oben angegebenen Beispiel geht das z.B. folgendermaßen:

Im ersten Schritt ziehen wir dazu die erste Zeile 4 mal von der zweiten ab, und 7 mal von der dritten:

${\begin{pmatrix}1&2&3&|&1&0&0\\0&-3&-6&|&-4&1&0\\0&-6&-14&|&-7&0&1\end{pmatrix}}$

Jetzt ziehen wir die zweite Zeile 2 mal von der dritten ab, und erhalten so links eine obere Dreiecksmatrix:

${\begin{pmatrix}1&2&3&|&1&0&0\\0&-3&-6&|&-4&1&0\\0&0&-2&|&1&-2&1\end{pmatrix}}$

Nun normieren wir die Zeilen um auf der Diagonalen 1er zu erhalten, indem wir die zweite mit -1/3 Multiplizieren, die dritte mit -1/2

${\begin{pmatrix}1&2&3&|&1&0&0\\0&1&2&|&{\frac {4}{3}}&-{\frac {1}{3}}&0\\0&0&1&|&-{\frac {1}{2}}&1&-{\frac {1}{2}}\end{pmatrix}}$

Von der ersten Zeile ziehen wir 2 mal die zweite ab und addieren die dritte 1 mal, anschließend ziehen wir die dritte Zeile 2 mal von der zweiten ab:

${\begin{pmatrix}1&0&0&|&-{\frac {13}{6}}&{\frac {5}{3}}&-{\frac {1}{2}}\\0&1&0&|&{\frac {7}{3}}&-{\frac {7}{3}}&1\\0&0&1&|&-{\frac {1}{2}}&1&-{\frac {1}{2}}\end{pmatrix}}$

Auf der rechten Seite steht jetzt die Inverse der ursprünglichen Matrix.

Welche elementaren Zeilenumformungen man verwendet ist hierbei egal, es empfielt sich jedoch ein zielgerichtetes Arbeiten, wie eben gezeigt (Spalte für Spalte auf obere Dreiecksmatrix bringen, anschließend die Diagonale auf 1en normieren. Anschließend ist es meistens am einfachsten, wenn man (anders als eben gezeigt) von unten anfangend die Einheitsmatrix herstellt. (Im eben gezeigten Beispiel hätte man dazu statt dem letzten Schritt erst die letzte Zeile 2 mal von der vorletzten abgezogen, anschließend die letzte 3 mal und die mittlere Zeile 2 mal von der ersten)

Vektor-Vektor-Produkte (Skalarprodukt)

Hat man zwei Spaltenvektoren v und w der Länge n, dann ist das Matrixprodukt v*w nicht definiert, aber die beiden Produkte v^T*w und v*w^T existieren.

Das erste Produkt ist eine 1×1-Matrix, die als Zahl interpretiert wird, sie wird das kanonische Skalarprodukt von v und w genannt und mit <v,w> bezeichnet.

${\begin{pmatrix}1\\2\\3\end{pmatrix}}^{T}\cdot {\begin{pmatrix}-2\\-1\\1\end{pmatrix}}=1\cdot (-2)+2\cdot (-1)+3\cdot 1=3$

Das zweite Produkt ist eine n×n-Matrix und heißt das dyadische Produkt von v und w.

${\begin{pmatrix}1\\2\\3\end{pmatrix}}\cdot {\begin{pmatrix}-2\\-1\\1\end{pmatrix}}^{T}={\begin{pmatrix}1\cdot (-2)&1\cdot (-1)&1\cdot 1\\2\cdot (-2)&2\cdot (-1)&2\cdot 1\\3\cdot (-2)&3\cdot (-1)&3\cdot 1\end{pmatrix}}={\begin{pmatrix}-2&-1&1\\-4&-2&2\\-6&-3&3\end{pmatrix}}$

Umformen von Matrizengleichungen

Speziell in den Multivariaten Verfahren werden häufig Beweisführungen, Herleitungen usw. im Matrizenkalkül durchgeführt.

Gleichungen werden im Prinzip wie algebraische Gleichungen umgeformt. Es wird von der Gleichung

$X=Y$

ausgegangen mit X,Y als n×m-Matrix.

Addieren und Subtrahieren

Die Gleichung kann von links und rechts mit einer n×m-Matrix A additiv erweitert werden zu

$X+A=Y+A$

bzw.

$-A+X=Y-A$ .

Multiplizieren mit einer Matrix

Die Gleichung kann multiplikativ von links durch die r×n-Matrix A oder von rechts durch die m×s-Matrix B erweitert werden:

$AX=AY$

bzw.

$XB=YB$

"Division" durch eine Matrix A

Die Gleichung wird mit der Inversen der Matrix A multipliziert, wobei A invertierbar sein muss.

Beispiele

Lineares Gleichungssystem als einfache Umformung

Gesucht ist der Lösungsvektor x eines linearen Gleichungssystems

$Ax=b$

mit A als n×m-Koeffizientenmatrix. Man erweitert von links

$A^{-1}\cdot A\cdot x=A^{-1}\cdot b\Leftrightarrow E\cdot x=A^{-1}\cdot b$

und erhält die Lösung

$x=A^{-1}\cdot b$.

Eine etwas aufwendigere Umformung erfordert ein

Orthogonalitätsbeweis im Multiplen Regressionsmodell

Im Multiplen Regressionsmodell geht man davon aus, dass eine abhängige Variable y durch p vorgegebene Variablen x_j (j=1,...,p) erklärt werden kann. Man schätzt mit Hilfe der Methode der kleinsten Quadrate die Regressionskoeffizienten als

$b=(X^{T}X)^{-1}X^{T}y$

und erhält mit der Schätzung das Gleichungssystem

$y=Xb+d$

mit y als (n×1)-Vektor der n y-Werte, d als (n×1)-Vektor der Residuen, b als ((p+1)×1)-Vektor der Regressionskoeffizienten und der (n× (p+1))-Datenmatrix X.

Es soll gezeigt werden, dass d^TXb = 0 ist: Es ist zunächst

$d^{T}Xb=(y-Xb)^{T}Xb=(y^{T}Xb-b^{T}X^{T}Xb)$

Wegen

$Xb=X(X^{T}X)^{-1}X^{T}y$

mit

$M:=X(X^{T}X)^{-1}X^{T}$

als idempotenter Matrix, d.h.

$M\cdot M=M$ .

erhält man

$d^{T}Xb=y^{T}My-y^{T}MMy$

was wegen MM=M das Skalar 0 ergibt.

Die Regressionshyperebene Xb steht auf d senkrecht, was man dahingehend interpretieren kann, dass in den Residuen keine verwertbare Information von X bzw. Xb mehr enthalten ist.

Verallgemeinerungen

Man kann auch Matrizen mit unendlich vielen Spalten oder Zeilen betrachten. Diese kann man immer noch addieren. Um sie jedoch multiplizieren zu können, muss man zusätzliche Bedingungen an ihre Komponenten stellen (da die auftretenden Summen unendliche Reihen sind und nicht konvergieren müssten).

Lässt man mehr als zwei Indices zu, erhält man Strukturen, die man sich als drei- oder höherdimensionale Tabellen denken kann. Diese Strukturen sind spezielle Tensoren.