Permutationsmatrix

Eine Permutationsmatrix oder auch Vertauschungsmatrix ist in der Mathematik eine Matrix, bei der in jeder Zeile und in jeder Spalte genau ein Eintrag eins ist und alle anderen Einträge null sind. Jede Permutationsmatrix entspricht genau einer Permutation einer endlichen Menge von Zahlen. Wird eine Permutationsmatrix mit einem Vektor multipliziert, dann werden die Komponenten des Vektors entsprechend dieser Permutation vertauscht. Permutationsmatrizen sind orthogonal, doppelt-stochastisch und ganzzahlig unimodular. Die Menge der Permutationsmatrizen fester Größe bildet mit der Matrizenmultiplikation eine zur symmetrische Gruppe isomorphe Untergruppe der allgemeinen linearen Gruppe. Permutationsmatrizen werden unter anderem in der linearen Algebra, der Kombinatorik und der Kryptographie verwendet.

Eine Permutationsmatrix ist eine quadratische Matrix, bei der genau ein Eintrag pro Zeile und Spalte gleich 1 {\displaystyle 1} ist und alle anderen Einträge gleich 0 {\displaystyle 0} sind.^[1] Hierbei sind im Allgemeinen 1 {\displaystyle 1} und 0 {\displaystyle 0} das Einselement und Nullelement eines zugrunde liegenden Rings R {\displaystyle R} . Jede Permutationsmatrix der Größe n × n {\displaystyle n\times n} entspricht genau einer Permutation ( π ( 1 ) , … , π ( n ) ) {\displaystyle (\pi (1),\ldots ,\pi (n))} der Zahlen von 1 {\displaystyle 1} bis n {\displaystyle n} . Die zu einer Permutation π ∈ S n {\displaystyle \pi \in S_{n}} zugehörige Permutationsmatrix

P π = ( p i j ) ∈ R n × n {\displaystyle P_{\pi }=(p_{ij})\in R^{n\times n}}

hat dann als Einträge

