Orthogonale Matrix

Eine orthogonale Matrix ist in der linearen Algebra eine quadratische, reelle Matrix, deren Zeilen- und Spaltenvektoren paarweise orthonormal zueinander sind. Damit ist die Inverse einer orthogonalen Matrix gleichzeitig ihre Transponierte. Die Menge der orthogonalen Matrizen fester Größe bildet mit der Matrizenmultiplikation als Verknüpfung die orthogonale Gruppe. Orthogonale Matrizen stellen Kongruenzabbildungen im euklidischen Raum, also Drehungen, Spiegelungen und Drehspiegelungen, dar. Sie werden beispielsweise bei der numerischen Lösung linearer Gleichungssysteme oder Eigenwertprobleme eingesetzt. Der analoge Begriff bei komplexen Matrizen ist die unitäre Matrix.

Definition

Eine reelle quadratische Matrix $Q\in \mathbb {R} ^{n\times n}$ heißt orthogonal, wenn das Produkt mit ihrer transponierten Matrix $Q^{T}$ die Einheitsmatrix $I$ ergibt, also

$Q\cdot Q^{T}=I$

gilt. Werden die Zeilenvektoren der Matrix $Q$ mit $q_{1},\ldots ,q_{n}$ bezeichnet, dann ist diese Bedingung gleichbedeutend damit, dass das Standardskalarprodukt zweier Zeilenvektoren

$q_{i}\cdot q_{j}^{T}=\delta _{ij}={\begin{cases}1&{\text{falls}}~i=j\\0&{\text{sonst}}\end{cases}}$

ergibt, wobei $\delta _{ij}$ das Kronecker-Delta ist. Die Zeilenvektoren einer orthogonalen Matrix bilden damit eine Orthonormalbasis des Koordinatenraums $\mathbb {R} ^{n}$. Dies trifft auch für die Spaltenvektoren einer orthogonalen Matrix zu, denn $Q^{T}$ ist aufgrund von

$Q^{T}\cdot Q=(Q\cdot Q^{T})^{T}=I^{T}=I$

ebenfalls orthogonal. Auch wenn die Bezeichnung „orthogonale Matrix“ so verstanden werden könnte, reicht es nicht aus, wenn die Zeilen- oder Spaltenvektoren lediglich paarweise orthogonal sind; sie müssen zusätzlich normiert sein, also die Länge eins aufweisen.

Beispiele

Konkrete Beispiele

• Die Matrix

$Q={\begin{pmatrix}0&1\\1&0\end{pmatrix}}$

ist orthogonal, denn es gilt

$Q\,Q^{T}={\begin{pmatrix}0&1\\1&0\end{pmatrix}}\cdot {\begin{pmatrix}0&1\\1&0\end{pmatrix}}={\begin{pmatrix}0\cdot 0+1\cdot 1&0\cdot 1+1\cdot 0\\1\cdot 0+0\cdot 1&1\cdot 1+0\cdot 0\end{pmatrix}}={\begin{pmatrix}1&0\\0&1\end{pmatrix}}=I$.

• Auch die Matrix

$Q={\frac {1}{5}}{\begin{pmatrix}3&-4\\4&3\end{pmatrix}}$

ist orthogonal, denn es gilt

$Q\,Q^{T}={\frac {1}{5}}{\begin{pmatrix}3&-4\\4&3\end{pmatrix}}\cdot {\frac {1}{5}}{\begin{pmatrix}3&4\\-4&3\end{pmatrix}}={\frac {1}{25}}{\begin{pmatrix}9+16&12-12\\12-12&16+9\end{pmatrix}}={\frac {1}{25}}{\begin{pmatrix}25&0\\0&25\end{pmatrix}}=I$.

Allgemeine Beispiele

• Permutationsmatrizen, also Matrizen, bei denen genau ein Eintrag pro Zeile und Spalte gleich $1$ ist und alle anderen Einträge gleich $0$ sind, sind orthogonal. Bezeichnet $P_{\pi }$ die zu einer Permutation $\pi$ zugehörige Permutationsmatrix, dann gilt

$P_{\pi }\,P_{\pi }^{T}=P_{\pi }\,P_{\pi ^{-1}}=P_{\pi \circ \pi ^{-1}}=P_{\mathrm {id} }=I$,

denn die transponierte Permutationsmatrix ist gleich der Permutationsmatrix der inversen Permutation, die alle Vertauschungen rückgängig macht, und das Produkt von Permutationsmatrizen entspricht der Hintereinanderausführung der Permutationen.

• Drehmatrizen, also Matrizen, die eine Drehung um den Koordinatenursprung in der euklidischen Ebene beschreiben, sind orthogonal. Bezeichnet

