Jordansche Normalform

Die jordansche Normalform ist ein Begriff aus dem mathematischen Teilgebiet der linearen Algebra. Benannt wurde sie nach Marie Ennemond Camille Jordan, der sie 1870 für endliche Körper und 1871 im Zusammenhang mit der Lösung komplexer Differentialgleichungssysteme für komplexe Matrizen herleitete, die aber auch schon 1868 Karl Weierstraß in seiner Behandlung bilinearer Formen im Komplexen bekannt war. Die jordansche Normalform ist ein einfacher Vertreter der Äquivalenzklasse der zu einer trigonalisierbaren Matrix ähnlichen Matrizen. Die Trigonalisierbarkeit ist gleichbedeutend damit, dass das charakteristische Polynom der Matrix vollständig in Linearfaktoren zerfällt. Matrizen über einem algebraisch abgeschlossenen Körper sind immer trigonalisierbar und daher immer ähnlich einer jordanschen Normalform.

Für jede lineare Abbildung eines endlichdimensionalen Vektorraums, deren charakteristisches Polynom vollständig in Linearfaktoren zerfällt, kann eine Vektorraumbasis gewählt werden, so dass die Abbildungsmatrix, die die Abbildung bezüglich dieser Basis beschreibt, jordansche Normalform hat.

Für jede beliebige, auch nicht trigonalisierbare Matrix liefert die rationale Normalform oder Frobenius-Normalform einen standardisierten Repräsentanten der Ähnlichkeitsklasse dieser Matrix.

Definition

Die jordansche Normalform zu einer quadratischen n\times n-Matrix A über den komplexen Zahlen \mathbb {C} ist eine Matrix J in der folgenden Blockdiagonalform:

J={\begin{pmatrix}J_{1}&&0\\&\ddots &\\0&&J_{k}\end{pmatrix}}=Q^{{-1}}AQ.

Die Matrix Q ist die Matrix der Eigenvektoren und Hauptvektoren, aus denen sie spaltenweise besteht. Q^{-1} bezeichnet dabei die inverse Matrix von Q. Die Matrizen J_i heißen Jordanblöcke und sie sind Bidiagonalmatrizen mit der folgenden Form:

J_{j}={\begin{pmatrix}\lambda _{j}&1&&&0\\&\lambda _{j}&1&&\\&&\ddots {}&\ddots {}\\&&&\lambda _{j}&1\\0&&&&\lambda _{j}\end{pmatrix}}\in \mathbb{C} ^{{s_{j}\times s_{j}}}.

Die \lambda _{j} sind dabei die Eigenwerte von A. Zu jedem Eigenwert \lambda _{j} gibt es seiner geometrischen Vielfachheit entsprechend viele Jordanblöcke. Die geometrische Vielfachheit ist dabei die Dimension des Eigenraums zum Eigenwert \lambda _{j}. Die Gesamtdimension aller Jordanblöcke eines Eigenwertes entspricht seiner algebraischen Vielfachheit, d.h. seiner Vielfachheit im charakteristischen Polynom.

In einem Jordanblock sind die sogenannten Jordanketten „gespeichert“ (Hauptvektor). Besteht A z.B. nur aus einem Jordanblock mit Eigenwert \lambda und bezeichne v_{l} einen Hauptvektor l-ter Stufe. Dabei ist v_{1} ein Eigenvektor zum Eigenwert \lambda und es gilt (A-\lambda E)v_{1}=0 und (A-\lambda E)v_{l}=v_{{l-1}} für l=2,\dots ,n. Es gilt dann Av_{1}=\lambda v_{1} und Av_{l}=v_{{l-1}}+\lambda v_{l} für l=2,\dots ,n, das heißt, die Abbildungsmatrix bezüglich der Basis (v_{1},\dotsc ,v_{n}) ist tatsächlich ein Jordanblock.

Es existiert noch die alternative Darstellung der Jordanblöcke mit 1 in der unteren Nebendiagonalen.

Im Spezialfall einer diagonalisierbaren Matrix ist die jordansche Normalform eine Diagonalmatrix.

Form der Transformationsmatrix

