matlab实现龙格古塔算法
1、文档下载:
本算法已经整理成文档如下,有需要的朋友可以点击进行下载
序号 | 文档(点击下载) |
---|---|
本项目文档 | 【老生谈算法】matlab实现龙格古塔算法.docx |
2、算法详解:
龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高,采取措施对误差进行抑制,所以其实现原理也较复杂。该算法是构建在数学支持的基础之上的。龙格库塔方法的理论基础来源于泰勒公式和使用斜率近似表达微分,它在积分区间多预计算出几个点的斜率,然后进行加权平均,用做下一点的依据,从而构造出了精度更高的数值积分计算方法。如果预先求两个点的斜率就是二阶龙格库塔法,如果预先取四个点就是四阶龙格库塔法。一阶常微分方程可以写作:y’=f(x,y),使用差分概念。
(Yn+1-Yn)/h= f(Xn,Yn)推出(近似等于,极限为Yn’)
Yn+1=Yn+hf(Xn,Yn)
另外根据微分中值定理,存在0<t<1,使得