Lookback Option Pricing with Fixed Proportional Transaction Costs under Fractional Brownian Motion

The pricing problem of lookback option with a fixed proportion of transaction costs is investigated when the underlying asset price follows a fractional Brownian motion process. Firstly, using Leland's hedging method a partial differential equation satisfied by the value of the lookback option is derived. Then we obtain its numerical solution by constructing a Crank-Nicolson format. Finally, the effectiveness of the proposed form is verified through a numerical example. Meanwhile, the impact of transaction cost rate and volatility on lookback option value is discussed.


Introduction
In 1985, the problems of option pricing and replication with transaction costs were firstly examined by Leland [1]. Because of infinite variance of geometric Brownian motion process, transaction costs will become infinite, thus the arbitrary argument used by Black-Scholes to price options can no longer be used. Barles and Soner [2] considered the difference between the maximum utility from final wealth when there is no option liability and when there is such a liability. They derived the nonlinear Black-Scholes equation satisfied by the option value by using the utility maximization theory. Considering both transaction costs and the risk from a volatile portfolio, in 2005, Jandačka andŠevčovič [3] got another more complicated model. The influence of paying transaction costs and discontinuous change of stock price on option value was studied by Amster et al. [4]. Under a stochastic volatility process, Mariani et al. [5] investigated the option pricing problem with transaction cost and derived the pricing model. Meanwhile they got numerical solutions by using the finite difference method.
As a kind of exotic options, the holder of lookback option could look back on the evolution process of underlying assets during the life of option at maturity. It is one kind of path-dependent options where the payoff is based on the maximum or the minimum of the underlying asset price during the drift of the option. In 1979, the lookback option pricing formula was firstly given by Goldman et al. [6]. Then the pricing formula was extended by Conze and Viswanathan [7]. They obtained explicit formulas of various European lookback options and also provided some results for the American counterparts by means of probability method. With underlying asset price following geometric Brownian motion, analytic expressions of discrete lookback option value were presented by Heynen and Kat [8]. Many empirical studies have shown that distribution of asset price is not entirely lognormal, and the probability density function of its logarithm yield tends to have the feature of "fat-tailedness, " so geometric Brownian motion processes do not accord with the realistic environment. For characteristics of self-similarity, fat tails, and long-term dependence, fractional Brownian motion (fBm) becomes an effective tool to describe the stock price process instead of geometric Brownian motion. fBm is neither a Markov process nor a semimartingale, so stochastic analysis theory for semimartingale cannot be applied directly. In incomplete markets, arbitrage opportunities will exist. A considerable number of arbitrage strategies for fBm have been provided by Rogers [9], Shiryaev [10], Salopek [11], and Cheridito [12]. In 2012, Gu et al. [13] dealt with the problem of discrete time option pricing in the presence of transaction costs by a time-changed geometric fractional Brownian motion model. Through a mean self-financing 2 International Scholarly Research Notices delta-hedging argument, they obtained the pricing formula for the European call option in discrete time setting. Feng [14] discussed the lookback option under fBm model. They derived the explicit solutions of lookback put option pricing formula.
In this paper, we will firstly give the price dynamics model under fBm process. Then a three-dimensional nonlinear mathematical pricing model for lookback put option value with fixed transaction costs will be established. We reduced it to a two-dimensional model through variable substitution. But the reduced model is still nonlinear, so it is not easy to get the analytic solution. So we get its numerical solution by constructing a Crank-Nicolson numerical scheme for the transformed model. Finally, we analyze the effectiveness of the numerical scheme and the influence of parameters on lookback put option value.
The fBm has long-term dependency for ∈ (0.5, 1) and is antisustainable for ∈ (0, 0.5). For ̸ = 0.5, it is neither a Markov process nor semimartingale, so the classical Itô's calculus cannot be used to define a fully stochastic calculus for fBm. Lin [15] and Decreusefond and ustünel [16] developed the stochastic calculus of variations under fBm model, but Shiryaev [10] investigated that there exists arbitrage chance in the market. Hu and Øksendal [17] presented the fractional Itô's calculus based on the Wick integration driven by fBm. Bender [18] constructed fractional Itô's formula with this type of calculus. For more details about fBm, interested readers may refer to Hu and Peng [19].
Consider a financial market with two primitive securities, namely, a risky asset and a risk-free bond . We will need the following assumptions.
(i) The securities trading is carried out continuously with short selling allowed; (ii) Suppose V shares of the traded stock are bought (V > 0) and sold (V < 0) at the price . Then a proportion of transaction cost is given by |V | in either buying or selling, where > 0 is the fixed transaction cost rate.
(iii) There exists no arbitrage opportunities and the portfolio's expected rate is the risk-free interest rate , where is a constant.
(iv) The portfolio is revised every time interval , where is a finite and fixed time step.
The price process of risk-free debt evolves over time as Suppose that the price process of the risky share evolves over time as where indicates the volatility of the stock price.

The Valuation of Lookback Put Option with Transaction Costs under fBm
In this section, we adopt Leland hedging strategy to derive a model for lookback option price. The following lemma gives Fractional Itô's differential rule which takes great value in deriving the formula.
International Scholarly Research Notices 3

Theorem 3. The pricing model for the value of lookback put option with transaction costs under fractional Brownian motion process is given as
Proof. For the option purchaser, a hedging portfolio Π is constructed: one long position of lookback put option and short Δ shares of stock; then the value of the portfolio at the current time is Π = − Δ . Assume that the trade occurs only at times and + where denotes the hedging time interval. Then over a small time interval , the change in value of the portfolio will be The stochastic partial differential equation (3) can be shown as By using discrete form of Lemma 2, we can get Substituting (9) into (7) yields In order to eliminate risk, let Δ = / . So in the same interval , the traded shares of stock are Since is not differentiable, we need to approximate it as We use ( ) instead of path dependence variable ; since ≤ , we can get Consequently, from (11) The transaction costs is as follows: According to the Definition 1, we know that Therefore From the assumption (iii), we have ( Π ) = Π . Thus Consequently The above equation can be simplified as 4

