Question

MATLAB四阶五级 Runge-Kutta-Felhberg 算法计算一阶微分方程，方程系数会随着上一步计算结果变化怎么办？

matlab 四阶五级 Runge-Kutta-Felhberg 算法，计算一阶微分方程，其中，微分方程中的一些系数会随着上一步计算结果变化怎么办？还可以用 ode45（） 吗？cp=1.676480*10^(-2)*(3.167349*1000+p)^(0.4666343)*(47.33943+T)^(0.2518024)*(T+2.001211*p)^(0.01549791);miu=1.4*10^(-15)*(6.426347*100+p)^(6.466196)*(-8.049315+T)^(2.867102)*(T+(-1.51141/100)*p)^(-4.537339); lou=7.007429*10^(6)*(7.367185*10+p)^(0.1159146)*(4208.254+T)^(-1.962483)*(T+1.308541*p)^(0.02299075);





Answer 1

题主给出的这种指茄类型微分方程，是可以用四阶五级 Runge-Kutta-Felhberg 算法和ode45（）函汪首数。
ode45（）函数实际上就是变步长四阶五级 Runge-Kutta（龙格困逗数-库塔）法。为了更好地回答，请题主给出具体的微分方程，以便说明问题。

追问

密度、粘度等那3个公式matlab代码：
miu=1.4*10^(-15)*(6.426347*100+p)^(6.466196)*(-8.049315+T)^(2.867102)*(T+(-1.51141/100)*p)^(-4.537339);

追答

根据实际给出的微分方程，属于刚性方程，不适合用四阶五级 Runge-Kutta-Felhberg 算法（ode45），而应该用ode23s函数来求解。