舒尔补/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]

在这里插入图片描述

  • 36
    点赞
  • 100
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值