Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 1393456, 9 pages
https://doi.org/10.1155/2020/1393456

Research Article

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

Juan He1,2 and Aiqing Zhang¹

¹Business School, Central University of Finance and Economics, Beijing 100081, China
²Accounting School, Guizhou University of Finance and Economics, Guiyang 550025, Guizhou, China

Correspondence should be addressed to Aiqing Zhang; zhangaiqing@cufe.edu.cn

Received 4 May 2020; Accepted 27 July 2020; Published 4 September 2020

Academic Editor: Hamdy Nabih Agiza

Copyright © 2020 Juan He and Aiqing Zhang. 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.

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 [1], which is based on the assumption that stocks and options are in an “ideal state” in the market and Samuelson’s model [2]:

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 [3]. 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 [6] and Liang et al. [7] considered fractional Brownian motion in Samuelson’s model equation:

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

Jumarie [8] 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. [9] 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. [14–16] solved a series of fractional equations by using homotopy analysis technique and Laplace transform algorithm.

As we all know, as an effective method, formula [17] has been widely used in the calculation of fractional differential equations. Langlands and Henry [18], Sun and Wu [19], and Lin and Xu [20] discussed the error estimate for scheme. De Staelena and Hendybc [21] constructed a numerical method of fourth-order finite difference in space and in time. Stability, uniqueness, and error estimates are analyzed. Zhang et al. [22] 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) as

where , 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 obtain

It is easy to verify that the following formula is correct:

Thus,

Giving up some positive terms, we have

Assume the following inequality holds:

Next, we will show is still valid. If , taking the inner product of (15) with , we derive

Note the fact that

Thus, we get

This 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 [25]:

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 that

Therefore, 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) as

Subtracting (38) from (10) at , we note that

Then, we have

Set , we have

Dropping some positive terms, we find

Assume

Next, we will prove that it holds also for . Subtracting (28) from a reformulation of (10) at , we find

Let , we have

Thus, we have

Note that , thus

Note that

This 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:

where

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

Table 1: Temporal convergence orders of various time steps for Example 1.

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

Figure 1: Error in time direction for different .

Figure 2: The and errors in space direction with .

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 3–5, 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.

Figure 3: The dynamic behavior of solution to FBSM equation at .

Figure 4: The dynamic behavior of solution to FBSM equation at .

Figure 5: The dynamic behavior of solution to FBSM equation at .

Figure 6: The influence of various on the option price.

Figure 7: The influence of various on the option price.

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

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,” The Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141–183, 1973.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

K. B. Oldham and J. Spanier, “The fractional calculus,” Mathematical Gazette, vol. 56, no. 247, pp. 396–400, 1974.

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.

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.

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.

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.

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.

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.

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.

A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer, Berlin, Germany, 1994.