大型稀疏法方程组极优化分解算法

摘要

针对大型法方程和误差方程解算问题,提出为矩阵寻求一个好的排序,使得它的Cholesky分解后的下三角矩阵有较少的非零元素个数。同时介绍稀疏对称正定矩阵的Cholesky结构分解,它可先确定下三角矩阵中非零元素的位置和个数,从而预先分配其存储空间,且在数值计算时,只需按照已安排好的非零元素的位置进行计算,即可得到数值解。最后给出实现方法和数值试验结果。