#### Abstract

We study the stability and stabilization of complex fractional Lorenz system. The fractional calculus are taken in sense of the Caputo derivatives. The technique is based on stability theory of fractional-order systems. Numerical solutions are imposed.

#### 1. Introduction

Chaos theory is a new branch of mathematics and physics. It provided a new method to understand the behavior of the processes in the world. Chaotic behaviors virtue in different areas of science and engineering such as mechanics, physics, electronics, medicine, biology, ecology, economy. In 1963, Lorenz discovered the chaotic system [1]. A dynamical system is called chaotic if it is deterministic, has long-term aperiodic behavior, and exhibits sensitive dependence on the initial conditions. If the system has one positive Lyapunov exponent, then it is called chaotic. The Lorenz system is a system of ordinary differential equations that imposing chaotic solutions for certain parameter values and initial conditions. Researches along this line ultimately led to the ruling of different classes of relevant chaotic systems. The Lorenz system has advantages, it presents a symmetry, and the equations are invariant under . Also the system is dissipative with a value of the divergence and therefore the volume in phase space always contracts under the flow. Moreover the Lorenz system is bounded.

This system was widely considered over the course of the 1970s, 1980s, and 1990s (see [2–5]) as a paradigmatic chaotic problem. In 1994, Sprott defined the diffusionless Lorenz system (a one-parameter version of Lorenz system) [6]. In 1999, Chen and Ueta established an interesting system [7], which is the dual system to the Lorenz system [1], in the sense defined by Vaněček and Čelikovský [8, 9]. Later, Lü and Chen [10] furthermore introduced a chaotic system, which refers to the transition between the Lorenz and the Chen systems. In 2005, Čelikovský and Chen described another generalized Lorenz canonical form [11]. More systems are posed in 2006-2007; a unified Lorenz-type system is recognized by Yang et al. [12, 13]. In 2009, Huang and Yang found a new chaotic system, which contains special cases as the modified Lorenz system and conjugate Chen system [14].

Recently, the Lorenz system has got many studies and attentions from researchers. Zhang et al. considered the synchronization of two identical hyperchaotic Lorenz systems with time delay only by a single linear controller [15], while Shi and Wang investigated a class of new synchronization (hybrid synchronization) coexisting in the same system [16]. Camargo et al. considered a system of two identical linearly coupled Lorenz oscillators presenting synchronization of chaotic motion for a specified range of the coupling strength. They verified the existence of global synchronization and antisynchronization attractors with intermingled basins of attraction, such that the basin of one attractor is riddled with holes belonging to the basin of the other attractor and vice versa [17]. Finally, Li et al. concerned with model-free control of the Lorenz chaotic system, where only the online input and output are available, while the mathematic model of the system is unknown [18].

Fractional-order chaotic systems have been studied widely in current years in many fields such as chaotic phenomena, chaotic control, chaotic synchronization, and other related studies. It has been provided that many fractional-order dynamical systems, as some well-known integer-order systems, can be displayed complex bifurcation and chaotic phenomena. For example, the fractional-order Lorenz system, the fractional-order Chen system, the fractional-order Lü system, and the fractional-order unified system also exhibit chaotic behavior.

Čelikovský and Chen (2002) generalized a new Lorenz canonical form of one parameter [19, 20]. The dynamics of the fractional-order Lorenz model was also investigated by I. Grigorenko and E. Grigorenko (2003) [21]. Li and Peng (2004) studied the time-domain algorithm to compute the chaotic attractor of the fractional differential system, which is completely different from the frequency-domain one [22]. Yan and Li (2007) [23] considered chaos synchronization of fractional-order Lorenz, Rossler, and Chen systems taking one as master and second one as slave. Zhou and Cheng (2008) [24] discussed a synchronization between different fractional-order chaotic systems, namely, Rossler and Chen systems and Chua and Chen systems. Based on the qualitative theory, the existence and uniqueness of solutions for a class of fractional-order Lorenz chaotic systems have been investigated theoretically by Yu et al. (2009) [25]. Fractional-order diffusion less Lorenz system has been investigated numerically by Sun and Sprott (2009) [26], and results have shown that the system has complex dynamics with interesting characteristics. Recently, Agrawal et al., synchronized fractional-order chaotic systems using active control method [27]. Xu et al., investigated the phenomenon of diffusionless Lorenz system with fractional-order [28]. Zhou and Ding studied the stability and stabilizations of equilibrium points of the fractional-order Lorenz chaotic system [29]. Finally, Si et al. employed an uncertain complex network with four fractional-order Lorenz systems to verify the effectiveness of the proposed approach [30].

