Matlab Code for Lyapunov Exponents of Fractional-Order Systems, Part II: The Noncommensurate Case.

In this paper, the Benettin–Wolf algorithm for determining all Lyapunov exponents of noncommensurate 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 systems, the Adams–Bashforth–Moulton scheme for fractional differential equations is utilized. Like 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 a function of orders, or a function of a bifurcation parameter. Special attention is paid to the periodicity of fractional-order systems and its influences. The case of noncommensurate Lorenz system is demonstrated. [ABSTRACT FROM AUTHOR]