AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 10.1155/2014/326354 326354 Research Article A New Fractional-Order Chaotic Complex System and Its Antisynchronization Jiang Cuimei 1 Liu Shutang 1 Luo Chao 2,3 Sadyrbaev Felix 1 College of Control Science and Engineering Shandong University Jinan 250061 China sdu.edu.cn 2 School of Information Science and Engineering Shandong Normal University Jinan 250014 China sdnu.edu.cn 3 Shandong Provincial Key Laboratory for Novel Distributed Computer Software Technology Jinan 250014 China 2014 19102014 2014 28 03 2014 28 07 2014 28 08 2014 20 10 2014 2014 Copyright © 2014 Cuimei Jiang 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.

We propose a new fractional-order chaotic complex system and study its dynamical properties including symmetry, equilibria and their stability, and chaotic attractors. Chaotic behavior is verified with phase portraits, bifurcation diagrams, the histories, and the largest Lyapunov exponents. And we find that chaos exists in this system with orders less than 5 by numerical simulation. Additionally, antisynchronization of different fractional-order chaotic complex systems is considered based on the stability theory of fractional-order systems. This new system and the fractional-order complex Lorenz system can achieve antisynchronization. Corresponding numerical simulations show the effectiveness and feasibility of the scheme.

1. Introduction

Chaotic behavior and synchronization of fractional-order dynamical systems have been extensively studied over the last decade. Many fractional-order systems can behave chaotically, such as the fractional-order Chua’s system , the fractional Rössler system , the fractional-order Lorenz system , the fractional-order Chen system , and the fractional-order Lü system . It has been shown that some fractional-order systems have chaotic behavior with orders less than 3. Meanwhile, chaos synchronization of fractional-order systems has attracted much attention, such as the complete synchronization (CS) , projective synchronization (PS) , and lag projective synchronization .

However, most of the studies about fractional-order systems had been based on real variables, and complex systems are rarely involved. Complex systems provide an excellent instrument to describe a variety of physical phenomena, such as detuned laser systems, amplitudes of electromagnetic fields, and thermal convection of liquid flows . And now complex systems have played an important role in many branches of physics, for example, superconductors, plasma physics, geophysical fluids, modulated optical waves, and electromagnetic fields . There are some new kinds of synchronization for complex dynamical systems, for example, complex complete synchronization (CCS) , complex projective synchronization (CPS) , complex modified projective synchronization (CMPS) [15, 16], and so forth. These new kinds of synchronization have been widely studied for applications in secure communication , because complex variables (doubling the number of variables) increase the contents and security of the transmitted information. Therefore, the dynamical behavior and synchronization of the fractional-order complex nonlinear systems are worth studying. Recently, Luo and Wang proposed the fractional-order complex Lorenz system  and the fractional-order complex Chen system  and studied their dynamical properties and chaos synchronization. To our best knowledge, there are few results on fractional-order chaotic complex systems until now.

Motivated by the above discussion, the aim of this paper is to investigate the chaotic phenomena in a newly proposed fractional-order complex Lü system, which may provide potential applications in secure communication. As will be shown below, this new system displays many interesting dynamical behaviors, such as fixed points, periodic motions, and chaotic motions. Besides, when the parameters of the system are fixed, the lowest order for chaos to exist is determined. Furthermore, antisynchronization between the new system and fractional-order complex Lorenz system is studied. More generally, we investigate antisynchronization of different fractional-order chaotic complex systems and give a usual scheme.

The remainder of this paper is organized as follows. In Section 2, The fractional-order complex Lü system is presented and its dynamics is discussed by phase portraits, bifurcation diagrams, the histories, and the largest Lyapunov exponents. In Section 3, the antisynchronization of different fractional-order chaotic complex systems is studied, and the proposed new system can antisynchronize the fractional-order complex Lorenz system. A concluding remark is given in Section 4.

2. The Fractional-Order Complex Lü System 2.1. The Proposal of the Fractional-Order Complex Lü System

