Matlab code for Lyapunov exponents of fractional order systems


International Journal of Bifurcation and Chaos, 28(5), 1850067 (2018) .


In this paper the Benettin-Wolf algorithm to determine all Lyapunov exponents for a class offractional-order systems modeled by Caputo's derivative and the corresponding Matlab codeare presented. First it is proved that the considered class of fractional-order systems admitsthe necessary variational system necessary to find the Lyapunov exponents. The underlying nu-merical method to solve the extended system of fractional order, composed of the initial valueproblem and the variational system, is the predictor-corrector Adams-Bashforth-Moulton forfractional differential equations. The Matlab program prints and plots the Lyapunov exponentsas function of time. Also, the programs to obtain Lyapunov exponents as function of the bifur-cation parameter and as function of the fractional order are described. The Matlab program forLyapunov exponents is developed from an existing Matlab program for Lyapunov exponents ofinteger order. To decrease the computing time, a fast Matlab program which implements theAdams-Bashforth-Moulton method, is utilized. Four representative examples are considered.

Add your rating and review

If all scientific publications that you have read were ranked according to their scientific quality and importance from 0% (worst) to 100% (best), where would you place this publication? Please rate by selecting a range.

0% - 100%

This publication ranks between % and % of publications that I have read in terms of scientific quality and importance.

Keep my rating and review anonymous
Show publicly that I gave the rating and I wrote the review