是否有一个python模块来解决/整合随机微分方程组？

Question

我有一个随机微分方程系统,我想解决.我希望这个问题已经解决了.我有点担心构造我自己的求解器,因为我担心我的解算器会太慢,并且可能存在数值稳定性的问题.

是否存在针对此类问题的python模块？

如果没有,是否有解决此类系统的标准方法.

Answer 1

有一个: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)