Matriz permutación

La matriz permutación es la matriz cuadrada con todos sus n×n elementos iguales a 0, excepto uno cualquiera por cada fila y columna, el cual debe ser igual a 1. De acuerdo a esta definición existen n! matrices de permutación distintas, de las cuales una mitad corresponde a matrices de permutación par (con el determinante igual a 1) y la otra mitad a matrices de permutación impar (con el determinante igual a -1).

Para n = 3 se tiene:

Matrices de permutación par:

( 1 0 0 0 1 0 0 0 1 ) ( 0 0 1 1 0 0 0 1 0 ) ( 0 1 0 0 0 1 1 0 0 ) {\displaystyle {\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\\\end{pmatrix}}{\begin{pmatrix}0&0&1\\1&0&0\\0&1&0\\\end{pmatrix}}{\begin{pmatrix}0&1&0\\0&0&1\\1&0&0\\\end{pmatrix}}}

Matrices de permutación impar:

( 0 0 1 0 1 0 1 0 0 ) ( 0 1 0 1 0 0 0 0 1 ) ( 1 0 0 0 0 1 0 1 0 ) {\displaystyle {\begin{pmatrix}0&0&1\\0&1&0\\1&0&0\\\end{pmatrix}}{\begin{pmatrix}0&1&0\\1&0&0\\0&0&1\\\end{pmatrix}}{\begin{pmatrix}1&0&0\\0&0&1\\0&1&0\\\end{pmatrix}}}

Puede notarse que las matrices de permutación conforman un grupo de orden n! respecto al producto.

Definición

Dada una permutación π {\displaystyle \pi } de m elementos,

π : { 1 , , m } { 1 , , m } {\displaystyle \pi :\lbrace 1,\ldots ,m\rbrace \to \lbrace 1,\ldots ,m\rbrace }

representada en forma de dos líneas por

( 1 2 m π ( 1 ) π ( 2 ) π ( m ) ) , {\displaystyle {\begin{pmatrix}1&2&\cdots &m\\\pi (1)&\pi (2)&\cdots &\pi (m)\end{pmatrix}},}

hay dos maneras naturales de asociar la permutación con una matriz de permutación; a es decir, comenzando con la matriz identidad de m × m, I m {\displaystyle I_{m}} , permuta las columnas o permuta las filas, según π {\displaystyle \pi } . Ambos métodos para definir las matrices de permutación aparecen en la literatura y las propiedades expresadas en una representación pueden ser fácilmente convertidas a la otra representación. Este artículo tratará principalmente de una sola de estas representaciones y la otra sólo se mencionará cuando haya una diferencia que haya que tener en cuenta.

La matriz permutación de m × m P π {\displaystyle P_{\pi }} = (pij) que se obtiene al permutar las columnas de la matriz identidad I m {\displaystyle I_{m}} es decir, por cada i, pij = 1 si j = π {\displaystyle \pi } (i) y 0 en caso contrario, se referirá a la representación de la columna en este artículo. Dado que las entradas en la fila i son todas 0 excepto que un 1 aparece en la columna π {\displaystyle \pi } (i), podemos escribir

P π = [ e π ( 1 ) e π ( 2 ) e π ( m ) ] , {\displaystyle P_{\pi }={\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\vdots \\\mathbf {e} _{\pi (m)}\end{bmatrix}},}

donde e j {\displaystyle \mathbf {e} _{j}} , un vector de base estándar, denota un vector de longitud m con un 1 en la posición j y 0 en cualquier otra posición.

Por ejemplo, la matriz permutación P π {\displaystyle P_{\pi }} correspondiente a la permutación: π = ( 1 2 3 4 5 1 4 2 5 3 ) , {\displaystyle \pi ={\begin{pmatrix}1&2&3&4&5\\1&4&2&5&3\end{pmatrix}},} es

P π = [ e π ( 1 ) e π ( 2 ) e π ( 3 ) e π ( 4 ) e π ( 5 ) ] = [ e 1 e 4 e 2 e 5 e 3 ] = [ 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 ] . {\displaystyle P_{\pi }={\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\mathbf {e} _{\pi (3)}\\\mathbf {e} _{\pi (4)}\\\mathbf {e} _{\pi (5)}\end{bmatrix}}={\begin{bmatrix}\mathbf {e} _{1}\\\mathbf {e} _{4}\\\mathbf {e} _{2}\\\mathbf {e} _{5}\\\mathbf {e} _{3}\end{bmatrix}}={\begin{bmatrix}1&0&0&0&0\\0&0&0&1&0\\0&1&0&0&0\\0&0&0&0&1\\0&0&1&0&0\end{bmatrix}}.}

