# Hidden Coexisting Attractors in a Fractional-Order System without Equilibrium: Analysis, Circuit Implementation, and Finite-Time Synchronization

Authors: Guangchao Zheng; Ling Liu; Chongxin Liu
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In this paper, a novel three-dimensional fractional-order chaotic system without equilibrium, which can present symmetric hidden coexisting chaotic attractors, is proposed. Dynamical characteristics of the fractional-order system are analyzed fully through numerical simulations, mainly including finite-time local Lyapunov exponents, bifurcation diagram, and the basins of attraction. In particular, the system can generate diverse coexisting attractors varying with different orders, which presents ample and complex dynamic characteristics. And there is great potential for secure communication. Then electronic circuit of the fractional-order system is designed to help verify its effectiveness. What is more, taking the disturbances into account, a finite-time synchronization of the fractional-order chaotic system without equilibrium is achieved and the improved controller is proven strictly by applying finite-time stable theorem. Eventually, simulation results verify the validity and rapidness of the proposed method. Therefore, the fractional-order chaotic system with hidden attractors can present better performance for practical applications, such as secure communication and image encryption, which deserve further investigation.


---

## Body

## 1. Introduction

Since Lorenz system [1] was put forward in the 1960s, chaos theory has made a great leap forward and researchers have never ceased to delve deeply into it. Among them, three-dimensional autonomous chaotic systems are investigated most extensively at home and abroad, such as Rössler system [2], Chen system [3], Lü system [4], and Liu system [5]. And these chaotic systems have achieved comprehensive application in real life, including secure communication, image encryption, medicine, and aerospace.In recent years, however, increasing attention has been paid to hidden attractors, which were discovered first in the Chua’circuit [6]. For those conventional chaotic systems above, their attractors are called self-excited attractors [7], whose basins of attraction can intersect unstable equilibrium points. In contrast, the attractors are defined as hidden attractors if its basins of attraction do not touch any unstable equilibrium point. The current literatures indicate that hidden attractors originate from several types of particular dynamical systems. And to be specific, the dynamical systems with one stable equilibrium [8, 9], a line equilibrium [10, 11], or no equilibrium [12–15] can all fall into this category. Even there are hidden attractors found in some special chaotic systems [16] which have both unstable equilibrium and stable equilibria. It is the presence of hidden attractors that pose a new challenge to Shilnikov criteria which illustrate that no less than one unstable equilibrium is needed to generate chaos. Meanwhile, with the lack of correlation between hidden attractors and equilibria, precise localization of these unpredictable attractors can be of great intractability. In addition, it also has been a focus of concern about coexistence of hidden attractors [17–20], which presents more complex and diverse dynamic behavior, and enhances the security ability in the field of cryptography application.However, the existing research findings about hidden attractors have been mostly focused on integer-order dynamical systems and relatively few efforts have been done for fractional-order systems with hidden attractors [16, 21–24]. The latter are worthy of further investigation owing to excellent characteristics of fractional calculus. To begin with, fractional-order systems, a further generalization, are more applicable by contrast with integer-order systems. Indeed, it can be said integer-order systems are the particular case of fractional-order systems and the idealization of real matter. Accordingly, fractional-order systems can describe the physical characteristics of the system engineering more clearly and reflect the essence and law of things in application domain. Then, fractional-order systems have special historical memory function, which are suitable for describing systems that rely on historical information. Furthermore, they can exhibit more complex and richer dynamic characteristics, with parameter selection more flexible. Thus, fractional calculus shows strong potential in practical engineering applications [25–28].Based on the above issues, this paper discusses the fractional-order form of a new 3D symmetric system without equilibrium and gives emphasis to the influence of the order on dynamic behavior of the system. The detailed dynamic analysis was carried out, mainly including existence of multiple attractors. And a finite-time synchronization of the fractional-order chaotic system with hidden attractors is realized. It has enormous value for the application of chaos technology in practical engineering field, especially in secure communication.This paper is organized as follows. In Section2, several definitions of fractional-order derivative and the fractional-order chaotic system model without equilibrium are introduced. In Section 3, dynamic analysis of the system is carried out at length. Section 4 gives the electronic circuit implementation results of the fractional-order system with hidden coexisting attractors. In Section 5, based on the finite-time stable theorem, a finite-time synchronization of the proposed system is implemented. Section 6 summarizes the conclusions.


