Abstract

We introduced a parameter which was related to ; then two numerical schemes for variable-order Caputo fractional derivatives were derived; the second-order numerical approximation to variable-order fractional derivatives and -order approximation for are established. For the given parameter , the error estimations of formulas were proven, which were higher than some recently derived schemes. Finally, some numerical examples with exact solutions were studied to demonstrate the theoretical analysis and verify the efficiency of the proposed methods.

1. Introduction

Fractional differential equations include constant-order and variable-order equations; a great quantity of natural phenomena can be modeled by variable-order fractional differential equations; the study of such problems has attracted much attention. In recent years, profound background of physical applications for the variable-order fractional calculus has been already established; the definition of variable-order operator has been investigated in [14]. Numerous problems in mathematical physics and engineering have been modeled by variable-order fractional differential equations, such as successful applications in mechanics [5], in the simulation of linear and nonlinear viscoelasticity oscillators [6], and in other cases where the order of the derivative varies with time [7]. In addition, a physical experimental investigation of variable-order operators has been considered in [8]. Variable-order fractional derivatives can be used to model anomalous diffusion, as they can describe the time dependent diffusion process more specifically than fractional derivatives of constant order, just as shown in [9]. The difference schemes of fractional derivatives with constant and variable order are investigated in [10].

Due to the existence of variable fractional derivative, it is usually difficult to obtain the analytical solution of such equations; therefore, it is particularly significant to give numerical solutions to these problems. Cao and Qiu [11] derived a high order numerical method for variable-order fractional ordinary differential equation by establishing a second-order numerical approximation to variable-order Riemann-Liouville fractional derivative. Fu et al. [12] adopted the method of approximate particular solutions for both constant-order and variable-order time fractional diffusion models. Several finite difference methods for variable-order fractional partial diffusion equations were proposed in [1318].

As what we can see above, there are many different definitions of the variable-order fractional derivatives. However, different from other definitions in mathematics, only initial condition is needed for the variable-order Caputo definition which can be easily used in physical field; this definition means that the memory rate of system is determined by the current time instant and changes with time. For simplicity, only variable-order Caputo definition was discussed in the whole paper.

A good approximation of the variable-order Caputo derivatives was observed in [19], it was defined at the points and ; then what will happen if it is redefined in a neighborhood at the point ? Approximation formula for the -order Caputo derivatives at the point was investigated in [20], where is a superconvergence point which ensures that the convergence order improves from the local second order to the overall order. In the current work, we applied the ideas in [20]; a parameter was given, which was changed with , and the selection of played a vital role in the numerical results in the paper. We proposed two new approximation formulas of second-order and accuracy for variable-order time fractional operator with orders and , respectively. Specifically, we adopted the following definition of variable-order Caputo fractional derivatives:

The paper is organized as follows. In the next section, we present two new formulas, select for and for , respectively, and provide the corresponding error analysis. In Section 3, the numerical verifications are presented; as can be seen from the numerical results, minor changes to the selected parameters will have a significant impact on the error estimates. Finally, in Section 4, the conclusion is drawn.

2. Approximation Formulas

In this section, we present the following lemma, which plays a vital role in the later analysis.

Lemma 1 (see [19]). Assuming that the derivatives of the function exist to the order of on and to the order of on is the th degree interpolation polynomial of based on the points , then, for , the following holds: where depends on and , and
For a given function and an integer , denote , where is the step. Introduce the following notation: and, for convenience, denote .

2.1. Second-Order Formula for

Denote

Estimating the derivative with order at the grid point , from (1), we directly obtain

Evaluate the integration on each subinterval, leading to

For each interval , denote the second-degree interpolation polynomial in the Lagrange form as follows: and for the last interval , denote the first-degree interpolation polynomial in the Lagrange form as

As an approximation formula, the following result is obtained: where where

Substituting (8) and (10) into (5) yields and the truncation error is

Here we denote the discrete approximation formula for the variable-order derivative with order as ; thus from (12), it is easy to get the following result: where

Next, the analysis for the approximation error of formula (12) was given in detail.

Theorem 2. Let ; the following holds: where is a positive constant.

