五线性代数方程组的直接解法.pptVIP

* §5 直接法的误差分析/*Error Analysis of Direct Method */ 一、扰动方程组的误差界/*Error Bound of Perturbed Systems */ 由实际问题得到的方程组的系数矩阵或者常数向量的元素，本身会存在一定的误差；这些初始数据的误差在计算过程中就会向前传播，从而影响到方程组的解。 初始数据误差和方程组的近似解的误差之间关系 例12 考察方程组： 精确解为 设方程组存在扰动 精确解为 上例说明该方程组的解对初始元素的扰动非常敏感。 设方程组为 系数矩阵 和常数向量 的扰动分别记为： 和 实际求解的方程组为 如果 和 很小，而 很大，则称 方程组 为病态(/*ill-conditioned*/)方程组， 称系数矩阵 为关于求解方程组或求逆的病态矩阵； 反之，如果 和 微小时， 也微小，则称 方程组 为良态(/*Well-conditioned*/)方程组， 称系数矩阵 为关于求解方程组或求逆的良态矩阵。 病态方程组对任何算法都将产生数值不稳定性 设 为可逆阵， 为一种从属矩阵范 数，则称 为矩阵 的条件数 ?若矩阵范数取2-范数，则得到谱条件数： ?若矩阵范数取1-范数，则得到1-条件数： ?若矩阵范数取 -范数，则得到 -条件数 ?如果 是正交矩阵，则 ?条件数的性质（ ） 设 ? 正交矩阵： ?如果 为正交矩阵，则 Lemma 设 ，在正交变换下，谱范数 和F-范数 保持不变，即设 则有： 证明： （1） 同理可证 （2） 对 设 ?设 和 是 按模最大和最小的特征值，则 取 2-范数等号成立 设 为可逆阵， 和 分别满 足方程组 和 其中 ，且 适当小，使得 则 其中所用的是任意向量范数和从属于它的矩阵范数 证明： 设 ， 如果 则 可逆，且 且 由条件 所以 可逆，且 （摄动定理） 方程组 有唯一解 Corollary(推论) ?如果 ，则 ?如果 ，则 在上述定理的条件下， 仅由 或者 引起的解的相对误差限可大致认为是 原始扰动的相对误差的倍数，倍数近似等于条件数。 ?如果 ，则 证明： 结论得证 实际的相对误差 与条件数的关系 例如： Hilbert矩阵就是一个著名的病态矩阵 3.54E+13 3.39E+10 2.9E+7 28375 1.60E+13 1.53E+10 1. 5E+7 15514 3.54E+13 3.39E+10 2.9E+7 28375 10 8 6 4 n 对称正定矩阵 *