如果您能分享一些 Ax = b 的共轭梯度迭代不收敛的原因,我将不胜感激?我的矩阵 A 是对称正定矩阵。
太感谢了!
编辑更多信息:我的矩阵是针对简单约束问题的优化算法中的简化 Hessian 矩阵。矩阵是 50% 稀疏的,我正在使用 matlab。
如果您能分享一些 Ax = b 的共轭梯度迭代不收敛的原因,我将不胜感激?我的矩阵 A 是对称正定矩阵。
太感谢了!
编辑更多信息:我的矩阵是针对简单约束问题的优化算法中的简化 Hessian 矩阵。矩阵是 50% 稀疏的,我正在使用 matlab。
如果您的矩阵是对称的、正定的,则 CG 方法可能收敛缓慢,但它收敛于. 它在计算机上不收敛的唯一原因是舍入误差,特别是如果矩阵的条件数,最大和最小特征值的商很大。
经验是,使用双精度算术,即使条件数适中,收敛也会停滞在到的原始错误。