Complexity
Volume 2018 (2018), Article ID 1872409, 10 pages
https://doi.org/10.1155/2018/1872409

Research Article

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

Lina Song

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

Correspondence should be addressed to Lina Song; l_n_song@163.com

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 is

where is the dimension of the fractal media. The diffusion kernel of the fractal media can be expressed as

where 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 equation

where is an easily or trivially invertible linear operator, is the remaining linear part, and is a nonlinear operator.

Equation (10) can be rewritten as

where 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 is

where 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.

(a)

(b)

Figure 1: The -curve.

Figure 2: .

Figure 3: .

Example 2. Considering the SSE 50 ETF put option 10000188 and the call option 10000615 as the second example, the daily closing prices (the data came from the trading software of Essence Securities) of options on 04/23/2015–12/16/2015 and 04/28/2016–12/27/2016 are taken as two groups of the sample data.
The first is to estimate parameters by the option 10000188. The strike price of the put option equals 3. Hurst exponent is estimated as 0.6567, which reflects that the time series has the long-term positive autocorrelation and indicates that the option market exhibits the fractal characters. The historical volatility for 50 ETF is 0.4385, with the number of the trading days being 242. Hurst exponent is 0.6468 and can be 1.3532. Based on the historical data from the people’s bank of China, the interest rate is made four adjustments in 2015 and the average value is around 0.02. So the value of the risk-free interest rate is considered as 0.02. and are and 0.001, successively.
is the first approximation; namely, . Repeating the steps of Example 1, the real price on 12/01/2015 is used to calculate and then evaluates the option value on 12/02/2015. In addition, the estimations from 12/03/2015 to 12/15/2015 are calculated and listed in Table 1 together with those from the Black-Scholes model and the actual values. Based on the study for the option 10000188, the option 10000615 with strike price 2.006 is estimated under the condition that the parameters , , , , , and have the same values as those of the put option. The estimations by the paper algorithm, the Black-Scholes model, and the actual values on 04/28/2016–05/12/2016 are listed in Table 2. Further, numerical comparisons from 05/13/2016 to 12/27/2016 are shown in Figures 4 and 5. From the figures, the solid lines, the diamonds, and the stars are close and the dots are over them at most of the trading time points. It can be seen that the curves of the real parts of our solutions are closer to the true values than those from the Black-Scholes formula. Namely, our errors are smaller. Through the comparative analysis, approximation (27) is stable and executable. It has embodied the superiority than the Black-Scholes model especially for the option 10000615.

Table 1: Numerical comparison for option 10000188.

Date
12/02/2015
12/03/2015
12/04/2015
12/07/2015
12/08/2015
12/09/2015
12/10/2015
12/11/2015
12/14/2015
12/15/2015

( )

()

Paper result ()

Paper result ()

Black-Scholes model
0.5385
0.5399
0.6001
0.6079
0.6317
0.6288
0.6429
0.6541
0.5834
0.6165

Actual value
0.5536
0.5507
0.6136
0.6182
0.6592
0.6422
0.6575
0.6509
0.5842
0.6179

Table 2: Numerical comparison for option 10000615.

Date
04/28/2016
04/29/2016
05/03/2016
05/04/2016
05/05/2016
05/06/2016
05/09/2016
05/10/2016
05/11/2016
05/12/2016

0
−0.05104, −0.5153
−0.06578, −0.5399
0.1661, −0.6499
0.1686, −0.6764
0.1847, −0.7002
−0.1134, −0.7773
−0.626, −1.087
−0.4913, −1.265
−0.3173, −1.131

0

0.1366, −0.5452
0.1383, −0.5656
0.1564, −0.5889

Paper result
0.166
0.1618, 0.1636
0.1862, 0.1791
0.1636, 0.1659
0.1614, 0.1622
0.1136, 0.1358
0.1013, 0.1204
0.118, 0.1187
0.1202, 0.1104
0.1208, 0.1135

Paper result ()
0.1622

