Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)



On quasi-Monte Carlo simulation
of Stochastic Differential Equations

Authors: Norbert Hofmann and Peter Mathé
Journal: Math. Comp. 66 (1997), 573-589
MSC (1991): Primary 65C05, 65C10; Secondary 60H10
MathSciNet review: 1397444
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: In a number of problems of mathematical physics and other fields stochastic differential equations are used to model certain phenomena. Often the solution of those problems can be obtained as a functional of the solution of some specific stochastic differential equation. Then we may use the idea of weak approximation to carry out numerical simulation. We analyze some complexity issues for a class of linear stochastic differential equations (Langevin type), which can be given by

\begin{displaymath}dX_{t}=-\alpha X_{t}dt+\beta (t)dW_{t}, \quad X_{0}:= 0, \end{displaymath}

where $\alpha >0$ and $\beta \! : [0,T]\to {\mathbb R}$. It turns out that for a class of input data which are not more than Lipschitz continuous the explicit Euler scheme gives rise to an optimal (by order) numerical method. Then we study numerical phenomena which occur when switching from (real) Monte Carlo simulation to quasi-Monte Carlo simulation, which is the case when we carry out the simulation on computers. It will easily be seen that completely uniformly distributed sequences yield good substitutes for random variates, while not all uniformly distributed (mod 1) sequences are suited. In fact we provide necessary conditions on a sequence in order to serve quasi-Monte Carlo purposes. This condition is expressed in terms of the measure of well-distributions. Numerical examples complement the theoretical analysis.

References [Enhancements On Off] (What's this?)

  • 1. N. N. Chentsov. Pseudo-random numbers for modeling Markov chains. Zh. Vychisl. Mat. i Mat. Fiz., 7:632 - 643, 1967.
  • 2. J. N. Franklin. Deterministic simulation of random processes. Math. of Computation, 17:28 - 59, 1963. MR 26:7125
  • 3. E. Hlawka. Lösung von Integralgleichungen mittels zahlentheoretischer Methoden. Österreich. Akad. Wiss. Math.-Natur. Kl. Sitzungsber. II, 171:103 - 123, 1962. MR 27:548
  • 4. N. Hofmann. Beiträge zur schwachen Approximation stochastischer Differentialgleichungen. Dissertation, HU Berlin, 1995.
  • 5. I. Karatzas and S. E. Shreve. Brownian Motion and Stochastic Calculus, volume 113 of Graduate Texts in Math., Springer, New York, Berlin, Heidelberg, 1988. MR 89c:60096
  • 6. P. E. Kloeden and E. Platen. Numerical Solution of Stochastic Differential Equations, volume 23 of Applications of Math., Springer, Berlin, Heidelberg, New York, 1992. MR 94b:60069
  • 7. D. E. Knuth. The Art of Computer Programming, Vol. 2/ Seminumerical Algorithms. Addison-Wesley Publ. Co., Reading, Mass., 1969. MR 44:3531
  • 8. L. Kuipers and H. Niederreiter. Uniform Distribution of Sequences. Wiley & Sons, New York, London, Sydney, Toronto, 1974. MR 54:7415
  • 9. P. Mathé. Approximation theory of Monte Carlo methods. Habilitation thesis, 1994.
  • 10. G. N. Milstein. Numerical Integration of Stochastic Differential Equations, volume 313 of Mathematics and its Appl., Kluwer Acad. Publ., Dordrecht, Boston, London, 1995. MR 96e:65003
  • 11. H. Niederreiter. Quasi-Monte Carlo methods and pseudo-random numbers. Bull. Amer. Math. Soc., 84(6):957 - 1041, 1978. MR 80d:65016
  • 12. H. Niederreiter. Random Number Generation and Quasi-Monte Carlo Methods, volume 63 of CBMS-NSF Region. Conf. Series in Appl. Math. SIAM, Philadelphia, 1992. MR 93h:65008

Similar Articles

Retrieve articles in Mathematics of Computation of the American Mathematical Society with MSC (1991): 65C05, 65C10, 60H10

Retrieve articles in all journals with MSC (1991): 65C05, 65C10, 60H10

Additional Information

Norbert Hofmann
Affiliation: Mathematisches Institut, Universität Erlangen–Nürnberg, Bismarckstr. 1 1/2, D–91054 Erlangen, Germany

Peter Mathé
Affiliation: Weierstraß Institute for Applied Analysis and Stochastics, Mohrenstraße 39, D–10117 Berlin, Germany

Received by editor(s): September 26, 1995
Received by editor(s) in revised form: March 27, 1996
Article copyright: © Copyright 1997 American Mathematical Society

American Mathematical Society