Optimal quadratic quantization for numerics: the Gaussian case

TitleOptimal quadratic quantization for numerics: the Gaussian case
Publication TypeJournal Article
Year of Publication2003
AuthorsGilles Pagès, and Jacques Printems
JournalMonte Carlo Methods and Applications
Volume9
Pagination135–166
Keywordsnumerical integration, optimal quantization, stochastic gradient methods
Abstract

Optimal quantization has been recently revisited in multi-dimensional numerical integration (see [1]), multi-asset American option pricing (see [2]), control theory (see [3]) and nonlinear filtering theory (see [4]). In this paper, we enlighten some numerical procedures in order to get some accurate optimal quadratic quantization of the Gaussian distribution in one and higher dimensions. We study in particular Newton method in the deterministic case (dimension $ d = 1 $) and stochastic gradient in higher dimensional case ($ d \geq 2 $). Some heuristics are provided which concern the step in the stochastic gradient method. Finally numerical examples borrowed from mathematical finance are used to test the accuracy of our Gaussian optimal quantizers.

References

  1. Gilles Pagès, "A space quantization method for numerical integration", J. Comput. Appl. Math., vol. 89, no. 1, Amsterdam, The Netherlands, Elsevier Science Publishers B. V., pp. 1–38, 1998.
  2. Vlad Bally, Gilles Pagès, and Jacques Printems, "A stochastic quantization method for nonlinear problems", Monte Carlo Methods and Applications, vol. 7(1), pp. 21-34, 2001.
  3. Gilles Pagès, Huyên Pham, and Jacques Printems, "An optimal Markovian quantization algorithm for multidimensional stochastic control problems", Stochastics and Dynamics, vol. 4(4), pp. 501-545, 2004.
  4. Gilles Pagès, and Huyên Pham, "Optimal quantization methods for nonlinear filtering with discrete-time observations", Bernoulli, vol. 11(5), 2005.
AttachmentSize
Gaussian Case.pdf563.34 KB