鸠魁 发表于 2020-6-8 15:48:23

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,应该如何编写,谢谢!
页: [1]
查看完整版本: 1stopt拟合常微分方程组