三对角线性方程组的解法.pptVIP

目 录 一、实验分析‥ ‥ ‥ ‥ ‥ ‥ ‥ 3 二、实验内容‥ ‥ ‥ ‥ ‥ ‥ ‥ 4 三、实验要求‥ ‥ ‥ ‥ ‥ ‥ ‥ 5 四、实验过程‥ ‥ ‥ ‥ ‥ ‥ ‥ 6 五、结果分析‥ ‥ ‥ ‥ ‥ ‥ ‥ 19 六、程序代码‥ ‥ ‥ ‥ ‥ ‥ ‥ 21 七、课程信息‥ ‥ ‥ ‥ ‥ ‥ ‥ 30 一、实验分析 关于线性方程组的数值解法一般有两类：直接法和迭代法。而对于系数矩阵较特殊的一类线性方程组，解法也相对有所不同。在实际问题中，例如解常微分方程边值问题，解热传导问题等，都会要求解系数矩阵为对角占优的三对角方程组。 如何求解三对角方程组就成为关键。三对角线性方程组是系数矩阵为三对角矩阵的一类方程组，其数值解法有：追赶法，Jacobi迭代法，Gauss迭代法和SOR法。 二、实验内容 考虑线性方程组： 其中A为三对角矩阵，编制程序求解该方程组 三、实验要求 （1）考虑不同的数值解法，对方程组进行编程，编写其通用代码； （2）对不同的程序代码用实例进行验证，取矩阵： 方程初值 ： 。 （3）比较不同方法对求解三对角线性方程组的优异性。 四、实验过程 (一)追赶法 编程思想： 矩阵的三对角线性方程组是由矩阵的直接三角分解法来推到的来的，求解等价于解两个三角形方程组： （1） ，求y；（2） ，求x 从而得到解三对角线方程组的追赶法公式。 1. 计算的递推公式： （i=2,3,……n-1）； 2. 解 （i=2,3,……n）； 3. 解 （i=n-1，n-2，……，2,1）。 运行结果：（点击查看源程序） a=-1*ones(1,4); c=a; b=2*ones(1,5); f=[1 0 1 0 0]; [x]=ZhuiGan(a,b,c,f)‘ n =5 x = 1.33333333333333 1.66666666666667 2.00000000000000 1.33333333333333 0.66666666666667 注：其中a为对角下向量，b为对角向量，c为对角上向量，f为方程常 系数。x为方程组的解。 （二）迭代法 对于 ，A为非奇异矩阵，除了可用追赶法计算之外，也可建立迭代法求解。迭代法的一般格式为： ，其中B为迭代矩阵，对由迭代格式产生迭代序列 ，若 ，则 即 为 方程 的解。 迭代法有以下几种方法：Jacobi迭代法，Gauss迭代法和SOR法。 Jacobi迭代法 1.迭代过程： 对于 ，设A可逆，且对于 令D为A的对角元素部分，则 ,继而对于方程组有 ，所以有解 ，令 为迭代矩阵， 。则Jacobi迭代法的迭代格式为： 。 2.算法设计： （1）将方程组系数矩阵用 表示，常数项用 表示。 （2）用 存储初值； （3）用库函数diag,tril,triu求取对角矩阵D,下三角矩阵L,上三角矩阵U； （4）利用得到的D,L,U，求取迭代矩阵J，和f; （5）循环判断， 存取循环值，用二范式迭代循环，当 时，终止循环。 3.运行结果：（点击查看源程序）