1.龙格-库塔法(单步法、精度高)
基于改进的欧拉法,改进欧拉法如下,使用两端的导数的算术平均值来近似两点间的差商(增量函数),从而计算下一点,其中端点的y通过前向欧拉法计算得到
增量函数可以进一步改写为两点间多个点的导数的加权和,从而进一步提升改进欧拉法的精度,这就是龙格库塔法。下面是一种推广方式
要使这种方法保持2阶精度,需要满足以下条件,可以通过局部截断误差计算
龙格库塔法的一般形式(k个点导数的加权平均):
其中,下面的符号为待定系数
可以通过下列表格的形式描述参数(注意这里a是上面c和a的积)
这里在求解待定系数或者验证阶次时需要注意,二元泰勒展开公式为:
此外,还需注意全导数与偏导数之间的区别:
龙格库塔法的导出基于泰勒展开,因此精度受函数的光滑性影响,如果光滑性不好,应使用低阶算法,h取小。
截图源自【《数值分析》| 华科 | 研究生基础课】https://www.bilibili.com/video/BV1AK4y1k7Px?p=34&vd_source=a53b34e44cbfd40d72a5b337c3e5a13d