分析09—常微方程数值解法.pptVIP

第八章;第八章目录;第八章 序;第八章 序;微分方程 数值解;§1 欧拉（Euler）法 ;1.1 Euler法及其简单改进 ;Euler公式的推导(续1）;Euler公式的推导（续2）;Euler公式的推导（续3）;Euler公式的推导（续4）;Euler公式的推导（续5）;;Euler公式注释; 显然，步长h越小，阶数P越高，局部截断误差越小，当 然计算精度越高； ;关于Euler法的整体截断误差注释; 式（8-7）表明，Euler方法的整体截断误差与h同阶，当h?0时，en?0。 ;结 论 ;3）计算结果与解法的阶数p，真解的导数y (p+1)有关，p越大，h p+1越小， |y(p+1)(ξ)|的上限越大，M 也越大，因此为保证精度当然应选阶数p较高的方法。但如果M 很大，当f (x,y)是分段连续的函数时，则应采用低阶的方法如用Euler法。;1.2 改进的Euler法 ;改进的Euler法(续）;预测──校正型公式 ;改进的Euler公式的局部截断误差分析;改进的Euler公式举例;表格8-1;表格8-2;??8-2计算结果说明（续）;例2 ;§2 龙格-库塔（Runge-kutta）方法 ;龙格-库塔方法的基本思想(续）;2.2 二阶龙格-库塔公式 ;二阶龙格-库塔公式（续1）;二阶龙格-库塔公式（续2）; 特别，当l=1即xn+l=xn+1时， c1=c2=1/2，二阶R-K公式就 是改进欧拉公式。 如果取l=1/2，则c1=0,c2=1， 这时二阶R-K公式称为变形的欧拉公式，其形式见左边： ;构造二阶R-K公式的主要步骤 ;2.3 高阶R-K公式 ;高阶R-K公式（续1）;高阶R-K公式（续2）; 若需再将精度提高至四阶，用上述类似处理方法，只 是必须在[xi,xi+1]上用四个点处的斜率加权平均作为k*的近 似值，构成一族四阶龙格-库塔公式，由于推导复杂，这 里从略，只将常用的公式介绍如下： ;Gill 公式;初值问题两种方法举例;对四阶标准R-K法， 取h=0.2，计算公 式如下：;例 3 （续2）;表8-6 计算f 的次数与精度阶数的关系 ;表8-6 的说明（续）;2.4 变步长R-K法 ;变步长R-K法（续1）;具体做法;变步长四阶标准R-K法解初值问题举例 ;步长折半次数;§3 线性多步法 ;3.1 阿当姆斯Adams公式 ;阿当姆斯Adams公式（续）;r+1步阿当姆斯显式公式 （又称Adams-Bashforth公式） ;r+1步阿当姆斯隐式公式 （又称Adams- Moulton公式） ; 其系数见下表8-9，显然也具有r+1阶精度。 ;四阶阿当姆斯显式公式;四阶阿当姆隐式公式内插公式及其局部截断误差为（取xn-2,xn-1,xn,xn+1为插值节点）： ;阿当姆斯公式举例;阿当姆斯预测—校正公式;于是得pn+1与cn+1的事后误差估计： ;阿当姆斯带误差修正的预测—校正公式;3.2 哈明(Hamming)公式 ;哈明(Hamming)公式（续1）;哈明(Hamming)公式（续2）;Hamming公式（预测一校正公式） ; 前面介绍的四套预测一校正公式（8-28）、（8-29）、（8-32）、（8-33）都是实际计算中常用的方法，它们的 共同特点是计算量省而精度高，但都必须用别的单步法 为其提供出发值，由于解题过程是递推进行的，所以出 发值的精度对以后的计算影响很大，故一般采用四阶龙 格一库塔法求其出发值。 ;§4 一阶方程组与高阶方程初值问题 ;一阶方程组（续1）;一阶方程组（续2）;4.2 高阶方程初值问题 ;高阶方程初值问题举例;例5（续）;§5 收敛性与稳定性 ;5.1 收敛性 ;5.2 稳定性 ;稳定性（续1）;表8-12列出了几个常用公式的绝对稳定区间。 从表中可以看到，隐式公式比同阶显式公式的稳定区间大，这就是为什么最好立足于隐式公式求解的原因。 ;稳定性（续3）;综上所述，当用数值方法解微分方程时，步长h的选择由两个因素确定。一是精度要求，由收敛性可知，h越小精度越高。 二是稳定性要求， 应落在所选公式的 绝对稳定域中。虽然h十分小时两个条件都满足，但h太小计算步数势必大增，每步的微小误差累积起来也很可观。因此，实际计算中在满足精度要求与稳定性的条件下，步长h应尽可能取得大些。 ;第八章