Simulating Stochastic Differential Equations with Conserved Quantities by Improved Explicit Stochastic Runge–Kutta Methods
Abstract
:1. Introduction
2. Preliminaries
- Denoteas the trees with only one vertex of color l,,
- Denoteas the tree formed by joining subtreeseach by a single branch to a common root.
- ,
- ,
- If, then.
3. Improved ESRK Methods
Algorithm 1 Improved ESRK method |
Step 1. Set the initial value and choose a constant . |
Step 2. |
if then |
end the procedure; |
else |
compute by the ESRK method; |
end if |
Step 3. |
if then |
let and go to Step 2; |
else |
go to Step 4; |
end if |
Step 4. Compute such that reaches the minimum. Compute
|
4. Numerical Examples
4.1. Stochastic Kubo Oscillator
- In order to show the proposed methods have the same mean-square convergence order as the original method. The mean-square convergence rate of methods (18) and (30) are displayed in Figure 1. Here, the mean-square errors are computed over 1000 simple paths with five different step sizes: , , , , and .
- In order to illustrate the advantages of the proposed methods in computational efficiency, we also apply the stochastic projection method and the stochastic discrete gradient method to solve (38) with . Figure 2 displays the mean-square sample errors and the corresponding CPU time of the stochastic projection method, the stochastic discrete gradient method, the Platen method, method (18), and method (30).
- In order to show the proposed methods’ ability to preserve the conserved quantity and get more accurate numerical solutions than the original method. We display the numerical sample paths and absolute errors in conserved quantity I of the Platen method, method (18), and method (30) in Figure 3 and Figure 4, respectively. Here, .
4.2. Stochastic Mathematical Pendulum
4.3. Stochastic Rigid Body
4.4. Stochastic Kepler Problem
5. Conclusions and Discussion
Author Contributions
Funding
Acknowledgments
Conflicts of Interest
References
- Mao, X. Stochastic Differential Equations and Applications; Elsevier Science: Amsterdam, The Netherlands, 2007. [Google Scholar]
- Cao, W.; Zhang, Z.; Karniadakis, G.E. Numerical methods for stochastic delay differential equations via the Wong–Zakai approximation. SIAM J. Sci. Comput. 2015, 37, A295–A318. [Google Scholar] [CrossRef]
- Hu, L.; Gan, S. Numerical analysis of the balanced implicit methods for stochastic pantograph equations with jumps. Appl. Math. Comput. 2013, 223, 281–297. [Google Scholar] [CrossRef]
- Li, X.; Cao, W. On mean-square stability of two-step Maruyama methods for nonlinear neutral stochastic delay differential equations. Appl. Math. Comput. 2015, 261, 373–381. [Google Scholar] [CrossRef]
- Wang, X.; Wu, J.; Dong, B. Mean-square convergence rates of stochastic theta methods for SDEs under a coupled monotonicity condition. BIT Numer. Math. 2020, 60, 1–32. [Google Scholar] [CrossRef]
- Wu, Q.; Hu, L.; Zhang, Z. Convergence and stability of balanced methods for stochastic delay integro-differential equations. Appl. Math. Comput. 2014, 237, 446–460. [Google Scholar] [CrossRef]
- Li, Q.; Gan, S. Mean-square exponential stability of stochastic theta methods for nonlinear stochastic delay integro-differential equations. J. Appl. Math. Comput. 2012, 39, 69–87. [Google Scholar] [CrossRef]
- Wang, X.; Gan, S. B-convergence of split-step one-leg theta methods for stochastic differential equations. J. Appl. Math. Comput. 2012, 38, 489–503. [Google Scholar] [CrossRef]
- Li, M.; Huang, C. Projected Euler-Maruyama method for stochastic delay differential equations under a global monotonicity condition. Appl. Math. Comput. 2020, 366, 124733. [Google Scholar] [CrossRef]
- Milstein, G.N. Numerical Integration of Stochastic Differential Equations; Kluwer Academic Publishers: Amsterdam, Poland, 1995. [Google Scholar]
- Kloeden, P.E.; Platen, E. Numerical Solution of Stochastic Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 23. [Google Scholar]
- Burrage, K.; Burrage, P.M. High strong order explicit Runge-Kutta methods for stochastic ordinary differential equations. Appl. Numer. Math. 1996, 22, 81–101. [Google Scholar] [CrossRef]
- Zeghdane, R.; Abbaoui, L.; Tocino, A. Higher-order semi-implicit Taylor schemes for Itô stochastic differential equations. J. Comput. Appl. Math. 2011, 236, 1009–1023. [Google Scholar] [CrossRef] [Green Version]
- Burrage, K.; Tian, T. Predictor-corrector methods of Runge–Kutta type for stochastic differential equations. SIAM J. Numer. Anal. 2002, 40, 1516–1537. [Google Scholar] [CrossRef] [Green Version]
- Milstein, G.N.; Platen, E.; Schurz, H. Balanced implicit methods for stiff stochastic systems. SIAM J. Numer. Anal. 1998, 35, 1010–1019. [Google Scholar] [CrossRef] [Green Version]
- Omar, M.; Aboul-Hassan, A.; Rabia, S.I. The composite Milstein methods for the numerical solution of Stratonovich stochastic differential equations. Appl. Math. Comput. 2009, 215, 727–745. [Google Scholar] [CrossRef]
- Burrage, K.; Burrage, P.M. Order conditions of stochastic Runge–Kutta methods by B-series. SIAM J. Numer. Anal. 2000, 38, 1626–1646. [Google Scholar] [CrossRef] [Green Version]
- Rößler, A. Runge–Kutta methods for the strong approximation of solutions of stochastic differential equations. SIAM J. Numer. Anal. 2010, 48, 922–952. [Google Scholar] [CrossRef]
- Li, X.; Zhang, C.; Ma, Q.; Ding, X. Discrete gradient methods and linear projection methods for preserving a conserved quantity of stochastic differential equations. Int. J. Comput. Math. 2018, 95, 2511–2524. [Google Scholar] [CrossRef]
- Hong, J.; Zhai, S.; Zhang, J. Discrete gradient approach to stochastic differential equations with a conserved quantity. SIAM J. Numer. Anal. 2011, 49, 2017–2038. [Google Scholar] [CrossRef]
- Cohen, D.; Dujardin, G. Energy-preserving integrators for stochastic Poisson systems. Commun. Math. Sci. 2014, 12, 1523–1539. [Google Scholar] [CrossRef]
- Malham, S.J.A.; Wiese, A. Stochastic Lie group integrators. SIAM J. Sci. Comput. 2008, 30, 597–617. [Google Scholar] [CrossRef] [Green Version]
- Wang, Z.; Ma, Q.; Yao, Z.; Ding, X. The Magnus expansion for stochastic differential equations. J. Nonlinear Sci. 2020, 30, 419–447. [Google Scholar] [CrossRef]
- Chen, C.; Cohen, D.; D’Ambrosio, R.; Lang, A. Drift-preserving numerical integrators for stochastic Hamiltonian systems. Adv. Comput. Math. 2020, 46, 1–22. [Google Scholar] [CrossRef] [Green Version]
- D’Ambrosio, R.; Scalone, C. On the numerical structure preservation of nonlinear damped stochastic oscillators. Numer. Algorithms 2020. [Google Scholar] [CrossRef]
- Misawa, T. A Lie algebraic approach to numerical integration of stochastic differential equations. SIAM J. Sci. Comput. 2001, 23, 866–890. [Google Scholar] [CrossRef]
- Misawa, T. Symplectic integrators to stochastic Hamiltonian dynamical systems derived from composition methods. Math. Probl. Eng. 2010, 2010, 384937. [Google Scholar] [CrossRef]
- Hairer, E.; Lubich, C.; Wanner, G. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006; Volume 31. [Google Scholar]
- Brent, R.P. Algorithms for Minimization without Derivatives; Courier Corporation: North Chelmsford, MA, USA, 2013. [Google Scholar]
- Milstein, G.N.; Repin, Y.M.; Tretyakov, M.V. Numerical methods for stochastic systems preserving symplectic structure. SIAM J. Numer. Anal. 2002, 40, 1583–1604. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Wang, Z.; Ma, Q.; Ding, X. Simulating Stochastic Differential Equations with Conserved Quantities by Improved Explicit Stochastic Runge–Kutta Methods. Mathematics 2020, 8, 2195. https://doi.org/10.3390/math8122195
Wang Z, Ma Q, Ding X. Simulating Stochastic Differential Equations with Conserved Quantities by Improved Explicit Stochastic Runge–Kutta Methods. Mathematics. 2020; 8(12):2195. https://doi.org/10.3390/math8122195
Chicago/Turabian StyleWang, Zhenyu, Qiang Ma, and Xiaohua Ding. 2020. "Simulating Stochastic Differential Equations with Conserved Quantities by Improved Explicit Stochastic Runge–Kutta Methods" Mathematics 8, no. 12: 2195. https://doi.org/10.3390/math8122195
APA StyleWang, Z., Ma, Q., & Ding, X. (2020). Simulating Stochastic Differential Equations with Conserved Quantities by Improved Explicit Stochastic Runge–Kutta Methods. Mathematics, 8(12), 2195. https://doi.org/10.3390/math8122195