1stopt拟合常微分方程组

来源:matlab  › 求助  › 倒序浏览
  • 回答数

    0

  • 浏览数

    2440

  • 收藏数

    0

作者:鸠魁 发表于 2020-6-8 15:48:23
跳转到指定楼层
求助,考虑用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,应该如何编写,谢谢!
分享:
回复

使用道具

成为第一个回答人

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录