chap7常微分方程及数值解法.ppt

第7章 常微分方程（组）的数值解法 刘东毅 天津大学理学院数学系 第7章 常微分方程（组）的数值解法 主要目标： 掌握常微分方程初值问题数值解法的基本理论 掌握计算机上的常用算法 主要内容： 初值问题计算格式的建立 Runge-Kutta 方法 一阶常微分方程组与高阶方程的数值解法 第7章 常微分方程（组）的数值解法 在科学研究和工程实践中会遇到很多微分方程，虽然从理论上可以证明其解的存在性，但其解的解析表达式往往是很难求解的，或者即使可以写出来，但也难于计算，此时，只能借助数值解来解决问题. 常微分方程（组）定解问题是自然科学和工程技术领域中常见的数学模型. 本章介绍求解此类问题的基本理论和数值解法。 下面，我们给出常微分方程初值问题解的存在惟一性定理。 定理 假设 f (x, y)∈C(D),且关于 y 满足 Lipschitz 条件, 则一阶常微分方程初值问题(7.0.1) 存在唯一解. 接下来，在此前提下, 我们讨论上述初值问题 (7.0.1) 的数值解法。 其基本思想是在区间 [a, b] 上引入一系列节点 7.1 初值问题计算格式的建立 略去余项, 并以数值解 yn , yn+1 代替 y (xn) 及 y (xn+1), 则得差分方程 如果利用下列数值微分公式 如果利用下列三点数值微分公式 2. Taylor 展开法 3. 数值积分方法 若对此式的积分采用数值积分梯形公式, 隐式梯形公式的迭代格式 定理 7.1.1 设 f (x , y) 在区域 D 上关于 y 满足 Lipschitz 条件, 即 例7.1.1：利用Euler公式与改进的Euler公式求解初值问题（步长h = 0.1） 计算结果见下表(见书P227表7.1) 解此问题的改进的Euler公式为 7.1.2 截断误差与方法的精度 整体截断误差 en+1 是在没有引进舍入误差的情况下, 纯粹因为不准确的计算格式造成的, 故又称为方法误差.它不仅与 x = xn+1 这一步的计算有关, 而且和 xn , xn-1 ,. . . , x1 这几步的计算都有关系. 如Euler公式 7.2 Runge-Kutta 方法 进一步假设 y∈Cp+1[a, b],由 Taylor 公式有 7.2 Runge-Kutta 方法 在(7.2.3)中略去余项,用yn, yn+1代替y(xn)及y(xn+1) 由于局部截断误差 7.2.1 Runge-Kutta方法的基本思想 其一般形式为: 7.2.2 二阶 Runge-Kutta 方法 得到 再根据 y(xn+1) 在点 xn的一元 3 阶 Taylor 展开式 即由 7.2.3 四阶Runge-Kutta方法 7.5一阶常微分方程组与高阶方程初值问题的数值解法 7.5.1 一阶常微分方程组初值问题 写成向量形式： 利用向量值函数的微积分理论, 很容易推导出一阶常微分方程组初值问题的数值解法. (7.5.3)的分量形式为 四阶标准的Runge-Kutta公式 四阶标准的Runge-Kutta公式的分量形式 例7.5.1：试写出用中点公式解下列初值问题的计算公式： 中点公式的向量形式 整理得分量计算格式 7.5.2 高阶常微分方程初值问题的数值解法 高阶常微分方程初值问题的一般形式为 等价的一阶常微分方程组初值问题： 例7.5.2写出用标准四阶Runge-Kutta公式求解 求此问题数值解的标准四阶Runge-Kutta公式为 7.5.3 高阶常微分方程组初值问题的数值解法 对于高阶常微分方程组的初值问题, 考察下面的二阶常微分方程组初值问题: 等价的一阶常微分方程组初值问题 本章小结 基本知识 建立初值问题数值解的计算格式的方法 数值微分方法、Taylor 展开法和数值积分方法 截断误差与方法的精度 Runge-Kutta方法 基本思想 一阶、二阶、四阶Runge-Kutta方法 一阶常微分方程组数值解法 高阶方程（组）的数值解法 考虑 r = 2 的情况, 此时有 利用二元函数的一阶Taylor 公式，即全微分公式 希望适当选择参数 , 使上式的局部截断 误差为 , 即为二阶方法. 下面将yn+1与y(xn+1)作比较 从而有 将上式代入 再由 在下面要将yn+1与y(xn+1)作比较，使它们的局部截断误差满足 为此考虑y(xn+1)。 由刚才已得到的 让它们满足 左式含有四个未知元三个 方程, 因此解不唯一. 参数 满足左式的一族公式统称 二阶 Runge-Kutt