第三章 线性代数方程组的数值解.ppt

一般地说，差分格式被写成如下的形式： 方程组（3.1）可被进一步写成矩阵的形式 : 解线性代数方程组的直接法。 大家知道，线性方程组 显然，按Cramer法则求解方程组（3.3）需要计算n+1个n阶行列式。每个n阶行列式按直接展开办法来算，需作（n-1）×n！次乘法和n!次加法运算。当n=30时，共约需完成 次乘法和加法运算，这是一个十分惊人的数字，即使在一台每秒作一亿次运算的计算机上完成这一计算也是不可能的。所以，尽管这种办法也是一种直接法，并且理论上可行，但实际上是无法进行求解的。即使采用其它办法来计算行列式，按Cramer法则求解的工作量也比通常的直接法大得多。因而，Cramer法则对于数值计算来说是没有什么用处的，仅在一些特殊场合才有用。 另外，矩阵求逆也是人们熟悉的求解线性方程的一种方法。 现在计算机上常用的直接解法大多数是以系数矩阵的三角形化为基础的。就是说，先对方程组进行变换，使其化为等价的（即具有相同解的）三角形方程组。由于三角形方程组的求解十分容易，原方程的求解问题即告解决。 §2 三角形方程组的解法 所谓三角形方程组是指下面两种形式的方程组： 若用矩阵符号可分别写为： 三角形方程组的求解是很简单的。方程组(3.6)的计算公式可归结为： 对于方程组(3.7)，其计算公式可归结为： §3 Gauss消去法 Gauss消去法（简称消去法）的提出已有相当长的时间了，是一种古老的方法。然而，近年来在计算机上求解线性代数方程组的实践表明，它仍是直接法中最常用的一种方法，也是最有效的方法之一。其基本思想是：用逐次消去一个未知数的办法把原来的方程组化为等价的（具有相同解）三角形方程组。这样，求解就很容易了。 假定把要求的n 阶线性的方程组(3.1)改写成如下形式： 其中 为 方阵， 为 向量，它们分别为： 这时方程组(3.10)就变为如下等价方程组： 表示成矩阵形式为： 若令 ，则系数 和 的计 算公式应为： 其中 上述的消去步骤还可以进行下去。如此继续之，重复上述步骤(n -1)次以后，我们即可得到如下等价三角形方程组： 表示为矩阵形式： 三角形方程组 很容易用前述的回代过程(3.9)求解，这样就完成了消去法求解n阶线性代数方程组的过程。从原来方程组(3.10)得出等价三角形方程组(3.11)的过程称之为消去过程。采用前面的记号，我们可将消去过程的计算公式归结为对于 ，递推地计算如下各量： 对照上节中差分格式，不难看到，由有限差分法得到的系数矩阵是能够满足以上两个条件的。因为每个节点方程（第i个方程）都代表着一个单元（第i个单元）与周围单元或外界环境的热量交换关系。在这些热量交换中，无疑都与该单元的温度(Ti)有关，Ti 的系数aii 当然不能为零。 而且在第i个方程中，Ti 的地位比其周围单元温度更为突出，表现在系数上，它的绝对值总是最大的。对于给定温度的节点方程而言，这种性质更明显。因为方程中除了该节点温度以外再也没有别的节点温度，它的系数当然也就最大了。 综上所述，由于对稳定导热问题用有限差分法得到的代数方程具有上述性质，因此在求解方程组时，可大胆放心使用Gauss消去法。（对于更一般的方程组，目前更多采用主元素消去法，而不用Gauss消去法。） §4 迭代法 前面介绍的解线性代数方程组的直接法对于阶数不是很高的问题是非常有效的，这种场合一般不使用下面介绍的迭代法。 迭代法的优点：由于不需要存储系数矩阵的零元素，所以占用的存储单元少。同时程序也比较简单，对于稳定导热用有限差分法所得到的方程组，求解收敛较快，因此广泛地被采用。 迭代法所要讨论的问题，仍是如下线性代数方程组： 迭代法的基本思想: 构造一个由 组成的向量序列，使其收敛于某个极限向量 ，并且 就是方程组(3.14)的精确解。 1．简单迭代法 最简单的迭代法称为简单迭代法，也称同步迭代法，Jakobi迭代法。迭代的最终目的是求解方程组(3.14)中的 。 则可将式(3.14)改写成： (1). 任意给定各节点上的温度值 作为解的第零次近似，把它们代入式(3.17)的 右端，由此算得 为解的第k+1次的似。这样得到的序列