实验用追赶法求解三对角方程组.docVIP

班级： 化院研07-2班 姓名：陶凯 学号：实验5.3　　用追赶法求解三对角方程组 实验目的 会使用Matlab语言编程使用追赶法求解线性方程组。 实验原理 设系数矩阵为三对角矩阵 　　　　 则方程组Ax=f称为三对角方程组。 设矩阵A非奇异，A有Crout分解A=LU，其中L为下三角矩阵，U为单位上三角矩阵，记 　　　 可先依次求出L，U中的元素后，令Ux=y，先求解下三角方程组Ly=f得出y，再求解上三角方程组Ux=y。 事实上，求解三对角方程组的2追赶法将矩阵三角分解的计算与求解两个三角方程组的计算放在一起，使算法吏为紧凑。其计算公式为： 实验程序 function x=chase (a,b,c,f) %定义函数chase n=length(b); if n-1==length(a) for i=n-1:-1:1 a(i+1)=a(i); end end %将a设置为n维向量 c(1)=c(1)/b(1); f(1)=f(1)/b(1); for i=2:n-1 b(i)=b(i)-a(i)*c(i-1); c(i)=c(i)/b(i); f(i)=(f(i)-a(i)*f(i-1))/b(i); end f(n)=(f(n)-a(n)*f(n-1))/(b(n)-a(n)*c(n-1)); for i=n-1:-1:1 f(i)=f(i)-c(i)*f(i+1); end x=f; 实验内容 设A =2 7 0 3 2 9 0 4 5 F=5 6 7 求解线性方程组Ax=F. 实验步骤 1 在matlab中输入上面的程序； 2 在commond窗口依次输入数据 a=[3 4] b=[2 2 5] c=[7 9] f=[5 6 7] 3 在commond窗口调用函数x=chase (a,b,c,f). 得到输出结果x = -0.6433 0.8981 0.6815即为原线性方程组的解。 将x代入验证A*x 得到 ans = （5 6 7）=f ，即结果正确，x是原线形方程组的解。 实验心得 1 要有耐心，编matlab时出现错误是常有的事。 2 数学并不枯燥，那种解决问题后的成就感会让你很兴奋不已。