This paper develops a fast Laplace transform method for solving the complex PDE system arising from Parisian and Parasian option pricing.The value functions of the options are governed by a system of partial different...This paper develops a fast Laplace transform method for solving the complex PDE system arising from Parisian and Parasian option pricing.The value functions of the options are governed by a system of partial differential equations(PDEs)of two and three dimensions.Applying the Laplace transform to the PDEs with respect to the calendar time to maturity leads to a coupled system consisting of an ordinary differential equation(ODE)and a 2-dimensional partial differential equation(2d-PDE).The solution to this ODE is found analytically on a specific parabola contour that is used in the fast Laplace inversion,whereas the solution to the 2d-PDE is approximated by solving 1-dimensional integro-differential equations.The Laplace inversion is realized by the fast contour integral methods.Numerical results confirm that the Laplace transform methods have the exponential convergence rates and are more efficient than the implicit finite difference methods,Monte Carlo methods and moving window methods.展开更多
基金supported by National Natural Science Foundation of China(Grant No.11671323)Program for New Century Excellent Talents in University of China(Grant No.NCET-12-0922)+1 种基金the Fundamental Research Funds for the Central Universities of China(Grant No.JBK1805001)Hunan Province Science Foundation of China(Grant No.2020JJ4562)。
文摘This paper develops a fast Laplace transform method for solving the complex PDE system arising from Parisian and Parasian option pricing.The value functions of the options are governed by a system of partial differential equations(PDEs)of two and three dimensions.Applying the Laplace transform to the PDEs with respect to the calendar time to maturity leads to a coupled system consisting of an ordinary differential equation(ODE)and a 2-dimensional partial differential equation(2d-PDE).The solution to this ODE is found analytically on a specific parabola contour that is used in the fast Laplace inversion,whereas the solution to the 2d-PDE is approximated by solving 1-dimensional integro-differential equations.The Laplace inversion is realized by the fast contour integral methods.Numerical results confirm that the Laplace transform methods have the exponential convergence rates and are more efficient than the implicit finite difference methods,Monte Carlo methods and moving window methods.