第5章---常微分方程初值问题.pptVIP

* 2.局部截断误差与整体截断误差之间的关系 若单步法(17)的局部截断误差是p+1阶的，即 c1独立于h(不依赖于h或与h无关)， 且函数Φ(x, u, v, h)在区域 定理4 上关于u,v 满足条件 则单步法(17)是 p 阶方法。 (20) * 推论: 当f 在D上满足基本条件时,单步法(17)的阶由局部截断误差的阶来确定。 结论： 一般情况下，若局部截断误差是 p+1阶的，则单步法是p 阶方法。 说明： 可用Taylor展开法估计 的阶，即 的主项 高阶项 则对应的单步法是p 阶方法。 由于 (21) * 3. 定理的应用（讨论各种方法的阶数） 显式Euler方法是一阶方法 事实上,当(1),(2)的解y(x)二阶连续可导时, f (x, y)关于y满足Lipschitz条件，Lipschitz常数为L，其局部截断误差为: 将y(xk+1)在xk点展开可得 局部截断误差有界，存在M，使得 记 于是 * 反复递推得： * 显式Euler方法的整体截断误差与h同阶。 * 隐式Euler方法是一阶方法 事实上，局部截断误差为: 将 f (x, y)在y用微分中值定理，有 则 将方程的解作Taylor展开 * 因此 * 梯形方法是二阶方法 将 f (xk+1, yk+1)在xkTaylor展开，有 则 将方程的解作Taylor展开 将 f (x, y)在y用微分中值定理，有 * 因此 * 预估－校正Euler方法是二阶方法 事实上，局部截断误差为： * 二阶方法 即得整体截断误差: 结论： 一般情况下，若局部截断误差是 p+1阶的，则单步法是p 阶方法。 方 法 ? ? 显式欧拉 隐式欧拉 梯形公式 简单 精度低 稳定性最好 精度低, 计算量大 精度提高 计算量大 设y(x)是初值问题的解 ,yn是单步法产生的近似解.如果对任意固定的点xn,均有lim yn =y(xn)则称单步法是收敛的. 可见,若方法是收敛的,则当h?0时,整体截断误差en=y(xn)-yn将趋于零. 定理 设单步方法是p?1阶方法, 增量函数?(x,y,h)在区域?a?x?b,-?y+?,0?h?h0?上连续,且关于y满足Lipschitz条件,初始近似y0=y(a)=?,则方法是收敛的,且存在与h无关的常数C,使 |y(xn)-yn|?Chp 五、 单步法解的收敛性和稳定性 * 唯一性 若f (x, y)在D上满足基本条件，f (x, y)关于y的Lipschitz常数为L时，只要hkL1，则由定理1可知(13)确定了唯一的yk+1。同理，只要hkL2，(14)确定了唯一的yk+1. 收敛性 (以(14)式为例说明) 当 时,以y 为变量的函数: 同样，从任意初始值出发,迭代(13)、 (14)，都收敛到 事实上， 单步隐式法解的唯一性、收敛性 * 在 上关于y 满足Lipschitz条件,且L－常数为 , 由压缩不动点定理得方程 有唯一不动点 , 而且从 出发,迭代 (15) 都收敛到 见非线性方程数值解的定理3 * 单步法稳定性 例：考察初值问题 在区间[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 * 定义 若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的。 讨论数值方法的稳定性, 通常仅限于典型的试验方程 其中?是复数。在复平面上,当方法稳定时要求变量?h的取值范围称为方法的绝对稳定域, 它与实轴的交集称为绝对稳定区间. * 将Euler方法应用于方程y? = ?y, 得到 设在计算yn时产生误差?n,计算值?yn = yn + ?n,