# A Novel Control Method for Integer Orders Chaos Systems via Fractional-Order Derivative

Authors: Ping Zhou; Fei Kuang
Journal: Discrete Dynamics in Nature and Society
 (2011)

---

## Abstract

A fractional-order control method is obtained to stabilize the point in chaos attractor of integer orders chaos systems. The control law has simple structure and is designed easily. Two examples are also given to illustrate the effectiveness of the theoretical result.


---

## Body

## 1. Introduction

Chaotic phenomena has been observed in many areas of science and engineering such as mechanics, electronics, physics, medicine, ecology, biology, and economy. To avoid troubles arising from unusual behaviors of a chaotic system, chaos control has received a great deal of interest among scientists from various research fields in the past few decades [1–8]. On the other hand, derivatives and integrals of fractional order have been found in many applications in recent years in physics and engineering. Many systems [4–13] are known to display fractional-order dynamics, such as viscoelastic systems, dielectric polarization, electrode-electrolyte polarization, and electromagnetic waves.In the recent years, emergence of effective methods in the differentiation and integration of non integer order equations makes fractional-order systems more and more attractive for the systems control community. It is verified that the fractional-order controllers can have better disturbance rejection ratios and less sensitivity to plant parameter variations compared to the traditional controllers [14]. However, few results on control the saddle point of index 2 in integer orders chaos system are presented via fractional-order derivative. A fractional-order control method [14] is presented to stabilize the unstable equilibrium of integer orders chaos systems, but the method cannot apply to the desired point in chaotic attractor. Motivated by that, in this paper, a very simple control method is presented for a class of integer orders chaotic systems via fractional-order derivative. The order of fractional-order derivative is only determined by the eigenvalues at the desired point in integer orders chaotic attractor. The control technique, based on stability theory of fractional-order systems, is simple and theoretically rigorous.


## 2. Control Chaos via Fractional-Order Derivative

Consider the following 3D integer orders chaos system:(2.1)Ẋ=F(X),
where X=(x1,x2,x3)T are real state variables and F:R3→R3 is differentiable function.LetX0=(x10x20x30)T be the any point in the chaos attractor of integer orders system (2.1), A0 the Jacobian matrix of F(X) at point X0, and λ±=a±jb, λ3=-c the eigenvalues of A0, respectively. In this paper, we always assume that a>0, c>0, and b≠0, especially X0  is called the saddle point of index  2 if it also is the equilibrium point of system (2.1). In order to stabilize the point X0  via fractional-order derivative, we introduce the following controller:(2.2)V(t)=-F(X0)-dqXdtq+Ẋ,
where 0<q<1 is the fractional order and will be determined later. Then, we can construct the control system for (2.1) as follows:(2.3)Ẋ=F(X)+V(t).Hence, the problem of stabilization ofX0 in chaos attractor of system (2.1) is shifted into the stability of X0 in system (2.3). So, the stabilization problem of X0 in system (2.1) can be solved if q is suitably designed such that X0 in system (2.3) is asymptotically stable.Now, we state the main result in this paper as follows.Theorem 2.1.
If0<q<q0=Δ(2/π)arctan|b/a| in V(t), then the point X0 is asymptotically stable in control system (2.3); that is, the point X0in chaos system (2.1) can be stabilized via fractional-order derivative.Proof.
According to system (2.3) and controller (2.2), we can obtain that
(2.4)dqXdtq=F(X)-F(X0).
It can be seen that the pointX0  is one equilibrium point of system (2.4) and the Jacobian matrix at this point for system (2.4) is the same as the Jacobian matrix at this point for system (2.1); therefore, the eigenvalues of the Jacobian matrix of (2.4) at point X0 are λ±=a±jb, λ3=-c.
Forλ3=-c<0, we have
(2.5)arg(λ3)=π>qπ2.
For λ±=a±jb, we also derive that
(2.6)arg(λ±)=arctan|ba|=q0π2>qπ2.
In term of (2.5) and (2.6), for any eigenvalue λ of A0, we obtain that
(2.7)|argλ|>πq2,
which implies that the equilibrium point X0 of system (2.4) is asymptotically stable [15, 16]; that is, the point X0 in the chaos attractor of integer orders chaos system (2.1) can be stabilized via fractional-order derivative; the proof is completed.Remark 2.2.
Especially, ifX0 is the saddle point of index 2 of system (2.1), then the control law obtained in this paper is similar with the controller obtained in [14]. So, Theorem 2.1 proposed by this paper extends the controller proposed by [14].Remark 2.3.
By Theorem2.1 proposed in this paper, we know that it is possible to stabilize the closed-loop system without applying any changes on the places of the poles.Remark 2.4.
IfX0 is the same saddle point for hyperchaotic systems, which implies that the eigenvalues of A0 are λ±=a±jb(a>0,b≠0), λ3=-c1(c1>0), and λ4=-c2(c2>0), respectively, then this point X0 in hyperchaotic systems can be stabilized via the proposed controller.


