Approximations for the Caputo Derivative (I)
For researchers solving fractional differential equations, this work provides a more accurate numerical method for approximating the Caputo derivative, though the improvement is incremental over existing L1-type schemes.
This paper constructs new approximations for the Caputo derivative of various orders, achieving higher accuracy by modifying initial and last weights. The approximation of order 2-α with weights k^{-1-α}/Γ(-α) yields numerical solutions for fractional differential equations that are more accurate than those using the L1 approximation in all experiments.
In this paper we construct approximations for the Caputo derivative of order $1-α,2-α,2$ and $3-α$. The approximations have weights $0.5\left((k+1)^{-α}-(k-1)^{-α}\right)/Γ(1-α)$ and $k^{-1-α}/Γ(-α)$, and the higher accuracy is achieved by modifying the initial and last weights using the expansion formulas for the left and right endpoints. The approximations are applied for computing the numerical solution of ordinary fractional differential equations. The properties of the weights of the approximations of order $2-α$ are similar to the properties of the $L1$ approximation. In all experiments presented in the paper the accuracy of the numerical solutions using the approximation of order $2-α$ which has weights $k^{-1-α}/Γ(-α)$ is higher than the accuracy of the numerical solutions using the $L1$ approximation for the Caputo derivative.