Abstract:The asymptotic stable region and long-time decay rate of solutions to linear homogeneous Caputo time fractional ordinary differential equations (F-ODEs) are known to be completely determined by the eigenvalues of the coefficient matrix. Very different from the exponential decay of solutions to classical ODEs, solutions of F-ODEs decay only polynomially, leading to the so-called Mittag-Leffler stability, which was already extended to semi-linear F-ODEs with small perturbations. This work is mainly devoted to the qualitative analysis of the long-time behavior of numerical solutions. By applying the singularity analysis of generating functions developed by Flajolet and Odlyzko (SIAM J. Disc. Math. 3 (1990), 216-240), we are able to prove that both $\mathcal{L}$1 scheme and strong $A$-stable fractional linear multistep methods (F-LMMs) can preserve the numerical Mittag-Leffler stability for linear homogeneous F-ODEs exactly as in the continuous case. Through an improved estimate of the discrete fractional resolvent operator, we show that strong $A$-stable F-LMMs are also Mittag-Leffler stable for semi-linear F-ODEs under small perturbations. For the numerical schemes based on $\alpha$-difference approximation to Caputo derivative, we establish the Mittag-Leffler stability for semi-linear problems by making use of properties of the Poisson transformation and the decay rate of the continuous fractional resolvent operator. Numerical experiments are presented for several typical time fractional evolutional equations, including time fractional sub-diffusion equations, fractional linear system and semi-linear F-ODEs. All the numerical results exhibit the typical long-time polynomial decay rate, which is fully consistent with our theoretical predictions.

Polynomial decay of a linear system of PDEs via Caputo fractional‐time derivative

The Stability and Stabilization of Infinite Dimensional Caputo-Time Fractional Differential Linear Systems

Solutions of systems with the Caputo-Fabrizio fractional delta derivative on time scales

Mittag--Leffler stability of numerical solutions to time fractional ODEs

Finite‐time stability and numerical approximations of fractional neutral delay systems involving proportional Caputo derivative

Existence and Stability for Fractional Differential Equations with a ψ–Hilfer Fractional Derivative in the Caputo Sense

Analysis of Caputo–Katugampola fractional differential system

Computational and Numerical Analysis of the Caputo-Type Fractional Nonlinear Dynamical Systems via Novel Transform

Efficient results on fractional Langevin-Sturm-Liouville problem via generalized Caputo-Atangana-Baleanu derivatives

Stability and Logarithmic Decay of the Solution to Hadamard-Type Fractional Differential Equation

Linearized Asymptotic Stability for Fractional Differential Equations

Energy stability of variable-step L1-type schemes for time-fractional Cahn-Hilliard model

Theory on Linear L-Fractional Differential Equations and a New Mittag–Leffler-Type Function

A generalized definition of Caputo derivatives and its application to fractional ODEs

A new approach in handling one-dimensional time-fractional Schrödinger equations

Generalized Proportional Caputo Fractional Differential Equations with Delay and Practical Stability by the Razumikhin Method

Numerical Study of Caputo Fractional-Order Differential Equations by Developing New Operational Matrices of Vieta–Lucas Polynomials

Theory on linear L-fractional differential equations and a new Mittag-Leffler-type function

Stability analysis of a second-order difference scheme for the time-fractional mixed sub-diffusion and diffusion-wave equation

Existence and Stability Results for Differential Equations with a Variable-Order Generalized Proportional Caputo Fractional Derivative

Stability analysis of fractional order time‐delay systems: constructing new Lyapunov functions from those of integer order counterparts