我有一个随机微分方程系统需要解决。希望已经有类似的现成解决方案。如果我自己构建求解器，担心速度太慢且数值稳定性存在问题。
是否有相应的Python模块可以解决这些问题？
如果没有，那么有没有标准方法来解决这些系统。
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)
- Kabie
1
谢谢！这很有帮助。 - CraigF
1
链接在被接受的答案中已经无法使用。还有sdeint：
这个答案发布几年后，看起来仍在半积极的发展中。文档中的第二个示例具有恒定系数的SDE系统。我不确定它们是否支持更复杂的SDE系统。
- jman
网页内容由stack overflow 提供, 点击上面的可以查看英文原文，
原文链接
原文链接