$R_{\alpha }={\begin{pmatrix}\cos \alpha &-\sin \alpha \\\sin \alpha &\cos \alpha \end{pmatrix}}$

die Drehmatrix einer Drehung um einen Winkel $\alpha$, die den Ursprung festlässt, dann gilt mit dem trigonometrischen Pythagoras

$R_{\alpha }\,R_{\alpha }^{T}={\begin{pmatrix}\cos ^{2}\alpha +\sin ^{2}\alpha &\cos \alpha \sin \alpha -\sin \alpha \cos \alpha \\\sin \alpha \cos \alpha -\cos \alpha \sin \alpha &\sin ^{2}\alpha +\cos ^{2}\alpha \end{pmatrix}}=I$.

Allgemeiner sind auch Drehmatrizen, die eine Drehung in einer beliebigen Ursprungsebene im $n$-dimensionalen Raum beschreiben, orthogonal.

• Spiegelungsmatrizen, also Matrizen, die eine (senkrechte) Spiegelung an einer Ursprungsgerade in der euklidischen Ebene beschreiben, sind orthogonal. Bezeichnet

$S_{n}=I-2nn^{T}$

die Spiegelungsmatrix einer Spiegelung an einer Ursprungsgerade mit Einheits-Normalenvektor $n$, dann gilt

$S_{n}\,S_{n}^{T}=S_{n}\,S_{n}=(I-2nn^{T})\,(I-2nn^{T})=I-4nn^{T}+4n(n^{T}n)n^{T}=I$,

denn Spiegelungsmatrizen sind nach Definition symmetrisch und für einen Einheitsvektor gilt $n^{T}n=1$. Allgemeiner sind auch Spiegelungsmatrizen an einer beliebigen linearen Mannigfaltigkeit im $n$-dimensionalen Raum, beispielsweise einer Hyperebene, die durch den Ursprung verläuft, orthogonal.

Eigenschaften

Inverse

Eine orthogonale Matrix $Q\in \mathbb {R} ^{n\times n}$ ist immer invertierbar und ihre Inverse ist gleich ihrer Transponierten, das heißt

$Q^{T}=Q^{-1}$.

Die Inverse einer Matrix ist nämlich gerade diejenige Matrix, die mit der Ausgangsmatrix multipliziert die Einheitsmatrix ergibt. Es gilt auch die Umkehrung und jede Matrix $Q$, deren Transponierte gleich ihrer Inversen ist, ist orthogonal, denn es gilt dann

$Q\,Q^{T}=Q\,Q^{-1}=I$.

Längen- und Winkeltreue

Wird eine orthogonale Matrix $Q\in \mathbb {R} ^{n\times n}$ mit einem Vektor $x\in \mathbb {R} ^{n}$ multipliziert, ändert sich die Länge (euklidische Norm) des Vektors nicht, das heißt

$\|Q\,x\|_{2}=\|x\|_{2}$.

Weiter ist Standardskalarprodukt zweier Vektoren $x,y\in \mathbb {R} ^{n}$ invariant bezüglich der Multiplikation mit einer orthogonalen Matrix $Q$, also

$\left\langle Q\,x,Q\,y\right\rangle =\left\langle x,y\right\rangle$.

Beide Eigenschaften folgen dabei direkt aus der Verschiebungseigenschaft des Standardskalarprodukts. Aufgrund dieser Längen- und Winkeltreue stellt die lineare Abbildung

$f\colon \mathbb {R} ^{n}\to \mathbb {R} ^{n},x\mapsto Q\,x$

eine Kongruenzabbildung im euklidischen Raum dar. Umgekehrt ist die Abbildungsmatrix jeder längen- oder winkeltreuen linearen Abbildung zwischen zwei endlichdimensionalen reellen Vektorräumen orthogonal, denn es gilt

$\left\langle q_{i},q_{j}\right\rangle =\left\langle Q\,e_{i},Q\,e_{j}\right\rangle =\left\langle e_{i},e_{j}\right\rangle =\delta _{ij}$,

wobei $e_{i}$ der $i$-te Standardbasisvektor ist.

Determinante

Für den Betrag der Determinante einer orthogonalen Matrix $Q\in \mathbb {R} ^{n\times n}$ gilt

$|\det Q|=1$,

was mit Hilfe des Determinantenproduktsatzes über

$\operatorname {det} ^{2}Q=\det Q\cdot \det Q=\det Q\cdot \det Q^{T}=\det(QQ^{T})=\det I=1$

folgt. Damit kann die Determinante einer orthogonalen Matrix nur die Werte eins oder minus eins annehmen. Es gibt allerdings auch nicht orthogonale Matrizen, deren Determinante plus oder minus eins ist, zum Beispiel unimodulare Matrizen. Orthogonale Matrizen, deren Determinante eins ist, entsprechen Drehungen, während orthogonale Matrizen, deren Determinante minus eins ist, Drehspiegelungen darstellen.

