2.10 分块矩阵求逆

分块矩阵求逆

分块矩阵求逆法，把大型矩阵变成小型矩阵，可以提高计算效率。

1、准对角矩阵 $A=\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ \mathbf{O}& A_{22}\end{array}\right]$ 的逆，可以待定系数法求，设逆矩阵为 $A^{-1}=\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]$ ，则
$$\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ \mathbf{O}& A_{22}\end{array}\right]=\left[\begin{array}{cc}{E}_n& \mathbf{O}\\ \mathbf{O}& E_m\end{array}\right]$$

计算

$$\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ \mathbf{O}& A_{22}\end{array}\right]=\left[\begin{array}{cc}X{A}_{11}& Y{A}_{22}\\ Z{A}_{11}& W{A}_{22}\end{array}\right]=\left[\begin{array}{cc}{E}_n& \mathbf{O}\\ \mathbf{O}& E_m\end{array}\right]$$

令对应元素相等，得

$$\begin{gathered} X{A}_{11}=E_{n}Y{A}_{22}=\mathbf{O} \\ Z{A}_{11}=\mathbf{O}W{A}_{22}=E_m \end{gathered}$$

由于矩阵 $A_{11},A_{22}$ ，可得 $X=A_{11}^{-1},Y=\mathbf{O},Z=\mathbf{O},W=A_{22}^{-1}$ ，所以 $A^{-1}=\left[\begin{array}{cc}{A}_{11}^{-1}& \mathbf{O}\\ \mathbf{O}& A_{22}^{-1}\end{array}\right]$ 。

2、准下角矩阵 $A=\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ A_{21}& A_{22}\end{array}\right]$ 的逆，可以待定系数法求，设逆矩阵为 $A^{-1}=\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]$ ，则
$$\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ A_{21}& A_{22}\end{array}\right]=\left[\begin{array}{cc}{E}_n& 0\\ 0& E_m\end{array}\right]$$

计算

$$\left[\begin{array}{cc}X& Y\\ Z& W\end{array}\right]\left[\begin{array}{cc}{A}_{11}& \mathbf{O}\\ A_{21}& A_{22}\end{array}\right]=\left[\begin{array}{cc}X{A}_{11}+Y{A}_{21}& Y{A}_{22}\\ Z{A}_{11}+W{A}_{21}& W{A}_{22}\end{array}\right]=\left[\begin{array}{cc}{E}_n& \mathbf{O}\\ \mathbf{O}& E_m\end{array}\right]$$

令对应元素相等，得

$$\begin{gathered} X{A}_{11}+Y{A}_{21}=E_{n}Y{A}_{22}=\mathbf{O} \\ Z{A}_{11}+W{A}_{21}=\mathbf{O}W{A}_{22}=E_m \end{gathered}$$

由于矩阵 $A_{11},A_{22}$ ，可得 $X=A_{11}^{-1},Y=\mathbf{O},Z=-A_{22}^{-1}{A}_{21}{A}_{11}^{-1},W=A_{22}^{-1}$ ，所以 $A^{-1}=
\left[ \begin{matrix}
A^{-1}{11} & \mathbf{O} \
-A^{{-1}_{22}A_{21}A}{-1}{11} & A^{-1}_{22} \
\end{matrix} \right] $ 。

3、准上角矩阵 $A=\left[\begin{array}{cc}{A}_{11}& A_{12}\\ \mathbf{O}& A_{22}\end{array}\right]$ 的逆 $A^{-1}=\left[\begin{array}{cc}{A}_{11}^{-1}& -A_{11}^{-1}{A}_{12}{A}_{22}^{-1}\\ \mathbf{O}& A_{22}^{-1}\end{array}\right]$ 。

上面三个例子显示，如果子矩阵是 $\mathbf{O}$ 矩阵，逆矩阵对应的子矩阵也是 $\mathbf{O}$ 矩阵，这些位置的元素不需要计算，直接是 $0$ ，在实际计算逆矩阵时，可以提高效率。
分块矩阵计算逆矩阵，如3阶准上角矩阵，$A=\left[\begin{array}{ccc}1& 3& 5\\ 2& 4& 6\\ 0& 0& 7\end{array}\right]$ ，4个子矩阵分别为 $A_{11}=\left[\begin{array}{cc}1& 3\\ 2& 4\end{array}\right],A_{12}=\left[\begin{array}{c}5\\ 6\end{array}\right],A_{21}=\left[\begin{array}{cc}0& 0\end{array}\right],A_{22}=\left[\begin{array}{c}7\end{array}\right]$ 。根据2阶逆矩阵公式，得 $A_{11}^{-1}=\left[\begin{array}{cc}4& -3\\ -2& 1\end{array}\right]/(-2)$ ， $A_{22}^{-1}=\left[\begin{array}{c}1/7\end{array}\right]$ ，$-A_{11}^{-1}{A}_{12}{A}_{22}^{-1}=\left[\begin{array}{c}1\\ -2\end{array}\right]/7$ , 所以逆矩阵为 $A^{-1}=\left[\begin{array}{ccc}-2& 3/2& 1/7\\ 1& -1/2& -2/7\\ 0& 0& 1/7\end{array}\right]$ 。 读者可以验证 $A{A}^{-1}=E,A^{-1}A=E$ 。比直接计算3阶矩阵的逆要快很多。