A Wiener-Hopf based approach to numerical computations in fluctuation theory for Lévy processes

摘要

This paper focuses on numerical evaluation techniques related to fluctuation theory for Lévy processes; they can be applied in various domains, e.g., in finance in the pricing of so-called barrier options. More specifically, with X_t:= sup_(0≤ s≤t) X_s denoting the running maximum of the Lévy process X_t, the aim is to evaluate ?X_t ε dx) for t,x>0. The starting point is the Wiener-Hopf factorization, which yields an expression for the transform E e~(-αXe(ν))of the running maximum at an exponential epoch (with ν~(-1) the mean of this exponential random variable). This expression is first rewritten in a more convenient form, and then it is pointed out how to use Laplace inversion techniques to numerically evaluate ?X _tε dx). In our experiments we rely on the efficient and accurate algorithm developed in den Iseger (Probab Eng Inf Sci 20:1-44, 2006). We illustrate the performance of the algorithm with various examples: Brownian motion (with drift), a compound Poisson process, and a jump diffusion process. In models with jumps, we are also able to compute the density of the first time a specific threshold is exceeded, jointly with the corresponding overshoot. The paper is concluded by pointing out how our algorithm can be used in order to analyze the Lévy process' concave majorant.