## 3. Applications

Now, we take Lü chaos system [17] and hyperchaotic Chen system [18] for numerical simulation, respectively. All the numerical simulation of fractional-order system in this paper is based on [19]. It is a direct time-domain approximation numerical simulation. Reference [20] has shown that using frequency-domain approximation in the numerical simulations of fractional systems may result in wrong consequences. This mistake has occurred in the recent literature that found the lowest-order chaotic systems among fractional-order systems.We introduce the numerical solution of fractional differential equations in [19]. There have been several definitions of fractional derivatives. In the following, we introduce the most common one of them(3.1)dqf(t)dtq=Jm-qf(m)(t),q>0,
where mis the first integer which is not less than q, f(m)(t) is the m-order derivative in the usual sense, and Jq(q>0) is the q-order Riemann-Liouville integral operator with the expression(3.2)Jqf(t)=1Γ(q)∫0t(t-τ)q-1f(τ)dτ.
Here, Γ(·) stands for gamma function, and dq/dtq is generally called “q-order Caputo differential operator.”Consider the following fractional-order system:(3.3)dq1xdtq1=f(x,y),dq2ydtq2=g(x,y),0<q1,q2<1,
with initial condition (x0,y0). Now, set h=T/N, tn=nh(n=0,1,2…,N). The above system can be discretized as follows:(3.4)xn+1=x0+hq1Γ(q1+2)[f(xn+1p,yn+1p)+∑j=0nα1,j,n+1f(xj,yj)],yn+1=y0+hq2Γ(q2+2)[g(xn+1p,yn+1p)+∑j=0nα2,j,n+1g(xj,yj)],
where(3.5)xn+1p=x0+1Γ(q1)∑j=0nβ1,j,n+1f(xj,yj),yn+1p=y0+1Γ(q2)∑j=0nβ2,j,n+1g(xj,yj),
and for i=1,2,(3.6)αi,j,n+1={nqi+1-(n-qi)(n+1)qi,j=0,(n-j+2)qi+1+(n-j)qi+1-2(n-j+1)qi+1,1≤j≤n,1,j=n+1,(3.7)βi,j,n+1=hqiqi[(n-j+1)qi-(n-j)qi],0≤j≤n.The error of this approximation is described as follows:(3.8)|x(tn)-xn|=o(hp1),p1=min⁡(2,1+q1),|y(tn)-yn|=o(hp2),p2=min⁡(2,1+q2).Now, we choose Lü chaos system for numerical simulation. The Lü chaos system [17] is(3.9)dx1dt=36(x2-x1),dx2dt=-x1x3+20x2,dx3dt=x1x2-3x3.Now, we choose(x10,x20,x30)=(10,9,15) in chaos attractor of chaos system (3.9). Since the eigenvalues of the Jacobian matrix at X0 for system (3.9) are λ±=4.5554±14.2605j, λ3=-28.1108, we have q0=Δ(2/π)arctan|b/a|=0.8032. According to Theorem 2.1, in order to stabilize this point, we choose q<q0 for controller (2.2).For example, the simulation results are shown in Figure1 for q=0.7, and the initial conditions are (15,14,20)T, where e1=x1-x10, e2=x2-x20, and e3=x3-x30.The control simulation result for (10, 9, 15) whenq=0.7.
(a)(b)(c)Finally, we choose hyperchaotic Chen system for numerical simulation. The hyperchaotic Chen system [18] is(3.10)dx1dt=35(x2-x1)+x4,dx2dt=7x1-x1x3+12x2,dx3dt=x1x2-3x3,dx4dt=x2x3+0.5x4.Now, we choose(x10,x20,x30,x40)=(0,0,19,0) in chaos attractor of the hyperchaotic Chen system (3.10). Since the eigenvalues of the Jacobian matrix at X0 for system (3.10) are λ±=0.4537±3.0878j,λ3=-23.4074, and λ4=-3, respectively, we have q0=Δ(2/π)arctan|b/a|=0.9071. According to Theorem 2.1, in order to stabilize this point, we choose q<q0 for controller (2.2).For example, the simulation results are shown in Figure2 for q=0.8 and the initial conditions are (10,14,49,12)T, where e1=x1-x10, e2=x2-x20, e3=x3-x30, and e4=x4-x40.The control simulation result for (0, 0, 19, 0) whenq=0.8.
(a)(b)(c)(d)


## 4. Conclusion

We construct a novel control law for integer orders chaos system via fractional-order derivative. Any desired point in chaos attractor of integer orders chaos system can be stabilized via the fractional-order derivative. The order of fractional-order derivative of chaos system is only determined by the eigenvalues at the point. By comparison with the traditional controllers, we know that this controller without causing any changes in the eigenvalues of the system at the desired fixed points. The proposed controller is employed to control Lü chaotic system and hyperchaotic Chen system, and some numerical simulation results are obtained. Theoretical analysis and simulation results show that the control method in this paper is effective.



---

*Source: 217843-2011-05-05.xml*