新-地下水流数值模拟讲义2007.pdf

《地下水流数值模拟》讲义，吉林大学，杜新强，2007 年 绪言 一、数值计算方法概述 1、数值计算的重要意义 17 世纪是数学发展史上一个划时代的时期，当我们今天享受着高科技成果 所带来的各种便利条件时，应该意识到笛卡尔（Descartes,1596-1650 ）和牛顿 （Newton,1642-1727 ）这两位卓越先驱人物所作出的具有划时代意义的贡献。 首先，笛卡尔创立了平面解析几何，使我们能够用数学形式描述动态变化着 的客观对象，接下来牛顿和莱布尼兹（Leibniz ）等人创立了微积分学，为我们研 究连续变量的变化规律给出了完整的方法体系。随后，又在这个基础上产生了更 多的数学分支以及相关的一些交叉学科分支。 由于数学研究的范围在不断扩大，而且研究的对象更为复杂，一个伴随的问 题就是相应的数值计算更加困难。如果不解决与理论方法平行的数值计算问题， 再好的数学理论也难以发挥应有的作用。 例 1：解线性方程组Ax ＝b ，x ∈Rn ，A 为 n 阶可逆方阵，用著名的 Crammer 法则求解，行列式的计 算按原始的方法计算，假设计算机每秒可算 1 亿次乘法运算，我们来估计机器所花的时间。 何为 Cramer 法则：该法则本身是十分重要的，它告诉我们线性代数方程组有解的充分必要条件以及 解如何求。 计算一个行列式所需要的乘法数：共 n ！项，每项n 个数相乘，故共需 n ！（n －1）次乘法；完成计算 的乘法数：共n ＋1 个行列式，故共需（n2 －1）n ！次乘法 行列式：将每一行、每一列在每一次都取一个数相乘，得到的结果再相加。 当n ＝10 时，需要 3592561200(次) ≈3.59251 （秒） 当n ＝100 时，需 9.33169 ×10161 （次）≈9.33169 ×10153 （秒）≈2.95906 ×10146 （年） 因此，当n 较大时，采用该方法进行计算是不可行的。 这表示在相应的计算机生产工艺条件下，一个好的数学方法（公式）未必是一个有效的计算方法。而 方法的有效性将直接影响到计算的效率，我们的任务就是要构造有效的数值计算方法。 例 2 ：设多项式为( p) x ( x2) − 9 ，我们来计算其在区间[1.92,2.08]上的值。 ( p) x ( x2) − 9 令： 1 《地下水流数值模拟》讲义，吉林大学，杜新强，2007 年 q x x −x +x − x + x − x ( +) x 18− x 144+ x 672− 9 8 7 6 2016 5 4032 4 5376 3 460823042 512 此例说明：即使数学上的恒等公式，用计算机来计算，结果也是不一样的。 科学的计算方法已成为第三种科学方法，而不仅仅是数学或计算机学科的一 个部分。数值计算的重要性一方面促进了计算方法的研究，另一方面也促进了计 算工具的发展。随着 20 世纪 40 年代中期人类第一台电子数字计算机的问世，数 值计算终于有了理想的支撑工具。 坚实的数学理论，科学的计算方法以及先进的计算工具的有机结合，可以为 我们打造了一个无比坚实的工作平台。完全可以这么讲，今天几乎所有的高科技