《计算方法》课件第七章 常微分方程数值方法.pptVIP

第七章 常微分方程数值方法 内容提要 §7.0 引言 §7.1 Euler 法 §7.2 Runge-Kutta法 习题七 第七章答案 §7.0 引言 §7.1 Euler 法 一、 Euler公式 二、梯形公式 三、改进的Euler公式 一、 Euler公式 例7.1 二、梯形公式 三、改进的Euler公式 例7.3 §7.2 Runge-Kutta法 一、Taylor级数法 二、Runge-Kutta法的基本思想 三、二阶Runge-Kutta公式 四、四阶Runge-Kutta公式 一、Taylor级数法 二、Runge-Kutta法的基本思想 三、二阶Runge-Kutta公式 习题七 第七章答案 这是含有4个未知数的、由3个方程所组成的方程组： 它有无穷多组解； 因而可以指定其中任意一个未知数。 恰好就是我们前面所得到的改进Euler公式。 返回引用 又称变形的Euler公式。 把系数满足式(7-19)的数值解公式，统称为二阶Runge-Kutta公式，它具有2阶精度。 返回节 按照上述的推导步骤，令 p = m = 4 ， 通过比较复杂的推导， 可以得到一族具有四阶精度的Runge-Kutta公式。 其中，最常用的标准四阶Runge-Kutta公式 四、四阶Runge-Kutta公式 由上面，可以看出： Runge-Kutta法是通过多计算 f 的函数值的办法， 避开了Taylor级数中的 f 的导数计算。 在这个意义上，可以说： Runge-Kutta法是Taylor级数法的一种变形， 而且四阶以下的Runge-Kutta公式，其级数与阶数(精度)是一致的(即m=p)。 但高于四阶的公式： 其级数要大于精度的阶数，即mp， 因而工作量大大增加。 因此，对于一般实际问题，四阶公式就够用了。 例7.4 用标准四阶Runge-Kutta公式，解初值问题 解： 先给出解本题目的具体标准四阶R-K公式。(下页) 取 h = 0.2 , 计算结果为 应当注意： 高阶R-K公式的推导是基于初值问题的解y(x)的Taylor展开， 因而要求y(x)具有较好的光滑性。 如果解 y(x) 的光滑性较差 虽然使用高精度的R-K公式， 有时甚至没有使用低精度的公式效果好。 返回章 返回节 1.用改进Euler公式和梯形公式求初值问题 的数值解，并求出其与精确解之间的误差限 ，计算到 。 2．用Euler公式和改进Euler公式求初值问题 的数值解，并求出其与精确解之间的误差限 ，计算到 。 3． 对 用下面方法求其数值解, (1)Euler公式；（2）改进Euler公式；（3）经典R-K公式。 返回章 1.095909 1.084096 1.266201 1.343360 1.416402 1.095656 1.083594 1.265441 1.342323 1.415059 1.095445 1.083216 1.264911 1.341641 1.414214 改进Euler公式 梯形公式 1. 2. 3. 返回章 * * （ 2学时 ） 《计算方法》课程讲义课件 大连海事大学 计算机科学与技术学院学院 在科学技术中，常常需要求解常微分方程的定解问题。这类问题中的最基本形式，就是本章讨论的一阶常微分方程的初值问题。 (7-1) 虽然求解常微分方程初值问题已经有了一些解析方法，但是： 这些解析方法只能求解某些特殊类型的方程。 因而，在实际中出现的许多有价值的问题，主要还得靠数值方法来求解。 返回引用 所谓数值方法： 不是寻求初值问题的解的解析表达式； 而是求出在一系列离散节点 a = x0 x1 ? ? ? ??? xn = b 的近似值 y0, y1, ….., yn 。 为了计算方便，可取 返回章 返回节 （7-2） 设初值问题（7-1）的解 把 返回引用 公式(7-3)叫做Euler公式