Abstract
Price of the arithmetic Asian options is not known in a closed-form solution, since arithmetic Asian option PDE is a degenerate partial differential equation in three dimensions. In this work we provide a new method for computing the continuous arithmetic Asian option price by means of partial differential equations. Using Fourier transform and changing some variables of the PDE we get a new direct method for solving the governing PDE without reducing the dimensionality of the PDE as most authors have done. We transform the second-order PDE with nonconstant coefficients to the first order with constant coefficients, which can be solved analytically.
1. Introduction
An Asian option is one whose payoff includes a time average of the underlying asset price. Asian options can be classified by their method of averaging, such as arithmetic or geometric. The geometric average Asian option is easy to price because a closed-form solution is available [1]. So in this work we only focus on arithmetic Asian option, which is the most commonly used, though an exact analytical solution for arithmetic average rate Asian options has not existed. This missing solution is primarily because the arithmetic average of a set of lognormal random variables is not lognormally distributed. There are several approaches to pricing arithmetic Asian options. The first approach is deriving approximations of closed-form solutions. Turnbull and Wakeman [2] and Levy [3] find approximate valuation formula by matching the first several moments of the arithmetic average. Geman and Yor [4] derive an analytical solution for pricing arithmetic Asian option in terms of the inverse Laplace transform. However, this transform is only applicable in some cases. A one-dimensional PDE is derived to price Asian options contingent on dividend paying stocks by Večeř [5]. Vecer and Xu [6] show that the price of arithmetic Asian option satisfies an integrodifferential equation in the case that the underlying asset is driven by special martingale processes, of which the Levy process is a special case. Dewynne and Shaw [7] provide a simplified means of pricing arithmetic Asian options by PDE approach, they derive an analytical formula for the Laplace transform in time of the Asian option, and then they obtained asymptotic solutions for Black-Scholes PDE for Asian options for low-volatility limit which is the big problem on using Laplace transform. Cruz-Báez and González-Rodrígues [8] obtain the same solution of Geman and Yor for arithmetic Asian options using partial differential equations, integral transforms, and mathematica programming, instead Bessel’s processes. Elshegmani et al. [9] derive a modified arithmetic Asian options PDE, together with its analytical solution. In addition, there are several numerical approaches to pricing arithmetic Asian options such as Monte Carlo’s simulation, Binomial Tree, and Finite Element Method.
In this paper, we derive the PDE for continuous arithmetic Asian option, and give a new method for solving this equation using Fourier transform.
2. Derivation of The PDE for Continuous Arithmetic Asian Option
We begin by assuming that the spot price of the underlying asset of the Asian option satisfies the stochastic differential equation where is a standard Brownian motion, and are constants.
We now consider continuous arithmetic Asian option with the average rate defined by the running sum of the underlying asset price or in differential form Suppose an Asian option has pay-off function at an expiration date . Then, the value of the Arithmetic Asian option at time is By multidimension of Ito’s Lemma, we have Set up a portfolio of one Asian option and a number of the underlying assets, the value of this portfolio is and the change in the value of this portfolio is To get rid of the stochastic term, choose , But we know that the change in the value of a portfolio with risk-less asset is From (2.9) and (2.8), yield its This is the Black-Scholes PDE of the continuous arithmetic Asian options.
There are four different types of the continuous arithmetic Asian options depending on the pay-off function as follows. (1)Arithmetic average fixed strike call option: (2)Arithmetic average fixed strike put option: (3)Arithmetic average floating strike call option: (4)Arithmetic average floating strike put option:
3. Analytical Solution of The PDE
To solve (2.10) we will make the first change of the variables Assume further , and taking into account that
then (3.1) becomes: Assume where is a complex number, then we get Note that all the coefficients in the above equation are constant except the coefficient with the term , so we can easily apply the Fourier transform.
Fourier’s transform for a function is defined by and the inverse Fourier transform is Some properties of the Fourier transform that we need in this work as follows: where is a constant, and is a Dirac delta function.
Applying the Fourier transform in on (3.4), Assume Then (3.9) is reduced to Applying another Fourier transform in , One solution of the above equation is From the initial condition, we have Applying the inverse Fourier transform in , To prove that expression (3.21) is a direct solution for (2.1), we will differentiate, (3.14) with respect to all variables, and then substituting into (2.1), Substituting (3.22) into (2.10) yields since we have , so
4. Conclusion
The valuation of the arithmetic Asian options with continuous sampling has been an outstanding issue in finance for several decades. Describing the distribution of the integral of lognormals is found to be challenging. In this paper, we have solved the problem with the PDE approach. We show that the governing PDE from the second order can be transformed to a simple PDE from the first order with constant coefficients, which can be easily solved. We have the solution for all types of the continuous arithmetic Asian options only by changing the pay-off function with respect to which one of the options we want to price. Our approach could be extended to the continuous arithmetic Asian options with constant dividend yield.