代数方程组的求解方法.PDFVIP

第八章 代数方程组的数值解法 王 娴 西安交通大学航天学院 内容 本章内容：对流动传热问题的控制方程进行离散后所 形成的代数方程组的求解方法，这是对物理过程进行 数值模拟的最后一个重要环节。  代数方程组求解方法概述  求解代数方程组的迭代法  迭代法的收敛性及加速收敛的方法 1. 代数方程组求解方法概述 1.1 多维流动换热问题离散方程组系数矩阵的特点 压力泊松方程，导热问题，离散方程总可以表示成以下形式： AP = b x方向L个节点，y方向M个节点，共有L x M 个节点 系数矩阵A中各主对角元素均大于零； 除了主对角线及其上下相邻几个位置上元素不为零外，只有 离开主元素 L 个元素的位置上才为非零元素，其余全为零； 1. 代数方程组求解方法概述 AP = b 1. 1 多维流动换热问题离散方程组系数矩阵的特点 0 L+1个元素 0 0 0 1. 代数方程组求解方法概述 1.2 求解代数方程组的直接解法与迭代法  直接解法：通过有限步的数值计算可以获得代数方 程真解的方法（设不考虑舍入误差）。  Cramer法则：只适于求解未知数个数极少时的情形。 若未知数个数为N ，则这种方法的计算次数近似地正 比于(N+1)! Cramer法则根本无法应用与数值计算。 1. 代数方程组求解方法概述 1.2 求解代数方程组的直接解法与迭代法 直接解法  Gauss消元法：先要把系数矩阵通过消元而化为上三 角阵然后逐一回代 若未知数个数为N ，该法所用乘法的次数近 3 似正比于N 。 即使是收敛很慢的 如二维正方形域内的Laplace方程，当采用 迭代方法也可能比 3200个节点时，采用Gauss消元法所需的计 直接法更有效！ 算次数是完成一轮Jacob 迭代所需计算次数 的168倍！ 1. 代数方程组求解方法概述 1.2 求解代数方程组的直接解法与迭代法 直接解法  求解三对角阵（一维问题）的直接解法：TDMA ， Thomas算法 TDMA 与Thomas 算法常用于五对角阵 （二维） 及七对角阵 （三维）迭代求解法中的直接求 解部分加快收敛速度 1. 代数方程组求解方法概述 1.2 求解代数方程组的直接解法与迭代法  迭代解法：假设解记作P = A-1b ，构造多维空间的 一个无限序列 P(n), 当n∞ 时，它收敛于A-1b 。一般 地，第n 次迭代所得之值取决于A ，b及上一次迭代值 P(n-1)，即：P(n) = f ( A, b, P(n-1) ) 研究： 1. 如何构造迭代方式； 2. 构造的迭代序列是否收敛； 3. 如果收敛，如何提高收敛速度。 1. 代数方程组求解方法概述 1.3 代数方程组迭代求解的终止判据 1. 规定迭代轮数 （全场的求解变量都作一次更新的计算过程为一轮次） 2. 相对偏差小于允许值 取决于计算问题、需要的精度及采用的迭代方法。如果采用的 迭