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. [1416] 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) 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 [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 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.

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.


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