第3章线性方程组的数值解法shopo.ppt

3.7 误差分析 3.7.1 方程组的性态 在建立方程组时，其系数往往含有误差（如观测误差或计算误差），就是说，所要求解的运算是有扰动的方程组，因此需要研究扰动对解的影响。 例3.14 考察方程组 和 上述两个方程组尽管只是右端项有微小扰动, 但解大不相同, 第1个方程组的解是 第2个方程组的解是 。这类方程组称为病态的。 例3.15考察方程组： 精确解为 设方程组存在扰动 精确解为 上例说明该方程组的解对初始元素的扰动非常敏感。 定义3.8 A或b的微小变化(又称扰动或摄动)引起方程组Ax=b解的巨大变化，则称方程组为病态方程组，矩阵A称为病态矩阵。否则方程组是良态方程组，矩阵A也是良态矩阵 为了定量地刻画方程组“病态”的程度，要对方程组Ax=b进行讨论，考察A（或b）微小误差对解的影响。为此先引入矩阵条件数的概念。 定义3.9 （矩阵条件数）设A为非奇异矩阵，称 为矩阵A条件数。 设方程组为 系数矩阵 和常数向量 的扰动分别记为： 和 实际求解的方程组为 对应例3.15即为 使解 变为 如果 和 很小，而解的变动 很大，则称方程组 为病态方程组，称系数矩阵 为关于求解方程组或求逆的病态矩阵； 反之，如果 和 微小时， 也微小，则称方程组 为良态方程组，称系数矩阵 为关于求解方程组或求逆的良态矩阵。 病态方程组对任何算法都将产生数值不稳定性 定理 3.12 (b的扰动对解的影响) 设A非奇异， Ax=b≠0，且 则有 证:设A精确且非奇异,b有扰动δb,使解x有扰动δx, 则 消去Ax=b，有 又 两边相乘再同除x和b的范数 定理 3.13 (A的扰动对解的影响)设A非奇异， Ax=b≠0，且 若， 则 证 ：见p66 我们还可证明更为一般的结论： 当方程组的系数矩阵A非奇异和常数项b为非零向量时，且同时有扰动δA，δb，满足 ，若x和x+δx分别是方程组 Ax=b 及 的解 则 例3.16 线性方程组 的系数矩阵带误差，成为如下方程组 求方程组系数矩阵的条件数, 并说明方程组的性态 解: 因为 所以 因此方程组是良态的 例3.16 Hilbert矩阵就是一个著名的病态矩阵 n 4 6 8 10 28375 2.9E+7 3.39E+10 3.54E+13 15514 1. 5E+7 1.53E+10 1.60E+13 28375 2.9E+7 3.39E+10 3.54E+13 对称正定矩阵 3.7.2 精度分析 求得方程组Ax=b的一个近似解以后,希望判断其精度，检验精度的一个简单办法是将近似解再回代到原方程组去求出余量r. r = b-A 如果r很小，就认为解是相当精确的。 定理3.14 设 是方程组Ax=b的一个近似解,其精确解记为 ,r为 的余量。则有 证明见P68 分析：由正交矩阵和条件数的定义便可推得 解：因为A是正交矩阵， 故ATA= AAT=I， A-1= AT，从而 例3.17 设A为正交矩阵，证明：cond2(A)=1 cond (AB) = || AB || · || (AB)-1 || ≤ ||A|| · ||B|| ||A-1|| · ||B-1|| = ||A|| · ||A-1|| ||B|| · ||B-1|| = cond (A) · cond (B) 例3.18 设A，B为n阶矩阵，证明： cond(AB)≤ cond(A) · cond(B) 分析: 由矩阵范数性质和条件数定义 便可证明 证： 例3.19 设A，B为n阶非奇异矩阵，||?||表示 矩阵的任一种范数，证明： || A-1-B-1 || ≤|| A-1 || || B-1 || || A-B || cond(AB)≤ cond(A) · cond(B) 分析: 由矩阵范数的基本性质即可推证 证：