第6章 常微分方程数值解法.pptVIP

四阶Runge-Kutta方法第27页，共76页，星期日，2025年，2月5日6.2.3R-K法的稳定性第28页，共76页，星期日，2025年，2月5日R-K法的稳定性第29页，共76页，星期日，2025年，2月5日R-K法的稳定性第30页，共76页，星期日，2025年，2月5日6.2.5隐式R-K法第31页，共76页，星期日，2025年，2月5日隐式R-K法第32页，共76页，星期日，2025年，2月5日隐式R-K法第33页，共76页，星期日，2025年，2月5日隐式R-K法第34页，共76页，星期日，2025年，2月5日隐式R-K法第35页，共76页，星期日，2025年，2月5日6.3线形多步法单步法主要依据yn的信息去计算yn+1。线性多步法是想依据yn，yn-1…，yn-r(r≥1)的信息去计算yn+1。考虑到线性组合较为方便，因此，线性多步法一般形式可设为第36页，共76页，星期日，2025年，2月5日6.3.1基于数值积分的方法第37页，共76页，星期日，2025年，2月5日基于数值积分的方法第38页，共76页，星期日，2025年，2月5日第39页，共76页，星期日，2025年，2月5日基于数值积分的方法第40页，共76页，星期日，2025年，2月5日基于数值积分的方法第41页，共76页，星期日，2025年，2月5日基于数值积分的方法第42页，共76页，星期日，2025年，2月5日基于数值积分的方法Adams预估—校正法预估校正并取第43页，共76页，星期日，2025年，2月5日6.3.2基于Taylar展开式的方法第44页，共76页，星期日，2025年，2月5日第45页，共76页，星期日，2025年，2月5日基于Taylar展开式的方法第46页，共76页，星期日，2025年，2月5日基于Taylar展开式的方法第47页，共76页，星期日，2025年，2月5日第48页，共76页，星期日，2025年，2月5日第49页，共76页，星期日，2025年，2月5日第50页，共76页，星期日，2025年，2月5日第51页，共76页，星期日，2025年，2月5日6.4一阶常微分方程组数值解法在许多实际问题中，常常出现高阶微分方程和高阶微分方程组，通过引入新的变量，总可化为一阶微分方程组。由此可知，讨论一阶常微分方程组的数值解法是很有意义的。第52页，共76页，星期日，2025年，2月5日6.4.1解一阶常微分方程组的R-K方法第53页，共76页，星期日，2025年，2月5日一阶常微分方程组的R-K方法第54页，共76页，星期日，2025年，2月5日第1页，共76页，星期日，2025年，2月5日绪论在工程和科学计算中，所建立的各种常微分方程的初值或边值问题，除很少几类的特殊方程能给出解析解，绝大多数的方程是很难甚至不可能给出解析解的，其主要原因在于积分工具的局限性。因此，人们转向用数值方法去解常微分方程，并获得相当大的成功，讨论和研究常微分方程的数值解法是有重要意义的。第2页，共76页，星期日，2025年，2月5日6.1初值问题的Euler方法第3页，共76页，星期日，2025年，2月5日初值问题的Euler方法第4页，共76页，星期日，2025年，2月5日初值问题的Euler方法第5页，共76页，星期日，2025年，2月5日初值问题的Euler方法第6页，共76页，星期日，2025年，2月5日初值问题的Euler方法第7页，共76页，星期日，2025年，2月5日初值问题的Euler方法第8页，共76页，星期日，2025年，2月5日初值问题的Euler方法第9页，共76页，星期日，2025年，2月5日初值问题的Euler方法第10页，共76页，星期日，2025年，2月5日初值问题的Euler方法第11页，共76页，星期日，2025年，2月5日初值问题的Euler方法第12页，共76页，星期日，2025年，2月5日xEuler法y改进的Euler法y精确解01.0000001.0000001.0000000.11.0000001.0959091.0954450.21.1918181.1840971.1832160.31.2774381.2662011.2649110.41.3582131.3433601.3416410.51.4351331.4164021.4142140.61.5089661.4