Abstract

We present a numerical scheme to calculate fluctuation identities for exponential Lévy processes in the continuous monitoring case. This includes the Spitzer identities for touching a single upper or lower barrier, and the more difficult case of the two-barriers exit problem. These identities are given in the Fourier-Laplace domain and require numerical inverse transforms. Thus we cover a gap in the literature that has mainly studied the discrete monitoring case; indeed, there are no existing numerical methods that deal with the continuous case. As a motivating application we price continuously monitored barrier options with the underlying asset modelled by an exponential Lévy process. We perform a detailed error analysis of the method and develop error bounds to show how the performance is limited by the truncation error of the sinc-based fast Hilbert transform used for the Wiener–Hopf factorisation. By comparing the results for our new technique with those for the discretely monitored case (which is in the Fourier-z domain) as the monitoring time step approaches zero, we show that the error convergence with continuous monitoring represents a limit for the discretely monitored scheme.

Highlights

  • Identities providing the Fourier-z transform of probability distribution functions of the extrema of a random path subject to monitoring at discrete intervals were first published by Spitzer (1956)

  • The identities for the minimum and maximum of a path, for use with a single upper or lower barrier and for the two-barrier exit problem, are comprehensively described in the discrete monitoring case by Fusai, Germano, and Marazzina (2016), who proposed numerical methods to compute them for exponential Lévy processes

  • We present results for the Spitzer–Laplace pricing procedure for continuously monitored single and double-barrier options for the normal inverse Gaussian (NIG), Kou and variance gamma (VG) processes

Read more

Summary

Introduction

Identities providing the Fourier-z transform of probability distribution functions of the extrema of a random path subject to monitoring at discrete intervals were first published by Spitzer (1956). They were extended to the continuous case by Baxter and Donsker (1957) and to double barriers by Kemperman (1963). The discretely and continuously monitored identities are in the Fourier-z and Fourier-Laplace domains respectively This means that, with the application of the inverse z or Laplace transform as appropriate, they can be used within Fourier transform option pricing methods, which we will use as an example in this paper. For example, the application to queuing systems, see the classical contributions by Cohen (1975, 1982) and Prabhu (1974) and more recent work by Bayer and Boxma (1996), Markov chains (Rogers, 1994), insurance (Chi & Lin, 2011), inventory systems (Cohen & Pekelman, 1978; Grassmann & Jain, 1989), and applied probability (Grassman, 1990), as well as in mathematical finance

Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call