## Abstract

Discovering unknown aspects of non-equilibrium systems with hidden strange attractors is an attractive research topic. A novel quadratic hyperjerk system is introduced in this paper. It is noteworthy that this non-equilibrium system can generate hidden chaotic attractors. The essential properties of such systems are investigated by means of equilibrium points, phase portrait, bifurcation diagram, and Lyapunov exponents. In addition, a fractional-order differential equation of this new system is presented. Moreover, an electronic circuit is also designed and implemented to verify the feasibility of the theoretical model.

## 1 Introduction

After the emergence of chaos as a field of study, various chaotic systems have been reported in the literature. Most of them are 3D differential equations with a limited number of equilibria [1], [2], [3]. Therefore Shilnikov criteria can be used to prove the chaos in those systems [4], [5]. It has always been of interest to find chaotic systems with special features like multi-scroll attractors [6], [7], [8], [9], multistability [10], [11], [12], [13], and extreme multistability [14], [15], [16].

Leonov and Kuznetsov categorized dynamical systems into two groups: dynamical systems with self-excited attractors, and dynamical systems with hidden attractors [17], [18], [19], [20]. If the basin of attraction for an attractor involves (even in its border) at least one equilibrium, then that attractor is self-excited [21], [22], [23]; otherwise, it is hidden. Hidden attractors have been observed in many real dynamical systems [24], [25]. Design, localization, and control of hidden attractors have been of great interest in recent years [26], [27], [28], [29], [30].

In 1994, the first non-equilibrium chaotic flow was found [31]. After about 20 years, some other chaotic systems with non-equilibrium were introduced [32], [33], [34], [35], [36]. It can be easily concluded that the chaotic attractor in such systems is hidden [37], [38], [39], [40]. They have attracted the interest of the scientific community rapidly because chaotic systems with hidden attractors allow unexpected responses to perturbations in dynamical systems [41]. Hidden attractors also may be seen in 4D systems. As noted earlier, chaotic behaviours can emerge in (at least) 3D differential equations. The first 4D chaotic system was found by Rossler [42], which was the first step in designing a 4D chaotic system.

There is another common type of chaotic systems, which is jerk systems. It is the lowest derivative for which an ordinary differential equation (ODE) with smooth, continuous functions can give chaos [43]. A dynamical system that is described by an *n*^{th}-order ODE with *n*>3 is a hyper jerk system, which can be shown to have an explicit form *J* is a simple quadratic or cubic polynomial. Such systems are surprisingly general and prototypical examples of complex dynamical systems in a high-dimensional phase space [44], [45].

In this paper, a new hyper jerk 4D chaotic system with hidden attractors is introduced. This system is the first example of 4D hyper jerk systems with no equilibrium. In the next section, we present the dynamical analysis of this system. Bifurcation and Lyapunov exponent analysis are studied in Section 3. In Section 4, fractional-order and time-delayed versions of this new 4D system are proposed. To use this new chaotic system in engineering applications, in Section 5 we report a circuital implementation of the theoretical model. Finally, we conclude the paper in Section 6.

## 2 Dynamics of a New 4D Hyper Jerk Chaotic Flow

Consider the following new four-dimensional hyper jerk chaotic flow:

If *a*×*e*>0, this system has no equilibrium. Here *x*, *y*, *z*, and *w* are state variables, while *a*, *b*, *c*, *d*, and *e* are positive parameters. Interestingly, when *a*=2, *b*=5.05, *c*=2.5, *d*=0.7, and *e*=2 and the initial condition [*x*(0), *y*(0), *z*(0), *w*(0)]=(0, −5, 0, 0) is selected, the new system (1) has a chaotic solution. In this case, the Lyapunov exponents of the system (2) are *λ*_{1}=0.0772, *λ*_{2}=0, *λ*_{3}=−0.8044, and *λ*_{4}=−1.8293. Different projections of the chaotic attractor are illustrated in Figure 1.

## 3 Bifurcation Analysis

In this part, the dynamical behaviors of system (1) with respect to changing parameters *d* are investigated (see Fig. 2). In Figure 2a, the bifurcation diagram of the system is shown, and in Figure 2b the Lyapunov exponents are shown. It can be seen that changing the parameter *d* shows a familiar period doubling route to chaos. Figure 3 shows the basin of attraction for the hidden attractor. It is hidden because there is no equilibrium in this new chaotic hyper jerk system.