## 2. Fractional-Order Form of the System without Equilibrium

### 2.1. Definitions of Fractional-Order Derivative

So far, a few different definitions [29, 30] of fractional-order derivative have been proposed, such as Grunwald–Letnikov (G-L), Riemann–Liouville (R-L), Caputo, and Conformable definitions. In general, the first three definitions are extensively employed and the Caputo definition is more suitable for engineering applications. In this paper, we conduct the study on the basis of the Caputo definition, which is given by(1)aCDtqft=1Γn-q×∫atfnτt-τq+1-ndτIn definition (1),n is the first integer which is not less thanq and n-1<q<n; Γ(·) is the Gamma function.


### 2.2. The Fractional-Order Chaotic System Model without Equilibrium

Sprott A system [31] is a three-dimensional conservative system and it has no equilibrium which depends on its third equation (1-y2). Based on Sprott A system, we construct a new 3D fractional-order chaotic system without equilibrium by imposing nonlinear functions associated with variablesx andy on the third equation of Sprott A system, which can generate diverse coexisting attractors, described as(2)dqxdtq=ydqydtq=-x-yzdqzdtq=ax2+by2+cxy+dwhere the parametersa,b,c,d are real constants,q is the fractional-order and 0<q<1. When setting the system parameters and the fractional-order as a=1.2, b=1, c=1, d=-1.3 and q=0.992, the fractional-order system (2) with initial conditions (x(0),y(0),z(0))=(1,1,0) presents hidden attractors as shown in Figure 1 and system (2) shows four-wing structure on thex-z plane. The Poincaré map on thex-zplane withy=0 shown in the three-dimensional space in Figure 2 also can verify it.Figure 1
The chaotic attractor of the fractional-order system (2): (a)x-y-z space; (b)x-y plane; (c)x-z plane; (d)y-z plane.
(a) (b) (c) (d)Figure 2
3D Poincaré map of the fractional-order system (2) on thex-z plane with y=0 (both directions through the cross-section in red and blue, respectively).To further confirm the chaotic motion, by the Benettin–Wolf algorithm, the finite-time local Lyapunov exponents on the time intervalt∈[0,1500] with the same initial conditions and parameters as above are shown in Figure 3. And it can be seen that the largest one of the finite-time local LEs on the time interval t∈[0,1500] always maintains to be positive.Figure 3
Dynamics of the finite-time local LEs of the fractional-order system (2) with initial conditions (x(0),y(0),z(0))=(1,1,0) on the time interval t∈[0,1500].Remark 1.
In numerical simulations, the Lyapunov exponents may differ significantly from different trajectories and only a finite-time interval can be considered for their computation. Therefore, we follow the concept of the finite-time local Lyapunov exponents in the article.


## 2.1. Definitions of Fractional-Order Derivative

So far, a few different definitions [29, 30] of fractional-order derivative have been proposed, such as Grunwald–Letnikov (G-L), Riemann–Liouville (R-L), Caputo, and Conformable definitions. In general, the first three definitions are extensively employed and the Caputo definition is more suitable for engineering applications. In this paper, we conduct the study on the basis of the Caputo definition, which is given by(1)aCDtqft=1Γn-q×∫atfnτt-τq+1-ndτIn definition (1),n is the first integer which is not less thanq and n-1<q<n; Γ(·) is the Gamma function.


## 2.2. The Fractional-Order Chaotic System Model without Equilibrium

