完整版本2.6追赶法.docxVIP

第 16 讲 追赶法、误差分析 在实际应用问题中， 经常会遇到解三对角线方程组。 例如： 用三次样条函数的插值问题中得 到的三转弯及三弯矩方程组， 当时说可用追赶法来求解。 还有用差分法解二阶线性常微分方 程边值问题， 若用三点插值格式也得到解三对角线方程组， 本节介绍该类方程组中的特例及 该种方程组的解法：追赶法。 优点： 1.计算量小。 2.方法简单，存贮量小。 3.数值稳定的（对舍入误差来说） 。 追赶法 三对角线方程组的一般表示方法： l1 a2 l2 L a3 l3 . . . . an l n b1 l1 ci li ui bi 1 ai 1ui li 1 l1 b1 ui ci / li li 1 bi 1 ai 1ui 可见，对 A 的分解只需求 l i , ui 且按 l1 u1 l 2 u2 ..... 的递推过程进行，形象地称为 “追 ”的过程 y1 f1 / l1 yi ( f i ai yi 1 ) / l (i 2,....n)  1 u1 1 u2 1 . U . . . un 1 l n (i 1,2,3,..., n 1) (i 1,2,3,..., n 1) l n 1 un 1 ln xn yn xi yi ui xi 1 (i n 1,.....2,1) 形象地称回代求解过程为 “赶 ”的过程 追赶法的计算量为 5n-4 次乘除法， 可用 4 个 一 维数组存放 ai , bi , ci , fi 。共占用 4n-2 个单元， 在计算过程中 li , ui , yi 依次覆盖掉 bi , ci , f i 最后， xi 覆盖掉 yi ，所以， 追赶法具有计算量小，占用内存单元少的特点。 2、误差分析 病态方程组与条件数 一个线性方程组 Ax=b 是由它的系数矩阵 A 和它的右端项 b 所确定， 在实际问题中， 由 于各种原因， A 或 b 往往有误差，从而使得解也产生误差。本节方程组的系数矩阵 A 或右 端项 b 的微小误差对解向量的影响问题。 定义 1 求解线性方 程 组 Ax b 时，若 A 或 b 有微小扰动 A 或 b 时，解 x 的扰 动 x 很大，则称此方程组为病态方程组，相应的系数矩阵 A 称为病态矩阵，反之，若此 时 x 很小，则称方程组为良态方程组，矩阵 A 称为良态矩阵。 定义 2 称 x ， A ， b 分别为解向量 x，矩阵 A 与右端向量 b 的相对扰动。 x A b 定义 3 设 A 为 n 阶可逆矩阵，则称数 Cond ( A) A 1 A 为矩阵 A 的条件数，其 中 ? 是矩阵的算子范数常用的条件数为： Cond ( A) A A Cond1 ( A) A 1 A Cond2 ( A) A 2 A  1 1 1 1 2 判断一个线性方程组 Ax=b 是否病态，需要算系数矩阵的条件数 Cond ( A) A A1， 而条件数的计算首先要计算逆矩阵的范数， 计算量很大。 在实际中常可通过求解过程直观地 判断方程组的病态特征。 (1)若在主元素消元过程中出现小 主元，则 A 可能是病态阵，但病态阵未必一定有这 种小主元。 (2)系数矩阵的行列式的值相对来说 很 小，则 A 有可能是病态阵。 (3)从矩阵本身来看， 若元素间数量级很大且无一定规律， 或者矩阵的某些行 (列 )近似线 性相关，这样的矩阵有可能是病态的. (4)如果 A 的最大特征值和最小特征值之比 (按绝对值 )较大，则 A 有可能是病态的。 用选主元的消去法不能解决病态问题，对于病态方程组可采用以下措施： 采用高精度的运算，减轻病态影响。 (2) 采用预处理方法改善 A 的条件数，例如可选择非奇异的对角阵或三角阵 P， Q，将 Ax=b 转化为等价形式  pAQy pb y Q 1 x