Option pricing of geometric Asian options in a subdi ﬀ usive Brownian motion regime

: In this paper, pricing problem of the geometric Asian option in a subdi ﬀ usive Brownian motion regime is discussed. The subdi ﬀ usive property is manifested by the random periods of time, during which the asset price does not change. Subdi ﬀ usive partial di ﬀ erential equations for geometric Asian option are derived by using delta-hedging strategy. Explicit formula for geometric Asian option is obtained by using partial di ﬀ erential equation method. Furthermore, numerical studies are performed to illustrate the performance of our proposed pricing model.


Introduction
For modeling of fluctuations in movement of underlying asset price, Brownian motion has been used traditionally as the driving force [1].However, based on some empirical studies it has been shown that financial data exhibiting periods of constant values which this property can not represent by Brownian motion [2].In recent years, many researchers attempt to fix this gap by using subdiffusive Brownian motion.Magdziarz [3] introduced subdiffusive geometric Brownian motion as the model of underlying asset price, and obtained the corresponding subdiffusive Black-Scholes (BS) formula for the fair price of European options.Liang et al. [4] extended the model of [3] into a fractional regime.Based on the fractional Fokker-Planck equation, they obtained the corresponding BS formula for European options.Wang et al. [5] considered the European option pricing in subdiffusive fractional Brownian motion regime with transaction costs.They obtained the pricing formula for European options in continuous time.One can refer to [6][7][8][9] to see more about the option pricing model in subdiffusive regime.
Asian options are financial derivatives whose payoff depends on the average of the prices of the underlying asset over a pre-fixed time interval.We will denote by {S t } t∈[0,T ] the risky asset price process.
According to the payoffs on the expiration date, Asian options can be differentiated in to two classes: fixed strike price Asian options and floating strike price Asian options.The payoff for a fixed strike price Asian option is (J T − K) + and (K − J T ) + for a call and put option respectively.The payoff for a floating strike price Asian option is (S T − J T ) + and (J T − S T ) + for a call and put option respectively.Where T is the expiration date, K is the strike price, J t is the average price of the underlying asset over the predetermined interval.According to the definition of J t , Asian options can again be divided into two types: the arithmetic average Asian option, where and the geometric average Asian option, where In recent years, scholars considered Asian option pricing under different models.Prakasa Rao [10] studied pricing model for geometric Asian power options under mixed fractional Brownian motion regime.They derived the pricing formula for European option when Hurst index H > 3 4 .Mao and Liang [11] discussed geometric Asian option under fractional Brownian motion framework.They derived a closed form solution for geometric Asian option.Zhang et al. [12] evaluated geometric Asian power option under fractional Brownian motion frame work.
To the best of our knowledge, pricing of geometric Asian option under subdiffusive regime has not been considered before.The main purpose of this paper is to evaluate the price of Asian power option under a subdiffusive regime.
The rest of the paper proceeds as follows: In section 2, the concept of subdiffusive Brownian motion and basic characteristics of inverse α−stable subordinator are introduced.In section 3, the subdiffusive partial differential equations for geometric Asian option and the explicit formula for geometric Asian option are derived.In section 4, some numerical results are given.

Subdiffusive Brownian motion
Let B(t) is a standard Brownian motion, then B(T α (t)) is called a subdiffusive Brownian motion.Where T α (t) is the inverse α-stable subordinator defined as below here U α (τ) τ≥0 is a strictly increasing α-stable Lévy process [13][14] with Laplace transform: ,where d = denotes "is identical in law to".The moments of the considered process can be found in [15] , where Γ(•) is the gamma function.Moreover, the Laplace transform of T α (t) equals where the function is the Mittag-Leffler function [16].Specially, when α ↑ 1, T α (t) reduces to the "objective time"t.

