Mo Logo [Home] [Lexikon] [Aufgaben] [Tests] [Kurse] [Begleitmaterial] [Hinweise] [Mitwirkende] [Publikationen] Englische Flagge
Mathematik-Online-Lexikon:

Implementierung der QR-Iteration

A B C D E F G H I J K L M N O P Q R S T U V W X Y Z Übersicht
Ein Schritt $ A \mapsto B$ der tridiagonalen $ QR$-Iteration mit Shift $ \lambda$ berechnet eine $ QR$-Faktorisierung,

$\displaystyle A-\lambda E = QR, $

die die tridiagonale Form erhält und bildet

$\displaystyle B=RQ +\lambda E =Q^{\operatorname t} A Q \,. $

Sie kann folgendermaßen implementiert werden. Zuerst wird eine Givens-Rotation $ Q_1^{\operatorname t} $, basierend auf dem Vektor $ (a_{1,1} - \lambda ,\, a_{2,1})^{\operatorname t}$, auf die ersten beiden Zeilen von links und $ Q_1$ auf die ersten beiden Spalten von rechts angewendet. Dadurch entsteht ein Eintrag in $ (3,1)$, der die tridiagonale Form von $ A$ zerstört. Nun werden Givens-Rotationen $ Q_i^{\operatorname t},\, i=2, \hdots ,n-1$ sukzessive auf die Zeilen $ i$ und $ i+1$ angewendet und ihre Transponierten $ Q_i$ auf die entsprechenden Spalten:

$\displaystyle B= Q_{n-1}^{\operatorname t} \cdots Q_{1}^{\operatorname t} A Q_1 \cdots Q_{n-1} \,. $

Die Transformation $ Q_i^{\operatorname t}$ basiert auf dem Vektor der aktuellen Einträge in den Positionen $ (i:i+1,i-1)$. Folglich stellt die Ähnlichkeitstransformation mit $ Q_i$ die tridiagonale Form bis zur Spalte $ i-1$ wieder her. Dieser Prozess wird nachstehend schematisch illustriert.

$\displaystyle \left( \begin{array}{c} {\mbox{\includegraphics[width=0.2\linewidth,bb=125 619 214 708,clip]{Num_18_2_Bild1}}} \end{array} \right)$ $\displaystyle \stackrel{Q_1}{\longrightarrow} \left( \begin{array}{c} {\mbox... ...8,clip]{Num_18_2_Bild2}}} \end{array} \right) \stackrel{Q_2}{\longrightarrow}$    
$\displaystyle \left( \begin{array}{c} {\mbox{\includegraphics[width=0.2\linewidth,bb=125 619 214 708,clip]{Num_18_2_Bild3}}} \end{array} \right)$ $\displaystyle \stackrel{Q_3}{\longrightarrow} \left( \begin{array}{c} {\mbox... ...214 708,clip]{Num_18_2_Bild4}}} \end{array} \right) \longrightarrow \cdots$    

Änderungen in den Einträgen sind durch verschiedene Symbole gekennzeichnet. Darüber hinaus werden die Einträge, die die folgende Givens-Rotation bestimmen, hervorgehoben. Diese Transformation annulliert den Eintrag in der eingekreisten Position.

Falls zufällig beide Einträge in Position $ (i:i+1,i-1)$ Null sind, kann die $ i$-te Transformation entfallen. In diesem Fall ist der Nebendiagonaleintrag $ b_{i,i-1}$ Null, und die Dimension des Eigenwertproblems kann reduziert werden.

Download:

(Dateityp: .m, 1.4K,  07.05.2007)
[Erläuterungen] [Verweise]
  automatisch erstellt am 19.  8. 2013