A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications

  • Authors:

  • Guang-Hua Gao;Zhi-Zhong Sun;Hong-Wei Zhang

  • Affiliations:

  • College of Science, Nanjing University of Posts and Telecommunications, Nanjing 210023, PR China and Department of Mathematics, Southeast University, Nanjing 210096, PR China;Department of Mathematics, Southeast University, Nanjing 210096, PR China;School of Mathematical Sciences, Tianjin Normal University, Tianjin 300387, PR China

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2014

Quantified Score

Hi-index 31.45

Visualization

Abstract

In the present work, first, a new fractional numerical differentiation formula (called the L1-2 formula) to approximate the Caputo fractional derivative of order @a (0=2), while the linear interpolation approximation is applied on the first small interval [t"0,t"1]. As a result, the new formula can be formally viewed as a modification of the classical L1 formula, which is obtained by the piecewise linear approximation for f(t). Both the computational efficiency and numerical accuracy of the new formula are superior to that of the L1 formula. The coefficients and truncation errors of this formula are discussed in detail. Two test examples show the numerical accuracy of L1-2 formula. Second, by the new formula, two improved finite difference schemes with high order accuracy in time for solving the time-fractional sub-diffusion equations on a bounded spatial domain and on an unbounded spatial domain are constructed, respectively. In addition, the application of the new formula into solving fractional ordinary differential equations is also presented. Several numerical examples are computed. The comparison with the corresponding results of finite difference methods by the L1 formula demonstrates that the new L1-2 formula is much more effective and more accurate than the L1 formula when solving time-fractional differential equations numerically.