计算声学第二章线性方程组的数值解法.pptVIP

第二章 线性方程组的数值解法 主要内容： 本章着重讨论了利用计算机求解系数行列式不为零的n 阶非齐次线性方程组 的两类主要解法——直接法和 迭代法。 要求： 熟练掌握高斯列主元消去法，雅可比（Jacobi）迭代 法，高斯-塞德尔（Gauss-Seidel）迭代法，了解逐次超松 弛（SOR）迭代法，分析迭代法收敛性。 第二章 线性方程组的数值解法 直接法： 又叫精确法，是指在没有舍入误差的假设下，经过有限 步运算即可求得方程组的精确解的方法。 目前较实用的直接法主要有高斯（Gauss）消去法及其 变形——选主元消去法及矩阵的三角分解法等。 适用于求解低阶稠密矩阵方程组及某些大型稀疏方程组 第二章 线性方程组的数值解法 迭代法： 用某种极限过程去逐步逼近线性方程组精确解的方法， 即从一个初始向量出发，按照一定的迭代公式产生一个向量 序列，使其收敛到方程组的解。 特点： 所需存储单元少，程序设计简单，但是存在收敛性与收 敛速度问题。 适用于求解大型稀疏矩阵方程组。 §1 线性方程组的直接解法 高斯消去法基本思想： 通过逐步消元，把线性方程组的系数矩阵约化为上三角 矩阵，得到一个与原方程组等价的三角方程组，再通过回代 过程逐一求出各未知数。消元过程是通过对增广矩阵的一系 列初等行变换得到的。 §1 线性方程组的直接解法 设三元线性方程组如下： 若 ，则将方程（2.1.1）分别乘以 和 ，然后 加到方程（2.1.2）、（2.1.3）中，消去这两个方程中的 得到同解方程组 §1 线性方程组的直接解法 同理，若 ，消去方程（2.1.5）中的 ，方程组又 可化为 通过逐步消元将方程组化为上三角方程组求解的方法称 为高斯消去法。 原方程组逐步化为同解的上三角方程组的过程称为消元 过程，按照方程相反顺序逐步求解上三角方程组的过程称为 回代过程。 §1 线性方程组的直接解法 一般形式： 设n元线性方程组 ，其中 如果约化的主元 ，则可通过高斯消去法 将方程组 约化为同解的三角方程组 ，即 §1 线性方程组的直接解法 并由回代公式 求得方程组的解。主对角线上的元素 称为主元。 §1 线性方程组的直接解法 特点： 每次都按照系数矩阵的主对角线上元素的顺序依次消去 对角线下方的元素。 存在的问题： （1）如果某个主元 ，则消元过程无法继续进行； （2）即使主元不为零，但是与该主元所在列的主对角线以 下的各元素相比，如果该主元的绝对值很小，用其做除数时 很小的舍入误差会引起计算结果的失真和误差的扩散。 §1 线性方程组的直接解法 例：如方程组 ，其准确到小数点后第9位 的解为 ， 用第一个方程消去第二个方程中的 得(计算中用4位浮点数) 解得 ，结果严重失真。 §1 线性方程组的直接解法 交换两个方程位置，改写为 对此方程进行消元得 解得 因此，在消元过程中，应避免选取绝对值较小的数做主 元，否则可能导致结果错误，计算失败。 §1 线性方程组的直接解法 列主元消去法： 列主元消去法能避免高斯消去法中出现的两个问题，是直接 法中最常用的一种方法。在消元过程中，依次按列选取绝对 值最大的元素作为主元，只进行行变换。 特点： 在系数矩阵中依次按列在主对角线及以下的元素中，选取绝 对值最大的元素作为主元，将它换至主对角线上，然后用它 消去主对角线以下的元素，最后方程组变为同解的三角方程 组进行求解。 §1 线性方程组的直接解法 例：求解线性方程组 解：由方程组的增广矩阵 最后利用回代可求得方程组的解为 。 §2 线性方程组的迭代解法 迭代法是一类逐次逼近方法，它按照一个适当的计算法则， 从某一初始向量 出发，逐次计算得到一个向量序列 ，使得此向量序列收敛于方程组的解 。 迭代法要解决的主要问题： （1）构造一种迭代格式，把所给方程组