使用python建立物理模型并解常微分方程,普通方程组

最近学校项目需要使用python重现并扩展一个物理模型,其主要方程包含5个常微分方程和7个普通方程组。
方程组如下:

其中S8、S4、S2、S、Sp、EH、EL、iH、iL、nH、nL、V为变量,共12个,且互相影响。

其他详细参数可以参考原文:https://doi.org/10.1039/C5CP05755H

之前尝试使用了scipy.intergrate.odeint函数求解,但是报错...

在文章中作者使用MATLAB对该方程组求解并且使用了jacobian矩阵(和dae函数?)结果如下:

所以想请教各位大佬有没有什么思路或者可用的python函数推荐?最好能给出一个大致的框架。谢谢!

查看全部
jpdyc
极品豆芽菜
2020/11/29 20:52
  • python
  • 点赞
  • 收藏
  • 回答
    私信
满意答案
查看全部

1个回复