Complexity

Volume 2018, Article ID 1872409, 10 pages

https://doi.org/10.1155/2018/1872409

## A Semianalytical Solution of the Fractional Derivative Model and Its Application in Financial Market

School of Mathematics, Dongbei University of Finance and Economics, Dalian 116025, China

Correspondence should be addressed to Lina Song; moc.361@gnos_n_l

Received 25 September 2017; Accepted 28 February 2018; Published 1 April 2018

Academic Editor: Thierry Floquet

Copyright © 2018 Lina Song. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Fractional differential equation has been introduced to the financial theory, which presents new ideas and tools for the theoretical researches and the practical applications. In the work, an approximate semianalytical solution of the time-fractional European option pricing model is derived using the method of combining the enhanced technique of Adomian decomposition method with the finite difference method. And then the result is introduced in China’s financial market. The work makes every effort to test the feasibility of the fractional derivative model in the actual financial market.

#### 1. Introduction

Since it appeared in 1973, the Black-Scholes model [1, 2] has become the most popular method for option pricing. This classical Black-Scholes equation is a linear parabolic partial differential equation, as follows:where is the option price, is the asset price, is the risk-free interest rate, and is the volatility for the asset price .

The Black-Scholes model provides a simple and effective technique to estimate the option price. With the researches developing further, empirical studies [3–6] proved that the financial market has the fractal characters both at home and abroad. So, the classical model established on a partial differential equation with integer-order derivatives is not enough to reflect the reality of the financial market. As everyone knows, the differential equations involving the fractional derivatives are powerful tools to study the fractal geometry and the fractal dynamics. Fractional differential equations show the superiority in modeling the important phenomena in various fields such as fluid flow, electromagnetic, acoustics, electrochemistry, cosmology, and material science. Then, whether the fractional differential equation can be applied into the financial market? The answer is Yes! Fractional derivative has the property of self-similarity and it replies to the long-range dependence better than the derivative with integer order. The superiorities are beneficial to depict the fractal structure in the financial market. In recent years, the researches on the application of fractional calculus in the financial theory show an increasing tendency. Cartea and del-Castillo-Negrete [7] presented the fractional diffusion models of the option prices in the derivatives market. In this article, the authors introduced Riemann-Liouville fractional derivative into the pricing problems of financial derivatives and established the space-fractional diffusion equations using three special processes of FMLS, CGMY, and KoBoL in markets with jumps, which provided new tools to deal with financial problems. Reference [7] made the numerical treatment of knock-out barrier options for FMLS processes and the analysis showed the effectiveness of fractional models. Finally, [7] pointed out that the models can be used to deal with American option. Based on the fractional models presented by Cartea and del-Castillo-Negrete, Marom and Momoniat [8] compared the numerical solutions of three models and analyzed European option. Subsequently, Xi and Cao [9] gave numerical solutions of the fractional American option pricing model by constructing linear complementarity problem. An explicit closed-form analytical solution of the FMLS model was found by Chen et al. [10] where the effectiveness of this model is proved by comparing with derivatives of integer order and the validity of the put-call parity is verified. Jumarie [11, 12] derived the time and space-time fractional Black-Scholes equations and gave optimal fractional Merton’s portfolio. Based on Jumarie’s theoretical derivation, Liang et al. [13] introduced a bifractional Black-Scholes model whose time-fractional derivative is produced using the heuristic arguments. Using the heuristic arguments, Chen et al. [14] used a modified Black-Scholes equation with a time-fractional derivative to govern the double barrier option and gave analytical solution. Song [15] presented a space-time fractional derivative model for European option pricing with transaction costs by the replicating portfolio.

The above works present new ideas and methods for the researches on derivatives pricing theory. This study is in the initial stage and lots of problems need to be solved. The aim of the paper is to solve the time-fractional pricing model and test the practicability of the results by the market data. The paper has been organized as follows. In Section 2, the time-fractional derivative model of option pricing is succinctly described. In Section 3, the semianalytical solution of the fractional model is derived by the enhanced technique and the finite difference method. In Section 4, the result is tested by the real data from China’s warrants and options market. Some problems and conclusions are presented in Section 5.

#### 2. Time-Fractional Black-Scholes Equation

Giona and Roman [16, 17] pointed out that the relationship between the total flow of the probability flow and the average probability density iswhere is the dimension of the fractal media. The diffusion kernel of the fractal media can be expressed aswhere is a transmission exponent and is a constant that can be determined [18]. Based on [13, 14, 19], is supposed as 1 in following work.

The change of option price is considered as a fractal transmission system. Let denote the option price with ; the relationship between and for the option problem is determined by (2), as follows:where

Then, one can establish the following time-fractional model:where is just the Riemann-Liouville fractional derivative defined as

The idea of the above transformation from a Black-Scholes equation with integer-order derivative to a fractional Black-Scholes equation has been adopted for the option problem in [13, 14, 19]. By the relationship presented in the mentioned documents, one can also obtain the time-fractional model. The analytical solution of the model containing the dividend yield and is given by [14]. Li [19] gave the numerical solutions of the time-fractional model with being 1 by the finite difference method. The main task of this paper is to present solutions of (6) with the condition and give the actual application. Under the situation, finding analytical solutions is difficult and so it is necessary to turn to semianalytical approximations for getting a solution.

The values of European call and put option in the work are supposed to be a solution of (6) with the following initial and boundary conditions, respectively:

#### 3. Semianalytical Solution

In this section, two aspects of work need to be complete. The one is to give an enhanced technique of Adomian decomposition method; the other is to solve a semianalytically solution of the pricing model by the enhanced technique and the finite difference method.

##### 3.1. Description of the Enhanced Method

