1stopt拟合常微分方程组
求助,考虑用1stopt软件拟合微分方程组参数,点了运行没有跑,不知道问题出在哪里,求大神们帮助一下,谢谢!
Title "pd1";
parameters k1,k2,k3,k4,k5;
Variable t,c1,c2,c3,c4,c5;
InitialODEValue c1=3.5,c2=3.5,c3=0,c4=0,c5=0;
ODEFunction
c1'=(-k1*c1*c2)*t;
c2'=(-k1*c1*c2-k2*c3*c2-k3*c4*c2)*t;
c3'=(k1*c2*c3-k2*c3*c2-k4*c3*c3-k5c4*c3)*t;
c4'=(k2*c3*c2+k4*c3*c3-k3*c4*c2-k5*c4*c3)*t;
c5'=(k3*c4*c2+k5*c4*c3)*t;
Data;
t c1 c2 c3 c4 c5
0 3.500 3.500 0.000 0.000 0.000
15 2.079 1.817 1.225 0.239 0.107
30 1.575 0.756 1.313 0.758 0.301
60 1.376 0.242 1.241 1.075 0.414
90 1.310 0.000 1.225 1.097 0.568
150 1.354 0.000 1.146 1.033 0.678
再就是,如果用matlab,应该如何编写,谢谢! 返回小木虫查看更多
如果不是正版的,是不能运行5个参数的
我用5.0计算了四参数,计算结果也不对。
Title "pdy123";
parameters k2,k3,k4,k5;
Variable t,c1,c2,c3,c4,c5;
InitialODEValue c1=3.5,c2=3.5,c3=0,c4=0,c5=0;
ODEFunction
c1'=(-0.005*c1*c2)*t;
c2'=(-0.005*c1*c2-k2*c3*c2-k3*c4*c2)*t;
c3'=(0.005*c2*c3-k2*c3*c2-k4*c3*c3-k5*c4*c3)*t;
c4'=(k2*c3*c2+k4*c3*c3-k3*c4*c2-k5*c4*c3)*t;
c5'=(k3*c4*c2+k5*c4*c3)*t;
Data;
t c1 c2 c3 c4 c5
0 3.500 3.500 0.000 0.000 0.000
15 2.079 1.817 1.225 0.239 0.107
30 1.575 0.756 1.313 0.758 0.301
60 1.376 0.242 1.241 1.075 0.414
90 1.310 0.000 1.225 1.097 0.568
150 1.354 0.000 1.146 1.033 0.678
这是4个参数的,但模拟出计算结果根本对不上
,
还有个问题,这个常微分方程,左边是dc/ct,我是把t直接放到右边,但感觉不对,正确应该如何处理,谢谢啦!
InitialODEValue t=0,c1=3.5,c2=3.5,c3=0,c4=0,c5=0;
Data;
//t c1 c2 c3 c4 c5
dc/(c*t) ? 还是 dc/dt
微分方程 中 将 dc/dt 写成 c'
直接写成c'就行?
左边是dc/dt
万分感谢
直接写成c',但我方程右边没有变量t了。实际是浓度随时间t的变化。这个不知道该怎么处理。