数值第5章解线性方程组的直接方法摘要.pptVIP

定义4 如果 的某个非负实函数N(A)，记作‖A‖， 满足条件： 则称 显然 满足定义中的四个条件, (3)，(4)两条均可由 Cauchy-Schwarz不等式证明，故 是一种矩阵范数. 除矩阵自身的运算外，在解方程中矩阵乘向量的运算即Ax，也是必不可少的.因此要求所引进的范数应满足条件： 上式称为相容性条件.为使引进的矩阵范数满足条件(4.5)，我们给出以下定义. (4.5) 定义6(矩阵的算子范数) 设 当给定向量范数 时可定义 称为矩阵的算子范数或从属范数. (4.6) 定理17 设 上的一种向量范数， 则由(4.6)定义的 是一种矩阵范数，且满足相容性条件 证明 因 中有界闭集 上的连续函数，故 在D上有最大值，即 使 而对 故 所以 从而当 成立，而x=0时 显然也成立. 定理17 设 则 这里 为矩阵的谱半径. 例7 已知 解 从定理可以看出，计算 较容易，而计算 时因为要求 的特征值，所以较为困难.但当A对称时，有 定理19 定理18 对任何 为任一种从属范数则 反之，对任意ε＞0，至少存在一种从属范数 使 证明:设 为A的特征值，则 由 得 非奇异，且 证明　用反证法.假定(I+B)奇异，则齐次方程 有非零解 而 与‖B‖＜1的假设矛盾，故(I+B)非奇异. 又 得 取范数 得 定理20　设 5.6.1矩阵条件数与扰动方程组误差界　　在解方程组Ax=b时，由于各种原因，A或b往往有误差，从而使得解也产生误差. 例8 方程组  的准确解为 ,当A与b有微小变化时，如变为方程 则准确解为 它表明A,b的微小扰动引起方程解x的很大 变化，这就是病态方程. §5.6 误差分析与病态方程组 定义7 求解线性方程组Ax=b时，若A或b有微小扰动 解x的误差 很大, ，则称此方程组为 病态方程组，相应的系数矩阵A称为病态矩阵. 反之，若此时 很小, ，则称此方程组为 良态方程组，相应的系数矩阵A称为良态矩阵. 注意方程组是否病态与用什么数值方法无关，它是由方程自身性质决定的. 在例8中因为行列式 因此出现病态.但有时A从表面上看性质很好，也可能是病态的. 　那么如何判断A是否病态?先给出如下定义. 例9 方程组Ax=b表示为 它的准确解 A对称正定且 表面看性质较好，但若对右端b作微小变化，如方程改为 则解变为 这里b的相对误差大约只有 但解的相对误差却很大， 故A也是病态矩阵. 定义8 设 非奇异，‖·‖v为矩阵的任一种从属范数，则 称为矩阵A的条件数. 从定义看到矩阵条件数依赖于范数的选取，如范数为2-范数， 则记为 同理有 等等. 条件数有以下性质： (1) (2) (3) U为正交矩阵，则 (4) 若 与 为A的按模最大与最小特征值，则 若A对称，则 下面给出扰动方程组解的误差分析.先考察b有扰动 则扰动方程为 由于Ax=b,故得 于是 再由Ax=b,有 即 故得 下面再研究方程Ax=b，当A有扰动 时，其解 的误差分析. 此时扰动方程为 因Ax=b,故有 因 存在, 若假定 则由定理20可知 非奇异, 并有: (5.6) 由(5.6)可得 因此 (5.7) 定理22 设A为非奇异矩阵,Ax=b≠0,且 如果 则(5.7)式成立. 从(5.7)看到，当A的条件数Cond(A)很大时，解的相对误差 也很大，故方程组为病态.在例9中 而 于是 条件数很大，故方程是严重病态的. 例10 Hibert矩阵是一个著名的病态矩阵，记作 它是一个对称正定矩阵，当n≥3时它是病态矩阵.例如 故 另外还有 等等. 因此Hn是严重的病态矩阵，且n越大 Cond(Hn)越大. 例11 在例10的方程组中可算出A的特征值 故 例中 实际相对误差是 而根据(5.6)的误差估计为 这与实际相差不大，即相对误差放大了将近3 000倍.故方程为病态方程组. 　定理23(事后误差估计)　设方程组 ?，则 若实际求得解为 证明　记剩余 则 它表明如果方程组病态，即使剩余‖r‖很小，解的相对误差仍可能很大. 5.6.2 病态方程组的解法 如果A的条件数Cond(A)1，则Ax=b为病态方程，但计算Cond(A)时需要求 A-1,计算量很大，相当于解方程组， 在实际中常可通过求解过程直观地判断方程组的病态性质，如果解方程时出现下述情况之一，则可能是病态方程组. (1) 在列主元消去法中出现小主元； (2) 在计算过程中行或列几乎线性相关或三角分解中对角元出现近似零的元素； (3) 矩阵A的元素数量级相差很大且无规律； (4) 剩余 很小，而解 很大，又达不