Sprott A system [31] is a three-dimensional conservative system and it has no equilibrium which depends on its third equation (1-y2). Based on Sprott A system, we construct a new 3D fractional-order chaotic system without equilibrium by imposing nonlinear functions associated with variablesx andy on the third equation of Sprott A system, which can generate diverse coexisting attractors, described as(2)dqxdtq=ydqydtq=-x-yzdqzdtq=ax2+by2+cxy+dwhere the parametersa,b,c,d are real constants,q is the fractional-order and 0<q<1. When setting the system parameters and the fractional-order as a=1.2, b=1, c=1, d=-1.3 and q=0.992, the fractional-order system (2) with initial conditions (x(0),y(0),z(0))=(1,1,0) presents hidden attractors as shown in Figure 1 and system (2) shows four-wing structure on thex-z plane. The Poincaré map on thex-zplane withy=0 shown in the three-dimensional space in Figure 2 also can verify it.Figure 1
The chaotic attractor of the fractional-order system (2): (a)x-y-z space; (b)x-y plane; (c)x-z plane; (d)y-z plane.
(a) (b) (c) (d)Figure 2
3D Poincaré map of the fractional-order system (2) on thex-z plane with y=0 (both directions through the cross-section in red and blue, respectively).To further confirm the chaotic motion, by the Benettin–Wolf algorithm, the finite-time local Lyapunov exponents on the time intervalt∈[0,1500] with the same initial conditions and parameters as above are shown in Figure 3. And it can be seen that the largest one of the finite-time local LEs on the time interval t∈[0,1500] always maintains to be positive.Figure 3
Dynamics of the finite-time local LEs of the fractional-order system (2) with initial conditions (x(0),y(0),z(0))=(1,1,0) on the time interval t∈[0,1500].Remark 1.
In numerical simulations, the Lyapunov exponents may differ significantly from different trajectories and only a finite-time interval can be considered for their computation. Therefore, we follow the concept of the finite-time local Lyapunov exponents in the article.


## 3. Dynamic Analysis of the Fractional-Order System without Equilibrium

### 3.1. Basic Dynamic Properties

For system (2), it is symmetric about the z-axis on account of the coordinate transformation (x,y,z)→(-x,-y,z). And if let the right hand side of system (2) to zero, then it can be written as(3)y=0-x-yz=0ax2+by2+cxy+d=0Whilex=y=0, it is obvious that there is no equilibrium due to d≠0. Even so, through the above phase portraits, finite-time local Lyapunov exponents, and Poincaré map, the fractional-order system can show obvious chaotic characteristics and be regarded as the sort of hidden attractors.


### 3.2. Bifurcation Analysis of System Parameters

Bifurcation diagram can help us observe the further dynamic behaviors of the fractional-order system (2). While setting b=1, c=1, d=-1.3 and q=0.99, the bifurcation diagram and the largest finite-time local LEs varying with the parametera are shown in Figures 4 and 5. The performance via period-doubling bifurcations and hidden coexisting attractors are observed in Figure 6. The fractional-order system presents several different symmetric pair of coexisting limit cycles for a≤1.11. As shown on Figures 6(a) and 6(b), there coexist a symmetric pair of limit cycles with period-1 orbits and another symmetric limit cycle with period-2 orbits, respectively. For a>1.11, the fractional system begins to exhibit obvious chaotic oscillations and hidden coexisting chaotic attractors while the parameter is not over 1.21 as shown in Figure 6(c).Figure 4
Bifurcation diagram with the variable parametera, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 5
Largest finite-time local LEs on the time intervalt∈[0,1500] with the variable parametera, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 6
Coexisting attractors with initial conditions(x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red) when the parametera is kept as (a) a=0.9; (b) a=1.0; (c) a=1.2.
(a) (b) (c)When the parameters and fractional-order are chosen asa=1.2, b=1, c=1, and q=0.99, the bifurcation diagram of state variablexvarying withd is provided in Figure 7. Similarly, the states of the fractional-order system evolve from periodic state to multiperiodic state and then to chaos, which also exists periodic window in the neighborhood of d=-1.2. For purpose of further study in system parameter analysis, when keeping b=1, c=1 and q=0.99, Figure 8 plots the largest finite-time local LEs diagram with the mutual variation of parametersaandd, and the right color bar denotes different levels of largest Lyapunov exponent. In Figure 8, different colors intuitively represent influences of parameters on dynamical behavior where the warmer colors mean the higher values of largest finite-time local LEs.Figure 7
Bifurcation diagram with the variable parameterd, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 8
Dynamic map of Largest finite-time local LEs on the time intervalt∈[0,1500] with parametersa andb by selecting initial conditions (x(0),y(0),z(0))=(1,1,0).


### 3.3. Influences of Fractional-Order q and Initial Conditions

