Special Issue

## Discrete Fractional-Order Systems with Applications in Engineering and Natural Sciences

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 1393456 | https://doi.org/10.1155/2020/1393456

Juan He, Aiqing Zhang, "Finite Difference/Fourier Spectral for a Time Fractional Black–Scholes Model with Option Pricing", Mathematical Problems in Engineering, vol. 2020, Article ID 1393456, 9 pages, 2020. https://doi.org/10.1155/2020/1393456

# Finite Difference/Fourier Spectral for a Time Fractional Black–Scholes Model with Option Pricing

Accepted27 Jul 2020
Published04 Sep 2020

#### Abstract

We study the fractional Black–Scholes model (FBSM) of option pricing in the fractal transmission system. In this work, we develop a full-discrete numerical scheme to investigate the dynamic behavior of FBSM. The proposed scheme implements a known formula for the -order fractional derivative and Fourier-spectral method for the discretization of spatial direction. Energy analysis indicates that the constructed discrete method is unconditionally stable. Error estimate indicates that the -order formula in time and the spectral approximation in space is convergent with order , where is the regularity of and and are step size of time and degree, respectively. Several numerical results are proposed to confirm the accuracy and stability of the numerical scheme. At last, the present method is used to investigate the dynamic behavior of FBSM as well as the impact of different parameters.

#### 1. Introduction

The classical option pricing model is proposed by Black and Scholes , which is based on the assumption that stocks and options are in an “ideal state” in the market and Samuelson’s model :where be the the stock value, is the Brownian motion with the unit variance, and and are two constants. But in many cases, fractional Brownian motion is more accurate than integer order . On the other hand, more and more diffusion processes were found to be non-Fickian [4, 5], and the fractional order stochastic differential equation is considered as an extension of the stochastic differential equation. One view is that fractional order option trading equation is regarded as nonrandom growth process caused by Brownian motion. Therefore, Jumarie  and Liang et al.  considered fractional Brownian motion in Samuelson’s model equation:

Combining Itô lemma and fractional Taylor expansion of the option price , Jumarie  obtained the following FBSM:

Jumarie  derived new families of the exact solution of the above equations. Moreover, traditional pricing models for double barrier options are often biased when price changes are considered as fractal transmission systems. Chen et al.  revealed that it would be better to use fractional order Black–Scholes equation to explain the pricing in fractal transmission systems.

There is a lot of work in the modeling and calculation of fractional equations. Yang et al. [10, 11] developed a new definition of fractional derivative. The advantage of this definition is that it does not contain singular kernel. Inc et al. [12, 13] studied the isolated solutions of a class of fractional equations with Kerr law nonlinearity by Riccati–Bernoulli method. The exact dark optical and periodic singular soliton solution is obtained. Singhet al.  solved a series of fractional equations by using homotopy analysis technique and Laplace transform algorithm.

As we all know, as an effective method, formula  has been widely used in the calculation of fractional differential equations. Langlands and Henry , Sun and Wu , and Lin and Xu  discussed the error estimate for scheme. De Staelena and Hendybc  constructed a numerical method of fourth-order finite difference in space and in time. Stability, uniqueness, and error estimates are analyzed. Zhang et al.  presented a discrete implicit finite scheme to solve time fractional Black–Scholes model. They discussed the stability and error estimation of numerical schemes by Fourier analysis. Unfortunately, their analysis methods are local approach. Due to the importance of FBSM, it is necessary to reconstruct an efficient numerical method and analyze global stability and error estimates.

In this work, we will develop an efficient full-discrete scheme to approximate the Black–Scholes model with -order fractional derivative. We apply method to discretize the direction of time and Fourier-spectral method to discretize the direction of space. Using the energy analysis method, we discuss the stability and error estimate of the fully discrete numerical method. The detailed analysis shows that the scheme is unconditionally energy stable, and the error estimates indicate that our full-discrete scheme can achieve -order accuracy in time and exponential accuracy in space direction. Finally, some numerical examples are conducted to support the theoretical claims. At the same time, the dynamic behavior of FBSM is studied by the proposed method.

