龙格-库塔(Runge-Kutta)方法
龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高,采取措施对误差进行抑制,所以其实现原理也较复杂。该算法是构建在数学支持的基础之上的。
对于一阶精度的欧拉公式有:
yi+1=yi+hki
其中 h 为步长,则
当用点 xi 处的斜率近似值 k1 与右端点 xi+1 处的斜率 k2 的算术平均值作为平均斜率 k∗ 的近似值,那么就会得到二阶精度的改进欧拉公式:
yi+1=yi+h(k1+k2)
其中 k1=f(xi,yi) ,