摘要: 非刚性常微分方程的数值解法通常会用四阶龙格库塔算法,其matlab函数对应ode45。 对于dy/dx = f(x,y),y(0)=y0。 其四阶龙格库塔公式如下: 对于通常计算,四阶已经够用,四阶以上函数f(x,y)计算工作量大大增加而精度提高较慢。 下面以龙格库塔法解洛伦兹方程为例: matla 阅读全文
posted @ 2020-01-28 14:35 Dsp Tian 阅读(21161) 评论(0) 推荐(0) 编辑