Lineare Algebra

Die lineare Algebra (auch Vektoralgebra) ist ein Teilgebiet der Mathematik, das sich mit Vektorräumen beschäftigt. Ähnlich wie in anderen Teilgebieten der Mathematik, sind die strukturerhaltenden Abbildungen, welche in der linearen Algebra die linearen Abbildungen sind, von besonderem Interesse. Diese können durch Matrizen repräsentiert werden. Die lineare Algebra schließt insbesondere auch die Betrachtung von linearen Gleichungssystemen mit ein.

Vektorräume und deren lineare Abbildungen sind ein wichtiges Hilfsmittel in vielen Bereichen der Mathematik. Außerhalb der reinen Mathematik finden sich Anwendungen unter anderem in den Naturwissenschaften, in der Informatik und in der Wirtschaftswissenschaft (zum Beispiel in der Optimierung).

Die lineare Algebra entstand aus zwei konkreten Anforderungen heraus: einerseits dem Lösen von linearen Gleichungssystemen, andererseits der rechnerischen Beschreibung geometrischer Objekte, der sogenannten analytischen Geometrie (daher bezeichnen manche Autoren lineare Algebra als lineare Geometrie).

Die Anfänge der Algebra und somit auch der Begriff selbst gehen weitestgehend auf den persisch-choresmischen Mathematiker, Astronomen, Geographen und Universalgelehrten Al-Chwarizmi zurück, der aufgrund der Islamisierung im Iran seine Werke ins Arabische übersetzen musste und so auf den Namen „al-jabr“ kam. Daraus leitet sich der Begriff der Algebra her.^[1]

Während die Entwicklung der Algebra bereits im alten Ägypten begann, entstand die lineare Algebra als eigenständiges Teilgebiet erst im 17. Jahrhundert mit der Theorie der Determinante. Die Entwicklung dieser Theorie wurde unabhängig voneinander von Gottfried Wilhelm Leibniz und Seki Takakazu gestartet. Im Jahr 1750 veröffentlichte dann Gabriel Cramer die nach ihm benannte cramersche Regel. Damit war man erstmals im Besitz einer Lösungsformel für viele lineare Gleichungssysteme.^[2]

Die Geschichte der modernen linearen Algebra reicht zurück bis in die Jahre 1843 und 1844. 1843 erdachte William Rowan Hamilton (von dem der Begriff Vektor stammt) mit den Quaternionen eine Erweiterung der komplexen Zahlen. 1844 veröffentlichte Hermann Graßmann sein Buch Die lineale Ausdehnungslehre. Arthur Cayley führte dann 1857 mit den ( 2 × 2 ) {\displaystyle (2\times 2)} -Matrizen eine der grundlegendsten algebraischen Ideen ein.

Ab dem 20. Jahrhundert befasste man sich dann mehrheitlich mit dem Begriff des Vektorraums. Insbesondere die Mathematiker August Ferdinand Möbius, Constantin Carathéodory und Hermann Weyl leisteten hierfür die Vorarbeit. So wurde beispielsweise festgestellt, dass lineare Abbildungen zwischen endlichdimensionalen Vektorräumen durch Matrizen beschrieben werden können. Auf dieser Erkenntnis basierend konnte Stefan Banach als Erster eine axiomatische Definition für reelle Vektorräume angeben.

Als lineares Gleichungssystem bezeichnet man eine Zusammenfassung von Gleichungen der Art

x 1 + x 2 = 1 {\displaystyle x_{1}+x_{2}=1}

3 x 1 + 6 x 2 = 4 {\displaystyle 3x_{1}+6x_{2}=4}

Derartige Gleichungssysteme erhält man aus vielen alltäglichen Fragestellungen, beispielsweise:

In welchem Verhältnis muss man eine 30%ige Lösung (entspricht x 1 {\displaystyle x_{1}} ) und eine 60%ige Lösung (entspricht x 2 {\displaystyle x_{2}} ) mischen, um eine 40%ige Lösung zu erhalten?

Der wesentliche Abstraktionsschritt der linearen Algebra besteht nun darin, die linken Seiten als eine Funktion A {\displaystyle A} der Unbekannten x = ( x 1 , x 2 ) {\displaystyle x=(x_{1},x_{2})} (in diesem Fall die Menge der jeweiligen Lösungen) aufzufassen:

