毕业论文 线性方程组的迭代法及程序实现

学校代码：11517学号：200810111217HENANINSTITUTEOFENGINEERING毕业论文题目线性方程组的迭代法及程序实现学生姓名专业班级学号系（部）数理科学系指导教师(职称)完成时间2012年5月20日河南工程学院论文版权使用授权书本人完全了解河南工程学院关于收集、保存、使用学位论文的规定，同意如下各项内容：按照学校要求提交论文的印刷本和电子版本；学校有权保存论文的印刷本和电子版，并采用影印、缩印、扫描、数字化或其它手段保存论文；学校有权提供目录检索以及提供本论文全文或者部分的阅览服务；学校有权按有关规定向国家有关部门或者机构送交论文的复印件和电子版；在不以赢利为目的的前提下，学校可以适当复制论文的部分或全部内容用于学术活动。论文作者签名：年月日河南工程学院毕业设计(论文)原创性声明本人郑重声明：所呈交的论文，是本人在指导教师指导下，进行研究工作所取得的成果。除文中已经注明引用的内容外，本论文的研究成果不包含任何他人创作的、已公开发表或者没有公开发表的作品的内容。对本论文所涉及的研究工作做出贡献的其他个人和集体，均已在文中以明确方式标明。本学位论文原创性声明的法律责任由本人承担。论文作者签名：年月日河南工程学院毕业设计（论文）任务书题目：线性方程组的迭代法及程序实现专业：信息与计算科学学号：姓名一、主要内容：通过本课题的研究，学会如何运用有限元方法来解决线性代数方程组问题，特别是Gaussie-Seidel迭代法和Jacobi迭代法来求解线性方程组。进一步学会迭代方法的数学思想，并对程序代码进行解析与改进，这对于我们以后学习和研究实际问题具有重要的意义。本课题运用所学的数学专业知识来研究，有助于我们进一步掌握大学数学方面的知识，特别是迭代方法。通过这个课题的研究，我进一步掌握了迭代方法的思想，以及程序的解析与改进，对于今后类似实际问题的解决具有重要的意义。二、基本要求：1．学会编写规范论文，独立自主完成。2．运用所学知识发现问题并分析、解决。3.通过对相关资料的收集、整理，最终形成一篇具有自己观点的学术论文，以期能对线性方程组迭代法的研究发展有一定的实践指导意义。4.在毕业论文工作中强化英语、计算机应用能力。完成期限：2012年月指导教师签名：专业负责人签名：年月日目录中文摘要…………………………………………………………………………Ⅰ英文摘要…………………………………………………………………………Ⅱ1综述............................................12经典迭代法概述..................................32.1Jacobi迭代法...............................32.2Gauss—Seidel迭代法........................42.3SOR（successiveoverrelaxation）迭代法....42.4SSOR迭代法.................................52.5收敛性分析.................................52.6数值试验...................................63matlab实现的两个例题............................83.1例1迭代法的收敛速度.......................83.2例2SOR迭代法松弛因子的选取..............12致谢...............................................16参考文献...........................................17附录............................................19线性方程组的迭代法及程序实现Ⅰ线性方程组的迭代法及程序实现摘要线性代数方程组的迭代方法是一种极限方法是解大型稀疏矩阵方程组的有效方法。它的基本思想是用某种极限过程去逐步逼近线性方程组的精确解,是一种逐步逼近的方法。迭代法将n阶线性方程组变形为某种迭代公式。对于任意给定的迭代初始值，由某一迭代格式便可生成一向量序列，我们的目的是求解方程组的解，因此我们会希望向量序列的极限逼近方程组的解。本文首先介绍了求解大型线性方程组的主要迭代算法，对一些经典迭代法(Jacobi方法、Gauss—Seidel方法、SOR方法和SSOR方法)进行了详细的讨论，其次着重讨论了经典迭代法的收敛性，详细总结并给出了各种迭代方法的收敛性定理，并通过举例及其Matlab程序实现进一步阐述了迭代法的收敛性。关键字线性方程组/Jacobi迭代法/Gauss-Seidel方法/SOR方法/收敛性线性方程组的迭代法及程序实现ⅡIterativemethodandproceduresforimplementationofthelinearequationsABSTRACTTheiterativemethodoflinearalgebraicequationsisanextrememethodisaneffectivemethodforthesolutionoflargesparsematrixequations.Thebasicideaistoacertainlimitprocesstograduallyapproachtheexactsolutionoflinearequations,astep-by-stepapproximationmethod.Theiterativemethodwillbeiterativeformulaforadeformationofnlinearequations.Foranygiveniterationoftheinitialvalue,byaniterativeschemecangenerateavectorsequence,ouraimisthesolutiontosolvingtheequations,sowewillwanttolimitapproximationthesolutionofequationsofvectorsequences.Thispaperfirstintroducesthemainiterativealgorithmforsolvinglargelinearequations,adetaileddiscussionofsomeclassicaliterativemethod(Jacobimethod,Gauss-Seidelmethod,SORandSSORmethods),followedfocusedontheclassicaliterativeconvergenceofthemethod,summarizedindetailandgivesavarietyofiterativemethodsconvergencetheorem,andfurtherelaboratedbyexampleandMatlabprogramiterativemethods.KEYWORDSlinearequations，Jacobiiterativemethod，Gauss-Seidelmethod，theSORmethod，convergence线性方程组的迭代法及程序实现11综述在科学研究和大型工程设计中出现了越来越多的数学问题，而这些问题往往需要求数值解。在进行数值求解时，经离散后，常常归结为求解行如Ax=b的大型线性方程组。20世纪50年代至70年代，由于电子计算机的发展，人们开始考虑和研究在计算机上用迭代法求线性方程组Ax=b的近似解，用某种极限过程去逐渐逼近精确解，并发现了许多非常有效的迭代方法。迭代法是按照某种规则构造一个向量序列｛x)(k｝，使其极限向量x是Ax=b的精确解。因此，对迭代法来说一般有下面几个问题：（1）如何构造迭代序列？（2）构造的迭代序列是否收敛？在什么情况下收敛？（3）如果收敛，收敛的速度如何？我们应该给予量的刻划，用以比较各种迭代法收敛的快慢。（4）因为计算总是有限次的，所以总要讨论近似解的误差估计和迭代过程的中端处理问题，这又和舍入误差的分析有关。迭代法具有需要计算机存储单元少、程序设计简单、原始系数矩阵在计算过程中始终不变等优点。例如Jacobi方法、Gauss-Seidel方法、SOR方法、SSOR方法，这几张迭代方法是最常用的一阶线性定常迭代法。大量偏微分方程的离散形式是大规模线性代数方程组，其数值计算是科学工程计算的核心，占有绝大部分的总体运算时间，解大规模稀疏线性方程组的Krylov子空间方法显示出与众不同的有效性。当矩阵是对称正定时，常用的方法是具有短递推的共轭度方法(CG)。系数矩阵不对称时，常用的方法中有完全正交化方法（FOM）和广义最小参量方法（GMRES）。还有很多迭代方法正在被人们发现和研究，新的有效的方法层出不穷，其中基于大型稀疏非Bermitian的正定阵的系数矩阵的Bermitian和Skew-Bermitian分裂的BSS方法，IBSS方法等具有非常好的实用性。但没有一种算法是通用的，对于具体问题必须根据所得到的线性方程组和算法的特点进行选择。MATLAB是Mathworks公司的产品，MATLAB的产生是与数学计算紧密联系在一起的，其基本数据结构是矩阵，它的表达式与数学工程计算中使用的形式十分相似，便于用户学习和使用．系统包括5个部分：MATLAB语言、MATLAB工作环境、MATLAB图形处理系统、MATLAB数学函数库和MATLAB应用程序接线性方程组的迭代法及程序实现2口．其主要功能包括：数值计算；符号计算；数据分析和可视化；文字处理；Simulink动态仿真．在数值计算中，线性方程组的求解是一个很重要的问题．用MATLAB来求解线性方程组，有几种方法，非常简单，通过对一些矩阵和函数的操作可以轻松地得到线性方程组的解，不需要使用者掌握任何程序设计语言，对迭代法只需编写简单的程序线性方程组的迭代法及程序实现32经典迭代法概述20世纪50年代至70年代，人们开始考虑和研究用迭代法求解线性方程组bAx（2-1）的近似解，发展了许多有效的方法，其中有Jacobi方法、Gauss-Seidel方法、SOR方法、SSOR方法，这几种迭代方法均属一阶线性定常迭代法，即若系数矩阵A的一个分裂：NMA；M为可逆矩阵，线性方程组（2-1）化为：bXNM)(；;bNXMX.11bMNXMX得到迭代方法的一般公式：gBXXkk)()1(（2-2）其中：NMB1,bMg1。2.1Jacobi迭代法若D为A的对角素构成的对角矩阵，且对角线元素全不为零。系数矩阵A的一个分解：A=D－(E+F);这里D为A的对角矩阵，E为严格下三角阵，F为严格上三角阵。其中：nnaaaD002211000001,21323121nnnnaaaaaaE0000,122311312nnnnaaaaaaFJacobi迭代的矩阵形式为：线性方程组的迭代法及程序实现4JkJkgXBX)()1(（2-3）(2-3)式中：bDgJ1；FEDBJ1，称JB为Jacobi迭代矩阵。其计算公式为：;1,,,,110010nijjkjijiiikiTnxabaxxxx（2-4）2.2Gauss—Seidel迭代法对于非奇异方程组，若D为A的对角素构成的对角矩阵，且对角线元素全不为零；系数矩阵A的一个分解：FEDA)((2-5)Gauss—Seidel迭代矩阵形式为bEDFxEDxkk111