Матричное представление перестановок

Материал из Викиконспекты
Версия от 19:51, 6 января 2017; Shersh (обсуждение | вклад) (Матрица перестановок)
(разн.) ← Предыдущая | Текущая версия (разн.) | Следующая → (разн.)
Перейти к: навигация, поиск

Матрица перестановок[править]

Определение:
Матрица перестановок (англ. Permutation matrix) — квадратная бинарная матрица, в каждой строке и в каждом столбце которой находится лишь одна единица.


Каждая матрица перестановки размера [math]n \times n[/math] является матричным представлением перестановки порядка [math]n[/math].

Пусть дана перестановка [math]\sigma[/math] порядка [math]n[/math]:

[math]\begin{pmatrix} 1 & 2& \ldots & n\\ \sigma(1)& \sigma(2) & \ldots & \sigma(n) \end{pmatrix}[/math]

Соответствующей матрицей перестановки является матрица [math]n \times n[/math] вида:

[math]P_\sigma = \begin{pmatrix} \mathbf{e}_{\sigma(1)}\\ \mathbf{e}_{\sigma(2)}\\ \vdots \\ \mathbf{e}_{\sigma(n)} \end{pmatrix}[/math], где [math]\mathbf{e}_{i}[/math] — двоичный вектор длины [math]n[/math], [math]i[/math]-й элемент которого равен единице, а остальные равны нулю.

Элементарная матрица перестановок[править]

Определение:
Если матрица перестановок [math]P[/math] получена из единичной матрицы [math]E[/math] перестановкой местами двух строк (или двух столбцов), то такая матрица называется элементарной матрицей перестановок (англ. Elementary permutation matrix).


Пример[править]

Пусть дана перестановка: [math]\pi = \begin{pmatrix} 1 & 2 & 3\\ 1 & 3 & 2 \end{pmatrix}[/math]. В соответствующей матрице в первом столбце единица будет стоять на первом месте, во втором столбе на третьем месте, в третьем на втором. Итого: [math]P = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & 1 & 0 \\ \end{pmatrix}[/math]. Также эта матрица является элементарной матрицей перестановок, так как получена из единичной, перестановкой второго и третьего столбцов.

Применение[править]

Благодаря своим свойствам, матрицам перестановок нашлось применение в линейной алгебре. Они используются в элементарных преобразованиях матриц, то есть домножение слева или справа на матрицу перестановок, есть перестановка любых строк или столбов соответственно.

Свойства[править]

Утверждение:
Для любых двух перестановок [math]\sigma, \pi[/math] их матрицы обладают свойством:
[math]P_\sigma P_\pi = P_{\pi \circ \sigma}[/math]
где [math]\circ[/math] — операция умножения перестановок.
[math]\triangleright[/math]
Рассмотрим [math]{P_\sigma P_\pi}_{i,j} = \sum\limits_{x = 1}^{n}{{P_\sigma}_{i,x} {P_\pi}_{x,j}}[/math]. Эта сумма может быть равна нулю или единице, причем единице в том случае, если в [math]i[/math]-той строчке на [math]k[/math]-том столбце матрицы [math]P_\sigma[/math] и в [math]k[/math]-той строчке на [math]j[/math]-том столбце матрицы [math]P_\pi[/math] стоят единицы. [math]{P_\sigma}_{i,k} = 1[/math] значит, что в перестановке [math]\sigma[/math] на [math]i[/math]-том месте стоит элемент [math]k[/math], и [math]{P_\pi}_{k,j} = 1[/math] означает что в перестановке [math]\pi[/math] на [math]k[/math]-том месте стоит элемент [math]j[/math], а [math]{P_\sigma P_\pi}_{i,j} = 1[/math] означает что в перестановке, которой соответствует эта матрица, так же на [math]i[/math]-том месте стоит элемент [math]j[/math]. Но также известно, что [math] (\pi \circ \sigma)(i) = \pi(\sigma(i)) = j [/math]. В результате если [math]{P_\sigma P_\pi}_{i,j} = 1[/math], то [math]{P_{\pi \circ \sigma}}_{i,j} = 1[/math]. Аналогичные рассуждения можно провести когда [math]{P_\sigma P_\pi}_{i,j} = 0[/math], и также получим, что [math]{P_{\pi \circ \sigma}}_{i,j} = 0[/math]. Поэтому для любых [math]i,j[/math] справедливо [math]{P_\sigma P_\pi}_{i,j} = {P_{\pi \circ \sigma}}_{i,j}[/math], а раз такое равентсво выполняется, то [math]P_\sigma P_\pi = P_{\pi \circ \sigma}[/math].
[math]\triangleleft[/math]
Утверждение:
Для любой матрицы перестановок [math]P[/math] справедливо:
[math]P^T P = P P^T = E[/math]
где [math]E[/math] — единичная матрица.
[math]\triangleright[/math]

