LU分解（考虑行交换）

$PA=LU$

对矩阵$A$做$LU$分解（考虑行交换）
$$A=\left[\begin{array}{cccc}{a}_{11}& a_{12}& a_{13}& a_{14}\\ a_{21}& a_{22}& a_{23}& a_{24}\\ a_{31}& a_{32}& a_{33}& a_{34}\\ a_{41}& a_{42}& a_{43}& a_{44}\end{array}\right]$$
第一步，将矩阵$A$的第一列元素$a_{11},a_{21},a_{31},a_{41}$中绝对值最大值所在行与第$1$行交换，即构造行交换矩阵$P_1$，用行交换矩阵$P_1$左乘矩阵$A$完成第一次行交换。如果$a_{41}$的绝对值最大，则将第$4$行和第$1$行交换，对应的行交换矩阵就是将$4$阶单位矩阵的第$1$行和第$4$行交换后得到的矩阵
$$P=\left[\begin{array}{cccc}0& 0& 0& 1\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 1& 0& 0& 0\end{array}\right]$$
用矩阵$P$左乘矩阵$A$相当于交换矩阵$A$的第$1$行和第$4$行
$$\begin{array}{rl}PA& =\left[\begin{array}{cccc}0& 0& 0& 1\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 1& 0& 0& 0\end{array}\right]\left[\begin{array}{cccc}{a}_{11}& a_{12}& a_{13}& a_{14}\\ a_{21}& a_{22}& a_{23}& a_{24}\\ a_{31}& a_{32}& a_{33}& a_{34}\\ a_{41}& a_{42}& a_{43}& a_{44}\end{array}\right]\\ & =\left[\begin{array}{cccc}{a}_{41}& a_{42}& a_{43}& a_{44}\\ a_{21}& a_{22}& a_{23}& a_{24}\\ a_{31}& a_{32}& a_{33}& a_{34}\\ a_{11}& a_{12}& a_{13}& a_{14}\end{array}\right]\end{array}$$
一般的，将第一次行交换后的矩阵记为
$$P_{1}A=\left[\begin{array}{cccc}{b}_{11}& b_{12}& b_{13}& b_{14}\\ b_{21}& b_{22}& b_{23}& b_{24}\\ b_{31}& b_{32}& b_{33}& b_{34}\\ b_{41}& b_{42}& b_{43}& b_{44}\end{array}\right]$$
第二步，构造矩阵$M_1$
$$M_1=\left[\begin{array}{cccc}1& 0& 0& 0\\ -l_{21}& 1& 0& 0\\ -l_{31}& 0& 1& 0\\ -l_{41}& 0& 0& 1\end{array}\right]$$
其中$l_{i1}=\frac{b_{i1}}{b_{11}},i=2,3,4.$ 用矩阵$M_1$左乘$(P_{1}A)$
$$M$$