Research Article  Open Access
Yaoyu Wang, Ling Liu, Xinshan Cai, Chongxin Liu, Yan Wang, Guangchao Zheng, "FractionalOrder Hidden Attractor Based on the Extended Liu System", Mathematical Problems in Engineering, vol. 2020, Article ID 1418272, 22 pages, 2020. https://doi.org/10.1155/2020/1418272
FractionalOrder Hidden Attractor Based on the Extended Liu System
Abstract
In this paper, a new commensurate fractionalorder chaotic oscillator is presented. The mathematical model with a weak feedback term, which is named hypogenetic flow, is proposed based on the Liu system. And with changing the parameters of the system, the hidden attractor can have no equilibrium points or line equilibrium. What is more interesting is that under the occasion that no equilibrium point can be obtained, the phase trajectory can converge to a minimal field under the lead of some initial conditions, similar to the fixed point. We call it the virtual equilibrium point. On the other hand, when the value of parameters can produce an infinite number of equilibrium points, the line equilibrium points are nonhyperbolic. Moreover than that, there are coexistence attractors, which can present hyperchaos, chaos, period, and virtual equilibrium point. The dynamic characteristics of the system are analyzed, and the parameter estimation is also studied. Then, an electronic circuit implementation of the system is built, which shows the feasibility of the system. At last, for the fractional system with hidden attractors, the finitetime synchronization control of the system is carried out based on the finitetime stability theory of the fractional system. And the effectiveness of the controller is verified by numerical simulation.
1. Introduction
In the past few decades, chaos has been established as an important branch of modern physics and mathematics due to its objectivity and universality as well as its wide application, for example, in cryptography [1], economics [2], electronic communication [3], and even in the popular field of neural network [4, 5] in recent years. As a pioneer and basic of research, the famous Lorenz system was proposed in 1963 [6]. Based on this, Chen and Ueta system [7] and Lü and Chen system [8] were proposed one after another, which are similar to the Lorenz system but have different topological structure. In 2004, Liu et al. put forward a Liu system which is also a butterfly attractor in [9]. Later, its circuit experimental confirmation was introduced in [10]. With the development of theory and technology, the research has been extended to fractional differential system. Facts show that fractional calculus with a history of more than 300 years is very helpful for interdisciplinary research and has applications in many aspects, such as physics, finance, and sociology. Naturally, the study of the fractional chaotic system has also been reported a lot, such as chaos control and synchronization [11] and fractional memristor [12]. In 2006, Lu and Liu analyzed the dynamics of the fractional Liu system and realized the circuit based on the approximation theory of the fractional operator [13]. In [14], the coupling synchronization of two similar fractional Liu systems is discussed. It is also based on the fractional Liu system, and Hegazi et al. studied the existence and uniqueness of solutions and the influence of the fractional order on chaos control and realized the chaos synchronization of the commensurate fractional Liu system and its integer order counterpart by function projection synchronization in [15]. In [16], a fourthorder commensurate fractional hyperchaotic system is proposed by adding a feedback variable to the Liu system, but with 11 terms including three quadratic nonlinearities. In this paper, we propose a new fourthorder commensurate fractional system based on an extension of the threeorder Liu system, which is simpler than the system in [16]. It is a particular system with many interesting features that are not available in other systems that have been similarly proposed: (1) it has nine terms, including two quadratic nonlinearities, a weak feedback term, and a constant. (2) It is a hidden hyperchaotic attractor with no equilibrium points. (3) It has small neighborhoods of the origin, and orbit from which is all attracted by the origin, the virtual equilibrium point, despite the fact that its steadystate equations have no real solutions. (4) Although it is fourdimensional, the rank of its Jacobian matrix is less than 4. That is, each row or column of the Jacobian matrix has at least one element that is zero. (5) It is expressed as quaternion fractional differential equations containing a signum function of additional state variable, through which amplitude and polarity of state variables can be transferred to each other, except for the additional state variable. Because the additional state is fed back to another in the form of signum function, whose return is an integer variable, indicating only the sign of the variable. If not all informationcontaining amplitudes and polarities can be transmitted between variables effectively, the hypogenetic chaotic is defined as a weak feedback chaotic flow with incomplete information transfer [17]. (6) With the constant term selected conditionally, the hidden attractor has nonhyperbolic line equilibrium on the coordinate axis, not no equilibrium points.
The rest of the paper is organized as follows: Section 2 describes some preliminaries of fractional calculus. In Section 3, the new commensurate fractional system is introduced and the dynamical behaviors are numerically analyzed. In Section 4, approximation of the fractionalorder integrator and circuit implementation of the new system is discussed. In Section 5, the DSP implementation is introduced. The chaos synchronization controller is designed for the fractionalorder system and numerical simulations are carried out in Section 6. Finally, some concluding remarks are drawn in Section 7.
2. Preliminaries of Fractional Calculus
Definition 1. Fractional integral of order α for the function can be expressed as follows:where Γ (.) is the Gamma function, a and t are the bounds of the operation, and α, a ∈ R, α > 0. There are many definitions of fractional calculus, and two widely used definitions are the Riemann–Liouville definition and Caputo definition. For the reason that the initial values of the Caputo fractionalorder differential equations are clear, it is more common in engineering applications.
Definition 2. The Caputo definition of fractional derivatives can be written aswhere n − 1 < α < n, and the C represents the Caputo definition.
3. System Model and Dynamics Analysis
In this section, we first present the Liu system given by Liu et al. [9]:
In which, x, y, and z ∈ R are the states and a, b, c, k, and h ∈ R are parameters. By simulation, it can be easily come out that system (3) is chaotic with a = 10, b = 40, c = 2.5, k = 1, and h = 4 [9]. Hence, a 4D commensurate fractionalorder system is generated from system (3) as follows:with the same values of (a, b, c, k, and h), and settle = 5 and and d = 20. When m = 20, the equilibrium point of system (4) is in the negative half of the axis, and the hidden attractor is shown in Figure 1 with final time 3600 and a time step of 0.005.
(a)
(b)
(c)
(d)
(e)
Comparing with the original, a weak feedback term about is the feedback to the second equation of the system in the form of a symbolic function. Then, the steadystate equations can be written as
If m = −20, the equilibrium point of system (4) is in the positive half of the axis. Then, a system with line equilibrium can be obtained. The line equilibrium is located on the axis except for the origin. To estimate the Jacobian matrix, the partial derivatives of the state variables are required. For calculation, the sign () is replaced by the tanh (nw), where n is a large positive number. It can be seen from Figure 2 that the larger the n is, the more similar the graph of the function tanh (nw) is to the graph of the signum function.
Via linearizing, the Jacobian matrix can be gotten as
And the eigenvalues are calculated as
for
Hence, regardless of the value of m, there is an eigenvalue whose real part is zero. A nonhyperbolic equilibrium point has one or more eigenvalues with a zero real part. Therefore, the system has a nonhyperbolic line equilibrium which is uncommon in chaotic systems with hidden attractors.
On the other hand, when m ≠ 20, there are three results in the steadystate equations x_{e} = y_{e} = z_{e} = 0. However, there is no result for , which signifies that system (4) has no equilibrium. And the hidden attractor can be seen in Figure 3 with m = 1 and q = 0.99.
(a)
(b)
(c)
(d)
(e)
When there is no equilibrium, the 3D Poincare mapping is also plotted in Figure 4, which also illustrates the chaotic characteristics of the system (4). The plane we chose is x + = 0. Then, in order to better analyze the trend of the trajectories, we mark them with different colors according to the different directions of the trajectories. If phase trajectories of system (4) cross the plane in a positive direction (x + > 0), the trace will be painted in blue. If phase trajectories of system (4) cross the plane negatively (x + < 0), the trace will be painted in orange. The scattered points shown in Figure 4 are consistent with the characteristics of chaos.
3.1. The Influence of Parameters on the System
Proverbially, bifurcation diagram and the LE are two dynamical tools to investigate the dynamics behaviors of nonlinear systems. Through investigating the bifurcation with respect to changing parameters, the dynamic behaviors of the system can be explored. While the LE not only can exhibit the projections of the system but also can display the changes for each state variable.
The calculated values of the LEs of system (4) considering various sampling and observation time with m = 1 and m = 20 are given in Table 1 and 2, respectively. Here, the hyperbolic tangent function is still utilized and the Jacobian matrix used in the calculation is shown in equation (6).