Eigenwerte

Die Eigenwerte einer orthogonalen Matrix $Q\in \mathbb {R} ^{n\times n}$ sind nicht notwendigerweise alle reell. Sie haben jedoch den komplexen Betrag eins, sind also von der Form

$\lambda =e^{it}$

mit $t\in \mathbb {R}$. Ist nämlich $x$ ein zu $\lambda$ gehöriger Eigenvektor, dann gilt aufgrund der Längentreue und der absoluten Homogenität einer Norm

$\|x\|_{2}=\|Q\,x\|_{2}=\|\lambda \,x\|_{2}=|\lambda |\,\|x\|_{2}$

und daher $|\lambda |=1$. Eine orthogonale Matrix besitzt demnach höchstens die reellen Eigenwerte $\pm 1$. Die komplexen Eigenwerte treten immer paarweise komplex konjugiert auf, das heißt mit $\lambda =e^{it}$ ist auch ${\bar {\lambda }}=e^{-it}$ ein Eigenwert.

Diagonalisierbarkeit

Eine orthogonale Matrix $Q\in \mathbb {R} ^{n\times n}$ ist normal, das heißt es gilt

$Q\,Q^{T}=Q^{T}\,Q$,

und damit über den komplexen Zahlen unitär diagonalisierbar. Nach dem Spektralsatz gibt es nämlich eine unitäre Matrix $U\in \mathbb {C} ^{n\times n}$, sodass

$U^{-1}\,Q\,U=D$

gilt, wobei $D\in \mathbb {C} ^{n\times n}$ eine Diagonalmatrix mit den Eigenwerten von $Q$ ist. Im Allgemeinen ist eine orthogonale Matrix $Q\in \mathbb {R} ^{n\times n}$ jedoch nicht reell diagonalisierbar. Es existiert allerdings eine orthogonale Matrix $V\in \mathbb {R} ^{n\times n}$, sodass

$V^{-1}\,Q\,V={\begin{pmatrix}D_{1}&&0\\&\ddots &\\0&&D_{s}\end{pmatrix}}$

eine Blockdiagonalmatrix ergibt, bei der die einzelnen Blöcke entweder Drehmatrizen der Größe $2\times 2$ sind oder aus der Zahl $+1$ oder $-1$ bestehen.

Normen

Die Spektralnorm einer orthogonalen Matrix $Q\in \mathbb {R} ^{n\times n}$ ist

$\|Q\|_{2}=\max _{\|x\|_{2}=1}\|Q\,x\|_{2}=\max _{\|x\|_{2}=1}\|x\|_{2}=1$.

Für die Frobeniusnorm gilt mit dem Frobenius-Skalarprodukt entsprechend

$\|Q\|_{F}={\sqrt {\langle Q,Q\rangle }}={\sqrt {\langle I,I\rangle }}={\sqrt {n}}$.

Das Produkt mit einer orthogonalen Matrix erhält sowohl die Spektralnorm, als auch die Frobeniusnorm einer gegebenen Matrix $A\in \mathbb {R} ^{n\times n}$, denn es gilt

$\|Q\,A\|_{2}=\max _{\|x\|_{2}=1}\|Q\,A\,x\|_{2}=\max _{\|x\|_{2}=1}\|A\,x\|_{2}=\|A\|_{2}$

und

$\|Q\,A\|_{F}={\sqrt {\langle Q\,A,Q\,A\rangle }}={\sqrt {\langle A,A\rangle }}=\|A\|_{F}$.

Damit bleibt auch die Kondition einer Matrix bezüglich dieser Normen nach Multiplikation mit einer orthogonalen Matrix erhalten.

Orthogonale Matrizen als Gruppe

→ Hauptartikel: Orthogonale Gruppe

Das Produkt zweier orthogonaler Matrizen $P,Q\in \mathbb {R} ^{n\times n}$ ist wieder orthogonal, denn es gilt

$(P\,Q)\,(P\,Q)^{T}=P\,(Q\,Q^{T})\,P^{T}=P\,P^{T}=I$.

Weiter ist die Inverse einer orthogonalen Matrix $Q\in \mathbb {R} ^{n\times n}$ ebenfalls orthogonal, denn es gilt

$Q^{-1}\,Q^{-T}=Q^{T}\,Q^{-T}=(Q^{-1}\,Q)^{T}=I^{T}=I$.

