jacobi迭代法和Gauss-Seidel迭代法

1、.數(shù)值計(jì)算方法實(shí)驗(yàn)報(bào)告（五）班級(jí)： 地信10801 序號(hào)： 姓名：一、 實(shí)驗(yàn)題目：jacobi迭代法和Gauss-Seidel迭代法二、 實(shí)驗(yàn)學(xué)時(shí): 2學(xué)時(shí)三、 實(shí)驗(yàn)?zāi)康暮鸵螅? 掌握迭代法的基礎(chǔ)原理。2 掌握jacobi迭代法和Gauss-Seidel迭代法的步驟。3 能用程序語言對(duì)jacobi迭代法和Gauss-Seidel迭代法進(jìn)行編程實(shí)現(xiàn)。四、實(shí)驗(yàn)過程代碼及結(jié)果1、代碼：#include<iostream.h>#include<math.h>float x100,xk100;float e;int N,M=1000;float a100101;void ini

2、tdata()cout<<"輸入方程階數(shù)："cin>>N;cout<<"輸入誤差限e："cin>>e;cout<<"輸入方程系數(shù)："<<endl;for(int i=1;i<=N;i+)for(int j=1;j<=N+1;j+)cin>>aij;cout<<"輸入初始解向量x0:"<<endl;for(i=1;i<=N;i+)cin>>xki;void jocobi()int

3、Nx=0,times=0;while(Nx<N)times+; Nx=0;if(times>=M)cout<<"發(fā)散"<<endl; break;for(int i=1;i<=N;i+)float sum=0;for(int j=1;j<=N;j+)if(i!=j)sum+=xkj*aij;xi=(aiN+1-sum)/aii;if(fabs(xi-xki)<e)Nx+;for(i=1;i<=N;i+)xki=xi;cout<<"times="<<times<<

4、;endl;for(int i=1;i<=N;i+)cout<<"x"<<i<<"="<<xi<<endl;void guass_seidel()int Nx=0,times=0;while(Nx<N)times+;Nx=0;if(times>=M)cout<<"發(fā)散"<<endl;break;for(int i=1;i<=N;i+)float sum1=0;float sum2=0;for( int j=i+1;j<=N

5、;j+)sum1+=xkj*aij;for( j=1;j<=i-1;j+)sum2+=aij*xkj;xi=(aiN+1-sum1-sum2)/aii;if(fabs(xi-xki)<e)Nx+;for(i=1;i<=N;i+)xki=xi;cout<<"times="<<times<<endl;for(int i=1;i<=N;i+)cout<<"x"<<i<<"="<<xi<<endl;void main()char ch;initdata();cout<<"請(qǐng)選擇解方程的方法：n"cout<<"A:jocobi B:guass_seidel n"ci