There are many definitions of fractional derivatives [20, 21], such as Riemann-Liouville, Grünwald-Letnikov, and Caputo definitions. In this paper, we use the Caputo definition which is defined as follows: (1) D * α x ( t ) = J n - α x ( n ) ( t ) . Here n is the first integer which is not less than α and α > 0 , x ( n ) is the n -order derivative in the usual sense, and J β ( β > 0 ) is the β -order Riemann-Liouville integral operator with expression (2) J β φ ( t ) = 1 Γ ( β ) 0 t ( t - τ ) β - 1 φ ( τ ) d τ . Here Γ stands for Gamma function, and the operator D * α is generally called α -order Caputo differential operator.

In 2007, the complex Lü system was proposed by Mahmoud et al. , which can be described as (3) x ˙ 1 = a 1 ( x 2 - x 1 ) , x ˙ 2 = - x 1 x 3 + a 2 x 2 , x ˙ 3 = 1 2 ( x ¯ 1 x 2 + x 1 x ¯ 2 ) - a 3 x 3 , where x = ( x 1 , x 2 , x 3 ) T is the vector of state variables, x 1 = m 1 + j m 2 and x 2 = m 3 + j m 4 are complex variables, and x 3 = m 5 is a real variable. j = - 1 , a i > 0    ( i = 1,2 , 3 ) is the system real parameter. When the parameters are chosen as a 1 = 40 ,   a 2 = 22 ,   a 3 = 5 , the system (3) is chaotic as shown in Figure 1.

Chaotic attractor of complex Lü system ( a 1 , a 2 , a 3 ) = ( 40,22,5 ) .

In this paper, we modify the derivative operator in (3) to be with respect to the fractional order α    ( 0 < α 1 ) . Thus, the fractional-order complex Lü system can be expressed by (4) D * α 1 x 1 = a 1 ( x 2 - x 1 ) , D * α 2 x 2 = - x 1 x 3 + a 2 x 2 , D * α 3 x 3 = 1 2 ( x ¯ 1 x 2 + x 1 x ¯ 2 ) - a 3 x 3 ; when the derivative order α i = 1    ( i = 1,2 , 3 ) , system (4) will be the common integer-order complex Lü system. The Caputo differential operator is a linear operator ; that is, D * α ( a y ( t ) + b z ( t ) ) = a D * α ( y ( t ) ) + b D * α ( z ( t ) ) , for arbitrary constants a and b . Therefore, D * α i ( x i ) = D * α i ( m 2 i - 1 + j m 2 i ) = D * α i m 2 i - 1 + j D * α i m 2 i ,   ( i = 1,2 , 3 ) .

Separating the real and imaginary parts of system (4), we can obtain the following system: (5) D * α 1 m 1 = a 1 ( m 3 - m 1 ) , D * α 1 m 2 = a 1 ( m 4 - m 2 ) , D * α 2 m 3 = - m 1 m 5 + a 2 m 3 , D * α 2 m 4 = - m 2 m 5 + a 2 m 4 , D * α 3 m 5 = m 1 m 3 + m 2 m 4 - a 3 m 5 .

2.2. Numerical Algorithm for the Fractional-Order Complex Lü System

In 2002, Diethelm et al. proposed the Adams-Bashforth-Moulton predictor-corrector scheme , which is numerically stable and can be used to both linear and nonlinear fractional differential equations. According to this algorithm, system (5) for initial condition ( m 10 , m 20 , m 30 , m 40 , m 50 ) can be discretized as (6) m 1 , n + 1 = m 10 + h α 1 Γ ( α 1 + 2 ) × [ a 1 ( m 3 , n + 1 p - m 1 , n + 1 p ) j = 0 n iiiiiiiii + j = 0 n γ 1 , j , n + 1 a 1 ( m 3 j - m 1 j ) ] , m 2 , n + 1 = m 20 + h α 1 Γ ( α 1 + 2 ) × [ a 1 ( m 4 , n + 1 p - m 2 , n + 1 p ) j = 0 n i i i i i i i i i + j = 0 n γ 1 , j , n + 1 a 1 ( m 4 j - m 2 j ) ] , m 3 , n + 1 = m 30 + h α 2 Γ ( α 2 + 2 ) × [ a 2 m 3 , n + 1 p - m 1 , n + 1 p m 5 , n + 1 p j = 0 n iiiiiiiiiii + j = 0 n γ 2 , j , n + 1 ( a 2 m 3 j - m 1 j m 5 j ) ] , m 4 , n + 1 = m 40 + h α 2 Γ ( α 2 + 2 ) × [ a 2 m 4 , n + 1 p - m 2 , n + 1 p m 5 , n + 1 p j = 0 n iiiiiiiiiiii + j = 0 n γ 2 , j , n + 1 ( a 2 m 4 j - m 2 j m 5 j ) ] , m 5 , n + 1 = m 50 + h α 3 Γ ( α 3 + 2 ) × [ m 1 , n + 1 p m 3 , n + 1 p + m 2 , n + 1 p m 4 , n + 1 p - a 3 m 5 , n + 1 p j = 0 n iiiiiiiiiiii + j = 0 n γ 3 , j , n + 1 ( m 1 j m 3 j + m 2 j m 4 j - a 3 m 5 j ) ] ,