## 4 Fractional Order Hidden Jerk System (FOHS)

Recently, many researchers have discussed fractional-order calculus and its applications [46], [47]. Fractional-order nonlinear systems with different control approaches have been investigated in [48], [49], [50]. Fractional-order memristors based on non-equilibrium chaotic and hyper chaotic systems have been proposed in [35], [51]. A novel fractional-order non-equilibrium chaotic system has been investigated in [52], and a fractional-order hyper chaotic system without equilibrium points has been investigated in [53]. A memristor-based fractional-order system with a capacitor and an inductor is discussed in [54]. Numerical analysis and methods for simulating fractional-order nonlinear system are proposed in [55], and MATLAB solutions for fractional-order chaotic systems are discussed in [56].

Time-delayed differential equations play an important role in most engineering applications [57], [58]. Stability analysis of differential equations with delays has also been discussed in [59]. Synchronization of such time-delayed systems is a major complex problem, and therefore many synchronization schemes are discussed in the literature [57], [58]. A time-delayed chaotic system was obtained from a logistic chaotic map [60]. A parameter identification problem for a general time-delayed chaotic system is considered and analysed in [61]. A novel time-delayed chaotic system with hidden attractors is discussed in [62].

### 4.1 Designing Fractional-Order Hidden Jerk System (FOHS)

In this section we derive the fractional-order hidden jerk system (FOHS). There are three commonly used definition of the fractional-order differential operator, namely Grunwald-Letnikov, Riemann-Liouville, and Caputo [46], [47]. The dimensionless mathematical model of the FOHS can be derived as follows:

where *q _{x}*,

*q*,

_{y}*q*, and

_{z}*q*are the fractional orders of the FOHS system. The three main approaches derived to solve fractional-order chaotic systems are the frequency-domain method [63], the Adomian decomposition method (ADM) [64], and the Adams-Bashforth-Moulton (ABM) algorithm [65]. The frequency-domain method is not always reliable in detecting the chaotic behaviour in nonlinear systems [66]. On the other hand, ABM and ADM are more accurate and convenient to analyse dynamical behaviours of a nonlinear system. Compared to ABM, ADM yields more accurate results and needs less computing and memory resources [67]. Hence, the proposed FOHS system is implemeted in a field-programmable gate array (FPGA) by applying the ADM scheme. As the ADM algorithm converges fast [68], the first six terms are used to get the solution of the FOHS system as in real cases, and it is impossible to find an accurate value of

_{w}*x*when

*t*takes larger values [69]. Hence we have to design a time discretization method. That is to say, for a time interval

*t*(initial time) to

_{i}*t*(final time), we divide the interval into (

_{f}*t*,

_{n}*t*

_{n}_{+1}) and we get the value of

*x*(

*n*+1) at time

*t*

_{n}_{+1}by applying

*x*(

*n*) at time

*t*

_{n}_{+1}using the relation

*x*(

*n*+1)=

*F*(

*x*(

*n*)) [69].

We use the ADM method [64], [69] to numerically solve the FOHS system. The fractional-order discrete form of the dimensionless state equations for the FOHS system can be given as

where *i*=1, 2, 3, 4, and *x*^{2}, *xy*, *xz*, and *xw*. The Adomian polynomials for the linear terms can be calculated as

Nonlinear term | Adomian polynomials |
---|---|

x^{2} | |

xy | |

xz | |

xw |

In Table 1, *h*=*t _{n}*

_{+1}−

*t*

_{n}and Γ(∗) is the gamma function. More generally, the discrete Adomian form [31] of the FOHS can be presented as

Using (3) and (4), the FOHS is numerically solved. Figure 4 shows the 2D phase portraits of the FOHS with *h*=0.001, *q*=0.995, the parameters *a*=2, *b*=5.05, *c*=2.5, *d*=0.7, *e*=2, and the initial conditions [0, −5, 0, 0].

### 4.2 Dynamic Analysis of the Fractional-Order Hidden Jerk System

#### 4.2.1 Bifurcation with Fractional Order

Dynamical behaviors of the FOHJS are illustrated in Figure 5. We first derive the bifurcation plots of the FOHJS with parameters *b* and *e*, and the plots are shown in Figures 5a and b. The FOHJS shows chaotic oscillations for 5≤*b*≤5.1 and 1.5≤*e*≤2 with routine period doubling route to chaos. For *b*<5 and *e*>2, the system goes unbounded and loses chaotic oscillations. The next most important analysis of interest when investigating a fractional-order system is the bifurcation with fractional order. As can be seen from Figure 5c, bifurcation of the FOHS system for change in fractional order shows that the system oscillations remain chaotic if *q _{i}*>0.992 and the system shows hyper chaotic behaviour for 0.993≤

