追赶法求解三对角方程及其算例(共4页).doc

精选优质文档-倾情为你奉上追赶法求解三对角方程组要求：对于给定的三对角系数矩阵和右端项，可以求解线性代数方程组一、 追赶法的数学理论设系数矩阵为三对角矩阵则方程组Ax=f称为三对角方程组。设矩阵A非奇异，A有Crout分解A=LU，其中L为下三角矩阵，U为单位上三角矩阵，记可先依次求出L，U中的元素后，令Ux=y，先求解下三角方程组Ly=f得出y，再求解上三角方程组Ux=y。事实上，求解三对角方程组的2追赶法将矩阵三角分解的计算与求解两个三角方程组的计算放在一起，使算法更为紧凑。其计算公式为：(*)二、 追赶法的算法和流程图1.预处理生成方程组的系数及其除数，事实上，按式(*)可交替生成与：其计算公式为2.追的过程顺序生成方程组右端：据式(*)的计算公式为 3.赶的过程逆序得出方程组的解:其计算公式按式为三、 追赶法的Matlab实现function x=chase(a,b,c,f)%求解线性方程组Ax=f,其