用追赶法解三对角线方程组.doc

云南大学数学与统计学实验教学中心实验报告课程名称：数值计算方法学期： 20112012学年第一学期成绩：指导教师： 学生姓名：学生学号：实验名称：用追赶法解三对角线方程组实验编号：No. 3实验日期： 2011/10/29实验学时： 3学院： 数学与统计学院专业： 数学与应用数学年级： 2010级一、实验目的加强追赶法在解方程组中的应用，掌握多种不同的解线性方程组的解法以及编程的应用2、 实验内容用追赶法解下列严格对角优势的三对角线方程组三、实验环境C语言，Turboc四.实验方法追赶法；其算法如下：（1） 计算： （i=2,3,4n-1）（2） 解 (i=2,3,4,n) （3） 解 （i=n-1,n-2,2,1） 五、实验过程1实验步骤 (1）编程: 根据所用算法及选用语言编出源程序（2）. 开机, 打开所用语言系统输入所编源程序.（3）. 调试程序, 修改错误直至能正确运行.（4）. 运行程序并输出计算结果.2 关键代码及其解释#include #include #include /*预处理指令*/#defineMAX_n100 /*宏定义，即最多可输入100阶方阵*/#define PRECISION 0.000001 /*精度*/void SulutionOutput(float x,int n) /*输出方程组的解*/ /*带分号的为声明*/void TriDiagonalMatrixInput(float a,float b,float c,float f,int n) /*三对角元素的输入*/int Z_G_method(float a,float b,float c,float f,int n) /*三对角方程组求解 追赶法 */void main() /*主函数*/int n;float aMAX_n,bMAX_n,cMAX_n, dMAX_n,fMAX_n; /*定义数组*/printf(nInput n=);scanf(%d,&n);TriDiagonalMatrixInput(a,b,c,d,f,n); /*以下三行均为调用函数声明*/Z_G_method(a,b,c,d,f,n);SulutionOutput(f,n);3 调试过程1.本打算用递归调用，但是对那个不太熟悉，就改为用循环了；2. 首先用dowhile循环，出现语法错误太多；无法调试；后改为了dowhile循环了；5.逐步调试修改程序，达到预期结果；六、实验总结1遇到的问题及解决过程（1）对C语言掌握不太好，编写程序时语法错误较多.解决过程：耐心修改，逐步排除！（2）首先对牛顿法掌握不太好.解决过程：认真看书，问同学。2产生的错误及原因分析（1）Floating point error: Domain. 原因分析： 编写程序时出现了语法错误 （2）编写程序时没去分大小写，造成多处语法错误；原因分析：都怪自己粗心马虎3体会和收获。无论干什么事都要认真，踏实！努力做好自己该做的。上课要认真听讲，专心学习，做作业才省力！各门课程之间相互联系，应搞好各科的学习。七：程序源代码：#include#include#include#defineMAX_n100#define PRECISION 0.000001/解输出void SulutionOutput(float x,int n)int i;for(i=1;i=n;+i)printf(nx%d=%f,i,xi);/三对角方程组元素输入void TriDiagonalMatrixInput(float a,float b,float c,float f,int n)int i;printf(Input b1,c1,f1:);scanf(%f %f %f,&b1,&c1,&f1);for(i=2;in;+i)printf(Input a%d,b%d,c%d,f%d:,i,i,i,i);scanf(%f %f %f %f,&ai,&bi,&ci,&fi);printf(Input a%d,b%d,f%d,f%d:,n,n,n,n);scanf(%f %f %f%f,&an,&bn,&cn,&fn);/三对角方程组求解 追赶法int Z_G_method(float a,float b,float c,float f,int n)int i;c1/=b1;for(i=2;in;+i)ci/=(bi-ai*ci-1);f1/=b1;for(i=2;i0;-i)fi-=ci*fi+1;void main()int n;float aMAX_n,bMAX_n,cMAX_n, dMAX_n,fMAX_n;printf(nInput n=);scanf(%d,&n);TriDiagonalMatrixInput(a,b,c,d,f,n);Z_G_method(a,b,c,d,f,n);SulutionOutput(f,n);运行结果：Input n=5Input b1,c1,f1:4 -1 100Input a2,b2,c2,f2:-1 4 -1 200Input a3,b3,c3,f3:-1 4 -1 200Input a4,