We show how spectral filtering techniques can improve the convergence of numerical schemes which use discrete Hilbert transforms based on a sinc function expansion, and thus ultimately on the fast Fourier transform. This is relevant, for example, for the computation of fluctuation identities, which give the distribution of the maximum or the minimum of a random path, or the joint distribution at maturity with the extrema staying below or above barriers. We use as examples the methods by Feng and Linetsky (2008) and Fusai, Germano and Marazzina (2016) to price discretely monitored barrier options where the underlying asset price is modelled by a Lévy process. Both methods show exponential convergence with respect to the number of grid points in most cases, but are limited to polynomial convergence under certain conditions. We relate these rates of convergence to the Gibbs phenomenon for Fourier transforms and achieve improved results with spectral filtering.

No 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 1 arXiv:1712.00077v1 [q-fin.CP] 30 Nov 2017 zero, we show that the error convergence with continuous monitoring represents a limit for the discretely monitored scheme.

We present new numerical schemes for pricing perpetual Bermudan and American options as well as α-quantile options. This includes a new direct calculation of the optimal exercise boundary for earlyexercise options. Our approach is based on the Spitzer identities for general Lévy processes and on the Wiener-Hopf method. Our direct calculation of the price of α-quantile options combines for the first time the Dassios-Port-Wendel identity and the Spitzer identities for the extrema of processes. Our results show that the new pricing methods provide excellent error convergence with respect to computational time when implemented with a range of Lévy processes.

No abstract