解线性方程组地迭代法

《解线性方程组地迭代法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案解线性方程组的迭代法Haha送给需要的学弟学妹摘要:因为理论的分析表明，求解病态的线性方程组是困难的，但是实际情况是否如此，需要我们来具体检验。系数矩阵H为Hilbert矩阵，是著名的病态问题。因而决定求解此线性方程组来验证上述问题。详细过程是通过用Gauss消去法、J迭代法、GS迭代法和SOR迭代法四种方法求解线性方程组。关键词：病态方程组、Gauss消去法、J迭代法、GS迭代法、SOR迭代法目录：一、问题背景介绍二、建立正确额数学模型三、求解模型的数学原理1、Gauss消去法求解原理2、Jacobi迭代法求解原理3、G-S迭代法求解原理4、SOR迭代法求解原理5、Jacobi和

2、G-S两种迭代法收敛的充要条件四、计算过程（一）Hilbert矩阵维数n=6时1、Gauss消去法求解2、Jacobi迭代法求解3、G-S迭代法求解4、SOR迭代法求解（二）Hilbert矩阵维数n=20、50和100时1、G-S迭代法求解图形2、SOR迭代法求解图形五、编写计算程序六、解释计算结果1、Gauss消去法误差分析2、G-S迭代法误差分析3、SOR迭代法误差分析G-S迭代法与SOR迭代法的误差比较七、心得体会正文：一、问题背景介绍。理论的分析表明，求解病态的线性方程组是困难的。实际情况是否如此，会出现怎样的现象呢？二、建立正确的数学模型。考虑方程组的求解，其中系数矩阵H为Hilbe

3、rt矩阵，这是一个著名的病态问题。通过首先给定解（为方便计算，笔者取x的各个分量等于1），再计算出右端这样的解就明确了，再用Gauss消去法、J迭代法、GS迭代法和SOR迭代法四种方法分别求解将求解结果与给定解比较，而后求出上述四种方法的误差，得出哪种方法比较好。三、求解模型的数学原理。1、Gauss消去法求解原理对于（A非奇异）求解时，可以先将A分解成一个下三角矩阵L和一个上三角矩阵U的乘积，即，就可以通过(1.1)文档大全实用标准文案求解出的值。接下来就具体讲讲如何将A分解成L和U，也就是Gauss消去法。欲把一个给定的矩阵A分解为一个下三角阵L与一个上三角阵U的乘积，最自然的做法便是通过

4、一系列的初等变换，逐步将A约化为一个上三角阵，而又能保证这些变换的乘积是一个下三角阵。这可归结为：对于一个任意给定的向量找一个尽可能简单的下三角阵，使经这一矩阵作用之后的第至第个分量均为零。能够完成这一任务的最简单的下三角阵便是如下形式的初等下三角阵：其中即这种类型的初等下三角阵称作Gauss变换，而称向量为Gauss向量。对于一个给定的向量我们有由此立即可知，只要取便有当然，这里我们要求而后经过多次变换可以得到(1.2)从而求出上三角阵U，而后通过求得下三角阵(1.3)将(1.2)和(1.3)带入到(1.1)式中求出的值即可。2、J迭代法求解原理考虑非奇异线性代数方程组令(1.4)其中文档大

5、全实用标准文案那么(1.4)式和合并后可以写成其中(1.5)若给定初始向量并代入(1.4)式右边，又可得到一个向量；一次类推，有(1.6)这就是所谓的Jacobi迭代法，其中叫做Jacobi迭代法的迭代矩阵，叫做Jacobi迭代法的常数项。3、GS迭代法求解原理注意到Jacobi迭代法中各分量的计算顺序是没有关系的，先算那个分量都一样。现在，假设不按Jacobi迭代格式，而是在计算的第一个分量用的各个分量计算，但当计算的第二个分量时，因已经算出，用它代替，其他分量仍用。类似地，计算时，因都已算出，用它们代替其他分量仍用的分量，于是有(1.7)我们称这种迭代格式为Gauss-Seidel迭代法，

6、简称为G-S迭代法。它的一个明显的好处是在编写程序是存储量减少了。如果存在，G-S迭代法可以改写成(1.8)我们把叫做G-S迭代法的迭代矩阵，而把叫做G-S迭代法的常数项。4、SOR迭代法求解原理我们知道，G-S迭代法的迭代格式为现在令则有(1.9)这就是说，对G-S迭代法来说，可以看作在向量上加上修正项而得到的。若修正项的前面加上一个参数，便得到松弛迭代法的迭代格式(1.10)用分量形式表示即为文档大全实用标准文案(1.11)其中叫做松弛因子。当时，相应的迭代法叫做超松弛迭代法；当时，叫做低松弛迭代法；当时，就是G-S迭代法。我们把超松弛迭代法简称为SOR迭代法。因为存在，所以(1.10)式

7、可以改写为其中叫做松弛迭代法矩阵。而SOR迭代法收敛的充要条件是(1.12)由(1.12)式知，SOR迭代法的谱半径依赖于，当然会问：能否适当选取使收敛速度最快？这就是选择最佳松弛因子的问题。经过相关计算可知，随着从0增加，减少，直至(1.13)时，达到极小(1.14)再增加时，开始增加。因此，称为最佳松弛因子。5、Jacobi和G-S两种迭代法收敛的充要条件Jacobi迭代法和G-S迭代法两种迭