第五单元解线性方程组的直接方法.pptVIP

例 解: 在一些实际问题中常有解三对角线性方程组Ax=f的问题，即： 对于具有条件（2）的方程组（1），我们介绍下面的追赶法求解。追赶法具有计算量少，方法简单，算法稳定的特点。 定理：设有三对角线性方程组Ax=f，且A满足条件（2），则A为非奇异矩阵。 证明：用归纳法证明。显然n=2时，有： 现假设定理对n-1阶的满足条件（2）的三对角矩阵成立，求证对满足条件（2）的n阶三对角矩阵定理亦成立。由条件 ，则利用消元法的第1步有： 证明：由于A是满足条件（2）的n阶三对角阵。因此A的任一个顺序主子式 亦满足条件（2）的n阶三对角矩阵。由上一个定理 根据这一结论以及三角分解定理知，这种矩阵A可进行三角分解： A=LU。在这里特别的有： 上面已经验证(5)对i=1成立。现设(5)对i-1成立。求证(5)对i成立。 由假设 ，再由(4)及假设条件(2)有: 由(4)可知: 可见由A的假设条件,我们完全求出了 实现了对A的LU分解。从而求解Ax=f 等价于求解两个三角方程组: 由此可得求解三对角线性方程组的追赶法: 1.计算 的递推公式: 2.求解 Ly=f 的公式: 3.求解 Ux=y 的公式: 我们将计算 及 的过程称为追的过程。计算方程组的解 的过程称为赶的过程。追赶法求解 Ax=f 仅需5n-4次乘除运算。工作量较小。 在计算机上，只需用3个一维数组分别存储A的系数 以及两个一维数组保存计算的中间结果 和 或 . 例:用追赶法求解: 解：（1）计算 ： （2）计算 ： 对于追赶法，由于已证明: （3）计算 ： 因此追赶法中不会出现中间结果迅速增长和全入误差产生积累现象。 需要对 （n维 列向量空间）中向量及 中的矩阵引进某种度 量，即引进向量或矩阵范数的概念。范数是长度概念的推广 为了对方程组的计算解进行误差分析，为了讨论迭代法的收敛性， 定义1：（向量范数） 第k步计算如下： 对于k=1，2……..，n-1做到（4） （1）按列选主元，即确定 使 （3）如果 ，则交换 第 行第k行元素 （2）如果 =0，则A为奇异矩阵，停止计算。 （4）消元计算： * * 解：精确解为（舍入值）： 例：用列主元消去法去解方程组 计算解 在常数项内 得到 （5）回代计算： 本例是具有舍入的4位浮点数进行运算，所得的计算解还是比较准确的。 下面是完全主元素消去法框图（图6-1） 图（6-1） stop 输入 n,A,b,c I(Z) I (I=1,2,…….n) k=1,2…….n 是 输出 detA=0 否 交换行 是 否 转下页 消元计算 交换列且 接上页 否 Stop ① ② ③ ②— ① ③— ① ②* ③* 第二步，选4.12为列主元，不需换行， ③* — ②* ③** 由①②*③**回代求解得： 解：第一步，选5.00为主列元，将 对凋位置， ③ ③ ① ① ② 与原方程组准确解 比较，可知，本题用3位有效数字计算的列元素法是相当准确的。 大量实践表明：列主元法为解线性方程组的精确方法。 完全主元素消去法在选主元素时要花费多机器时间。下面介绍另一种常用的方法即列主元素消去法，它仅考虑依次按列选主元素，然后换行使之变到主元素位置上，再进行消元计算。设用列主元素消去法解Ax=b已完成k-1步计算。即有： 列主元素消去法步骤： 设Ax=b。对于具有行交换的列主元素消去法，消元结果