A ( x ) = ( x 1 + x 2 3 x 1 + 6 x 2 ) {\displaystyle A(x)={\begin{pmatrix}x_{1}+x_{2}\\3x_{1}+6x_{2}\end{pmatrix}}}

Dann wird die Lösung des Gleichungssystems zu der Aufgabe: Finde ein x {\displaystyle x} , sodass

A ( x ) = ( 1 4 ) {\displaystyle A(x)={\begin{pmatrix}1\\4\end{pmatrix}}}

gilt. Das Übereinanderschreiben ist dabei lediglich ein Formalismus, um mit mehr als einer Zahl gleichzeitig umgehen zu können.

Statt A {\displaystyle A} schreibt man auch einfach die relevanten Zahlen in Form eines Rechtecks auf und nennt das Objekt eine Matrix:

A = ( 1 1 3 6 ) {\displaystyle A={\begin{pmatrix}1&1\\3&6\end{pmatrix}}}

Man stellt fest, dass die Funktion A {\displaystyle A} spezielle Eigenschaften hat, sie ist eine lineare Abbildung. Ist x {\displaystyle x} eine Lösung für das Gleichungssystem A ( x ) = b {\displaystyle A(x)=b} , und y {\displaystyle y} eine Lösung des Gleichungssystems A ( y ) = c {\displaystyle A(y)=c} , so ist

z = x + y = ( x 1 + y 1 x 2 + y 2 ) {\displaystyle z=x+y={\begin{pmatrix}x_{1}+y_{1}\\x_{2}+y_{2}\end{pmatrix}}}

eine Lösung von A ( z ) = b + c {\displaystyle A(z)=b+c} . Man kann das auch in der Form A ( x + y ) = A ( x ) + A ( y ) {\displaystyle A(x+y)=A(x)+A(y)} schreiben. Ist weiter λ {\displaystyle \lambda } irgendeine reelle Zahl, so ist A ( λ x ) = λ ⋅ A ( x ) {\displaystyle A(\lambda x)=\lambda \cdot A(x)} ; dabei ist

λ x = ( λ x 1 λ x 2 ) {\displaystyle \lambda x={\begin{pmatrix}\lambda x_{1}\\\lambda x_{2}\end{pmatrix}}} .

Der andere Ursprung der linearen Algebra findet sich in der rechnerischen Beschreibung des 2- und 3-dimensionalen (euklidischen) Raumes, auch „Anschauungsraum“ genannt. Mit Hilfe eines Koordinatensystems können Punkte im Raum durch Tripel ( x 1 , x 2 , x 3 ) {\displaystyle (x_{1},x_{2},x_{3})} von Zahlen beschrieben werden. Der Abbildungstyp der Verschiebung führt zum Begriff des Vektors, der Richtung und Betrag der Verschiebung angibt. Viele physikalische Größen, beispielsweise Kräfte, haben stets diesen Richtungsaspekt.

Da man auch Vektoren durch Zahlentripel ( a 1 , a 2 , a 3 ) {\displaystyle (a_{1},a_{2},a_{3})} beschreiben kann, verschwimmt die Trennung zwischen Vektoren und Punkten: Einem Punkt P {\displaystyle P} entspricht sein Ortsvektor, der vom Koordinatenursprung nach P {\displaystyle P} zeigt.

Viele der in der klassischen Geometrie betrachteten Abbildungstypen, beispielsweise Drehungen um Achsen durch den Ursprung oder Spiegelungen an Ebenen durch den Ursprung, gehören zur Klasse der linearen Abbildungen, die schon oben erwähnt wurde.

Der Begriff des Vektorraumes entsteht als Abstraktion der obigen Beispiele: Ein Vektorraum ist eine Menge, deren Elemente Vektoren genannt werden, zusammen mit

• einer Addition von Vektoren
• einer Multiplikation von Vektoren mit Elementen eines fixierten Körpers, Skalarmultiplikation (äußere Multiplikation) genannt.

Diese Addition und die Skalarmultiplikation müssen noch einige einfache Eigenschaften erfüllen, die auch für die Vektoren im Anschauungsraum gelten.

Man könnte sagen, dass Vektorräume gerade so definiert sind, dass man von linearen Abbildungen zwischen ihnen sprechen kann.

