舒尔补/schur补

在这里插入图片描述

本文地址:https://blog.csdn.net/itnerd/article/details/83385817

n x n 的矩阵可以写成分块形式: M = [ A B C D ] n × n M = \left[\begin{matrix} A & B \\ C & D \end{matrix} \right]_{n \times n} M=[ACBD]n×n其中 A 和 D 为方阵。

  • 若 A 是非奇异的,则 A 在 M 中的舒尔补为: D − C A − 1 B D - CA^{-1}B DCA1B
  • 若 D 是非奇异的,则 D 在 M 中的舒尔补为: A − B D − 1 C A - BD^{-1}C ABD1C

要记住上面的形式很容易,只要记住字母顺序 DCAB、ABDC 都是在 M 中顺时针排列的。


若 A 非奇异,则有

[ I 0 − C A − 1 I ] [ A B C D ] [ I − A − 1 B 0 I ] = [ A 0 0 D − C A − 1 B ] (1) \left[\begin{matrix} I & 0 \\ -CA^{-1} & I \end{matrix} \right] \left[\begin{matrix} A & B \\ C & D \end{matrix} \right] \left[\begin{matrix} I & -A^{-1}B \\ 0 & I \end{matrix} \right] =\left[\begin{matrix} A & 0 \\ 0 & D-CA^{-1}B \end{matrix} \right]\tag{1} [ICA10I][ACBD][I0A1BI]=[A00DCA1B](1)
上式可以看作对矩阵 M 实施分块矩阵的初等行列变换。由 (1) 可得 ∣ A B C D ∣ = ∣ A 0 0 D − C A − 1 B ∣ = ∣ A ∣ ∣ D − C A − 1 B ∣ \left|\begin{matrix} A & B \\ C & D \end{matrix} \right| = \left|\begin{matrix} A & 0 \\ 0 & D-CA^{-1}B \end{matrix} \right| = |A||D-CA^{-1}B| ACBD=A00DCA1B=ADCA1B此时 M 非 奇 异 ⟺ D − C A − 1 B 非 奇 异 M 非奇异 \Longleftrightarrow D-CA^{-1}B 非奇异 MDCA1B
同理,当 D 非奇异时 M 非 奇 异 ⟺ A − B D − 1 C 非 奇 异 。 M 非奇异 \Longleftrightarrow A-BD^{-1}C 非奇异。 MABD1C


分块矩阵的逆

由式 (1) 还可以方便地求得分块矩阵的逆。
[ I − A − 1 B 0 I ] − 1 [ A B C D ] − 1 [ I 0 − C A − 1 I ] − 1 = [ A − 1 0 0 ( D − C A − 1 B ) − 1 ] \left[\begin{matrix} I & -A^{-1}B \\ 0 & I \end{matrix} \right]^{-1} \left[\begin{matrix} A & B \\ C & D \end{matrix} \right]^{-1} \left[\begin{matrix} I & 0 \\ -CA^{-1} & I \end{matrix} \right]^{-1} =\left[\begin{matrix} A^{-1} & 0 \\ 0 & (D-CA^{-1}B)^{-1} \end{matrix} \right] [I0A1BI]1[ACBD]1[ICA10I]1=[A100(DCA1B)1]因而 [ A B C D ] − 1 = [ I − A − 1 B 0 I ] [ A − 1 0 0 ( D − C A − 1 B ) − 1 ] [ I 0 − C A − 1 I ] \left[\begin{matrix} A & B \\ C & D \end{matrix} \right]^{-1} =\left[\begin{matrix} I & -A^{-1}B \\ 0 & I \end{matrix} \right] \left[\begin{matrix} A^{-1} & 0 \\ 0 & (D-CA^{-1}B)^{-1} \end{matrix} \right] \left[\begin{matrix} I & 0 \\ -CA^{-1} & I \end{matrix} \right] [ACBD]1=[I0A1BI][A100(DCA1B)1][ICA10I]

在这里插入图片描述

  • 28
    点赞
  • 5
    评论
  • 55
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

打赏
文章很值,打赏犒劳作者一下
相关推荐
©️2020 CSDN 皮肤主题: 程序猿惹谁了 设计师:白松林 返回首页

打赏

颹蕭蕭

白嫖?

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者