问题来源:1stopt 7参数拟合求助 - 计算模拟 - 程序代码 - 小木虫论坛-学术科研互动平台 (muchong.com)
原来理解有误,正确的解法参考:初值未知且需要传递中间变量的微分方程参数拟合_wlfc的专栏-CSDN博客
但本文下面的解法可作为 隐函数式微分方程参数拟合 例子。
根据两个式子(一个隐函数式微分方程,一个代数方程)拟合7个未知参数fac1~fac7。式中lambv'=d(lambv)/d(lamb),lambv为中间变量。
lambv'=((1.0/3/fac7)*(fac1*((lamb/lambv')^fac2-(lambv'/lamb)^(0.5*fac2))+fac3*((lamb/lambv')^fac4-(lambv'/lamb)^(0.5*fac4))+fac5*((lamb/lambv')^fac6-(lambv'/lamb)^(0.5*fac6))));
p=(pa+fac1*((lamb/lambv)^fac2/lamb-(lambv/lamb)^(fac2*0.5)/lamb)+fac3*((lamb/lambv)^fac4/lamb-(lambv/lamb)^(fac4*0.5)/lamb)+fac5*((lamb/lambv)^fac6/lamb-(lambv/lamb)^(fac6*0.5)/lamb));
实验数据:
t |
lamb |
pa |
p |
0 |
0.959926667 |
-0.218715994 |
-0.060966222 |
0.25 |
0.963063333 |
-0.200840559 |
-0.046836278 |
0.5 |
0.96612 |
-0.183553472 |