In current years, researchers introduced and studied several types of chaotic nonlinear systems with complex variables. These systems which involving complex variables are employed to describe the physics of a detuned laser, rotating fluids, disk dynamos, electronic circuits, and particle beam dynamics in high-energy accelerators. As special model, the chaotic complex Lorenz system which is used to describe and simulate the physics of detuned lasers and thermal convection of liquid flows. The complex Lorenz model is, in fact, the outcome of the study of Baroclinic models. This model corresponds to the equilibrium state of the atmosphere in which surfaces of constant density are not parallel to the surface of constant gravitational potential. In 1982, Fowler et al. defined and studied a complex Lorenz as a generalization of the real Lorenz system [31, 32]. The complex variables of Lorenz system are related, respectively, to electric field, the atomic polarization amplitudes, and the population inversion in a ring laser system of two-level atoms; for more details, see [33, 34]. Later Rauh et al. (1996) [35] introduced a set of equations to model such a system. This has the same structure as that of Lorenz equations, but it occurs in the complex domain. Panchev and Vitanov (2005) [36] studied asymptotic properties of some complex Lorenz-like systems, while other properties and chaotic synchronization of the complex Lorenz model are studied by Mahmoud et al. (2007) [37]. Further studies are introduced due to Mahmoud et al.; they constructed the new hyperchaotic complex Lorenz systems by extending the ideas of adding state feedback controls and introducing the complex periodic forcing [38, 39], while G. M. Mahmoud and E. E. Mahmoud (2010) [40] investigated complete synchronization of n-dimensional chaotic complex systems with uncertain parameters. Li et al. (2012) [41] considered the distributed robust control problems of uncertain linear multiagent systems with undirected communication topologies. It was assumed that the agents have identical nominal dynamics while subject to different norm-bounded parameter uncertainties, leading to weakly heterogeneous multiagent systems. Mahmoud (2012) [42] introduced a new hyperchaotic complex Lorenz system. This hyperchaotic complex system was constructed by adding a linear controller to the second equation of the chaotic complex Lorenz system. Finally, El-Sayed et al. [43] studied the dynamic properties of the continuous dynamical system of the fractional equations of complex variables. The five-dimensional complex Lorenz model is often used to explain and simulate the physics of tuned lasers. These studies have their origin in the shaping expression by Haken [44] of the isomorphism between the three equations of the real Lorenz model and the three real equations for a single-mode laser operating with its resonant cavity tuned to resonance with the material transition. Furthermore, analytic studies of model systems can be found in [45–47].

In this paper, we study the stability and stabilization of fractional complex Lorenz model. The fractional calculus are taken in sense of the Caputo derivatives. Caputo initial value problem holds for both homogeneous and nonhomogeneous conditions. For this reason choice Caputo derivative is better than other fractional derivatives. The technique based on stability theory of fractional-order systems. Numerical solutions are imposed.

#### 2. Fractional Calculus

The idea of the fractional calculus (i.e., calculus of integrals and derivatives of any arbitrary real or complex order) was planted over 300 years ago. Abel in 1823 investigated the generalized tautochrone problem and for the first time applied fractional calculus techniques in a physical problem. Later Liouville applied fractional calculus to problems in potential theory. Since that time the fractional calculus has haggard the attention of many researchers in all area of sciences. Recently, the theory of fractional calculus has found interesting applications in the theory of analytic functions. The classical definitions of fractional operators and their generalizations have fruitfully been applied in obtaining, for example, the characterization properties, coefficient estimates, distortion inequalities, and convolution structures for various subclasses of analytic functions and the works in the research monographs [48–50]. This section concerns with some preliminaries and notations regarding the fractional calculus. Let us start with the Riemann-Liouville fractional operators.

