jacobi迭代法和Gauss-Seidel迭代法.doc

数值计算方法实验报告（五） 班级： 地信10801 序号： 姓名： 实验题目：jacobi迭代法和Gauss-Seidel迭代法 实验学时: 2学时 实验目的和要求： 掌握迭代法的基础原理。 掌握jacobi迭代法和Gauss-Seidel迭代法的步骤。 能用程序语言对jacobi迭代法和Gauss-Seidel迭代法进行编程实现。 四、实验过程代码及结果 1、代码： #includeiostream.h #includemath.h float x[100],xk[100]; float e; int N,M=1000; float a[100][101]; void initdata() { cout输入方程阶数：; cinN; cout输入误差限e：; cine; cout输入方程系数：endl; for(int i=1;i=N;i++) for(int j=1;j=N+1;j++) cina[i][j]; cout输入初始解向量x0:endl; for(i=1;i=N;i++)cinxk[i]; } void jocobi() { int Nx=0,times=0; while(NxN) { times++; Nx=0; if(times=M){cout发散endl; break;} for(int i=1;i=N;i++) { float sum=0; for(int j=1;j=N;j++) if(i!=j)sum+=xk[j]*a[i][j]; x[i]=(a[i][N+1]-sum)/a[i][i]; if(fabs(x[i]-xk[i])e)Nx++; } for(i=1;i=N;i++)xk[i]=x[i]; } couttimes=timesendl; for(int i=1;i=N;i++) coutx[i]=x[i]endl; } void guass_seidel() { int Nx=0,times=0; while(NxN) { times++; Nx=0; if(times=M) { cout发散endl; break; } for(int i=1;i=N;i++) { float sum1=0; float sum2=0; for( int j=i+1;j=N;j++) {sum1+=xk[j]*a[i][j];} for( j=1;j=i-1;j++) {sum2+=a[i][j]*xk[j];} x[i]=(a[i][N+1]-sum1-sum2)/a[i][i]; if(fabs(x[i]-xk[i])e) Nx++; } for(i=1;i=N;i++)xk[i]=x[i]; } couttimes=timesendl; for(int i=1;i=N;i++) coutx[i]=x[i]endl; } void main() { char ch; initdata(); cout请选择解方程的方法：\n; coutA:jocobi B:guass_seidel \n; cinch; if(ch==A) jocobi(); else if(ch==B) guass_seidel(); } 2．结果： 3