In Figure9, keeping the system parameters as a=1.2, b=1, c=1, and d=-1.3, the bifurcation of system (2) with the change of fractional-orderq illuminates the system could undergo the period-doubling bifurcation and end up in chaos as the order is set on q>0.982. And there exists diverse coexisting attractors under different orders. For q=0.99, a symmetric pair of strange attractors coexist with different initial conditions as shown in Figures 6(c) and 11, and its basins of attractions in the z=0 plane are shown in Figure 10, where the one kind of chaotic attractor is distinguished from another by red and orange and the yellow region indicates initial conditions in it can reach unbounded orbits. For q=0.98, there coexists a symmetric pair of limit cycles with period-4, as shown in Figure 12, and its basins in the z=0 plane are provided in Figure 13. Corresponding to the above analysis of basic dynamic properties, the basins have the similar symmetry about thez-axis.Figure 9
Bifurcation diagram with the variable fractional-orderq, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 10
Cross-section forz=0 of the basins of attraction for the symmetric pair of strange attractors (red and orange) with fractional-order q=0.99 (the yellow region can reach unbounded orbits).Figure 11
Coexisting chaotic attractors of the fractional-order system (2) with q=0.99, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red): (a)x-y plane; (b)x-z plane; (c)y-z plane.
(a) (b) (c)Figure 12
Coexisting limit cycles of the fractional-order system (2) with q=0.98, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 13
Cross-section forz=0 of the basins of attraction for the symmetric pair of limit cycles (blue and cyan) with fractional-order q=0.98 (the yellow region can reach unbounded orbits).Nevertheless, there are other cases of coexisting attractors. With order chosen asq=0.999, a strange attractor coexists with a limit cycle as shown in Figures 14 and 16, whose basins of attractions can be seen in Figure 15. The two attractors are respectively symmetric with respect to thez-axis and so do their basins. From the presented numerical results, the yellow region where initial values will lead to unbounded orbits largely remains invariant with different fractional orders, and the region of oscillation displays multitudinous and diverse dynamical behaviors.Figure 14
Coexistence of a chaotic attractor and a limit cycle of system (2) with q=0.999, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-0.5,5,0) (red).Figure 15
Cross-section forz=0 of the basins of attraction for the chaotic attractor (red) and the limit cycle (blue) with fractional-order q=0.999 (the yellow region can reach unbounded orbits).Figure 16
Coexistence of a chaotic attractor and a limit cycle of system (2) with q=0.999, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-0.5,5,0) (red): (a)x-y plane; (b)x-z plane; (c)y-z plane.
(a) (b) (c)


## 3.1. Basic Dynamic Properties

For system (2), it is symmetric about the z-axis on account of the coordinate transformation (x,y,z)→(-x,-y,z). And if let the right hand side of system (2) to zero, then it can be written as(3)y=0-x-yz=0ax2+by2+cxy+d=0Whilex=y=0, it is obvious that there is no equilibrium due to d≠0. Even so, through the above phase portraits, finite-time local Lyapunov exponents, and Poincaré map, the fractional-order system can show obvious chaotic characteristics and be regarded as the sort of hidden attractors.


## 3.2. Bifurcation Analysis of System Parameters

Bifurcation diagram can help us observe the further dynamic behaviors of the fractional-order system (2). While setting b=1, c=1, d=-1.3 and q=0.99, the bifurcation diagram and the largest finite-time local LEs varying with the parametera are shown in Figures 4 and 5. The performance via period-doubling bifurcations and hidden coexisting attractors are observed in Figure 6. The fractional-order system presents several different symmetric pair of coexisting limit cycles for a≤1.11. As shown on Figures 6(a) and 6(b), there coexist a symmetric pair of limit cycles with period-1 orbits and another symmetric limit cycle with period-2 orbits, respectively. For a>1.11, the fractional system begins to exhibit obvious chaotic oscillations and hidden coexisting chaotic attractors while the parameter is not over 1.21 as shown in Figure 6(c).Figure 4
Bifurcation diagram with the variable parametera, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 5
Largest finite-time local LEs on the time intervalt∈[0,1500] with the variable parametera, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 6
Coexisting attractors with initial conditions(x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red) when the parametera is kept as (a) a=0.9; (b) a=1.0; (c) a=1.2.
(a) (b) (c)When the parameters and fractional-order are chosen asa=1.2, b=1, c=1, and q=0.99, the bifurcation diagram of state variablexvarying withd is provided in Figure 7. Similarly, the states of the fractional-order system evolve from periodic state to multiperiodic state and then to chaos, which also exists periodic window in the neighborhood of d=-1.2. For purpose of further study in system parameter analysis, when keeping b=1, c=1 and q=0.99, Figure 8 plots the largest finite-time local LEs diagram with the mutual variation of parametersaandd, and the right color bar denotes different levels of largest Lyapunov exponent. In Figure 8, different colors intuitively represent influences of parameters on dynamical behavior where the warmer colors mean the higher values of largest finite-time local LEs.Figure 7
Bifurcation diagram with the variable parameterd, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 8
Dynamic map of Largest finite-time local LEs on the time intervalt∈[0,1500] with parametersa andb by selecting initial conditions (x(0),y(0),z(0))=(1,1,0).


