解线性方程组地迭代法.docVIP

解线性方程组地迭代法 解线性方程组地迭代法 解线性方程组地迭代法 实用标准文案 解线性方程组的迭代法 Haha 送给需要的学弟学妹 摘要 :因为理论的分析表明，求解病态的线性方程组是困难的，但是实际情况是否如此，需要我们来具体检 验。系数矩阵H 为 Hilbert矩阵，是著名的病态问题。因而决定求解Hx b此线性方程组来验证上述问 题。 详细过程是通过用Gauss 消去法、 J 迭代法、 GS 迭代法和 SOR 迭代法四种方法求解Hxb 线性方程组。 关键词 ：病态方程组、 Gauss 消去法、 J 迭代法、 GS 迭代法、 SOR 迭代法 目录： 一、问题背景介绍 二、建立正确额数学模型 三、求解模型的数学原理 1 、 Gauss 消去法求解原理 2 、 Jacobi 迭代法求解原理 3 、 G-S 迭代法求解原理 4 、 SOR 迭代法求解原理 5 、 Jacobi 和 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 迭代法的误差比较 七、心得体会 正文 ： 一、问题背景介绍。 理论的分析表明，求解病态的线性方程组是困难的。实际情况是否如此，会出现怎样的现象呢？ 文档大全 实用标准文案 二、建立正确的数学模型。 考虑方程组 Hx b 的求解，其中系数矩阵 H 为 Hilbert 矩阵， H ( hi , j ) n n , 1 , i , j 1,2,L , n hi, j i j 1 这是一个著名的病态问题。通过首先给定解（为方便计算，笔者取 x 的各个分量等于 1 ），再计算出右端 b Hx, 这样Hx b的解就明确了，再用Gauss消去法、J迭代法、GS迭代法和SOR迭代法四种方法分别求解 Hx b, 将求解结果与给定解比较，而后求出上述四种方法的误差，得出哪种方法比较好。 三、求解模型的数学原理。 1 、 Gauss 消去法求解原理 对于 Ax b（ A 非奇异）求解时， 可以先将 A 分解成一个下三角矩阵 L 和一个上三角矩阵 U 的乘积， 即 A LU ，就可以通过 Ly b, Ux y, (1.1) 求解出 x 的值。 接下来就具体讲讲如何将 A 分解成 L 和 U ，也就是 Gauss 消去法。 欲把一个给定的矩阵 A 分解为一个下三角阵 L 与一个上三角阵 U 的乘积，最自然的做法便是通过一 系列的初等变换，逐步将 A 约化为一个上三角阵，而又能保证这些变换的乘积是一个下三角阵。这可归结 为：对于一个任意给定的向量 x Rn ,找一个尽可能简单的下三角阵， 使 x 经这一矩阵作用之后的第 k 1 至 第 n 个分量均为零。能够完成这一任务的最简单的下三角阵便是如下形式的初等下三角阵： Lk I lk ekT , 其中 lk 0,L ,0, lk T , 1,k ,L ,lnk 即 1 O Lk 1 lk 1,k 1 M O ln ,k 1 这种类型的初等下三角阵称作 Gauss 变换，而称向量 lk 为 Gauss 向量。 x x1,L T Rn , 我们有 对于一个给定的向量 , xn Lk x x1,L , xk , xk xklk 1,k ,L , xn xk lnk T 1 . 由此立即可知，只要取 lik xi , i k 1,L , n, xk 便有 Lk x x1 ,L , xk ,0, L T ,0 . 当然，这里我们要求 xk 0. 而后经过多次变换可以得到 文档大全 实用标准文案 Ln Ln 1 L L1A U, (1.2) 从而求出上三角阵 U ，而后通过 LA U , L LnLn 1L L1 求得下三角阵 L UA 1. (1.3) 将 (1.2) 和(1.3) 带入到 (1.1) 式中求出x的值即可。 2 、 J 迭代法求解原理 考虑非奇异线性代数方程组 Ax b. 令 A D L U , (1.4) 其中 A aij , D diag a11, a22 ,L , ann , 0 a21 0 L a31 a32 0 , M M O O an1 an 2 L an ,n 10 0 a12 a13 L a1n 0 a23 L a2n U O