大型矩阵迭代.docx

苏州大学本科生毕业设计（论文）18目录摘要31前言52迭代法的思想62.1迭代法的基本概念62.2迭代法的收敛条件63常见迭代法73.1Jacobi迭代法73.1.1Jacobi迭代法的计算公式73.1.2Jacobi迭代法的收敛条件83.2Gauss迭代法93.2.1Gauss迭代法的计算公式93.2.2Gauss迭代法的收敛条件103.3SOR迭代法103.3.1SOR迭代法的计算公式103.3.2SOR迭代法的收敛条件113.4CG迭代法113.5PCG迭代法123.6GMRES迭代法134数值算例分析13算例4.113算例4.2145总结16参考文献17致谢18摘要这篇论文介绍了大型矩阵运算时要用到的迭代法和其思想，并简单介绍了常用的几种迭代方法：如基于矩阵分解原理的Jacobi迭代法，Gauss迭代法和SOR迭代法；属于共轭方向法的CG法和PCG法；最后介绍了适用于解决非线性方程组的GMRES迭代法及他们其中一些迭代法的收敛条件.并通过对不同矩阵的具体运算比较了他们的迭代次数，简单分析了几种迭代方法的优劣.关键字:迭代法思想，Jacobi迭代法，Gauss迭代法，SOR迭代法,CG迭代法,PCG迭代法，GMRES迭代法，迭代次数，收敛体条件.AbstractThisessayintroducestheiterativemethodwhichisneededduringlargematrixoperation.Meanwhile,someofcommonmethodsofiterationarealsosuggestedinthisessay,forinstance:Jacobianinteration,GaussinterationandsuccessiveoverrelaxationmethodwhichwereallbasedontheMatrixdecompositiontechnique;conjugategradientalgorithmandPreprocessingconjugategradientalgorithmwhichbothbelongtoconjugatedirectionmethod;GMRESalgorithmwhichisusedtosolvenon-linearequationandsomeoftheirownconditionofconvergence.Inaddition,bycomparingtheiriterationsteps,Iroughlyanalysedtheirmeritsanddemerits.Keywords:iterativemethod，Jacobianinteration,Gaussinteration,successiveoverrelaxationmethod,conjugategradientalgorithm,Preprocessingconjugategradientalgorithm,GMRESalgorithm,iterationstepsandconditionofconvergence.1前言矩阵自19世纪被英国数学家凯利提出来，历经了200多年的发展，在现代生活中有及其广泛的应用，下面我将列举几项矩阵在现实生活中的应用：矩阵可用于预测水质量状况以及用于处理水污染等问题;矩阵还可以用于模拟飞机飞行所需要的飞行环境;在物理学领域的磁电场的探查测量方面也有很广泛的应用;许多图片在电脑上的存储方式就是以矩阵的格式存储的,以便对图片进行调整和美化的时候可以利用数学公式对图片进行处理等。除此以外，矩阵还在许多领域有着极为广泛的应用，这里就不一一举例了。而这些矩阵通常都是大型矩阵，如何高效的处理大型矩阵和与其相关的线性方程组成了一个值得深思的问题。而迭代法在处理零元素较多的大型稀疏矩阵时具有运算和存储两方面的优势。比较典型的有基于矩阵分解原理的Jacobi迭代法，这是由伟大的普鲁士数学家雅可比在19世纪提出来的，在当时Jacobi迭代法具有许多优点，其计算公式简单，每一次的迭代中只需计算一次矩阵和向量的乘法。这也是最早期的迭代法。同雅可比迭代法一样，高斯-塞德尔和超松弛迭代法也是基于矩阵分解原理，这些都是传统的迭代法，出现时间较久。此外，还有其他一些优异的矩阵迭代算法，如以极小化的方法来讨论方程组的解的著名的共轭梯度（CG）算法，CG法出现在20世纪中期，它有许多优异的性质：比如相较于传统迭代法，CG法收敛速度较快；此外CG法只有极少量的非向量运算。然而由于实际运算中舍入误差的积累以及正交性的逐渐丧失，CG法长时间并未得到广泛应用。直到20世纪70年代，预处理共轭梯度算法（PCG算法）出现了，它通过对CG算法中未加工过的矩阵进行某种形式的加工改造，使