QR-Zerlegung
# QR-Zerlegung

Die QR-Zerlegung oder QR-Faktorisierung ist ein Begriff aus den mathematischen Teilgebieten der linearen Algebra und Numerik. Man bezeichnet damit die Zerlegung einer Matrix $A$ in das Produkt

    $A=Q\cdot R$

zweier anderer Matrizen, wobei $Q$ eine orthogonale bzw. unitäre Matrix und $R$ eine obere Dreiecksmatrix ist. Die QR-Zerlegung ist ein Spezialfall der Iwasawa-Zerlegung.

Eine solche Zerlegung existiert stets und kann mit verschiedenen Algorithmen berechnet werden. Die bekanntesten davon sind

* Householdertransformationen
* Givens-Rotationen
* Gram-Schmidtsches Orthogonalisierungsverfahren.

Das letztere wird üblicherweise in der linearen Algebra benutzt, ist aber in seiner Standardform numerisch instabil. Man kann das Verfahren aber erweitern und numerisch stabilisieren.

## Definition

Jede reelle Matrix $A\in \mathbb {R} ^{m\times n}$, $m\geq n$ besitzt eine (fast – siehe weiter unten) eindeutige reduzierte QR-Zerlegung

    $A={\hat {Q}}\cdot {\hat {R}}$

als Produkt einer in den Spalten orthogonalen Matrix ${\hat {Q}}\in \mathbb {R} ^{m\times n}$ und einer oberen Dreiecksmatrix ${\hat {R}}\in \mathbb {R} ^{n\times n}$.

Diese Lösung ist erweiterbar zu einer vollständigen QR-Zerlegung

    $A=Q\cdot R$,

indem man ${\hat {Q}}$ mit weiteren orthogonalen Spalten ${\tilde {Q}}$ zu einer quadratischen $m\times m$-Matrix erweitert, und an ${\hat {R}}$ unten Nullzeilen anfügt, so dass als Matrixprodukt eine $m\times n$-Matrix entsteht:

    $Q\cdot R={\begin{pmatrix}{\hat {Q}}&{\tilde {Q}}\end{pmatrix}}\cdot {\begin{pmatrix}{\hat {R}}\\0\end{pmatrix}}={\hat {Q}}\cdot {\hat {R}}.$

Die QR-Zerlegung ist eindeutig für $m\geq n$ und $\operatorname {Rang} (A)=n$, wenn man die Vorzeichen der Diagonalelemente von $R,{\hat {R}}$ vorgibt – üblicherweise wählt man alle positiv.

Anmerkungen: Eine reelle Matrix $Q$ heißt orthogonal, wenn das Matrixprodukt mit ihrer transponierten Matrix $Q^{T}$ die Einheitsmatrix $I$ ergibt, also $Q^{T}Q=I$ ist. Im allgemeineren Fall einer komplexen Matrix $A$ ist die Matrix $Q$ eine unitäre Matrix, das heißt, das Matrixprodukt mit ihrer adjungierten Matrix $Q^{H}={\overline {Q}}^{T}={\overline {Q^{T}}}$ ergibt die Einheitsmatrix $I$, also $Q^{H}Q=I$.

### Beispiel

Für die reelle Matrix $A={\begin{pmatrix}12&-51&4\\6&167&-68\\-4&24&-41\end{pmatrix}}$ ergibt sich die QR-Zerlegung

    $Q={\begin{pmatrix}6/7&-69/175&-58/175\\3/7&158/175&6/175\\-2/7&6/35&-33/35\end{pmatrix}}$, $R={\begin{pmatrix}14&21&-14\\0&175&-70\\0&0&35\end{pmatrix}}$,

denn es gilt:

    $QR={\begin{pmatrix}6/7&-69/175&-58/175\\3/7&158/175&6/175\\-2/7&6/35&-33/35\end{pmatrix}}{\begin{pmatrix}14&21&-14\\0&175&-70\\0&0&35\end{pmatrix}}={\begin{pmatrix}12&-51&4\\6&167&-68\\-4&24&-41\end{pmatrix}}=A$

## Anwendung

Die QR-Zerlegung spielt in vielen Verfahren der numerischen Mathematik eine wichtige Rolle, beispielsweise um eine orthogonale oder unitäre Basis zu bestimmen oder um lineare Ausgleichsprobleme zu behandeln. Sie ist integraler Bestandteil des QR-Algorithmus und der Unterraumiteration zur Berechnung von Eigenwerten einer Matrix.

### Lösung regulärer oder überbestimmter Gleichungssysteme

Um die Lösung $x\in \mathbb {R} ^{n}$ eines linearen Gleichungssystems $Ax=b$ mit Matrix $A\in \mathbb {R} ^{m\times n},\ m\geq n$ von vollem Rang zu bestimmen, sind folgende drei Schritte durchzuführen:

* Bestimme eine QR-Zerlegung der Matrix $A$.
* Berechne $z=Q^{T}b\in \mathbb {R} ^{n}$, üblicherweise unter Benutzung der Faktorisierung von $Q$ aus Schritt 1.
* Löse $Rx=z$ durch Rückwärtseinsetzen.

Für $m=n$ ist dies eine Alternative zur LR-Zerlegung, sie hat den doppelten Aufwand der LR-Zerlegung, ist aber möglicherweise numerisch stabiler.

Im Fall $m>n$ gibt es im Gleichungssystem mehr Gleichungen als Variablen und es liegt ein überbestimmtes Gleichungssystem vor. Hier wird $x$ durch Lösung des Ausgleichproblems nach der Methode der kleinsten Quadrate (s. auch Regressionsanalyse) bestimmt:

    Minimiere $\|Ax-b\|^{2}=\sum _{j=1}^{m}\left(\sum _{k=1}^{n}a_{jk}x_{k}-b_{j}\right)^{2}$.

In diesem Fall ist $A^{+}=R^{+}Q^{T}$ die Moore-Penrose-Pseudoinverse von $A$ und für die berechnete Kleinste-Quadrate-Lösung $x$ gilt die Beziehung $x=A^{+}b$, die die übliche Darstellung $x=A^{-1}b$ des regulären Falls $m=n$ verallgemeinert.

### Lösung unterbestimmter Gleichungssysteme

Für $m<n$ hat die Matrix $A$ einen nichttrivialen Kern. Bei vollem Rang von $A$ bilden die Lösungen des Gleichungssystems $Ax=b$ daher einen affinen Unterraum. Diejenige Lösung mit kleinster Norm liegt im orthogonalen Komplement des Kerns und man bekommt sie mit Hilfe einer QR-Zerlegung von $A^{T}$:

* Bestimme eine QR-Zerlegung der Matrix $A^{T}=QR$.
* Löse $R^{T}z=b\in \mathbb {R} ^{m}$ durch Vorwärtseinsetzen.
* Berechne $x=Qz\in \mathbb {R} ^{n}$.

Auch hier ist wieder $A^{+}=Q(R^{T})^{+}$ die Moore-Penrose-Pseudoinverse von $A$ und für die berechnete Lösung kleinster Norm gilt die Beziehung $x=A^{+}b$.

