第二章线性代数方程组的数值方法.ppt

对称矩阵的三角分解 （2-78） 对称矩阵的三角分解 现将A阵作如下的三角分解 这里 ，由式（2-78）得 于是对称阵A可分解为 （2-81） 其中 是下三角阵。 因此可定义对角阵 这样 （2-83） 以上说明，对称正定矩阵可分解为两个互为转置的三角阵之积。 对称矩阵的三角分解 若A为对称正定矩阵，可以证明D亦正定，从而di0。 A阵不仅对称而且正定 乔列斯基分解的一般公式 由矩阵乘法，可直接获得S的计算公式 对称矩阵的三角分解 （2-85） 式（2-85）就是乔列斯基分解的一般公式。 对称矩阵的三角分解 对称矩阵的三角分解 求S阵 第一列 第二列 第三列 依次逐列计算。 对称系数阵方程组的乔列斯基解法的步骤： 第一步：将A阵按式（2-85）分解，得到S阵； 第二步：由 求解向量Y，有 第三步：由 求解向量X，有 （2-86） （2-87） 对称矩阵的三角分解 对称矩阵的三角分解 用平方根法求解 举例 解 对称矩阵的三角分解 由 求解向量Y，有 由 求解向量X，有 X=（2，1，-1） T Y=（3，0.5，-1）T 对称矩阵的三角分解 题6 用平方根法求解 对称矩阵的三角分解 六、三对角阵的三角分解——追赶法 追赶法是适用于三对角线方程组的有效解法。 三对角方程组形式： AX=D 三对角阵的三角分解 其系数阵是三对角阵 将A阵分解为L阵和U阵。 三对角阵的三角分解 验证两端必有 三对角阵的三角分解 三对角阵的三角分解算式为 解三对角线方程组的追赶法的步骤 （1）将A阵按式（2-89）分解，得到L、U阵。 （2）向前回代求解LY=D，得 （3）向后回代求解UX=Y，得到 （2-89） 追 赶 三对角阵的三角分解 用追赶法解三对角方程组 举例 解: 三对角阵的三角分解 三对角阵的三角分解 由LY=D，得 由UX=Y，得 三对角阵的三角分解 题7 用追赶法解三对角方程组 三对角阵的三角分解 七、方程组的逆矩阵解法与矩阵求逆 （一）方程组的逆矩阵解法 对方程 的两端左乘以逆矩阵A-1得 解为 方程组的逆矩阵解法 （二）矩阵求逆 用计算机求n×n阶非奇异方阵A的逆矩阵A-1，等价于求解具有相同系数阵（被求逆的矩阵）且右端项分别为 的n个方程组，即求解下述方程组 或改写为 以A-1左乘上式两端，则有 （2-94） 列向量Xi是用前面讲过的任一种解法针对AX=B而得到的解，它们构成A-1的相应的列，所以计算机求逆的一个实用方法就是求解方程组系式（2-94）。 求解（2-94）的一个简便方法是考虑增广矩阵 应用高斯消去法、主元消去法或高斯——约当消去法解n×2n阶的方阵[AI]。 方程组的逆矩阵解法 举例 用高斯——约当消去法计算矩阵 的逆矩阵。 解 所求的逆矩阵为 方程组的逆矩阵解法 解 举例：全主元法 解为 x2和x3次序变化 全主元法的精度略优于列主元法，这是由于全主元法是在全体系数中选主元，故它对控制舍入误差比较有效。但全主元法在计算过程中，需同时作行与列的互换，因而程序比较复杂，计算时间较长。列主元法的精度虽稍低于全主元法，但其计算简单，工作量大为减少，实践表明，它与全主元法同样具有良好的数值稳定性，故列主元法是求解中小型稠密线性方程组的最好方法之一。一般采用列主元法。 说明： 主 元 消 去 法 题3 分别用列主元消去法和全主元消去法求解线性方程组。 精确解舍入到4位有效数字为 主 元 消 去 法 关于主元法的几点讨论 （1）有两类系数阵可以避免选主元，直接使用高斯消去法便可得到满意的结果。它们是：严格对角占优矩阵和对称正定矩阵。 严格对角占优矩阵是指其元素满足 或以列为准，其元素满足 严格对角占优矩阵必定是非奇异阵，且每步消元后，仍保持对角占优的性质，因此，不必选主元。 对称正定矩阵其行列式大于零，方程组可解，且在每步消元中，剩下的子阵仍保持对角占优和对称性，因此不必选主元。 行 列 （2）关于病态方程 其具有5位有效数字的正确解为 如把其系数阵也取成5位有效数字 给出的错误解答为 第一步消元后 该方程为一个病态方程。 判断是否为病态方程: 行列式很大或很小（如某些行列式近似相关）； 元素间相差大数量级，且无规则； 主元消去过程中出现小主元； 特征值相差大数量级。 解决方