数值分析-雅可比迭代课程设计.docVIP

郑州轻工业学院 《数 值 分 析》 课 程 设 计 报 告 题 目： 雅可比迭代 姓 名： 院（系）： 数学与信息科学系 专业班级： 信息与计算科学10-01班 学 号： 指导教师： 时 间： 2013年1月07日至2013年1月14日 摘要 随着工程技术中产生的大型稀疏矩阵方程组目录 目录 《数 值 分 析》 1 课 程 设 计 报 告 1 题 目： 1 雅可比迭代法 1 摘要 3 目录 4 1 理论基础 5 2 算法分析 5 2.1 功能分析 6 2.2 算法分析 6 3 程序设计 6 3.1 具体程序 6 3.2 具体例子的引入 11 3.3 模块设计 12 3.3.1为所用的二维和一维数组申请空间： 12 3.3.2迭代矩阵的求解： 13 3.3.3迭代计算过程设计： 14 3.3.4迭代计算过程设计： 14 4 总结 15 5 参考文献 16 1 理论基础 这次报告的理论基础主要有两个方面：一是基于对C语言编程知识的了解，在计算中能够把矩阵转换成数组的形式，利用循环在C的环境下进行迭代；二是对雅可比迭代法的掌握，在数值分析中，雅可比迭代法是求解大型方程组的比较适合的方法，迭代公式非常简单，算法也是简单的矩阵和向量的乘法。 2 算法分析 算法方面主要是两个方面的结合：首先是雅可比的理论迭代分析，线性方程组的一般形式为： 以矩阵的形式来分析Jacobi迭代法： ， =， 其中，。这就是Jacobi迭代法的矩阵表示形式。 然后是在C语言中先把上述分析的矩阵都装入数组，用一个for循环计算出迭代矩阵，进而再用循环迭代求解。 2.1 功能分析 题目要求用雅可比迭代法解大型方程组的程序，并用数值例子计算。所以在功能上，这次的报告主要实现的是用雅可比迭代法解大型线性方程组，并且用C语言实现迭代的求解过程。 2.2 算法分析 程序上主要的设计在于两个模块：一是计算雅可比迭代的迭代方程；二是根据结果精度的需要进行迭代求出结果。 3 程序设计 程序的设计部分包括了整体程序，具体例子的引入，具体分模块分析。 3.1 具体程序 #includestdio.h #includemath.h #includemalloc.h void main() { int i,j,k,N; printf(请输入所要求解方程组系数矩阵的阶数N:\n); scanf(%d,N); double t,r,h=0; double **a,**b,**c,**p,*e,*f,*d; a=(double **)malloc(sizeof(double)*N); b=(double **)malloc(sizeof(double)*N); c=(double **)malloc(sizeof(double)*N); p=(double **)malloc(sizeof(double)*N); d=(double *)malloc(N*sizeof(double)); e=(double *)malloc(N*sizeof(double)); f=(double *)malloc(N*sizeof(double)); for(i=0;iN;i++) { a[i]=(double *)malloc(sizeof(double)*N); b[i]=(double *)malloc(sizeof(double)*N); c[i]=(double *)malloc(sizeof(double)*N); p[i]=(double *)malloc(sizeof(double)*N); } printf(从第一个方程开始依次输入未知数的系数,每输入一个系数按回车键输入下一个:\n); for(i=0;iN;i++) { for(j=0;jN;j++) scanf(%lf,p[i][j]); /*系数矩阵*/ } printf(依次输入每个方程右边的常数并按回车键结束:\n); for(i=0;iN;i++) scanf(%lf,d[i]); printf(输入结果需要的精度:\n); scanf(%lf,r); printf(常数项矩阵b:\n); for(i=0;iN;i++) /*显示等号右端常数*/ printf(%lf\n,d[i]); for(i=0;iN;i++) { for(j=0;jN;j++) { if(i==j) /*为了