OrthogonalisierungsverfahrenOrthogonalisierung zum Lösen linearer Gleichungssysteme
Beschreibung
Man möchte ein Gleichungssystem Ax=b lösen. Dazu berechnet man eine Zerlegung der Matrix A in das Produkt zweier Matrizen A=QR. Wobei diese Matrizen folgende schöne Eigenschaften haben:
![$ Q^T=Q^{-1} $ $ Q^T=Q^{-1} $](/teximg/5/9/00050995.png)
![$ R=\pmat{r_{11} & \cdots & r_{1n} \\ \vdots & \ddots &\vdots &\\ 0 & \cdots & r_{nn}} $ $ R=\pmat{r_{11} & \cdots & r_{1n} \\ \vdots & \ddots &\vdots &\\ 0 & \cdots & r_{nn}} $](/teximg/1/6/00413561.png)
Die Inverse von Q ist also gleich ihrer Transponierten und R ist eine obere Dreiecksmatrix. Deren Lösung man durch Rückwärtseinsetzen bestimmen kann.
Algorithmus
Gegeben (A,b)
Gesucht die Lösung x von Ax=b
1. Berechne die QR-Zerlegung von A
A->(Q,R)
2. Berechne ![$ c=Q^Tb $ $ c=Q^Tb $](/teximg/9/7/00044479.png)
3. Löse das Dreieckssystem Rx=c
Verfahren zur praktischen Umsetzung
Householder Spiegelungen
Givens Drehungen
Gram-Schmidt Orthogonalisierung
Matheraum Links
Fragen zum Artikel
Fragen zu diesem Artikel können im Matheforum gestellt werden.
|