*q*≤0.998, and the largest positive Lyapunov exponents (

*L*

_{1}=0.3166,

*L*

_{2}=0.08217) of the FOHS system appear when

*q*=0.998, as can be seen from Figure 5d against its largest integer-order Lyapunov exponents (

*L*

_{1}=0.2991,

*L*

_{2}=0.07634). Figure 5e and f show the 2D phase portraits in the

*XY*-plane for various fractional orders.

#### 4.2.2 Time-Delayed Hidden Jerk System (TDHJS)

In this section we derive the time-delayed, dimensionless model of the hidden jerk system (HS) by introducing multiple time delays in the second state variable (*y*). TDHJS can be defined as

where *a*=2, *b*=5.05, *c*=2.5, *d*=0.7, and *e*=2. Similar to HS, the TDHJS shows chaotic oscillations for the initial conditions [0, −5, 0, 0], and the time delays are *τ*_{1}=0.05 and *τ*_{2}=0.04. Figure 6 shows the 2D phase portraits of TDHJS.

#### 4.2.3 Dynamic Properties of TDHJS

##### 4.2.3.1 Lyapunov Exponents

There are various algorithms proposed based on chaos synchronization for the estimation of the Lyapunov exponent of time-delayed dynamical systems [70], [71]. In this paper we adopt the technique employing the synchronization of identical systems coupled by a linear negative feedback mechanism [72] for finding the exact Lyapunov exponents of TDHJS. The calculated Lyapunov exponents are *L*_{1}=0.1684, *L*_{2}=0, *L*_{3}=−0.9974, and *L*_{4}=−2.078.

##### 4.2.3.2 Bifurcation

To understand the chaotic behaviour of TDHJS and the impact of time delays (*τ*_{1}, *τ*_{2}) and parameters (*a _{τ}*,

*b*,

_{τ}*c*,

_{τ}*d*,

_{τ}*e*), we derive the bifurcation plots and investigate the various regions of bifurcation. Figure 7 shows the bifurcation of TDHJS with time delays, and the time delays

_{τ}*τ*

_{1}and

*τ*

_{2}show chaotic oscillations for 0.03≤

*τ*

_{1}≤0.0785; for

*τ*

_{1}>0.785, the TDHJS system goes unbounded. For 0.001≤

*τ*

_{1}≤0.06, TDHJS shows chaotic oscillations, and for

*τ*

_{2}>0.06 the system spirals out and goes into periodic orbits.

In addition, the bifurcation diagrams of the TDHJS system with different parameters are reported in Figure 8. To investigate the bifurcation of the TDHJS system in the parameter space, the time delays are fixed at *τ*_{1}=0.05, *τ*_{2}=0.04. First we investigate the bifurcation of TDHJS with the parameter *a _{τ}*, with all the other parameter fixed to their respective values. TDHJS shows a chaotic region for 1.9≤

*a*≤2.8 and shows period-2 oscillations for

_{τ}*a*<1.7 and period-3 oscillations for

_{τ}*a*>2.85 with period-doubling route to chaos. Figure 8a shows the bifurcation of TDHJS with the parameter

_{τ}*a*and Figure 8b shows the bifurcation of TDHJS with the parameter

_{τ}*b*. The plots show two distinct chaotic regions for 5≤

_{τ}*b*≤5.09 and 5.1≤

_{τ}*b*≤5.2 separated by a symmetry-breaking and symmetry-restoring region for 5.09≤

_{τ}*b*≤5.1 with period-4 oscillations for

_{τ}*b*≥5.25 and enters and exits the chaotic regions by period-doubling and period-halving, respectively. The third bifurcation analysis is for the parameter

_{τ}*c*, as seen from Figure 8c. The bifurcation plot of

_{τ}*c*shows a chaotic region for 2.4≤

_{τ}*c*≤2.6, period-2 oscillations for

_{τ}*c*<2.25, and period-4 oscillations for 2.25≤

_{τ}*c*≤2.32 with the routine period doubling route to chaos. For

_{τ}*c*>2.6, the TDHJS system goes unbounded. Figure 8d shows the bifurcation of TDHJS with