In gewisser Weise ist der Begriff des Vektorraums für die lineare Algebra bereits zu allgemein. Jedem Vektorraum ist eine Dimension zugeordnet, beispielsweise hat die Ebene Dimension 2 {\displaystyle 2} und der Anschauungsraum die Dimension 3 {\displaystyle 3} . Es gibt aber Vektorräume, deren Dimension nicht endlich ist, wodurch viele der bekannten Eigenschaften verloren gehen. Es hat sich aber als sehr erfolgreich erwiesen, unendlichdimensionale Vektorräume mit einer zusätzlichen topologischen Struktur auszustatten; die Untersuchung topologischer Vektorräume ist Gegenstand der Funktionalanalysis.

(Der Rest dieses Artikels beschäftigt sich mit dem Fall endlicher Dimensionen.)

Jeder Vektorraum hat (unter der Annahme, dass das Auswahlaxiom gilt,) mindestens eine Basis. Die Basis kann endlich oder unendlich viele Elemente enthalten.

Falls ein Vektorraum eine Basis aus endlich vielen Elementen hat, haben alle Basen dieses Vektorraumes endlich viele Elemente und die Anzahl der Elemente ist für alle Basen gleich. Falls ein Vektorraum eine Basis aus unendlich vielen Elementen hat, haben alle Basen dieses Vektorraumes unendlich viele Elemente. Deshalb ist es sinnvoll, von der Dimension eines Vektorraumes als die Anzahl der Elemente einer Basis sowie von endlich- und unendlich-dimensionalen Vektorräumen zu sprechen. Für Summen und Durchschnitte von Untervektorräumen gilt die Dimensionsformel und für die Dimensionen von Faktorräumen eines endlich-dimensionalen Vektorraumes V {\displaystyle V} die Formel dim ⁡ V / U = dim ⁡ V − dim ⁡ U {\displaystyle \dim V/U=\dim V-\dim U} .

Jede lineare Abbildung f : V → W {\displaystyle f\colon V\to W} ist durch die Angabe der Bilder einer Basis von V {\displaystyle V} eindeutig festgelegt. Für lineare Abbildungen gelten der Homomorphiesatz und der Rangsatz. Lineare Abbildungen können bezüglich fest gewählter Basen durch Matrizen dargestellt werden. Dabei entspricht der Hintereinanderausführung von linearen Abbildungen die Multiplikation ihrer Darstellungsmatrizen.

Ein lineares Gleichungssystem A ⋅ x = b {\displaystyle A\cdot x=b} mit A ∈ K m × n {\displaystyle A\in \mathbb {K} ^{{m}\times {n}}} , b ∈ K m {\displaystyle b\in \mathbb {K} ^{m}} und x ∈ K n {\displaystyle x\in \mathbb {K} ^{n}} ist genau dann lösbar, wenn der Rang der Matrix A {\displaystyle A} gleich dem Rang der erweiterten Koeffizientenmatrix ( A b ) {\displaystyle {\begin{pmatrix}A&b\end{pmatrix}}} ist. In diesem Fall ist die Lösungsmenge des Systems ein affiner Unterraum von K n {\displaystyle \mathbb {K} ^{n}} der Dimension n − r a n g ( A ) {\displaystyle n-\mathrm {rang} (A)} . Für nicht zu große Gleichungssysteme können die Rangbestimmung und die Berechnung des Lösungsraumes mit dem Gaußschen Eliminationsverfahren durchgeführt werden.

Eine lineare Abbildung f : V → V {\displaystyle f\colon V\to V} (also ein Endomorphismus) eines endlichdimensionalen Vektorraumes V {\displaystyle V} ist bereits invertierbar, wenn sie injektiv oder surjektiv ist. Dies ist wiederum genau dann der Fall, wenn ihre Determinante ungleich null ist. Hieraus folgt, dass die Eigenwerte eines Endomorphismus genau die Nullstellen seines charakteristischen Polynoms sind. Eine weitere wichtige Aussage über das charakteristische Polynom ist der Satz von Cayley-Hamilton.

Ein Endomorphismus (beziehungsweise eine quadratische Matrix) ist genau dann diagonalisierbar, wenn das charakteristische Polynom in Linearfaktoren zerfällt und für jeden Eigenwert dessen algebraische Vielfachheit gleich der geometrischen Vielfachheit, also die Nullstellenordnung des Eigenwerts im charakteristischen Polynom gleich der Dimension des zugehörigen Eigenraumes ist. Äquivalent dazu ist die Existenz einer Basis des Vektorraumes, die aus Eigenvektoren der linearen Abbildung besteht. Endomorphismen, deren charakteristisches Polynom in Linearfaktoren zerfällt, sind immerhin noch trigonalisierbar, können also durch eine Dreiecksmatrix dargestellt werden. Ein etwas tiefer liegendes Ergebnis ist, dass die darstellende Matrix dabei sogar in jordansche Normalform gebracht werden kann.

