常微分方程初值问题数值解法.ppt

{注意到第30页，共69页，星期日，2025年，2月5日对于初值问题，如果关于满足（向前Euler方法的整体截断误差）Lipschitz条件，为对应的Lipschitz常数，当时，向前Euler方法的数值解一致收敛于初值问题的精确解，且整体截断误差满足估计式如果，Euler方法的整体截断误差为第31页，共69页，星期日，2025年，2月5日一、Runge-Kutta方法的基本思想§3龙格-库塔（Runge-Kutta）方法显式单步法的一般形式：R-K方法是利用一些点的线性组合构造增量函数，使得相应方法的局部截断误差的阶数尽可能高。?二阶Runge-Kutta方法确定参数，使得与在点的Taylor展开式有尽可能多的相同项。第32页，共69页，星期日，2025年，2月5日比较两式的相同项得方程组有无穷多解第33页，共69页，星期日，2025年，2月5日若取其一组解则得到改进的Euler公式（二阶方法）若取其另一组解则得到二阶的Heun（休恩）公式。第34页，共69页，星期日，2025年，2月5日二、显式Runge-Kutta方法及其稳定性和设是一个正整数，代表使用函数值的个数，是一些特定的权因子（均为实数），则称下列方法（公式）为初值问题的m级显式Runge–Kutta公式，其中第35页，共69页，星期日，2025年，2月5日类似前面的处理方法，可以得到四级方法：m=4局部截断误差最常用的一种四阶方法：经典显式Runge-Kutta公式第36页，共69页，星期日，2025年，2月5日解：例2：用经典的四阶Runge-Kutta方法求解下列初值问题。经典的四阶Runge-Kutta公式：第37页，共69页，星期日，2025年，2月5日第38页，共69页，星期日，2025年，2月5日第39页，共69页，星期日，2025年，2月5日注：?对于显式N级R-K方法，最多只能得到N阶方法。?上述方法的缺陷：计算非常复杂。?可通过积分方法确定参数。例2：确定如下三级三阶显式Runge-Kutta公式中的参数：解：对微分方程两边积分得第40页，共69页，星期日，2025年，2月5日采用Simpson公式计算上式右端积分项可设参数则有选择剩余参数，使得第41页，共69页，星期日，2025年，2月5日取第42页，共69页，星期日，2025年，2月5日第43页，共69页，星期日，2025年，2月5日取利用Taylor展开式第44页，共69页，星期日，2025年，2月5日代入当时，第45页，共69页，星期日，2025年，2月5日例3：求经典四阶的R-K方法的绝对稳定域。解：第46页，共69页，星期日，2025年，2月5日其绝对稳定域为三、隐式Runge-Kutta方法m级隐式R–K方法的一般形式其中系数的确定方法同显式R–K方法完全类似第47页，共69页，星期日，2025年，2月5日(1)一级二阶的隐式中点方法：(2)二级四阶的隐式R-K方法：N级隐式R-K法可以达到2N阶缺陷：需要求解非线性方程（组）第48页，共69页，星期日，2025年，2月5日一、k步线性多步法§4线性多步法/*LinearMutistepMethodandPredictor-CorrectorFormat*/所谓的线性多步法，指的是某一步解的公式不仅与前一步的值有关，而且与前面若干步解的值有关的方法。对初值问题两边积分得第49页，共69页，星期日，2025年，2月5日数值分析第9章常微分方程初值问题数值解法*第1页，共69页，星期日，2025年，2月5日第九章常微分方程初值问题数值解法/*NumericalMethodforOrdinaryDifferentialEquations*/许多实际问题的数学模型是微分方程或微分方程的初值问题,如物体运动,电路震荡,化学反映及生物群体的变化等。能用解析方法求出精确解的微分方程为数不多,而且有的方程即使有解析解,也可能由于解的表达式非常复杂而不易计算,因此有必要研究