实验4 Jcobi迭代法和GS迭代2.doc

《实验4 Jcobi迭代法和GS迭代2.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验题目1．Jacobi迭代法用Jacobi迭代法求解线性方程组，保留四位有效数字(err＝1e-4)，其中A=[8-11;210-1;11-5];b=[1;4;3]。2．Gauss-Seidel迭代法用Gauss-Seidel迭代法求解线性方程组，保留四位有效数字(err＝1e-4)，其中A=[8-11;210-1;11-5];b=[1;4;3]。3．分别用Jacobi迭代法和Gauss-Seidel迭代法求解方程组要求精度，初始，最大迭代次数N=25，试比较这几种迭代法的迭代次数和收敛速度。1.程序：#include#includeintma

2、in(){intn=3,i,k,j,mm=1000;//最大迭代次数;floatt,x[mm][n],dx[n],dx_max=1,err=1e-4;//精度for(i=0;i

3、d]=",n);for(k=0;dx_max>1e-8;k++){dx_max=0;for(i=0;i

4、20.22620.3005-0.487030.22340.3061-0.494740.22510.3058-0.494150.22500.3056-0.493860.22490.3056-0.493970.22490.3056-0.493982.程序：#include#includeintmain(){intn=3,i,k,j,mm=1000;//最大迭代次数;floatt,x[mm][n],dx[n],dx_max=1,err=1e-4;//精度for(i=0;i

5、0,-1,1,1,-5};floatb[3]={1,4,3};for(i=0;i1e-8;k++){dx_max=0;for(i=0;i

6、else{t=0;for(j=0;j<=i-1;j++)t+=a[i][j]*x[k+1][j];for(j=i+1;j

7、-0.49390.22490.3056-0.49393.Jacobi迭代法程序：#include#includeintmain(){intn=6,i,k,j,mm=1000;//最大迭代次数;floatt,x[mm][n],dx[n],dx_max=1,err=1e-5;//精度for(i=0;i