Numerische Lineare Algebra - TU Chemnitz
Numerische Lineare Algebra - TU Chemnitz Numerische Lineare Algebra - TU Chemnitz
Unterraumiteration Simultane Iteration Aufgrund der oberen Dreiecksstruktur von R k bilden die Spalten q (k) 1 , . . . , q (k) d von Q k nicht nur eine Orthonormalbasis von p(A)S k−1 , es gilt vielmehr span{p(A)q (k−1) 1 , . . . , p(A)q (k−1) } = span{q (k) 1 , . . . , q (k) j }, j = 1, . . . , d. j Die simultane Iteration erzeugt somit eine Folge geschachtelter Unterräume der Dimensionen 1 bis d. Ist A diagonalisierbar und sind v 1 , . . . , v n die Eigenvektoren zu den Eigenwerten λ 1 , . . . , λ n von A mit so gilt |p(λ 1 )| ≤ |p(λ 2 )| ≤ · · · ≤ |p(λ n )| span{q (k) 1 , . . . , q (k) j } → span{v n−j+1 , . . . , v n }, j = 1, . . . , d, sofern |p(λ n−j+1 )| > |p(λ n−j )|. Oliver Ernst (Numerische Mathematik) Numerische Lineare Algebra Wintersemester 2013/14 41 / 177
Unterraumiteration QR-Iteration Das QR-Verfahren entspricht simultaner Iteration mit d = n beginnend mit Q 0 = I . Hierbei wird also eine Folge geschachtelter Unterräume erzeugt, die jeweils gegen invariante Unterräume von A konvergieren. Alternativ kann das QR-Verfahren als Folge von Orthonormalbasen betrachtet werden, bezüglich derer im Grenzfall die Matrix A Schur-Form annimmt. Es wird sozusagen das kartesische Koordinatensystem im C n geeignet „zurechtrotiert“. Neben QR-Zerlegung sind auch andere Transformationen gebrauchlich, etwa LR-Zerlegungen. Diese sollten nach Moglichkeit die vorhandene Struktur in A erhalten. Oliver Ernst (Numerische Mathematik) Numerische Lineare Algebra Wintersemester 2013/14 42 / 177
- Seite 1 und 2: Numerische Lineare Algebra Oliver E
- Seite 3 und 4: Inhalt II 5.2 Störungstheorie 5.3
- Seite 5 und 6: Einleitung Nach dem Satz von Abel-R
- Seite 7 und 8: Vektoriteration Der Prototyp aller
- 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: Unterraumiteration Simultane Iterat
- 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
Unterraumiteration<br />
Simultane Iteration<br />
Aufgrund der oberen Dreiecksstruktur von R k bilden die Spalten q (k)<br />
1 , . . . , q (k)<br />
d<br />
von<br />
Q k nicht nur eine Orthonormalbasis von p(A)S k−1 , es gilt vielmehr<br />
span{p(A)q (k−1)<br />
1 , . . . , p(A)q (k−1) } = span{q (k)<br />
1 , . . . , q (k)<br />
j }, j = 1, . . . , d.<br />
j<br />
Die simultane Iteration erzeugt somit eine Folge geschachtelter Unterräume der<br />
Dimensionen 1 bis d.<br />
Ist A diagonalisierbar und sind v 1 , . . . , v n die Eigenvektoren zu den Eigenwerten<br />
λ 1 , . . . , λ n von A mit<br />
so gilt<br />
|p(λ 1 )| ≤ |p(λ 2 )| ≤ · · · ≤ |p(λ n )|<br />
span{q (k)<br />
1 , . . . , q (k)<br />
j } → span{v n−j+1 , . . . , v n }, j = 1, . . . , d,<br />
sofern |p(λ n−j+1 )| > |p(λ n−j )|.<br />
Oliver Ernst (<strong>Numerische</strong> Mathematik) <strong>Numerische</strong> <strong>Lineare</strong> <strong>Algebra</strong> Wintersemester 2013/14 41 / 177