In Vektorräumen, auf denen zusätzlich ein Skalarprodukt ⟨ ⋅ , ⋅ ⟩ {\displaystyle \langle \cdot ,\cdot \rangle } gegeben ist, wird durch ‖ x ‖ := ⟨ x , x ⟩ {\displaystyle \|x\|:={\sqrt {\langle x,x\rangle }}} eine Norm definiert. In diesen Skalarprodukträumen existieren stets Orthonormalbasen, die etwa durch das Gram-Schmidtsche Orthonormalisierungsverfahren konstruiert werden können. Nach dem Projektionssatz kann man in diesen Räumen die Bestapproximation aus einem Untervektorraum durch orthogonale Projektion bestimmen.

Bezüglich der Diagonalisierbarkeit von Endomorphismen in Skalarprodukträumen stellt sich die Frage, ob eine Orthonormalbasis aus Eigenvektoren existiert. Das zentrale Resultat hierzu ist der Spektralsatz. Insbesondere gilt im reellen Fall: Zu jeder symmetrischen Matrix A ∈ R n × n {\displaystyle A\in \mathbb {\mathbb {R} } ^{{n}\times {n}}} gibt es eine orthogonale Matrix Q {\displaystyle Q} , sodass Q T A Q {\displaystyle Q^{T}AQ} eine Diagonalmatrix ist. Wendet man dieses Ergebnis auf quadratische Formen an, ergibt sich der Satz von der Hauptachsentransformation.

Auch Bilinearformen und Sesquilinearformen können bei fest gewählten Basen durch Matrizen dargestellt werden. Eine Bilinearform ist genau dann symmetrisch und positiv definit, also ein Skalarprodukt, wenn ihre darstellende Matrix symmetrisch und positiv definit ist. Eine symmetrische Matrix ist genau dann positiv definit, wenn alle ihre Eigenwerte positiv sind. Allgemein gilt für symmetrische Bilinearformen und hermitesche Sesquilinearformen der Trägheitssatz von Sylvester, der besagt, dass die Anzahl der positiven und negativen Eigenwerte der darstellenden Matrizen nicht von der Wahl der Basis abhängen.

Vektoren endlichdimensionaler Räume können durch ihre Komponenten beschrieben werden, die (je nach Anwendung) als Spaltenvektor

a = ( 3 7 2 ) {\displaystyle \mathbf {a} ={\begin{pmatrix}3\\7\\2\end{pmatrix}}}

oder Zeilenvektor

b = ( 4 6 3 7 ) {\displaystyle \mathbf {b} ={\begin{pmatrix}4&6&3&7\end{pmatrix}}}

geschrieben werden. Häufig werden Zeilenvektoren mit einem hochgestellten T für transponiert, wie b T {\displaystyle b^{T}} , gekennzeichnet.

In der Literatur werden Vektoren auf unterschiedliche Weise von anderen Größen unterschieden: Es werden Kleinbuchstaben, fettgedruckte Kleinbuchstaben, unterstrichene Kleinbuchstaben, Kleinbuchstaben mit einem Pfeil darüber oder kleine Frakturbuchstaben benutzt. Dieser Artikel verwendet Kleinbuchstaben.

Eine Matrix wird durch ein „Raster“ von Zahlen angegeben. Hier ist eine Matrix mit vier Zeilen und drei Spalten:

M = ( 8 2 9 4 8 2 8 3 7 5 9 1 ) {\displaystyle \mathbf {M} ={\begin{pmatrix}8&2&9\\4&8&2\\8&3&7\\5&9&1\end{pmatrix}}}

Matrizen werden meistens mit Großbuchstaben bezeichnet.