Es seien v_{{j,1}},\ldots ,v_{{j,l}},\ldots ,v_{{j,s_{j}}} Hauptvektoren der jeweils l-ten Stufe, wobei s_j die Dimension des j-ten Jordanblocks ist, j=1,\dotsc ,k.

Dann ist Q, definiert durch

Q:=(v_{{1,1}}|\ldots |v_{{1,s_{1}}}|\ldots |v_{{k,1}}|\ldots |v_{{k,s_{k}}})

eine Transformationsmatrix, die mittels Q^{{-1}}AQ=J die Jordan-Normalform J von A herstellt.

In Worten: Die Spalten vonQ sind die Eigenvektoren mit den dazugehörigen Hauptvektoren in der Reihenfolge der dazugehörigen Jordanblöcke. Allerdings ist  Q nicht eindeutig bestimmt.

Algorithmus zur Bestimmung einer komplexen jordanschen Normalform

Für die jordansche Normalform eines Endomorphismus u:V\to V eines n-dimensionalen \mathbb {C} -Vektorraums V wählt man eine Basis B=(b_{{1}},\ldots ,b_{{n}}) des Vektorraums V und berechnet die jordansche Normalform der Abbildungsmatrix A=M_{{B}}(u) von u bezüglich der Basis B.

Im Folgenden wird daher V={\mathbb  {C}}^{n} gesetzt und die komplexe jordansche Normalform einer quadratischen Matrix A\in \mathbb{C} ^{{n\times n}} bestimmt. Die Einheitsmatrix wird mit E_{n} bezeichnet.

Bestimmung der Eigenwerte

Mit Hilfe des charakteristischen Polynoms

\chi _{A}=\det \left(\lambda E_{n}-A\right)

errechnet man aus seinen Nullstellen die paarweise verschiedenen Eigenwerte

\lambda _{1},\ldots ,\lambda _{k}\in \mathbb{C} \,.

Die Eigenwerte werden hier also nicht ihrer Vielfachheit entsprechend aufgeführt.

Bestimmung der Größe der Jordanblöcke

Hierfür müssen zunächst die Dimensionen der verallgemeinerten Eigenräume bestimmt werden. Das heißt, man berechnet für alle 1\leq j\leq k die Zahlen

a_s := \dim \ker\,(A - \lambda_{j}E_n)^{s}, \quad s \in \mathbb{N}_0\,.

Insbesondere ist stets a_{0}=0 und a_{1} ist gerade die geometrische Vielfachheit des Eigenwerts \lambda _{j}. Die Dimension des Kerns kann mit Hilfe des Dimensionssatzes aus dem Rang berechnet werden, der beispielsweise mit dem gaußschen Algorithmus bestimmt werden kann.

Die Folge der a_{s} ist monoton wachsend und wird ab einem bestimmten Wert für s stationär, spätestens bei der algebraischen Vielfachheit des Eigenwertes im charakteristischen Polynom. Die Anzahl der Jordanblöcke der Größe s zum Eigenwert \lambda _{{j}} lässt sich dann mit Hilfe der Formel

2a_{{s}}-a_{{s-1}}-a_{{s+1}}

berechnen. Außerdem gibt a_{1} die Gesamtzahl der zu diesem Eigenwert gehörigen Jordanblöcke an.

Komplexe jordansche Normalform

Die erhaltenen Jordanblöcke schreibt man in eine Matrix und erhält die komplexe jordansche Normalform einer Matrix. Haben alle Blöcke die Größe 1, liegt der Spezialfall einer Diagonalmatrix vor, und A ist somit diagonalisierbar.

Das Minimalpolynom g\in {\mathbb  C}[X] von A erhält man aus g=\prod_{j=1}^{k}(X-\lambda_{j})^{m_{j}}, worin m_{j} die Größe des größten Jordanblocks zum Eigenwert \lambda _{{j}} bezeichnet.

Die jordansche Normalform ist bis auf die Reihenfolge der Jordanblöcke eindeutig bestimmt. Sofern alle Eigenwerte in \mathbb {K} liegen, sind zwei Matrizen, welche dieselbe jordansche Normalform haben, zueinander ähnlich.

Beispiel

Man betrachte die Matrix A\in \mathbb{C} ^{{5\times 5}}, die wie folgt definiert ist