_{τ}*d*, which shows a chaotic region for 0.55≤

_{τ}*d*≤0.875 with period-doubling and period-halving entry and exit to the chaotic region. The fifth bifurcation analysis is for the parameter

_{τ}*e*, which shows chaotic region for 0.8≤

_{τ}*e*≤3, as can be seen from Figure 8e.

_{τ}## 5 The Electronic Circuit Implementation of 4D New Chaotic System

In this section, the new chaotic system in (1) was implemented as an electronic circuit like in Figure 11. The electronic circuit of the new chaotic system was executed in the electronics lab. The circuit includes basic electronic materials like resistors, multipliers, and capacitors.

### 5.1 Scaling of the 4D New Chaotic System

To keep the amplitude values in the linear range of operational amplifiers, the chaotic system should be scaled. The *x*, *y*, *z*, and *w* signal values in the interval (−40, 60) are shown in Figure 9. The *z* and *w* values are not in the interval (−15, 15) for electronic circuit implementation. Thus, the *z* and *w* values need to be scaled for real-time applications. For the scaling process, let *x*, *y*, *z*, *w* instead of the variables *X*, *Y*, *Z*, *W*, the scaled system first becomes

According to new values (*X*, *Y*, *Z*, *W*), we regenerate of all the values as

With the new values, the scaled system is as shown below.

Finally, new scaled chaotic system is given by

In Figure 10, the new time series of the scaled system with *a*=2, *b*=5.05, *c*=2.5, *d*=0.7, and *e*=2 is shown. The amplitudes of *z* and *w* are decreased according to Figure 1. After these processes, we can execute the electronic circuit application.

### 5.2 Circuit Implementation of the New Scaled Chaotic System

A circuit is designed as shown in Figure 11 for the scaled chaotic system. The experimental electronic circuit of the scaled chaotic system was designed for parameters *a*=2, *b*=5.05, *c*=2.5, *d*=0.7, and *e*=2 with initial conditions *x*(0)=0, *y*(0)=−5, *z*(0)=0, and *w*(0)=0.

We select *C*_{1}=*C*_{2}=*C*_{3}=*C*_{4}=1 nF, *R*_{1}=400 kΩ, *R*_{2}=200 kΩ, *R*_{3}=*R*_{4}=*R*_{6}=*R*_{7}=*R*_{9}=*R*_{13}=*R*_{14}=100 kΩ, *R*_{5}=160 kΩ, *R*_{8}=15 MΩ, *R*_{10}=39.6 kΩ, *R*_{11}=40 kΩ, and *R*_{12}=57 kΩ; the corresponding phase portraits on the oscilloscope are shown in Figure 12.

## 6 Conclusion

A novel quadratic four-dimensional hyperjerk system with hidden chaotic attractors was proposed in this paper. Dynamical properties of this new system were investigated by means of equilibrium points, phase portrait, bifurcation diagram, and Lyapunov exponents through numerical simulation and circuit implementation. In addition, a fractional-order differential equation of this new system, along with analysis of its dynamical properties, was presented.

## References

[1] E. N. Lorenz, J. Atmos. Sci. **20**, 130 (1963).10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2Search in Google Scholar

[2] O. E. Rössler, Phy. Lett. A **57**, 397 (1976).10.1016/0375-9601(76)90101-8Search in Google Scholar

[3] G. Chen and T. Ueta, Int. J. Bifurcation Chaos **9**, 1465 (1999).10.1142/S0218127499001024Search in Google Scholar

[4] L. Shilnikov, Sov. Math. **6**, 163 (1965).Search in Google Scholar

[5] L. P. Shil’nikov, Methods of Qualitative Theory in Nonlinear Dynamics, World Scientific, Singapore, Vol. 5, 2001.10.1142/4221Search in Google Scholar

[6] J. Muñoz-Pacheco, E. Tlelo-Cuautle, I. Toxqui-Toxqui, C. Sánchez-López, and R. Trejo-Guerra, Int. J. Electron. **101**, 1559 (2014).10.1080/00207217.2014.880999Search in Google Scholar

[7] E. Tlelo-Cuautle, J. J. Rangel-Magdaleno, A. D. Pano-Azucena, P. J. Obeso-Rodelo, and J. C. Nunez-Perez, Commun. Nonlinear Sci. Numer. Simulat. **27**, 66 (2015).10.1016/j.cnsns.2015.03.003Search in Google Scholar

