工程数值方法上(2015-3-17讲稿).doc

工程数值方法学习内容： Chapter 1 线性代数方程组的数值解法 Chapter 2 插值问题与数值微分 Chapter 3 数值积分方法 Chapter 4 常微分方程(组)初值问题的数值方法 Chapter 5 常微分方程(组)边值问题的数值方法 Chapter 6 椭圆型偏微分方程的数值方法 Chapter 7 加权残值方法 参考书目： 武大山大编计算方法，高教版 林成森，数值计算方法（上下），科学出版社，2000 中科院研究生数学丛书，工程中的数值方法，科学出版社，2000 曾绍林，工程数学基础（研究生数学丛书），科学出版社，2001 李庆扬，数值分析基础教程，高等教育出版社 李庆扬，数值分析（第4版），清华版 关治，数值计算方法，清华版 李岳生黄有谦，数值逼近，高教版 李荣华编，微分方程数值解法，人教版 邱吉宝，加权残值法的理论与应用，宇航版 Chapter 1 线性代数方程组的数值解法 线性代数方程组的求解是工程实践中最常遇到的问题。据不完全统计，在工程实践中提出的计算问题中，有近一半涉及到求解线性方程组。例如结构有限元分析问题，大地测量问题，气象预报问题，电力传输网分析问题，各种电路分析问题，数据拟合问题，以及非线性方程组与微分方程的数值求解问题等等。因此，学习掌握线性代数方程组求解的基本理论方法是十分必的。 本章将介绍目前一些利用计算机求解线性代数方程组常用的且简单有效的数值方法。 求解线性方程组的数值方法尽管很多，但归并起来可为两类： 直接法（精确法）经有限次的四则运算若运算中没有舍入误差即可求得方程组精确解的方法。如克莱姆（Cramer）则方法、消元法、分解法、分解法等等。 迭代法（）将求解方程组的问题化为构造一个无限迭代的序列，在实现该序列过程中的每一步计算结果是把前一步所得结果施行相同的计算步骤而获得，这一无限序列的极限就是原方程组的精确解答。如：简单迭代法、赛德尔迭代法、牛顿法、共轭斜量法等。 需要指出的，在一般情况下，我们使用直接法和迭代法两类方法都不可能完全获得原方程组的精确解答。显然（1）实际中在使用直接法时不可能没有数值计算的舍入误差，故所谓精确方法的解并不是绝对精确的（2）实际中在使用迭代法时，不可能极限过程进行到底，而只能进行有限次的迭代，获得满足精度要求的近似解答。 关于这两类方法求解的误差分析，我们将在每类方法的介绍之后进行简要讨论。 §1.1 直接法—Cramer法则与求逆方法 n元n非齐次线性代数方程组为： (1.1) 利用矩阵和向量符号，这个方程组可表为: (1.2) 其中,——方程组（1.1）系数矩阵（阶方阵）； ——方程组（1.1）的右端向量（阶列阵）； ——方程组（1.1）解向量（阶列阵）。§1.1.1 Cramer法则 由线性代数中线性方程组解的定理可知： 若A的行列式，则方程（1.1）有唯一解此时，根据Cramer法则，其解的表达式为： 1.3) 其中， 即将中的第j列元素依次换为右端向量的元素所构成的阶方阵的行列式。 易见，利用Cramer法则给出的（1.3）式求解n阶线性方程组（1.2）式，需要计算个阶行列式，而每个n阶行列式将有项，其中每一项又含有n个因子，故展开每一n阶行列式仅乘法运算就需次（忽略加法运算次数）而对个n阶行列式其乘法运算次数为： 对式（1.3）其除法次数为：n 解方程（1.2）乘除法运算次数为：次 例如，阶，次 这是一个十分惊人的数字，即使利用超高速的电子计算机能够胜任此计算次数，仅由于多个数的连乘亦有可能造成溢出而无法继续运算。因此，Cramer法则这个在理论上完善且精确的求解方法，仅在理论上和一些特殊情况下可以发挥作用，而高阶线性代数方程组的实际求解中几乎没有多少实用价值。为此，人们不得不研究其它一些计算简单且行之有效的求解方法。 §1.1求逆方法 若，则A非奇异，存在，则方程组（1.2）的解向量可表为： （1.4） 矩阵求逆方法有：（1）伴随矩阵法；（2）初等变换法； 当A的阶数n较大时，求逆麻烦且计算量大。高阶线性代数方程组求逆方法 §1.2 直接法—Gauss（高斯）消去法 尽管是一种古老的方法，但仍不失为最常用和最有效的方法之一。 基本思想：逐次消元将原方程组化为等价的三角形方程组进行求解。 §1.三角形方程组 所谓三角形方程组是以下两种形式的方程组(1)下三角形式 (2.1)矩阵记为 (2.1') 其中系数