第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越来越小时，