A:={\begin{pmatrix}25&-16&30&-44&-12\\13&-7&18&-26&-6\\-18&12&-21&36&12\\-9&6&-12&21&6\\11&-8&15&-22&-3\end{pmatrix}}.

Ihr charakteristisches Polynom lautet \chi (A)=(X-3)^{{5}}. Somit besitzt diese Matrix genau einen Eigenwert, nämlich 3. Mit der Abkürzung B:=A-3E_{5} werden nun die a_{{s}} bestimmt:

Es gilt {\mathrm  {rg}}(B)=2. Somit ist a_{{1}}=\dim(V)-{\mathrm  {rg}}(B)=5-2=3.

Weiterhin ist B^{2} die Nullmatrix, also gilt {\mathrm  {rg}}(B^{2})=0 und somit a_{2}=5-0=5 und die Folge a_{s} wird ab dieser Stelle stationär.

Damit folgt:

Es gibt a_{1}=3 Jordanblöcke, und davon

2a_{{1}}-a_{{0}}-a_{{2}}=6-0-5=1 Jordanblock mit Größe 1 und

2a_{{2}}-a_{{1}}-a_{{3}}=10-3-5=2 Jordanblöcke mit Größe 2.

Somit ist {\begin{pmatrix}3&1&0&0&0\\0&3&0&0&0\\0&0&3&1&0\\0&0&0&3&0\\0&0&0&0&3\end{pmatrix}} die jordansche Normalform von A. Das Minimalpolynom von A ist (X-3)^{{2}}.

Bestimmung einer Basistransformation zur komplexen jordanschen Normalform

Nun soll eine Basistransformationsmatrix P\in {{\rm {GL}}}(n;{\mathbb  {C}}) bestimmt werden, die

J=P^{{-1}}AP

erfüllt. Sie ist durch diese Gleichung bekanntlich nicht eindeutig bestimmt. Das Standard-Verfahren verwendet die vorherige Kenntnis der komplexen jordanschen Normalform J.

Ein Standard-Verfahren

Ein gängiges Verfahren, um eine Basistransformation zu erhalten, ist das folgende: Man bestimme (wie auch bei obigem naiven Ansatz) zunächst die Jordannormalform J. Dann hat man insbesondere schon alle Eigenwerte \lambda berechnet sowie die Kerne {{\rm {Kern}}}(A-\lambda I)^{k} für alle 1\leq k\leq m(\lambda ), worin m(\lambda )\in {\mathbb  N} die Dimension des größten Jordanblocks zum Eigenwert \lambda bezeichnet. Anschließend arbeite man zur Bestimmung einer regulären Matrix P mit J=P^{{-1}}AP die Blöcke nacheinander ab. Dabei ist zu beachten, dass man bei Jordanblöcken zum selben Eigenwert stets vom größten Block zum kleinsten Block vorgeht.

Zu jedem Block der Größe s und Eigenwert \lambda werden s Spalten der Basistransformationsmatrix v^{1},\ldots ,v^{s} nach einem bestimmten Schema bestimmt. Wenn der Block in J die Spalten m,\ldots ,m+s-1 belegt, so werden die Vektoren v^{1},\ldots ,v^{s} in P ebenso (von links nach rechts) in die Spalten m,\ldots ,m+s-1 eingefügt. Die Vektoren v^{1},\ldots ,v^{s} werden nun wie folgt bestimmt:

Nachdem man auf obige Weise alle Jordanblöcke abgearbeitet hat, wurden am Ende alle Spalten von P aufgefüllt. Es gilt: P ist regulär und erfüllt P^{{-1}}AP=J, und ihre Spalten bilden eine Basis, bezüglich derer A die Darstellung J besitzt.

Wird die alternative Darstellung der Jordanblöcke gewählt, d.h. mit 1 in der unteren Nebendiagonalen, muss lediglich die Reihenfolge der Basisvektoren pro Jordanblock umgekehrt werden.

Beispiel

Als erläuterndes Beispiel betrachte man hierzu die Matrix

A:={\begin{pmatrix}25&-16&30&-44&-12\\13&-7&18&-26&-6\\-18&12&-21&36&12\\-9&6&-12&21&6\\11&-8&15&-22&-3\end{pmatrix}}

wie oben. Es gilt

