#### Abstract

A 3D fractional-order nonlinear system with coexisting chaotic attractors is proposed in this paper. The necessary condition of the existence chaos is . The fractional-order system exhibits chaotic attractors with the order as low as 2.5431. The largest Lyapunov exponent varying as fractional order is given. Furthermore, there are the coexisting “positive attractor” and “negative attractor” in this fractional-order chaotic system, and the necessary condition for “positive attractor” and “negative attractor” is obtained. Meanwhile, a control scheme for the stabilization of the unstable equilibrium is suggested via a single state variable linear controller. Numerical results show that the control scheme is valid.

#### 1. Introduction

Chaotic behaviors in nonlinear is a very interesting phenomenon. The high irregularity, unpredictability, and complexity in chaotic systems [1, 2] have been widely used in the field of engineering and technology such as secure communications, image steganography, authenticated encryption, motor control, and power system protection. Recently, coexisting chaotic attractors have been found in chaotic systems [3–9]. For example, the coexisting chaotic attractors in a 3D no-equilibrium system were reported by Pham et al. [3], the coexisting multiple attractors in Hopfield neural network were found by Bao et al. [4], the coexisting chaotic attractors in a hyperchaotic hyperjerk system were given by Wang et al. [5], the coexisting “positive attractor” and “negative attractor” in a 3D autonomous continuous chaotic system were found by Zhou and Ke [6], and so on [7–9]. Therefore, more and more attention has been focused on the coexisting chaotic attractors in nonlinear chaotic systems.

On the other hand, the fractional-order differential equations [10–12] can be accurately described in the real-world physical systems such as viscoelasticity, dielectric polarization, electrode-electrolyte polarization, electromagnetic waves, heat conduction, diffusion-wave, and superdiffusion. Chaotic behaviors have been found in many real-world physical fractional-order nonlinear systems, for example, the fractional-order chaotic brushless DC motor [13], the fractional-order electronic circuits [14], the fractional-order microelectromechanical system [15], and the fractional-order gyroscopes [16]. Therefore, more and more attention has been paid to the chaotic behaviors in fractional-order nonlinear systems.

Motivated by the above discussions, based on a 3D autonomous continuous chaotic system reported by Zhou and Ke [6], we suggested a 3D autonomous continuous fractional-order system. We have shown that the chaotic system reported by Zhou and Ke [6] can be extended to its fractional-order version where the coexisting “positive attractor” and “negative attractor” can be observed. We obtained that the fractional-order system with the order as low as 2.5431 exhibits chaotic attractors. Moreover, we obtained the largest Lyapunov exponent varying as fractional order. Finally, for the stabilization of the unstable equilibrium, one control scheme is proposed via a single state variable linear controller.

The outline of this paper is organized as follows. In Section 2, based on a 3D autonomous continuous chaotic system reported by Zhou and Ke [6], the fractional-order version nonlinear system is given, and some basic dynamical properties of this fractional-order version nonlinear system are obtained including the necessary condition of the existence chaos, the largest Lyapunov exponent varying as fractional order, and the coexisting “positive attractor” and “negative attractor.” In Section 3, by a single state variable, stabilization of the unstable equilibrium points of the fractional-order chaotic system is discussed. Finally, the conclusions are given in Section 4.

#### 2. System Model and Basic Characteristics

In this paper, the Caputo definition of the fractional derivative will be used in next. The Caputo definition of the fractional derivative is described aswhere is the Caputo operator, is the first integer which is not less than , and is the -order derivative in usual sense of .

Next, based on the 3D autonomous continuous chaotic system reported by Zhou and Ke [6], a fractional-order system is addressed aswhere fractional order . Fractional-order system (2) has five equilibrium points. They are , , , , and , respectively.

Now, we can obtain the eigenvalues of the five equilibrium points. The eigenvalues of equilibrium point are (−6, −1, 4). Thus, equilibrium point is an unstable saddle point of index one. The eigenvalues of equilibrium points and are (−5.3957, *j*, *j*). Thus, equilibrium points and are unstable saddle points of index two. The eigenvalues of equilibrium points and are (−6.7558, *j*, *j*). Thus, equilibrium points and are unstable saddle points of index two.

