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

由于 ，有 ∴ 证毕 第三十一页，共八十一页，2022年，8月28日 §9.3 龙格-库塔（Runge-Kutta）法 9.3.1 龙格-库塔(Runge-Kutta)法的基本思想 Euler公式可改写成 则yi+1的表达式y(xi+1)与的Taylor展开式的前两项完全相同,即局部截断误差为 。 改进的Euler公式又可改写成 第三十二页，共八十一页，2022年，8月28日 上述两组公式在形式上有一个共同点:都是用f(x,y)在某些点上值的线性组合得出y(xi+1)的近似值yi+1,而且增加计算的次数f(x,y)的次数,可提高截断误差的阶。如欧拉公式:每步计算一次f(x,y)的值,为一阶方法。改进欧拉公式需计算两次f(x,y)的值，它是二阶方法。它的局部截断误差为 。 第三十三页，共八十一页，2022年，8月28日 于是可考虑用函数f(x,y)在若干点上的函数值的线性组合来构造近似公式，构造时要求近似公式在(xi,yi)处的Taylor展开式与解y(x)在xi处的Taylor展开式的前面几项重合，从而使近似公式达到所需要的阶数。既避免求偏导,又提高了计算方法精度的阶数。或者说,在 这一步内多预报几个点的斜率值，然后将其加权平均作为平均斜率，则可构造出更高精度的计算格式，这就是龙格—库