我有一个随机微分方程系统，我想解决。我希望这个问题已经解决了。我有点担心构造我自己的求解器，因为我担心我的解算器会太慢，并且可能存在数值稳定性的问题。
是否存在针对此类问题的python模块？
如果没有，是否有解决此类系统的标准方法。
答案 0 :(得分：7)
有一个：http://diffusion.cgu.edu.tw/ftp/sde/
网站示例：
""" add required Python packages """
from pysde import *
from sympy import *
""" Variables acclaimed """
x,dx=symbols('x dx')
r,G,e,d=symbols('r G epsilon delta')
""" Solve Kolmogorov Forward Equation """
l=sde.KolmogorovFE_Spdf(r*(G-x),e*x*(1-x),0,1)
sol=l.subs({e:r*d})
pprint(sol)
答案 1 :(得分：0)
已接受答案中的链接不再起作用。还有sdeint：
https://pypi.org/project/sdeint/#description
在接受答案几年后发布，看起来处于半主动开发中。文档中的第二个示例有一个具有常数系数的 SDE 系统。我不确定他们是否支持更复杂的 SDE 系统。