扰动误差分析.docVIP

第2章 线性代数方程组数值解法I：直接法 考虑方程组 (非奇异， 且 ) (2.6.1) 设有误差 有误差，则因此引起解有误差，即有扰动方程组 现在来研究如何通过和对的影响作出估计。 定理2.6.1 设 方程组(2.6.1)中分别有扰动，，因而解向量有误差;又足够小，使得 ，则有误差估计式 证明 由 两边取范数有 得到 得到 又注意到 有 从而得到 ，故上述不等式左边乘以，右边圆括号第一项乘以，第二项乘以，并从括号中提出，则得(2.6.3) 定理的结果实际包含两种特殊情形： A精确，即 ，有扰动，从而 (2) 有扰动，精确，即，这时 当很小时，上式可近似表示为 2．条件数与病态方程组 定义 2.6.1 设为非奇异矩阵，称数 即 为矩阵的条件数。 矩阵的条件数的一些基本性质： 任何非奇异矩阵，对任一算子范数均有 根据定义 ,可得 若 为正交阵，即，则 又非奇异，则 (4)设 与为按模最大和最小的特征值，则 特别地，若（即A对称），则 若对称正定，则 证明 略 定理 2.6.2 （事后误差估计）设方程组 ,非奇异，，是精确解， 是近似解，剩余向量 ,则有估计式 证明： 因，得 ，从而，于是 ，又由 ，于是得估计式右端 类似地，由上述 ，得，或 ，由 得 ，综合两式得估计式两端 例 2.6.1 3．事后误差估计 定理 2.6.2 设 方程组，非奇异，非奇异，，是精确解， 是近似解，剩余向量 ,则 有估计式 例题讲解2 例题2.1 对方程组Ax=b, A非奇异不一定能作顺序Gauss消去过程，或者说，A非奇异不一定有LU分解。 证 这类命题只需举一个反例即可。反例要尽可能简单，这里可考虑2阶、元素为1或0。构造反例一般要经过多次“失败-修正”过程。 考虑A= ，易见A非奇异。显然，顺序Gauss消去过程的第1步就不能进行。从LU分解来看，设有A=LU分解，则有 == 于是有b=0和ab=1同时成立，这就自相矛盾了。 例题2.2 设方程组 = 试手算（或辅以计算器）分别用 （1）顺序Guass消去法 （2）列主元Gauss消去法 （3）直接三角分解法 求解，要求计算中取4位有效数字，最后结果舍入成3位有效数字。上述方程组用4位浮点数进行计算而舍入为3位有效数字的精确解为： 解 （1）顺序Gauss消去过程计算有 第1步消元 第2步消元 回代过程求解并舍入得 （2）列主元Gauss消去过程计算有 选主元换行 第1步消元 选主元，仍为 第2步消元 回代过程求解并舍入 （3）令A=LU，用Doolittle分解计算得 由Ly=b 解得 有Ux=y 解得 与精确解比较可见，列主元Gauss消去法的解精确度最高；其他两种方法的解精确度较差，但彼此接近（这正好符合两种方法实质是一样的情况）。 例题2.3 Gauss消去法的一种自然而又简单的改进是所谓Gauss-Jordan消去法。先考虑顺序Gauss消去法过程的情形。它只是把a这一列中a下面的元素消为0，而Gauss-Jordan消去过程则把a这一列元素的a以外全部消去为0，并且约化a=1。为此，需进行n步消元，第n列也消为只剩下一个元素1，其余均为0（因此，△≠0在这里也是必要的）。这样一来，不用回代过程，方程组的解就在b的位置上。现在依据上述导引，做 试用Gauss-Jordan消去法解方程组 = 试给出Gauss-Jordan算法的核心部分。 由上述两点能得出什么思考？ 解（1）用增广矩阵的演变来描述求解过程 → →→ 于是有x=（2，2，3）。 （2）Gauss-Jordan消去法算法的核心部分为： 对k=1,2,…，n，做 a←(j=k，k+1，…，n，n+1) 对k=1，2，…，n∧i≠k,做 a← a-aa(j=k+1,…,n,n+1) 输出解x=a(i=1,2,…,n) 从上述做法可引发如下思考：Gauss-Jordan消去过程自然也可考虑加上