Pricing the geometric average asian options under subdiffusive Brownian motion model
In this section, under the framework of PDE method and delta-hedging strategy, we will discuss the pricing problem of the geometric Asian call options in a subdiffusive environment.
Consider a subdiffusive version of the Black-Scholes model, i.e., a simple financial market model consists of a risk-less bond and a stock,whose price dynamics are respectively given by with constant interest rate r > 0. The stock price S t = X (T α (t)), in which X(τ) follows where µ, σ are constants.
In addition, we assume that the following assumptions holds: (i) There are no transaction costs, margin requirements, and taxes; all securities are perfectly divisible; there are no penalties to short selling; the stock pays no dividends or other distributions; and all investors can borrow or lend at the same short rate.
(ii) The option can be exercised only at the time of maturity.The value of a geometric average Asian call option at time t is function of time and of S t and J t , that is, Then we can obtain the following result.Theorem 3.1.Suppose the stock price S follows the model given by Eq (3.2), the price of the geometric average Asian call option V(t, S , J) satisfies the following PDE.
with the terminal condition where Proof of Theorem 3.1 Consider a replicating portfolio Π consists with one unit option V t = V(t, S t , J t ) and ∆ units of stock.At time t the value of this portfolio is where to simplify the notation we omit t.
Suppose that ∆ does not change over the time interval (t, t + dt), then we will select appropriate ∆ and make Π is risk-free over the time interval (t, t + dt).
It follows from [5] that the differential of portfolio Π can be expressed as From Eq (1.1) we know Substituting Eq (3.7) into Eq (3.6) we have This equation is written on the pair (S t , J t ), but since it has a continuous distribution with support on (0, +∞) × (0, +∞), it follows that the function V(t, S , J) solves the PDE (3.3) and this completes the proof.
Proof is completed.Solving the terminal value problem of partial differential Eqs (3.3) and (3.4), we can obtain Theorem 3.2.If the stock price S follows the model given by Eq (3.2),then the price of a fixed strike geometric average Asian call option V(t, S , J) is given by where the function Φ(x) is the cumulative probability function for a standard normal distribution and Proof of theorem 3.2 We make the same transformation of variables as [12] and V(t, S , J) = U(t, ξ). (3.11) Though calculating we can get the following Cauchy problem (3.12) Furthermore, we apply the following transformation of variables Substituting Eqs (3.13) and (3.15) into Eq (3.12) we can obtain By calculating we have Then from Eqs (3.13) and (3.15), the cauchy problem Eq (3.12) changes into (3.16) According to heat equation theory, the solution to Eq (3.16) is given by By utilizing the inverse transformation of variables and algebraic operation to Eq (3.17), we can obtain the pricing formula of a fixed strike geometric average Asian call option Eq (3.9).Remark 3.1.It is not to obtain that γ (t) < 0, thus the maximum value of γ(t) is γ(0) = We can obtain the pricing formula for a floating strike geometric average Asian call option by using Fourier transform.Please refer to the Appendix for further details.
Furthermore, we can obtain the following Theorem.
Theorem 3.3.The put-call parity relationship for the fixed strike geometric average Asian option can be given by V(t, S , J) − P(t, S , J) = a(t)J where P(t, S , J) is the price of the fixed strike geometric average Asian put option and Proof of Theorem 3.3 Letting H(t, S , J) = V(t, S , J) − P(t, S , J).

Numerical analysis
In this section, we will give some numerical results.Theorem 3.2 leads the following results.
Corollary 4.1.When t = 0, the price of a fixed strike geometric average Asian call option V(K, T ) is given by where , Corollary 4.2.Letting α ↑ 1, then from Theorem 3.2 we can obtain the price of a fixed strike geometric average Asian call option V(t, S , J) is given by V(t, S , J) = J t S (T −t) 1 this is consistent with the result in [10,12].par From figure 1, we can see that the value of V α (K, T ) decreases with K increases and increases with T increases.
From Figure 2, we can see that the price of a fixed strike geometric average Asian call option in subdiffusive regime (V α (K, T )) is decrease with the increase of α.Furthermore, when α ↑ 1, reduces to the "objective time"t, then the price of a fixed strike geometric average Asian call option in subdiffusive regime is larger than the price of a fixed strike geometric average Asian call option in Brownian motion regime.
From Figure 3, it is obviously to see that the price of a fixed strike geometric average Asian call option in subdiffusive regime (V α (K, T )) is usually larger than that in Brownian motion regime V B (K, T ).Specially, the value of V α (K, T ) − V B (K, T ) decrease as K increases.

Conclusion
The Asian options have been traded in major capital markets, and offer great flexibility to the market participants.Therefore, it is important to price them accurately and efficiently both in theory and practice.In order to capture the periods of constant values property in the dynamics of underlying asset price, this paper discuss the pricing problem of the fixed strike geometric Asian option under a subdiffusive regime.We derive both subdiffusive partial differential equations and explicit formula for geometric Asian option by using delta-hedging strategy and partial differential equation method.Furthermore, numerical studies are performed to illustrate the performance of our proposed pricing model.

Figure 1 .
Figure 1.The price of a fixed strike geometric average Asian call option in subdiffusive regime V α (K, T ), according to the exercise date T and strike price K. Here, S 0 = 20,r = 0.06,σ = 0.3,α = 0.8.

Figure 3 .
Figure 3.The difference between the price of a fixed strike geometric average Asian call option in subdiffusive regime V α (K, T ) and the price of a fixed strike geometric average Asian call option in Brownian motion regime V B (K, T ), according to the exercise date T and strike price K. Here, S 0 = 20,r = 0.06,σ = 0.3,α = 0.7.