## 3.3. Influences of Fractional-Order q and Initial Conditions

In Figure9, keeping the system parameters as a=1.2, b=1, c=1, and d=-1.3, the bifurcation of system (2) with the change of fractional-orderq illuminates the system could undergo the period-doubling bifurcation and end up in chaos as the order is set on q>0.982. And there exists diverse coexisting attractors under different orders. For q=0.99, a symmetric pair of strange attractors coexist with different initial conditions as shown in Figures 6(c) and 11, and its basins of attractions in the z=0 plane are shown in Figure 10, where the one kind of chaotic attractor is distinguished from another by red and orange and the yellow region indicates initial conditions in it can reach unbounded orbits. For q=0.98, there coexists a symmetric pair of limit cycles with period-4, as shown in Figure 12, and its basins in the z=0 plane are provided in Figure 13. Corresponding to the above analysis of basic dynamic properties, the basins have the similar symmetry about thez-axis.Figure 9
Bifurcation diagram with the variable fractional-orderq, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 10
Cross-section forz=0 of the basins of attraction for the symmetric pair of strange attractors (red and orange) with fractional-order q=0.99 (the yellow region can reach unbounded orbits).Figure 11
Coexisting chaotic attractors of the fractional-order system (2) with q=0.99, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red): (a)x-y plane; (b)x-z plane; (c)y-z plane.
(a) (b) (c)Figure 12
Coexisting limit cycles of the fractional-order system (2) with q=0.98, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-1,-1,0) (red).Figure 13
Cross-section forz=0 of the basins of attraction for the symmetric pair of limit cycles (blue and cyan) with fractional-order q=0.98 (the yellow region can reach unbounded orbits).Nevertheless, there are other cases of coexisting attractors. With order chosen asq=0.999, a strange attractor coexists with a limit cycle as shown in Figures 14 and 16, whose basins of attractions can be seen in Figure 15. The two attractors are respectively symmetric with respect to thez-axis and so do their basins. From the presented numerical results, the yellow region where initial values will lead to unbounded orbits largely remains invariant with different fractional orders, and the region of oscillation displays multitudinous and diverse dynamical behaviors.Figure 14
Coexistence of a chaotic attractor and a limit cycle of system (2) with q=0.999, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-0.5,5,0) (red).Figure 15
Cross-section forz=0 of the basins of attraction for the chaotic attractor (red) and the limit cycle (blue) with fractional-order q=0.999 (the yellow region can reach unbounded orbits).Figure 16
Coexistence of a chaotic attractor and a limit cycle of system (2) with q=0.999, initial conditions (x(0),y(0),z(0))=(1,1,0) (blue) and (x(0),y(0),z(0))=(-0.5,5,0) (red): (a)x-y plane; (b)x-z plane; (c)y-z plane.
(a) (b) (c)


## 4. Circuit Implementation Results of the Fractional-Order System