[8] F. R. Tahir, S. Jafari, V.-T. Pham, C. Volos, and X. Wang, Int. J. Bifurcation Chaos **25**, 1550056 (2015).10.1142/S021812741550056XSearch in Google Scholar

[9] S. Jafari, V.-T. Pham, and T. Kapitaniak, Int. J. Bifurcation Chaos **26**, 1650031 (2016).10.1142/S0218127416500310Search in Google Scholar

[10] Q. Lai and S. Chen, Int. J. Bifurcation Chaos **26**, 1650177 (2016).10.1142/S0218127416501777Search in Google Scholar

[11] J. Kengne, A. N. Negou, and D. Tchiotsop, Nonlinear Dynam. **87**, 2281 (2017).10.1007/s11071-016-3189-zSearch in Google Scholar

[12] P. Sharma, M. D. Shrimali, A. Prasad, N. V. Kuznetsov, and G. A. Leonov, Eur. Phys. J. Spec. Top. **224**, 1485 (2015).10.1140/epjst/e2015-02474-ySearch in Google Scholar

[13] J. C. Sprott, S. Jafari, A. J. M. Khalaf, and T. Kapitaniak, Eur. Phys. J. Spec. Top. **226**, 1979 (2017).10.1140/epjst/e2017-70037-1Search in Google Scholar

[14] B. Bao, T. Jiang, Q. Xu, M. Chen, H. Wu, et al., Nonlinear Dynam. **86**, 1711 (2016).10.1007/s11071-016-2988-6Search in Google Scholar

[15] B.-C. Bao, Q. Xu, H. Bao, and M. Chen, Electron. Lett. **52**, 1008 (2016).10.1049/el.2016.0563Search in Google Scholar

[16] B. Bao, H. Bao, N. Wang, M. Chen, and Q. Xu, Chaos, Soliton. Fract. **94**, 102 (2017).10.1016/j.chaos.2016.11.016Search in Google Scholar

[17] G. A. Leonov and N. V. Kuznetsov, Int. J. Bifurcation Chaos **23**, 1330002 (2013).10.1142/S0218127413300024Search in Google Scholar

[18] G. A. Leonov, N. V. Kuznetsov, and T. N. Mokaev, Eur. Phys. J. Spec. Top. **224**, 1421 (2015).10.1140/epjst/e2015-02470-3Search in Google Scholar

[19] Z. Wei, M. Chen, B. Bao, and X. Quan, Int. J. Bifurcation Chaos **26**, 1650125 (2016).10.1142/S021812741650125XSearch in Google Scholar

[20] Z. Wei, I. Moroz, J. C. Sprott, A. Akgul, and W. Zhang, Chaos: Interdis. J. Nonlinear Sci. **27**, 033101 (2017).10.1063/1.4977417Search in Google Scholar
PubMed

[21] D. Dudkowski, O. A. Kuznetsova, G. A. Leonov, T. N. Mokaev, and N. V. Stankevich, Phys. Rep. **637**, 1 (2016).10.1016/j.physrep.2016.05.002Search in Google Scholar

[22] Z. Wei, P. Yu, W. Zhang, and M. Yao, Nonlinear Dynam. **82**, 131 (2015).10.1007/s11071-015-2144-8Search in Google Scholar

[23] Z. Wei, W. Zhang, Z. Wang, and M. Yao, Int. J. Bifurcation Chaos **25**, 1550028 (2015).10.1142/S0218127415500285Search in Google Scholar

[24] G. Leonov, N. Kuznetsov, and V. Vagaitsev, Phys. D: Nonlinear Phenomena **241**, 1482 (2012).10.1016/j.physd.2012.05.016Search in Google Scholar

[25] G. Leonov, N. V. Kuznetsov, M. A. Kiseleva, E. P. Solovyeva, and A. M. Zaretskiy, Nonlinear Dynam. **77**, 277 (2014).10.1007/s11071-014-1292-6Search in Google Scholar

[26] P. R. Sharma, M. D. Shrimali, A. Prasad, N. V. Kuznetsov, and G. A. Leonov, Int. J. Bifurcation Chaos **25**, 1550061 (2015).10.1142/S0218127415500613Search in Google Scholar

[27] U. Chaudhuri and A. Prasad, Phys. Lett. A **378**, 713 (2014).10.1016/j.physleta.2014.01.003Search in Google Scholar

