第七章常微方程数值解.docVIP

PAGE PAGE 1 第七章 常微分方程数值解 §1 引言 一 一阶初值问题解的存在唯一性 一阶常微分方程初值问题 SKIPIF 1 0 （*） 其中 SKIPIF 1 0 是 SKIPIF 1 0 平面某一区域D上的连续函数，如果 SKIPIF 1 0 ，满足 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 存在，并满足方程 SKIPIF 1 0 那么 SKIPIF 1 0 是初值问题（*），在 SKIPIF 1 0 上的解。 对于（*）是否有唯一解？对 SKIPIF 1 0 还要附加一些条件。 定义 如果存在正常数 SKIPIF 1 0 ，使得对任意 SKIPIF 1 0 有 SKIPIF 1 0 则称 SKIPIF 1 0 满足Lipschitz条件，L称为Lipschitz常数。 如果 SKIPIF 1 0 ，那么有 SKIPIF 1 0 导数有界T满足Lipschitz条件 如果存在常数 SKIPIF 1 0 ，使得对一切 SKIPIF 1 0 及 SKIPIF 1 0 有 SKIPIF 1 0 则称 SKIPIF 1 0 对 SKIPIF 1 0 满足Lipschitz条件。同理，只要 f(x,y)对 y的偏导数有界，则f(x,y) 满足对y的Lipschitz条件。 定理（存在唯一性），设 SKIPIF 1 0 是在 SKIPIF 1 0 上的连续函数，而且 SKIPIF 1 0 对 SKIPIF 1 0 满足Lipschitz条件，则对任意 SKIPIF 1 0 ， 初值问题（*）在 SKIPIF 1 0 上存在唯一的连续可微解 SKIPIF 1 0 。 二 本章研究的问题 例1 y y(x)=e- 为了对初值问题进行求解，一些简单问题有解析解，大量非线性问题没有解析表达式，就是线性问题也不一定有解析解，因此，近似求解和数值求解常微分方程是非常必要的。 1 初值问题数值解基本概念 已知 y 满足（1.1）的解是过点（x0 首先对连续区间 SKIPIF 1 0 离散化 SKIPIF 1 0 为常数。 离散点是等矩的也可以是不等距的，下面仅讨论等距情况。 为方便起见： 把 SKIPIF 1 0 在xn处的精确值记为 SKIPIF 1 0 ，其近似值用 SKIPIF 1 0 表示 如何求xn处解的近似值yn？ 从X0 ,y x 写出计算公式有： y1 y2 ? yn+1 n=0,1,? 方程（1.2）叫差分方程 这种求微分方程近似解的过程称为步进式的方法，计算yn+1 x y y 局部截断误差，即xi 综上所述，求微分方程数值解需要处理以下几个问题： 1 把微分方程（连续的）离散化为差分方程（离散的） 2 用差分方程和初始条件计算出微分方程数值解 3 有关理论 误差，局部截断误差，主局部截断误差，阶 收敛性 稳定性（绝对稳定性） 本章讲授的内容： 1 单步法：显示Euler方法，隐式Euler方法，梯形方法，改进Euler方法，Rung-kutta方法 2 单步法的收敛性，稳定性，相容性 3 线性多步法 三 预备知识 1 一元Taylor多项式 y yx y y y 2 二元Taylor多项式 fx+h,y+k=f 1 其中： x+an= 3 数值积分 4 Lagrang插值 §2 简单数值方法 已知 y x （I）显式Euler方法（也称为Euler方法） Euler方法是求常微分方程初值问题的最简单办法。 上节的公式(1.2)就是显示Euler