Circuit implementation of the fractional-order chaotic system plays a crucial role in actual application, which can also verify whether the results obtained by theory analysis and numerical computation are consistent and correct or not. Based on PSpice software, the schematic of the fractional-order chaotic system with orderq=0.992 and parameters a=1.2, b=1, c=1, and d=-1.3 is designed and implemented as shown in Figure 17. In Figure 17, AD633 and OP-07 are adopt as all analogue multipliers and operational amplifiers.Figure 17
The electronic circuit schematic of the fractional-order system (2) in PSpice.Therein, with the proposed transfer function approximations method [32] utilized and ωmax=100, PT=0.01 assumed, the fractional integrator is derived by (4). In engineering practice it can be realized by three chain fractances composed of three resistor-capacitor pairs, as shown by (5). With comparing (4) with (5), the values of these circuit components are determined as Ra=95957.169kΩ, Ca=1.002μF, Rb=1.153kΩ, Cb=13.893μF and Rc=0.205Ω, Cc=12.959μF.(4)1s0.99=1.1471s+58.2265s+3.5105×105s+0.0104s+62.4256s+3.7636×105(5)Fs=1/Cas+1/RaCa+1/Cbs+1/RbCb+1/Ccs+1/RcCc=1C0C0/Ca+C0/Cb+C0/Ccs2+Cb+Cc/Ra+Ca+Cc/Rb+Ca+Cb/Rc/CaCb+CbCc+CaCcs+Ra+Rb+Rc/RaRbRc/CaCb+CbCc+CaCcs+1/RaCas+1/RbCbs+1/RcCcWith the Kirchhoff’s circuit laws, the mathematical model of system (2) is described by the circuit equations as follow:(6)dqxdtq=R7R1R6C0ydqydtq=-1R3C0x-1ηR2C0yzdqzdtq=R5ηR8R4C0x2+R7ηR9R6C0y2+R7ηR10R6C0xy+R12R13R11C0VCC-where C0=1μF and 1/η=1/10 hinged on the characteristic of AD633. The parameters of electronic components are calculated as: VCC+=-VCC-=15VDC, R1=R3=R11=1000kΩ, R2=R4=R5=R6=R7=R9=R10=100kΩ, R8=83.33kΩ, R12=10kΩ, R13=115.38kΩ. Figure 18 demonstrates circuit implementation results in which hidden strange attractors can been observed distinctly. And it also indicates that the circuital results are comparatively consistent with the numerical simulations.Figure 18
Hidden chaotic attractors of the fractional-order system (2) in PSpice: (a) on thex-y plane; (b) on thex-z plane; (c) on they-z plane.
(a) (b) (c)


## 5. Finite-Time Synchronization of the Fractional-Order System with Hidden Coexisting Attractors