Observemos que la columna j de la matriz identidad I 5 {\displaystyle I_{5}} ahora aparece como la columna π ( j ) {\displaystyle \pi (j)} de P π {\displaystyle P_{\pi }}

La otra representación, obtenida por permutar las filas de la matriz identidad I m {\displaystyle I_{m}} , es decir, por cada j, pij = 1 si i = π ( j ) {\displaystyle \pi (j)} y 0 en caso contrario, se denominará representación de la fila.

Propiedades

La representación de las columnas de una matriz permutación se utiliza a lo largo de esta sección, salvo que se indique lo contrario.

Multiplicando P π {\displaystyle P_{\pi }} veces un vector columna g permutará las filas del vector:

P π g = [ e π ( 1 ) e π ( 2 ) e π ( n ) ] [ g 1 g 2 g n ] = [ g π ( 1 ) g π ( 2 ) g π ( n ) ] . {\displaystyle P_{\pi }\mathbf {g} ={\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\vdots \\\mathbf {e} _{\pi (n)}\end{bmatrix}}{\begin{bmatrix}g_{1}\\g_{2}\\\vdots \\g_{n}\end{bmatrix}}={\begin{bmatrix}g_{\pi (1)}\\g_{\pi (2)}\\\vdots \\g_{\pi (n)}\end{bmatrix}}.}

El uso repetido de este resultado muestra que si M {\displaystyle M} es una matriz de tamaño apropiado, el producto, P π {\displaystyle P_{\pi }} M {\displaystyle M} es sólo una permutación de las filas de M {\displaystyle M} . Sin embargo, observando que

P π e k T = e π 1 ( k ) T {\displaystyle P_{\pi }\mathbf {e} _{k}^{\mathsf {T}}=\mathbf {e} _{\pi ^{-1}(k)}^{\mathsf {T}}}

para cada k {\displaystyle k} se muestra que la permutación de las filas viene dada por π 1 {\displaystyle \pi ^{-1}} . ( M T {\displaystyle M^{T}} es la transpuesta de M {\displaystyle M} ).

Como las matrices de permutación son matrices ortogonales (es decir, P π P π T = I {\displaystyle P_{\pi }P_{\pi }^{\mathsf {T}}=I} ), existe la matriz inversa y se puede escribir como

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

Multiplicar un vector fila h veces P π {\displaystyle P_{\pi }} permutará las columnas del vector:

h P π = [ h 1 h 2 h n ] [ e π ( 1 ) e π ( 2 ) e π ( n ) ] = [ h π 1 ( 1 ) h π 1 ( 2 ) h π 1 ( n ) ] {\displaystyle \mathbf {h} P_{\pi }={\begin{bmatrix}h_{1}\;h_{2}\;\dots \;h_{n}\end{bmatrix}}{\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\vdots \\\mathbf {e} _{\pi (n)}\end{bmatrix}}={\begin{bmatrix}h_{\pi ^{-1}(1)}\;h_{\pi ^{-1}(2)}\;\dots \;h_{\pi ^{-1}(n)}\end{bmatrix}}}

Una vez más, la aplicación repetida de este resultado muestra que la post-multiplicación de una matriz M {\displaystyle M} por la matriz permutación P π {\displaystyle P_{\pi }} o sea, M {\displaystyle M} P π {\displaystyle P_{\pi }} resulta en la permutación de las columnas de M {\displaystyle M} . Nótese también que

e k P π = e π ( k ) . {\displaystyle \mathbf {e} _{k}P_{\pi }=\mathbf {e} _{\pi (k)}.}

Dadas dos permutaciones π {\displaystyle \pi } y σ de m elementos, las correspondientes matrices de permutación P π {\displaystyle P_{\pi }} y P σ {\displaystyle P_{\sigma }} actúan sobre vectores columna y están compuestos por

P σ P π g = P π σ g . {\displaystyle P_{\sigma }P_{\pi }\,\mathbf {g} =P_{\pi \,\circ \,\sigma }\,\mathbf {g} .}

Las mismas matrices que actúan sobre los vectores fila (es decir, post-multiplicación) se componen según la misma regla

h P σ P π = h P π σ . {\displaystyle \mathbf {h} P_{\sigma }P_{\pi }=\mathbf {h} P_{\pi \,\circ \,\sigma }.}

Para ser claros, las fórmulas anteriores utilizan la notación polaca para la composición de permutación, es decir,

π σ ( k ) = π ( σ ( k ) ) . {\displaystyle \pi \,\circ \,\sigma (k)=\pi \left(\sigma (k)\right).}

Dejemos que Q π {\displaystyle Q_{\pi }} sea la matriz permutación correspondiente a π {\displaystyle \pi } en su representación en filas. Las propiedades de esta representación se pueden determinar a partir de las de la representación de la columna, ya que Q π = P π T = P π 1 . {\displaystyle Q_{\pi }=P_{\pi }^{\mathsf {T}}=P_{{\pi }^{-1}}.} En particular,