where (7) m 1 , n + 1 p = m 10 + 1 Γ ( α 1 ) j = 0 n ω 1 , j , n + 1 a 1 ( m 3 j - m 1 j ) , m 2 , n + 1 p = m 20 + 1 Γ ( α 1 ) j = 0 n ω 1 , j , n + 1 a 2 ( m 4 j - m 2 j ) , m 3 , n + 1 p = m 30 + 1 Γ ( α 2 ) j = 0 n ω 2 , j , n + 1 ( a 2 m 3 j - m 1 j m 5 j ) , m 4 , n + 1 p = m 40 + 1 Γ ( α 2 ) j = 0 n ω 2 , j , n + 1 ( a 2 m 4 j - m 2 j m 5 j ) , m 5 , n + 1 p = m 50 + 1 Γ ( α 3 ) j = 0 n ω 3 , j , n + 1 ( m 1 j m 3 j + m 2 j m 4 j - a 3 m 5 j ) , γ i , j , n + 1 = { n α i + 1 - ( n - α i ) ( n + 1 ) α i , j = 0 , ( n - j + 2 ) α i + 1 + ( n - j ) α i + 1 - 2 ( n - j + 1 ) α i + 1 , 1 j n , 1 ,    j = n + 1 , ω i , j , n + 1 = h α i α i ( ( n + 1 - j ) α i - ( n - j ) α i ) .

2.3. Dynamics of the Fractional-Order Complex Lü System 2.3.1. Symmetry and Invariance

Note that the symmetry of system (5) is symmetric about m 5 -axis, which means it is invariant for the coordinate transformation of ( m 1 , m 2 , m 3 , m 4 , m 5 ) ( - m 1 , - m 2 , - m 3 , - m 4 , m 5 ) .

2.3.2. Equilibria and Stability

The equilibria of system (5) can be calculated by solving the equations D * α j m i = 0 ( j = 1,2 , 3 ; i = 1,2 , , 5 ) , and this system has an isolated equilibria E 0 = ( 0,0 , 0,0 , 0 ) and nontrivial equilibria E θ = ( r cos θ , r sin θ , r cos θ , r sin θ , a 2 ) , where r = a 2 a 3 , θ [ 0,2 π ] .

As to the equilibrium E 0 , it is stable when a 2 < 0 and unstable when a 2 > 0 . For E θ , the characteristic polynomial of Jacobian matrix is λ ( λ + a 1 - a 2 ) ( λ 3 + ( a 1 + a 3 - a 2 ) λ 2 + a 1 a 3 λ + 2 a 1 a 2 a 3 ) = 0 when a 2 > 0 . According to the fractional-order Routh-Hurwitz conditions , when ( a 1 + a 3 - a 2 ) · ( a 1 a 3 ) > 2 a 1 a 2 a 3 , E θ will be stable.

2.3.3. Chaotic Behavior and Attractors

Using the above discretization scheme (6), we find that chaotic behaviors exist in the fractional-order complex Lü system. In the numerical simulations, the system parameters are chosen as ( a 1 , a 2 , a 3 ) = ( 42,22,5 ) , and an initial value is ( x 1 , x 2 , x 3 ) T = ( 1 + 2 j , 3 + 4 j , 5 ) T . When varying the fractional derivative order α i    ( i = 1,2 , 3 ) , system (5) will display diverse motions. The existence of chaos is demonstrated with the time histories, phase diagrams, bifurcation diagrams, and the largest Lyapunov exponents. It is well known that there are many effective algorithms for the calculation of the Lyapunov exponents . In this paper, the largest Lyapunov exponents are calculated by Wolf algorithm .

( 1) Commensurate Order α 1 = α 2 = α 3 = α . The bifurcation diagram is calculated numerically against α [ 0.9,1 ] , while the incremental value of α is 0.0002. From the bifurcation diagram shown in Figure 2, it is found that chaotic range is α [ 0.928,1 ] . To identify the route to chaos, the time history of m 5 is shown in Figures 3(a)3(d). It is clearly shown that the state variables are stable at the fixed point at α = 0.922 , which can be seen in Figure 3(a). When α increases, intermittent dynamical behavior is observed in Figures 3(b)-3(c). As α is further increased, the motion become chaotic as shown for α = 0.928 , where the largest Lyapunov exponent is λ = 0.0154 . In Figures 4(a)-4(b), phase portraits are shown at α = 0.927 and 0.928 , respectively. Numerical evidence displays that the lowest order to yield chaos is 4.64, where α = 0.928 .

Bifurcation diagram of system (5) with α 1 = α 2 = α 3 = α [ 0.9,1 ] .

Time histories showing the rout to chaos via intermittency for system (5) at α 1 = α 2 = α 3 = α .

α = 0.922

α = 0.926

α = 0.927

α = 0.928

Phase portraits of system (5) at α 1 = α 2 = α 3 = α .

α = 0.927

α = 0.928

( 2) α 2 = α 3 = 1 , and Let α 1 Vary. The bifurcation diagram is calculated numerically against α 1 [ 0.78,1 ] , while the incremental value of α 1 is 0.0002. Figure 5 shows that chaotic motions exist in the range α 1 [ 0.812,1 ] . To identify the route to chaos, the time history of m 5 is shown in Figures 6(a)6(d). At α 1 = 0.807 , the state variables are stable at the fixed point as depicted in Figure 6(a). When α 1 increases, intermittent dynamical behavior is observed in Figures 6(b)-6(c). As α 1 is further increased, the motion become chaotic as shown for α 1 = 0.812 , where the largest Lyapunov exponent is λ = 0.0899 . In this case, the lowest order for system (5) to be chaotic is 4.624, where α 1 = 0.812 .

Bifurcation diagram of system (5) with α 1 [ 0.78,1 ] .

Time histories showing the rout to chaos via intermittency for system (5) at α 2 = α 3 = 1 with α 1 varying.

α 1 = 0.807

α 1 = 0.810

α 1 = 0.811

α 1 = 0.812

( 3) α 1 = α 3 = 1 , and Let α 2 Vary. The bifurcation diagram is calculated numerically against α 2 [ 0.75,1 ] , while the incremental value of α is 0.0005. Figure 7(a) shows that the chaotic zone covers most of the range of α 2 [ 0.805,1 ] . To observe the dynamical behavior of system, the region of α 2 [ 0.75,0.85 ] is expanded in step size of 0.0002 as shown in Figure 7(b). The period-doubling bifurcations can be seen in Figure 7(b). Phase diagrams shown in Figures 8(a)8(d) exhibit period-1, period-2, period-4, and chaotic behaviors. Thus, Figure 8 identifies a period doubling route to chaos. In this case, the lowest order for system (5) to be chaotic is 4.61, where α 2 = 0.805 and the largest Lyapunov exponent is λ = 0.0461 .

Bifurcation diagram of system (5) at α 1 = α 3 = 1 with α 2 varying.

α 2 ( 0.75,1 )

α 2 ( 0.75,0.85 )

Phase portraits of system (6) at α 1 = α 3 = 1 with α 2 varying.

α 2 = 0.801

α 2 = 0.802

α 2 = 0.804

α 2 = 0.805

( 4) α 1 = α 2 = 1 , and Let α 3 Vary. The system is calculated numerically against α 3 [ 0.8,1 ] with the step size of 0.001. Figures 9(a)9(d) displays the phase portraits at α 3 = 0.830,0.831,0.90 and 0.950 , respectively. Results show that chaos exists in the range α 3 [ 0.831,1 ] . To identify the route to chaos, the time history of m 5 is shown in Figures 10(a)10(d). At α 1 = 0.820 , the state variables are stable at the fixed point as depicted in Figure 10(a). When α increases, intermittent dynamical behavior is observed in Figures 10(b)-10(c). As α is further increased, the motion become chaotic as shown for α = 0.831 , where the largest Lyapunov exponent is λ = 0.0516 . Numerical evidence displays that the lowest order for system (5) to be chaotic is 4.831, where α 3 = 0.831 .

Phase portraits of system (5) at α 1 = α 2 = 1 with α 3 varying.

α 3 = 0.830

α 3 = 0.831

α 3 = 0.90

α 3 = 0.950

Time histories showing the rout to chaos via intermittency for system (5) at α 1 = α 2 = 1 with α 3 varying.

α 3 = 0.820

α 3 = 0.825

α 3 = 0.830

α 3 = 0.831

3. Antisynchronization between Different Fractional-Order Complex Systems

In this section, we give a general method to achieve antisynchronization of different fractional-order complex systems firstly. Consequently, antisynchronization between fractional-order complex Lü and Lorenz system can be achieved. Without loss of generality, we assume that the derivative order is α    ( α < 1 ) in both master system and slave system.

3.1. A General Method for Antisynchronization of Fractional-Order Complex Systems

Consider the following fractional-order complex system: (8) D * α x = A x + f ( x ) , where x = ( x 1 , x 2 , , x n ) T is the state complex vector, x = x r + j x i , and define x 1 = m 1 + j m 2 , x 2 = m 3 + j m 4 , , x n = m 2 n - 1 + j m 2 n . f = ( f 1 , f 2 , , f n ) T is a vector of nonlinear complex functions and A R n × n is the matrix of system parameters. Superscripts r and i stand for the real and imaginary parts of the state complex vector. System (8) is considered as the master system and the slave system is given by (9) D * α y = B y + g ( y ) + u , where y = ( y 1 , y 2 , , y n ) T is the state complex vector, y = y r + j y i , and define y 1 = s 1 + j s 2 , y 2 = s 3 + j s 4 , , y n = s 2 n - 1 + j s 2 n . u = u r + j u i is designed controller, where u r = ( u 1 , u 3 , , u 2 n - 1 ) T ,   u i = ( u 2 , u 4 , , u 2 n ) T .

Remark 1.

Some fractional-order chaotic complex systems can be described by (8), such as the fractional-order complex Lorenz, Lü, and Chen systems.

Now we give the stability results for linear fractional-order systems.

Lemma 2 (see [<xref ref-type="bibr" rid="B28">28</xref>]).

Autonomous linear system of the fractional-order D * α x = A x , with x ( 0 ) = x 0 is asymptotically stable if and only if | arg ( λ i ( A ) ) | > α π / 2 , ( i = 1,2 , 3 ) . In this case, the component of the state decay towards 0 like t - α . Also, this system is stable if and only if either it is asymptotically stable or those critical eigenvalues which satisfy | arg ( λ i ( A ) ) | = α π / 2 have geometric multiplicity one, where arg ( λ i ( A ) ) denotes the argument of the eigenvalue λ i of A .

Theorem 3.

Antisynchronization between (8) and (9) will be achieved, if the controller is designed as follows: (10) u ( t ) = u r + j u i = - f ( x ) - g ( y ) + ( B - A ) x - K e . The real and imaginary parts of (10) are (11) u r = - f r ( x ) - g r ( y ) + ( B - A ) x r - K e r , u i = - f i ( x ) - g i ( y ) + ( B - A ) x i - K e i ; here K is the control gain matrix, which satisfies | arg ( λ i ( B - K ) ) | > α π / 2 for all the eigenvalues of B - K .

Proof.

From the definition of antisynchronization, we obtain the error vector between (8) and (9) as follows: (12) e ( t ) = y ( t ) + x ( t ) , The derivative of the error vector (12) can be expressed as (13) D * α e ( t ) = D * α y ( t ) + D * α x ( t ) = B y + g ( y ) + u + A x + f ( x ) , Substituting (10) into (13), the error dynamical system (13) can be written as (14) D * α e ( t ) = ( B - K ) e ( t ) . Since | arg ( λ i ( B - K ) ) | > α π / 2 , according to Lemma 2, the error vector e ( t ) asymptotically converges to zero as t . So antisynchronization between different fractional-order complex systems is achieved by using the controller (10). This completes the proof.

Remark 4.

If A = B and f ( · ) = g ( · ) , systems (8) and (9) become identical. Therefore, our scheme is also applicable to achieve antisynchronization of two identical fractional-order chaotic complex systems.

3.2. Antisynchronization between Fractional-Order Complex Lü and Lorenz System

In this section, the antisynchronization behavior between the fractional-order complex Lü and Lorenz systems is made. It is assumed that the fractional-order complex Lü system drives the fractional-order complex Lorenz system . Thus the master system is described by (15) D * α 1 x 1 m = a 1 ( x 2 m - x 1 m ) , D * α 2 x 2 m = - x 1 m x 3 m + a 2 x 2 m , D * α 3 x 3 m = 1 2 ( x ¯ 1 m x 2 m + x 1 m x ¯ 2 m ) - a 3 x 3 m , where (16) A = ( - a 1 a 1 0 0 a 2 0 0 0 - a 3 ) , f ( x ) = ( 0 - x 1 m x 3 m 1 2 ( x ¯ 1 m x 2 m + x 1 m x ¯ 2 m ) ) . The slave system is (17) D * α 1 x 1 s = b 1 ( x 2 s - x 1 s ) + u 1 + j u 2 , D * α 2 x 2 s = b 2 x 1 s - x 2 s - x 1 s x 3 s + u 3 + j u 4 , D * α 3 x 3 s = 1 2 ( x ¯ 1 s x 2 s + x 1 s x ¯ 2 s )    - b 3 x 3 s + u 5 , where (18) B = ( - b 1 b 1 0 b 2 - 1 0 0 0 - b 3 ) , g ( y ) = ( 0 - x 1 s x 3 s 1 2 ( x ¯ 1 s x 2 s + x 1 s x ¯ 2 s ) ) . In the numerical simulations, the initial values of the master and slave systems are ( x 1 m , x 2 m , x 3 m ) T = ( 1 + 2 j , 3 + 4 j , 5 ) T and ( x 1 s , x 2 s , x 3 s ) T = ( - 1 + 9 j , 8 - 5 j , 1 ) T , respectively. Choose the parameters of the master and slave system as ( a 1 , a 2 , a 3 ) = ( 42,22,5 ) , ( b 1 , b 2 , b 3 ) = ( 10,180,1 ) , α 1 = α 2 = α 3 = 0.95 . In order to satisfy | arg ( λ i ( B - K ) ) |    > α π / 2 , we choose the gain control matrix K as follows: (19) K = ( 0 0 0 180 0 0 0 0 0 ) .

According to Theorem 3, the controller is designed as (20) u = - ( 0 - x 1 m x 3 m 1 2 ( x ¯ 1 m x 2 m + x 1 m x ¯ 2 m ) ) - ( 0 - x 1 s x 3 s 1 2 ( x ¯ 1 s x 2 s + x 1 s x ¯ 2 s ) ) + ( B - A ) ( x 1 m x 2 m x 3 m ) - K ( x 1 s + x 1 m x 2 s + x 2 m x 3 s + x 3 m ) = ( ( b 1 - a 1 ) ( x 2 m - x 1 m ) x 1 m x 3 m + x 1 s x 3 s + b 2 x 1 m - ( a 2 + 1 ) x 2 m - 180 ( x 1 s + x 1 m ) - 1 2 ( x ¯ 1 m x 2 m + x 1 m x ¯ 2 m ) - 1 2 ( x ¯ 1 s x 2 s + x 1 s x ¯ 2 s ) + ( a 3 - b 3 ) x 3 m ) .

The errors of antisynchronization converge asymptotically to zero in a quite short period as depicted in Figure 11. Figure 12 shows state variables of drive system and response system, in which the state complex variables are demonstrated by real and imaginary part, respectively. The above results verify that antisynchronization between fractional-order complex Lü system and Lorenz system has been achieved.

The time evolution of the synchronization errors between systems (15) and (17) (time/s).

State variables of master system (15) and slave (17) (time/s).

4. Conclusions

In this paper, a new fractional-order chaotic complex system is proposed. By means of phase portraits, bifurcation diagrams, the histories, and the largest Lyapunov exponents, we investigate chaotic behavior of this new system. Our results show that the new system displays many interesting dynamical behaviors, such as fixed points, periodic motions, and chaotic motions. Two typical routes to chaos—period doubling and intermittency—are found in this system. Besides, when the parameters of the system are fixed, the lowest order for chaos to exist is determined. Moreover, antisynchronization of different fractional-order chaotic complex systems has been studied. Meanwhile, the new system and the fractional-order complex Lorenz system can achieve antisynchronization.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The research is supported by the National Nature Science Foundation of China (nos. 61273088 and 10971120) and the Nature Science Foundation of Shandong province (no. ZR2010FM010). The authors would like to thank the editors and anonymous referees for their constructive comments and suggestions.

Hartley T. T. Lorenzo C. F. Qammer H. K. Chaos in a fractional order Chua's system IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications 1995 42 8 485 490 10.1109/81.404062 2-s2.0-0029359880 Li C. Chen G. Chaos and hyperchaos in the fractional-order Rössler equations Physica A: Statistical Mechanics and its Applications 2004 341 1–4 55 61 10.1016/j.physa.2004.04.113 MR2092674 2-s2.0-3342927052 Grigorenko I. Grigorenko E. Chaotic dynamics of the fractional Lorenz system Physical Review Letters 2003 91 3 034101 2-s2.0-0041384356 Li C. Chen G. Chaos in the fractional order Chen system and its control Chaos, Solitons and Fractals 2004 22 3 549 554 10.1016/j.chaos.2004.02.035 ZBL1069.37025 2-s2.0-2042470796 Lu J. G. Chaotic dynamics of the fractional-order Lü system and its synchronization Physics Letters A 2006 354 4 305 311 10.1016/j.physleta.2006.01.068 2-s2.0-33646366905 Pan L. Zhou W. Fang J. Li D. Synchronization and anti-synchronization of new uncertain fractional-order modified unified chaotic systems via novel active pinning control Communications in Nonlinear Science and Numerical Simulation 2010 15 12 3754 3762 10.1016/j.cnsns.2010.01.025 ZBL1222.34063 2-s2.0-77952883806 Si G. Sun Z. Zhang Y. Chen W. Projective synchronization of different fractional-order chaotic systems with non-identical orders Nonlinear Analysis: Real World Applications 2012 13 4 1761 1771 10.1016/j.nonrwa.2011.12.006 MR2891007 2-s2.0-84856545857 Chen L. Chai Y. Wu R. Lag projective synchronization in fractional-order chaotic (hyperchaotic) systems Physics Letters A: General, Atomic and Solid State Physics 2011 375 21 2099 2110 10.1016/j.physleta.2011.04.015 2-s2.0-79955688348 Roldán E. de Valcárcel G. J. Vilaseca R. Mandel P. Single-mode-laser phase dynamics Physical Review A 1993 48 1 591 598 10.1103/PhysRevA.48.591 2-s2.0-0000367731 Ning C.-Z. Haken H. Detuned lasers and the complex Lorenz equations: subcritical and supercritical Hopf bifurcations Physical Review A 1990 41 7 3826 3837 10.1103/PhysRevA.41.3826 2-s2.0-0013044957 Toronov V. Y. Derbov V. L. Boundedness of attractors in the complex Lorenz model Physical Review E 1997 55 3 3689 3692 10.1103/PhysRevE.55.3689 MR1438730 2-s2.0-0000288953 Mahmoud G. M. Bountis T. The dynamics of systems of complex nonlinear oscillators: a review International Journal of Bifurcation and Chaos in Applied Sciences and Engineering 2004 14 11 3821 3846 10.1142/S0218127404011624 MR2114882 ZBL1091.34524 2-s2.0-17044382505 Mahmoud E. E. Complex complete synchronization of two nonidentical hyperchaotic complex nonlinear systems Mathematical Methods in the Applied Sciences 2014 37 3 321 328 10.1002/mma.2793 MR3158634 Wu Z. Duan J. Fu X. Complex projective synchronization in coupled chaotic complex dynamical systems Nonlinear Dynamics 2012 69 3 771 779 10.1007/s11071-011-0303-0 MR2943344 2-s2.0-84863983340 Zhang F. Liu S. Full state hybrid projective synchronization and parameters identification for uncertain chaotic (hyperchaotic) complex systems Journal of Computational and Nonlinear Dynamics 2014 9 2 021009 10.1115/1.4025475 2-s2.0-84887939971 Mahmoud G. M. Mahmoud E. E. Complex modified projective synchronization of two chaotic complex nonlinear systems Nonlinear Dynamics 2013 73 4 2231 2240 10.1007/s11071-013-0937-1 MR3094790 2-s2.0-84882449836 Liu S. Zhang F. Complex function projective synchronization of complex chaotic system and its applications in secure communication Nonlinear Dynamics 2014 76 2 1087 1097 10.1007/s11071-013-1192-1 MR3192201 2-s2.0-84899632205 Luo C. Wang X. Chaos in the fractional-order complex Lorenz system and its synchronization Nonlinear Dynamics 2013 71 1-2 241 257 10.1007/s11071-012-0656-z MR3010577 2-s2.0-84867977717 Luo C. Wang X. Chaos generated from the fractional-order complex Chen system and its application to digital secure communication International Journal of Modern Physics C 2013 24 23 1350025 10.1142/S0129183113500253 MR3047663 2-s2.0-84876258326 Oldham K. B. Spanier J. The Fractional Calculus 1974 San Diego, Calif, USA Academic Press Podlubny I. Fractional Differential Equations 1999 New York, NY, USA Academic Press Mahmoud G. M. Bountis T. Mahmoud E. E. Active control and global synchronization of the complex Chen and Lü systems International Journal of Bifurcation and Chaos 2007 17 12 4295 4308 10.1142/S0218127407019962 MR2394229 2-s2.0-40549096015 Diethelm K. Ford N. J. Freed A. D. A predictor-corrector approach for the numerical solution of fractional differential equations Nonlinear Dynamics 2002 29 1–4 3 22 10.1023/A:1016592219341 2-s2.0-0036650479 Wang X. He Y. Wang M. Chaos control of a fractional order modified coupled dynamos system Nonlinear Analysis: Theory, Methods &; Applications 2009 71 12 6126 6134 10.1016/j.na.2009.06.065 MR2566518 2-s2.0-72149089014 Wolf A. Swift J. B. Swinney H. L. Vastano J. A. Determining Lyapunov exponents from a time series Physica D: Nonlinear Phenomena 1985 16 3 285 317 10.1016/0167-2789(85)90011-9 2-s2.0-0008494528 Pardalos P. M. Yatsenko V. A. Optimization approach to the estimation and control of Lyapunov exponents Journal of Optimization Theory and Its Applications 2006 128 1 29 48 10.1007/s10957-005-7554-1 MR2201888 2-s2.0-31144456498 Nair S. P. Shiau D.-S. Principe J. C. Iasemidis L. D. Pardalos P. M. Norman W. M. Carney P. R. Kelly K. M. Sackellares J. C. An investigation of EEG dynamics in an animal model of temporal lobe epilepsy using the maximum Lyapunov exponent Experimental Neurology 2009 216 1 115 121 10.1016/j.expneurol.2008.11.009 2-s2.0-59749085335 Matignon D. Stability results for fractional differential equations with applications to control processing Proceedings of the International IMACS IEEE-SMC Multiconference on Computational Engineering in Systems Applications 1996 Lille, France