矩阵三角分解法.ppt

5.3.2 高斯-若当消去法 5.4 矩阵三角分解法 5.4.2 平方根法 5.4.3 追赶法 定理11 如果 为 阶对称正定矩阵，则存在一个实的非奇异下三角阵 当限定 的对角元素为正时，这种 可以用直接分解方法来确定计算 元素的递推公式. (对称正定矩阵的三角分解或Cholesky分解) 使 分解是唯一的. 因为 其中 于是得到解对称正定方程组 的平方根法计算公式： 对于 l. （4.7） 由矩阵乘法及 按等式两边对应元素相等，得 2. 求解 即求解两个三角形方程组： 4. 3. （4.8） 由计算公式1知 所以 于是 这个结果说明，分解过程中元素 的数量级不会增长 且对角元素 恒为正数. 当求出 的第 列元素时， 的第 行元素亦算出. 所以平方根法约需 次乘除法，大约为一般直接分解法计算量的一半. 于是不选主元素的平方根法是一个数值稳定的方法. 由于 为对称阵，因此在计算机实现时只需存储 的 下三角部分. 下三角部分共需存储 个元素，可按行主序用一维数组存放， 即 矩阵元素 在一维数组中表