Reguläre Matrix

quadratische Matrix, die eine Inverse besitzt

Eine reguläre, invertierbare oder nichtsinguläre Matrix ist in der Mathematik eine quadratische Matrix, die eine Inverse besitzt. Reguläre Matrizen können auf mehrere äquivalente Weisen charakterisiert werden. Zum Beispiel zeichnen sich reguläre Matrizen dadurch aus, dass die durch sie beschriebene lineare Abbildung bijektiv ist. Daher ist ein lineares Gleichungssystem mit einer regulären Koeffizientenmatrix stets eindeutig lösbar. Die Menge der regulären Matrizen fester Größe mit Einträgen aus einem Ring oder Körper bildet mit der Matrizenmultiplikation als Verknüpfung die allgemeine lineare Gruppe.

Nicht zu jeder quadratischen Matrix existiert eine Inverse. Eine quadratische Matrix, die keine Inverse besitzt, wird singuläre Matrix genannt.

Definition

Bearbeiten

Eine quadratische Matrix   mit Einträgen aus einem unitären Ring   (in der Praxis meist dem Körper der reellen Zahlen) heißt regulär, wenn eine weitere Matrix   existiert, sodass

 

gilt, wobei   die Einheitsmatrix bezeichnet. Die Matrix   ist hierbei eindeutig bestimmt und heißt inverse Matrix zu  . Die Inverse einer Matrix   wird üblicherweise mit   bezeichnet. Bei einer singulären Matrix existiert keine solche Matrix  .

Ist   ein kommutativer Ring, Körper oder Schiefkörper, so sind die beiden Bedingungen   und   äquivalent, das heißt, eine linksinverse Matrix ist dann auch rechtsinvers und umgekehrt, sprich, die obige Bedingung lässt sich durch   beziehungsweise   abschwächen.

Beispiele

Bearbeiten

Die reelle Matrix

 

ist regulär, denn sie besitzt die Inverse

 ,

mit

 .

Die reelle Matrix

 

ist singulär, denn für eine beliebige Matrix

 

gilt

 .

Äquivalente Charakterisierungen

Bearbeiten

Reguläre Matrizen über einem Körper

Bearbeiten

Eine  -Matrix   mit Einträgen aus einem Körper  , zum Beispiel die reellen oder komplexen Zahlen, ist genau dann invertierbar, wenn eine der folgenden äquivalenten Bedingungen erfüllt ist:

  • Es gibt eine Matrix   mit  .
  • Die Determinante von   ist ungleich null:  .
  • Die Null ist kein Eigenwert von  .
  • Das lineare Gleichungssystem   besitzt nur die triviale Lösung  .
  • Für jedes   existiert mindestens eine Lösung   des linearen Gleichungssystems  .
  • Für jedes   existiert höchstens eine Lösung   des linearen Gleichungssystems  .
  • Die Zeilenvektoren sind linear unabhängig.
  • Die Zeilenvektoren erzeugen  .
  • Die Spaltenvektoren sind linear unabhängig.
  • Die Spaltenvektoren erzeugen  .
  • Die durch   beschriebene lineare Abbildung  ,  , ist bijektiv.
  • Die transponierte Matrix   ist invertierbar.
  • Der Rang der Matrix   ist gleich  .

Bei einer singulären  -Matrix   mit Einträgen aus einem Körper   ist keine der obigen Bedingungen erfüllt.

Reguläre Matrizen über einem unitären kommutativen Ring

Bearbeiten

Allgemeiner ist eine  -Matrix   mit Einträgen aus einem kommutativen Ring mit Eins   genau dann invertierbar, wenn eine der folgenden äquivalenten Bedingungen erfüllt ist:

  • Es gibt eine Matrix   mit  .
  • Die Determinante von   ist eine Einheit in   (man spricht auch von einer unimodularen Matrix).
  • Für alle   existiert genau eine Lösung   des linearen Gleichungssystems  .
  • Für alle   existiert mindestens eine Lösung   des linearen Gleichungssystems  .
  • Die Zeilenvektoren bilden eine Basis von  .
  • Die Zeilenvektoren erzeugen  .
  • Die Spaltenvektoren bilden eine Basis von  .
  • Die Spaltenvektoren erzeugen  .
  • Die durch   beschriebene lineare Abbildung  ,  , ist surjektiv (oder gar bijektiv).
  • Die transponierte Matrix   ist invertierbar.

Bei einer singulären  -Matrix   mit Einträgen aus einem kommutativen Ring mit Eins   ist keine der obigen Bedingungen erfüllt.

Der wesentliche Unterschied zum Fall eines Körpers ist hier also, dass im Allgemeinen aus der Injektivität einer linearen Abbildung nicht mehr ihre Surjektivität (und damit ihre Bijektivität) folgt, wie bereits das einfache Beispiel  ,   zeigt.

Weitere Beispiele

Bearbeiten

Die Matrix

 