We organize the rest of the paper as follows. Section 2 will briefly introduce the FBSM. In Section 3, we develop a time-discrete method for FBSM and then present its discrete energy law. In Section 4, we will study the error estimate of the full-discrete scheme. In Section 5, we present accuracy/stability tests and numerous numerical examples to demonstrate the validity of the full-discrete method. In addition, we will discuss the properties of the solution of the FBSM. Some concluding remarks are given in Section 6.

#### 2. Black–Scholes Model

In this work, we will consider the following time fractional Black–Scholes model:where is the price of the option, the price of the underlying asset, the interest rate, and the volatility of the stock price, ; the time fractional derivative is defined by

This is a linear parabolic partial differential equation which has been studied extensively.

We transform the problem to an initial value problem by using the time to mature , and we then set ; we can rewrite (4) aswhere , and with the following boundary (barrier) and initial conditions,

In order to solve the above model by numerical method, it is necessary to truncate the original unbounded region into a finite interval. Therefore, we will consider problem (8) in bounded interval . Then, we will study the following problem:

Remark 1. In fact, one can choose homogeneous or inhomogeneous boundary conditions. It all depends on the actual option price. We have tested it, and it does not make any difference in actual numerical examples.

#### 3. Order Numerical Method

Here, we will develop the time-discrete method for equation (10). First, given a positive integer , set be the time step size, and denote as the mesh point. Then, we introduce an method to discrete the Caputo fractional derivative of order :where .

Lemma 1 (see [20, 23, 24]). The coefficients of formula (13) satisfy the following properties:

Then, we can obtain the following time-discrete scheme:where . It should be noted that if , we can rewrite the above equation as

First of all, we have the following energy stability results for time-discrete (15).

Theorem 1. The time-discrete scheme (15) is unconditionally stable. It satisfies the following energy dissipation law:

Proof. When , computing the inner product of (16) with , we obtainIt is easy to verify that the following formula is correct:Thus,Giving up some positive terms, we haveAssume the following inequality holds:Next, we will show is still valid. If , taking the inner product of (15) with , we deriveNote the fact thatThus, we getThis yields (17).

#### 4. Error Estimate for Full Discretization

In this part, we will study the Fourier-spectral method for the time-discrete method (15). First, we define as the polynomial space. Define be the -projection operator which satisfies

We have the following estimate :

Then, we can develop the following full-discrete scheme:

We now present the stability results of the fully discrete scheme (28).

Theorem 2. Let be the solution of (28), then we derive

Next, we begin to analyze the error estimates of the full-discrete scheme (28). Define the following error function:

From [20, 19, 24, 23], we know that satisfies

We also define the following error functions:

Lemma 2. For and , we have the following results:

Proof. Note thatTherefore, we obtain

Remark 2. In [23, 24], readers will also find similar parameter estimation. This result is very useful for us to analyze error estimate.

Theorem 3. For the constructed numerical scheme (28), we have the following error estimate:

Proof. For , we can write equation (28) asSubtracting (38) from (10) at , we note thatThen, we haveSet , we haveDropping some positive terms, we findAssumeNext, we will prove that it holds also for . Subtracting (28) from a reformulation of (10) at , we findLet , we haveThus, we haveNote that , thusNote thatThis ends the proof.

#### 5. Numerical Examples

In this section, several numerical examples will be present to confirm the accuracy and applicability of the full-discrete scheme (28). We consider a rectangular computed domain of . In order to better simulate the periodic boundary conditions, the Fourier-spectral method will be used to discretize space direction.

##### 5.1. Verification of Convergence of Numerical Method

First, in order to conduct a time accuracy test, an exact solution will be constructed to evaluate the convergence of the full-discrete scheme (28).

