华科计算方法课件第六章.ppt

基本思路 三阶龙格-库塔公式 三阶龙格-库塔公式 为此，考虑格式的局部截断误差. 比较系数可知，若要格式具有二阶精度，只需 特别，若取 有 格式就是改进的Euler格式. 若取 有 格式称为变形的Euler格式（也称为中点格式）.其形式为 ? 三阶Runge-Kutta格式 问题：适当选取7个待定参数，使得格式具有三阶精度. 类似地，利用Taylor公式分析格式的局部截断误差.比较系数可知，若要格式具有三阶精度，只需 常用的取 * 计 算 方 法 华中科技大学数学与统计学院 第六章 常微分方程初值 问题的数值解法 计算方法课程组 §6.1 基本离散方法 §6.2 线性多步方法 §6.3 Runge-Kutta 方法 §6.4 收敛性与稳定性 §6 常微分方程数值解法 考虑一阶常微分方程的初值问题 : 例如： 其解析解为： §6.1 基本离散方法 但是, 只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。 因此，只能依赖于数值方法去获得微分方程的数值解。 例如： 其解析解为： 很难得到其解析解 例如： 其解析解为 只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。 因此，只能依赖于数值方法去获得微分方程的数值解。 要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 通常取节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。它适合计算机求解,应用广泛,具有应用价值。 6.1.2 Euler公式 利用等距分割，数值微分来代替导数项，建立差分格式。 1、向前差商公式 所以，可以构造差分方程 称为局部截断误差。显然，这个误差在逐步计算过程中会传播，积累。因此还要估计这种积累 定义:　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差。 定义: 　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 收敛性： 考察局部误差的传播和积累 2、向后差商公式 是隐格式，要迭代求解 可以由向前差商公式求出 3、中心差商公式 是多步，2阶格式，该格式不稳定 对微分方程 积分有: 类似，可以算出其误差估计式： 2阶的方法 所以，有 是个隐式的方法， 要用迭代法求解 局部截断误差 4、梯形公式 5、欧拉公式的改进： 隐式欧拉法 向后差商近似导数 x0 x1 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? ) 1 , ... , 0 ( ) , ( 1 1 1 - = + = + + + n i y x f h y y i i i i 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 欧拉公式，而前者称为显式 欧拉公式。 中点欧拉公式（欧拉二步法） 中心差商近似导数 x0 x2 x1 假设 ，则可以导出 即中点公式具有 2 阶精度。 需要2个初值 y0和 y1来启动递推 过程，这样的算法称为双步法 /* double-step method */，而前面的三种算法都是单步法 /* single-step method */。 方 法 ? ? 显式欧拉 隐式欧拉 梯形公式 中点公式 简单 精度低 稳定性最好 精度低, 计算量大 精度提高 计算量大 精度提高, 显式 多一个初值, 可能影响精度 改进欧拉法 Step 1: 先用显式欧拉公式作预测，算出 Step 2: 再将 代入隐式梯形公式的右边作校正，得到 Euler方法、隐式Euler方法、梯形方法与单步法计算公式的 显式单步法 对应关系 隐式单步法 显式 Euler方法 隐式Euler方法 梯形方法(隐式) 6.1.3 总结 算例：分别用Euler公式和改进的Euler公式求解： 取步长 ，计算y(0.5)的近似值 解：欧拉公式： 改进的Euler公式： 算例 分别用显式Euler方法，梯形方法和预估－校正Euler方法 初值问题 解： 取 h =0.1， (1) Euler方法为: 续 算例 分别用显式Euler方法，梯形方法和预估－校正Euler方法解初值