0.1636, 0.1651
0.1614, 0.1618
0.113, 0.1287

Black-Scholes model
0.3826
0.3771
0.3919
0.3870
0.3847
0.3526
0.3316
0.3300
0.3334
0.3361

Actual value
0.1688
0.1625
0.1695
0.1634
0.1598
0.1328
0.1188
0.1151
0.1156
0.1105

(a) 05/13/2016–08/31/2016

(b) 09/01/2016–12/27/2016

Figure 4: .

(a) 05/13/2016–08/31/2016

(b) 09/01/2016–12/27/2016

Figure 5: .

5. Problems and Conclusions

In the work, the approximate semianalytical solution of the fractional option pricing model is derived by the enhanced technique and the finite difference method. Further, the result is introduced into the real market. By the two examples in Section 4, one can find that the application of the fractional derivative model in the actual financial market is reliable and valuable. But the existing problems should be pointed out.

(i) The first problem is the parameter estimations. Except for some conventional parameters, the fractional derivative model adds the fractal dimension and the order of fractional derivative. The two parameters make the model better display the fractal features of the market. The semianalytical approximate solution presented in this work contains a convergence-control parameter that can improve the accuracy of the approximate solution. However, the existence of these parameters increases the difficulties of the empirical analysis. In the work, equals 2-H. The relationship between and on the previous day is used to calculate option price on that day. These are not accurate. If these parameters are made more accurate, one can obtain much better valuation. Their optimal determinations are the key to settle the problems.

(ii) The second problem is the computational difficulties and the actual applications are limited. Example 2 tells us that even the first approximation is also effective. But the point influences judging the rational values of by the -curve. If , the enhanced technique is reduced to the standard Adomian decomposition method. Research experience illustrates that the values of lie in the range close to zero and can be taken as the reasonable estimations. But it is undeniable that the computational difficulties cause the limitations in real application.

Although the fractional partial differential equation has been introduced into the financial theory, the literatures on the practical application are very few. This work makes every effort to apply the fractional partial differential equation to model the actual market. For existing problems, we continue the study in further work. The connections between the fractal dimension, the order of fractional derivative, and Hurst exponent are worthy of studying. In short, the effectiveness and advantage of the fractional derivative model are indubitable, but more problems need further researches. We hope and believe that the study of the fractional derivative model in the financial theory can advance the development of mathematics and finance.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work is supported by the National Natural Science Foundation of China (no. 71501031) and the Project of the Educational Department of Liaoning Province, China (no. LN2017FW010).

References

F. Black and M. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–654, 1973.

R. C. Merton, “Theory of rational option pricing,” Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141–183, 1973.

B. Mandelbrot, “The variation of certain speculative prices,” The Journal of Business, vol. 36, no. 4, pp. 394–413, 1963.

E. E. Peters, “Fractal structure in the capital market,” Financial Analyst, vol. 7, pp. 434–453, 1989.

H. Q. Li and C. Q. Ma, “An empirical study of long-term memory of return and volatility in chinese stock markrt,” Journal of Finance and Economics, vol. 31, no. 8, pp. 29–37, 2005.

T. F. Huang, B. Y. Li, and J. X. Xiong, “Test on the chaotic characteristic of chinese futures market,” Systems Engineering, vol. 30, no.1, pp. 43–53, 2012.

Á. Cartea and D. del-Castillo-Negrete, “Fractional diffusion models of option prices in markets with jumps,” Physica A: Statistical Mechanics and its Applications, vol. 374, no. 2, pp. 749–763, 2007.

O. Marom and E. Momoniat, “A comparison of numerical solutions of fractional diffusion models in finance,” Nonlinear Analysis: Real World Applications, vol. 10, no. 6, pp. 3435–3442, 2009.

J. Xi and J. W. Cao, “Fractional partial differential equations and numerical discretization method for pricing american option,” Journal on Numerical Methods and Computer Applications, vol. 35, no. 3, pp. 229–240, 2014.

