Cramersche Regel
# Cramersche Regel

Die Cramersche Regel oder Determinantenmethode ist eine mathematische Formel für die Lösung eines linearen Gleichungssystems. Sie ist bei der theoretischen Betrachtung linearer Gleichungssysteme hilfreich. Für die Berechnung einer Lösung ist der Rechenaufwand jedoch in der Regel zu hoch, da dabei verhältnismäßig viele Determinanten auftreten. Deshalb kommen dazu andere Verfahren aus der numerischen Mathematik zum Einsatz. Die Cramersche Regel ist nach Gabriel Cramer benannt, der sie im Jahr 1750 veröffentlichte, jedoch wurde sie bereits vorher von Leibniz entdeckt.

## Regel

Gegeben sei ein lineares Gleichungssystem mit gleich vielen Gleichungen wie Unbekannten in der Form
    ${\begin{matrix}a_{11}x_{1}+a_{12}x_{2}\,+&\cdots &+\,a_{1n}x_{n}&=&b_{1}\\a_{21}x_{1}+a_{22}x_{2}\,+&\cdots &+\,a_{2n}x_{n}&=&b_{2}\\&&&\vdots &\\a_{n1}x_{1}+a_{n2}x_{2}\,+&\cdots &+\,a_{nn}x_{n}&=&b_{n}\end{matrix}}$
