In this paper, a numerical algorithm to solve Caputo differential equations is proposed. The proposed algorithm utilizes the R2 algorithm for fractional integration based on the fact that the Caputo derivative of a function f(t) is defined as the Riemann–Liouville integral of the derivative . The discretized equations are integer order differential equations, in which the contribution of from the past behaves as a time-dependent inhomogeneous term. Therefore, numerical techniques for integer order differential equations can be used to solve these equations. The accuracy of this algorithm is examined by solving linear and nonlinear Caputo differential equations. When large time-steps are necessary to solve fractional differential equations, the high-speed algorithm (HSA) proposed by the present authors (Fukunaga, M., and Shimizu, N., 2013, “A High Speed Algorithm for Computation of Fractional Differentiation and Integration,” Philos. Trans. R. Soc., A, 371(1990), p. 20120152) is employed to reduce the computing time. The introduction of this algorithm does not degrade the accuracy of numerical solutions, if the parameters of HSA are appropriately chosen. Furthermore, it reduces the truncation errors in calculating fractional derivatives by the conventional trapezoidal rule. Thus, the proposed algorithm for Caputo differential equations together with the HSA enables fractional differential equations to be solved with high accuracy and high speed.
A Numerical Method for Caputo Differential Equations and Application of High-Speed Algorithm
TSK Bld. 302, 9-1, Aza-Sunada,
Iwaki, Fukushima 971-8135, Japan
Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received August 15, 2018; final manuscript received May 4, 2019; published online July 17, 2019. Assoc. Editor: Dumitru Baleanu.
- Views Icon Views
- Share Icon Share
- Search Site
Fukunaga, M., and Shimizu, N. (July 17, 2019). "A Numerical Method for Caputo Differential Equations and Application of High-Speed Algorithm." ASME. J. Comput. Nonlinear Dynam. September 2019; 14(9): 091007. https://doi.org/10.1115/1.4043794
Download citation file: