Matlab code for Lyapunov exponents of fractional order

版本 3.0.3 (3.2 KB) 作者: Marius-F. Danca
Program to compute LEs of autonomous continuous systems of commensurate Caputo's Fractional Order
4.0
(2)
428.0 次下载
更新时间 2022/12/21

查看许可证

August 2022: plot modified to overcome the problems with plot function in the last matlab variants
***********************************
**slight changes**
The program, determines the Lyapunov Exponents (LEs) of a commensurate Fractional Order (FO) autonomous continuous-time system modeled by Capto's derivative [1] (for non-commensurate order see [2])
Because there are continuous citations of the paper presenting this program for Lyapunov Exponents (LEs) and due to some (minor) errors existing in the codes, here there are the slightly improved variants of related codes
Thus, the names of the several used codes are synchronized so that the use of programs is easier now
For commensurate order, all programs uses FDE12.m to integrate the FDEs:
https://goo.gl/XScYmi
The function containing the extended system (see e.g. LE_RF.m or LE_RF_p.m), has to be in the same folder with FDE12.m and FO_Lyapunov.m
Note that there are the following possibilities:
  1. The program, FO_Lyapunov.m, determines LEs of the (extended) system (see e.g. LE_RF.m) as function of time.
2.1 A modified variant of FO_Lyapunov (FO_Lyapunov_q.m) determines the LEs of a system (see e.g. LE_RF.m) for a given value of q.
2.2 FO_Lyapunov_q.m can be also used to obtain the variation of LEs as function of q. For this purpose it is used by run_FO_LE_q.m
3.1 Another modified variant FO_Lyapunov (FO_Lyapunov_p.m) determines the LEs of a system (see e.g. LE_RF_p.m) for a given value of p.
3.2 FO_Lyapunov_p.m can be used to obtain the variation of LEs as function of p. For this purpose it is used by run_LE_FO_p.m
As mentioned in [1], 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 (e.g. h=0.002 and h_norm=0.2, but to the detriment of computational time)
For non-commensurate order see FO_NC_Lyapunov.m which uses another program to integrate non-commensurate FDEs (see [2]) .
[1] 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
[2] Marius.-F. Danca, Matlab code for Lyapunov exponents of fractional-order systems, Part II: The non-commensurate case, International Journal of Bifurcation and Chaos, 31(12), 2150187, (2021)

引用格式

Marius-F. Danca (2024). Matlab code for Lyapunov exponents of fractional order (https://www.mathworks.com/matlabcentral/fileexchange/114595-matlab-code-for-lyapunov-exponents-of-fractional-order), MATLAB Central File Exchange. 检索来源 .

[1] 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 版本兼容性
创建方式 R2022a
兼容任何版本
平台兼容性
Windows macOS Linux
标签 添加标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

版本 已发布 发行说明
3.0.3

FDE12.m location
3.0.2

FDE12.m address
3.0.1

FDE12 link modified
3.0.0

plot function replaced
2.0.4

description
2.0.3

Description
2.0.2

Description
2.0.1

Citation
2.0.0