Matlab code for Lyapunov exponents of fractional-order systems, Part II: The non-commensurate case

Marius-F. Danca
DOI: https://doi.org/10.48550/arXiv.2104.11323
2021-04-08
Computational Physics
Abstract:In this paper the Benettin-Wolf algorithm to determine all Lyapunov exponents adapted to a class of non-commensurate fractional-order systems modeled by Caputo's derivative and the corresponding Matlab code are presented. The paper continues the work started in [Danca & Kuznetsov, 2018], where the Matlab code of commensurate fractional-order systems is given. To integrate the extended system, the Adams-Bashforth-Moulton for fractional differential equations is utilized. As the the Matlab program for commensurate-order systems, the program presented in this paper prints and plots all Lyapunov exponents as function of time. The program can be simply adapted to plot the evolution of the Lyapunov exponents as function of orders, or a function of a bifurcation parameter. A special attention is payed to the periodicity of fractional-order systems and its influences. The case of the Lorenz system is considered.
What problem does this paper attempt to address?