主题
Search

共轭梯度平方方法


双共轭梯度方法 中,残差向量 r^((i)) 可以被视为 r^((0)) 和关于 Ai 次多项式的乘积,即:

 r^((i))=P_i(A)r^((0)).
(1)

这个相同的多项式满足

 r^~^((i))=P_i(A^(T))r^~^((0))
(2)

因此

rho_i=(r^~^((i)),r^((i)))
(3)
=(P_i(A^(T))r^~^((0)),P_i(A)r^((0)))
(4)
=(r^~^((0)),P_i^2(A)r^((0))).
(5)

这表明,如果 P_i(A)r^((0)) 缩小为一个较小的向量 r^((i)),那么将这个“收缩”算子应用两次并计算 P_i^2(A)r^((0)) 可能是有利的。迭代系数仍然可以从这些向量中恢复(如上所示),并且事实证明很容易找到 x 的相应近似值。这种方法是共轭梯度平方(CGS)方法(Sonneveld 1989)。

通常人们观察到 CGS 的收敛速度大约是 双共轭梯度方法 的两倍,这与同一个“收缩”算子应用两次的观察结果相符。然而,没有理由认为收缩算子,即使它真的减少了初始残差 r^((0)),也应该减少一次减少后的向量 r^((k))=P_k(A)r^((0))。CGS 经常出现高度不规则的收敛行为就证明了这一点。应该意识到,对当前解的局部校正可能非常大,以至于会发生抵消效应。这可能导致解的精度不如更新后的残差所暗示的那么高(van der Vorst 1992)。如果初始猜测接近解,则该方法倾向于发散。

CGS 每次迭代所需的操作次数与 双共轭梯度方法 大致相同,但不涉及与 A^(T) 的计算。因此,在与 A^(T) 的计算不切实际的情况下,CGS 可能很有吸引力。


另请参阅

双共轭梯度方法, 切比雪夫迭代, 正规方程上的共轭梯度法 共轭梯度法, 灵活的广义最小残差法, 广义最小残差法, 线性方程组, 最小残差法, 非平稳迭代法, 预处理器, 拟最小残差法 平稳迭代法, 对称 LQ 方法, 无转置拟最小残差法

此条目由 Noel Black 和 Shirley Moore 贡献,改编自 Barrett 等人 (1994) (作者链接)

使用 Wolfram|Alpha 探索

参考文献

Barrett, R.; Berry, M.; Chan, T. F.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; 和 van der Vorst, H. 线性系统求解模板:迭代方法的构建块,第二版 Philadelphia, PA: SIAM, 1994. http://www.netlib.org/linalg/html_templates/Templates.html.Sonneveld, P. "CGS:非对称线性系统的快速 Lanczos 型求解器。" SIAM 科学与统计计算杂志 10, 36-52, 1989.van der Vorst, H. "Bi-CGSTAB:Bi-CG 的快速平滑收敛变体,用于求解非对称线性系统。" SIAM 科学与统计计算杂志 13, 631-644, 1992.

在 Wolfram|Alpha 中被引用

共轭梯度平方方法

请引用为

Black, NoelMoore, Shirley. "共轭梯度平方方法." 来自 MathWorld--Wolfram Web 资源,由 Eric W. Weisstein 创建. https://mathworld.net.cn/ConjugateGradientSquaredMethod.html

主题分类