计算方法9.3.ppt

9.3 误差估计与稳定性 * * 常微分方程初值问题的求解，是将微分方程转化为差分方程来求解，并用数值解yi来近似替代精确解y(xi)，这种近似替代是否合理，还与分割区间[xi-1,xi]的长度有关，当步长h越来越小，即h=xi-xi-1→0时，判断yi→y(xi)是否成立。若成立，则称该数值方法是收敛的，否则称为不收敛的。 * * 在推导欧拉公式或梯形公式时，利用数值微分代替导数项，得到差分方程。 称为局部截断误差。显然，这个误差在逐步计算过程中会传播，积累。因此还要估计这种积累 一、误差与收敛性 　　　 若某数值方法的局部截断误差为O(hp+1)，则称这种数值方法的阶数是p ，或称该数值方法具有p 阶精度。 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，用某种数值方法计算yi+1，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为该数值方法计算yi+1的局部截断误差。 * * 定义9.1 定义9.2 讨论常微分方程的数值解法的误差，主要分析求解公式局部截断误差和整体截断误差，并引入阶数的概念。 步长(h1) 越小,p越高, 则局部截断误差越小,计算精度越高。 * * 1. 欧拉公式的局部截断误差 设f(x,y)充分光滑,将y(xi+1)在xi点作Taylor展开: 欧拉法具有 1 阶精度。 * * 2. 梯