Numerical Methods for Stochastic Systems Preserving Symplectic Structure

Abstract

Stochastic Hamiltonian systems with multiplicative noise, phase flows of which preserve symplectic structure, are considered. To construct symplectic methods for such systems, sufficiently general fully implicit schemes, i.e., schemes with implicitness both in deterministic and stochastic terms, are needed. A new class of fully implicit methods for stochastic systems is proposed. Increments of Wiener processes in these fully implicit schemes are substituted by some truncated random variables. A number of symplectic integrators is constructed. Special attention is paid to systems with separable Hamiltonians. Some results of numerical experiments are presented. They demonstrate superiority of the proposed symplectic methods over very long times in comparison with nonsymplectic ones.