第二章1 解线性方程组的直接法.ppt

* 如果 此时有 如果 ，此时有 上面定理给出了原方程组扰动后解得相对误差的上界。 结论：解的相对误差是在扰动相对误差 和 的基础上通过 量 所刻划的，其大小反映了解的相对误差的大小 。 * 为矩阵A的条件数。 　　如果矩阵范数取2范数，则记 　　定义 设A∈Rn×n为可逆矩阵，按算子范数，称 Cond（A） 同样可以定义cond∞（A）和cond1（A）。 利用条件数的概念，前面的定理结论 可写成 * 可见，条件数比较大时，A和b的小扰动可能会引起解的较大误差，所以条件数刻画了方程组Ax b的性态。一般地，条件数较大时方程组是“病态”的；条件数较小时方程组是“良态”的。当然，病态与良态是相对而言的。 * 例 下列Hilbert矩阵是一族著名的病态矩阵： 它是一个n×n的对称矩阵，可以证明是正定的。计算条件数有cond2（H4） 1.5514 × 104 ， cond2（ H6） 1.4951 × 107，cond2（ H8） 1.525 × 1010。由此可见，随着n的增加， Hn的病态可能越严重。 * 　　对于实际问题，条件数一般是很难计算的。下列现象可能表示方程组Ax b是病态的。 （1）如果矩阵A的按绝对值最大特征值和最小特征值之比很大，则A是病态的。 （2）如果系数矩阵A的元素间数量级差别很大，并且无一定规则，则A可能病态。 　　（3）如果系数矩阵A的某些行或列是近似相关的，或系数矩阵的行列式值相对说很小，则A可能病态。 （4）如果在A的三角化过程中出现绝对值很小的主元或采用选主元技术，主元素数量级相差悬殊时，则A可能病态。 对于给定的非齐次方程组Ax b，要判断它是否病态一般 并不容易。这是因为条件数涉及到系数矩阵的逆矩阵。 zhaojk制作 * 一般病态方程组的求解是比较困难的。 方程组给定后，其系数矩阵的条件数就随之确定，所以 方程组的性态是方程组的固有性质，与求解方法无关，一般来说， 在计算机上求解的方程组都是所给的方程组的扰动方程，这是因为 计算机存储数据的舍入误差所致。对于良态方程组，只要求解方法 稳定，即可得到比较满意的计算结果。但是，对于病态的方程组， 即使用稳定性很好的算法求解也未必得到理想的结果。 对于病态的方程组的求解，可在计算实践中考虑下述的方法： 1.采用高精度的算术运算，如采用双精度的运算，使由于舍入 误差的放大损失若干有效位数之后，还能保留一些有效数字，从 而，改善和减轻“病态”的影响。 2.在求解方程组之前，对原方程组作一些预处理，如进行矩阵 * 平衡，以降低矩阵的条件数。矩阵平衡就是A的每一行（每一列） 分别乘以适当的常数，即找可逆的对角阵 和 ，使方程组 转化为： 理论上应该选择对角阵 和 ，使 其中min是对所有可逆的对角阵集合来说的。 但实际上这很难实现，一般矩阵的平衡只是针对具体问题进行 具体的处理。 一种简单的处理方法是令 ，这称为行平衡。 的选 择可以使 每一行的 范数大体相当，以避免消元过程中 * “大数吃掉小数”。具体的处理方法如下： 设有n阶方程组Ax b，A非奇异，计算 取对角， 则得到与 等价的方程组 ，其系数矩阵 的条件数可能大大低于A的条件数。 * 例：对方程组 系数矩阵的条件数 ，很大，故此方程组是病态 的。 按上述方法进行平衡，取 ，则得到同 解方程组 此时系数矩阵的条件数 ，很小，故 该方程组是良态的。 * * * * 三、平方根法的数值稳定性 用平方根法求解对称正定方程组时不需选取主元 由 可知 因此 平方根法是数值稳定的 事实上,对称正定方程组也可以用顺序Gauss消去法求解 而不必加入选主元步骤 舍入误差的增长可以控制 第2章 线性代数方程组 优点 1.数值稳定. 2. 计算量小，大约为 次乘除法，是一般矩阵A的LU分解 改进的平方根法 缺点 计算lii时要开平方. 设n阶对称正定矩阵A有分解 . 其中L为单位下三角阵， D为对角阵.即 计算量的一半. 第2章 线性代数方程组 n阶正定矩阵 A的 分解公式： 求解对称正定方程组Ax b的改进平方根法 计算公式 ： 1. 分解计算 优点 1.计算量小，大约为 次乘除法 ，同平方根方法，是一般 程组的有效方法. 矩阵A的LU分解 消元法 计算量的一半，是目前解对称正定矩阵方 元素的储存及计算顺序 计算出T LD第i行元素tij j 1,1,…,i-1 ，存放在A第i行位置， 3.精度较高,不用选主元,数值稳定. 元素的储存比消去法节省一半，但比平方根法多用一个单元内存， 2. 计算简单 没有开方运算 . 缺点 元素的储存量大. 计算出L第i行元素lij j 1,1,…,i-1 仍存放在A第i行位置，同时计 算出di .以n 3为例说明 因此透明度很