置换矩阵就是重新排列后的单位矩阵
对一个矩阵进行行交换，需要通过置换矩阵（permutation matrix）来完成。
在对一个Ax=b的方程组进行高斯消元的过程中，我们常常会遇到一种情况，也就是消元消不下去的情况。下面，我列出了3×3矩阵的两次消元过程，第一行，是正常的消元流程，对于方阵而言，先令主对角线线上的第一个元素A11为主元（pivot），同时当前列也被称为主元列。用A11所在的行去乘以一个合适的系数减去下面各行（第二，第三行），使得A11这一列下面的元素都为0。之后，再令A22为主元，A22所在的列为主元列，用A22所在的行乘以一个相应的系数减去下面各行（第三行），使得A22所在列下面的元素都为0。一直进行到，对角线上的最后一个元素。
但有时候，我们会遇到主对角线上为0的情况，正如我下面写的第二个消元过程。当消元进行到主对角线上的A22时，消元无法进行下去了，因为，0无论乘以任何继续都无法消除0下面的所有元素。这时，我们要找到该列中A22下面的非零元素所在的行，和A22所在的行，进行行交换。最终使得主对角线上的A22不为0,并完成后续的消元过程。
而，行交换的矩阵表达形式就是用置换矩阵P乘以该矩阵。例如，我们下面这个例子：
消元一开始就遇到了问题，主对角线上的第一个元素A11=0，消元进行不下去了。这时我们只需把上面的两行对调一下就可以得到下面的方程组。
而如果我们把这一过程用矩阵的形式来表现的话就是：
重点：现在我要重复一遍我这篇文章开头的标题，置换矩阵就是重新排列后的单位矩阵。实际上，最简单的置换矩阵P就是单位矩阵I本身，他的作用就是对矩阵不做任何处理。
P = I（no exchanges）
为了更加直观的描述一个置换矩阵的功能，我们可以用矩阵来表示置换矩阵。用去左乘另一个矩阵，会交换这个矩阵I的i行和j行。如果我们仔细比较置换矩阵中的i，j行和相同尺寸的单位矩阵中的i，j行，我们会发现，他就是交换了单位矩阵中ij行后的矩阵。
它有如下的一些特性：
1，对任何一个置换矩阵P而言，每一行，梅一列都只有一个“1”。
2，同时，置换矩阵的转置也是一个置换矩阵，可能是同一个矩阵P，可能是另一个矩阵。
3，任意两个置换矩阵的乘积也是一个置换矩阵。
置换矩阵也可以是一次交换多行的矩阵P，这需要通过两个置换矩阵相乘得到：
例子：
4，对于任何n阶矩阵，共有n的阶乘种置换矩阵。其中，n！=1x2x3…xn。因此，若n=4就对应有24种置换矩阵，若n=5，就会有120种置换矩阵！
5，置换矩阵的逆矩阵就是他自己。已知，置换矩阵的作用就是交换另一个矩阵的i行和j行，要想恢复这个操作，那就是把另一个矩阵的i行和j行再交换一次，所以：
注意，第五条所说的性质，只争对进行一次行交换的置换矩阵，但，如果是一次执行多次行交换的置换矩阵，例如P=P21P32，他的矩阵和他的逆不同，如下：
但他的逆确实实现了对p操作的还原，如下：
6，置换矩阵被广泛的用于对非奇异矩阵（nonsingular matrix）进行高斯消元时，遇到主元为零时需要进行行交换的情况。即：PA=LU（或者PA=LDU）。对于奇异矩阵（singular matrix）,行交换也无法彻底杜绝主元为0的情形。
LU分解是一个伟大的矩阵分解！
以下是个人笔记：
置换矩阵：
置换矩阵的乘法：
置换矩阵的性质：
注：上文中提到的，置换矩阵P可用于列交换，需要进一步考证，待定！
（全文完）
作者 — 松下J27
格言摘抄：人就是三个口！（无名氏）
鸣谢：
文中截图均来自于：《Introduction to Linear Algebra》,5th Edition – Gilbert Strang
（配图与本文无关）
版权声明：所有的笔记，可能来自很多不同的网站和说明，在此没法一一列出，如有侵权，请告知，立即删除。欢迎大家转载，但是，如果有人引用或者COPY我的文章，必须在你的文章中注明你所使用的图片或者文字来自于我的文章，否则，侵权必究。 —-松下J27
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/217095.html原文链接：https://javaforall.net