{{\rm {Kern}}}(A-3I)={{\rm {Span}}}\left(\left\{{\begin{pmatrix}-2\\1\\2\\0\\0\\\end{pmatrix}},{\begin{pmatrix}2\\0\\0\\1\\0\\\end{pmatrix}},{\begin{pmatrix}2\\2\\0\\0\\1\\\end{pmatrix}}\right\}\right) und {{\rm {Kern}}}(A-3I)^{2}={\mathbb  {C}}^{5}.

Ihre Jordannormalform lautet

J={\begin{pmatrix}3&1&0&0&0\\0&3&0&0&0\\0&0&3&1&0\\0&0&0&3&0\\0&0&0&0&3\end{pmatrix}}.

Man beginne mit dem ersten Jordanblock der Dimension 2. Dazu wähle man

v^{2}\in {{\rm {Kern}}}(A-3I)^{2}\setminus {{\rm {Kern}}}(A-3I)^{1}

beliebig, beispielsweise v^{2}:={\begin{pmatrix}1\\0\\0\\0\\0\\\end{pmatrix}}. Dann ist v^{1}:=(A-3I)v^{2}={\begin{pmatrix}22\\13\\-18\\-9\\11\\\end{pmatrix}} zu wählen. Daraus erhält man P:={\begin{pmatrix}22&1&*&*&*\\13&0&*&*&*\\-18&0&*&*&*\\-9&0&*&*&*\\11&0&*&*&*\end{pmatrix}}. Nun gehe man zum zweiten Jordanblock der Größe 2 über. Man wähle nun

w^{2}\in {{\rm {Kern}}}(A-3I)^{2}\setminus {{\rm {Span}}}({{\rm {Kern}}}(A-3I)^{1}\cup \{v^{2}\})

beliebig, beispielsweise w^{2}:={\begin{pmatrix}0\\1\\0\\0\\0\\\end{pmatrix}}. Dann ist w^{1}:=(A-3I)w^{2}={\begin{pmatrix}-16\\-10\\12\\6\\-8\\\end{pmatrix}}, und man landet bei P={\begin{pmatrix}22&1&-16&0&*\\13&0&-10&1&*\\-18&0&12&0&*\\-9&0&6&0&*\\11&0&-8&0&*\end{pmatrix}}. Schließlich ist der letzte Jordanblock (der Größe 1) an der Reihe. Man wähle hierzu

x^{1}\in {{\rm {Kern}}}(A-3I)^{1}\setminus {{\rm {Span}}}(\{v^{1},w^{1}\})

beliebig, beispielsweise x^{1}:={\begin{pmatrix}2\\0\\0\\1\\0\\\end{pmatrix}}. Dann ist P={\begin{pmatrix}22&1&-16&0&2\\13&0&-10&1&0\\-18&0&12&0&0\\-9&0&6&0&1\\11&0&-8&0&0\end{pmatrix}} eine reguläre Matrix mit J=P^{{-1}}AP.

Reelle jordansche Normalform

Betrachtet man reelle Matrizen, so zerfällt deren charakteristisches Polynom im Allgemeinen nicht mehr vollständig in Linearfaktoren, sondern nur noch in irreduzible Faktoren, die in diesem Fall stets lineare oder quadratische Faktoren sind. Es stellt sich nun die Frage nach einer Normalform, wenn man ausschließlich reelle Basistransformationen zulässt.

Zu einem quadratischen irreduziblen Faktor (\lambda -a_{j})^{2}+b_{j}^{2} mit b_{j}>0 definiert man als Jordanblock

J_{j}={\begin{pmatrix}a_{j}&b_{j}&&&&&&0\\-b_{j}&a_{j}&1&&&&&\\&&a_{j}&b_{j}&&&&\\&&-b_{j}&a_{j}&1&&&\\&&&\ddots {}&\ddots {}&\ddots {}&&\\&&&&\ddots {}&\ddots {}&1&\\&&&&&\ddots {}&a_{j}&b_{j}\\0&&&&&&-b_{j}&a_{j}\\\end{pmatrix}}.

Wir nennen die Anzahl der Zeilen (bzw. Spalten) die Größe dieses Blocks. Dann bezeichnet man

J={\begin{pmatrix}J_{1}&&0\\&\ddots {}&\\0&&J_{k}\end{pmatrix}}=P^{{-1}}AP

