Verallgemeinertes Eigenwertproblem

Das verallgemeinerte Eigenwertproblem ist eine Problemstellung der linearen Algebra.

Definition

Das Problem, zu vorgegebenen Matrizen $A,B\in {ℂ}^{n\times n}$ gewisse Zahlen $\lambda \in ℂ$ und Vektoren $x\in {ℂ}^n$ mit $x\ne 0$ zu bestimmen, sodass

$Ax=\lambda Bx$

gilt, wird in Abgrenzung zum Eigenwertproblem als verallgemeinertes Eigenwertproblem bezeichnet.

Lösungsverfahren

Ist $B$ regulär, so lässt sich das verallgemeinerte Eigenwertproblem auf das gewöhnliche Eigenwertproblem

$B^{-1}Ax=\lambda x$

zurückführen. Dieser Lösungsansatz ist aber i. A. nur von theoretischer Bedeutung, da die Berechnung einer inversen Matrix numerisch oft nicht möglich oder sehr unpraktisch ist. Oftmals lassen sich aus der Aufgabenstellung schon gewisse Informationen über die betrachteten Matrizen sammeln, welche die Berechnung dann vereinfachen können. Sind z. B. $A,B$ symmetrisch und $B$ außerdem positiv definit, so lässt sich die Berechnung wesentlich vereinfachen: Die Matrix $B$ lässt sich mittels der Cholesky-Zerlegung in $B=L{L}^T$ zerlegen. Dann ist $B^{-1}A$ ähnlich zu einer Matrix $H:=L^{-1}A(L^{-1}{)}^T$. Die Inverse von $L$ lässt sich sehr effizient berechnen, da $L$ eine Dreiecksmatrix ist. Bestimmt man nun die Eigenwerte von $H$, so sind dies auch die Eigenwerte von $B^{-1}A$.

Für beliebige Matrizen $A,B$ kann auch der QZ-Algorithmus genutzt werden.

Beispiel

Betrachte das verallgemeinerte Eigenwertproblem

$\left[\begin{array}{ccc}4& 4& 8\\ 4& 0& 4\\ 8& 4& 4\end{array}\right]x=\lambda \left[\begin{array}{ccc}2& -1& 0\\ -1& 2& -1\\ 0& -1& 2\end{array}\right]x$.

Naiver Ansatz

Die Berechnung der Inversen von $B$ ergibt

$B^{-1}=\frac{1}{4}\left[\begin{array}{ccc}3& 2& 1\\ 2& 4& 2\\ 1& 2& 3\end{array}\right]$

und damit

$B^{-1}A=\left[\begin{array}{ccc}7& 4& 9\\ 10& 4& 10\\ 9& 4& 7\end{array}\right]$.

Die Eigenwerte dieser Matrix sind 20,7703 sowie -2 und - 0,7703.

Mittels der Cholesky-Zerlegung

$A,B$ sind symmetrisch und $B$ außerdem positiv definit. Die Cholesky-Zerlegung liefert die Matrix

$L^T=\left[\begin{array}{ccc}1,4142& -0,7071& 0\\ 0& 1,2247& -0,8165\\ 0& 0& 1,1547\end{array}\right]$.

Dann ist $H:=L^{-1}A(L^{-1}{)}^T=\left[\begin{array}{ccc}2,0000& 3,4641& 7,3485\\ 3,4641& 3,3333& 8,0139\\ 7,3485& 8,0139& 12,6667\end{array}\right]$.

Die Eigenwerte dieser Matrix sind wie zu erwarten mit den oben berechneten Eigenwerten identisch.

Literatur

• Vorlage:Literatur
• Josef Stoer, Roland Bulirsch: Numerische Mathematik 2. 5. Auflage, Springer, Berlin/Heidelberg/New York 2005, ISBN 978-3-540-23777-8.