Einzelne Elemente eines Vektors werden bei Spaltenvektoren in der Regel durch einen Index angegeben: Das zweite Element des oben angegebenen Vektors a {\displaystyle a} wäre dann a 2 = 7 {\displaystyle a_{2}=7} . In Zeilenvektoren wird manchmal eine Hochzahl verwendet, wobei man aufpassen muss, ob eine Vektorindizierung oder ein Exponent vorliegt: Mit dem obigen Beispiel b {\displaystyle b} hat man etwa b 4 = 7 {\displaystyle b^{4}=7} . Matrixelemente werden durch zwei Indizes angegeben. Dabei werden die Elemente durch Kleinbuchstaben dargestellt: m 2 , 3 = 2 {\displaystyle m_{2,3}=2} ist das Element in der zweiten Zeile der dritten Spalte (statt „in der dritten Spalte der zweiten Zeile“, denn so lässt sich m 2 , 3 {\displaystyle m_{2,3}} leichter lesen).

Der verallgemeinerte Begriff dieser Gebilde ist Tensor, Skalare sind Tensoren nullter Stufe, Vektoren Tensoren erster Stufe, Matrizen Tensoren zweiter Stufe. Ein Tensor n {\displaystyle n} -ter Stufe kann durch einen n {\displaystyle n} -dimensionalen Zahlenwürfel repräsentiert werden.

Oftmals ist es erforderlich, Matrizen mittels elementarer Zeilenumformungen oder Basiswechsel auf eine spezielle Form zu bringen. Wichtig sind dabei insbesondere die Dreiecksform, die Diagonalform und die jordansche Normalform.

Bei der Darstellung einer linearen Abbildung – wie unter Matrix beschrieben – gibt es den Sonderfall einer linearen Abbildung f {\displaystyle f} eines endlichdimensionalen Vektorraums auf sich selbst (eines sog. Endomorphismus). Man kann dann dieselbe Basis v {\displaystyle v} für Urbild- und Bildkoordinaten verwenden und erhält eine quadratische Matrix A {\displaystyle A} , sodass die Anwendung der linearen Abbildung der Linksmultiplikation mit A {\displaystyle A} entspricht. Um die Abhängigkeit von f {\displaystyle f} und v {\displaystyle v} zum Ausdruck zu bringen, verwendet man Schreibweisen wie A = M v ( f ) {\displaystyle A=M_{v}(f)} oder A = v f v {\displaystyle A={}_{v}f_{v}} . Die zweimalige Hintereinanderausführung dieser Abbildung entspricht dann der Multiplikation mit A 2 {\displaystyle A^{2}} usw., und man kann alle polynomialen Ausdrücke mit A {\displaystyle A} (Summen von Vielfachen von Potenzen von A {\displaystyle A} ) als lineare Abbildungen des Vektorraums auffassen.

Zu einer invertierbaren Matrix A {\displaystyle A} existiert eine inverse Matrix A − 1 {\displaystyle A^{-1}} mit A − 1 A = A A − 1 = E {\displaystyle A^{-1}A=AA^{-1}=E} . Analog zur Rechenregel x 0 = 1 {\displaystyle x^{0}=1} bei Zahlen ist die nullte Potenz einer quadratischen Matrix die Diagonalmatrix E {\displaystyle E} (Einheitsmatrix) mit Einsen auf der Diagonalen und in der alle restlichen Elemente Null sind, sie entspricht der Identitätsabbildung jedes Vektors auf sich selbst. Negative Potenzen einer quadratischen Matrix A {\displaystyle A} lassen sich nur berechnen, wenn die durch A {\displaystyle A} gegebene lineare Abbildung invertierbar ist, also keine zwei unterschiedlichen Vektoren u 1 {\displaystyle u_{1}} und u 2 {\displaystyle u_{2}} auf denselben Vektor A u 1 = A u 2 {\displaystyle Au_{1}=Au_{2}} abbildet. Anders ausgedrückt, muss für eine invertierbare Matrix A {\displaystyle A} aus u 1 − u 2 ≠ 0 {\displaystyle u_{1}-u_{2}\neq 0} stets A ( u 1 − u 2 ) ≠ 0 {\displaystyle A(u_{1}-u_{2})\neq 0} folgen, das lineare Gleichungssystem A u = 0 {\displaystyle Au=0} darf also nur die Lösung 0 {\displaystyle 0} haben.

Eine Determinante ist eine spezielle Funktion, die einer quadratischen Matrix eine Zahl zuordnet. Diese Zahl gibt Auskunft über einige Eigenschaften der Matrix. Beispielsweise lässt sich an ihr erkennen, ob eine Matrix invertierbar ist. Eine weitere wichtige Anwendung ist die Berechnung des charakteristischen Polynoms und damit der Eigenwerte der Matrix.