Consider a general differential equationwhere is an easily or trivially invertible linear operator, is the remaining linear part, and is a nonlinear operator.

Equation (10) can be rewritten aswhere is determined by the given conditions.

The standard Adomian decomposition method [20, 21] suggests the solution can be decomposed into the infinite series and the nonlinear term is decomposed; namely,where are Adomian polynomials, defined by

Substituting decomposition series (12) into both sides of (11) gives the following relationship:

From the above equation, the standard Adomian decomposition method defines the components by the following recursive relationship:

Further, [15] introduced the convergence-control parameter and new components are constructed, defined by

A solution of (10) is an infinite series with the computable components which contains a new parameter ; namely,

Comparing (15) with (16), it is knowable that the standard Adomian decomposition method is a case of the enhanced technique. The existence of the parameter can adjust the convergence region and it is beneficial to improve the accuracy of the truncated series. It is can be identified by the analysis presented in [22–24], where Professor Liao gives the efficient way to obtain the optimal convergence-control parameters.

##### 3.2. Derivation of the Solution

The finite difference method has been extended successfully to deal with various fractional differential equations [25–27]. Under the technique, the differential equation is transformed into a difference equation and numerical solutions are finally obtained.

The space variable is discretized by the difference approximation, where and is a realistic and practical approximation to infinity. Let the step and ; denotes a solution of (6) in the space point .

For spatial derivative, the following difference approximations are used:

Substituting (18) into (6) gives

Equation (8) is discretized and gives

From (9), one can get

According to Adomian decomposition method, the linear operator is chosen and the treatment of nonlinear term is not involved for this problem. Equation (19) can be rewritten as

The inverse operator of is assumed to be a Riemann-Liouville fractional integral operator, given by

Then, the recursive relationship of the call option is defined by

For the put option, we have

The components are determined by the following expressions:

Eventually, the semianalytical approximate solution of the fractional model (6) is

We can also establish other modes according to the recurrence formulas presented by [28]. But, the focus of the work is to explain the application of the time-fractional model, so we do not consider other forms. Further, the initial and boundary conditions that are satisfied by series solution (27) are checked out. For the boundary condition, it is evident via Maclaurin series expansion of and . The following equality tells us that the initial condition can be proved under condition (23). Therefore, series solution (27) is reasonably structured.

#### 4. Application and Analysis

With the help of the mathematical calculation software* MATLAB*, the result of Section 3 is tested in China’s financial market.

*Example 1. *The call warrant of Angang Steel 030001 is taken as an example to explain the application of results (27).

The data (the data came from the database of the Tai’an CSMAR series) from 12/05/2005 to 10/25/2006 are used to estimate parameters. The first step is to determine the parameter of the fractional derivative model. Theoretically speaking, it is a Hausdorff dimension. For the self-affine processes, [29] pointed out the celebrated relationship for a self-affine surface in -dimensional space, where is a fractal dimension and is Hurst coefficient. For this problem, Mandelbrot gave some discussions in [29], where he arranged a lot of his own and his partners work. Hurst is a statistic that characterizes the correlation effect of time series. For a time series, the fractal dimension can be considered as 2-H. In practical application, the calculation results are generally approximate. So, the traditional R/S analysis [30, 31] is employed to give Hurst exponent. The formula is applied to calculate the logarithmic return rate of the call warrant, where is the daily closing price. Hurst exponent for option is estimated as 0.6207. The result is greater than 0.5 and proves that the warrants market has the fractal character. This conclusion basically tallies with those studies [5, 6] on Chinese financial market.

Secondly, is estimated by the historical volatility and then is annualized. The formula iswhere stands for the daily closing price of the stock of Angang Steel and the value 242 is chosen as a general number of trading days. Then the volatilities are 0.3948. Hurst exponent is 0.6006 and takes 1.3994.

The risk-free interest rate takes 0.0252 (the data came from the people’s bank of China) that is the rate for one-year term deposit after adjustment in 2006. The strike price is 3.6. The space length and the step length are and 0.01, respectively.

After completing the above parameters estimation, the data , , and on 10/25/2006 are substituted into the second approximation . Let be (1-H), the parameter is solved as and . If is (Hurst exponent), the corresponding is and .

The next step is to calculate the option price on 10/26/2006, when , . Based on the convergence-control parameter pointed out by [22–24], the -curves of the solution with being and are drawn in Figure 1. The interval where the solution curves are near horizontal is valid and the homologous solutions can be convergent. Although Figure 1 gives the curves for the real values of , it is knowable that the two sets of values for are acceptable.

Finally, substituting and the first set of values for into , the obtained warrant values in turn are 3.304 and under the condition that is 0.3793. Similarly, one can get 3.304 and for being 0.6207. The classical Black-Scholes solution is 2.9486. The actual value is 3.256. Enlarging the analyzing scope, the estimations by the paper algorithm, the Black-Scholes model, and the actual values of the call warrant on 09/18/2006–11/27/2016 are listed in Figures 2 and 3 under the determined parameters.

The solid line and the diamonds in Figures 2 and 3 are obtained under enhanced technique (26). Actually, our results are complex-valued. The real parts of solutions are cut out and described by the solid line and the diamonds. is one of important parameters in the algorithm and it is determined by the actual values of the previous trading day. That is to say, the function of is not only a control parameter of the series solution, but also a reference point for the fluctuation of the option price on the considered trading day. From the concrete example, the valuations by substituting these into the expressions are close. From Figures 2 and 3, we can know that the solid line and the diamonds are very intimate. The trajectory motions of the two curves are along the actual values which are labeled with stars and came from the actual market. In contrast, the dots that describe the values from the classical Black-Scholes model are far from the true values. It is evident that the estimates determined by results (27) are dependable.