mit Einträgen aus dem Polynomring   hat die Determinante   und   ist invertierbar in  . Somit ist   regulär in  . Die inverse Matrix ist

 .

Die Matrix

 

mit Einträgen aus dem Restklassenkörper   hat die Determinante   und   ist invertierbar in  . Somit ist   regulär in  . Die inverse Matrix ist

 .

Die Matrix

 

mit Einträgen aus dem Restklassenring   hat die Determinante  . Da   und   nicht teilerfremd sind, ist   in   nicht invertierbar. Daher ist   nicht regulär.

Eigenschaften

Bearbeiten

Ist die Matrix   regulär, so ist auch   regulär mit der Inversen

 .

Sind die beiden Matrizen   und   regulär, so ist auch ihr Produkt   regulär mit der Inversen

 .

Die Menge der regulären Matrizen fester Größe bildet demnach mit der Matrizenmultiplikation als Verknüpfung eine (im Allgemeinen nichtkommutative) Gruppe, die allgemeine lineare Gruppe  . In dieser Gruppe ist die Einheitsmatrix das neutrale Element und die inverse Matrix das inverse Element. Für eine reguläre Matrix   gelten damit auch die Kürzungsregeln

 

und

 ,

wobei   und   beliebige Matrizen passender Größe sind.

Eine singuläre Matrix besitzt den Eigenwert null, d. h., es gibt einen vom Nullvektor verschiedenen Vektor, der von der Matrix auf ersteren abgebildet wird. Alle Vektoren, die von der Matrix auf den Nullvektor abgebildet werden, erzeugen den Eigenraum zum Eigenwert null. Die Dimension dieses Eigenraumes ist die geometrische Vielfachheit des Eigenwerts null, siehe Jänich (2008), S. 197 ff.

Blockmatrizen

Bearbeiten

Ist eine quadratische Blockmatrix   gegeben, wobei   und das Schur-Komplement   von   in   eine reguläre Matrix ist, dann ist auch   eine reguläre Matrix und es gilt

 

Daraus folgt für die inverse Matrix

 

Wenn   und das Schur-Komplement   von   in   eine reguläre Matrix ist, gilt entsprechend

 

und für die inverse Matrix[1]

 

Mithilfe dieser Formel kann die inverse Matrix einer quadratischen ( )-Blockmatrix  mit Blöcken der Dimension   effizient berechnet werden. Es ist also  . Die Laufzeit für die Inversion beträgt  . Im Vergleich dazu beträgt die Laufzeit für den Gauß-Jordan-Algorithmus  .[2]

Reguläre Matrizen über einem Restklassenkörper

Bearbeiten

Eine Matrix mit Einträgen aus einem Restklassenkörper   mit einer Primzahl   ist genau dann regulär, wenn die Zeilenvektoren linear unabhängig sind.

Für den Restklassenkörper   kann die Anzahl der regulären  -Matrixen wie folgt berechnet werden:

  • Jedes der   Elemente der 1. Zeile kann unabhängig voneinander 2 Werte annehmen. Der Nullvektor ist ausgeschlossen. Für die 1. Zeile gibt es also   Möglichkeiten.
  • Für die 2. Zeile sind alle Vektoren ausgeschlossen, die eine Linearkombination der 1. Zeile sind, also   Vektoren. Für die 2. Zeile gibt es also   Möglichkeiten.
  • Für die 3. Zeile sind alle Vektoren ausgeschlossen, die eine Linearkombination der 1. Zeile und 2. Zeile sind, also   Vektoren. Für die 3. Zeile gibt es also   Möglichkeiten.
  • Allgemein gibt es für die Zeile mit dem Index   also   mögliche Werte. Für alle Zeilen der Matrix ergeben sich daher insgesamt   Möglichkeiten.

Daraus lässt sich der Anteil der regulären  -Matrixen an allen  -Matrixen bestimmen. Es gibt   verschiedene  -Matrixen, weil jedes der   Elemente unabhängig voneinander 2 Werte annehmen kann. Der Anteil der regulären  -Matrixen beträgt daher

 

Für   gegen unendlich konvergiert dieses Produkt nach dem Pentagonalzahlensatz wegen   gegen einen endlichen Grenzwert. Dieser beträgt etwa 0,289.

Dieses Ergebnis lässt sich für beliebige Primzahlen   auf den Restklassenkörper   verallgemeinern. Es gibt   verschiedene  -Matrixen, von denen   reguläre  -Matrixen sind. Der Anteil der regulären  -Matrixen beträgt  .[3]

Literatur

Bearbeiten
Bearbeiten

Einzelnachweise

Bearbeiten
  1. Stephen M. Watt, University of Western Ontario: Pivot-Free Block Matrix Inversion
  2. Iria C. S. Cosme, Isaac F. Fernandes, Joao L. de Carvalho, Samuel Xavier-de-Souza: Memory-Usage Advantageous Block Recursive Matrix Inverse
  3. StackExchange: Number of non singular matrices over a finite field of order 2