36第六节 误差分析.pptVIP

数学学院 信息与计算科学系 第六节、误差分析 　　对线性方程组 Ax=b, 其中设A为非奇异矩阵, x为方程组的精确解. 由于A和b元素是实验测量得到的, 或者是计算的结果, 因此实际是A+?A和b+?b, 方程组变成了 (A+?A)(x+?x)=(b+?b) 下面我们来研究?A或?b与?x的关系，即误差分析. 首先考察一个例子. 例1 设有方程组Ax=b，其中 求的精确解(无任何误差)为 对两组不同的常数项 它们的差只有 但所得解的误差却是 即两组不同的常数项分量误差不过是 可是解的分量误差却高达 对于这样的方程组，不管用什么样的数值方法，我们总难(甚至不可能)算出合理的(与真正精确解相差不大的)解，像这样的方程组称为病态方程组. 定义1 如果矩阵A或常数项b 的微小变化(小扰动), 引起方程组Ax=b解的巨大变化, 则称此方程组为“病态”方程组, 其系数矩阵 A 称为“病态”矩阵(相对于方程组而言), 否则称方程组为“良态”方程组, A称为“良态”矩阵. 应该注意, 矩阵的“病态”性质是矩阵本身的特性,下面我们希望找出刻画矩阵“病态”性质的量. 设有方程组 Ax=b, 其中A为非奇异矩阵, x为方程组的精确解. 以下我们分别研究方程组的常数项b (和A)(即先假设?A＝0，再假设?b＝0)的微小误差(小扰动)时对解的影响. (1) 现设A是精确的, x为Ax=b的精确解，当方程组右端有误差?b, 受扰解为 x+?x, 则 A(x+?x)=b+?b, Ax=b, ?x=A-1?b, ||?x||≤||A-1||·||?b||. （1） 由Ax=b有 ||b||≤||A|| ||x||. 于是得 即解x的相对误差的上界是b的相对误差的||A-1|| ||A||倍. （2） (2) 现设b是精确的, x为Ax=b的精确解，当A有微小误差(小扰动)?A, 受扰解为 x+?x, 则 　　　(A+?A)(x+?x)=b,　有(A+?A)?x= -(?A)x. 而　(A+?A)=A(I+A-1?A). 由3.6节定理7知, ||A-1?A||1时, (I+A-1?A)-1存在. 有 ?x= -(I+A-1?A)-1A-1(?A)x. 因此 （P65定理7） 设||A-1||·||?A||1, 即得 如果?A充分小, 且在条件||A-1||·||?A||1下, 那么此式说明矩阵A的相对误差||?A||/||A||在解x中也可能放大||A-1||·||A||倍. (3) 现设x为Ax=b的精确解，当A有微小误差(小扰动)?A, 而b同时也有微小误差?b(小扰动)时, 受扰解为 x+?x, 则还可以推出相对误差估计式为 总之, 量||A-1||·||A||越小, 由A(或b或两者)的相对误差引起的解的相对误差就越小; 量||A-1||·||A||越大, 解的相对误差就可能越大. 所以量||A-1||·||A||事实上刻画了解对原始数据变化的灵敏程度, 即刻画了方程组的“病态”程度, 于是引进下述定义: 定义2 设A是非奇异矩阵, 称数 condv(A) =||A-1||v||A||v (v=1,2或?) 为矩阵A的条件数 . 由此看出矩阵的条件数与范数有关. 矩阵的条件数是一个十分重要的概念. 由上面讨论知,当A的条件数相对的大, 即cond(A)1时,则方程组是“病态”的(即A是“病态”矩阵, 或者说A是坏条件的, 相对于方程组), 当A的条件数相对的小, 则方程组是“良态”的(或者说A是好条件的). 注意, 方程组病态性质是方程组本身的特性. A的条件数越大, 方程组的病态程度越严重, 也就越难用一般的计算方法求得比较准确的解. 通常使用的条件数, 有 (1) cond∞(A)=||A-1||?||A||?; (2) A的谱条件数; 当A为对称矩阵时 其中?1, ?n为A的绝对值最大和绝对值最小的特征值. 条件数的性质: (1). 对任何非奇异矩阵, 都有condv(A) ≥1. 事实上 condv(cA) =condv(A). (2). 设A为非奇异矩阵且c≠0(常数), 则 (3). 如果A为正交矩阵, 则cond2(A) =1; 如果B为非奇异矩阵, A为