International Scholarly Research Notices
Here The boundary condition is whose financial significance is that the value of option is not sensitive to the maximum when the underlying asset price reaches the maximum [20]. According to the definition of the lookback put option with floating strike price, the formula can be given as Hence the result is proved.
Since the parameter̃is nonlinear in the model (6), it's difficult to get its analytic solution. Then we will give its numerical scheme in the next section.

Crank-Nicolson Scheme
Firstly we should reduce the three-dimensional nonlinear mathematical model (6) into a corresponding twodimensional model so that it can be solved more easily. Let With the time step = / and spatial step ℎ = max / .
Time derivatives and spatial derivatives are, respectively, approximated by difference method Substituting the above equations into the model (26) gives standard Crank-Nicolson scheme. Since the corrected volatility (27) has first-and second-order derivatives, its solution can be obtained through nonlinear iteration. But it will consume too much time. We will use standard central difference scheme to approximate it. Hence we get the new numerical scheme as follows: with So the numerical scheme equals The format (33) is a six-point implicit scheme, whose local truncation error is (ℎ 2 + 2 ). Only a little computation can achieve satisfactory accuracy. For numerical convenience, we define a vector = [ 1 , 2 , . . . , −1 ] . Then the format (33) can be given as below: ( ) +1 = ( ) + , = 0, 1, 2, . . . , − 1. (35) Here We also need to know the value of the boundary points and 0 . So discrete the boundary conditions by second-order Gear formula. Hence When the price of underlying asset reaches 0, the lookback put option must be exercised. Therefore another boundary condition can be given as According to variable substitution (25), we have Substituting the boundary conditions (37) and (39) into (35) and using chasing method gives this format's solution. The value of the lookback put option at the starting time is 0 . It can be seen that it does not contain the variable in the new numerical scheme so that it avoids the difficulty of determining the variable in the iterative process.

Numerical Example
In this section, we will firstly study the convergence of the numerical scheme to explain its validity and discuss the influence of various parameters on the lookback option's value. For simplicity, we will only consider valuation of lookback put option for ∈ (0.5, 1). Example 1. Consider a European put option in half a year when the underlying asset is a stock, whose price follows Fractional Brown Motion (3), and the option value is satisfied by the model (6). Suppose that the stock price is $30, the volatility is 20% per annum, the risk-free interest rate is 5%, the Hurst parameter is 0.6, the hedging time interval is 0.005, and the number of time steps is taken as 100, spatial steps taken as 300 in the numerical scheme (35). With our usual notation, this means that For different values of Hurst parameter, Table 1 shows the price for lookback put option with transaction cost through Crank-Nicolson numerical approach and proves its convergence of this method. From Table 1, we can also see that, with the increase of Hurst parameter, the magnitude of decrease for lookback put option value is increasing. This implies the self-similarity of fBm model for ∈ (0.5, 1).
Then we study the influence of transaction cost rates on lookback put option value. Other values of parameters remain unchanged. Figure 1 displays that the value of the option is decreasing with the increase of transaction costs. This is mainly because in this paper, for the long position of options, we only consider the purchase price. Transaction rates' increase makes a unit of an asset's cost increase, directly leading to the option hedging costs' raising. So option value will decrease. Besides, with the initial stock price increased, the lookback put option value increases. Figure 2 shows the effect of expiration time and stock price on lookback put option value. From it we can see that longer maturity time can produce larger option value. Moreover, when the expiration time is shorter, it grows more obviously. This is mainly because that the option value increases over time.
Then we discuss the impact of risk-free rate and volatility change on lookback put option value. From Figure 3, it can be seen that the option value decreases with risk-free rate increasing. One factor is that the expected rate of return rises as risk-free rate grows. Another factor is that risk-free rate's growth leads to the present value of future cash flows 6 International Scholarly Research Notices  St o ck p ri ce T r a n s a c t io n r a t e Lookback put option value ×10 −3   dropping. On the other hand, the larger the volatility, the larger the option value. It is mainly because that volatility's increase amplifies the probability of making large gains in the future.
Finally, we illustrate the advantages of the model in this paper from the data in financial market. Taking China Guodian warrants for example, we will use the new model presented in this paper and Black-scholes model to reprice the warrants and obtain the error between theoretical and   Table 2 below shows the statistical characteristics of the logarithmic rate of return for China Guodian warrants. It can be observed that the value of kurtosis for the warrants is greater than that for the standard normal distribution and the probability to obtain corresponding statistic -is less than 0.05, which indicates share earnings is distributed the feature aiguille large remaining part.
By comparison between theoretical and empirical prices, Figure 4 shows pricing errors under the pricing model of lookback option with transaction costs based on fractional Brownian motion process ( = 0.7) and Black-Scholes model, respectively. As can be seen from Figure 4, the pricing errors under FBM model are significantly less than BS model. This indicates that the new model involved in this paper is more realistic.

Conclusions
In this paper, the lookback option pricing problem, with proportional transaction costs under the Fractional Brown Motion model, was studied. By using hedging principle, the nonlinear partial differential equations, satisfied by International Scholarly Research Notices 7  the option value, have been derived. As for the equations, we have constructed the Crank-Nicolson scheme and verified its effectiveness by using the Matlab software. It can be seen that the model we have built up is more in line with the reality of the market environment through Matlab software simulation.