高斯-赛德尔算法.docVIP

数值计算方法实验报告 班级：09医软1班 指导教师：金力 工作组成员：谢清泉（负责程序主体部分的编写） 魏月圆（负责流程图和实验报告的书写） 项波波（负责程序的调试和修改） 王兆德（负责材料的收集和整理） 一 实验名称 高斯—塞德尔迭代法 二 实验目的 理解高斯塞德尔算法的基本思想，及公式的推导过程，会用此公式解简单线性方程组。 三 实验原理 在雅可比迭代中，总是用前次近似分量x1（k），x2（k），.. ….xi-1(k),xi(k),……xn(k)去计算当前分量xi（k+1）（i=1,2，……,n）。实际上此刻前i-1个分量x1,x2,…..xi-1的新近似值x1（k+1），x2(k+1),……xi-1(k+1)代替x1(k),x2(k),…… Xi-1(k)去计算xi（k+1）可能会得到更满意的效果。据此得到迭代公式就是高斯—塞德尔公式，即 Xi（k+1）=（bi - - ）/aii (i=1,2,……n) 四 流程图 K++开始输入线性方程行数n，输入系数矩阵a，常数矩阵b，精度mK=1i=0 K++ 开始 输入线性方程行数n，输入系数矩阵a，常数矩阵b，精度m K=1 i=0 ｜xi+1-xi︳ ｜xi+1-xi︳ m i++ Xi（k+1）=（bi - - ）/aii Xi in Y N Y NY输出结果结束输出此次迭代失败K10001 N Y 输出结果 结束 输出此次迭代失败 K10001 N N 五 程序代码 #include string.h #include math.h #define N 100 float Table(int n,float a[N][N],float b[N]) /*定义数组并分配值*/ { int i,j; printf(请输入系数矩阵a：\n); label:for(i=0;in;i++) { printf(行 %d:,i+1); for(j=0;jn;j++) scanf(%f,a[i][j]); /*二维数组分配值*/ } printf(请输入常量向量 b:); for(i=0;in;i++) scanf(%f,b[i]); /*一维数组分配值*/ printf(联合矩阵为:\n); for(i=0;in;i++) /*将数组结合，检验输入是否有误*/ { for(j=0;jn;j++) printf(%10.5f,a[i][j]); printf(%10.5f,b[i]); printf(\n); } } float init_vec(int n,float x[N]) /*获取初始值*/ { int i; printf(请输入初始解 x:); for(i=0;in;i++) scanf(%f,x[i]); } float gs(int n,float a[N][N],float b[N],float x[N]) /*迭代的实现*/ { int i,j,k; float m; float tmp1,tmp2,x2[N]; printf(请输入精度m： ); scanf(%f,m); for(k=0;k10001;k++) /*设定最大迭代次数不能超过1000*/ { for(i=0;in;i++) x2[i]=x[i]; for(i=0;in;i++) { tmp1=0.0; tmp2=0.0; for(j=0;ji;j++) tmp1+=a[i][j]*x[j]; for(j=i+1;jn;j++) tmp2+=a[i][j]*x2[j]; x[i]=(b[i]-tmp1-tmp2)/a[i][i]; } for(i=0,