数值计算方法解三对角方程组的追赶法实 验 报 告4.18.docVIP

实 验 报 告 实验课程名称 数值计算方法 实验项目名称 解三对角方程组的追赶法 年 级 09 专 业 信息与计算方法 学生姓名 学 号 理 学 院 实验时间： 2011 年 4 月 18 日 学生所在学院：理学院 专业：信息与计算科学 班级：091 姓 名 学 号 实验组 实验时间 2011.4.18 指导教师 刘轶中 成 绩 实验项目名称 解三对角方程组的追赶法 实验目的及要求： 由于系数矩阵只有三条对角线上有元素，其余元素为零，所以算法实现只用4个一维数组，求学生掌握算法设计技巧，画出流程图，在计算机上完成整个算法。 实验硬件及软件平台：PC机，tc2.0，Internet网 实验步骤： 根据算法事先写出相应程序。 启动PC机，进入tc集成环境，输入代码。 编译调试。 调试通过，计算出正确结果后。 实验内容（包括实验具体内容、算法分析、源代码等等）：#includestdio.h 用追赶法解三对角方程组Ax=b， #includestdio.h #includeconio.h #includemath.h #define MAX_n 100 #define PRECISION 0.00001 void SulutionOutput(float x[],int n) { int i; for(i=1;i=n;i++) printf(\nx[%d]=%f,i,x[i]); } void TriDiagonalMatrixInput( float a[],float b[],float c[],float f[],int n ) { int i; printf(Input b[1],c[1],f[1]:); scanf(%f%f%f,b[1],c[1],f[1]); for(i=2;in;++i) { printf(Input a[%d],b[%d],c[%d],f[%d]:\n,i,i,i,i); scanf(\n%f%f%f%f,a[i],b[i],c[i],f[i]); } { printf(Input a[%d],b[%d],f[%d]:\n,n,n,n); scanf(\n%f%f%f,a[n],b[n],f[n]); } } int Z_G_method (float a[],float b[],float c[],float f[],int n) { int i ; c[1]/=b[1]; for(i=2;in;++i) c[i]/=(b[i]-a[i]*c[i-1]); f[1]/=b[1]; for(i=2;i=n;++i) f[i]=(f[i]-a[i]*f[i-1])/(b[i]-a[i]*c[i-1]); for(i=n-1;i0;--i) f[i]-=c[i]*f[i+1]; return f[i]; } void main() { int n; float a[MAX_n],b[MAX_n]; float c[MAX_n],f[MAX_n]; printf (\nInput n=); scanf(%d,n); TriDiagonalMatrixInput(a,b,c,f,n); Z_G_method(a,b,c,f,n); SulutionOutput(f,n); } 实验结果与讨论： 指导教师意见： 签名： 年 月 日