bzw. in Matrixschreibweise $Ax=b$ mit
    $A={\begin{pmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{n1}&a_{n2}&\cdots &a_{nn}\end{pmatrix}}\,,\quad x={\begin{pmatrix}x_{1}\\x_{2}\\\vdots \\x_{n}\end{pmatrix}}\,,\quad b={\begin{pmatrix}b_{1}\\b_{2}\\\vdots \\b_{n}\end{pmatrix}}\,.$
Vorausgesetzt sei außerdem, dass die quadratische Koeffizientenmatrix $A$ regulär (invertierbar) ist. Dies ist genau dann der Fall, wenn $\det(A)\neq 0$ ist.
Dann ist das Gleichungssystem eindeutig lösbar und die Komponenten $x_{i}$ des eindeutig bestimmten Lösungsvektors $x$ sind gegeben durch
    $x_{i}={\frac {\det(A_{i})}{\det(A)}}$ für alle $i$.
Hierbei ist $A_{i}$ die Matrix, die gebildet wird, indem die $i$-te Spalte der Koeffizientenmatrix $A$ durch die rechte Seite des Gleichungssystems $b$ ersetzt wird:
    $A_{i}={\begin{pmatrix}a_{1,1}&\cdots &a_{1,i-1}&b_{1}&a_{1,i+1}&\cdots &a_{1,n}\\a_{2,1}&\cdots &a_{2,i-1}&b_{2}&a_{2,i+1}&\cdots &a_{2,n}\\\vdots &&\vdots &\vdots &\vdots &&\vdots \\a_{n,1}&\cdots &a_{n,i-1}&b_{n}&a_{n,i+1}&\cdots &a_{n,n}\end{pmatrix}}$

### Begründung

Die Regel ergibt sich aus den zwei Regeln, denen die Determinante gehorcht:
* Sie ist multilinear in den Spalten (und Zeilen) einer Matrix, das heißt linear in jeder einzelnen Spalte (und Zeile).
* Sie ist alternierend in den Spalten (und Zeilen) einer Matrix, das heißt: Sind irgend zwei Spalten (oder Zeilen) gleich, so verschwindet die Determinante.
Beachtet man nun, dass der Vektor $b$ gemäß Gleichung $Ax=b$ gerade eine Linearkombination der Spalten der Matrix $A$ (jeweils mit Skalarfaktor $x_{i}$ versehen) ist, so wird deutlich, dass bei der Berechnung von $\det A_{i}$ gemäß den beiden Regeln allein diejenige Spalte von $A$ einen nicht verschwindenden Beitrag liefert, an deren Stelle der Vektor $b$ getreten ist, und dieser Beitrag besteht gerade aus dem $x_{i}$-Fachen der $i$-ten Spalte, so dass $\det A_{i}=x_{i}\cdot \det A$, wie die Cramersche Regel besagt.

## Beispiele
### Lineares Gleichungssystem 2. Ordnung

Diesem Beispiel liegt das folgende lineare Gleichungssystem zu Grunde:
    ${\begin{matrix}\color {blue}{1}\,\color {black}x_{1}+\color {blue}{2}\,\color {black}x_{2}=\color {green}{3}\\\color {blue}{4}\,\color {black}x_{1}+\color {blue}{5}\,\color {black}x_{2}=\color {green}{6}\end{matrix}}$
Die erweiterte Koeffizientenmatrix des Gleichungssystems ist dann:
    ${\begin{pmatrix}\color {blue}{A}&\color {green}b\end{pmatrix}}=\left({\begin{array}{cc|c}\color {blue}{1}&\color {blue}{2}&\color {green}{3}\\\color {blue}{4}&\color {blue}{5}&\color {green}{6}\end{array}}\right)$
Nach der Cramerschen Regel berechnet sich dessen Lösung wie folgt:
    $x_{1}={\frac {\det(A_{1})}{\det(A)}}={\frac {\det {\begin{pmatrix}\color {green}{3}&\color {blue}{2}\\\color {green}{6}&\color {blue}{5}\end{pmatrix}}}{\det {\begin{pmatrix}\color {blue}{1}&\color {blue}{2}\\\color {blue}{4}&\color {blue}{5}\end{pmatrix}}}}={\frac {3}{-3}}=-1\qquad$
    $x_{2}={\frac {\det(A_{2})}{\det(A)}}={\frac {\det {\begin{pmatrix}\color {blue}{1}&\color {green}{3}\\\color {blue}{4}&\color {green}{6}\end{pmatrix}}}{\det {\begin{pmatrix}\color {blue}{1}&\color {blue}{2}\\\color {blue}{4}&\color {blue}{5}\end{pmatrix}}}}={\frac {-6}{-3}}=2$

### Lineares Gleichungssystem 3. Ordnung

Diesem Beispiel liegt das folgende lineare Gleichungssystem zu Grunde:
    ${\begin{matrix}\color {blue}{82}\,\color {black}x_{1}+\color {blue}{45}\,\color {black}x_{2}+\color {blue}{9}\,\color {black}x_{3}=\color {green}{1}\\\color {blue}{27}\,\color {black}x_{1}+\color {blue}{16}\,\color {black}x_{2}+\color {blue}{3}\,\color {black}x_{3}=\color {green}{1}\\\color {blue}{9}\,\color {black}x_{1}+\color {blue}{5}\,\color {black}x_{2}+\color {blue}{1}\,\color {black}x_{3}=\color {green}{0}\\\end{matrix}}$
Die erweiterte Koeffizientenmatrix des Gleichungssystems ist dann:
    ${\begin{pmatrix}\color {blue}{A}&\color {green}b\end{pmatrix}}=\left({\begin{array}{ccc|c}\color {blue}{82}&\color {blue}{45}&\color {blue}{9}&\color {green}{1}\\\color {blue}{27}&\color {blue}{16}&\color {blue}{3}&\color {green}{1}\\\color {blue}{9}&\color {blue}{5}&\color {blue}{1}&\color {green}{0}\end{array}}\right)$
Nach der Cramerschen Regel berechnet sich die Lösung des Gleichungssystems wie folgt:
    $x_{1}={\frac {\det(A_{1})}{\det(A)}}={\frac {\det {\begin{pmatrix}\color {green}{1}&\color {blue}{45}&\color {blue}{9}\\\color {green}{1}&\color {blue}{16}&\color {blue}{3}\\\color {green}{0}&\color {blue}{5}&\color {blue}{1}\end{pmatrix}}}{\det {\begin{pmatrix}\color {blue}{82}&\color {blue}{45}&\color {blue}{9}\\\color {blue}{27}&\color {blue}{16}&\color {blue}{3}\\\color {blue}{9}&\color {blue}{5}&\color {blue}{1}\end{pmatrix}}}}={\frac {1}{1}}=1\qquad$
    $x_{2}={\frac {\det(A_{2})}{\det(A)}}={\frac {\det {\begin{pmatrix}\color {blue}{82}&\color {green}{1}&\color {blue}{9}\\\color {blue}{27}&\color {green}{1}&\color {blue}{3}\\\color {blue}{9}&\color {green}{0}&\color {blue}{1}\end{pmatrix}}}{\det {\begin{pmatrix}\color {blue}{82}&\color {blue}{45}&\color {blue}{9}\\\color {blue}{27}&\color {blue}{16}&\color {blue}{3}\\\color {blue}{9}&\color {blue}{5}&\color {blue}{1}\end{pmatrix}}}}={\frac {1}{1}}=1\qquad$
    $x_{3}={\frac {\det(A_{3})}{\det(A)}}={\frac {\det {\begin{pmatrix}\color {blue}{82}&\color {blue}{45}&\color {green}{1}\\\color {blue}{27}&\color {blue}{16}&\color {green}{1}\\\color {blue}{9}&\color {blue}{5}&\color {green}{0}\end{pmatrix}}}{\det {\begin{pmatrix}\color {blue}{82}&\color {blue}{45}&\color {blue}{9}\\\color {blue}{27}&\color {blue}{16}&\color {blue}{3}\\\color {blue}{9}&\color {blue}{5}&\color {blue}{1}\end{pmatrix}}}}={\frac {-14}{1}}=-14$

## Geschichte

Die Cramersche Regel wurde 1750 von Gabriel Cramer in seinem Buch Introduction à l′analyse des lignes courbes algébriques veröffentlicht. Er gab darin explizit die Formeln für lineare Gleichungssysteme mit bis zu drei Gleichungen an und beschrieb, wie man die Lösungsformeln für Gleichungssysteme mit mehr Gleichungen erstellen kann. Da die Determinante noch nicht eingeführt war, verwendete er Brüche mit je einem Polynom im Zähler und Nenner. Wie der folgende Auszug aus der Originalarbeit zeigt, sind diese mit den Polynomen der Leibniz-Formel identisch.
An diesem Beispiel sieht man auch, dass Cramer noch nicht die heutige Notation linearer Gleichungssysteme verwendete. Mit dieser lautet die Formel wie folgt:
    $x_{1}={\frac {b_{1}a_{22}a_{33}-b_{1}a_{32}a_{23}-b_{2}a_{12}a_{33}+b_{2}a_{32}a_{13}+b_{3}a_{12}a_{23}-b_{3}a_{22}a_{13}}{a_{11}a_{22}a_{33}-a_{11}a_{32}a_{23}-a_{21}a_{12}a_{33}+a_{21}a_{32}a_{13}+a_{31}a_{12}a_{23}-a_{31}a_{22}a_{13}}}$
Cramer selbst war bewusst, dass lineare Gleichungssysteme nicht immer eindeutig lösbar sind. Étienne Bézout zeigte dann 1764, dass der Nenner null wird, wenn das Gleichungssystem nicht eindeutig lösbar ist. Des Weiteren gab Cramer keinen Beweis für seine Formel an. Diesen lieferte erst Augustin Louis Cauchy im Jahr 1815. Dabei führte er auch die heutzutage verwendete Notation der Cramerschen Regel ein.
Gottfried Wilhelm Leibniz brachte die Cramersche Regel schon 1678 in einem Manuskript zu Papier. Dieses wurde allerdings erst später entdeckt und hatte somit keine Auswirkung auf die Entwicklung von Lösungsverfahren für lineare Gleichungssysteme. Colin Maclaurin beschrieb in seinem Werk Treatise of Algebra, das 1748 veröffentlicht wurde, die Spezialfälle der Cramerschen Regel für lineare Gleichungssysteme aus zwei oder drei Gleichungen. Er hatte zwar die Idee, diese Formeln auf Gleichungssysteme mit mehreren Gleichungen zu erweitern, doch im Gegensatz zu Cramer fand er keine Regel, wie man die Vorzeichen in den dabei verwendeten Polynomen richtig setzt. Im 20. Jahrhundert entfachte Carl Benjamin Boyer einen Streit unter Mathematik-Historikern, ob Maclaurin oder Cramer der Entdecker der Formel war. Er empfahl dabei auch eine Umbenennung in Maclaurin-Cramer-Regel.

## Rechenaufwand

Um mit der Cramerschen Regel ein lineares Gleichungssystem mit $n$ Unbekannten zu lösen, müssen $n+1$ Determinanten berechnet werden. Die Anzahl der auszuführenden arithmetischen Operationen hängt damit allein vom Algorithmus zur Berechnung der Determinanten ab.
Werden die Determinanten wie von Cramer mit Hilfe der Leibniz-Formel berechnet, so sind für jede Determinante $(n-1)\cdot n!$ Multiplikationen und $n!-1$ Additionen notwendig. Schon bei einem System aus vier Gleichungen sind 360 Multiplikationen, 115 Additionen und 4 Divisionen notwendig. Im Vergleich zu anderen Verfahren sind dies sehr viele Operationen. Auch unter Verwendung effizienter Algorithmen zur Determinantenberechnung ist der Rechenaufwand für die Lösung eines linearen Gleichungssystems mit der Cramerschen Regel wesentlich höher als beispielsweise beim gaußschen Eliminationsverfahren.
Bei der Berechnung einer $n\times n$-Matrix auf einem Rechner mit $10^{8}$ Gleitkommaoperationen pro Sekunde (100 Mflops) würden sich die folgenden Rechenzeiten ergeben:
| n | 10 | 12 | 14 | 16 | 18 | 20 |
| Rechenzeit | 0,4 s | 1 min | 3,6 h | 41 Tage | 38 Jahre | 16 000 Jahre |

## Beweis

Für diesen Beweis verwendet man eine Matrix $X_{i}$, die entsteht, indem man die $i$-te Spalte der Einheitsmatrix durch den Lösungsvektor $x$ des Gleichungssystems $Ax=b$ ersetzt. So sieht $X_{2}$ für ein Gleichungssystem mit vier Gleichungen folgendermaßen aus:
    $X_{2}={\begin{pmatrix}1&x_{1}&0&0\\0&x_{2}&0&0\\0&x_{3}&1&0\\0&x_{4}&0&1\end{pmatrix}}$
Anhand dieses Beispiels lässt sich auch sehen, dass $AX_{i}=A_{i}$ gilt.
    ${\begin{aligned}A\cdot X_{2}&={\begin{pmatrix}a_{11}&a_{12}&a_{13}&a_{14}\\a_{21}&a_{22}&a_{23}&a_{24}\\a_{31}&a_{32}&a_{33}&a_{34}\\a_{41}&a_{42}&a_{43}&a_{44}\end{pmatrix}}\cdot {\begin{pmatrix}1&x_{1}&0&0\\0&x_{2}&0&0\\0&x_{3}&1&0\\0&x_{4}&0&1\end{pmatrix}}\\&={\begin{pmatrix}a_{11}&a_{11}x_{1}+a_{12}x_{2}+a_{13}x_{3}+a_{14}x_{4}&a_{13}&a_{14}\\a_{21}&a_{21}x_{1}+a_{22}x_{2}+a_{23}x_{3}+a_{24}x_{4}&a_{23}&a_{24}\\a_{31}&a_{31}x_{1}+a_{32}x_{2}+a_{33}x_{3}+a_{34}x_{4}&a_{33}&a_{34}\\a_{41}&a_{41}x_{1}+a_{42}x_{2}+a_{43}x_{3}+a_{44}x_{4}&a_{43}&a_{44}\end{pmatrix}}\\&={\begin{pmatrix}a_{11}&b_{1}&a_{13}&a_{14}\\a_{21}&b_{2}&a_{23}&a_{24}\\a_{31}&b_{3}&a_{33}&a_{34}\\a_{41}&b_{4}&a_{43}&a_{44}\end{pmatrix}}\\&=A_{2}\end{aligned}}$
Da zusätzlich $\det(X_{i})=x_{i}$ gilt, folgt mit der Produktregel für Determinanten
    ${\begin{matrix}\det(A)&\cdot &\det(X_{i})&=&\det(A_{i})\\\det(A)&\cdot &x_{i}&=&\det(A_{i})\\&&x_{i}&=&\det(A_{i})\cdot \det(A)^{-1}\end{matrix}}$
Da die Determinante $\det(A)$ nach Voraussetzung nicht das Null-Element ist, existiert $\det(A)^{-1}$.

## Verallgemeinerung

Eine Verallgemeinerung – und gleichzeitig einen Schritt im Beweis – der Cramerschen Regel stellt der folgende Satz dar: Gegeben sei ein quadratisches lineares Gleichungssystem der Form
    $Ax=b$.
Ist $x=(x_{1},x_{2},\ldots ,x_{n})$ eine Lösung dieses linearen Gleichungssystems, dann gilt
    $\det(A)\cdot x_{i}=\det(A_{i})$ für jedes $i$.
Die Matrix $A_{i}$ entsteht aus der Koeffizientenmatrix $A$, indem die $i$-te Spalte von $A$ durch die rechte Seite des Gleichungssystems $b$ ersetzt wird.
Es entfällt die Einschränkung auf ein eindeutig lösbares Gleichungssystem. Da zudem keine Division mehr auftritt, gilt der Satz für alle Gleichungssysteme mit Koeffizienten aus einem kommutativen Ring. Diese Verallgemeinerung wird nicht mehr Cramersche Regel genannt.

## Folgerungen aus der Cramerschen Regel
### Die Inverse einer Matrix

Die einzelnen Spalten der Inversen einer Matrix $A$ werden jeweils von der Lösung des Gleichungssystems $Ax=e_{j}$ mit dem $j$-ten Einheitsvektor auf der rechten Seite gebildet. Berechnet man diese mit der Cramerschen Regel, so erhält man unter Verwendung der Adjunkten $\operatorname {adj} (A)$ bzw. der Kofaktormatrix $\operatorname {cof} (A)=\operatorname {adj} (A)^{T}$ die Formel
    $A^{-1}={\frac {1}{\det(A)}}\operatorname {adj} (A)={\frac {1}{\det(A)}}\operatorname {cof} (A)^{T}$
Diese Formel zeigt auch eine Eigenschaft von Matrizen mit Einträgen aus einem kommutativen Ring anstatt eines Körpers. Diese sind demnach genau dann invertierbar, wenn $\det(A)$ invertierbar ist.

### Lösung eines homogenen linearen Gleichungssystems

Mit Hilfe der Cramerschen Regel lässt sich einfach zeigen, dass die triviale Lösung $x_{1}=x_{2}=\dotsb =x_{n}=0$ die einzige Lösung eines jeden homogenen linearen Gleichungssystems mit $\det(A)\neq 0$ ist. Da bei allen Matrizen $A_{i}$ in der $i$-ten Spalte nur Nullen stehen, sind deren Spalten nicht mehr linear unabhängig, und es gilt deshalb $\det(A_{i})=0$. Damit berechnen sich alle $x_{i}$ zu null.
Aus der obigen Eigenschaft folgt direkt, dass der Kern eines linearen Gleichungssystems $Ax=b$ mit $\det(A)\neq 0$ der Nullvektor ist und dieses deshalb eindeutig lösbar ist.