Example 1. We consider the following FBSM with -order Caputo derivative:whereIt is easy to verify that the exact solution will be .
We set . The default values for the parameters are set as . In Table 1, we show the temporal convergence orders of various time steps. As can be seen from Table 1, our full-discrete scheme is close to -order accuracy in time, which is confirm with the result in Theorem 3.

 = 1/200 = 1/400 = 1/800 = 1/1600 = 1/3200 1.8166 1.8263 1.8345 1.8415 1.8475 1.6680 1.6746 1.6797 1.6837 1.6869 1.4892 1.4925 1.4947 1.4963 1.4974 1.3939 1.3960 1.3974 1.3983 1.3989 1.2964 1.2979 1.2987 1.2992 1.2995 1.0980 1.0990 1.0995 1.0997 1.0999

Fix , in Figure 1, we give error for different . It is obvious that our numerical scheme has good convergence in time direction. Let , and . Figure 2 shows that the full-discrete scheme (28) has excellent convergence behavior in space direction.

##### 5.2. Effect of Various Parameters

This section is devoted to investigate the dynamic behavior of FBSM equation with different . In the following numerical experiments, we fix . From Figures 35, we know that has certain influence on the solution, as increases, and the solution becomes smoother. In order to test the influence of on the option price, we set and let change at the same time. Figure 6 shows that the parameter has a significant effect on the price of options, and there will be an inflection point around . Finally, we investigated the influence of on the option price, and the results are shown in Figure 7, it can be seen that when increases, the option price also increases.

#### 6. Conclusion

In this paper, a new full-discrete numerical method is developed to solve the FBSM. An efficient -order and unconditionally energy stable method is constructed by combining the approach in time and Fourier method in space direction. It is proved that the full-discrete converges to the order globally. Numerical examples demonstrate the robustness and accuracy of the developed full-discrete method, numerically. Finally, we also study the properties of the solution of the FBSM.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

Juan He carried out an efficient numerical approach to time fractional Black–Scholes model. Aiqing Zhang helped to draft the manuscript. All authors read and approved the final manuscript.

#### Acknowledgments

The work of Juan He was supported by the China Scholarship Council (no. 202008520027). The work of Aiqing Zhang was supported by the Cultivation Project of Major Scientic Research Projects of Central University of Finance and Economics (no. 14ZZD007).

