jpdyc
极品豆芽菜
2020-11-29 20:52
采纳率: 100%
浏览 82

使用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函数推荐?最好能给出一个大致的框架。谢谢!

  • 点赞
  • 写回答
  • 关注问题
  • 收藏
  • 邀请回答

2条回答 默认 最新

  • weixin_43162240
    Bryce__LY 2020-11-30 10:35
    已采纳
    from scipy.integrate import odeint  这个可以  但是方程太多了不知道行不行  之前解过一个三个等式的可以用 
    点赞 评论
  • jpdyc
    极品豆芽菜 2020-11-30 16:10

    关键是这里面既有微分方程又有普通方程影响微分方程里的参数。。。不知道你之前有没有解过类似的方程?

    点赞 评论

相关推荐