Proof. From (12), obviously Using Lemma 1, it produces Next, the error on each interval is analyzed.
When , from (19), it follows thatWhen , from (18), this leads to the following estimation: where is a positive constant.
From (19), it can be obtained in the same method. Substituting (20)–(22) into (17) leads to the theorem.

2.2. Second-Order Formula for

Denote , . Computing the variable-order derivative with order at ,

For the interval , the cubic interpolation polynomial in the Hermite form is as and it follows thatthe truncation error is

For the interval , the cubic interpolation polynomial at points , and in the Lagrange form is as leading to the truncation error is using the cubic interpolation polynomial at points , and in the Lagrange form to approximate on the last interval , we get the following holds: the truncation error is

Substituting (25), (28), and (31) into (23), there is

Denote the approximation formula of the variable-order derivative with order as and obtain where

Substituting the above formula into (34), thus

For the first term in the right hand side of (36), this leads to

For the second term in the right hand side of (36), the following holds:

For the third term in the right hand side of (36), this yields

Substituting (37)–(39) into (36), we obtain the approximate difference scheme of the : where

Similar to the proof of the Theorem 2, the truncation error of formula (40) is given by following theorem.

Theorem 3. For , the following holds: where is a positive constant.

Proof. From (33), we obtain the truncation error From Lemma 1, the result below is natural. Substituting (45) into (44), it yields the estimation where is a positive constant.
Substituting (46) into (44) leads to the theorem.

Remark 4. From the second term of the right hand of the truncation error estimate in Theorem 2, we can easily obtain the fact that formula (14) does not have exact second-order accuracy; it is related to the selection of and . This can be seen from Table 2 in the numerical example below; its accuracy is much higher than the second order.

3. Numerical Verification

In this section, the validity and numerical accuracy of the new presented formula (14) and (40) are demonstrated, respectively. Meanwhile, the corresponding computational results with the formula of [19] are given for contrast.

Take a positive integer , let , and denote

Example 1 (accuracy of formula (14) and [19]). Take . Compute the Caputo fractional derivative of at numerically.
The exact solution is given by From the results presented in Table 1, taking , we found that although both formula (14) and [19] have second-order approximation accuracy, the numerical accuracy by the formula (14) is significantly higher than that by the formula of [19]; furthermore, the computational errors are also obviously smaller than the formula of [19]. Thus, the new formula (14) is valid for solving complex problems that require high accuracy, and better computational results can be obtained from the formula.
In Table 2, taking , it shows that the convergence order is higher as compared with the theoretical result, and the error is smaller than the theoretical result, but why this happens is not clear yet; it requires further investigation.
Moreover, on the basis of Table 1, if a slight disturbance to is added to Tables 3 and 4, the experimental results will be greatly different. It can be easily found that there is a very close relationship between the convergence order and the selection of . The experimental results show that the convergence order is the ideal result only when . Tables 2, 5, and 6 also illustrate this fact.

Example 2 (accuracy of formula (40)). Take . Compute the Caputo fractional derivative of at numerically.
The exact solution is given by Table 7 lists the maximum errors and convergence rates for the variable-order . In particular, when , , and , the numerical results demonstrate that formula (40) leads to accuracy.
Similar to Example 1, we also give Tables 8 and 9 as a comparison; the results show that is the best choice.

Remark 3. Due to the fact that the selection of is very important, it is worth noting that different selected leads to the different numerical accuracy. From considerable trials, we observed that only if we take for and for , will numerical results be closer to theoretical analysis. In addition, carefully observing Table 9, although the convergence order is not higher than in Table 7, the error is lower than in Table 7; further investigation is needed on this phenomenon.

4. Conclusion

In this paper, two approximations to variable-order Caputo fractional derivatives were developed, and the analysis for the truncation errors of new formulas was made. In addition, numerical examples support theoretical results. In further work, we also want to give a strict proof of the selection of and apply these new methods to solve more complex cases within a reasonable accuracy.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by Science and Technology Department Key Projects of Fujian Province, China (Grant 2014H0034), and the Natural Science Foundation of Fujian Province, China (Grants 2017J01402, 2017J01557, and 2016J01310).