Numerical simulation of fractional-order Duffing system with extended Mittag-Leffler derivatives

In this paper, we studied the dynamics of a nonlinear fractional-order Duffing system combined with Mittag-Leffler derivatives in order to provide dynamic behaviors different from existing ones. The Mittag-Leffler derivative is a generalized version of the exponential kernel derivative. To achieve this goal, we introduced a modified extension to higher-order Mittag-Leffler derivatives to overcome the initialization problem. Moreover, we discussed some properties and relationships of the studied derivatives. Then we presented numerical schemes to handle fractional extensions of the considered oscillatory system including the Mittag-Leffler and the Caputo derivatives. Numerical simulations are carried out and the resulting simulation dynamics of the studied fractional oscillatory system are compared.