数值分析作业-matlab上机作课件.docVIP

数值分析 ———Matlab上机作业 学院： 班级： 老师： 姓名： 学号： 第二章 解线性方程组的直接解法 第14题 【解】 1、编写一个追赶法的函数 输入a，b，c，d输出结果x，均为数组形式 function x=Zhuiganfa(a,b,c,d) %首先说明：追赶法是适用于三对角矩阵的线性方程组求解的方法，并不适用于其他类型矩阵。 %定义三对角矩阵A的各组成单元。方程为Ax=d %b为A的对角线元素(1~n)，a为-1对角线元素(2~n)，c为+1对角线元素(1~n-1)。 % A=[2 -1 0 0 % -1 3 -2 0 % 0 -2 4 -3 % 0 0 -3 5] % a=[-1 -2 -3];c=[-1 -2 -3];b=[2 3 4 5];d=[6 1 -2 1]; n=length(b); u(1)=b(1); y(1)=d(1); for i=2:n l(i)=a(i-1)/u(i-1);%先求l(i) u(i)=b(i)-c(i-1)*l(i);%再求u(i) %A=LU,Ax=LUx=d,y=Ux, %Ly=d,由于L是下三角矩阵，对角线均为1，所以可求y(i) y(i)=d(i)-l(i)*y(i-1); e