Q π e k T = P π 1 e k T = e ( π 1 ) 1 ( k ) T = e π ( k ) T . {\displaystyle Q_{\pi }\mathbf {e} _{k}^{\mathsf {T}}=P_{{\pi }^{-1}}\mathbf {e} _{k}^{\mathsf {T}}=\mathbf {e} _{(\pi ^{-1})^{-1}(k)}^{\mathsf {T}}=\mathbf {e} _{\pi (k)}^{\mathsf {T}}.}

De ello se deduce que

Q σ Q π g = Q σ π g . {\displaystyle Q_{\sigma }Q_{\pi }\,\mathbf {g} =Q_{\sigma \,\circ \,\pi }\,\mathbf {g} .}

Del mismo modo,

h Q σ Q π = h Q σ π . {\displaystyle \mathbf {h} \,Q_{\sigma }Q_{\pi }=\mathbf {h} \,Q_{\sigma \,\circ \,\pi }.}


Otras propiedades:

  • El elemento neutro del grupo es la matriz identidad.
  • El elemento inverso de cada elemento del grupo de matrices de permutación es la matriz traspuesta correspondiente.
  • Cada elemento del grupo de matrices de permutación es una matriz ortogonal.
  • El producto de matrices de permutación par siempre genera una matriz de permutación par.
  • El producto de matrices de permutación impar siempre genera una matriz de permutación par.
  • El producto de matrices de permutación de paridad distinta siempre genera una matriz de permutación impar.
  • Observe que las matrices de permutación par conforman un semigrupo y que además el grupo de matrices de permutación no es conmutativo.
  • Cada elemento del grupo de matrices de permutación fuera del semigrupo es una matriz simétrica.

Grupo matricial

Si (1) denota la permutación identidad, entonces P 1 {\displaystyle P_{1}} es la matriz identidad.

Dejemos que S n {\displaystyle S_{n}} denote denote el grupo simétrico, o grupo de permutaciones, en {1,2,...,n}. Puesto que hay n permutaciones, hay n! matrices permutación. Por las fórmulas anteriores, las matrices permutación de n × n forman un grupo bajo la matriz multiplicación con la matriz identidad como elemento neutro.

El mapa SnA ⊂ GL(n, Z2) es una representación fiel. Por lo tanto, |A| = n!

Matrices dobles estocásticas

Una matriz permutación es en sí misma una matriz doble estocástica, pero también desempeña un papel especial en la teoría de dichas matrices. El teorema de Birkhoff-von Neumann dice que cada matriz real doble estocástica es una combinación convexa de matrices permutación del mismo orden y que las matrices permutación son precisamente los puntos extremos del conjunto de matrices dobles estocásticas. Es decir, el politopo de Birkhoff, el conjunto de matrices dobles estocásticas, es la envolvente convexa del conjunto de matrices permutación.

Propiedades algebraicas lineales