W. Chen, X. Xu, and S.-P. Zhu, “Analytically pricing European-style options under the modified Black-Scholes equation with a spatial-fractional derivative,” Quarterly of Applied Mathematics, vol. 72, no. 3, pp. 597–611, 2014.

G. Jumarie, “Stock exchange fractional dynamics defined as fractional exponential growth driven by (usual) Gaussian white noise. Application to fractional Black-Scholes equations,” Insurance: Mathematics and Economics, vol. 42, no. 1, pp. 271–287, 2008.

G. Jumarie, “Derivation and solutions of some fractional Black-Scholes equations in coarse-grained space and time: application to Merton's optimal portfolio,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1142–1164, 2010.

J.-R. Liang, J. Wang, W.-J. Zhang, W.-Y. Qiu, and F.-Y. Ren, “Option pricing of a bi-fractional Black-Merton-Scholes model with the Hurst exponent H in ,” Applied Mathematics Letters, vol. 23, no. 8, pp. 859–863, 2010.

W. Chen, X. Xu, and S.-P. Zhu, “Analytically pricing double barrier options based on a time-fractional Black-Scholes equation,” Computers & Mathematics with Applications, vol. 69, no. 12, pp. 1407–1419, 2015.

L. Song, “A space-time fractional derivative model for European option pricing with transaction costs in fractal market,” Chaos, Solitons and Fractals, vol. 103, pp. 123–130, 2017.

M. Giona and H. E. Roman, “Fractional diffusion equation on fractals: one-dimensional case and asymptotic behaviour,” Journal of Physics A: Mathematical and General, vol. 25, no. 8, pp. 2093–2105, 1992.

M. Giona and H. E. Roman, “Fractional diffusion equation for transport phenomena in random media,” Physica A: Statistical Mechanics and Its Applications, vol. 185, no. 1-4, pp. 87–97, 1992.

F.-Y. Ren, Z.-G. Yu, and F. Su, “Fractional integral associated to the self-similar set or the generalized self-similar set and its physical interpretation,” Physics Letters A, vol. 219, no. 1-2, pp. 59–68, 1996.

W. Y. Li, Application and numerical solution of fractional order equation in financial model [M.S. thesis], Hangzhou Dianzi University, Hangzhou, China, 2009.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

G. Adomian, Solving frontier problems of physics: the decomposition method, Kluwer Academic Publishers, Massachusetts, Mass, USA, 1994.

S. J. Liao, beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman and Hall/CRC Press, Boca Raton, Fla, USA, 2003.

S. Liao, “Notes on the homotopy analysis method: some definitions and theorems,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 983–997, 2009.

S. Liao, “An optimal homotopy-analysis approach for strongly nonlinear differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 2003–2016, 2010.

S. Chen, F. Liu, and V. Anh, “A novel implicit finite difference method for the one-dimensional fractional percolation equation,” Numerical Algorithms, vol. 56, no. 4, pp. 517–535, 2011.

E. Sousa, “A second order explicit finite difference method for the fractional advection diffusion equation,” Computers & Mathematics with Applications, vol. 64, no. 1, pp. 3141–3152, 2012.

Q. Zhang, M. Ran, and D. Xu, “Analysis of the compact difference scheme for the semilinear fractional partial differential equation with time delay,” Applicable Analysis: An International Journal, vol. 96, no. 11, pp. 1867–1884, 2017.

L. Song and W. Wang, “A new improved Adomian decomposition method and its application to fractional differential equations,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1590–1598, 2013.

B. B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman and Company, New York, NY, USA, 1982.

B. B. Mandelbrot, “Statistical methdology for nonperiodic cycles: from the covariance to R/S analysis,” National Bureau of Economic Research, vol. 1, no. 3, pp. 259–290, 1972.

H. E. Hurst, “The long-term storage capacity of reservoirs,” Transcactions of the American Society of Civil Engineers, vol. 116, pp. 87–92, 1951.