北京交通大学数值分析第七章 常微分方程数值解法（新）.ppt

三. R-K公式的具体构造——基于数值积分的角度 二级Runge-Kutta方法 取n=2 由此得 另一方面 为使局部截断误差为 ,应取 改进的Euler方法 取 中点方法 取 二阶Heun方法 取 二级Runge-Kutta方法不超过二阶 记 则 因此局部截断误差只能达到 三级Runge-Kutta方法 取n=3 记 又由于 因此要使局部截断误差为O(h4)，必须 Kutta方法 取 三阶Heun方法 取 三级Runge-Kutta方法不超过三阶 完全类似于二级Runge-Kutta方法的分析 将 和 都展开到 项 易证 三级Runge-Kutta方法的局部截断误差只能达到 四级R-K方法 取n=4 经典R-K方法 四. 变步长R-K方法 具体做法 变步长四阶标准R-K法解初值问题举例 五.再论单步法的稳定性 考虑一般的显式单步法． 若值y i有一个扰动?i，那么用显式单步法计算，得到的值yi+1就会产生一个偏差?i+1．若记 或 由此可知，单步法(35)式绝对稳定的条件是 （1）首先考察Euler方法的稳定性． 此时增量函数? (x, y,h) =f (x, y)=y，因而有 （2）讨论改进Euler方法的稳定性．此时增量函数 （3）经典Runge-Kutta方法的稳定性． 此时增量函数 代入后得 于是有 如果仅限于讨论?是实数的情形，则上述几个单步法的稳定性条件可分别简化为 Euler法稳定性条件：?2 ? h? 0， 改进Euler法稳定性条件：?2 ? h ? 0， 经典Runge-Kutta法稳定性条件： ?2.785 ? h ? 0. 例 对初值问题 取h = 0.1和0.2，用经典Runge-Kutta方法求解. 隐式单步法的稳定性讨论. （1）考察向后Euler方法. 对于试验方程(46)，其向后Euler法的公式为 yi+1 = y i + h? yi+1 解出yi+1，有 ，从而得到误差(扰动)公式为 由此得到绝对稳定的条件为 其绝对稳定区域是以1为半径、以1为中心的圆外部. （2）讨论梯形法的稳定性． 对于试验方程(46)，相应的梯形法公式为 解出yi+1，有 从以上的分析讨论可以看到，隐式方法的稳定性比显式方法好，这也是隐式方法的主要优点！ §4 线性多步法 求解初值问题的数值方法都是“步进式”的，即求解过程从初值y0开始，顺着节点的排列次序，一步一步地向前推进．所以，在计算yi+1 时，前面的i + 1 个值y0, y1, ?, y i 都是已知的．如果在计算yi+1 时能充分利用这些已有的信息，而不是像单步法中那样，只用其前一步的值y i，则可望构造出精度高，但计算量小的求解公式．线性多步法就是基于这一思想发展起来的，其计算公式可表示为 其中， ，而?j ,?j 都是常数. 线性多步法(50)式的实质是用若干节点处的函数值及导数值的线性组合来逼近y(xi+1)的值．由于在计算yi+1时需要用到其前k + 1个值：y i , yi?1, yi?2,?, y i ?k，故多步法(50)式又称为k + 1步法，且当??1 = 0时，此k + 1步法是显式的，而在??1? 0时，此k + 1步法是隐式的. 原则上，一切形如(50)式的多步法都可用Taylor级数展开的方法来导出(即确定其中的系数?j ,?j )，但有些多步法也可用数值积分法来构造．下面讨论实际中较为常用的几种线性多步法. 一、Adams外插法 1．公式的推导 对方程y0 = f (x, y)的两边从xi到xi+1积分，得 为了近似计算(51)式中的积分，我们以xi ?k , xi?k+1,?,xi?1, xi 为插值节点，作函数f (x, y (x)) 的k 次插值多项式 p k (x)，从而有 f (x, y (x) ) = p k (x) + R (x), 其中，R (x)为插值余项．将上式代入(51)式，得 略去积分余项R (x) d x，并用y i代替y (xi )，可得到计算公略去积分余项R (x) d x，并用y i代替y (xi )，