Tavazoei and Haeri [17] have obtained that a necessary condition for a fractional-order nonlinear system to exist chaotic iswhere is the eigenvalues of saddle equilibrium point of index two in fractional-order nonlinear system.

Now, we can obtain the necessary condition of the existence chaos in fractional-order system (2). According to (3), we have the following:Thus, the necessary condition of existence chaos in fractional-order system (2) is . This result indicates that fractional-order system (2) with the order as low as 2.5431 can exhibit chaotic attractors.

In this paper, the improved version of Adams-Bashforth-Moulton [18] (denoted by IVABM) numerical algorithm is used to deal with fractional-order system (2). The IVABM numerical algorithm will be introduced next. Now, consider the fractional-order systemwith initial condition . Let , and . By IVABM numerical algorithm, system (5) can be discretized as follows:where

The error of this IVABM numerical algorithm is

We can yield the largest Lyapunov exponent varying as via IVABM numerical algorithm, which is shown in Figure 1, where and the initial condition is (−5, 2, 5).

According to Figure 1, we can obtain that the largest Lyapunov exponent is larger zero for . This result indicates that the chaotic attractor is emerged in fractional-order system (2) for . For example, let , the largest Lyapunov exponent is 0.9853, and the chaotic attractor in fractional-order system (2) with initial condition (−5, 2, 5) is shown as Figure 2.

**(a)**

**(b)**

