Numerical Picard iteration methods for simulation of non-Lipschitz stochastic differential equations
- In this paper, we present splitting approaches for stochastic/deterministic coupled differential equations, which play an important role in many applications for modelling stochastic phenomena, e.g., finance, dynamics in physical applications, population dynamics, biology and mechanics. We are motivated to deal with non-Lipschitz stochastic differential equations, which have functions of growth at infinity and satisfy the one-sided Lipschitz condition. Such problems studied for example in stochastic lubrication equations, while we deal with rational or polynomial functions. Numerically, we propose an approximation, which is based on Picard iterations and applies the Doléans-Dade exponential formula. Such a method allows us to approximate the non-Lipschitzian SDEs with iterative exponential methods. Further, we could apply symmetries with respect to decomposition of the related matrix-operators to reduce the computational time. We discuss the different operator splitting approaches for a nonlinear SDE with multiplicative noise and compare this to standard numerical methods.
Creative Commons - CC BY 4.0 - Attribution 4.0 International