p i j = δ π ( i ) , j = { 1 , falls  π ( i ) = j 0 , sonst, {\displaystyle p_{ij}=\delta _{\pi (i),j}={\begin{cases}1,&{\text{falls }}\pi (i)=j\\0,&{\text{sonst,}}\end{cases}}}

wobei δ i j {\displaystyle \delta _{ij}} das Kronecker-Delta bezeichne. Werden durch die Permutation π {\displaystyle \pi } genau zwei Zahlen miteinander vertauscht, so bezeichnet man P π {\displaystyle P_{\pi }} auch als Vertauschungsmatrix. Ist e i {\displaystyle e_{i}} der i {\displaystyle i} -te kanonische Einheitsvektor als Zeilenvektor, dann lässt sich die Permutationsmatrix P π {\displaystyle P_{\pi }} auch durch

P π = ( e π ( 1 ) ⋮ e π ( n ) ) {\displaystyle P_{\pi }={\begin{pmatrix}e_{\pi (1)}\\\vdots \\e_{\pi (n)}\end{pmatrix}}}

darstellen. Gelegentlich findet sich allerdings in der Literatur auch die umgekehrte Variante, bei der die Einheitsvektoren spaltenweise zusammengesetzt werden, wodurch die Permutationsmatrizen entsprechend transponiert werden.^[2] Im Folgenden wird jedoch die gebräuchlichere erste Variante verwendet.

Die zu der Permutation

π = ( 1 2 3 4 5 4 2 1 5 3 ) ∈ S 5 {\displaystyle \pi ={\begin{pmatrix}1&2&3&4&5\\4&2&1&5&3\\\end{pmatrix}}\in S_{5}}

zugehörige Permutationsmatrix ist

P π = ( e 4 e 2 e 1 e 5 e 3 ) = ( 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 1 0 0 ) ∈ R 5 × 5 {\displaystyle P_{\pi }={\begin{pmatrix}e_{4}\\e_{2}\\e_{1}\\e_{5}\\e_{3}\end{pmatrix}}={\begin{pmatrix}0&0&0&1&0\\0&1&0&0&0\\1&0&0&0&0\\0&0&0&0&1\\0&0&1&0&0\end{pmatrix}}\in R^{5\times 5}} .

Nachdem durch die Permutation π {\displaystyle \pi } beispielsweise die Zahl 5 {\displaystyle 5} auf die Zahl 3 {\displaystyle 3} abgebildet wird, findet sich in der fünften Zeile von P π {\displaystyle P_{\pi }} die 1 {\displaystyle 1} in der dritten Spalte.

Wird eine Permutationsmatrix mit einem gegebenen Spaltenvektor v = ( v 1 , … , v n ) T {\displaystyle v=(v_{1},\ldots ,v_{n})^{T}} multipliziert, dann ergibt das Matrix-Vektor-Produkt

P π ⋅ v = ( e π ( 1 ) ⋮ e π ( n ) ) ⋅ ( v 1 ⋮ v n ) = ( v π ( 1 ) ⋮ v π ( n ) ) {\displaystyle P_{\pi }\cdot v={\begin{pmatrix}e_{\pi (1)}\\\vdots \\e_{\pi (n)}\end{pmatrix}}\cdot {\begin{pmatrix}v_{1}\\\vdots \\v_{n}\end{pmatrix}}={\begin{pmatrix}v_{\pi (1)}\\\vdots \\v_{\pi (n)}\end{pmatrix}}}

einen neuen Spaltenvektor, dessen Einträge entsprechend der Permutation π {\displaystyle \pi } vertauscht wurden. Ist beispielsweise v = ( v 1 , v 2 , v 3 , v 4 , v 5 ) T {\displaystyle v=(v_{1},v_{2},v_{3},v_{4},v_{5})^{T}} , dann ergibt das Matrix-Vektor-Produkt mit der obigen Beispiel-Permutationsmatrix den Spaltenvektor

P π ⋅ v = ( 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 1 0 0 ) ⋅ ( v 1 v 2 v 3 v 4 v 5 ) = ( v 4 v 2 v 1 v 5 v 3 ) . {\displaystyle P_{\pi }\cdot v={\begin{pmatrix}0&0&0&1&0\\0&1&0&0&0\\1&0&0&0&0\\0&0&0&0&1\\0&0&1&0&0\end{pmatrix}}\cdot {\begin{pmatrix}v_{1}\\v_{2}\\v_{3}\\v_{4}\\v_{5}\end{pmatrix}}={\begin{pmatrix}v_{4}\\v_{2}\\v_{1}\\v_{5}\\v_{3}\end{pmatrix}}.}

Wird eine Matrix von links mit einer Permutationsmatrix multipliziert, dann werden die Zeilen der Matrix gemäß der Permutation vertauscht. Umgekehrt ergibt die Multiplikation eines Zeilenvektors mit der transponierten Permutationsmatrix wieder einen Zeilenvektor mit entsprechend der Permutation π {\displaystyle \pi } vertauschten Elementen, also

v T ⋅ P π T = ( v 1 … v n ) ⋅ ( e π ( 1 ) T … e π ( n ) T ) = ( v π ( 1 ) … v π ( n ) ) {\displaystyle v^{T}\cdot P_{\pi }^{T}={\begin{pmatrix}v_{1}&\ldots &v_{n}\end{pmatrix}}\cdot {\begin{pmatrix}e_{\pi (1)}^{T}&\ldots &e_{\pi (n)}^{T}\end{pmatrix}}={\begin{pmatrix}v_{\pi (1)}&\ldots &v_{\pi (n)}\end{pmatrix}}} .

Für obiges Beispiel erhält man somit

v T ⋅ P π T = ( v 1 v 2 v 3 v 4 v 5 ) ⋅ ( 0 0 1 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 ) = ( v 4 v 2 v 1 v 5 v 3 ) . {\displaystyle v^{T}\cdot P_{\pi }^{T}={\begin{pmatrix}v_{1}&v_{2}&v_{3}&v_{4}&v_{5}\end{pmatrix}}\cdot {\begin{pmatrix}0&0&1&0&0\\0&1&0&0&0\\0&0&0&0&1\\1&0&0&0&0\\0&0&0&1&0\end{pmatrix}}={\begin{pmatrix}v_{4}&v_{2}&v_{1}&v_{5}&v_{3}\end{pmatrix}}.}

Wird eine Matrix von rechts mit der transponierten Permutationsmatrix multipliziert, werden entsprechend die Spalten der Matrix gemäß der Permutation vertauscht.

Permutationsmatrizen sind stets invertierbar, wobei die Inverse einer Permutationsmatrix gerade ihre Transponierte ist. Die transponierte Matrix ist dabei die Permutationsmatrix der inversen Permutation, es gilt also

P π − 1 = P π T = P π − 1 {\displaystyle P_{\pi }^{-1}=P_{\pi }^{T}=P_{\pi ^{-1}}} .

Reelle Permutationsmatrizen sind demnach stets orthogonal und haben vollen Rang n {\displaystyle n} .

Das Produkt zweier Permutationsmatrizen ist wieder eine Permutationsmatrix, die der Hintereinanderausführung der zugehörigen Permutationen entspricht. Die Permutationsmatrix der Hintereinanderausführung zweier Permutationen π , σ ∈ S n {\displaystyle \pi ,\sigma \in S_{n}} ergibt sich zu

P π ∘ σ = P σ ⋅ P π {\displaystyle P_{\pi \circ \sigma }=P_{\sigma }\cdot P_{\pi }} .

Die Abbildung π ↦ P π {\displaystyle \pi \mapsto P_{\pi }} stellt somit einen Antihomomorphismus dar. Die Menge der Permutationsmatrizen bildet zusammen mit der Matrizenmultiplikation eine Gruppe, und zwar eine Untergruppe der allgemeinen linearen Gruppe G L ( n , R ) {\displaystyle \mathrm {GL} (n,R)} . Jede Permutationsmatrix kann dabei als Produkt von elementaren zeilenvertauschenden Matrizen dargestellt werden.

Ganzzahlige Potenzen von Permutationsmatrizen sind wieder Permutationsmatrizen. Für jede Permutationsmatrix P π {\displaystyle P_{\pi }} gibt es dabei eine Potenz k {\displaystyle k} , sodass

P π k = I {\displaystyle P_{\pi }^{k}=I}

ergibt, wobei I {\displaystyle I} die Einheitsmatrix ist. Das kleinste positive k {\displaystyle k} mit dieser Eigenschaft ist gleich der Ordnung von P π {\displaystyle P_{\pi }} in der allgemeinen linearen Gruppe. Diese Ordnung ist gleich dem kleinsten gemeinsamen Vielfachen der Längen der disjunkten Zyklen von π {\displaystyle \pi } .

Die Determinante einer Permutationsmatrix ist entweder + 1 {\displaystyle +1} oder − 1 {\displaystyle -1} und entspricht dem Vorzeichen der zugehörigen Permutation:

det ( P π ) = sgn ⁡ ( π ) {\displaystyle \det(P_{\pi })=\operatorname {sgn} (\pi )} .

Eine Permutationsmatrix über den ganzen Zahlen ist damit eine ganzzahlige unimodulare Matrix. Die Spur einer ganzzahligen Permutationsmatrix entspricht der Anzahl der Fixpunkte der Permutation.

Die Determinante kann mit dem folgenden Schema ermittelt werden. Dazu wird für jede Spalte der Matrix die Zeilennummer, in der die eins steht, nebeneinander in eine Tabelle geschrieben. Darunter wird für jede Zahl z der ersten Zeile die Anzahl der Zahlen geschrieben, die größer als z sind und in der Tabelle links von z stehen; diese Anzahl heißt Kennmarke. Bei der eingangs angegebenen 8×8-Permutationsmatrix ergibt sich

Zeilennummer z	4	7	1	6	2	8	5	3
Kennmarke a	0	0	2	1	3	0	3	5

Ist die Summe der Kennmarken, wie hier, eine gerade Zahl, dann ist die Determinante 1, sonst −1; als Formel bei einer Permutationsmatrix der Größe n × n {\displaystyle n\times n} :^[3]

d e t ( P π ) = ( − 1 ) ν {\displaystyle \mathrm {det} (P_{\pi })=(-1)^{\nu }}   mit   ν = ∑ j = 1 n a j {\displaystyle \nu =\sum _{j=1}^{n}{\mathsf {a}}_{j}}

Die Eigenwerte einer reellen Permutationsmatrix sind nicht notwendigerweise alle reell, sie liegen aber auf dem komplexen Einheitskreis. Sind l 1 , … , l s {\displaystyle l_{1},\ldots ,l_{s}} die Längen der Zyklen einer Permutation π {\displaystyle \pi } , dann sind die Eigenwerte der zugehörigen Permutationsmatrix P π {\displaystyle P_{\pi }} die komplexen Einheitswurzeln

λ j k = e 2 π i k / l j {\displaystyle \lambda _{jk}=e^{2\pi ik/l_{j}}}

für j = 1 , … , s {\displaystyle j=1,\ldots ,s} und k = 1 , … , l j {\displaystyle k=1,\ldots ,l_{j}} . Eine reelle Permutationsmatrix besitzt demnach genau dann den Eigenwert e 2 π i k / m {\displaystyle e^{2\pi ik/m}} , wobei k {\displaystyle k} und m {\displaystyle m} teilerfremd seien, wenn die zugrunde liegende Permutation mindestens einen Zyklus aufweist, dessen Länge durch m {\displaystyle m} teilbar ist. Die Vielfachheit dieses Eigenwerts entspricht dann der Anzahl solcher Zyklen. Eine reelle Permutationsmatrix besitzt daher stets den Eigenwert 1 {\displaystyle 1} mit Vielfachheit gleich der Gesamtzahl der Zyklen s {\displaystyle s} der zugrunde liegenden Permutation.

Da reelle Permutationsmatrizen orthogonal sind, gilt für ihre Spektralnorm

‖ P π ‖ 2 = 1 {\displaystyle \|P_{\pi }\|_{2}=1} .

Für die Spalten- und Zeilensummennorm einer reellen Permutationsmatrix ergibt sich ebenfalls

‖ P π ‖ 1 = ‖ P π ‖ ∞ = 1 {\displaystyle \|P_{\pi }\|_{1}=\|P_{\pi }\|_{\infty }=1} .

Eine reelle Permutationsmatrix ist damit eine doppelt-stochastische Matrix. Nach dem Satz von Birkhoff und von Neumann ist eine quadratische Matrix genau dann doppelt-stochastisch, wenn sie eine Konvexkombination von Permutationsmatrizen ist.

• Die Permutationsmatrix der identischen Permutation ist die Einheitsmatrix.
• Eine Permutationsmatrix ist genau dann symmetrisch, wenn die zugrunde liegende Permutation selbstinvers ist.
• Weist eine Permutationsmatrix eine Blockstruktur auf, dann lässt sich die zugrunde liegende Permutation als Summe von Permutationen darstellen.

Permutationsmatrizen werden unter anderem verwendet:

• in der linearen Algebra als Elementarmatrizen bei der Gauß-Elimination zur Lösung linearer Gleichungssysteme
• in der Kombinatorik bei der Matrixdarstellung von Permutationsgruppen
• in der Kryptographie als Komponenten von Blockverschlüsselungsverfahren

In der Schachmathematik bilden die Permutationsmatrizen gerade die Lösungen des Problems, n {\displaystyle n} Türme auf ein Schachbrett der Größe n × n {\displaystyle n\times n} so zu verteilen, dass sich keine Türme gegenseitig angreifen. Schwieriger zu lösen ist das Damenproblem, bei dem die Türme durch Damen ersetzt werden, die auch diagonal angreifen können. Die Lösungen des Damenproblems sind ebenfalls Permutationsmatrizen.

Eine verallgemeinerte Permutationsmatrix oder monomiale Matrix ist eine quadratische Matrix G ∈ R n × n {\displaystyle G\in R^{n\times n}} , bei der genau ein Eintrag pro Zeile und Spalte ungleich 0 {\displaystyle 0} ist. Monomiale Matrizen haben die Darstellung

G = P ⋅ D {\displaystyle G=P\cdot D} ,

wobei P ∈ R n × n {\displaystyle P\in R^{n\times n}} eine gewöhnliche Permutationsmatrix und D ∈ R n × n {\displaystyle D\in R^{n\times n}} eine Diagonalmatrix ist, deren Diagonaleinträge alle ungleich 0 {\displaystyle 0} sind. Die regulären monomialen Matrizen bilden mit der Matrizenmultiplikation als Verknüpfung die monomiale Gruppe M ⁡ ( n , R ) {\displaystyle \operatorname {M} (n,R)} , eine weitere Untergruppe der allgemeinen linearen Gruppe GL ⁡ ( n , R ) {\displaystyle \operatorname {GL} (n,R)} . Spezielle monomiale Matrizen sind vorzeichenbehaftete Permutationsmatrizen, bei denen in jeder Zeile und jeder Spalte genau ein Eintrag + 1 {\displaystyle +1} oder − 1 {\displaystyle -1} ist und alle übrigen Einträge 0 {\displaystyle 0} sind.

• Levi-Civita-Symbol
• Rothe-Diagramm
• Irreduzible Matrix
• Zyklische Matrix

• Siegfried Bosch: Lineare Algebra. Springer, 2006, ISBN 3-540-29884-3. 
• Jörg Liesen, Volker Mehrmann: Lineare Algebra. 3. Auflage. Springer, Berlin, Heidelberg 2021, ISBN 978-3-662-62741-9, doi:10.1007/978-3-662-62742-6. 

1. ↑ Jörg Liesen, Volker Mehrmann: Lineare Algebra. Springer, 2011, S. 45. 
2. ↑ Siegfried Bosch: Lineare Algebra. Springer, 2006, S. 275. 
3. ↑ R. Zurmühl, S. Falk: Matrizen und ihre Anwendungen 1. Grundlagen, Für Ingenieure, Physiker und Angewandte Mathematiker. Springer, Berlin u. a. 1997, ISBN 3-540-61436-2, S. 57. 

• Eric W. Weisstein: Permutation Matrix. In: MathWorld (englisch).
• Wkbj79: Permutation Matrix. In: PlanetMath. (englisch)