Same as the results reported by Zhou and Ke [6], there are coexisting “positive attractor” and “negative attractor” in fractional-order system (2). The chaotic attractor in fractional-order system (2) can be decided by the initial conditions; that is, the chaotic attractor depends on the distance from the initial point (initial condition to the unstable points [6]. Next, the distances from the initial point to points , , , and are denoted by , , , and , respectively. Same as [6], the following results can be obtained:

(1) If the initial point (initial condition) is near the unstable or , there is the same chaotic attractor in fractional-order chaotic system (2), which is called “positive attractor” by Zhou and Ke [6], where the “positive attractor” means and a necessary condition for “positive attractor” is [6].

(2) If the initial point (initial condition) is near the unstable or , there is the same chaotic attractor in fractional-order chaotic system (2), which is called “negative attractor” by Zhou and Ke [6], where the “positive attractor” means and a necessary condition for “negative attractor” is [6].

Next, some numerical simulations are given for . Here, the largest Lyapunov exponent is 0.763 for .

*Case 1 (the initial point is near unstable saddle point or ). *For example, choose the initial conditions as (). Therefore, , , , and are 0.7636, 8.8563, 6.1735, and 9.7172, respectively. Therefore, the initial point () is near unstable saddle point . Therefore, the fractional-order system (2) has the “positive attractor,” which is shown as Figure 3.

**(a)**

**(b)**

For example, choose the initial conditions as (). Therefore, , , , and are 8.8563, 0.7636, 9.7172, and 6.1735, respectively. Therefore, the initial point () is near unstable saddle point . Thus, fractional-order system (2) has the “negative attractor,” which is shown as Figure 3.

*Case 2 (the initial point is near unstable saddle point or ). *For example, choose the initial conditions as (). Therefore, , , , and are 5.5875, 9.1586, 1.2594, and 10.9026, respectively. Therefore, the initial point () is close to unstable saddle point . Thus, fractional-order system (2) has the “positive attractor,” which is shown as Figure 4.

**(a)**

**(b)**

For example, choose the initial conditions as (−4, −3, −2). Therefore, , , , and are 9.1586, 5.5875, 10.9026, and 1.2594, respectively. Therefore, the initial point (−4, −3, −2) is close to unstable saddle point . Thus, fractional-order system (2) has the “negative attractor,” which is shown as Figure 4.

According to Figures 3 and 4, the coexisting chaotic attractors are found in fractional-order chaotic system (2). These results indicate that the chaotic system reported by Zhou and Ke [6] can be extended to its fractional-order version where the coexisting “positive attractor” and “negative attractor” can be observed.

*Remark 1. *There are overlaps between the coexisting chaotic attractors in [3–5, 7–9]. However, there are two isolated chaotic attractors in fractional-order chaotic system (2); that is, there are no overlaps between the coexisting chaotic attractors in fractional-order chaotic system (2).

#### 3. Stabilization of the Unstable Equilibrium Points of the Fractional-Order Chaotic System (2)

First, a result on stability of fractional-order nonlinear system is recalled. Consider the fractional-order nonlinear system as follows:where , denotes the state vector, is a nonlinear vector field, is a constant matrix, and and denote the nonlinear and linear parts of .

Lemma 2 (for more details, see [19]). *Fractional-order nonlinear system (9) is said to be asymptotically stable if the next two conditions hold:*(i)*, and ,*(ii)* , and ,** where and denote the eigenvalues and the -norm with respect to matrix , respectively.*

Next, we discuss how to stable the unstable equilibrium point in fractional-order chaotic system (2) via single state variable linear controller. Now, let the unstable equilibrium points in system (2) be , and the controlled fractional-order system is shown as follows:where and is a linear controller determined by a single state variable . are feedback gains determined later.

Theorem 3. *Let real matrix beand choose suit feedback gains ; if , then the unstable equilibrium points of fractional-order chaotic system (2) in controlled system (10) are asymptotically stable.*

*Proof. *Let , , and ; then controlled system (10) can be changed asThus, the unstable equilibrium points in fractional-order chaotic system (2) are changed as the point in system (12).

One can easily obtain that point is the origin of system (12). Therefore, the problem of stabilization of the unstable equilibrium points in fractional-order chaotic system (2) is turned to the problem of stabilization of fractional-order system (12). If the origin of system (12) is asymptotically stable, then the unstable equilibrium points of fractional-order chaotic system (2) in controlled system (10) are asymptotically stable.

Now, controlled system (12) can be rewritten aswhere .

First, it is easy to obtain that , andThus, the first condition in Lemma 2 is satisfied.

Second, according to, and using , we can yield that the second condition in Lemma 2 is satisfied.

Therefore, according to Lemma 2, the origin of system (12) is asymptotically stable. That is, the unstable equilibrium points of fractional-order chaotic system (2) in controlled system (10) are asymptotically stable. The proof is finished.

*Remark 4. *In our control scheme, the linear controller is only determined by one single state variable, so our control scheme is different from many previous control schemes.

#### 4. Numerical Simulations Results

Next, in order to show the effectiveness of the proposed control approach, the numerical simulations are performed for .

For the unstable equilibrium point , if we choose , , then . Thus, we obtain , , , and . According to Theorem 3, the unstable equilibrium point in controlled system (10) is asymptotically stable. Figure 5 displays the time evolution of state variables. Here, the initial conditions be ().

For the unstable equilibrium point , if we choose , , then . Thus, we obtain *j*, and . According to Theorem 3, the unstable equilibrium point in controlled system (10) is asymptotically stable. Figure 6 shows the time evolution of state variables. Here, the initial conditions are ().

For the unstable equilibrium point , if we choose , , then . Therefore, we obtain , *j*, and . According to Theorem 3, the unstable equilibrium point in controlled system (10) is asymptotically stable. Figure 7 shows the time evolution of state variables. Here, the initial conditions are ().

For the unstable equilibrium point , if we choose , , , then . Thus, we obtain *j*, and . According to Theorem 3, the unstable equilibrium point in controlled system (10) is asymptotically stable. Figure 8 shows the time evolution of state variables. Here, the initial conditions are ().

For the unstable equilibrium point , if we choose , , , then . Thus, we obtain *j*, and . According to Theorem 3, the unstable equilibrium point in controlled system (10) is asymptotically stable. Figure 9 shows the time evolution of state variables. Here, the initial conditions are ().

The simulative results in Figures 5–9 show the effectiveness of our control scheme.

#### 5. Conclusions

In this paper, a fractional-order chaotic system is proposed. The necessary condition of existence chaos in this fractional-order chaotic system is . The largest Lyapunov exponent varying as fractional order is given. The coexisting “positive attractor” and “negative attractor” can be observed, and the necessary conditions for “positive attractor” and “negative attractor” are obtained. Meanwhile, by a single state variable, a linear controller is used for the stabilization of the unstable equilibrium points of the fractional-order chaotic system. The numerical results show that the control approach is effective.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare no conflicts of interest.