[理学]第8章线性方程组的直接解法.ppt

西安电子科技大学理学院 主讲: 王卫卫 由实际问题建立起来的线性方程组Ax=b本身存在模型误差和观测误差，或者是由计算得到的，存在舍入误差等。总之，A,b都会有一定扰动?A ,?b, 因此实际处理的是A+?A或b+?b ，我们需要分析A或b的扰动对解的影响。 西安电子科技大学理学院 主讲: 王卫卫 x — exact solution x+?x — exact solution 当方程组的系数矩阵A和齐次项b受到扰动?A ,?b后，其解x会受到怎样的扰动？分三种情况讨论 1. b有扰动?b, A无扰动(?A=0) 解的相对误差 齐次项b的相对误差 当b受到扰动?b时，引起的解的相对误差不超过b的相对误差的??A??·??A-1??倍 西安电子科技大学理学院 主讲: 王卫卫 2. A有扰动?A, b无扰动(?b =0) 解的相对误差 A的相对误差 设??A-1 ?? ?? ?A??1 Th1 A+?A=A(I+ A-1?A)可逆 ??A-1?A??1 ?? ?A??充分小 西安电子科技大学理学院 主讲: 王卫卫 3. A有扰动?A, b有扰动?b,类似可得 解的相对误差 A的相对误差 b的相对误差 当A受到扰动?A时，引起的解的相对误差不超过A的相对误差的??A??·??A-1??倍 当方程组的系数矩阵A或齐次项b受到扰动?A ,?b时，引起的解的相对误差完全由??A??·??A-1??来决定，它刻画了方程组的解对原始数据的敏感程度。 西安电子科技大学理学院 主讲: 王卫卫 Def 7 矩阵的条件数 A非奇异， ??A??p·??A-1??p =cond(A)p称为矩阵的条件数(p=1,2,?) cond(A)p相对较大，称方程组是病态的；否则，若cond(A)p相对较小，称方程组是良态的。 常用的矩阵条件数 1. ??A??? ·??A-1??? =cond(A)? 西安电子科技大学理学院 主讲: 王卫卫 条件数的性质 西安电子科技大学理学院 主讲: 王卫卫 西安电子科技大学理学院 主讲: 王卫卫 example Hilbert matrix 西安电子科技大学理学院 主讲: 王卫卫 How to verify whether a given linear system is ill posed? Cond(A)p 1(不易计算) 系数矩阵A 的三角约化中出现小主元 |?max(A)|/| ?min(A)|1 Det(A) 很小，或A的某些行近似线性相关 A的元素间数量级相差很大，且无一定规则。 西安电子科技大学理学院 主讲: 王卫卫 三角方程组，易于求解 改进的平方根法计算量与平方根法类似，但不需开方，只需存储下三角部分 西安电子科技大学理学院 主讲: 王卫卫 三对角方程组 追赶法—求解三对角方程组 三对角矩阵 西安电子科技大学理学院 主讲: 王卫卫 Theorem 4 若三对角矩阵A为弱对角占优，即满足 则它可唯一分解成单位下二对角阵L和上二对角阵U的乘积 追赶法的代数基础 西安电子科技大学理学院 主讲: 王卫卫 三角方程组，易于求解 追 赶 追赶法的求解过程就是将系数矩阵分解两个简单的二对角矩阵，从而归结为求解两个简单方程组的过程。追赶法的原理和高斯消去法相同，但考虑到方程组的特点，计算时会把大量零元素撇开，从而大大节省计算量：5n-4次乘除法.用4个一维数组存放ai, bi , ci , di , li 占ai , ui 占bi , yi 占di , xi 占yi 西安电子科技大学理学院 主讲: 王卫卫 summary Cramer rule Gauss elimination LU factorization Gauss-Jordan elimination Square root /improved square root 追赶法 (n+1)! n3/3 n3/3 n3/2 n3/6 5n-4 Row/column/complete Pivoting Only eliminating elements in a column below the diagonal one No pivoting Directly factorization column pivoting Eliminating elements in a row except for the diagonal element A symmetric and positive definite No pivoting A三对角，弱