实验2追赶法算法设计及MATLAB实现.doc

数 值 计算方法实验报告实验序号：实验二 实验名称：追赶法算法设计及MATLAB实现实 验 人：专业年级：教 学 班：学 号：实验时间：实验二　追赶法算法设计及MATLAB实现一、实验目的１．初步掌握算法设计规则；２．初步掌握MATLAB程序设计规则．二、实验内容1．构造利用追赶法求解三对角线性方程组的算法；2．在MATLAB环境下编写追赶法的程序〔函数；3．自由选择若干个三对角线性方程组求解三、实验步骤1．追赶法算法：算法名称：thomas输入参数：向量a,b,c,f输出参数：输出解信息x算法的自然语言：Step1：u=b,y=b;Step2：对于 i=2,3,….n;Step2.1：当u,否则转step5 l=a/u; u=b-l*c； y=f-l*y；Step3：当u时,x=y/u,否则转step5Step4：对于：i=n-1,n-2,…..,2,1,转step6 x=/uStep5：无解信息,转step7Step6：输出xStep7：关机2．MATLAB程序function [x,L,U]=thomasn=length;% 对A进行分解u<1>=b<1>;for i=2:n if~=0> l=a/u; u=b-l*c; else break; endendL=eye+diag;U=diag+diag;x=zeros;y=x;% 求解Ly=by<1>=f<1>;for i=2:n y=f-l*y;end% 求解Ux=yif~=0> x=y/u;endfor i=n-1:-1:1 x=-c*x>/u;end3．求解实例例1．方程组例2．方程组例3．方程组四、实验结论对于追赶法我最先写的是如下的程序：3 / 4.但是出现了如上截图中的错误,后来与同学讨论还是没能解决我的问题,最后借鉴了她的算法得到了正确的结果。