[28] D. Dudkowski, A. Prasad, and T. Kapitaniak, Phys. Lett. A **379**, 2591 (2015).10.1016/j.physleta.2015.06.002Search in Google Scholar

[29] F. Nazarimehr, S. Jafari, S. Mohammad Reza Hashemi Golpayegani, and J. C. Sprott, Int. J. Bifurcation Chaos **27**, 1750023 (2017).10.1142/S0218127417500237Search in Google Scholar

[30] F. Nazarimehr, B. Saedi, S. Jafari, and J. C. Sprott, Int. J. Bifurcation Chaos **28**, 1750037 (2017).10.1142/S0218127417500377Search in Google Scholar

[31] J. C. Sprott, Phys. Rev. E **50**, R647 (1994).10.1103/PhysRevE.50.R647Search in Google Scholar

[32] Z. Wei, Phys. Lett. A **376**, 102 (2011).10.1016/j.physleta.2011.10.040Search in Google Scholar

[33] S. Jafari, J. Sprott, and S. M. R. H. Golpayegani, Phys. Lett. A **377**, 699 (2013).10.1016/j.physleta.2013.01.009Search in Google Scholar

[34] V.-T. Pham, C. Volos, S. Jafari, and T. Kapitaniak, Nonlinear Dynam. **87**, 2001 (2017).10.1007/s11071-016-3170-xSearch in Google Scholar

[35] K. Rajagopal, A. Karthikeyan, and P. Duraisamy, Chaos, Soliton. Fract. **103**, 476 (2017).10.1016/j.chaos.2017.07.007Search in Google Scholar

[36] K. Rajagopal, S. Jafari, and G. Laarem, Pramana **89**, 92 (2017).10.1007/s12043-017-1487-8Search in Google Scholar

[37] V.-T. Pham, A. Akgul, C. Volos, S. Jafari, and T. Kapitaniak, AEU – Int. J. Electron. Commun. **78**, 134 (2017).10.1016/j.aeue.2017.05.034Search in Google Scholar

[38] V.-T. Pham, S. Jafari, C. Volos, T. Gotthans, X. Wang, et al., Optik – Int. J. Light Electron Opt. **130**, 365 (2017).10.1016/j.ijleo.2016.10.100Search in Google Scholar

[39] V.-T. Pham, S. Takougang, K. Christos, V. S. Jafari, and T. Kapitaniak, AEU – Int. J. Electron. Commun. **78**, 220 (2017).10.1016/j.aeue.2017.04.012Search in Google Scholar

[40] V.-T. Pham, X. Wang, S. Jafari, C. Volos, and T. Kapitaniak, Int. J. Bifurcation Chaos **27**, 1750097 (2017).10.1142/S0218127417500973Search in Google Scholar

[41] S. Jafari and J. Sprott, Chaos, Soliton. Fract. **57**, 79 (2013).10.1016/j.chaos.2013.08.018Search in Google Scholar

[42] O. Rossler, Phys. Lett. A **71**, 155 (1979).10.1016/0375-9601(79)90150-6Search in Google Scholar

[43] J. C. Sprott, Elegant Chaos: Algebraically Simple Chaotic Flows, World Scientific, Singapore 2010.10.1142/7183Search in Google Scholar

[44] J. L. Sprott, Int. J. Chaos Theory Appl. **5**, 1 (2000).Search in Google Scholar

[45] K. E. Chlouverakis and J. Sprott, Chaos, Soliton. Fract. **28**, 739 (2006).10.1016/j.chaos.2005.08.019Search in Google Scholar

[46] V. Lakshmikantham and A. Vatsala, Nonlinear Anal.: Theor. Methods Appl. **69**, 2677 (2008).10.1016/j.na.2007.08.042Search in Google Scholar

[47] K. Diethelm, The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type, Springer, Germany, 2010.10.1007/978-3-642-14574-2Search in Google Scholar

[48] M. P. Aghababa, J. Comput. Nonlinear Dynam. **7**, 021010 (2012).10.1115/1.4005323Search in Google Scholar

[49] E. A. Boroujeni and H. R. Momeni, Signal Process. **92**, 2365 (2012).10.1016/j.sigpro.2012.02.009Search in Google Scholar

[50] R. Zhang and J. Gong, Syst. Sci. Control Eng.: Open Access J. **2**, 751 (2014).10.1080/21642583.2014.891955Search in Google Scholar

