第9章常微分方程初值问题数值解法导论.ppt

当在某节点上xi的yi值有大小为δ的扰动时，如果在其后的各节点 上的值yi产生的偏差都不大于δ，则称这种方法是稳定的。 稳定性不仅与算法有关，而且与方程中函数f(x,y)也有关，讨论起来比较复杂。为简单起见，通常只针对模型方程 来讨论。一般方程若局部线性化，也可化为上述形式。模型方程相对比较简单，若一个数值方法对模型方程是稳定的，并不能保证该方法对任何方程都稳定，但若某方法对模型方程都不稳定，也就很难用于其他方程的求解。 先考察显式Euler方法的稳定性。模型方程 的Euler公式为 将上式反复递推后，可得 或 式中 要使yi有界，其充要条件为 即 由于 ，故有 （9.26） 可见，如欲保证算法的稳定，显式Euler格式的步长h的选取要受到式（9.26）的限制。 的绝对值越大，则限制的h值就越小。 用隐式Euler格式，对模型方程 的计算公式为，可化为 由于 ,则恒有 ,故恒有 因此，隐式Euler格式是绝对稳定的（无条件稳定）的（对任何h0）。 9.4.2 收敛性 常微分方程初值问题的求解，是将微分方程转化为差分方程来求解，并用计算值yi来近似替代y(xi)，这种近似替代是否合理，还须看分割区间 的长度h越来越小时，即 时， 是否成立。若成立，则称该方法是收敛的，否则称为不收敛。 这里仍以Euler方法为例，来分析其收敛性。 Euler格式 设 表示取 时, 按Euler公式的计算结果, 即 Euler方法局部截断误差为： 设有常数 ，则 （9.27） 总体截断误差 （9.28） 又 由于f(x,y)关于y满足李普希茨条件,即 代入(9.28)上式，有 再利用式（9.27），式（9.28） 即 上式反复递推后，可得 （9.29） 设 （T为常数） 因为 所以 把上式代入式（9.29），得 若不计初值误差，即 ，则有 （9.30） 式(9.30)说明,当时 , ,即 ,所以Euler方法是收敛的，且其收敛速度为 ，即具有一阶收敛速度。同时还说明Euler方法的整体截断误差为 ，因此算法的精度为一阶。 龙格-库塔方法是一类重要算法，但这类算法在每一步都需要先预报几个点上的斜率值，计算量比较大。考虑到计算yi+1之前已得出一系列节点 上的斜率值，能否利用这些已知的斜率值来减少计算量呢？ 这就是亚当姆斯（Adams）方法的设计思想。 §9.5 亚当姆斯方法 9.5.1 亚当姆斯格式 设用xi,xi+1两点的斜率值加权平均作为区间 上的平均斜率，有计算格式 （9.21） 选取参数λ，使格式（9.21）具有二阶精度。 将 在xi处Taylor展开 代入计算格式 化简, 因此有 与y(xi+1)在xi处的Taylor展开式 相比较, 需取 才使格式(9.21)具有二阶精度。这样导出的计算格式 称之为二阶亚当姆斯格式。 四阶亚当姆斯格式。 (9.22) 上述几种亚当姆斯格式都是显式的，算法比较简单，但用节点 的斜率值来预报区间 上的平均斜率是个外推过程，效果不够理想。为了进一步改善精度，变外推为内插，即增加节点xi+1的斜率值来得出 上的平均斜率。譬如考察形如 类似地可以导出三阶亚当姆斯格式。 （9.23） 的隐式格式,设(9.23)右端的 Taylor展开有 可见要使格式(9.23)具有二阶精度,需令 , 这样就可构造二阶隐式亚当姆斯格式 其实是梯形格式。类似可导出三阶隐式亚当姆斯格式 ( 9.10 ) 预测 校正 左矩形公式 梯形公式 改进的欧拉公式 右矩形公式 可以证明,改进的欧拉公式(9.10)的精度为二阶。这是一种一步显式格式,它可以表示为嵌套形式。 或者表示成下列平均化形式 改进的欧拉公式 左矩形公式 右矩形公式 左、右矩形公式的平均值，即梯形公式. 9.2.6 改进欧拉法算法实现 （1）计算步骤 ① 输入 , h , N