线性代数方程组.pptVIP

Jacobi迭代法的计算步骤：浙江大学研究生学位课程《实用数值计算方法》*Gauss--Seidel迭代法。01-S迭代法的基本思想来自Jacobi02迭代法。只是迭代矩阵B不同。03也可写为：若用（3-29）式记号，则G--S迭代式可表示为：或其中这就是Gauss--Seidel迭代格式，它0201030506041有关收敛性和收敛速度的讨论完全与Jacobi方法相同，只是迭代阵为L。23.3.3逐次超松弛（SOR）法SOR法定义为：三种迭代法的比较：01一般情况下，J方法与G--S方法比较02并无优劣。收敛情况与速度均不一定。03例：042.但是具有相容次序的矩阵，在相同精度要求下，迭代次数分别为：Jacobi方法：1154G--S方法：578SOR方法：59可见对这类矩阵，G--S法比J方法快一倍，而SOR法的收敛速度可提高一个数量级。最后介绍一类对于三种方法都收敛的矩阵。先定义可约矩阵和对角占优矩阵：排列阵P，使得其中A11为r阶方阵，A22为（n-r）阶方阵不是可约矩阵就是不可约矩阵(1).称n阶方阵为可约矩阵，如果存在n阶若上式严格不等号成立，则称A为严格对角占优矩阵。(2).称n阶方阵A是对角占优矩阵，如果(3).称A是不可约对角占优矩阵。如果A是不可约的且对角占优，而（3-42）式中至少对一个i有严格不等号成立。我们有如下结论：1.若A是严格对角占优或不可约对角占优矩阵，则A非奇异。2.若线性方程组系数A是上述矩阵，则：1）J方程和G--S方法都收敛。2）当SOR方法收敛。3.4稀疏矩阵前面介绍了直接法和迭代法。但是实践中如何选择计算方法是一个很复杂的问题。需要考虑的因素很多，主要有：1.方程组的性质；（病态、对角占优、正定等）2.相类似问题出现的次数；3.方程组的阶数；4.计算机的容量、字长、运算速度等；5.系数矩阵的结构；（稀疏、三对角等）6.对结果的精度要求。3.4浙江大学研究生学位课程《实用数值计算方法》*迭代法是解超高阶线性代数方程组01的重要方法之一，但是它也有缺点：02不适合解病态问题；03精度不高；04收敛速度依赖于加速因子的选取；05收敛性不能保证；06所需乘除运算量大。07在这些方面，相对来说，直接法运算08步骤一定，运算次数有限，精度高。对解09大型稀疏矩阵问题可以利用和保持系数矩10阵的稀疏性来降低存储容量和缩短计算时11间。1201但是，为了保持系数矩阵的稀疏性。必须02使用特殊的计算方法和高级程序设计技术。03近年来在这方面有较大的发展。其优越性04大大超过了迭代法。05在这一节里，不详细介绍计算方法，06只给出一些常用的标准稀疏矩阵形式。同07时给出一种一般稀疏矩阵的存储结构。这08一节里还给出一个有效的计算公式。3.4.1稀疏矩阵的结构和存储从本质上讲，解稀疏问题的方法与一般Gausss消去法LU分解法并没有什么不同，只不过前者更注重于运算过程中对零元素填入量的控制。稀疏系数矩阵问题的直接法必定依赖于矩阵的稀疏形状，下面列举的是常见的标准的稀疏矩阵。1.稀疏矩阵的结构关于L,U元素的存放从（3-22），（3-23）可以看出。A的元素aij在计算出或以后就不再有用了。故L,U的非零元素便可以存放在矩阵A中的相应位置上了。最后，A所存放的元素是：