Numerical Computation of First-Passage Times of Increasing Lévy Processes |
| |
Authors: | Mark Veillette Murad S Taqqu |
| |
Institution: | (1) Division of Applied Mathematics, Brown University, Providence, RI 02912, USA;(2) Center for Nonlinear Analysis, Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA 15213-3890, USA |
| |
Abstract: | Let {D(s), s ≥ 0} be a non-decreasing Lévy process. The first-hitting time process {E(t), t ≥ 0} (which is sometimes referred to as an inverse subordinator) defined by $E(t) = \inf \{s: D(s) > t \}$E(t) = \inf \{s: D(s) > t \} is a process which has arisen in many applications. Of particular interest is the mean first-hitting time
U(t)=\mathbbEE(t)U(t)=\mathbb{E}E(t). This function characterizes all finite-dimensional distributions of the process E. The function U can be calculated by inverting the Laplace transform of the function (U)\tilde](l) = (lf(l))-1\widetilde{U}(\lambda) = (\lambda \phi(\lambda))^{-1}, where ϕ is the Lévy exponent of the subordinator D. In this paper, we give two methods for computing numerically the inverse of this Laplace transform. The first is based on
the Bromwich integral and the second is based on the Post-Widder inversion formula. The software written to support this work
is available from the authors and we illustrate its use at the end of the paper. |
| |
Keywords: | |
本文献已被 SpringerLink 等数据库收录! |
|