Asian Option Pricing with Transaction Costs and Dividends under the Fractional Brownian Motion Model

The pricing problem of geometric average Asian option under fractional Brownian motion is studied in this paper. The partial differential equation satisfied by the option’s value is presented on the basis of no-arbitrage principle and fractional formula.Then by solving the partial differential equation, the pricing formula and call-put parity of the geometric average Asian optionwith dividend payment and transaction costs are obtained. At last, the influences of Hurst index and maturity on option value are discussed by numerical examples.


Introduction
Option pricing theory has been an unprecedented development since the classic Black-Scholes option pricing model [1] was proposed.Asian options are a kind of common strong path-dependent options, whose value depends on the average price of the underlying asset during the life of the option.Fusai and Meucci [2] have discretely studied Asian option pricing problem under the Levy process.Vecer [3] has got the unified algorithm of the Asian option value based on the basic theory of stochastic analysis.Večeř and Xu [4] extended the method of removing path correlation to the case in a semimartingale model and obtained the partial differential equation of the option value under the standard Brownian motion.
However, the empirical analysis shows that there is a longterm correlation between the underlying asset prices, so that the geometric Brownian motion is not considered as an ideal tool to describe the process of asset price.Since fractional Brownian motion has the properties of self-similarity, thick tail, and long-term correlation, that fractional Brownian motion has become a good tool to depict the process of underlying asset price.According to the standard Brownian motion, Mandelbrot and Van Ness [5] obtained a stochastic integral form of fractional Brownian motion.Based on the wick product, Duncan et al. [6] introduced fractional Itô integral.Elliott and van der Hoek and others [7] have studied fractional Brownian motion with Hurst parameter belonging to the interval (1/2, 1), and they obtained fractional Girsanov theorem and fractional Itô formula by using wick product.Because an Itô formula for generalized functionals of a fractional Brownian motion with arbitrary Hurst parameter [8] was obtained by Christian Bender, it brought great convenience to option pricing.
In the reality of the securities market, investors were faced with considerable and nonignorable transaction costs and Leland [8] firstly examined the problems of option pricing and hedging with transaction costs.Due to infinite variation of geometric Brownian motion, transaction costs would become infinite in the continuous time completely hedging strategy.So Leland suggested that no-arbitrage assumption is replaced by Delta hedging strategy under the condition of discrete time occasions and transaction costs.The model was then extended by Hoggard et al. and others [9].Guasoni [10] studied the standard option with transaction costs under the fractional Brownian motion, but he did not obtain option pricing formula.Then Liu and Chang and others [11] extended the option pricing with transaction costs under fractional Brownian motion and provide an approximate solution of the nonlinear Hoggard-Whalley-Wilmott equation.have systematically discussed the European option pricing problems with transaction costs and long-range dependence.But these studies are usually aimed at European standard options.To the authors' knowledge, there does not exist systematic research about Asian options under time-varying fractional Brownian motion.
In this paper, Asian option pricing problems with transaction costs and dividends under fractional Brownian motion are studied.Firstly, the partial differential equation satisfied by geometric average Asian option value is obtained on the basis of no-arbitrage principle.Then the analytic expressions of option value and parity formula are presented by solving the partial differential equation.At last, the influences of Hurst exponent and maturity on option value are discussed by numerical examples.

Geometric Average Asian Options Pricing Model under Fractional Brownian Motion
Definition 1 (see [17]).Let (Ω, , ) be a complete probability space on which a standard fractional Brownian motion with Hurst exponent  (0 <  < 1) is continuous, centered Gaussian processes {  (),  ≥ 0} with covariance functions Cov(  (),   ( In 2003, Bender has obtained an Itô formula for generalized functionals of a fractional Brownian motion with arbitrary Hurst parameter [18].The following lemma is obtained by using the integral Itô formula.Lemma 2. Suppose that stochastic process   satisfied the following equation: where   and   are, respectively, drift coefficient and diffusion coefficient.Suppose that stochastic process  = (,  ,   ); then, for any  ∈ [0, ], one has where   =  (1/) ∫  0 ln    is geometric average of   between the time period of [0, ].
In this paper, the following basic assumptions were needed.
(i) Underlying asset price,   , satisfied the stochastic differential equations where   is the expected return,   denotes dividend yield,  is volatility, and   () is a fractional Brownian motion.
(ii) Risk-free interest rate   is a certain function of time .
(iii) Transaction costs are proportional to the value of the transaction in the underlying.Let  denote the transaction cost per unit dollar of transaction, where  is a constant.To buy or sell ]  shares of the underlying asset need pay proportional transaction costs (|]  |  ); note that ]  > 0 denotes buying the underlying asset and ]  < 0 denotes selling.
(iv) The expected return of the hedge portfolio equals the risk-free rate   .
Let  = (,   ,   ) denote the value of the geometric average Asian call at time , where   =  (1/) ∫  0 ln    is geometric average of underlying asset in [0, ].Construct a portfolio Π: long one position of the geometric average Asian call, and sell Δ shares of the underlying asset.Then the value of the portfolio at time  is After the time interval , the change in the value of the portfolio Π is as follows: where   denotes the change in the underlying asset price and ]  = Δ + −Δ  is the change of the underlying asset share in [,  + ].Choose Δ  = /  ; then, (5) becomes where The mathematical expectation of transaction costs is obtained in the following form: By assumption (iv), the following relation holds: By ( 6) and ( 8), one has Substituting ( 10) into ( 9), the following partial differential equation is obtained: Substituting   =  (1/) ∫  0 ln    and   / =   ln(  /  )/ into (11), the following equation is obtained: Substituting ( 13) into (12), the following result is obtained.
Theorem 3. Suppose that the underlying asset price   satisfied (3); then, the value of the geometric average Asian call at time  (0 ≤  ≤ ), (,   ,   ), satisfies the following mathematical model: Remark 4. Theorem 3 is obtained for the long position of the option.If the short position of option is considered, similarly, we can also get the mathematical model ( 14) and only the corresponding modified volatility is given by the following form: Let Le() = √2/(/)() −1 , which is called fractional Leland number [8].
Remark 5.For the long position of a single European Asian option, its final payoff is (  − ) + or ( −   ) + and they are both convex function, so   > 0, and noticing   =  (1/) ∫  0 ln    , thus   > 0. However, for the short position of a single European Asian option, its final payoff at maturity is −(  − ) + or −( −   ) + and they are both concave function, so that   < 0,   < 0. So for a single European Asian option, ( 13) and ( 15) can be represented as