Рассмотрим [math]{(P P^T)}_{ij} = \sum\limits_{k = 1}^{n}{(P)}_{ik} {(P^T)}_{kj} = \sum\limits_{k=1}^{n} {(P)}_{ik} {(P)}_{jk} = {\delta}_{ij} = {E} [/math]

Теперь в обратную сторону [math]{(P^T P)}_{ij} = \sum\limits_{k = 1}^{n}{(P^T)}_{ik} {(P)}_{kj} = \sum\limits_{k=1}^{n} {(P)}_{ki} {(P)}_{kj} = {\delta}_{ij} = {E} [/math] где [math] {\delta}_{ij} [/math] — символ Кронекера.

Отсюда следует, что [math] P^T=P^{-1} [/math], так как по определению обратной матрицы [math] PP^{-1}=P^{-1}P = E [/math].
[math]\triangleleft[/math]
Утверждение:
При умножение слева матрицы перестановок [math] {P}_{ij} [/math] на матрицу [math]A[/math] происходит перестановка [math] {i} [/math]-й и [math] {j} [/math]-й строк матрицы [math]A[/math]. Умножение справа матрицы перестановок [math] {P}_{ij} [/math] на матрицу [math]A[/math] приводит к перестановке [math] {i} [/math]-го и [math] {j} [/math]-го столбцов матрицы [math]A[/math].
[math]\triangleright[/math]

Рассмотрим сначала умножение слева, т.е. матрицу [math] {P}_{ij}{A} [/math], которую обозначим [math] {B} = {b}_{kl} [/math]. Посчитаем чему равны элементы этой матрицы:

[math] {b}_{kl} = {(\ 0\ \ldots\ 0\ 1\ 0\ \ldots\ 0\ )} \begin {pmatrix} {a}_{1l}\\ {a}_{2l}\\ \vdots\\ {a}_{ml} \end {pmatrix} = \begin {cases} {a}_{kl}, & k \ne i,j,\\ {a}_{jl}, & k = i,\\ {a}_{il}, & k = j. \end {cases} [/math]

Действительно, по определению матрицы перестановок единица в строке стоит на [math] {k} [/math]-м месте, если , [math] {k \ne i,j} [/math], на [math] {j} [/math]-м месте, если [math] {k = i} [/math], и на [math] {i} [/math]-м месте, если [math] {k = j} [/math]. Итак, если [math] {k \ne i,j} [/math], то [math] {k} [/math]-я строка матрицы [math]B[/math] просто совпадает с [math] {k} [/math]-й строкой матрицы [math]A[/math]. Далее, [math] {i} [/math]-я строка матрицы [math]B[/math] совпадает с [math] {j} [/math]-й строкой матрицы [math]A[/math], и наоборот. Поэтому [math]B[/math] получается из [math]A[/math] перестановкой [math] {i} [/math]-й и [math] {j} [/math]-й строк.

Теперь рассмотрим умножение справа. Пусть [math] {B} = {A}{P}_{ij} [/math].

[math] {b}_{kl} = {(\ {a}_{k1}\ {a}_{k2}\ \ldots\ {a}_{kn}\ )} \begin {pmatrix} 0\\ \vdots\\ 0\\ 1\\ 0\\ \vdots\\ 0 \end {pmatrix} = \begin {cases} {a}_{kl}, & l \ne i,j,\\ {a}_{kj}, & l = i,\\ {a}_{ki}, & l = j. \end {cases} [/math]

