A new symplectic chaos synchronization of chaotic systems with uncertain chaotic parameters is studied. The traditional chaos synchronizations are special cases of the symplectic chaos synchronization. A sufficient condition is given for the asymptotical stability of the null solution of error dynamics and a parameter difference. The symplectic chaos synchronization with uncertain chaotic parameters may be applied to the design of secure communication systems. Finally, numerical results are studied for symplectic chaos synchronized from two identical Lorenz-Stenflo systems in three different cases.
1. Introduction
Chaos has been detected in a large number of nonlinear dynamic systems of physical characteristics. In addition to the control and stabilization of chaos, chaos synchronization systems are a fascinating concept which has received considerable interest among nonlinear scientists in recent times. However, chaos is desirable in some systems, such as convective heat transfer, liquid mixing, encryption, power converters, secure communications, biological systems and chemical reactions. There are a chaotic master system (or driver) and either an identical or a different slave system (or responser). Our goal is the synchronization of the chaotic master system and the chaotic slave system by coupling or by other methods. In practice, some or all of the parameters of chaotic system parameters are uncertain. A lot of works have been preceded to solve this problem using adaptive control concept.
Among many kinds of chaos synchronization, generalized chaos synchronization is investigated [1–6]. There exists a functional relationship between the states of the master system and those of the slave system. The symplectic chaos synchronization concept [7]:
(1)y=H(t,x,y)+F(t),
is studied, where x, y are the state vectors of the master system and of the slave system, respectively, F(t) is a given function of time in different form. The F(t) may be a regular function or a chaotic function. When H(t,x,y)+F(t)=x and H(t,x,y)=x of (1) reduces to the traditional generalized chaos synchronization and the traditional chaos synchronization given in [8–10], respectively.
As numerical examples, in 1996, Stenflo originally used a four-dimensional autonomous chaotic system to describe the low-frequency short-wavelength gravity wave disturbance in the atmosphere [11]. This system is similar to the celebrated Lorenz equation but more complex than it due to the introduction of a new feedback control and a new state variable and thus is called a Lorenz-Stenflo system after the names of Lorenz and Stenflo [12]. The nonlinear dynamical behaviors of the Lorenz-Stenflo system have been investigated in [13–16].
This paper is organized as follows. In Section 2, by the Lyapunov asymptotic stability theorem, the symplectic chaos synchronization with uncertain chaotic parameters by adaptive control scheme is given. In Section 3, various adaptive controllers and update laws are designed for the symplectic chaos synchronization of the identical Lorenz-Stenflo systems. Numerical simulations are also given in Section 3. Finally, some concluding remarks are given in Section 4.
2. Symplectic Chaos Synchronization with Chaotic Parameters by Adaptive Control Scheme
There are two identical nonlinear dynamical systems, and the partner A controls the partner B. The partner A is given by
(2)x˙=f(t,x,A(t)),
where x=[x1,x2,…,xn]T∈Rn is a state vector, A(t)=[A1(t),A2(t),…,Am(t)]T∈Rm is a vector of uncertain coefficients in f, and f is a vector function.
The partner B is given by(3a)y˙=f(t,y,A^(t)),
where y=[y1,y2,…,yn]T∈Rn is a state vector, A^(t)=[A^1(t),A^2(t),…,A^m(t)]T∈Rm is an estimated vector of uncertain coefficients in f.
So a controller u(t) is added on partner B, and the partner B becomes
(3b)y˙=f(t,y,A^(t))+u(t),where u(t)=[u1(t),u2(t),…,un(t)]T∈Rn is the control vector function.
Our goal is to design the controller u(t) so that the state vector y of the partner B asymptotically approaches H(t,x,y)+F(t), a given function H(t,x,y) plus a given vector function F(t)=[F1(t),F2(t),…,Fn(t)]T which is the regular function or the chaotic function.
To define error vector e(t)=[e1,e2,…,en]T:
(4)e=H(t,x,y)-y+F(t),(5)limt→∞e=0
is demanded.
From (4), it is obtained that
(6)e˙=∂H∂t+∇HTΨ˙-y˙+F˙(t),
where Ψ˙T=[x˙y˙].
By (2), (3a), and (3b), (6) can be rewritten as
(7)e˙=∂H∂t+∂H∂xf(t,x,A(t))+∂H∂yf(t,y,A^(t))-f(t,y,A^(t))-u(t)+F˙(t).
A positive definite Lyapunov function V(e,A~(t)) is chosen:
(8)V(e,A~(t))=12eTe+12A~T(t)A~(t),
where A~(t)=A(t)-A^(t).
Its derivative along any solution of (7) is
(9)V˙(e,A~(t))=eT{∂H∂t+∂H∂xf(t,x,A(t))+∂H∂yf(t,y,A^(t))-f(t,y,A^(t))+F˙(t)-u(t)∂H∂t+∂H∂x}+A~T(t)A~˙(t).
In (9), the u(t) and A~˙(t) are designed so that V˙=eTCn×ne+A~T(t)Dm×mA~(t) where Cn×n and Dm×m are two diagonal negative definite matrices. The V˙ is a negative definite function of e and A~(t). By Lyapunov theorem of asymptotical stability
(10)limt→∞e=0,limt→∞A~(t)=0.
The symplectic chaos synchronization with uncertain chaotic parameters is obtained [3–7, 11, 17–23].
3. Numerical Results for the Symplectic Chaos Synchronization of Lorenz-Stenflo System with Uncertain Chaotic Parameters via Adaptive Control
The master Lorenz-Stenflo system can be described as [11]
(11)x˙1=a(x2-x1)+rx4,x˙2=dx1-x2-x1x3,x˙3=-bx3+x1x2,x˙4=-x4-cx1.
And the slave Lorenz-Stenflo system can be described as
(12)y˙1=a^(y2-y1)+r^y4,y˙2=d^y1-y2-y1y3,y˙3=-b^y3+y1y2,y˙4=-y4-c^y1,
where a, b, c, d, and r are called the first Prandtl number, geometric parameter, second Prandtl number, Rayleigh number, and rotation number, respectively. The parameters of Lorenz-Stenflo system are chosen as a=11, b=2.9, c=5, d=23, and r=1.9.
The controllers, u1, u2, u3, and u4, are added to the four equations of (12), respectively:
(13)y˙1=a^(y2-y1)+r^y4+u1,y˙2=d^y1-y2-y1y3+u2,y˙3=-b^y3+y1y2+u3,y˙4=-y4-c^y1+u4.
The initial values of the states of the master system and of the slave system are taken as x1(0)=10, x2(0)=-13, x3(0)=12, x4(0)=-5, y1(0)=-3, y2(0)=-4, y3(0)=3, and y4(0)=-5.
Case 1.
A time delay symplectic synchronization.
We take F1(t)=x1(t-T), F2(t)=x2(t-T), F3(t)=x3(t-T), and F4(t)=x4(t-T). They are chaotic functions of time, where time delay T=0.5 sec. The Hi(x,y,t)=(xi2+yi)(sinwt+3)(i=1,2,3,4) are given. By (5) we have
(14)limt→∞ei=limt→∞((xi2+yi)(sinωt+3)-yi+xi(t-T))=0,i=1,2,3,4.
From (14) we have
(15)e˙i=(2xix˙i+y˙i)(sinωt+3)+ωcosωt(xi2+yi)-y˙i+x˙i(t-T),i=1,2,3,4,
where ω=19.
Equation (15) can be expressed as
(16)e˙1=[2x1(a(x2-x1)+rx4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+a(x2(t-T)-x1(t-T))+rx4(t-T),e˙2=[2x2(dx1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3+u2+dx1(t-T)-x2(t-T)-x1(t-T)x3(t-T),e˙3=[2x3(-bx3+x1x2)-b^y3+y1y2](sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-bx3(t-T)+x1(t-T)x2(t-T),e˙4=[2x4(-x4-cx1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T)-cx1(t-T),
where
(17)e1=(x12+y1)(sinωt+3)-y1+x1(t-T),e2=(x22+y2)(sinωt+3)-y2+x2(t-T),e3=(x32+y3)(sinωt+3)-y3+x3(t-T),e4=(x42+y4)(sinωt+3)-y4+x4(t-T).
Choose a positive definite Lyapunov function:
(18)V(e1,e2,e3,e4,a~,b~,c~,d~,r~)=12(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2),
where a~=a-a^, b~=b-b^, c~=c-c^, d~=d-d^, r~=r-r^, and a^, b^, c^, d^, and r^ are estimates of uncertain parameters a, b, c, d, and r, respectively.
Its time derivative along any solution of (16) is
(19)V˙=e1{(x2(t-T)-x1(t-T))[2x1(a(x2-x1)+rx4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+a(x2(t-T)-x1(t-T))+rx4(t-T)}+e2{[2x2(dx1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3-u2+dx1(t-T)-x2(t-T)-x1(t-T)x3(t-T)[2x2(dx1-x2-x1x3)+(d^y1-y2-y1y3)]}+e3{[2x3(-bx3+x1x2)-b^y3+y1y2]×(sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-bx3(t-T)+x1(t-T)x2(t-T)[2x3(-bx3+x1x2)-b^y3+y1y2]}+e4{[2x4(-x4-cx1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T)-cx1(t-T)[2x4(-x4-cx1)-y4-c^y1]}+a~a~˙+b~b~˙+c~c~˙+d~d~˙+r~r~˙.
The adaptive controllers are chosen as
(20)u1=[2x1(a(x2-x1)+rx4)+(a(y2-y1)+ry4)]×(sinωt+3)+ωcosωt(x12+y1)-a(y2-y1)-ry4+a(x2(t-T)-x1(t-T))+rx4(t-T)+e1,u2=[2x2(dx1-x2-x1x3)+(dy1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-dy1+y2+y1y3+dx1(t-T)-x2(t-T)-x1(t-T)x3(t-T)+e2,u3=[2x3(-bx3+x1x2)-by3+y1y2](sinωt+3)+ωcosωt(x32+y3)+by3-y1y2-bx3(t-T)+x1(t-T)x2(t-T)+e3,u4=[2x4(-x4-cx1)-y4-cy1](sinωt+3)+ωcosωt(x42+y4)+y4+cy1-x4(t-T)-cx1(t-T)+e4,
and the update laws are chosen as
(21)a~˙=-a^˙=2x1(y2-y1)(sinωt+3)e1+(y2-y1)e1-a~,b~˙=-b^˙=-2x3y3e3(sinωt+3)-y3e3-b~,c~˙=-c^˙=-2x4y1e4(sinωt+3)-y1e4-c~,d~˙=-d^˙=2x2y1e2(sinωt+3)+y1e2-d~,r~˙=-r^˙=2x1y4e1(sinωt+3)+y4e1-r~.
The initial values of estimate for uncertain parameters are a^(0)=b^(0)=c^(0)=d^(0)=r^(0)=0.
Equation (19) becomes
(22)V˙=-(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2)<0,
which is negative definite. The Lyapunov asymptotical stability theorem is satisfied. The time delay symplectic synchronization of the identical Lorenz-Stenflo systems is achieved. The numerical results of the phase portrait of master system, chaotic system (13), the time series of states, state errors, and parameter differences are shown in Figures 1, 2, and 3, respectively. The symplectic chaos synchronization is accomplished using adaptive control method.
Projections of phase portrait for master system (11).
Projections of the phase portrait for chaotic system (13) of Case 1.
Time histories of states, state errors, parameter differences, F1, F2, F3, F4, H1, H2, H3, and H4 for Case 1.
Case 2.
A time delay symplectic synchronization with uncertain chaotic parameters.
The Lorenz-Stenflo master system with uncertain chaotic parameters is
(23)x˙1=A(t)(x2-x1)+R(t)x4,x˙2=D(t)x1-x2-x1x3,x˙3=-B(t)x3+x1x2,x˙4=-x4-C(t)x1,
where A(t), B(t), C(t), D(t), and R(t) are uncertain chaotic parameters. In simulation, we take
(24)A(t)=a(1+k1z1),B(t)=b(1+k2z3),C(t)=c(1+k3z4),D(t)=d(1+k4z2),R(t)=r(1+k5z1),
where k1, k2, k3, k4, and k5 are positive constants. Chosen are k1=0.07, k2=0.08, k3=0.06, k4=0.07, and k5=0.08. The system (23) is chaotic dynamic motion, shown in Figure 4.
Projections of the phase portrait for chaotic system (23).
The F(t) is chaotic system, and the chaotic signal of goal system can be described as
(25)z˙1=a(z2-z1)+rz4,z˙2=dz1-z2-z1z3,z˙3=-bz3+z1z2,z˙4=-z4-cz1,
where initial conditions of the chaotic signal of system are z1(0) = 2, z2(0) = 5, z3(0) = −4, and z4(0) = −6. The Hi(x,y,t)=(xi2+yi)(sinwt+3)(i=1,2,3,4) are given. By (5) we have
(26)limt→∞ei=limt→∞((xi2+yi)(sinωt+3)-yi+xi(t-T))=0,i=1,2,3,4.
From (26) we have
(27)e˙i=(2xix˙i+y˙i)(sinωt+3)+ωcosωt(xi2+yi)-y˙i+x˙i(t-T),i=1,2,3,4.
Equation (27) can be expressed as
(28)e˙1=[2x1(A(t)(x2-x1)+R(t)x4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+A(t)(x2(t-T)-x1(t-T))+R(t)x4(t-T),e˙2=[2x2(D(t)x1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3+u2+D(t)x1(t-T)-x2(t-T)-x1(t-T)x3(t-T),e˙3=[2x3(-B(t)x3+x1x2)-b^y3+y1y2](sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-B(t)x3(t-T)+x1(t-T)x2(t-T),e˙4=[2x4(-x4-C(t)x1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T)-C(t)x1(t-T),
where
(29)e1=(x12+y1)(sinωt+3)-y1+x1(t-T),e2=(x22+y2)(sinωt+3)-y2+x2(t-T),e3=(x32+y3)(sinωt+3)-y3+x3(t-T),e4=(x42+y4)(sinωt+3)-y4+x4(t-T).
Choose a positive definite Lyapunov function:
(30)V(e1,e2,e3,e4,a~,b~,c~,d~,r~)=12(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2),
where a~=A(t)-a^, b~=B(t)-b^, c~=C(t)-c^, d~=D(t)-d^, r~=R(t)-r^, and a^, b^, c^, d^, and r^ are estimates of uncertain parameters A(t), B(t), C(t), D(t), and R(t), respectively.
Its time derivative along any solution of (28) is
(31)V˙=e1{[2x1(A(t)(x2-x1)+R(t)x4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+A(t)(x2(t-T)-x1(t-T))+R(t)x4(t-T)}+e2{[2x2(D(t)x1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3-u2+D(t)x1(t-T)-x2(t-T)-x1(t-T)x3(t-T)[2x2(D(t)x1-x2-x1x3)+(d^y1-y2-y1y3)]}+e3{[2x3(-B(t)x3+x1x2)-b^y3+y1y2]×(sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-B(t)x3(t-T)+x1(t-T)x2(t-T)[2x3(-B(t)x3+x1x2)-b^y3+y1y2]}+e4{[2x4(-x4-C(t)x1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T)-C(t)x1(t-T)}+a~a~˙+b~b~˙+c~c~˙+d~d~˙+r~r~˙.
The adaptive controllers are chosen as
(32)u1=[2x1(A(t)(x2-x1)+R(t)x4)+(a(y2-y1)+ry4)]×(sinωt+3)+ωcosωt(x12+y1)-a(y2-y1)-ry4+A(t)(x2(t-T)-x1(t-T))+R(t)x4(t-T)+e1,u2=[2x2(D(t)x1-x2-x1x3)+(dy1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-dy1+y2+y1y3+D(t)x1(t-T)-x2(t-T)-x1(t-T)x3(t-T)+e2,u3=[2x3(-B(t)x3+x1x2)-by3+y1y2](sinωt+3)+ωcosωt(x32+y3)+by3-y1y2-B(t)x3(t-T)+x1(t-T)x2(t-T)+e3,u4=[2x4(-x4-C(t)x1)-y4-cy1](sinωt+3)+ωcosωt(x42+y4)+y4+cy1-x4(t-T)-C(t)x1(t-T)+e4,
and the update laws are chosen as
(33)a~˙=2x1(y2-y1)(sinωt+3)e1+(y2-y1)e1-a~,b~˙=-2x3y3e3(sinωt+3)-y3e3-b~,c~˙=-2x4y1e4(sinωt+3)-y1e4-c~,d~˙=2x2y1e2(sinωt+3)+y1e2-d~,r~˙=2x1y4e1(sinωt+3)+y4e1-r~.
The initial values of estimate for uncertain parameters are a^(0)=b^(0)=c^(0)=d^(0)=r^(0)=0.
Equation (31) becomes
(34)V˙=-(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2)<0,
which is negative definite. The Lyapunov asymptotical stability theorem is satisfied. The time delay symplectic synchronization with uncertain chaotic parameters of the identical Lorenz-Stenflo systems is achieved. The numerical results of the phase portrait of chaotic system (13), the time series of states, state errors, parameters, and parameter differences are shown in Figures 5, 6, and 7, respectively. The symplectic chaos synchronization is accomplished using adaptive control method.
Projections of the phase portrait for chaotic system (13) of Case 2.
Time histories of states, state errors, F1, F2, F3, F4, H1,H2, H3, and H4 for Case 2.
Time histories of A(t), B(t), C(t), D(t), and R(t) and parameter differences for Case 2.
Case 3.
A multitime delay symplectic synchronization with uncertain chaotic parameters.
We take F1(t)=x1(t-T1), F2(t)=x2(t-T2), F3(t)=x3(t-T3), and F4(t)=x4(t-T4). They are chaotic functions of time, where multi time delay T1, T2, T3, and T4 are positive constants, T1 = 0.5 sec, T2=0.7 sec, T3=0.8 sec, and T4=0.9 sec. The Hi(x,y,t)=(xi2+yi)(sinwt+3)(i=1,2,3,4) are given. By (5) we have
(35)limt→∞ei=limt→∞((xi2+yi)(sinωt+3)-yi+xi(t-Ti))=0,i=1,2,3,4.
From (35) we have
(36)e˙i=(2xix˙i+y˙i)(sinωt+3)+ωcosωt(xi2+yi)-y˙i+x˙i(t-Ti),i=1,2,3,4.
Equation (36) can be expressed as
(37)e˙1=[2x1(A(t)(x2-x1)+R(t)x4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+A(t)(x2(t-T2)-x1(t-T1))+R(t)x4(t-T4),e˙2=[2x2(D(t)x1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3+u2+D(t)x1(t-T1)-x2(t-T2)-x1(t-T1)x3(t-T3),e˙3=[2x3(-B(t)x3+x1x2)-b^y3+y1y2]×(sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-B(t)x3(t-T3)+x1(t-T1)x2(t-T2),e˙4=[2x4(-x4-C(t)x1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T4)-C(t)x1(t-T1),
where
(38)e1=(x12+y1)(sinωt+3)-y1+x1(t-T1),e2=(x22+y2)(sinωt+3)-y2+x2(t-T2),e3=(x32+y3)(sinωt+3)-y3+x3(t-T3),e4=(x42+y4)(sinωt+3)-y4+x4(t-T4).
Choose a positive definite Lyapunov function:
(39)V(e1,e2,e3,e4,a~,b~,c~,d~,r~)=12(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2),
where a~=A(t)-a^, b~=B(t)-b^, c~=C(t)-c^, d~=D(t)-d^, r~=R(t)-r^, and a^, b^, c^, d^, and r^ are estimates of uncertain parameters A(t), B(t), C(t), D(t), and R(t), respectively.
Its time derivative along any solution of (37) is
(40)V˙=e1{[2x1(A(t)(x2-x1)+R(t)x4)+(a^(y2-y1)+r^y4)]×(sinωt+3)+ωcosωt(x12+y1)-a^(y2-y1)-r^y4-u1+A(t)(x2(t-T2)-x1(t-T1))+R(t)x4(t-T4)}+e2{[2x2(D(t)x1-x2-x1x3)+(d^y1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-d^y1+y2+y1y3-u2+D(t)x1(t-T1)-x2(t-T2)-x1(t-T1)x3(t-T3)[2x3(-B(t)x3+x1x2)-b^y3+y1y2]}+e3{[2x3(-B(t)x3+x1x2)-b^y3+y1y2](sinωt+3)+ωcosωt(x32+y3)+b^y3-y1y2-u3-B(t)x3(t-T3)+x1(t-T1)x2(t-T2)[2x3(-B(t)x3+x1x2)-b^y3+y1y2]}+e4{[2x4(-x4-C(t)x1)-y4-c^y1](sinωt+3)+ωcosωt(x42+y4)+y4+c^y1-u4-x4(t-T4)-C(t)x1(t-T1)}+a~a~˙+b~b~˙+c~c~˙+d~d~˙+r~r~˙.
The adaptive controllers are chosen as
(41)u1=[2x1(A(t)(x2-x1)+R(t)x4)+(a(y2-y1)+ry4)]×(sinωt+3)+ωcosωt(x12+y1)-a(y2-y1)-ry4+A(t)(x2(t-T2)-x1(t-T1))+R(t)x4(t-T4)+e1,u2=[2x2(D(t)x1-x2-x1x3)+(dy1-y2-y1y3)]×(sinωt+3)+ωcosωt(x22+y2)-dy1+y2+y1y3+D(t)x1(t-T1)-x2(t-T2)-x1(t-T1)x3(t-T3)+e2,u3=[2x3(-B(t)x3+x1x2)-by3+y1y2](sinωt+3)+ωcosωt(x32+y3)+by3-y1y2-B(t)x3(t-T3)+x1(t-T1)x2(t-T2)+e3,u4=[2x4(-x4-C(t)x1)-y4-cy1](sinωt+3)+ωcosωt(x42+y4)+y4+cy1-x4(t-T4)-C(t)x1(t-T1)+e4,
and the update laws are chosen as
(42)a~˙=2x1(y2-y1)(sinωt+3)e1+(y2-y1)e1-a~,b~˙=-2x3y3e3(sinωt+3)-y3e3-b~,c~˙=-2x4y1e4(sinωt+3)-y1e4-c~,d~˙=2x2y1e2(sinωt+3)+y1e2-d~,r~˙=2x1y4e1(sinωt+3)+y4e1-r~.
The initial values of estimate for uncertain parameters are a^(0)=b^(0)=c^(0)=d^(0)=r^(0)=0.
Equation (40) becomes
(43)V˙=-(e12+e22+e32+e42+a~2+b~2+c~2+d~2+r~2)<0,
which is negative definite. The Lyapunov asymptotical stability theorem is satisfied. The multi time delay symplectic synchronization with uncertain chaotic parameters of the identical Lorenz-Stenflo systems is achieved. The numerical results of the phase portrait of chaotic system (13), the time series of states, state errors, and parameter differences are shown in Figures 8 and 9, respectively. The symplectic chaos synchronization is accomplished using adaptive control method.
Projections of the phase portrait for chaotic system (13) of Case 3.
Time histories of states, state errors, parameter differences, F1, F2, F3, F4, H1, H2, H3, and H4 for Case 3.
4. Conclusions
A novel symplectic synchronization of a Lorenz-Stenflo system with uncertain chaotic parameters is obtained by the Lyapunov asymptotical stability theorem. The simulation results of three cases are shown in corresponding figures which imply that the adaptive controllers and update laws we designed are feasible and effective. The symplectic synchronization of chaotic systems with uncertain chaotic parameters via adaptive control concept can be used to increase the security of secret communication system.
Acknowledgment
This research was supported by the National Science Council, Republic of China, under Grant no. 98-2218-E-011-010.
FematR.Solís-PeralesG.On the chaos synchronization phenomena19992621506010.1016/S0375-9601(99)00667-2MR1732081ZBL0936.37010GeZ.-M.YangC.-H.Pragmatical generalized synchronization of chaotic systems with uncertain parameters by adaptive control20072312879410.1016/j.physd.2007.03.019MR2345768ZBL1167.34357YangS. S.DuanC. K.Generalized synchronization in chaotic systems19989101703170710.1016/S0960-0779(97)00149-5MR1663088ZBL0946.34040KrawieckiA.SukiennickiA.Generalizations of the concept of marginal synchronization of chaos20001191445145810.1016/S0960-0779(99)00062-4MR1745739ZBL0982.37022GeZ. M.YangC. H.ChenH. H.LeeS. C.Non-linear dynamics and chaos control of a physical pendulum with vibrating and rotating support200124222472642-s2.0-003595353010.1006/jsvi.2000.3353ChenM. Y.HanZ. Z.ShangY.General synchronization of Genesio-Tesi systems20041413473542-s2.0-864429258810.1142/S0218127404009077GeZ. M.YangC. H.Symplectic synchronization of different chaotic systems200940525322543ZBL1198.932002-s2.0-6734922803010.1016/j.chaos.2007.10.055GeZ. M.YangC. H.The generalized synchronization of a Quantum-CNN chaotic oscillator with different order systems2008355980990ZBL1141.370172-s2.0-3524886840410.1016/j.chaos.2006.05.090GeZ.-M.YangC.-H.Synchronization of complex chaotic systems in series expansion form20073451649165810.1016/j.chaos.2006.04.072MR2335412ZBL1152.37314PecoraL. M.CarrollT. L.Synchronization in chaotic systems199064882182410.1103/PhysRevLett.64.821MR1038263ZBL0938.37019StenfloL.Generalized Lorenz equations for acoustic-gravity waves in the atmosphere1996531838410.1088/0031-8949/53/1/0152-s2.0-0002804382VincentU. E.Synchronization of identical and non-identical 4-D chaotic systems using active control20083741065107510.1016/j.chaos.2006.10.005MR2411534ZBL1153.37359LiuZ.The first integrals of nonlinear acoustic gravity wave equations200061552610.1238/Physica.Regular.061a00526ZhouC.LaiC. H.YuM. Y.Bifurcation behavior of the generalized Lorenz equations at large rotation numbers199738105225523910.1063/1.531938MR1471924ZBL0897.76038BanerjeeS.SahaP.ChowdhuryA. R.Chaotic scenario in the Stenflo equations2001633177180ZBL1115.373662-s2.0-000175754810.1238/Physica.Regular.063a00177ChenY.WuX.GuiZ.Global synchronization criteria for two Lorenz-Stenflo systems via single-variable substitution control2010621-236136910.1007/s11071-010-9723-5MR2736999ZBL1217.34086GeZ.-M.ChenY.-S.Synchronization of unidirectional coupled chaotic systems via partial stability200421110111110.1016/j.chaos.2003.10.004MR2033687ZBL1048.37027YangC. H.LiS. Y.TsenP. C.Synchronization of chaotic system with uncertain variable parameters by linear coupling and pragmatical adaptive tracking20127032187220210.1007/s11071-012-0609-6GeZ. M.ChenC. C.Phase synchronization of coupled chaotic multiple time scales systems2004203639647ZBL1069.340562-s2.0-024227095810.1016/j.chaos.2003.08.001TamL. M.ChenJ. H.ChenH. K.Si TouW. M.Generation of hyperchaos from the Chen-Lee system via sinusoidal perturbation20083838268392-s2.0-4304913039210.1016/j.chaos.2007.01.039Sebastian SudheerK.SabirM.Modified function projective synchronization of hyperchaotic systems through Open-Plus-Closed-Loop coupling201037419-2020172023ZBL1236.3407210.1016/j.physleta.2010.02.068SunM.TianL.JiaQ.Adaptive control and synchronization of a four-dimensional energy resources system with unknown parameters200939419431949ZBL1197.931002-s2.0-6344908610310.1016/j.chaos.2007.06.117LiZ.ZhaoX.Generalized function projective synchronization of two different hyperchaotic systems with unknown parameters20111252607261510.1016/j.nonrwa.2011.03.009MR2813207ZBL1225.34061