Существует ли модуль python для решения / интеграции системы стохастических дифференциальных уравнений?

У меня есть система стохастических дифференциальных уравнений, которые я хотел бы решить. Я надеялся, что этот вопрос уже решен. Я немного обеспокоен построением моего собственного решателя, потому что боюсь, что мой решатель будет слишком медленным, и могут возникнуть проблемы с численной стабильностью.

есть ли модуль python для таких проблем?

Если нет, существует ли стандартный подход для решения таких систем.

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)