Numerische Lineare Algebra - TU Chemnitz
Numerische Lineare Algebra - TU Chemnitz Numerische Lineare Algebra - TU Chemnitz
Inhalt 1 Organisatorisches 2 Einleitung 3 Das QR-Verfahren 3.1 Vektoriteration 3.2 Unterraumiteration 3.3 Reduktion auf Hessenberg-Gestalt 3.4 QR Iteration mit impliziten Shifts 3.5 Konvergenz der Unterraumiteration 3.6 Konvergenz des QR-Verfahrens 4 Hermitesche Eigenwertaufgaben 5 Anhang: Grundlagen aus der Linearen Algebra Oliver Ernst (Numerische Mathematik) Numerische Lineare Algebra Wintersemester 2013/14 27 / 177
Vektoriteration Der Prototyp aller numerischer Eigenwertverfahren ist die sog. Vektoriteration nach von Mises, auch Potenzmethode (engl. power method) genannt. Es beruht auf der Tatsache, dass für (nahezu) jeden Vektor q ∈ C n die Vektorfolge q, Aq, A 2 q, A 3 q, . . . zunehmend in Richtung eines Eigenvektors von A zum Eigenwert ρ(A) zeigt. In seiner einfachsten Form lautet das Verfahren wie folgt: Algorithmus 1 : Vektoriteration nach von Mises. Gegeben : A ∈ C n×n , q 0 ∈ C n . 1 for k = 1 to . . . do 2 z k ← Aq k−1 3 q k ← z k /‖z k ‖ 2 4 µ k ← q H k Aq k Wie man sieht, wird, allein schon um Gleitpunktüber-/unterlauf zu vermeiden, der Vektor zusätzlich in jedem Schritt normiert . Oliver Ernst (Numerische Mathematik) Numerische Lineare Algebra Wintersemester 2013/14 28 / 177
- Seite 1 und 2: Numerische Lineare Algebra Oliver E
- Seite 3 und 4: Inhalt II 5.2 Störungstheorie 5.3
- Seite 5: Einleitung Nach dem Satz von Abel-R
- Seite 9 und 10: Vektoriteration Beispiel Wir betrac
- Seite 11 und 12: Vektoriteration Beispiel q 0 = v 1
- Seite 13 und 14: Unterraumiteration Vektoriteration
- Seite 15 und 16: Unterraumiteration Invariante Unter
- Seite 17 und 18: Unterraumiteration Verschieben des
- Seite 19 und 20: Unterraumiteration Simultane Iterat
- Seite 21 und 22: Unterraumiteration QR-Iteration Das
- Seite 23 und 24: Unterraumiteration Invarianzerkennu
- Seite 25 und 26: Unterraumiteration Rekursive Darste
- Seite 27 und 28: Inhalt 1 Organisatorisches 2 Einlei
- Seite 29 und 30: Reduktion auf Hessenberg-Gestalt Al
- Seite 31 und 32: Reduktion auf Hessenberg-Gestalt Mi
- Seite 33 und 34: Inhalt 1 Organisatorisches 2 Einlei
- Seite 35 und 36: QR-Schritt mit impliziten Shifts QR
- Seite 37 und 38: QR-Schritt mit impliziten Shifts Un
- Seite 39 und 40: QR-Schritt mit impliziten Shifts Im
- Seite 41 und 42: Konvergenz Praktische Koordinaten S
- Seite 43 und 44: Konvergenz Praktische Koordinaten D
- Seite 45 und 46: Konvergenz Konvergenzsatz II Satz 3
- Seite 47 und 48: Inhalt 1 Organisatorisches 2 Einlei
- Seite 49 und 50: QR Konvergenz Satz 3.12 Beim QR-Ver
- Seite 51 und 52: QR Konvergenz In der Praxis besitze
Vektoriteration<br />
Der Prototyp aller numerischer Eigenwertverfahren ist die sog. Vektoriteration nach<br />
von Mises, auch Potenzmethode (engl. power method) genannt. Es beruht auf der<br />
Tatsache, dass für (nahezu) jeden Vektor q ∈ C n die Vektorfolge<br />
q, Aq, A 2 q, A 3 q, . . .<br />
zunehmend in Richtung eines Eigenvektors von A zum Eigenwert ρ(A) zeigt. In<br />
seiner einfachsten Form lautet das Verfahren wie folgt:<br />
Algorithmus 1 : Vektoriteration nach von Mises.<br />
Gegeben : A ∈ C n×n , q 0 ∈ C n .<br />
1 for k = 1 to . . . do<br />
2 z k ← Aq k−1<br />
3 q k ← z k /‖z k ‖ 2<br />
4 µ k ← q H k Aq k<br />
Wie man sieht, wird, allein schon um Gleitpunktüber-/unterlauf zu vermeiden, der<br />
Vektor zusätzlich in jedem Schritt normiert .<br />
Oliver Ernst (<strong>Numerische</strong> Mathematik) <strong>Numerische</strong> <strong>Lineare</strong> <strong>Algebra</strong> Wintersemester 2013/14 28 / 177