als reelle jordansche Normalform. Um sie und eine geeignete reelle Matrix P\in {\mathbb  {R}}^{{n\times n}} zu bestimmen, kann man folgendermaßen vorgehen:

\chi (\lambda )=\prod _{{j=1}}^{k}\left(\lambda -\lambda _{j}\right)^{{\mu _{j}}}\cdot \prod _{{j=1}}^{l}\left(\left(\lambda -a_{j}\right)^{2}+b_{j}^{2}\right)^{{\nu _{j}}},
wobei \lambda _{1},\ldots ,\lambda _{k}\in {\mathbb  {R}} paarweise verschiedene Eigenwerte mit Vielfachheit \mu _{j}\in {\mathbb  {N}} bezeichnen. Weiter seien darin a_{1},\ldots ,a_{l}\in {\mathbb  {R}}, b_{1},\ldots ,b_{l}>0, \nu _{1},\ldots ,\nu _{l}\in {\mathbb  {N}} und (a_{1},b_{1}),\ldots ,(a_{l},b_{l}) paarweise verschieden.
K_{{j,m}}:={{\rm {Kern}}}_{{{\mathbb  {R}}}}(A-\lambda _{j}E)^{m} für m=1,2,\ldots ,m_{j},
worin m_{j}\leq \mu _{j} die kleinste natürliche Zahl ist mit \dim _{{{\mathbb  {R}}}}K_{{j,m_{j}}}=\mu _{j}. Analog bestimme man für jedes j\in \{1,\ldots ,l\}
K^{{j,m}}:={{\rm {Kern}}}_{{{\mathbb  {R}}}}\left(\left(A-a_{j}E\right)^{2}+b_{j}^{2}E\right)^{m} für m=1,2,\ldots ,n_{j},
worin n_{j}\leq \nu _{j} die kleinste natürliche Zahl ist mit \dim _{{{\mathbb  {R}}}}K^{{j,n_{j}}}=2\nu _{j}.
Zudem setzen wir K_{{j,0}}:=K^{{j,0}}:=\{0\}.
Außerdem ist \mu _{j} die Summe der Jordanblockgrößen zum Eigenwert \lambda _{j} und 2\nu _{j} die Summe der Jordanblockgrößen zum Faktor (\lambda -a_{j})^{2}+b_{j}^{2}. Aus diesen Angaben kann man eindeutig die jordansche Normalform J bestimmen.

Ein Verfahren, um eine Basistransformation zu erhalten, ist das folgende:

Dann setze man für t=2m,\ldots ,2 sukzessiv v^{{t-1}}:=\left\{{\begin{array}{ll}{\frac  {1}{b_{j}}}(A-a_{j}E)v^{t}\ ,&{\textrm  {falls}}\ t\ {\textrm  {gerade}}\ ,\\(A-a_{j}E)v^{t}+b_{j}v^{{t+1}}\ ,&{\textrm  {falls}}\ t\ {\textrm  {ungerade}}\ .\\\end{array}}\right.
Schließlich setzt man P wie gehabt aus den Vektoren v^{1},\ldots ,v^{{2m}} zusammen.

Beispiel

Man betrachte die Matrix B\in M_{{5}}({\mathbb  R}), die wie folgt definiert ist

B:={\begin{pmatrix}6&-2&6&1&1\\1&-1&2&1&-2\\-2&0&-1&0&-1\\-1&0&-2&2&-1\\-4&4&-6&-2&3\end{pmatrix}}.

Ihr charakteristisches Polynom lautet \chi (\lambda )=((\lambda -2)^{2}+1)^{{2}}(\lambda -1), wobei (\lambda -2)^{2}+1 irreduzibel über \mathbb {R} ist. Nun berechnen wir die jordansche Normalform:

{{\rm {Kern}}}_{{{\mathbb  {R}}}}(B-1E)={{\rm {Span}}}_{{{\mathbb  {R}}}}\left(\left\{{\begin{pmatrix}1\\-1\\-1\\-1\\0\\\end{pmatrix}}\right\}\right).

Dieser Kern hat die Dimension 1. Also gibt es nur einen Jordanblock der Größe \geq 1. Andererseits muss die Summe der Jordanblockgrößen 1 sein (die Potenz von \lambda -1), so dass es genau einen Jordanblock zum Eigenwert 1 gibt, und er hat die Größe 1. Weiter hat

{{\rm {Kern}}}_{{{\mathbb  {R}}}}\left(\left(B-2E\right)^{2}+1^{2}E\right)={{\rm {Span}}}_{{{\mathbb  {R}}}}\left(\left\{{\begin{pmatrix}2\\0\\-1\\-1\\-1\\\end{pmatrix}},{\begin{pmatrix}0\\1\\0\\1\\-1\\\end{pmatrix}}\right\}\right)

die Dimension 2, so dass es demzufolge nur {\tfrac  {1}{2}}\cdot 2=1 Jordanblock der Größe \geq 2 gibt. Da die Summe der Jordanblockgrößen 4 sein muss (das Doppelte der Potenz von (\lambda -2)^{2}+1), ergibt sich, dass dieser eine Jordanblock die Größe 4 besitzt. Außerdem errechnen wir

{{\rm {Kern}}}_{{{\mathbb  {R}}}}\left(\left(B-2E\right)^{2}+1^{2}E\right)^{2}={{\rm {Span}}}_{{{\mathbb  {R}}}}\left(\left\{{\begin{pmatrix}1\\0\\0\\0\\-1\\\end{pmatrix}},{\begin{pmatrix}0\\1\\0\\0\\-1\\\end{pmatrix}},{\begin{pmatrix}0\\0\\1\\0\\-1\\\end{pmatrix}},{\begin{pmatrix}0\\0\\0\\1\\0\\\end{pmatrix}}\right\}\right).

Somit ist J:={\begin{pmatrix}1&&&&\\&2&1&&\\&-1&2&1&\\&&&2&1\\&&&-1&2\end{pmatrix}} die reelle jordansche Normalform von B.

Zum Vergleich, die komplexe jordansche Normalform lautet J_{{{\mathbb  {C}}}}:={\begin{pmatrix}1&&&&\\&2+i&1&&\\&&2+i&&\\&&&2-i&1\\&&&&2-i\end{pmatrix}}.

Zum Berechnen einer Basistransformationsmatrix beginne man mit dem ersten reellen Eigenwert und dann mit dem (ersten) Jordanblock der Dimension 1. Man wähle

u^{1}\in {{\rm {Kern}}}_{{{\mathbb  {R}}}}(B-1E)^{1}\setminus {{\rm {Kern}}}_{{{\mathbb  {R}}}}(B-1E)^{0}

beliebig, also beispielsweise u^{1}:={\begin{pmatrix}1\\-1\\-1\\-1\\0\\\end{pmatrix}}. Daraus erhält man P:={\begin{pmatrix}1&*&*&*&*\\-1&*&*&*&*\\-1&*&*&*&*\\-1&*&*&*&*\\0&*&*&*&*\end{pmatrix}}.

Nun gehe man zum ersten irreduziblen Faktor (komplexen Eigenwert) und dann zum Jordanblock der Größe 4 über. Dazu wähle man

v^{4}\in {{\rm {Kern}}}_{{{\mathbb  {R}}}}((B-2E)^{2}+1^{2}E)^{2}\setminus {{\rm {Kern}}}_{{{\mathbb  {R}}}}((B-2E)^{2}+1^{2}E)^{1}

beliebig, beispielsweise v^{4}:={\begin{pmatrix}0\\0\\0\\1\\0\\\end{pmatrix}}. Dann ist v^{3}:={\frac  {1}{1}}(B-2E)v^{4}={\begin{pmatrix}1\\1\\0\\0\\-2\\\end{pmatrix}}, v^{2}:=(B-2E)v^{3}+bv^{4}={\begin{pmatrix}0\\2\\0\\2\\-2\\\end{pmatrix}} und v^{1}:={\frac  {1}{1}}(B-2E)v^{2}={\begin{pmatrix}-4\\0\\2\\2\\2\\\end{pmatrix}} zu wählen. Daraus erhält man: P:={\begin{pmatrix}1&-4&0&1&0\\-1&0&2&1&0\\-1&2&0&0&0\\-1&2&2&0&1\\0&2&-2&-2&0\end{pmatrix}}. P ist eine reguläre Matrix mit J=P^{{-1}}BP.

Jordansche Normalform in allgemeinen Körpern

Die jordansche Normalform kann noch weiter verallgemeinert werden auf allgemeine Körper. In diesem Zusammenhang wird sie häufig auch als Weierstraß-Normalform (bzw. Frobenius-Normalform) bezeichnet. Dies erlaubt eine eindeutige Matrixdarstellung von Endomorphismen von endlichdimensionalen Vektorräumen, bei der sich alle ähnlichen Endomorphismen durch eine eindeutige Matrix darstellen lassen. So können ähnliche lineare Abbildungen identifiziert werden. Das Lemma von Frobenius charakterisiert zueinander ähnliche Matrizen durch die Elementarteiler ihrer charakteristischen Matrizen und liefert die Frobenius-Normalform als Normalform des Vektorraums unter der Operation eines Polynomrings.

Durch die Darstellung in der Weierstraß-Normalform ist der Aufbau des Minimalpolynoms sofort erkennbar und das charakteristische Polynom leicht zu berechnen.

Anwendung bei linearen Differentialgleichungssystemen erster Ordnung mit konstanten Koeffizienten

Gegeben sei ein lineares Differentialgleichungssystem (von n Gleichungen) erster Ordnung mit konstanten Koeffizienten

y'=A\cdot y+g(x)

durch eine Matrix A\in {\mathbb  {C}}^{{n\times n}} und eine stetige Funktion g:{\mathbb  {R}}\rightarrow {\mathbb  {C}}^{n}. Es ist bekannt, dass die eindeutige Lösung des Anfangswertproblems

y(x_{0})=y_{0}\in {\mathbb  {C}}^{n}

gegeben ist durch

y(x)=e^{{(x-x_{0})A}}\cdot y_{0}+\int _{{x_{0}}}^{x}e^{{(x-t)A}}g(t){{\rm {d}}}t\ ,

worin

\exp(B):=e^{B}:=\sum _{{k=0}}^{\infty }{\frac  {1}{k!}}B^{k} für B\in {\mathbb  {C}}^{{n\times n}}

die Matrixexponentialfunktion bezeichnet. Man beachte:

\exp \left(t\cdot {\begin{pmatrix}\lambda &1&0&\cdots &0\\0&\lambda &1&\ddots &0\\\vdots &\ddots &\ddots &\ddots &\vdots \\0&\cdots &0&\lambda &1\\0&\cdots &\cdots &0&\lambda \end{pmatrix}}\right)=e^{{t\lambda }}\cdot {\begin{pmatrix}1&{\frac  {t^{1}}{1!}}&{\frac  {t^{2}}{2!}}&\cdots &{\frac  {t^{{n-1}}}{(n-1)!}}\\0&1&{\frac  {t^{1}}{1!}}&\cdots &{\frac  {t^{{n-2}}}{(n-2)!}}\\\vdots &\ddots &\ddots &\ddots &\vdots \\0&\cdots &0&1&{\frac  {t^{1}}{1!}}\\0&\cdots &\cdots &0&1\end{pmatrix}}.
\exp(t\cdot {{\rm {diag}}}(J_{1},\ldots ,J_{m}))={{\rm {diag}}}(\exp(tJ_{1}),\ldots ,\exp(tJ_{m})).
\exp(tA)=\exp(t\cdot PJP^{{-1}})=P\cdot \exp(tJ)\cdot P^{{-1}}.

Mit anderen Worten: Kennt man eine Darstellung A=PJP^{{-1}} mit der komplexen jordanschen Normalform J, so kann man \exp(tA) für jedes t \in \mathbb{R} explizit ausrechnen, so dass zum Bestimmen von

y(x)=e^{{(x-x_{0})A}}\cdot y_{0}+\int _{{x_{0}}}^{x}e^{{(x-t)A}}g(t){{\rm {d}}}t\

nur noch das Integrationsproblem zu lösen ist, welches im homogenen Fall g = 0 völlig entfällt.

Siehe auch

Trenner
Basierend auf einem Artikel in: externer Link Wikipedia.de
Seitenende
Seite zurück
©  biancahoegel.de
Datum der letzten Änderung: Jena, den: 27.12. 2020