计算机数值分析课件4.pptVIP

华长生制作 * 2.5 平方根法 一、对称正定矩阵的三角分解(Cholesky分解) 记为 因此 Diagonal:对角 为非奇异下三角阵 为非奇异上三角阵 因此 所以 综合以上分析, 则有 定理1. (Cholesky分解) 且该分解式唯一 这种关于对称正定矩阵的分解称为Cholesky分解 二、对称正定线性方程组的解法 线性方程组 则线性方程组可化为两个三角形方程组 对称正定方程 组的平方根法 例1. 用平方根法解对称正定方程组 解: ? - = - = 1 1 2 r k rk rr rr l a l 11 1 1 l a l i i = rr r k rk ik ir ir l l l a l ? - = × - = 1 1 11 11 a l = 分解 T LL 6 29 174 13 29 25 L = 6 29 174 13 29 25 = ) , ( b L 即 = ) , ( y L T 所以原方程组的解为 思考 本例中出现了大量的根式运算 原因为 考虑改变分解方式 三、平方根法的数值稳定性 用平方根法求解对称正定方程组时不需选取主元 由 可知 因此 平方根法是数值稳定的 2.6 追赶法(Thomas算法) 对角占优矩阵: 补充 有一类方程组,在今后要学习的插值问题和边值问题中 有着重要的作用,即三对角线方程组,其形式为: 其中 以