Statistics for stochastic differential equations and approximations of resolvent
arXiv:2207.13831 · doi:10.7566/JPSJ.92.074006
Abstract
The numerical evaluation of statistics plays a crucial role in statistical physics and its applied fields. It is possible to evaluate the statistics for a stochastic differential equation with Gaussian white noise via the corresponding backward Kolmogorov equation. The important notice is that there is no need to obtain the solution of the backward Kolmogorov equation on the whole domain; it is enough to evaluate a value of the solution at a certain point that corresponds to the initial coordinate for the stochastic differential equation. For this aim, an algorithm based on combinatorics has recently been developed. In this paper, we discuss a higher-order approximation of resolvent, and an algorithm based on a second-order approximation is proposed. The proposed algorithm shows a second-order convergence. Furthermore, the convergence property of the naive algorithms naturally leads to extrapolation methods; they work well to calculate a more accurate value with fewer computational costs. The proposed method is demonstrated with the Ornstein-Uhlenbeck process and the noisy van der Pol system.
13 pages, 5 figures
References in corpus (6)
- Duality and hidden symmetries in interacting particle systems
- Algorithms for Solving High Dimensional PDEs: From Nonlinear Monte Carlo to Machine Learning
- A mass, momentum, and energy conservative dynamical low-rank scheme for the Vlasov equation
- Dynamic tensor approximation of high-dimensional nonlinear PDEs
- Numerical methods to evaluate Koopman matrix from system equations
- Connection among stochastic Hamilton-Jacobi-Bellman equation, path-integral, and Koopman operator on nonlinear stochastic optimal control