《雅克比迭代法》由会员分享，可在线阅读，更多相关《雅克比迭代法（3页珍藏版）》请在金锄头文库上搜索。
1、 雅克比迭代法上机题目： 用雅克比迭代法解线性方程组上机程序：#include#includemain()double A33=5,2,1,-1,4,2,2,-3,10,b3=-12,20,3; / 输入系数矩阵A 和 右端向量bdouble n=3, tol=1.0e-3,x3=-3,1,1; / 输入方程大小n，误差限tol， 和初始向量 xdouble y3; / 记录每次迭代产生的新的近似解double maxerr; / maxerr 记录相邻两个数值解想x 与 y 的差的绝对值的最大的那个分量。 int k,i,j; printf(此线性方程组的增广矩阵如下： n );/ 下面循环
2、是输出增广矩阵（A, b）for (k=0;kn;k+)for (i=0;in;i+) printf(%f ,Aki); printf(%f n, bk);/ 输出增广矩阵（A, b）完 printf(此方程组的精确解为 x=(-4, 3, 2),n); / 显示精确解以便于对比/ 输出初始迭代向量printf(n Jacobi迭代的初始值x0=();for (i=0;in;i+) printf( %f ,xi); printf() n);/ 输出初始迭代向量完printf(误差限是%f n,tol); / 输出误差限printf( Jacobi迭代解序列X(k) max|x(k+1)-x(k
3、)| n );printf(x%d = ,k=0);for(i=0;in;i+)printf(%f ,xi); printf(n);/ Jacobi迭代 k=0; do for(i=0;in;i+) double T=0.0; for(j=0;jn;j+) if(j=i)continue; T=T+Aij*xj; yi=(bi-T)/Aii;/ Jacobi迭代完/求相邻两个数值解想x 与 y 的差的绝对值的最大的那个分量maxerr=fabs(y0-x0);for (j=1;jn;j+) if (maxerrfabs(yj-xj)maxerr=fabs(yj-xj); /求相邻两个数值解想x 与 y 的差的绝对值的最大的那个分量完for(i=0;in;i+)xi=yi; k=k+1; printf(x%d=,k);for(i=0;itol);运行结果：
