课程设计追赶法解三对角方程组.docVIP

追赶法解三对角方程组 摘 要 本文主要通过运用追赶法来求解三对角方程组的问题.追赶法是用来求解三对角方程组的专用方法,对于三对角方程组,追赶法比Gauss消去法的计算量要小的多,本文主要介绍了追赶法的原理,并用Matlab编写求解程序,以实现对三对角方程组的求解,进一步解决实际中的问题.并且根据所得出的结果分析追赶法算法和高斯消去的法的计算精度. 关键词：追赶法,三对角方程组,追赶法的Matlab程序 目 录 1.前 言 1 2.解题思想和方法 1 2.1 追赶法解题思想 1 2.2追赶法解题原理 4 3.对追赶法的MATLAB求解 4 3.1实验程序 4 3.2 应用举例 5 4.与高斯消去法的精度比较 6 课程设计总结 8 参考文献 8 前 言 当今很多科学与工程计算问题大都可以化为线性方程组的形式,所以有效的求解线性方程组在科学和工程计算中是非常重要的.虽然线性代数方程的求解方法和数值计算软件包均很成熟,但随着并行计算机的发展,问题的求解速度和解题规模都大大提高,因而使数值计算方法和响应的数学软件包都产生了变化,相应的线性方程组的有效并行求解也引起了人们的普遍关注. 追赶法是用来求解三对角方程组的专用方法,生活中很多实际问题,都归结为求解线性方程组.例如解常微分方程边值问题,解热传导方程以及船体数学放样中建立三次样条函数等,都会要求就解系数矩阵成三对角线性的方程组.而解三对角方程组的最简单方法是追赶法,公式简单,计算量小,所占用的存储单元少,所以在小机器上也能求解.追赶大事用来求解三对角方程组的专用方法,对于三对角方程组,追赶法比Gauss消去法的计算量要小得多.应用追赶法求解三对角线性方程,追赶法仍然保持LU分解特性,它是一种特殊的LU分解.充分利用了系数矩阵的特点,而且使之分解更简单,得到对三对角线性方程组的快速解法. 本文讨论使用追赶法解线性方程组.介绍追赶法的理论,求解线性方程组的追赶法的实现以及追赶法的应用. 2.解题思想和方法 三对角矩阵是一种具有特殊意义的带状矩阵.用差分法求解二阶常微分方程边值问题时,最后常规解为求解具有三对角系数矩阵的线性方程组.对三对角矩阵实行Doolittle（或Crout）分解,便得到求解三对角方程组的最有效方法---追赶法. 设矩阵A非奇异,A有Crout分解A=LU,其中ＬＵUx=y,先求解下三角方程组Ly=f得出y,再求解上三角方程组Ux=y,事实上求解三对角方程组的2追赶法将矩阵三角分解的计算与求解两个三角方程组的计算放在一起,使算法更为紧凑. 本文主要介绍追赶法的原理,以及编写Matlab程序实现在计算机上的应用并分析他们的计算精度,比较解三对角方程组的最优解的问题. 2.1 追赶法解题思想 在实际问题中,经常遇到以下形式的方程组 （2.11） 这种方程组的系数矩阵A为三对角矩阵,即 以下针对这种方程组的特点提供一种简便有效的算法—追赶法.追赶法实际上是高斯消去法的一种简化形式,它同样分消元与回代两个过程. 先将（2.11）第一个方程中x1的系数化为1 （2.12） 有注意到剩下的方程中,实际上只有第二个方程中含有变量x1,因此消元手续可以简化.利用（2.12）可将第二个方程化为,这样一步一步地顺序加工（2.11）的每个方程,设第k – 1个方程已经变成 （2.13） 再利用（2.13）从第k个方程中消去xk-1,得： 同除,得 记 则有 这样做n – 1步以后,便得到： 将上式与（2.11）中第11个方程联立,即可解出xn = yn ,这里 于是,通过消元过程,所给方程组（2.11）可归结为以下更为简单的形式： （2.14） 这种方程组称作二对角型方程组,其系数矩阵中的非零元素集中分步在主对角线和一条次主对角线上 对加工得到的方程组（2.14）自下而上逐步回代,即可依次求出xn,xn-1,…,x1,计算公式为： （2.15） 上述算法就是追赶法,它的消元过程与回代过程分别称作“追”过程与“赶”过程.综合追与赶的过程,得如下计算公式： （2.16） （2.17） 2.2追赶法解题原理 (1) 计算的递推公式 ,. (2) 解Ly=f , (3) 解Ux=y . 我们将计算系数及的过程称为追的过程.将计算方程组的解的过程称为追赶的过程. 3.对追赶法的MATLAB求解 3.1实验程序 function x=chase(a,b,c,f) %定义函数chase n=length(b