*Definition 1. *The fractional (arbitrary) order integral of the function of order is defined by
When , we write , where denoted the convolution product, , and and , as where is the delta function.

*Definition 2. *The fractional (arbitrary) order derivative of the function of order is defined by

*Remark 3. *From Definitions 1 and 2, , we have

*Definition 4. *The Caputo fractional derivative of order is defined, for a smooth function by
where , (the notation stands for the largest integer not greater than .

*Remark 5. *The following relations hold the following:(i)Representation
(ii)The Caputo fractional derivative of the power function
(iii)(iv)Linearity
(v)Noncommutation

#### 3. Stability of Fractional Complex System

A complex Lorenz-like system has been found in Laser Physics and while analyzing baroclinic instability of the geophysical flows in the atmosphere which takes the form where and are complex variables; that is, and , but is real variable, and , such that , and are real parameters. It assumed that is a function in .

Here, we consider the following fractional complex Lorenz system, in sense of the Caputo derivative: where ; in this note we shall take . For real variables, system (13) reduces to the system of the form

In this section, we shall discuss the stability of the system (14). For this purpose we need the following preliminaries which can be found in [51].

*Definition 6. *The zero solution of the equation is said to be stable if, for any initial values , there exists such that . The zero is said to be asymptotically stable if it is stable, and .

Lemma 7. *Assume the system of the form
**
where . If , denotes the eigenvalues of , and denotes the norm; and , then the system (15) is asymptotically stable.*

Lemma 8. *Assume the controller system with the linear feedback control input
**
where is a feedback, , . If , and , then the system (16) is asymptotically stable.*

Theorem 9. *If , then the system (14) is asymptotically stable at the equilibrium point .*

* Proof. *System (14) can be written in the form (15), where
It is obvious that satisfies
Moreover the characteristic equation of the system is
and when , the eigenvalues are
where are of multiplicity 2. Since , therefore
and for suitable choice of the parameters , and we have . According to Lemma 7, it implies that the equilibrium point of system (14) is asymptotically stable. Now for and , the eigenvalues are
hence relation (16) holds, and consequently the system (14) is asymptotically stable. Finally, the case and such that , and implies the eigenvalues of multiplicity 2
which yields that the system is asymptotically stable at . Note that the other cases impose complex . This completes the proof.

The system (14) can be assumed in the form (16), we derive the following result.

Theorem 10. *If , then the controlled system (14) is asymptotically stable at the equilibrium point .*

* Proof. *Applying control input on (14), where and are taken in which , where the characteristic equation of the system is
Obviously, for the eigenvalues of the system are
where and of multiplicity 2. Also when , the eigenvalues are
Finally for , and, the system admits the following eigenvalues:
Furthermore, ; hence in view of Lemma 8, controlled system (14) is asymptotically stable at .

#### 4. Stabilizing

In this section, we design an evaluation controller for fractional-order Lorenz chaotic system (14) via fractional-order derivative. For this puropose we need the following result which can be found in [52].

Lemma 11. *The fixed points of the following nonlinear commensurate fractional-order autonomous system,
**
are asymptotically stable if all eigenvalues of the Jacobian matrix evaluated at the fixed points satisfy , where are continuous nonlinear vector functions, and the fixed points of this nonlinear commensurate fractional-order system are calculated by solving equation .*

We can obtain the following results.

Theorem 12. *Assume the controlled fractional-order Lorenz chaotic system
**
where is the fractional-order controller, and is the feedback coefficient. For the case , if ,
**
then system (30) will asymptotically converge to the unstable equilibrium point . Furthermore, for the case , and , if and , then system (30) will asymptotically converge to the unstable equilibrium point .*

* Proof. *The Jacobi matrix of the controlled fractional-order Lorenz chaotic system (30) at is
Since , we have
Thus the characteristic equation takes the form
which leads to the eigenvalues
when . And for , we have
Since , and , therefore,
where . According to Lemma 11, it yields that the equilibrium point of system (30) is asymptotically stable; that is, the unstable equilibrium point in fractional-order Lorenz system (30) can be stabilized via fractional-order derivative. For the case , and , since , then the characteristic equation
has negative eigenvalues, and hence the system is asymptotically stable. The proof is completed.

Theorem 13. *Assume the controlled fractional-order Lorenz chaotic system
**
where is the fractional-order controller, and is the feedback coefficient. For the case , if ,
**
then system (39) will asymptotically converge to the unstable equilibrium point . Furthermore, for the case , and , if and , then system (39) will asymptotically converge to the unstable equilibrium point .*

* Proof. *The Jacobi matrix of the controlled fractional-order Lorenz chaotic system (39) at is

Since , we have

Therefore, the characteristic equation is
which yields the eigenvalues
when . while for , we impose
Since , and , therefore,
where . According to Lemma 11, it yields that the equilibrium point of system (39) is asymptotically stable; that is, the unstable equilibrium point in fractional-order Lorenz system (39) can be stabilized via fractional-order derivative. For the case , since , then the characteristic equation
has negative eigenvalues, and hence the system is asymptotically stable. This completes the proof.

Theorem 14. *Assume the controlled fractional-order Lorenz chaotic system
**
where is the fractional-order controller, and is the feedback coefficient. If , and , then system (48) will asymptotically converge to the unstable equilibrium point .*

* Proof. *The Jacobi matrix of the controlled fractional-order Lorenz chaotic system (48) at is
Because , so we have
The characteristic equation is
which implies the eigenvalues
By the assumptions of the theorem and in virtue of Lemma 11, we obtain that system (48) will asymptotically converge to the unstable equilibrium point .

In the same manner of Theorem 14, we have the following result.

Theorem 15. *Assume the controlled fractional-order Lorenz chaotic system
**
where is the fractional-order controller, and is the feedback coefficient. If , and , then system (53) will asymptotically converge to the unstable equilibrium point .*

#### 5. Numerical Solution

In this section, we discuss the numerical solution of fractional differential equations. All the numerical simulations of fractional-order system in this paper are based on [22]. Set with the initial condition , therefore, the system (14) can be described as follows: where and for The error of this approximation can be computed as follows:

#### 6. Synchronizing

we consider a feedback controller for the fractional-order Lorenz chaotic system (14) via fractional-order derivative and obtain the controlled response system where is the fractional-order controller, and is the feedback coefficient. We need the following result which can be found in [53].

Lemma 16. *The following linear commensurate fractional-order autonomous system,
**
is asymptotically stable if and only if is satisfied for all eigenvalues () of matrix . Also, this system is stable if and only if is satisfied for all eigenvalues () of matrix , and those critical eigenvalues which satisfy have geometric multiplicity one, where , and .*

We have the following theorem.

Theorem 17. *If , and , then the fractional-order Lorenz chaotic system (14) and the controlled fractional-order Lorenz chaotic system (58) achieved synchronization via fractional-order derivative.*

* Proof. *Define the synchronization error variables as follows:
Therefore, we obtain the system
where
Consider and with , such that
Thus, the characteristic equation is of the form
Because yields that . According to Lemma 16, it implies that the equilibrium point of error system (60) is asymptotically stable,and hence the fractional-order Lorenz chaotic system (14) and the controlled fractional-order Lorenz chaotic system (58) achieved synchronization via fractional-order derivative. The proof is completed.

#### 7. Conclusion and Discussion

The five-dimensional complex models are habitually used to clarify and simulate the physics of tuned lasers. By using fractional-order derivative (in sense of the Caputo derivatives), we stabilized the unstable equilibrium points of the complex fractional-order Lorenz chaotic system and included chaos synchronization for the fractional-order Lorenz chaotic system. Numerical solutions are computed for these systems. Conditions are imposed to study the stability, and several cases are discussed for these systems. Figure 1 shows the asymptotically stable at , such that , and the values of are taken positive in (a) with the initial values , negative in (b) with the initial values , and for in (c), while in (d) the parameters are valued as with the initial values . Figures 2 and 3 confirm the time evolution at . Figures 4–6 show the affectedness of , where in Figure 4, in Figure 5, and in Figure 6. Figure 7 shows the phase plane of . Here we discuss the system for . Note that and are congruent on and , respectively. Finally, Figure 8 represents to the chaotic attractors of the proposed system with the parameters and with the initial values , respectively.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(d)**