Es gibt geschlossene Formeln zur Berechnung der Determinanten, wie den Laplace’schen Entwicklungssatz oder die Leibniz-Formel. Diese Formeln sind jedoch eher von theoretischem Wert, da ihr Aufwand bei größeren Matrizen stark ansteigt. In der Praxis kann man Determinanten am leichtesten berechnen, indem man die Matrix mit Hilfe des Gauß-Algorithmus in obere oder untere Dreiecksform bringt, die Determinante ist dann einfach das Produkt der Hauptdiagonalelemente.

Obige Begriffe sollen an einem durch die Fibonacci-Folge motivierten Beispiel verdeutlicht werden.

Die Fibonacci-Folge f n {\displaystyle f_{n}} ist rekursiv durch die Gleichungen f 0 = 0 {\displaystyle f_{0}=0} , f 1 = 1 {\displaystyle f_{1}=1} und f n + 1 = f n + f n − 1 {\displaystyle f_{n+1}=f_{n}+f_{n-1}} für n ≥ 1 {\displaystyle n\geq 1} definiert, was gleichbedeutend ist mit

( f 1 f 0 ) = ( 1 0 ) {\displaystyle {\begin{pmatrix}f_{1}\\f_{0}\end{pmatrix}}={\begin{pmatrix}1\\0\end{pmatrix}}}

und

( f n + 1 f n ) = ( 1 1 1 0 ) ⋅ ( f n f n − 1 ) für n ≥ 1 {\displaystyle {\begin{pmatrix}f_{n+1}\\f_{n}\end{pmatrix}}={\begin{pmatrix}1&1\\1&0\end{pmatrix}}\cdot {\begin{pmatrix}f_{n}\\f_{n-1}\end{pmatrix}}\quad {\text{für}}\quad n\geq 1} ,

woraus durch Iteration die nichtrekursive Formel

( f n + 1 f n ) = ( 1 1 1 0 ) n ⋅ ( 1 0 ) für n ≥ 0 {\displaystyle {\begin{pmatrix}f_{n+1}\\f_{n}\end{pmatrix}}={\begin{pmatrix}1&1\\1&0\end{pmatrix}}^{n}\cdot {\begin{pmatrix}1\\0\end{pmatrix}}\quad {\text{für}}\quad n\geq 0}

folgt, in der die n {\displaystyle n} -te Potenz einer Matrix A {\displaystyle A} vorkommt.

Das Verhalten einer solchen Matrix bei Potenzierung ist nicht leicht zu erkennen; hingegen wird die n {\displaystyle n} -te Potenz einer Diagonalmatrix einfach durch Potenzierung jedes einzelnen Diagonaleintrags berechnet. Wenn es nun eine invertierbare Matrix T {\displaystyle T} gibt, sodass T − 1 A T {\displaystyle T^{-1}AT} Diagonalform hat, lässt sich die Potenzierung von A {\displaystyle A} auf die Potenzierung einer Diagonalmatrix zurückführen gemäß der Gleichung ( T − 1 A T ) n = T − 1 A n T {\displaystyle (T^{-1}AT)^{n}=T^{-1}A^{n}T} (die linke Seite dieser Gleichung ist dann die n {\displaystyle n} -te Potenz einer Diagonalmatrix). Allgemein lässt sich durch Diagonalisierung einer Matrix ihr Verhalten (bei Potenzierung, aber auch bei anderen Operationen) leichter erkennen.

Fasst man A = v f v {\displaystyle A={}_{v}f_{v}} als Matrix einer linearen Abbildung auf, so ist die Transformationsmatrix T {\displaystyle T} die Basiswechselmatrix zu einer anderen Basis v ′ {\displaystyle v'} , also T = v e v ′ {\displaystyle T={}_{v}e_{v'}} (wobei die Identitätsabbildung e {\displaystyle e} jeden Vektor auf sich selbst abbildet). Dann ist nämlich T − 1 A T = v ′ f v ′ {\displaystyle T^{-1}AT={}_{v'}f_{v'}} .

Im oben genannten Beispiel lässt sich eine Transformationsmatrix T {\displaystyle T} finden, sodass

T − 1 ⋅ A ⋅ T = ( ϕ 0 0 1 − ϕ ) {\displaystyle T^{-1}\cdot A\cdot T={\begin{pmatrix}\phi &0\\0&1-\phi \end{pmatrix}}}