По определению матрицы перестановок единица в столбце стоит на [math] {l} [/math]-м месте, если [math] {l \ne i,j} [/math], на [math] {j} [/math]-м месте, если [math] {l = i} [/math], и на [math] {i} [/math]-м месте, если [math] {l = j} [/math]. Итак, если [math] {l \ne i,j} [/math], то [math] {l} [/math]-й столбец матрицы [math]B[/math] просто совпадает с [math] {l} [/math]-м столбцом матрицы [math]A[/math]. Далее, [math] {i} [/math]-й столбец матрицы [math]B[/math] совпадает с [math] {j} [/math]-м столбцом матрицы [math]A[/math], и

наоборот. Поэтому [math]B[/math] получается из [math]A[/math] перестановкой [math] {i} [/math]-го и [math] {j} [/math]-го столбцов.
[math]\triangleleft[/math]

Пример

Пусть задана матрица перестановки [math]P = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & 1 & 0 \\ \end{pmatrix}[/math], которая соответствует перестановке [math]\pi = \begin{pmatrix} 1 & 2 & 3 \\ 1 & 3 & 2 \end{pmatrix}[/math], и матрица [math]A = \begin{pmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \\ \end{pmatrix}[/math],

тогда перемножив получим:

  • [math]PA = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & 1 & 0 \\ \end{pmatrix} \times \begin{pmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \\ \end{pmatrix} = \begin{pmatrix} 1 & 2 & 3 \\ 7 & 8 & 9 \\ 4 & 5 & 6 \\ \end{pmatrix}[/math], видно, что вторая и третья строки поменялись местами.
  • [math]AP = \begin{pmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \\ \end{pmatrix} \times \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & 1 & 0 \\ \end{pmatrix} = \begin{pmatrix} 1 & 3 & 2 \\ 4 & 6 & 5 \\ 7 & 9 & 8 \\ \end{pmatrix}[/math], видно, что второй и третий столбец поменялись местами.


Утверждение:
Квадрат элементарной матрицы перестановок есть единичная матрица.
[math]\triangleright[/math]

Любая элементарная матрица перестановок является симметричной матрицей, следовательно [math] \forall{i,j} : {a}_{ij} = {a}_{ji} [/math]. Отсюда следует, что

[math] {P} = {P^T} [/math], а [math] {P P^T} = {E} [/math].
[math]\triangleleft[/math]
Утверждение:
Матрица перестановок [math]n[/math]-го порядка может быть представлена в виде произведения [math](n - 1)[/math] элементарных матриц перестановок [math]{(n \gt 2)}[/math].
[math]\triangleright[/math]

Обозначим [math]{t}_{ij}[/math] — элементарную матрицу, полученную из единичной путем изменения [math]i[/math]-й и [math]j[/math]-й строк. Рассмотрим матрицу перестановок [math] P = \begin {pmatrix} {a}_{11} & {a}_{12} & \ldots & {a}_{1n}\\ {a}_{21} & {a}_{22} & \ldots & {a}_{2n}\\ \vdots & \vdots & \ddots & \vdots\\ {a}_{n1} & {a}_{n2} & \ldots & {a}_{nn} \end {pmatrix}[/math]

Возьмем [math] {{a}_{ij} \ne 0} [/math] и перестановками строк (домножением соответствующей элементарной матрицей слева) или столбцов (домножением соответствующей элементарной матрицей справа) перемещаем его на первое место. Так как в каждой строке или столбце только одна единица, то получим: [math] \begin {pmatrix} 1 & 0 & \ldots & 0\\ 0 & {a}_{22}' & \ldots & {a}_{2n}'\\ \vdots & \vdots & \ddots & \vdots\\ 0 & {a}_{n2}' & \ldots & {a}_{nn}' \end {pmatrix}[/math] и так далее, пока не получится единичной матрицы.

В итоге: [math] t_1 \ldots t_kAt_{k+1} \ldots t_{k+l} = E [/math].

Все элементарные матрицы обратимы и обратная к элементарной матрице — это тоже элементарная матрица, следовательно: [math] A = t_k^{-1} \ldots t_1^{-1}Et_{k+l}^{-1} \ldots t_{k+1}^{-1} = t_k^{-1} \ldots t_1^{-1}t_{k+l}^{-1} \ldots t_{k+1}^{-1} [/math].

Заметим, что с каждым шагом мы домнажаем на одну элементарную матрицу перестановок, следовательно всего будет [math] (n - 1) [/math] таких матриц.
[math]\triangleleft[/math]

См. также[править]

Источники информации[править]