数值分析-第八章-常微分方程数值解法.pptVIP

ISCM 2007,Beijing China 数值分析Numerical Analysis 第八章 常微分方程数值解法 记数值误差为： 引进试验方程： 向前欧拉公式的稳定性 试验方程的欧拉公式 若每步计算有舍入误差，则 当λ为复数时 当λ为实数时 本章小结 例：考察初值问题 在区间[0, 0.5]上的解。 分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。 0.0 0.1 0.2 0.3 0.4 0.5 精确解 改进欧拉法 欧拉隐式 欧拉显式 节点 xi 1.0000 ?2.0000 4.0000 ?8.0000 1.6000?101 ?3.2000?101 1.0000 2.5000?10?1 6.2500?10?2 1.5625?10?2 3.9063?10?3 9.7656?10?4 1.0000 2.5000 6.2500 1.5626?101 3.9063?101 9.7656?101 1.0000 4.9787?10?2 2.4788?10?3 1.2341?10?4 6.1442?10?6 3.0590?10?7 §8.4 单步法的稳定性 §8.4 单步法的稳定性 定义 　　　若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的 /*absolutely stable */。 一般分析某算法的稳定性时，为简单起见，只考虑模型方程或试验方程 /* test equation */ §8.4 单步法的稳定性 当h=0.25时 0 0.5 0.75 1.0 1 0.25 h = 0.5 h = 0.25 §8.2 欧拉(Euler)法 欧拉方法的收敛性 §8.2 欧拉(Euler)法 局部截断误差 称为局部截断误差 §8.2 欧拉(Euler)法 欧拉方法的收敛性 定义 若给定方法的局部截断误差满足 则称该方法是 P 阶的，或称为具有 P 阶精度。 §8.2 欧拉(Euler)法 整体截断误差 §8.2 欧拉(Euler)法 欧拉方法的收敛性 由此知，当 §8.2 欧拉(Euler)法 注 整体截断误差与局部截断误差的关系： §8.2 欧拉(Euler)法 向后欧拉公式 隐式欧拉法或向后欧拉法 /* implicit Euler method or backward Euler method*/ xn+1点向后差商近似导数 隐式或后退欧拉公式 §8.2 欧拉(Euler)法 向后欧拉公式 由于未知数 yn+1 同时出现在等式的两边，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好（后面分析）。 隐式欧拉公式中的未知数 yn+1 可通过以下迭代法求解： §8.2 欧拉(Euler)法 向后欧拉公式 迭代法求隐式欧拉格式中yn+1的收敛性 见上图， 显然，这种近似也有一定误差， 如何估计这种误差y(xn+1) ? yn+1 ？ 方法同上，基于Taylor展开估计局部截断误差。 但是注意，隐式公式中右边含有f(xn+1 , yn +1 ) , 由于yn +1不准确，所以不能直接用 y (xn+1)代替f(xn+1 , yn +1 ) 设已知曲线上一点 Pn (xn , yn ),过该点作弦线，斜率为(xn+1 , yn +1 ) 点的方向场f(x,y)。若步长h充分小，可用弦线和垂线x=xn+1的交点近似曲线与垂线的交点。 几何意义 xn xn+1 Pn Pn+1 x y y(x) ? 隐式欧拉法的局部截断误差： ? 隐式欧拉法的局部截断误差： 即隐式欧拉公式具有 1 阶精度。 ? 隐式欧拉法的局部截断误差： §8.2 欧拉(Euler)法 向后欧拉公式 比较欧拉显式公式和隐式公式及其局部截断误差 显式公式 隐式公式 若将这两种方法进行算术平均，即可消除误差 的主要部分/*leading term*/ 而获得更高的精度,称为梯形法 ? 梯形公式 /* trapezoid formula */ — 显、隐式两种算法的平均 注：的确有局部截断误差 ， 即梯形公式具有2 阶精度，比欧拉方法