eine Diagonalmatrix ist, in der der goldene Schnitt ϕ = 1 + 5 2 {\displaystyle \phi ={\frac {1+{\sqrt {5}}}{2}}} vorkommt. Hieraus erhält man schließlich die Formel von Binet:

f n = 1 5 ⋅ [ ( 1 + 5 2 ) n − ( 1 − 5 2 ) n ] {\displaystyle f_{n}={\frac {1}{\sqrt {5}}}\cdot \left[\left({\frac {1+{\sqrt {5}}}{2}}\right)^{n}-\left({\frac {1-{\sqrt {5}}}{2}}\right)^{n}\right]}

Wie kommt man von der Matrix A {\displaystyle A} auf die Zahl ϕ {\displaystyle \phi } ? An der Diagonalmatrix erkennt man sofort

( ϕ 0 0 1 − ϕ ) ⋅ ( 1 0 ) = ( ϕ 0 ) {\displaystyle {\begin{pmatrix}\phi &0\\0&1-\phi \end{pmatrix}}\cdot {\begin{pmatrix}1\\0\end{pmatrix}}={\begin{pmatrix}\phi \\0\end{pmatrix}}} ,

dass es also einen Vektor u {\displaystyle u} ungleich Null gibt, der durch Multiplikation mit der Diagonalmatrix komponentenweise vervielfacht (genauer: ver- ϕ {\displaystyle \phi } -facht) wird: ( T − 1 A T ) u = ϕ u {\displaystyle (T^{-1}AT)u=\phi u} . Die Zahl ϕ {\displaystyle \phi } heißt wegen dieser Eigenschaft ein Eigenwert der Matrix T − 1 A T {\displaystyle T^{-1}AT} (mit Eigenvektor u {\displaystyle u} ). Im Fall von Diagonalmatrizen sind die Eigenwerte gleich den Diagonaleinträgen.

ϕ {\displaystyle \phi } ist aber auch zugleich Eigenwert der ursprünglichen Matrix A {\displaystyle A} (mit Eigenvektor T u {\displaystyle Tu} , denn A ( T u ) = ϕ ( T u ) {\displaystyle A(Tu)=\phi (Tu)} ), die Eigenwerte bleiben bei Transformation der Matrix also unverändert. Die Diagonalform der Matrix A {\displaystyle A} ergibt sich demnach aus deren Eigenwerten, und um die Eigenwerte von A {\displaystyle A} zu finden, muss man untersuchen, für welche Zahlen x {\displaystyle x} das lineare Gleichungssystem A u = x u {\displaystyle Au=xu} eine von Null verschiedene Lösung u {\displaystyle u} hat (oder, anders ausgedrückt, die Matrix x E − A {\displaystyle xE-A} nicht invertierbar ist).

Die gesuchten Zahlen x {\displaystyle x} sind genau diejenigen, die die Determinante der Matrix x E − A {\displaystyle xE-A} zu Null machen. Diese Determinante ist ein polynomialer Ausdruck in x {\displaystyle x} (das sogenannte charakteristische Polynom von A {\displaystyle A} ); im Falle der oben genannten 2×2-Matrix A {\displaystyle A} ergibt dies die quadratische Gleichung x 2 − x − 1 = 0 {\displaystyle x^{2}-x-1=0} mit den beiden Lösungen x = ϕ {\displaystyle x=\phi } und x = 1 − ϕ {\displaystyle x=1-\phi } . Die zugehörigen Eigenvektoren sind Lösungen der linearen Gleichungssysteme A u = ϕ u {\displaystyle Au=\phi u} beziehungsweise A u = ( 1 − ϕ ) u {\displaystyle Au=(1-\phi )u} , sie bilden dann die Spalten der Transformationsmatrix T {\displaystyle T} .

Ob eine Matrix diagonalisierbar ist, hängt vom verwendeten Zahlbereich ab. A {\displaystyle A} ist zum Beispiel über den rationalen Zahlen nicht diagonalisierbar, weil die Eigenwerte ϕ {\displaystyle \phi } und 1 − ϕ {\displaystyle 1-\phi } irrationale Zahlen sind. Die Diagonalisierbarkeit kann aber auch unabhängig vom Zahlbereich scheitern, wenn nicht „genügend“ Eigenwerte vorhanden sind; so hat etwa die Jordanform-Matrix

( 1 1 0 1 ) {\displaystyle {\begin{pmatrix}1&1\\0&1\end{pmatrix}}}

