第七章(分方程-3).docVIP

PAGE PAGE 16 第七章 常微分方程的数值解法 §1　引言 1、一阶常微分方程初值问题（微分方程加初值条件） 　　　　　　　　　(1) 　　特征：为已知的二元函数，只有一个变元，一阶导数，要求满足微分方程且过点的曲线，一个未知函数（微分方程的解函数）。 2、解的存在性 若连续且对满足Lipschitz条件，即存在常数，使对，有，则初值问题(1)的解存在唯一，本章均假设对满足Lipschitz条件。 3、数值解法 问题：求精确解（解析表达式）极其困难，实际应用中只要求数值解。 数值解：在取一系列等距离散节点，步长，求的近似值。 … … … … ＝ ≈ ≈ … ≈ ≈ … … … 方法：建立的递推公式，，从而按节点排列顺序（步进式）。§2　Euler方法 显示、隐式Euler法与梯形法 对（1）的方程两端由到定积分得： 从而有：（*） （积分中值定理） （1）若取，得（左矩形公式） 由（*）得： ， ,则有： 称式为解初值问题的Euler法.（显格式）  　　几何意义：初值问题(1)的解曲线过点，从出发，以为斜率作一段直线，与直线交点于，显然有，再从出发，以为斜率作直线推进到上一点，其余类推，这样得到解曲线的一条近似曲线，它就是折线（某点导数值的几何意义为该点切线的斜率）。 （2）若取，得（右矩形公式） 类似可得：　 称式为隐式Euler法（向后Euler法，隐格式） （3）取（梯形公式） 则得：称式为梯形法.（隐格式）　 例1 用Euler法、隐式Euler法、梯形法及改进Euler法求解 ，取，计算到，并与精确解比较.解：由于,, Euler法： 　　　　　　 时，. 隐式Euler法：　　　　　解出　　 　当时， 梯形法： 　 　　 解出 　　当时， 改进Euler法: 即  当时， 　　精确解：，(梯形法效果最好，改进Euler较好：具有相同的误差数量级，其它不好!) Euler法 隐式Euler法 梯形法 改进Euler法 精确解 0.1 1.000 000 1.009 091 1.004 762 1.005 000 1.004 837 0.2 1.010 000 1.026 446 1.018 594 1.019 025 1.019 731 0.3 1.029 000 1.051 315 1.040 633 1.041 218 1.040 818 0.4 1.056 100 1.083 014 1.070 097 1.070 802 1.070 320 0.5 1.090 490 1.120 922 1.106 278 1.107 076 1.106 531 2、隐式法的计算（Euler法及梯形法） 方法1：显示化（对线性时，以为未知量的一元线性方程，见上例） 方法2：迭代法（对非线性时，可看作一个关于方程，利用迭代法求解） Euler法： （*1） （以为未知量的一元非线性方程） ,n （*2） 梯形法：（*3） ，（*4） 计算步骤： （1）初值 ； 迭代 （当） （2）；迭代 收敛性：只要步长足够小，就可保证迭代收敛。 （a）当（*2）收敛 （b）当（*4）收敛. 证明：（b）迭代序列收敛当 有；由（*3）、（*4）式知： （满足李氏条件） ，则当，（）时，有，从而收敛得证。 3、 改进Euler法 （梯形法隐式迭代式仅限迭代一次）,为避免迭代，可先用Euler法计算出的近似 　　　　显式预报，隐式校正 或　 （6） 称以上两式为改进Euler法（校正就是迭代一次） 右端已不含，显式方法。 例2：用梯形法的迭代格式求的数值解，，计算到 解：梯形公式 迭代式，解为：3、 Euler法的局部截断误差　　定义1 在的前提下，称为在的局部截断误差。 定义2　若一种数值方法的局部截断误差：，则称这种方法是精度为阶的。（含的项称为局部截断误差主项.） 注：按某种方法由算出这一步的误差。 （1）Euler法 （*） （**） 由(*)（**）式有： 根据定义2，Euler法中的故此方法为一阶方法. （2）隐式Euler法（分析时隐式公式右边的） （***） 由（*）及（***）式：　故局部截断误差主项为，也是一阶方法. （3）梯形法 类似有，方法是二阶的. （4）改进Euler法 则 = ,方法是二阶的.§3 Runge-Kutta方法 一、 显式 Runge-Kutta法的一般形式 一般形式： 二、 2级显式R-K方法（r=2） 设想构造R—K公式：（*）