Version 3.0.0 (3.77 KB) by Marius-F. Danca
Program to compute LEs as function of q, of autonomous systems of commensurate Caputo's Fractional Order
Updated 28 Aug 2022
August 2022: plot modified to overcome the problems with plot function in the last matlab variants
The program, FO_Lyapunov_q, can be used either alone, to determined the LEs of a FO system for a fixed fractional order q (see e.g. LE_RF.m which contains the extended system), or can be used to obtain the variation of LEs as function of q, case when the code run_FO_LE_q must be used .
- To obtain the LEs for a given q, for, e.g., RF system, one uses
For example, for the RF system 
2. If one intends to obtain the evolution of LEs as function of q one uses
E.g., for the same system RF
Note that FO_Lyapunov_q.m, LE_RF.m, run_FO_LE_q and FDE12.m (necessary to integrate the system) must be in the same folder.
As mentioned in , the relation between h_norm and h is essential. Here both are chosen equal (0.02), but multiple of h for h_norm should be tried
 Marius-F. Danca and N. Kuznetsov, Matlab code for Lyapunov exponents of fractional order systems, International Journal of Bifurcation and Chaos, 28(05)(2018), 1850067
Marius-F. Danca (2023). FO_Lyapunov_q (https://www.mathworks.com/matlabcentral/fileexchange/114605-fo_lyapunov_q), MATLAB Central File Exchange. Retrieved .
Marius-F. Danca and N. Kuznetsov, Matlab code for Lyapunov exponents of fractional order systems, International Journal of Bifurcation and Chaos, 28(05)(2018), 1850067
MATLAB Release Compatibility
Created with R2022a
Compatible with any release
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.