数值线性代数课程设的计.pdfVIP

Cholesky 分解法即平方根法，是求解对称正定 数值课程设计 线性方程组最常用也是最有效的方法之一，并且在 Cholesky 分解法及改进的Cholesky 分 许多实际问题中，需要求解的线性方程组的系数是 解法解对称正定方程组的 学原理及实例分 实对称正定的。因此，对于这一类方程组，可以把 析 直接三角分解法简化，即对系数矩阵A 的LU分解, 简化为A = LL ，其中A 是n 阶实对称正定阵，L 为 信息与计算科学班 非奇异的下三角矩阵，且当L 的主对角元均为正时， 摘要：线性方程组的解法是一个古老的 学问题， 这种分解方法是唯一的。这样，对于线性方程组 而Gauss 消去法是适用于求解中小规模的系 矩阵 Ax = b，当把A 分解为A = LL 时，线性方程组变为 稠密的线性方程组的一种有效方法，但是对于大型 LLx = b，令y = Lx ，有Ly = b，再利用前代法解 稀疏稀疏矩阵的线性方程组，则有更为有效的方法 Ly = b，得y ，用回代法解L x = y，得 x ，方程组 求解，而对于工程实际当中常遇到的系 矩阵是对 Ax = b的解x 。 称正定矩阵的线性方程组，Cholesky 分解法 （平方 从后文的问题实例分析中可以看出，利用 根法）是解这类方程组的最常用和最有效的方法。 Cholesky 分解法解对称正定方程组时，计算L的对 关键词：线性方程组 对称正定 Cholesky 分解法 角元素 , 时需要用到开方运算，使最终得到的解误 差增大，为了避免开方，我们可以通过改进 Cholesky 分解法，使A 具有如下形式的分解： 目录 A = LDL ，其中A是对称正定阵，D是对角元素均为 第一章 概述……………………………………1 正 的对角矩阵，这一分解也称为LDL 分解或改进 第 章 实例分析………………………………1 的平方根法，通过后文的误差图像可以看出，改进 2.1 问题介绍………………………1 Cholesky 分解法后，可以减少解的误差。用改进的 2.2 问题解答……………………… 2 Cholesky 分解法解Ax = b时，A = LDL ，此时线性 2.2.1 问题（1）………………2 方程组变为LDLx = b，令y = DLx ，有Ly = b，用 2.2.2 问题（2）………………2 前代法解Ly = b，得y ，所以DL 仍是上三角阵，可 2.2.3 方法比较……………...3 再用回代法解DLx = y，得x ，即得Ax = b的解x 。 第三章 Cholesky 分解法的 学原理………3 3.1 Cholesky 分解法介绍…………3 3.2 Cholesky 分解法的 学原理… 3 3.3 改进的Cholesky 分解法与Cholesky 第 章 实例分析 分解法的比较…………………………3 2.1 问题介绍 第四章 结论……………………………………4 利用 MathCAD 将 Cholesky 分解法和改进的