舒尔补加速方程求解

(1) [ H p p H p l H l p H l l ] [ Δ X p Δ X l ] = [ b p b l ] [ I − H p l H l l − 1 0 I ] [ H p p H p l H l p H l l ] [ Δ X p Δ X l ] = [ I − H p l H l l − 1 0 I ] [ b p b l ] [ H p p − H p l H l l − 1 H l p 0 H l p H l l ] [ Δ X p Δ X l ] = [ b p − H p l H l l − 1 b l b l ] \begin{aligned} \begin{bmatrix} H_{pp} & H_{pl} \\ H_{lp} & H_{ll} \end{bmatrix} \begin{bmatrix} \Delta \Chi_{p} \\ \Delta \Chi_{l} \end{bmatrix} &= \begin{bmatrix} b_{p} \\ b_{l} \end{bmatrix} \\ \begin{bmatrix} I & -H_{pl}H^{-1}_{ll} \\ 0 & I \end{bmatrix} \begin{bmatrix} H_{pp} & H_{pl} \\ H_{lp} & H_{ll} \end{bmatrix} \begin{bmatrix} \Delta \Chi_{p} \\ \Delta \Chi_{l} \end{bmatrix} &= \begin{bmatrix} I & -H_{pl}H^{-1}_{ll} \\ 0 & I \end{bmatrix} \begin{bmatrix} b_{p} \\ b_{l} \end{bmatrix} \\ \begin{bmatrix} H_{pp}-H_{pl}H^{-1}_{ll}H_{lp} & 0 \\ H_{lp} & H_{ll} \end{bmatrix} \begin{bmatrix} \Delta \Chi_{p} \\ \Delta \Chi_{l} \end{bmatrix} &= \begin{bmatrix} b_{p}-H_{pl}H^{-1}_{ll}b_l \\ b_{l} \end{bmatrix} \\ \end{aligned} \tag{1} [HppHlpHplHll][ΔXpΔXl][I0HplHll1I][HppHlpHplHll][ΔXpΔXl][HppHplHll1HlpHlp0Hll][ΔXpΔXl]=[bpbl]=[I0HplHll1I][bpbl]=[bpHplHll1blbl](1)

以上均是一般行变换(以左乘矩阵的形式表示),故方程组的解不会发生变化。
【注】:右侧变量b的定义是将负号包含在内的!!
将公式(1)中最后一个等式,重新写成下面两个等式:
(2) [ H p p − H p l H l l − 1 H l p ] Δ X p = b p − H p l H l l − 1 b l H l p Δ X p + H l l Δ X l = b l \begin{aligned} [H_{pp}-H_{pl}H^{-1}_{ll}H_{lp} ] \Delta \Chi_{p}&= b_{p}-H_{pl}H^{-1}_{ll}b_l \\ H_{lp}\Delta \Chi_p+H_{ll}\Delta \Chi_l &= b_l \end{aligned} \tag{2} [HppHplHll1Hlp]ΔXpHlpΔXp+HllΔXl=bpHplHll1bl=bl(2)

故:
(3) Δ X p = [ H p p − H p l H l l − 1 H l p ] − 1 [ b p − H p l H l l − 1 b l ] Δ X l = H l l − 1 [ b l − H l p Δ X p ] \begin{aligned} \Delta \Chi_{p}&= [H_{pp}-H_{pl}H^{-1}_{ll}H_{lp} ]^{-1} [ b_{p}-H_{pl}H^{-1}_{ll}b_l ] \\ \Delta \Chi_l &= H^{-1}_{ll} [ b_l-H_{lp}\Delta \Chi_p] \end{aligned} \tag{3} ΔXpΔXl=[HppHplHll1Hlp]1[bpHplHll1bl]=Hll1[blHlpΔXp](3)

可以加速的原因

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值