《雅可比迭代法》word版

《《雅可比迭代法》word版》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、雅可比迭代法班级：计算062姓名：王保翔3060811028目的意义：设方程组Ax=b的系数矩阵A非奇异而且，将A分裂为A=D+L+U,可以使计算简便。其中，思想：A=D+L+U，其中，将方程组乘以，得到等价的方程组，i=1，2，…n算法：Step1：给定一组x，即初值。Step2：用for循环计算：x[k+1]=(b[i]-)/a[k][k].Step3:当abs(x[k+1]-x[k])#definep0.1#defineN5voidmain(){inti,j,k,s,n;doublex[N

2、],a[N][N],b[N],q[N],c[N],d[N];doublet,u;printf("Inputn:");scanf("%d",&n);for(i=0;i<=n;i++){printf("Inputx[%d],b[%d]:",i,i);scanf("%lf,%lf",&x[i],&b[i]);}for(i=0;i<=n;i++)for(j=0;j<=n;j++){printf("Inputa[%d][%d]:",i,j);scanf("%lf",&a[i][j]);}for(i=0;i<=n;i++)for(j=0;j<=n;j++)pri

3、ntf("a=%lf",a[i][j]);for(k=0;;k++){t=0;for(s=0;s<=n;s++)d[s]=x[s];for(i=0;i<=n;i++){u=0;d[i]=x[i];for(j=0;j<=n;j++){if(j!=i)u=u+a[i][j]*d[j];elsecontinue;}printf("u=%lf",u);x[i]=(b[i]-u)/a[i][i];q[i]=x[i]-d[i];printf("x[%d]=%lf",i,x[i]);}for(i=0;i<=n;i++){if(q[i]>=0)c[i]=

4、q[i];elsec[i]=-q[i];t=t+c[i];}if(t