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.
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:Ẋ=F(X),
where X=(x1,x2,x3)T are real state variables and F:R3→R3 is differentiable function.
Let X0=(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:V(t)=-F(X0)-dqXdtq+Ẋ,
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:Ẋ=F(X)+V(t).
Hence, the problem of stabilization of X0 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.
If 0<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
dqXdtq=F(X)-F(X0).
It can be seen that the point X0 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
arg(λ3)=π>qπ2.
For λ±=a±jb, we also derive that
arg(λ±)=arctan|ba|=q0π2>qπ2.
In term of (2.5) and (2.6), for any eigenvalue λ of A0, we obtain that
|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, if X0 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 Theorem 2.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.
If X0 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 themdqf(t)dtq=Jm-qf(m)(t),q>0,
where m is 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 expressionJqf(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: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: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)],
wherexn+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,α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,βi,j,n+1=hqiqi[(n-j+1)qi-(n-j)qi],0≤j≤n.
The error of this approximation is described as follows: |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] isdx1dt=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 Figure 1 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) when q=0.7.
Finally, we choose hyperchaotic Chen system for numerical simulation. The hyperchaotic Chen system [18] isdx1dt=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 Figure 2 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) when q=0.8.
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.
Acknowledgments
The work is supported jointly by National Natural Science Foundation of China under Grant no. 61004042 and Natural Science Foundation Project of CQ CSTC 2009BB2417.
PecoraL. M.CarrollT. L.Synchronization in chaotic systems1990648821824103826310.1103/PhysRevLett.64.821ZBL0938.37019OttE.GrebogiC.YorkeJ. A.Controlling chaos1990641111961199104152310.1103/PhysRevLett.64.1196ZBL0964.37501ChenG.YuX.2003292Berlin, GermanySpringerx+369Lecture Notes in Control and Information Sciences2014603ZBL1167.37335DengW.Generating 3-D scroll grid attractors of fractional differential systems via stair function2007171139653983238439310.1142/S0218127407019810DengW.dengwh@lzu.edu.cnLüJ.jhlu@iss.ac.cnGenerating multi-directional multi-scroll chaotic attractors via a fractional differential hysteresis system20073695-643844310.1016/j.physleta.2007.04.112AhmadW. M.wajdi@sharjah.ac.aeSprottJ. C.sprott@physics.wisc.eduChaos in fractional-order autonomous nonlinear systems200316233935110.1016/S0960-0779(02)00438-1ZBL1033.37019TikhonovD. A.MalchowH.malchow@uos.deChaos and fractals in fish school motion, II200316228728910.1016/S0960-0779(02)00319-3ZBL1032.92034LiC.cgli@uestc.edu.cnChenG.gchen@ee.cityu.edu.hkChaos in the fractional order Chen system and its control200422354955410.1016/j.chaos.2004.02.035ZBL1069.37025LiC. G.LiaoX. F.YuJ . B.Synchronization of fractional order chaotic systems2003686306720310.1103/PhysRevE.68.067203GaoX.gxin7993@163.comYuJ.Synchronization of two coupled fractional-order chaotic oscillators200526114114510.1016/j.chaos.2004.12.030ZBL1077.70013DengW. H.LiC. P.leecp@online.sh.cnChaos synchronization of the fractional Lü system20053531–4617210.1016/j.physa.2005.01.021LuJ. G.jglu@sjtu.edu.cnSynchronization of a class of fractional-order chaotic systems via a scalar transmitted signal200627251952510.1016/j.chaos.2005.04.032ZBL1086.94007YanJ. P.LiC. P.On chaos synchronization of fractional differential equations2007322725735228011510.1016/j.chaos.2005.11.062ZBL1132.37308TavazoeiM. S.HaeriM.haeri@sina.sharif.eduChaos control via a simple fractional-order controller2008372679880710.1016/j.physleta.2007.08.040MatignonD.Stability results of fractional differential equations with applications to control processing963Proceedings of the IMACS-IEEE Multiconference on Computational Engineering in Systems Applications (CESA '96)July 1996Lille, FranceOgataK.1990Englewood Cliffs, NJ, USAPrentice-HallLüJ.ChenG.ZhangS.Dynamical analysis of a new chaotic attractor200212510011015191398410.1142/S0218127402004851ZBL1044.37021LiY.TangW. K. S.ChenG.Generating hyperchaotic attractor via state feedback control200515103367337510.1142/S0218127405013988LiC. P.leecp@online.sh.cnPengG. J.Chaos in Chen's system with a fractional order200422244345010.1016/j.chaos.2004.02.013ZBL1060.37026TavazoeiM. S.HaeriM.haeri@sina.sharif.eduA necessary condition for double scroll attractor existence in fractional-order systems20073671-210211310.1016/j.physleta.2007.05.081