Die orthogonalen Matrizen fester Größe bilden daher mit der Matrizenmultiplikation als Verknüpfung eine Gruppe, die orthogonale Gruppe $\mathrm {O} (n)$. Dabei handelt es sich um eine Untergruppe der allgemeinen linearen Gruppe $\mathrm {GL} (n,\mathbb {R} )$. Die orthogonalen Matrizen mit Determinante eins, also die Drehmatrizen, bilden wiederum eine Untergruppe der orthogonalen Gruppe, die Drehgruppe (oder spezielle orthogonale Gruppe) $\mathrm {SO} (n)$. Die orthogonalen Matrizen mit Determinante minus eins, also die Drehspiegelungen, bilden keine Untergruppe der orthogonalen Gruppe, denn ihnen fehlt ein neutrales Element, sondern lediglich eine Nebenklasse.

Verwendung

Lineare Gleichungssysteme

Die Lösung linearer Gleichungssysteme der Form

$Q\,x=b$

mit einer orthogonalen Matrix $Q\in \mathbb {R} ^{n\times n}$ und einer rechten Seite $b\in \mathbb {R} ^{n}$ lässt sich numerisch effizient durch

$x=Q^{T}\,b$

berechnen. Die Ermittlung der Lösung $x\in \mathbb {R} ^{n}$ erfordert also lediglich eine Matrix-Vektor-Multiplikation, die mit einem Aufwand der Ordnung $O(n^{2})$ durchgeführt werden kann. Im Vergleich dazu benötigt die Lösung allgemeiner linearer Gleichungssysteme beispielsweise mit Hilfe der Gauß-Elimination einen Aufwand $O(n^{3})$. Dieser Vorteil wird beispielsweise bei der (reellen) diskreten Fourier-Transformation und der diskreten Kosinus-Transformation genutzt.

Matrixzerlegungen

Eine weitere Anwendung orthogonaler Matrizen ist die QR-Zerlegung einer gegebenen Matrix $A\in \mathbb {R} ^{m\times n}$ als Produkt

$A=Q\,R$

einer orthogonalen Matrix $Q\in \mathbb {R} ^{m\times m}$ und einer oberen Dreiecksmatrix $R\in \mathbb {R} ^{m\times n}$. Die Konstruktion der Matrix $Q$ kann dabei mit Givens-Rotationen, die Drehungen entsprechen, oder Householdertransformationen, die Spiegelungen entsprechen, durchgeführt werden. QR-Zerlegungen werden in der Numerik bei der Lösung schlecht konditionierter, überbestimmter oder unterbestimmter linearer Gleichungssysteme eingesetzt. Ein weiteres Anwendungsfeld besteht in der Berechnung von Eigenwertproblemen mit dem QR-Algorithmus.

Mit Hilfe der Singulärwertzerlegung lässt sich jede reelle Matrix $A\in \mathbb {R} ^{m\times n}$ auch als Produkt

$A=U\,\Sigma \,V^{T}$

einer orthogonalen Matrix $U\in \mathbb {R} ^{m\times m}$, einer Diagonalmatrix $\Sigma \in \mathbb {R} ^{m\times n}$ und der Transponierten einer weiteren orthogonalen Matrix $V\in \mathbb {R} ^{n\times n}$ darstellen. Die Diagonaleinträge der Matrix $\Sigma$ sind dann die Singulärwerte von $A$. Die Singulärwertzerlegung wird beispielsweise in der Geometrie bei der Hauptachsentransformation von Quadriken und in der Statistik bei der Hauptkomponentenanalyse multivariater Datensätze eingesetzt.

Siehe auch

• Orthogonale Abbildung
• Orthogonalisierungsverfahren

Literatur

• Gerd Fischer: Lineare Algebra. (Eine Einführung für Studienanfänger). 14. durchgesehene Auflage. Vieweg, 2003, ISBN 3-528-03217-0. 
• Jörg Liesen, Volker Mehrmann: Lineare Algebra. Springer, 2011, ISBN 978-3-8348-8290-5. 
• Hans Rudolf Schwarz, Norbert Köckler: Numerische Mathematik. Springer, 2009, ISBN 978-3-8348-0683-3. 
• Eberhard Zeidler, Wolfgang Hackbusch (Hrsg.): Taschenbuch der Mathematik. Band 1. Springer, 2012, ISBN 978-3-8351-0123-4. 

Weblinks

• D.A. Suprunenko: Orthogonal matrix. In: Michiel Hazewinkel (Hrsg.): Encyclopedia of Mathematics. Springer-Verlag und EMS Press, Berlin 2002, ISBN 1-55608-010-7 (englisch, encyclopediaofmath.org). 
• Todd Rowland: Orthogonal Matrix. In: MathWorld (englisch).
• akrowne: Orthogonal matrices. In: PlanetMath. (englisch)