nur den Eigenwert 1 {\displaystyle 1} (als Lösung der quadratischen Gleichung ( x − 1 ) 2 = 0 {\displaystyle (x-1)^{2}=0} ) und ist nicht diagonalisierbar. Bei genügend großem Zahlbereich (zum Beispiel über den komplexen Zahlen) lässt sich aber jede Matrix diagonalisieren oder in Jordansche Normalform transformieren.

Da die Transformation einer Matrix dem Basiswechsel einer linearen Abbildung entspricht, besagt diese letzte Aussage, dass man zu einer linearen Abbildung bei genügend großem Zahlbereich stets eine Basis wählen kann, die „auf einfache Weise“ abgebildet wird: Im Fall der Diagonalisierbarkeit wird jeder Basisvektor auf ein Vielfaches von sich abgebildet (ist also ein Eigenvektor); im Fall der Jordanform auf ein Vielfaches von sich plus evtl. den vorigen Basisvektor. Diese Theorie der linearen Abbildung lässt sich auf Körper verallgemeinern, die nicht „genügend groß“ sind; in ihnen müssen neben der Jordanform andere Normalformen betrachtet werden (zum Beispiel die Frobenius-Normalform).

• Howard Anton: Lineare Algebra. Spektrum Akademischer Verlag, Heidelberg, ISBN 978-3-8274-0324-7.
• Albrecht Beutelspacher: Lineare Algebra. Vieweg-Verlag, ISBN 978-3-658-02412-3.
• Siegfried Bosch: Lineare Algebra. Springer-Lehrbuch, ISBN 978-3-540-76437-3.
• Egbert Brieskorn: Lineare Algebra und analytische Geometrie. Band 1, Vieweg-Verlag, 2012, ISBN 978-3-322-83175-0.
• Egbert Brieskorn: Lineare Algebra und analytische Geometrie. Band 2, Vieweg-Verlag, 1985, ISBN 978-3-528-08562-9.
• Theodor Bröcker: Lineare Algebra und Analytische Geometrie. Birkhäuser Verlag, ISBN 978-3-7643-7144-9.
• Gerd Fischer: Lineare Algebra. Vieweg-Verlag, ISBN 978-3-658-03944-8.
• Günter Gramlich: Lineare Algebra. Carl Hanser Verlag, ISBN 978-3-446-44140-8.
• Günter Gramlich: Anwendungen der Linearen Algebra. Carl Hanser Verlag, ISBN 978-3-446-22655-5.
• Klaus Jänich: Lineare Algebra. Springer-Lehrbuch, ISBN 978-3-540-75501-2.
• Hans-Joachim Kowalsky: Lineare Algebra. de Gruyter Lehrbuch, ISBN 978-3-11-017963-7.
• Burkhard Lenze: Basiswissen Lineare Algebra. Springer-Vieweg, 2020, ISBN 978-3-658-29968-2.
• Jörg Liesen, Volker Mehrmann: Lineare Algebra. 3. Auflage. Springer, Berlin, Heidelberg 2021, ISBN 978-3-662-62741-9, doi:10.1007/978-3-662-62742-6. 
• Falko Lorenz: Lineare Algebra. 2 Bände, BI/Spektrum, 2003, ISBN 3-8274-1406-7.
• Gilbert Strang: Lineare Algebra. Springer-Lehrbuch, ISBN 978-0-9802327-7-6.
• Shafarevich Igor, Remizov Alexey: Linear Algebra and Geometry. Springer, 2012, ISBN 978-3-642-30993-9.

• MIT OpenCourseWare: Video-Lektionen zur linearen Algebra gehalten von Professor Gilbert Strang; aufgenommen 1999 (englisch).
• 17 Kapitel lineare Algebra bei: mathproject.de in Deutsch und Englisch.
• Mathematik-Online-Kurs – kurze Lektionen zu vielen Themen der linearen Algebra.
• Howard Anton, Lineare Algebra, Buch in deutscher Sprache: Download

1. ↑ John Stillwell: Mathematics and Its History. Springer, New York, NY 2010, ISBN 978-1-4419-6052-8, S. 88–89, doi:10.1007/978-1-4419-6053-5_6. 
2. ↑ Heinz-Wilhelm Alten: 4000 Jahre Algebra. Geschichte, Kulturen, Menschen. Springer, Berlin u. a. 2003, ISBN 3-540-43554-9, S. 335–339.