Option Pricing Formula
Theorem 6. Suppose that the underlying asset price   satisfied (3); then, the value, (,   ,   ), of the geometric average Asian call with strike price , maturity , and transaction fee rate  at time  is where Proof.By Theorem 3, the value, (,   ,   ), of the geometric average Asian call satisfies the following model: Combined with the boundary conditions of the call option, (,   ,   ) = (  − ) + , the model ( 19) can be converted to Let  = (),   =   + (), (,   ) = (,   ) () which satisfied the conditions and then we have Substituting ( 23) into (21), we can get Set Combining with the terminal conditions () = () = () = 0, we have where Thus the model ( 21) is converted into the classic heat conduction equation Its solution is ) . ( After variable reduction, we have where So the value of geometric average Asian call option at time  is obtained Theorem 7. Suppose the underlying asset price   satisfies (3); then the relationship between   (,   ,   ), the value of geometric average Asian call option, and   (,   ,   ), the value of put option with strike price , maturity , and transaction fee rate  at time , is where  * ,  * are the same as above. Proof.Let Then  is suitable for the following terminal question in {0 ≤  < ∞, 0 ≤  < ∞, 0 ≤  ≤ }: We let   = (1/)[ ln   + ( − ) ln   ]; then  is suitable for the following in {  ∈ , 0 ≤  ≤ }: Set the form solution of problem of (36) is Substituting (37) into (36) and comparing coefficients, one has Taking () = 1, () = −, the solutions of (38) are Thus The parity formula between call option and put option is Pricing formula of the geometric average Asian put options can be obtained by Theorems 6 and 7.
Theorem 8. Suppose   satisfies (3); then the value,   (,   ,   ), of the geometric average Asian put option with strike price , maturity , and transaction fee rate  at time  is In particular, if  = 1/2,  = 0, , ,  are all constant and price formula (17) is reduced to the following formula [19]: where The following results are consequences of the above theorems. where The rest of symbols are the same as Theorem 6.
Noticing that if  = 1/2,   () is standard Brownian motion (), the corresponding underlying asset price follows geometric Brownian motion, one has the following results. where ). (48)

Numerical Example
We discuss the impact of Hurst index and transaction rates on the Asian option value by numerical examples.Assume that the parameter selection is as follows: We calculate the value of the option by using the price formula of (45).The relationships between the value of call option or put option and the underlying asset price with different Hurst index are given in Figures 1 and 2, respectively.From Figures 1 and 2, the relationship between Hurst index and Asian option value is negative.Furthermore, the impact on the call option value decreases with the increase of the underlying asset price, but the impact on the put option value decreases with the decrease of the underlying asset price.By Figures 3 and 4, we can get the change trend of the value of Asian call and put options with the change of maturity and Hurst index at the same time.The option value increases with the maturity increases, but the value of the call option increases faster than the value of a put option increase.

Conclusions
Asian options are popular financial derivatives that play an essential role in financial market.Pricing them efficiently and accurately is very important both in theory and practice.We have investigated geometric average Asian option valuation problems with transaction costs under the fractional Brownian motion.By no-arbitrage principle and fractional Itô's formula, this paper has deduced PDE satisfied by the option's value.Meanwhile, the pricing formula and call-put parity of the geometric average Asian option with transaction costs are derived by solving PDE.At last, the influences of Hurst exponent and maturity on option value are discussed through numerical examples.

Figure 1 :Figure 2 :Figure 3 :Figure 4 :
Figure 1: The values of the call option with different .

Corollary 9 .
If risk-free interest rate , dividend yield , and volatility  are all constant, then the price formulas of geometric average Asian call and put option with strike price , maturity , and transaction fee rate  under fractional Brownian motion at time  are, respectively,

Corollary 10 .
If risk-free interest rate   and dividend yield   are the functions of time  and volatility  is constant, then the price formula (,   ,   ) with respect to geometric average Asian call option with strike price , maturity , and transaction fee rate  under standard Brownian motion at time  is  (,   ,