第七章节 常微分方程数值解.doc

第七章 常微分方程数值解 §1 引言 一 一阶初值问题解的存在唯一性 一阶常微分方程初值问题 （*） 其中是平面某一区域D上的连续函数，如果，满足 存在，并满足方程 那么是初值问题（*），在上的解。 对于（*）是否有唯一解？对还要附加一些条件。 定义 如果存在正常数，使得对任意有 则称满足Lipschitz条件，L称为Lipschitz常数。 如果 ，那么有 导数有界(满足Lipschitz条件 如果存在常数，使得对一切及有 则称对满足Lipschitz条件。同理，只要 f(x,y)对 y的偏导数有界，则f(x,y) 满足对y的Lipschitz条件。 定理（存在唯一性），设是在 上的连续函数，而且对满足Lipschitz条件，则对任意， 初值问题（*）在上存在唯一的连续可微解。 二 本章研究的问题 例1 ，满足微分方程和初始条件的解是 y(x)=，无法给出具体表达式。 为了对初值问题进行求解，一些简单问题有解析解，大量非线性问题没有解析表达式，就是线性问题也不一定有解析解，因此，近似求解和数值求解常微分方程是非常必要的。 1 初值问题数值解基本概念 已知 (1.1) 满足（1.1）的解是过点（的一条曲线y(x) 首先对连续区间离散化 为常数。 离散点是等矩的也可以是不等距的，下面仅讨论等距情况。 为方便起见： 把在处的精确值记为，其近似值用表示 (1.2) 方程（1.2）叫差分方程 这种求微分方程近似解的过程称为步进式的方法，计算若用到前面不止一个信息量，叫多步法。 ，即这一步的公式误差。 综上所述，求微分方程数值解需要处理以下几个问题： 1 把微分方程（连续的）离散化为差分方程（离散的） 2 用差分方程和初始条件计算出微分方程数值解 3 有关理论 误差，局部截断误差，主局部截断误差，阶 收敛性 稳定性（绝对稳定性） 本章讲授的内容： 1 单步法：显示Euler方法，隐式Euler方法，梯形方法，改进Euler方法，Rung-kutta方法 2 单步法的收敛性，稳定性，相容性 3 线性多步法 三 预备知识 1 一元Taylor多项式 y 2 二元Taylor多项式 k)f(x,y)+ , = 3 数值积分 4 Lagrang插值 §2 简单数值方法 (2.1) （I）显式Euler方法 Euler方法是求常微分方程初值问题的最简单办法。 上节的公式(1.2)就是显示Euler方法 (2.2) 问题： 要从(2.1) 离散化得到(2.2)还有什么办法？ 1数值积分方法 把（2.1）写成如下形式： (**) 用左矩形公式近似左边积分,用得到结果 (2.2) 2 Taylor展开方法 略去高阶项 即 设是初值问题的解，那么有 从而有 用近似值由此得出公式（2.2） 3数值微分方法 称为差商，即用差商近似微商， 也得出： 由此得出公式（2.2） 当已知时，可由公式（2.2）简单地求出，方法为显式的，由上的近似值可求出上的近似值，称为单步公式， （2.2）称为显式Euler方法也称为Euler方法。 （II）隐式Euler方法（也叫后退Euler方法） 数值积分方法 对(**) 公式右边积分用右矩形积分公式得到 (2.3) 也可用Taylor展开方法 略去高阶项，并设是初值问题的解，则有 即写为 同样用 的近似值代入有 (2.3)的右端含有，一般，不能直接由（2.3）得出，这种方法是隐式的。（2 .3）称为隐式Euler方法。 该方法也可用数值微分方法得出。 （III）梯形方法 在上对上式进行积分有 等式右边积分采用梯形公式近似有 用来代替就得到 （2.4） （2.4）称为梯形公式 (IV) 用公式做计算 显示Euler公式由，但是，对于（2.3）、（2.4）是隐式公式，不能由直接计算出，而是要解方程，一般用迭代方法，以（2.3）为例， 取，或用显式公式求出作为 ， 即 。 当 时，取 这样方法称为迭代法 下面考虑迭代收敛性： 当收敛；称为迭代收敛条件 对于梯形公式（2.4），由于等式右边含有，因而是隐式方法。 用它们来求时必须解方程，一般用迭代求解。 取， 迭代公式为 同样，当时，取；仿隐式Euler