我有一个我想解决的随机微分方程系统。我希望这个问题已经得到解决。我有点担心构建自己的求解器，因为我担心我的求解器会太慢，并且可能存在数值稳定性问题。
是否有针对此类问题的 python 模块？
如果没有，是否有解决此类系统的标准方法。
问问题
3304 次
3 回答
8
有一个：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)
于 2013-01-30T01:01:50.317 回答
1
已接受答案中的链接不再起作用。还有sdeint：
https://pypi.org/project/sdeint/#description
它是在接受答案几年后发布的，看起来处于半主动开发阶段。文档中的第二个示例有一个具有恒定系数的 SDE 系统。我不确定他们是否支持更复杂的 SDE 系统。
于 2021-07-05T19:23:14.060 回答
0
diffeqpy包将 Julia 的DifferentialEquations.jl引入了 Python。这可以做很多事情，包括随机微分方程。
于 2022-01-17T03:29:37.893 回答