实验七G-S迭.docVIP

数学与软件科学学院 实验报告 学期： 至 第 学期 年 月 日 课程名称：计算机数值计算方法 专业： 年级 班 实验编号：7 实验项目： Gauss-Seidel迭代格式 指导教师：张莉 姓名： 学号： 实验成绩： 一、实验目的及要求 实验目的： (1) 掌握Gauss-Seidel迭代算法的本质，比较其与Jacobi迭代格式的区别。 (2) 用一门自己熟悉的语言实现此算法。 实验要求： (1) 根据算法的特征写出算法流程图； (2) 测试给出的例子； (3) 分析什么时候此算法是收敛的。 二、实验内容 1.用高斯-赛德尔方法解方程组：P136例6.2。 2. 用高斯-赛德尔方法解方程组 自动取初值，当 三、实验步骤(该部分不够填写.请填写附页) #includestdio.h #includemath.h #define MAX_N 20 #define MAXREPT 100 #define epsilon 0.00001 int main() { int n; int i,j,k; double err; double xmax,oxmax; static double a[MAX_N][MAX_N]; static double l[MAX_N][MAX_N],u[MAX_N][MAX_N]; static double x[MAX_N],nx[MAX_N]; printf(\nInput n value(dim of AX=C):); scanf(%d,n); if(nMAX_N) { printf(The input n is larger than MAX_N,please redefine the MAX_N.\n); return 1; } if(n=0) { printf(Please input a number between 1 and %d.\n,MAX_N); return 1; } printf(Now input the matrix a(j,i),i,j=0,...,%d:\n,n-1); for(i=0;in;i++) for(j=0;jn;j++) scanf(%lf,a[i][j]); for(i=0;in;i++) x[i]=1; oxmax=0; for(i=0;iMAXREPT;i++) { for(j=0;jn;j++) { nx[j]=0; for(k=0;kn;k++) nx[j]+=a[j][k]*x[k]; } xmax=0.0; for(j=0;jn;j++) if(fabs(nx[j]xmax)) xmax=fabs(nx[j]); for(j=0;jn;j++) nx[j]/=xmax; for(j=0;jn;j++) x[j]=nx[j]; if(fabs(xmax - oxmax) epsilon) { printf(Solve... max lamda=%lf\n,xmax); printf(The vector is :\n); for(i=0;in;i++) printf(%lf\n,nx[i]); break; /* return 0; */ } oxmax=xmax; } printf(After %d repeat ,max no result...\n,MAXREPT); for(i=0;in;i++) u[i][i]=1; for(k=0;kn;k++) { for(i=k;in;i++) { l[i][k]=a[i][k]; for(j=0;j=k-1;j++) l[i][k]-=(l[i][j]*u[j][k]); } for(j=k+1;jn;j++) { u[k][j]=a[k][j]; for (i=0;i=k-1;i++) u[k][j]-=(l[k][i]*u[i][j]); u[k][j]/=l[k][k]; } } for(i=0;in;i++)