4.5 置换矩阵

4.5 置换矩阵

是不是任意可逆矩阵都可进行 $LDU$ 分解呢？其实不能，消元操作需要除以对角元素 $a_{ii}$ ，当其为 $0$ 时，则会失败。这时可在下面行中选择任一对角元素不为 $0$ 的行，对调这两行，则可继续消元。例如

$$A=\left[\begin{array}{ccc}0& 0& 2\\ 1& 2& 3\\ 0& 1& 2\end{array}\right]$$

第一行第一个元素 $a_{11}$ 为 $0$ ，无法消除第二行第一列的非零元素。矩阵后面两行中，第二行第一个元素 $a_{21}$ 非零，则对调这两行，矩阵变换为
$$\left[\begin{array}{ccc}1& 2& 3\\ 0& 0& 2\\ 0& 1& 2\end{array}\right]$$

此时第一列元素除对角线外已经都是 $0$ 。同理消除第二列时，第二行对角线元素为 $0$ ，此时也需要对调后两行，矩阵变换为

$$\left[\begin{array}{ccc}1& 2& 3\\ 0& 1& 2\\ 0& 0& 2\end{array}\right]$$

成为上三角阵。

重要性质 对任意可逆矩阵，经过适当的行对调操作，可以分解为 $LDU$ 。

类似消元操作，行对调操作也可以用矩阵乘法实现，该矩阵称为置换矩阵。

定义 置换矩阵 矩阵 $P_{ij}$ 是单位矩阵 $E$ 对调 $i,j$ 两行所得。

例如
$$P_{12}=\left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 0\\ 0& 0& 1\end{array}\right]P_{13}=\left[\begin{array}{ccc}0& 0& 1\\ 0& 1& 0\\ 1& 0& 0\end{array}\right]P_{23}=\left[\begin{array}{ccc}1& 0& 0\\ 0& 0& 1\\ 0& 1& 0\end{array}\right]$$

矩阵 $A$ 列向量左乘置换矩阵 $P_{ij}$ 就是对调向量的 $i,j$ 两个分量。
$$P_{ij}{\mathbf{a}}_k=P_{ij}(a_{1k},\cdots ,a_{ik},\cdots ,a_{jk},\cdots ,a_{mk})=(a_{1k},\cdots ,a_{jk},\cdots ,a_{ik},\cdots ,a_{mk})$$

$P_{ij}A$ 就是对调矩阵 $A$ 的 $(i,j)$ 两行。

置换矩阵是正交矩阵， $P^{T}P=E$ 。对矩阵进行多次行对调操作，就是多个置换矩阵连乘，记为 $P$ ， $P$ 是单位矩阵 $E$ 进行相应的多次行对调结果。

$$P=P_{21}{P}_{32}=\left[\begin{array}{ccc}0& 0& 1\\ 1& 0& 0\\ 0& 1& 0\end{array}\right]$$

重要性质 对任意可逆矩阵 $A$，经过适当的行对调操作 $P$，可以分解为 $PA=LDU$ 。

我们还可以换个角度看待 $PA=LDU$ ，由于各矩阵均可逆，得 $(LDU{)}^{-1}PA=E$ ，令 $P^{\prime }=(LDU{)}^{-1}P$ 则 $P^{\prime }A=E$ ，这说明 $P^{\prime }$ 是逆矩阵 $A^{-1}$。通过高斯消元法可得到逆矩阵 $A^{-1}=U^{-1}{D}^{-1}{L}^{-1}P$ ，对角阵 $D$ 可逆，需对角元素均不为零，故矩阵 $A$ 主元均不为零时，矩阵 $A$ 可逆。

当矩阵 $A$ 是对称矩阵时，假设没有行对调，则 $S=LDU$ ，取转置，$S^T=(LDU{)}^T=U^T{D}^T{L}^T=U^{T}D{L}^T=S=LDU$ ，所以有 $L^T=U$ 成立。

重要性质 对称矩阵，假设没有行对调，则可以分解为 $S=LD{L}^T$ 。