解三对交线方程组追赶法.pptVIP

* §7 解三对交线方程组的追赶法 必要性（即为什么讲追赶法）： 在实际应用问题中，经常会遇到解三对角线方程组。例如：用三次样条函数的插值问题中得到的三转弯（8.9、8.10或8.11）及三弯矩方程组（8.19或8.20），当时说可用追赶法来求解。还有用差分法解二阶线性常微分方程边值问题，若用三点插值格式也得到解三对角线方程组，本节介绍该类方程组中的特例及该种方程组的解法：追赶法。 优点： 1.计算量小（仅5n-4次乘除法运算）。 2.方法简单，存贮量小。 3.数值稳定的（对舍入误差来说）。 一、理论基础 三对角线方程组的一般表示方法： 或 其中，A满足条件： 以下研究在（7.2）成立的条件下（7.1）的解法。 定理15 满足（7.2）方程组(7.1) 的系数矩阵A为非奇异矩阵。 证明：用归纳法证明，显然对 时有 假设定理对 成立，下证对 定理亦成立。 ，则由高斯消去法： 显然， 下验证B为 阶三对角阵且满足 （7.2）: 则由归纳法假设 故有 ，即A为非奇异矩阵。 定理16 设 其中A为满足（7.2）的三对角阵，则A的所有顺序主子式都不为零，即 证明 由A是满足（7. 2）的n阶三对角阵，因此A的任一个顺序主子阵 亦是满足（7. 2）的k阶三对角阵，由定理15，则有 二、追赶法 设A为满足（7.2）n阶三对角阵 , 则有唯一三角分解A=LU A=LU ——杜里特而分解——一般三角分解； A=LU——Cront分解——追赶分解法； A=LDR——平方根法。 若L为下三角阵，U为单位上三角阵。 即有： 由矩阵乘法， 得: （1） （2） （3） 于是，得到解（7.1）的追赶法计算公式。 （1）分解计算公式： 次数: 1次 2(n-2)次 （2）求解 逆推公式 1次 3(n-1)次 （3）求解 逆推公式 n-1次 将计算 的过程称为追的过程，计算方程组解 过程称为赶的过程。乘法的次数为6n-6次，若令 ，则计算量为(6n-6)-(n-2)次乘除法，即5n-4次,但此时需贮存 ，尽管计算量小，但增加了计算机的内存，则两者各有利弊。这实际上也是一种计算方法好坏的衡量标准：两者兼顾。 稳定性可由定理17反映： 定理17 设有三对角线方程组 ，其中A满足（7.2）式，则由追赶法计算公式得到 满足： （1） （3） （2） 分析: (1)要证 ，只要证 ,而此式中含有 ，因此用归纳法证明。(2)、（3）只要用三角不等式即可证得。 于是，解 化为解两个方程组： 注： 实际上，条件（7.2）是充分的，但条件并非完全必要，三对角线上不能有零元素，条件太苛刻，于是条件可以作适当改变，使追赶法还可行 。 *