第六章 常微分方程的数值解法.pptVIP

第六章 常微分方程的数值解法 引子 常微分方程初值问题 解函数y(x), x0?x?xf. 解析法：常系数线性微分方程可用特征根法，一些特殊的微分方程可以用分离变量法等求解. 绝大多数常微分方程都难以求得精确解. 第六章 常微分方程的数值解法 §6.1 Euler法及其改进 §6.2 Runge-Kutta法 §6.3 收敛性与稳定性（略） §6.4 RKF法与Adams法（略） §6.5 微分方程组与高阶微分方程 §6.1 Euler法及其改进 1 数值解 取步长h，记xn=x0+nh，n=1,2, ?,N, xN=xf求得各节点xn解函数值y(xn)的近似值yn ，称y0，y1, …，yn, … yN为常微分方程的数值解。 2 Euler格式 Euler格式 ，n=0,1,2, ? 隐式Euler格式 ，n=0,1,2, ? 二步Euler格式 ，n=0,1,2, ? 原理(1)微分方程?数值微分 原理(2)微分方程?积分方程?数值积分 局部截断误差 设（*）为解常微分方程的差分格式，h为步长。假设y1, …，yn准确，称 为（*）的局部截断误差。当 ，称（*）为p阶格式或（*）有p阶精度。 注：当局部截断误差 时，“整体”截断误差为 Euler格式的局部截断误差 Euler格式 1阶格式 隐式Euler格式 1阶格式 二步Euler格式 2阶格式 证明方法：Taylor公式 隐式Euler格式 3 Euler法的改进 原理：微分方程?积分方程?数值积分?梯形公式 梯形格式 2阶精度 改进Euler格式 2阶精度 例6.2 h=0.2, x0=0, x1=0.2, x2=0.4, x3=0.6, x4=0.8, x5=1.0, Euler格式 ， n=0,1,2,3,4 隐式Euler格式 改进Euler格式 §6.2 Runge-Kutta法 原理：Lagrange微分中值定理 问题已转化为如何对K*进行近似计算。 Runge-Kutta法 Runge-Kutta法 2阶Runge-Kutta格式 求参数?1, ?2及p使 改进Euler格式 ( ) 中点格式 ( ) 4阶经典Runge-Kutta格式 例6.3 h=0.4 §6.3 收敛性与稳定性 （略） 收敛性定义: 称求解常微分方程的差分格式收敛, 如果对任意固定的x=x0+Nh, 当N?? (同时h?0), 数值解yN?y(x). P146例6.4. 定理6.1 设差分格式 p阶格式, 即局部截断误差为O(hp+1). 若增量函数 关于y满足Lipschitz连续性条件, 则数值解的(整体)截断误差为 。 证明 绝对稳定性 定义 设由某差分格式求试验方程y=?y(?0)数值解，若当yn有扰动?(数据误差或舍入误差)时，因此yn+1产生偏差的绝对值不超过|?|，则称该差分格式是绝对稳定的。 Euler格式绝对稳定域： 隐式Euler格式：任意h 改进Euler格式 : h-2/? 四阶R-K格式 : h-2.78/? ?§6.4 RKF法与Adams法(略) 怎样根据需要的精度来选定步长？ 变步长Runge-Kutta法 (i) 用上一步步长作为初始步长h； (ii) 计算y(h)和y(h/2)， 如果?=| y(h/2)? y(h)|/10 ?，继续折半步长直到???, 并将y(h/2)作为结果; 反之若???, 将步长加倍直到??，这时再将步长折半一次，就得到结果。 2 RKF(Runge-Kutta-Felhberg )格式 (Matlab ode45) 4阶Runge-Kutta格式和一个5阶Runge-Kutta格式组合