There is a wide variety of literature [24, 33–37] on synchronization of fractional-order chaotic systems by reason of comprehensive practical application in many fields, such as aerospace, power system, secure communication and so on. Research about finite-time synchronization in allusion to fractional-order chaotic system with hidden attractors, are quite limited. In this section, based on the finite-time stable theorem [38], in consideration of the disturbances, the modified controller is devised for synchronization of the fractional-order system with hidden chaotic attractors and achieves the anticipated effect.Theorem 2 (see [38]).
For the fractional-order system in general, if it satisfies(7)xaCDtqxT=Γ2Γ2+qaCDtqxxqT≤-kxxTα,α<q+q22,k>0,where x=[x1,x2,⋯xn], xq=[x1q,x2q,⋯xnq]. Then the state variables x of the given fractional-order system will converge to zero within a finite-time:(8)t=v0q-2α/1+qΓ1+q-2α/1+qΓ1+qΓ1+2q-2α/1+qkΓ2+q1/q,where v=xxqT.Lemma 3.
∀ a , b > 0, 0<c<1: (9)a+bc≤ac+bc
Making allowances for the internal disturbancesΔfim(X,t) and Δfis(X1,t) (i=1,2,3 and x,y,z∈X, x1,y1,z1∈X1) in the fractional-order chaotic system (2), the master system and the slave system can be obtained as(10)dqxdtq=y+Δf1mX,tdqydtq=-x-yz+Δf2mX,tdqzdtq=ax2+by2+cxy+d+Δf3mX,t(11)dqx1dtq=y1+Δf1sX1,t-u1dqy1dtq=-x1-y1z1+Δf2sX1,t-u2dqz1dtq=ax12+by12+cx1y1+d+Δf3sX1,t-u3where a=1.2, b=1, c=1, d=-1.3, and q=0.99; u=[u1,u2,u3] is the designed controller.Assumption 4.
The disturbancesΔfim(X,t), Δfis(X1,t) are bounded. In other words, positive constants Fim and Fis can be found to satisfy(12)ΔfimX,t≤Fim,ΔfisX1,t≤Fis
Afterwards, we define the error system ase1=x1-x, e2=y1-y, e3=z1-z, which can ultimately be derived as(13)dqe1dtq=e2+Δf1sX1,t-Δf1mX,t-u1dqe2dtq=-e1-y1e3-ze2+Δf2sX1,t-Δf2mX,t-u2dqe3dtq=1.2x1e1+1.2xe1+y1e2+ye2+x1e2+ye1+Δf3sX1,t-Δf3mX,t-u3Theorem 5.
It can realize the synchronization of the master system (10) and the slave system (11) within a finite-time, while the control input is provided as(14)u1=1.2x1e3+ye3+ke1β+k1signe1u2=-ze2+x1e3+ke2β+k2signe2u3=1.2xe1+ye2+ke3β+k3signe3where ki=Fis+Fim, i=1,2,3.Proof.
According to (7) and (12), the corresponding function with the error system (13) and the control input (14) can be constructed as(15)Γ2Γ2+qdqdtqe1e2e3e1qe2qe3qT=e1e2e3dqdtqe1dqdtqe2dqdtqe3T=e1e2-e1e2-y1e2e3-ze22+1.2x1e1e3+1.2xe1e3+y1e2e3+ye2e3+x1e2e3+ye1e3+Δf1sX1,t-Δf1mX,te1+Δf2sX1,t-Δf2mX,te2+Δf3sX1,t-Δf3mX,te3-u1e1-u2e2-u3e3=-ke1β+1-ke2β+1-ke3β+1+∑i=13ΔfisX1,t-ΔfimX,tei-kieisignei≤-ke1β+1-ke2β+1-ke3β+1+∑i=13Fis+Fimei-kieisignei=-ke1β+1-ke2β+1-ke3β+1+∑i=13kiei-eisignei≤-ke1β+1-ke2β+1-ke3β+1where ei-eisign(ei)=ei-ei≤0.
By reference to Lemma3, (15) can be derived as(16)-ke1β+1-ke2β+1-ke3β+1=-ke12β+1/2-ke22β+1/2-ke32β+1/2≤-ke12+e22+e32β+1/2=-keeTβ+1/2
Therefore, from Theorem2, the error system (13) states will converge to zero within a finite-time.The numerical simulations are carried out to verify the effectiveness of the proposed scheme. The disturbancesΔfim(X,t), Δfis(X1,t) are chosen as (17). The parameters for controllers are set as k=1.5, β=0.5 and k1=1.0, k2=1.4, k3=4.4. Initial states of the master system and the slave system are selected as (x(0),y(0),z(0))=(-0.5,-0.5,2) and (x1(0),y1(0),z1(0))=(1,1,0). As shown in Figure 19, the error system (13) can stabilize to zero within 1s. Meanwhile, time response of the master system (10) and the slave system (11) under control input (14) is presented in Figure 20 and the system states of the slave system will track the trajectories of the master system within 1s. In comparison with the original controllers which do not take into account disturbance factors in [38], the improved controllers in this paper can synchronize the master system and the slave system under disturbances. The simulation results in Figure 20 also indicate that the proposed method is fast and can achieve chaos synchronization within finite-time. All the results verify the effectiveness of the proposed scheme in synchronization of fractional-order chaotic system with hidden coexisting attractors.(17)Δf1mX,t=0.15xcos⁡2t,Δf1sX1,t=-0.2x1sin⁡2tΔf2mX,t=-0.2ysin⁡3t,Δf2sX1,t=0.1y1sin⁡2tΔf3mX,t=0.1zcos⁡t,Δf3sX1,t=0.15z1cos⁡tFigure 19
Synchronization errors between the master and slave systems.Figure 20
Synchronization results of the master and slave systems: (a)x-x1; (b)y-y1; (c)z-z1.
(a) (b) (c)


## 6. Conclusions

In this paper, a novel three-dimensional symmetric fractional-order system without equilibrium, which can generate various hidden coexisting attractors, has been introduced. We have analyzed dynamical characteristics of the fractional-order system in detail by numerical simulations and found the system presents abundant complex dynamic characteristics. Especially, through bifurcation diagram, the basins of attraction, and phase diagram, the influence of the order on dynamic behavior of the fractional-order system has been explored in depth. With the approximations of fractional-order integratorq=0.992 derived, circuit implementation of the fractional-order system with hidden coexisting attractors has been achieved in PSpice and it has a good agreement with the numerical simulations. Furthermore, the modified controller has been devised and proved, which can realize the finite-time synchronization of the fractional-order system without equilibrium under the disturbances. Finally, simulation results also illustrate the robustness and speediness of the proposed controller. In order to simplify control inputs, fractional-order finite-time backstepping method and dynamic surface method will be considered in the future. Furthermore, secure communication and image encryption based on the novel fractional-order chaotic system with hidden coexisting attractors will also be further investigated.



---

*Source: 6908607-2019-07-22.xml*