高等工程数学第五讲,常微分方程数值解法.pptVIP

（4）四阶经典(Runge--Kutta) 公式 算法9-2（见教材212） 用四阶经典(Runge--Kutta)法求解下面初值问题的数值解。（取h=0.2,N=5） 依题意四阶经典(Runge--Kutta)公式应为 解： 例3、 由初值y0=0逐次算出下表 特点：用四阶经典(Runge--Kutta)法计算的结果，其精度高于改进Euler法. 例1. 解: 利用经典的四阶龙格-库塔公式有 * * 主要内容 10.1 Euler方法 10.2 改进Euler法 10.3 收敛性与稳定性 10.4龙格—库塔(Runge--Kutta)法 第五讲 常微分方程数值解法 , 许多实际问题的数学模型是微分方程或微分方程组的定解问题,如物体运动,电路振荡,化学反应及生物种群的变化等. 能用解析方法求出精确解的微分方程为数不多,有的甚至方程即使有解析解,也可能由于解得表达式非常复杂而不易计算.因此十分有必要研究微分方程数值解法. 目的：用数值解逼近解析解 我们考虑一阶常微分方程初值问题 （10.1） 设f(x,y)连续，且关于y满足Lipschitz条件，存在常数L,使得 由常微分方程理论知，初值问题（10.1）有解且唯一。 处得近似值 的方法， 称为问题（9-1）的数值解。 为步长。 常用几种方法:用差商近似导数；用数值积分方法；用Taylor多项式近似。 数值解法------求问题（9-1）的解y(x)在若干点 现将微分方程离散化----建立数值解法 (1)?? 用差商近似导数 可由初值y0逐次算出 .称离散化问题（9-3）为差分方程初值问题。 Euler方法----用差分方程初值问题的解近似微分方程初值问题（9-1）。 可由初值y0逐次算出. 用Euler方法求下面初值问题的数值解 10.1 Euler方法 例1、 解： 若取h=0.1,N=10,可由初值y0=0逐次算出.见下表 二,隐式欧拉公式 三,欧拉中点公式 若用梯形公式代替微分方程（10.3）右端积分，有 梯形公式的局部截断误差 10.2 改进的Euler方法 一,梯形公式 由于梯形公式仍然是隐式方法，一般用迭代法求解，若每一步只代一次，相当将Euler公式与梯形公式结合使用。 为上机编程方便，改写为 二, 改进Euler法(梯形法) 改进Euler法求解下面初值问题的数值解。（取h=0.1,N=10） 用公式（9-7） 若取h=0.1,N=10,可由初值y0=0逐次算出.见下表 例2、 解： 改进Euler法优于Euler法.(二阶方法) 10.3 收敛性与稳定性 截断误差 收敛性: （2）截断整体误差（不含舍入误差） 配项 结合前面有 稳定性 绝对稳定区域---- 绝对稳定区域包含左半平面----该方法称为A稳定。 模型 （1）Euler方法稳定性 由Euler方法知，计算公式为 10.3 收敛性与稳定性 绝对稳定区域---- 绝对稳定区域包含左半平面----该方法称为A稳定。 收敛性 则Euler方法的绝对稳定区域为 特点： o -1 -2 （2）梯形公式的稳定性 由梯形公式知 类似于前面梯形公式的方法 绝对稳定区域为左半平面。 o 3-2 (Runge--Kutta)法的构造 一般有 确定相应参数的原则：近似公式在（Xn,Yn) 处的泰勒展式与y(x) 在Xn处的泰勒展式的前面的项尽可能重合，以提高精度。 当P=2时的模型 上式在（Xn,Yn)处展开 比较系数有： ? 改进Euler公式 （二阶）中点公式 （3）三阶(Runge--Kutta) 公式