jacobi迭代法的.doc

实 验 报 告 实验课程名称 计算机数值方法 实验项目名称 年 级 专 业 学生姓名 学 号 理 学 院 实验时间： 2012 年 月 日 学生所在学院： 专业： 班级 姓 名 学 号 实验组 1 实验时间 指导教师 成 绩 实验项目名称 迭代法 实验目的及要求： 实验（或算法）原理： 迭代法的基本思想是对元线性方程组，将其变形为等价方程组 ，其中成为迭代矩阵。从某一取定的初始向量出发,按照一个适当的迭代公式 ,逐次计算出向量 （）,使得向量序列收敛于方程组的精确解.设元线性方程组为： 设矩阵非奇异，且，上式变形为 ，该式称为迭代法。写成矩阵形式： 设， ， 则 。 令 则有 其迭代格式为 这就是迭代法的矩阵形式。 根据计算公式编程即可求解线性方程组。 实验硬件及软件平台： PC机，vc++6.0，Internet网 实验步骤： 1．根据算法事先写出相应程序。 2．启动PC机，进入vc集成环境，输入代码。 3．编译调试。 4．调试通过，计算出正确结果。 实验内容（包括实验具体内容、算法分析、源代码等等）： 算法设计 (1). (2) () (3),若是，输出，若否，置 。 编写相应的程序上机调试。 //jacobi迭代法程序： #includemath.h #includestdio.h #define N 3//方程阶数 void jacobi(double a[N][N],double b[N],double x[N],double e) {int i,j; double p,x1[N],ep,sub;//最大范数 do {for(i=0;iN;i++) x1[i]=x[i]; for(i=0;iN;i++) {p=0; for(j=0;jN;j++) if(j!=i) p=p+a[i][j]*x1[j]; x[i]=(b[i]-p)/a[i][i]; } ep=0; for(i=0;iN;i++) {sub=fabs(x[i]-x1[i]); if(epsub) ep=sub; } } while(epe); } void main() {double a[N][N]={{4,1,1},{1,4,1},{1,1,4}}; double b[N]={7,10,7}; double x[N]={0,0,0}; double e; int i; e=1e-4; jacobi(a,b,x,e); for(i=0;iN;i++) printf(x[%d]=%f\n,i,x[i]); } 3、用上述程序验证线性方程组=的解为 3运行结果： 4 、用上述程序解线性方程组。 4题运行结果： 实验结果与讨论： 解线性方程组有直接法和迭代法，直接法计算量较小，但需较大的储存量，且程序复杂，，而迭代法恰与直接法相反，对于迭代法我们需要讨论其收敛性及误差估计： 若所求的迭代序列{}收敛于解线性方程组的解，即，则迭代法收敛，否则称为发散。迭代法收敛的充要条件是：jacobi迭代法收敛??;若‖B‖1?则jacobi代法收敛;称为迭代法的第步误差向量，若迭代法收敛，当充分大时，如时，取作为解的近似向量，即线性方程组的近似解。 指导教师意见： 签名： 年 月 日 1