[51] K. Rajagopal, A. Karthikeyan, and A. K. Srinivasan, Nonlinear Dynam, **87**, 2281 (2017).10.1007/s11071-016-3189-zSearch in Google Scholar

[52] R. Li and W. Chen, Chinese Phys. B **22**, 040503 (2013).10.1088/1674-1056/22/4/040503Search in Google Scholar

[53] D. Cafagna and G. Grassi, Chinese Phys. B **24**, 080502 (2015).10.1088/1674-1056/24/8/080502Search in Google Scholar

[54] M.-F. Danca, W. K. Tang, and G. Chen, Chaos, Soliton. Fract. **84**, 31 (2016).10.1016/j.chaos.2015.12.018Search in Google Scholar

[55] I. Petráš, Acta Montanistica Slovaca **11**, 273 (2006).Search in Google Scholar

[56] Z. Trzaska, in: Engineering Education and Research Using MATLAB (Ed. A. H. Assi), Intech, Rijeka, Chapter 19, 2011.Search in Google Scholar

[57] Y. Xiao, W. X. Sufang, and T. XiuchunLi, App. Math. Comput. **213**, 538 (2009).10.1016/j.amc.2009.03.049Search in Google Scholar

[58] W. He, F. Qian, J. Cao, and Q.-L. Han, Phys. Lett. A **375**, 498 (2011).10.1016/j.physleta.2010.11.052Search in Google Scholar

[59] W. Deng, Y. Wu, and C. Li, Int. J. Bifurcation Chaos **16**, 465 (2006).10.1142/S0218127406014939Search in Google Scholar

[60] L. Acho, A continuous-time delay chaotic system obtained from a chaotic logistic map, IASTED International Conference Modelling, Identification and Control. “Modelling, Identification and Control (MIC 2017)”, ACTA Press, Innsbruck, 2017, p. 147.10.2316/P.2017.848-001Search in Google Scholar

[61] C. Qin-Qin, Acta Phys. Sin. **64**, 240506 (2015).10.7498/aps.64.240506Search in Google Scholar

[62] V.-T. Pham, S. Vaidyanathan, C. K. Volos, S. Jafari, N. V. Kuznetsov, et al., Eur. Phys. J. Spec. Top. **225**, 127 (2016).10.1140/epjst/e2016-02625-8Search in Google Scholar

[63] A. Charef, H. H. Sun, Y. Y. Tsao, and B. Onaral, IEEE Trans. Automat. Contr. **37**, 1465 (1992).10.1109/9.159595Search in Google Scholar

[64] G. Adomian, Math. Computer Model. **13**, 17 (1990).10.1016/0895-7177(90)90125-7Search in Google Scholar

[65] H. Sun, A. Abdelwahab, and B. Onaral, IEEE Trans. Automat. Contr. **29**, 441 (1984).10.1109/TAC.1984.1103551Search in Google Scholar

[66] M. Tavazoei and M. Haeri, IET Signal Process. **1**, 171 (2007).10.1049/iet-spr:20070053Search in Google Scholar

[67] H. Shao-Bo, S. Ke-Hui, and W. Hui-Hai, Acta Phys. **63**, 030502 (2014).10.7498/aps.63.030502Search in Google Scholar

[68] R. Caponetto and S. Fazzino, Int. J. Bifurcation Chaos **23**, 1350050 (2013).10.1142/S0218127413500508Search in Google Scholar

[69] S. He, K. Sun, and H. Wang, Entropy **17**, 8299 (2015).10.3390/e17127882Search in Google Scholar

[70] A. Stefański, T. Kapitaniak, and A. Dąbrowski, in: IUTAM Symposium on Chaotic Dynamics and Control of Systems and Processes in Mechanics (Eds. G. Rega and F. Vestroni), Springer, Netherlands, 2005.Search in Google Scholar

[71] A. Dabrowski, Nonlinear Dynam. **78**, 1601 (2014).10.1007/s11071-014-1542-7Search in Google Scholar

[72] A. Stefanski, A. Dabrowski, and T. Kapitaniak, Chaos, Soliton. Fract. **23**, 1651 (2005).10.1016/S0960-0779(04)00428-XSearch in Google Scholar

**Received:**2017-11-12

**Accepted:**2017-12-28

**Published Online:**2018-02-15

**Published in Print:**2018-02-23

©2018 Walter de Gruyter GmbH, Berlin/Boston