第八章 常分方程的数值解法.ppt

§1 欧拉（Euler）方法 列表分别进行计算，结果如下 欧拉公式 欧拉预估-校正公式 0 0.0 1 1 1 0.1 0.9 0.91 2 0.2 0.82 0.8368 3 0.3 0.75676 0.77858 4 0.4 0.70849 0.73435 5 0.5 0.67430 0.70364 §2 龙格-库塔（Runge-Kutta）方法 龙格-库塔方法基本思想 利用微分中值定理，存在一点 使得差商 又 ，所以 其中 称为 上的平均斜率。显然，只要 提供平均斜率的一种近似，便可导出相应的计算格式。 §2 龙格-库塔（Runge-Kutta）方法 欧拉方法：应用点 的斜率 作为平均斜率的近似值 改进欧拉方法：应用 、 两点斜率 、 的算术平均近似平均斜率 基本思想：在区间 内设法多预报几个点的斜率值，并把它们进行加权平均作为平均斜率的近似值，则可能构造精度更高的计算格式。 §2 龙格-库塔（Runge-Kutta）方法 二阶龙格-库塔公式 在区间 内任取一点 ，用 和 两点斜率 的加权平均得到平均斜率，所以令 其中 和 为待定常数。 §2 龙格-库塔（Runge-Kutta）方法 用欧拉公式预报 的值，则 则计算公式为 希望通过选择 和 使计算格式具有较高精度。 §2 龙格-库塔（Runge-Kutta）方法 将 在点 泰勒展开有 §2 龙格-库塔（Runge-Kutta）方法 所以有 由 在点 泰勒展开有 比较 和 在 的泰勒展开式可知，当 时，前 面公式具有二阶精度，称为二阶龙格-库塔公式。 第八章 常微分方程的数值解法 实际问题中，经常需要求解常微分方程。只有一些特殊 类型的常微分方程有解析解，实际中的大多数常微分方程无 法求得解析解。在多数的工程应用中，往往只需要获得若干 点上的函数值，即求数值解。 包含自变量、未知函数和未知函数的导数或微分的方程 称为微分方程。在求解微分方程时，必须附加某种定解条件。 微分方程和定解条件一起组成定解问题。定解条件通常有两 种，一种是给出积分曲线在初始时刻的性态，这类条件称为 初始条件，相应的定解问题就是初值问题；另一种是给出了 积分曲线首末两端的性态，这类条件称为边界条件，相应的 定解问题称为边值问题。 第八章 常微分方程的数值解法 未知函数为一元函数的微分方程称为常微分方程，未知 函数为多元函数，从而有多元函数偏导数的方程称为偏微分 方程。微分方程中各阶导数的最高阶数称为微分方程的阶。 本章讨论一阶常微分方程的初值问题。 由常微分方程基本理论知道：当 在矩形区域 上连续，并且关于变量 满足李 普希兹（Lipschitz）条件，即对任意的 ， 都有 其中 为常数，则一阶常微分方程初值问题的解存在且唯一。 第八章 常微分方程的数值解法 求解常微分方程有各种各样的解析方法，但是解析方法 只能求解一些特殊类型的方程，求解实际应用中得到的微分 方程主要利用数值解法。 对于一阶常微分方程初值问题 求其数值解，就是计算出解函数 在离散节点 处的近似值 。 第八章 常微分方程的数值解法 相邻两个节点的距离 称为步长，假设节点是 等距的，即 为常数，这时 对应的函数值为 这样把一个连续问题转化为一个离散型问题的过程称为 离散化过程。离散化过程是把连续的微分方程初值问题转化 为一个离散的差分方程问题，然后将求得的差分方程初值问 题的解 ，作为微分方程的解 在 处的值 的近 似值。这种离散变量法又称为差分法，对应的离散方程称为 差