La traza de una matriz permutación es el número de puntos fijos de la permutación. Si la permutación tiene puntos fijos, se puede escribir en forma de ciclo como π = (a1)(a2)...(ak donde σ no tiene puntos fijos, Entonces ea1,ea2,...,eak son vectores propios de la matriz permutación.

Para calcular los vectores propios de una matriz permutación P σ {\displaystyle P_{\sigma }} anotamos σ {\displaystyle \sigma } como un producto de permutaciónes cíclicas, es decir, σ = C 1 C 2 C t {\displaystyle \sigma =C_{1}C_{2}\cdots C_{t}} . Dejemos que las longitudes correspondientes de estos ciclos sean l 1 , l 2 . . . l t {\displaystyle l_{1},l_{2}...l_{t}} , y dejemos que R i ( 1 i t ) {\displaystyle R_{i}(1\leq i\leq t)} sea el conjunto de soluciones complejas de x l i = 1 {\displaystyle x^{l_{i}}=1} . La unión de todos los R i {\displaystyle R_{i}} forman el conjunto de autovalores de la matriz permutación correspondiente. La multiplicidad geométrica de cada autovalor es igual al número de R i {\displaystyle R_{i}} que lo contengan.

De la teoría de grupos sabemos que cualquier permutación puede ser expresada como producto de transposiciones. Por lo tanto, cualquier matriz permutación P {\displaystyle P} los factores como producto de un intercambio de matrices elementales, cada una de las cuales tiene como determinante -1. Por lo tanto, el determinante de una matriz permutación P {\displaystyle P} es sólo la paridad de la permutación correspondiente.

Ejemplos

Permutación de filas y columnas

Cuando una matriz permutación P es multiplicada por la izquierda por una matriz M para generar PM permutará las filas de M (aquí los elementos de un vector columna),

cuando P es multiplicada por la derecha por M para generar MP permutará las columnas de M (aquí los elementos de un vector fila):

P * (1,2,3,4)T = (4,1,3,2)T



(1,2,3,4) * P = (2,4,3,1)











Las permutaciones de filas y columnas son por ejemplo reflexiones (ver más abajo) y permutaciones cícicas (ver matriz circulante).

Permutación de filas

La matriz permutación Pπ corresponde a la permutación: π = ( 1 2 3 4 5 1 4 2 5 3 ) , {\displaystyle \pi ={\begin{pmatrix}1&2&3&4&5\\1&4&2&5&3\end{pmatrix}},} hay

P π = [ e π ( 1 ) e π ( 2 ) e π ( 3 ) e π ( 4 ) e π ( 5 ) ] = [ e 1 e 4 e 2 e 5 e 3 ] = [ 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 ] . {\displaystyle P_{\pi }={\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\mathbf {e} _{\pi (3)}\\\mathbf {e} _{\pi (4)}\\\mathbf {e} _{\pi (5)}\end{bmatrix}}={\begin{bmatrix}\mathbf {e} _{1}\\\mathbf {e} _{4}\\\mathbf {e} _{2}\\\mathbf {e} _{5}\\\mathbf {e} _{3}\end{bmatrix}}={\begin{bmatrix}1&0&0&0&0\\0&0&0&1&0\\0&1&0&0&0\\0&0&0&0&1\\0&0&1&0&0\end{bmatrix}}.}


Dado un vector g,

P π g = [ e π ( 1 ) e π ( 2 ) e π ( 3 ) e π ( 4 ) e π ( 5 ) ] [ g 1 g 2 g 3 g 4 g 5 ] = [ g 1 g 4 g 2 g 5 g 3 ] . {\displaystyle P_{\pi }\mathbf {g} ={\begin{bmatrix}\mathbf {e} _{\pi (1)}\\\mathbf {e} _{\pi (2)}\\\mathbf {e} _{\pi (3)}\\\mathbf {e} _{\pi (4)}\\\mathbf {e} _{\pi (5)}\end{bmatrix}}{\begin{bmatrix}g_{1}\\g_{2}\\g_{3}\\g_{4}\\g_{5}\end{bmatrix}}={\begin{bmatrix}g_{1}\\g_{4}\\g_{2}\\g_{5}\\g_{3}\end{bmatrix}}.}

Explicación

Una matriz permutación siempre será de la forma

[ e a 1 e a 2 e a j ] {\displaystyle {\begin{bmatrix}\mathbf {e} _{a_{1}}\\\mathbf {e} _{a_{2}}\\\vdots \\\mathbf {e} _{a_{j}}\\\end{bmatrix}}}

donde eai representa el i vector base (como fila) para Rj, y donde
[ 1 2 j a 1 a 2 a j ] {\displaystyle {\begin{bmatrix}1&2&\ldots &j\\a_{1}&a_{2}&\ldots &a_{j}\end{bmatrix}}}

es la forma permutación de la matriz permutación.

Ahora, al realizar la matriz multiplicación, uno forma esencialmente el producto escalar de cada fila de la primera matriz con cada columna de la segunda. En este caso, estaremos formando el producto escalar de cada fila de esta matriz con el vector de elementos que queremos permutar. Es decir, por ejemplo, v= (g0,...,g5)T,

eai·v=gai

Así, el producto de la matriz permutación con el vector v anterior, será un vector en la forma (ga1, ga2, ..., gaj) , y ésta es pues una permutación de v ya que hemos dicho que la forma de permutación es

( 1 2 j a 1 a 2 a j ) . {\displaystyle {\begin{pmatrix}1&2&\ldots &j\\a_{1}&a_{2}&\ldots &a_{j}\end{pmatrix}}.}

Por lo tanto, las matrices permutación sí permuta el orden de los elementos en vectores multiplicados por ellos.

Véase también

Referencias

  • Esta obra contiene una traducción total derivada de «Permutation matrix» de Wikipedia en inglés, concretamente de esta versión del 5 de abril de 2019, publicada por sus editores bajo la Licencia de documentación libre de GNU y la Licencia Creative Commons Atribución-CompartirIgual 4.0 Internacional.
Control de autoridades
  • Proyectos Wikimedia
  • Wd Datos: Q851512
  • Commonscat Multimedia: Permutation matrices / Q851512

  • Identificadores
  • GND: 4811820-5
  • Wd Datos: Q851512
  • Commonscat Multimedia: Permutation matrices / Q851512