1. F. Black and M. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–654, 1973. View at: Publisher Site | Google Scholar
2. R. C. Merton, “Theory of rational option pricing,” The Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141–183, 1973. View at: Publisher Site | Google Scholar
3. R. Cioczek-Georges and B. B. Mandelbrot, “Alternative micropulses and fractional brownian motion,” Stochastic Processes and Their Applications, vol. 64, no. 2, pp. 143–152, 1996. View at: Publisher Site | Google Scholar
4. W. Min, B. P. English, G. Luo, B. J. Cherayil, S. C. Kou, and X. S. Xie, “Fluctuating enzymes: lessons from single-molecule studies,” Accounts of Chemical Research, vol. 38, no. 12, pp. 923–931, 2005. View at: Publisher Site | Google Scholar
5. T. A. M. Langlands, B. I. Henry, and S. L. Wearne, “Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions,” Journal of Mathematical Biology, vol. 59, no. 6, pp. 761–808, 2009. View at: Publisher Site | Google Scholar
6. 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, 2007. View at: Google Scholar
7. J. R. Liang, J. Wang, W. J. Zhang, W. Y. Qiu, and F. Y. Ren, “The solution to a bifractional black-scholes-merton differential equation,” International Journal of Pure and Applied Mathematics, vol. 58, no. 1, pp. 99–112, 2010. View at: Google Scholar
8. 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. View at: Publisher Site | Google Scholar
9. 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. View at: Publisher Site | Google Scholar
10. X. J. Yang, F. Gao, J. A. Machado, and D. Baleanu, “A new fractional derivative involving the normalized sinc function without singular kernel,” European Physical Journal Special Topics, vol. 226, no. 16-18, pp. 3567–3575, 2017. View at: Publisher Site | Google Scholar
11. X.-J. Yang, F. Gao, Y. Ju, and H.-W. Zhou, “Fundamental solutions of the general fractional-order diffusion equations,” Mathematical Methods in the Applied Sciences, vol. 41, no. 18, pp. 9312–9320, 2018. View at: Publisher Site | Google Scholar
12. M. Inc, A. Yusuf, A. I. Aliyu, and D. Baleanu, “Dark and singular optical solitons for the conformable space-time nonlinear Schrödinger equation with Kerr and power law nonlinearity,” Optik, vol. 162, pp. 65–75, 2018. View at: Publisher Site | Google Scholar
13. M. Inc, A. I. Aliyu, and A. Yusuf, “Dark optical, singular solitons and conservation laws to the nonlinear Schrödinger's equation with spatio-temporal dispersion,” Modern Physics Letters B, vol. 31, no. 14, p. 1750163, 2017. View at: Publisher Site | Google Scholar
14. J. Singh, D. Kumar, and D. Baleanu, “On the analysis of fractional diabetes model with exponential law,” Advances in Difference Equations, vol. 231, no. 1, 2018. View at: Google Scholar
15. J. Singh, D. Kumar, and D. Baleanu, “New aspects of fractional biswas-milovic model with mittag-leffler law,” Mathematical Modelling of Natural Phenomena, vol. 14, no. 3, p. 303, 2019. View at: Publisher Site | Google Scholar
16. J. Singh, D. Kumar, D. Baleanu, and S. Rathore, “On the local fractional wave equation in fractal strings,” Mathematical Methods in the Applied Sciences, vol. 42, no. 5, pp. 1588–1595, 2019. View at: Publisher Site | Google Scholar
17. K. B. Oldham and J. Spanier, “The fractional calculus,” Mathematical Gazette, vol. 56, no. 247, pp. 396–400, 1974. View at: Google Scholar
18. T. A. M. Langlands and B. I. Henry, “The accuracy and stability of an implicit solution method for the fractional diffusion equation,” Journal of Computational Physics, vol. 205, no. 2, pp. 719–736, 2005. View at: Publisher Site | Google Scholar
19. Z.-z. Sun and X. Wu, “A fully discrete difference scheme for a diffusion-wave system,” Applied Numerical Mathematics, vol. 56, no. 2, pp. 193–209, 2006. View at: Publisher Site | Google Scholar
20. Y. Lin and C. Xu, “Finite difference/spectral approximations for the time-fractional diffusion equation,” Journal of Computational Physics, vol. 225, no. 2, pp. 1533–1552, 2007. View at: Publisher Site | Google Scholar
21. R. H. De Staelena and A. S. Hendybc, “Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions,” Computers and Mathematics with Applications, vol. 74, no. 6, pp. 1166–1175, 2017. View at: Google Scholar
22. H. Zhang, F. Liu, I. Turner, and Q. Yang, “Numerical solution of the time fractional black-scholes model governing european options,” Computers & Mathematics with Applications, vol. 71, no. 9, pp. 1772–1783, 2016. View at: Publisher Site | Google Scholar
23. C. Li, T. Zhao, W. Deng, and Y. Wu, “Orthogonal spline collocation methods for the subdiffusion equation,” Journal of Computational and Applied Mathematics, vol. 255, pp. 517–528, 2014. View at: Publisher Site | Google Scholar
24. C.-M. Chen, F. Liu, I. Turner, and V. Anh, “A fourier method for the fractional diffusion equation describing sub-diffusion,” Journal of Computational Physics, vol. 227, no. 2, pp. 886–897, 2007. View at: Publisher Site | Google Scholar
25. A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer, Berlin, Germany, 1994.

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.