Sinh-acceleration: efficient evaluation of probability distributions, option pricing, and Monte Carlo simulations

From MaRDI portal
Publication:5377002

DOI10.1142/S0219024919500110zbMATH Open1411.91615arXiv1808.05295OpenAlexW2885264721WikidataQ128432404 ScholiaQ128432404MaRDI QIDQ5377002FDOQ5377002


Authors: Svetlana Boyarchenko, Sergei Levendorskiĭ Edit this on Wikidata


Publication date: 21 May 2019

Published in: International Journal of Theoretical and Applied Finance (Search for Journal in Brave)

Abstract: Characteristic functions of several popular classes of distributions and processes admit analytic continuation into unions of strips and open coni around mathbbRsubsetmathbbC. The Fourier transform techniques reduces calculation of probability distributions and option prices to evaluation of integrals whose integrands are analytic in domains enjoying these properties. In the paper, we suggest to use changes of variables of the form xi=sqrt1omega1+bsinh(sqrt1omega+y) and the simplified trapezoid rule to evaluate the integrals accurately and fast. We formulate the general scheme, and apply the scheme for calculation probability distributions and pricing European options in L'evy models, the Heston model, the CIR model, and a L'evy model with the CIR-subordinator. We outline applications to fast and accurate calibration procedures and Monte Carlo simulations in L'evy models, regime switching L'evy models that can account for stochastic drift, volatility and skewness, and the Heston model. For calculation of quantiles in the tails using the Newton or bisection method, it suffices to precalculate several hundred of values of the characteristic exponent at points of an appropriate grid ({em conformal principal components}) and use these values in formulas for cpdf and pdf.


Full work available at URL: https://arxiv.org/abs/1808.05295




Recommendations




Cites Work


Cited In (6)





This page was built for publication: Sinh-acceleration: efficient evaluation of probability distributions, option pricing, and Monte Carlo simulations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5377002)