It is clear from Table 1 and Table 2 that under different parameter values, observation time, and sampling time, the hidden attractor has a (+, +, 0, −) sign of Lyapunov exponent demonstrating a 4D hyperchaotic system [18]. Moreover, the Lyapunov exponents for m = 1 and m = 20 are shown in Figures 5 and 6, respectively. It is obvious from Figures 5 and 6 that system (4) with has (+, +, 0, −) nature of the Lyapunov exponents and thus confirm hyperchaotic nature of system (4).
(a)
(b)
(a)
(b)
With the change of the parameter m, the bifurcation and LE of system (4) are shown in Figures 7(a) and 7(b). System (4) displays abundant dynamic behavior when changing the value of the bifurcation parameter m from 45 to 70. As can be seen in Figure 7(b), there is a reverse perioddoubling route to chaos, when increasing the value of the parameter m from 45 to 53.5, except for 47.5 < m < 48.25, and the chaotic behavior is shown. In addition, when the parameter m > 53.5, system (4) exhibits a series of evident perioddoubling bifurcations. When the control parameter m is increased above the value of 68.75, the attractor enters periodic behavior eventually.
(a)
(b)
For fractional system models, the behavior of the system varies greatly at different orders. The bifurcation with q is shown in Figure 8, from which chaos and periodic behavior can be clearly observed. And the chaotic behavior can still be observed when q = 0.835, illustrating that the chaotic characteristics of the proposed system are significant. Besides, when in the section of (0.877, 0.887) and (0.937, 0.96), the system is periodic. However, when q > 0.975, chaos appeared again. These characteristics are confirmed again by the phase plane plot. In Figure 9 with changing the q, the period (as Figure 9(a)), perioddoubling (as Figure 9(b)), and chaotic (as Figure 9(c)) trajectories are shown, which are consistent with the conclusion and analysis in Figure 8.
(a)
(b)
(a)
(b)
(c)
(d)
(e)
(f)
Furthermore, the largest Lyapunov exponents (LEmax) for two parameters are calculated in Figure 10. In Figure 10(a), the maximum of the LEmax can exceed 2, which indicates the chaotic system. However, in some parameter values, the LEmax can be less than 0, and there may be periodic trajectories or a stable case. Figure 10(b) is about parameters c and h. The LEmax shows obvious symmetry about the value of h. In Figure 10, as the LEmax increases, the color changes from dark blue to yellow. Hence, for the selection of parameters, the parameters in the dark blue region corresponding to the LEmax equal to or less than 0 are more likely to lead to limit cycle or stability. In the light green region corresponding to the positive LEmax, it can be considered as weak chaos. The yellow region corresponding to a large positive LEmax denotes chaotic motion. Therefore, choosing the right parameters is crucial for studying the behavior of different states of the system.
(a)
(b)
3.2. Hidden Attractor Coexistence
In this occasion, the interesting feature, the virtual equilibrium point, will be analyzed. As we all know that nonlinear dynamical systems basically have a variety of attractive characteristic behaviors, and so does system (4). Under the same parameter variable value, different initial values can lead the trajectory to different dynamic behaviors. For some initial values, the phase trajectory can be hyperchaos. For some special initial values, the phase trajectory can converge to a very small region, similar to the fixed point, which we call the virtual equilibrium point. Of interesting is the fact that although it is a hidden attractor without any equilibrium points, the system still can converge to the virtual equilibrium point.
By means of numerical simulation, the time series of the trajectory starting from different initial values can be obtained, through which we can intuitively know that whether the system is chaotic or convergent. Previously analyzed, the trajectories from most initial conditions will eventually lead to a chaotic path. And only the trajectories from some initial values in a particular region will converge to the virtual equilibrium point. So two representative initial values are selected under different circumstances. With m = 1 and q = 0.96 and the others parameters are the same with above, the time series are shown in Figure 11, in which the blue track exhibits chaos with the initial conditions X (0) = (2.2, 2.4, 0.8, 0) and the track in orange was from X (0) = (0, 0, 0.3, 0.1) tending to the virtual equilibrium point.
(a)
(b)
(c)
(d)
Besides, the trajectories for the system are plotted in Figure 12. Under the same parameters, the system exhibits chaos behavior from X (0) = (2.2, 2.4, 0.8, 0) (in blue), while the system converges to a minimal region around (0, 0, 0, 0) corresponding to the initial conditions X (0) = (0, 0, 0.3, 0.1) (marked with the pentacle in orange). Figures 12(a) and 12(b) are the phase plots without transient process. And Figures 12(c) and 12(d) exhibit the transient trajectories from different initial values. It is obvious that trajectories with (0, 0, 0.3, 0.1) as initial conditions eventually converged to the minimum area near the origin which is consistent with Figure 11.
(a)
(b)
(c)
(d)
3.3. Parameter Estimation
In this section, the parameter estimations will be studied employing the new cost function for parameter estimation of chaotic models, which was proposed by Sajad Jafari in 2014 [19]. The main idea of this method is to get the regression map of the real system and the estimated system by the time series of the real system and the estimated system. Then, according to the minimum distance between the regression map of the estimation system and the regression map of the real system, the cost function is calculated to judge the estimation effect of the unknown parameters, so as to realize the estimation of the unknown parameters of the chaotic system. The new cost function can effectively overcome the sensitivity of initial conditions of the chaotic system in time domain, so it can complete the estimation of parameters under the condition that the initial conditions are not very accurate.
System (4) is the actual system to be estimated, with the parameters a = 10, b = 40, c = 2.5, = 5, h = 4, d = 20, and m = 16 and initial conditions (2.2, 2.4, 0.8, 0). We make parameters c and h both unknown for estimate, and the simulation result of the cost function with parameters c and h is obtained as shown in Figure 13.
(a)
(b)
When the value of the cost function is the minimum, the corresponding parameter value is the estimate (where the white square is located). As can be seen in Figure 13(b), the parameter estimation results are very accurate. Different colors can generally reflect the value of the cost function and dynamical features of the system. Commonly, the parameter values of cold color region corresponding to a low cost function are unlikely to cause chaos oscillation of the system, while the parameter values of warm color region corresponding to a high cost function are more likely to cause the display of chaos of the system.
4. Circuit Implementation of the FractionalOrder System
Circuit designs have been used to confirm the realization ability of mathematical models [20–24]. Based on PSIM software, a circuit realization scheme of the commensurate fractional chaotic system with order q = 0.99 and parameters a = 10, b = 40, c = 2.5, k = 1, h = 4, = 5, d = 20, and m = 1 is designed and implemented.
The realization of the fractionalorder integrator is the most important for the circuit of the fractional chaotic system. With its zeros z_{i} and poles p_{i}, the transfer function 1/sq is transformed into its zeropolegain form as shown in equation (7). With the transfer function approximation method proposed in [25] and ω_{max} = 100, P_{T} = 0.01 assumed, the zeros z_{i} and poles p_{i} are calculated as equation (8) and the expression of the fractional integrator is obtained as equation (9) under the approximate error of ε = 0.3 dB:
It can be realized by three chain fractals composed of three resistorcapacitor pairs as shown in Figure 14, and its transfer function is derived as equation (10). Comparing equations (9) and (10), these parameters of electronic components are calculated as follows: Ca = 1.011 μF, Cb = 13.982 μF, Cc = 13.050 μF, Ra = 95081.958 kΩ, Rb = 6.441 kΩ, and Rc = 6.436 Ω:
Based on the operationalamplifier approach [26], the circuit of the commensurate fractional chaotic system is shown in Figure 15. To make sure to observe the hidden attractor completely, we reduce the output of chaotic signal to 5% of the original and set
Hence, we have 20k = 20, d/20 = 1, m/20 = 0.8, and 20h = 80. With Kirchhoff’s circuit laws, the circuital equations corresponding to system (4) are described as follows:where C_{0} = 1 μF.
The parameters of electronic components are calculated as follows: R1 = R3 = R5 = R7 = R12 = R14 = R15 = 10kΩ, R2 = R6 = 20kΩ, R4 = R10 = R13 = 100kΩ, R8 = R9 = R16 = 40kΩ, R11 = 5kΩ, R17 = 200kΩ, R18 = 400kΩ, and R19 = 500kΩ. And Figure 14 demonstrates the circuit implementation results.
Comparing with the numerical simulations in Figure 3, the hidden attractors observed in Figure 16 shrunk by around 95%. Intuitively, the axis range of z in Figure 16 is (0, 6), which is 5% of the axis range of z in Figure 3 (0, 120). And the axis ranges of x, y, and changed in the same proportion, indicating that the circuital results are comparatively consistent with the numerical simulations.
(a)
(b)
(c)
(d)
5. DSP Implementation of the FractionalOrder System
For fractionalorder chaotic systems, the design of the fractionalorder integrator is crucial to its realization on the DSP platform. At present, the Adomian decomposition method (ADM for short) which is proposed in [27] is most frequently used to solve the fractionalorder chaotic system [28, 29].
Firstly, we separate the fractionalorder system into two parts:where m = ceil (q), b_{k} is a specified constant relating to the initial values, Lx (t) and Nx (t) are the linear and nonlinear terms of the fractional differential equations, respectively, and is the constant of the system. And the nonlinear terms can be decomposed bywhere i = 0, 1, …, ∞. Then, the nonlinear terms are expressed as
Thus, the solution of equation (13) can be described as
The method to get each decomposition item is derived by
According to equations (13)–(17), the discrete solution of the fractionalorder chaotic system can be denoted as x (n + 1) = F (x (n)). Because of its high accuracy, fast convergence, and less computer resource consumption, the ADM becomes one of the best candidate methods to decompose the nonlinear differential equations [30].
The fractionalorder chaotic system is
By employing the ADM and setting step size h = t_{n} _{+} _{1} − t_{n}, the discrete iterative formulas of proposed system (18) can be obtained. Considering the fast convergence performance of the ADM, the first six terms of the iteration are selected:
According to the given conditions, the following equation can be obtained:where Γ (·) is the Gamma function.
Make , and the second term of each state variable is shown as follows:where
Similarly, the component x^{2} iswhere
By applying the same approach to the components x^{3} and x^{4}, the component x^{5} iswherewherewhere
According to equations (21)–(32), the discrete iterative equations can be obtained which are necessary for the fractionalorder system implemented on the DSP platform. For hardware design, the block diagram of the working principle is shown in Figure 17.
In the experiment, the YXDSPF28335 CORE is employed for its powerful computing ability, and the key chip in this board is the floatpoint DSP TMS320F28335. A lowpower 4 channel and 12bit parallel D/A converter DAC7724 is used to convert the time series generated by DSP. Then, the converted data are sent to an oscilloscope (Agilent DSOX 3052A) which is used to capture the phase portraits of the system. The flow diagram for DSP implementation is shown in Figure 18.
Firstly, the DSP is initialized, and initial value settings and some preparation before calculation are made. Then, the system is iterated according to equations (19)–(30) and the result is pushed in. In the result processing the data is converted into analog signals and sent to an oscilloscope which is used to record phase portraits. After that, the result is pop out from the DSP and replaces the initial value for the next iteration. Iterations continue until the phrase portrait is complete.
Through the DSP platform, the experimental results were obtained as shown in Figure 19. It is consistent with the numerical simulations and the circuit implementation results.
(a)
(b)
(c)
6. FiniteTime Synchronization of the FractionalOrder System
6.1. FiniteTime Stability Theory of FractionalOrder System
The synchronization of fractional chaotic systems has great potential in the field of secure communication. With the development of the research, many synchronization methods have been proposed. In [31], a theory of finitetime stability of the fractionalorder system is proposed, which can realize the stability of approximation error in finite time, with fast synchronization speed and strong robustness. In this paper, the finitetime stability theory of the fractionalorder system is used to realize system synchronization control.
Theorem 1. (see[31]). Consider a fractionalorder system satisfyingwhere 0 < q < 1, and the states and . Then, the state variable of the fractionalorder system tends to zero in finitetime t, and the system is stable:where .
Lemma 1. When satisfying a, b < 0 and 0 < c < 1, we can get the inequality as follows:Let system (4) be the master system, where q = 0.99, then its slave system is as follows:where are the controllers based on Theorem 1.
Then, let be the errors between master system (4) and slave system (35). Thus, get the error systemThe controllers can be designed aswhere are the parameters related to the controllers.
Theorem 2. Adding the controller of equation (36), the fractionalorder error system can be stable in a finitetime.
Proof. Based on error system (35) and designed controllers (36), the approximation error can be gotten:According to equation (31), we can getFor the inequality in Lemma 1, thenHence, we can getwhere the control parameters are set as λ = 1.6 and β = 0.8. Therefore, according to equation (31),satisfying Theorem 1. Hence, the masterslave system can realize finitetime synchronization and the error system can also be stable in a finitetime.
6.2. Simulation Results
With q = 0.99, the time step t_{Δ} = 0.001, and final time T_{sim} = 8s, the initial conditions of the master system (as equation (4)) and its slave system (as equation (34)) are (2.2, 2.4, 0.8, 0.1) and (3.5, 3.4, 1.4, −0.5), respectively. Under the control parameters λ = 1.6 and β = 0.8, the simulation result of the chaos synchronization is shown in Figure 20.
If with the initial conditions as (2.2, 2.4, 0.8, 0.1) and (3, 2, 1.2, 0.5), and the control parameters λ = 1.6, β = 0.8, the simulation result of the chaos synchronization is shown in Figure 21.
In Figure 20, the synchronization errors e_{1}, e_{2}, e_{3}, and e_{4} can convergence to zero in less than 2.5 seconds with the controllers. And the synchronization errors can convergence to zero in less than 2 seconds in Figure 21, for the smaller errors between the initial conditions. Both the results of the two conditions indicate that the state variables of the masterslave system can synchronize in a finite time weather the initial conditions for have different signs. Therefore, the designed controller is effective.
7. Conclusion and Perspectives
This paper has introduced a fractionalorder chaotic oscillator based on the Liu system. Changing the parameters, the hidden attractor can have no equilibrium point, in which situation the virtual equilibrium point was studied, and have line equilibrium point. It is worth noting that this infinite number of equilibrium points is not the usual hyperbolic equilibrium point, but the uncommon nonhyperbolic line equilibrium point. In addition, through threedimensional Poincare, bifurcation, and Lyapunov exponential spectrum and the parameter estimation of the system, the influence of the parameters are analyzed. Besides, the feasibility of the new 4D nonlinear system is confirmed by the help of electronic circuit implementation by means of PSIM and DSP. Finally, based on the theory of fractional finitetime stability, a finitetime synchronization controller is designed to realize the synchronization control of the fractional chaotic system with hidden attractors. And the validity of the method is verified by numerical simulations. Fractional hidden attractors have important research significance and application value in chaos and cryptography. In the future, we will further explore the application of this research in image encryption technology and secure communication.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
There are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (grant nos. 51877162 and 51977173).
References
 Z. Hua, F. Jin, B. Xu, and H. Huang, “2D LogisticSinecoupling map for image encryption,” Signal Processing, vol. 149, pp. 148–161, 2018. View at: Publisher Site  Google Scholar
 J. Ma, T. Xu, and Y. Hong, “Impact research on a nonlinear cold chain evolutionary game under three various contracts,” International Journal of Bifurcation and Chaos, vol. 29, no. 15, Article ID 1950058, 2019. View at: Publisher Site  Google Scholar
 J. Petrzela, “Fractionalorder chaotic memory with wideband constant phase elements,” Entropy, vol. 22, no. 4, p. 422, 2020. View at: Google Scholar
 Z. Zeng and W. Zheng, “Multistability of neural networks with timevarying delays and concaveconvex characteristic,” IEEE Transactions on Neural Networks and Learning Systems, vol. 23, no. 2, pp. 293–305, 2012. View at: Google Scholar
 Z. Zeng, T. Huang, and W. Zheng, “Multistability of recurrent networks with timevarying delays and the piecewise linear activation function,” IEEE Transactions on Neural Networks, vol. 21, no. 8, pp. 1371–1377, 2010. View at: Google Scholar
 E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, vol. 20, no. 2, pp. 130–141, 1963. View at: Publisher Site  Google Scholar
 G. Chen and T. Ueta, “Yet another chaotic attractor,” International Journal of Bifurcation and Chaos, vol. 9, no. 7, pp. 14651466, 1999. View at: Publisher Site  Google Scholar
 J. Lü and G. Chen, “A new chaotic attractor coined,” International Journal of Bifurcation and Chaos, vol. 12, no. 3, pp. 659–661, 2002. View at: Publisher Site  Google Scholar
 C. Liu, T. Liu, L. Liu, and K. Liu, “A new chaotic attractor,” Chaos, Solitons & Fractals, vol. 22, no. 5, pp. 1031–1038, 2004. View at: Publisher Site  Google Scholar
 F. Wang and C. Liu, “Studies on Liu chaotic system and its experimental confirmation,” Acta Physics Sinica, vol. 55, no. 10, pp. 5061–5069, 2006. View at: Google Scholar
 K. Li, J. Cao, and J. He, “Hidden hyperchaotic attractors in a new 4D fractional order system and its synchronization,” Chaos, vol. 30, no. 3, Article ID 033129, 2020. View at: Publisher Site  Google Scholar
 Y. Yu, M. Shi, H. Kang, M. Chen, and B. Bao, “Hidden dynamics in a fractionalorder memristive HindmarshRose model,” Nonlinear Dynamics, vol. 100, no. 1, pp. 891–906, 2020. View at: Publisher Site  Google Scholar
 J. Lu and C. Liu, “Realization of fractionalorder Liu chaotic system by circuit,” Chinese Physics, vol. 16, no. 6, pp. 1586–1590, 2007. View at: Google Scholar
 X. Wang and M. Wang, “Dynamic analysis of the fractionalorder Liu system and its synchronization,” Chaos, vol. 17, no. 3, Article ID 033106, 2007. View at: Google Scholar
 A. S. Hegazi, E. Ahmed, and A. E. Matouk, “On chaos control and synchronization of the commensurate fractional order Liu system,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 5, pp. 1193–1202, 2013. View at: Publisher Site  Google Scholar
 Q. Han, C. Liu, L. Sun, and D. Zhu, “A fractional order hyperchaotic system derived from a Liu system and its circuit realization,” Chinese Physics, vol. 22, no. 2, Article ID 020502, 2013. View at: Publisher Site  Google Scholar
 C. Li, J. Sprott, and H. Xing, “Hypogenetic chaotic jerk flows,” Physics Letters A, vol. 380, no. 1112, pp. 1172–1177, 2016. View at: Publisher Site  Google Scholar
 J. P. Singh and B. K. Roy, “The nature of Lyapunov exponents is (+, +, −, −). Is it a hyperchaotic system?” Chaos, Solitons & Fractals, vol. 92, pp. 73–85, 2016. View at: Publisher Site  Google Scholar
 S. Jafari, J. Sprott, V. Pham, S. Golpayegani, and A. Jafari, “A new cost function for parameter estimation of chaotic systems using return maps as fingerprints,” International Journal of Bifurcation and Chaos, vol. 24, no. 10, Article ID 1450134, 2014. View at: Google Scholar
 W. Xue, Y. Li, S. Cang, H. Jia, and Z. Wang, “Chaotic behavior and circuit implementation of a fractionalorder permanent magnet synchronous motor model,” Journal of the Franklin Institute, vol. 352, no. 7, pp. 2887–2898, 2015. View at: Publisher Site  Google Scholar
 Z. Wang, H. Tang, and Z. Chen, “The design and implementation of a multiwing chaotic attractor based on a fiveterm threedimension system,” International Journal of Circuit Theory and Applications, vol. 44, no. 5, pp. 1186–1201, 2016. View at: Publisher Site  Google Scholar
 S. Wang, S. Lao, H. Chen, J. Chen, and S. Chen, “Implementation of the fractionalorder ChenLee system by electronic circuit,” International Journal of Bifurcation and Chaos, vol. 23, no. 2, Article ID 1350030, 2013. View at: Google Scholar
 S. Zhang, Y. Zeng, Z. Li, and C. Zhou, “Hidden extreme multistability, antimonotonicity and offset boosting control in a novel fractionalorder hyperchaotic system wthout equilibrium,” International Journal of Bifurcation and Chaos, vol. 28, no. 13, Article ID 1850167, 2018. View at: Google Scholar
 M. Wang, X. Liao, Y. Deng, Z. Li, Y. Zeng, and M. Ma, “Bursting, dynamics, and circuit implementation of a new fractionalorder chaotic system with coexisting hidden attractors,” Journal of Computational and Nonlinear Dynamics, vol. 14, no. 7, Article ID 071002, 2019. View at: Publisher Site  Google Scholar
 A. Charef, H. H. Sun, Y. Y. Tsao, and B. Onaral, “Fractal System as represented by singularity function,” IEEE Transactions on Automatic Control, vol. 37, no. 9, pp. 1465–1470, 1992. View at: Publisher Site  Google Scholar
 A. Akgul, H. Calgan, I. Koyuncu, I. Pehlivan, and A. Istanbullu, “Chaosbased engineering applications with a 3D chaotic system without equilibrium points,” Nonlinear Dynamics, vol. 84, no. 2, pp. 481–495, 2016. View at: Publisher Site  Google Scholar
 G. Adomian, “A new approach to nonlinear partial differential equations,” Journal of Mathematical Analysis and Applications, vol. 102, no. 2, pp. 420–434, 1984. View at: Publisher Site  Google Scholar
 D. Cafagna and G. Grassi, “Bifurcation and chaos in the fractionalorder Chen system via a timedomain approach,” International Journal of Bifurcation and Chaos, vol. 18, no. 7, pp. 1845–1863, 2018. View at: Google Scholar
 N. Shawagfeh, “Analytical approximate solutions for nonlinear fractional differential equations,” Applied Mathematics and Computation, vol. 131, no. 23, pp. 517–529, 2002. View at: Publisher Site  Google Scholar
 S. He, K. Sun, and S. Banerjee, “Dynamical properties and complexity in fractionalorder diffusionless Lorenz system,” European Physical Journal Plus, vol. 131, no. 8, pp. 254–2016. View at: Google Scholar
 L. Zhao, J. Hu, and X. Liu, “Adaptive tracking control and synchronization of fractional hyperchaotic Lorenz system with unknown parameters,” Acta Physics Sinica, vol. 59, no. 4, pp. 2305–2309, 2010. View at: Google Scholar
Copyright
Copyright © 2020 Yaoyu Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.