华北电力大学在职研究生数学_教学案例设计_教学研究_教育专区.ppt

Mathematical Methods its Applications J. G. Liu School of Math. Phys. North China Elec. P.U. 常微分方程数值解法 §1 问题的提出和基本概念 1、问题的提出; 2、常微分方程初值问题的数值解法; 3、一些基本概念 §2 单步法 1、Euler法，梯形公式，改进Euler法 A）设计思想： （1） 差商代替导数； （2） 数值积分法； （3） 泰勒展开法。 B）改进的Euler公式 2、龙格-库塔（Runge-Kutta）法 A) 龙格-库塔法的设计思想； B) 龙格-库塔法的设计方法； C) 常用龙格-库塔法。 C ）总结 1、问题的提出 （1） 设初值问题（1）的解存在且唯一，求解（1）会遇到如下问题： A）（1）的解析解无法表出； B）（1）的解析解可以求出，但是计算其函数值很复杂。 如： 其解析解为 返回 2、常微分方程初值问题的数值解法 给定步长h0，取节点 通过数值方法计算问题(1)的解y(x)在各个节点处的近似值 该数值方法称为初值问题(1)的数值解法，近似解yk称为数值解。 常微分方程数值解法的特点 从初值出发，顺着节点排列次序一步步向前推进，即利用已知 信息 ，来计算 。 ——步进式 返回 3、一些基本概念 （1） 单步法和多步法 单步法：计算yk+1时，只用到前一节点的信息；（自开始） 多步法：计算yk+1时，用到了前面一个以上节点的信息； （2） 显式格式和隐式格式 显式格式：若 yk+1的计算表达式只含有已知信息；（便于计算） 隐式格式：yk+1的计算表达式含有未知信息。（数值稳定性好） 一般形式: 增量函数 （3） 截断误差 若计算数值解yk+1时用到的信息都是准确的，称 为该方法的局部截断误差， 返回 并称该方法为p阶方法。 （1） 差商代替导数 1） 向前差商 ——Euler公式（显式格式） （2） 2） 向后差商 （3） ——隐式Euler公式 3） 中心差商 ——Euler两步法, 显式格式 （4） 返回 （2） 数值积分法 1） 左矩形公式 ——Euler公式 2） 右矩形公式 ——隐式Euler公式 注：梯形公式即为Euler公式和隐式Euler公式的算术平均！ 3） 梯形公式 ——梯形公式 （5） 4） 中矩形公式 ——Euler两步法 返回 （3） 泰勒展开法 1）将y(x)在x=xk处展开 令x=xk+1，则 只保留h的线性项，得 —Euler公式 局部截断误差为（假设yk准确） 因此，为一阶方法！ （6） 主项 2）将y(x)在x=xk+1处展开 令x=xk，则 只保留h的线性项，得 —隐式Euler公式 为一阶方法！ （7） 局部截断误差为（假设yk准确） 3）（6）-（7）可得 保留h的线性项，得 ——梯形公式，隐式格式 （8） 二阶方法！ 局部截断误差（假设yk准确）： 4）将y(x)在x=xk处展开，得 保留h的线性项，得 ——Euler两步法，显示方法，二阶方法！ （9） （6）-（9），得 返回 改进的Euler公式 由Euler公式和梯形公式结合构成预报校正系统： 返回 总结 显式欧拉法 1、显式欧拉法是一阶收敛的算法； 2、显式欧拉法计算简单； 3、显式欧拉法对函数的光滑性要求不高； 4、显式欧拉法也形象地称为折线法。 隐式欧拉法 1、隐式欧拉法是一阶收敛的算法； 2、隐式欧拉法每步需解非线性方程； 3、隐式欧拉法对函数的光滑性要求不高。 收敛阶的提高需要更多函数值计算！ 返回 1、龙格-库塔法的设计思想 ——平均斜率 Euler公式 隐式Euler公式 龙格-库塔法是一类高精度的单步法! 由牛顿-莱布尼兹公式，我们知道 通过对K*取不同的近似，可以得到不同的计算公式！例如 龙格-库塔法的基本设计思想： 在[xk,xk+1]内多预测几个点的斜率，然后取其加权平均作为K*的近似，以构成具有较高精度的计算公式！ 改进的Euler公式 返回 2、龙格-库塔法的设计方法 （1）一级龙格-库塔法(1阶方法) 其一般形式为 确定参数?使得其收敛阶最高！ 由y(xk+1)在x=xk处的泰勒展开式 可知方法(10)的局部截断误差为 (10) 可见?=1时，(10)的局部截断误差最小，为O(h2)。 所得一阶龙格-库塔公式为 ——