# Synchronization and Antisynchronization of a Planar Oscillation of Satellite in an Elliptic Orbit via Active Control

Authors: Mohammad Shahzad
Journal: Journal of Control Science and Engineering
 (2011)

---

## Abstract

We have investigated the synchronization and antisynchronization behaviour of two identical planar oscillation of a satellite in elliptic orbit evolving from different initial conditions using the active control technique based on the Lyapunov stability theory and the Routh-Hurwitz criteria. The designed controller, with our own choice of the coefficient matrix of the error dynamics that satisfy the Lyapunov stability theory and the Routh-Hurwitz criteria, is found to be effective in the stabilization of the error states at the origin, thereby, achieving synchronization and antisynchronization between the states variables of two nonlinear dynamical systems under consideration. The results are validated by numerical simulations usingmathematica.


---

## Body

## 1. Introduction

In the last two decades, considerable research has been done in nonlinear dynamical systems and their various properties. One of the most important property of nonlinear dynamical systems is synchronization, which classically represents the entrainment of frequencies of oscillators due to weak interactions [1–3]. Studies in this field are partly motivated by experimental realization in lasers, electronic circuits, plasma discharge, and chemical reactions [2–4]. Synchronization techniques have been improved in recent years, and many different methods are applied theoretically as well as experimentally to synchronize the chaotic systems. Some of them are nonlinear feedback method [5], adaptive control method [6], antisynchronization method [7], sliding mode control method [8], and chaos synchronization using active control, which is introduced in [9] is one of these methods. Notable among these methods, chaos synchronization using active control scheme has recently been widely accepted as an efficient technique to synchronize the chaotic systems. The reason is because it can be used to synchronize identical as well as nonidentical systems; a feature that gives it an advantage over other synchronizing methods. This method based on the Lyapunov stability theory and the Routh-Hurwitz criteria to active control in order to achieve stable synchronization has been applied to many practical systems such as the electronic circuits, which model a third-order “jerk” equation [10], Lorenz, Chen, and Lü system [11], geophysical systems [12], nonlinear equations of acoustic gravity waves (Lorenz-Stenflo system) [13], Van-der Pol-Duffing oscillator [14], forced damped pendulum [15], nuclear magnetic resonance modeled by the Bloch equations [16], parametrically excited oscillators [17], permanent magnet reluctance machine [18], inertial ratchets [19], RCL-shunted Josephson junction [20], and modified projective synchronization [21].In this article, we have applied the active control technique based on the Lyapunov stability theory and the Routh-Hurwitz criteria to study the synchronization and antisynchronization behavior of two identical planar oscillation of a satellite in elliptic orbit evolving from different initial conditions. The system under consideration is chaotic for some values of parameter involved in the system. In synchronization, the two systems (master and slave) are synchronized and start with different initial conditions. The same problem may be treated as the design of control laws for full chaotic slave system using the known information of the master system so as to ensure that the controlled receiver synchronizes with the master system. Hence, the slave chaotic system completely traces the dynamics of the master system in the course of time. The aim of this study is to trace the chaotic dynamics of the planar oscillation of a satellite in elliptic orbit based on synchronization and antisynchronization. To the best of my knowledge nobody studied this before.


## 2. Description of the Model

Elliptically orbiting planar oscillations of satellites in the solar system make an interesting study, and significant contributions to this end can be found [22–30], all of which have studied the influence of certain perturbative forces, such as solar radiation pressure, tidal force, and air resistance. In the present work, we consider the planar oscillation of a satellite in elliptic orbit with the spin axis fixed perpendicular to the orbital plane. Let the long axis of the satellite makes an angle θ with a reference axis that is fixed in inertial space, the long axis of the satellite makes an angle ψ with satellites planet centre line [31], and the satellite to be a triaxial ellipsoid with principal moments of inertia A<B<C, where C is the moment about the spin axis. The orbit is taken to be a fixed ellipse with semimajor axis a, eccentricity e, true anomaly f, and instantaneous radius r. The equation of motion of satellite planar oscillation in an elliptic orbit around the earth is
(1)2Cθ̈-3(B-A)Gmpr3sin⁡2ψ=0.
Taking 3(B-A)/C=ω02, Gmp=μ, ψ=f-θ and l/r=1+ecosf, (1) can be written as
(2)d2θdt2=μω02l3(1+ecos⁡⁡f)3sin⁡2(f-θ).
In order to reduce three variables ψ,θ, and f we are taking fas independent variable.Now,(3)dθdt=dθdfdfdt=ḟdθdf=hr2dθdf,(usingr2ḟ=h),d2θdt2=-2hṙr3dθdf+hr2d2θdf2dfdt,d2θdt2=-2hṙr3dθdf+h2r4d2θdf2(usingr2ḟ=h).
Now, from  l/r=1+ecosf, we have
(4)ṙ=-l(-esin⁡f)ḟ(1+ecos⁡⁡f)2,ṙ=elsin⁡f(1+ecos⁡⁡f)2hr2.Using (3) and (4) in (2), we have
(5)h2r4d2θdf2-2elh2sin⁡fr5(1+ecos⁡⁡f)2dθdf=μω02r3sin⁡2(f-θ)h2r4d2θdf2-2eh2sin⁡fr4(1+ecos⁡⁡f)dθdf=h2ω02lr3sin⁡2(f-θ)(usingh2=μl),d2θdf2=2esin⁡f(1+ecos⁡⁡f)dθdf+ω02sin⁡2(f-θ)(1+ecos⁡⁡f).Using the binomial expansion for(1+ecosf)-1 and for very small value of e ignoring the higher order terms in e, we have
(6)d2θdf2=2esin⁡fdθdf+ω022(1-ecos⁡⁡f)sin⁡2(f-θ).


## 3. Synchronization via Active Control

For a system of two coupled chaotic oscillators, the master system (ẋ=f(x,y)) and the slave system (ẏ=g(x,y)), where x(t) and y(t) are the phase space (state variables), and f(x,y) and g(x,y) are the corresponding nonlinear functions, synchronization in a direct sense implies limt→∞|x(t)-y(t)|→0. When this occurs the coupled systems are said to be completely synchronized. Chaos synchronization is related to the observer problem in control theory [32]. The problem may be treated as the design of control laws for full chaotic slave system using the known information of the master system so as to ensure that the controlled receiver synchronizes with the master system. Hence, the slave chaotic system completely traces the dynamics of the master in the course of time.In order to formulate the active controllers, we write the system (6) in two first-order differential equations as shown below.Letdθ/df=x1′ and d2θ/df2=x2′, then we have
(7)x1′=x2,x2′=2ex1′sin⁡f+ω022(1-ecos⁡⁡f)sin⁡2(f-x1).
Let us define another system
(8)y1′=y2+u1(f),y2′=2ey1′sin⁡f+ω022(1-ecos⁡⁡f)sin⁡2(f-y1)+u2(f),
where (7) and (8) are called the master system and the slave system, respectively, and in slave system, u1(f) and u2(f) are control functions to be determined. Let e1(f)=y1(f)-x1(f) and e2(f)=y2(f)-x2(f) be the synchronization errors such that limf→∞ei(f)→0 for i=1,2. From (7) and (8), we have
(9)e1′(f)=e2(f)+u1(f),e2′(f)=2ee1′(f)sin⁡f+ω022(1-ecos⁡⁡f)×[sin⁡2(f-y1)-sin⁡2(f-x1)]+u2(f).
In order to express (9) as only linear terms in e1(f) and e2(f), we redefine the control functions as follows:
(10)u1(f)=v1(f),u2(f)=-2ee1′(f)sin⁡f-ω022(1-ecos⁡⁡f)×[sin⁡2(f-y1)-sin⁡2(f-x1)]+v2(f).
From (9) and (10), we have
(11)e1′(f)=e2(f)+v1(f),e2′(f)=v2(f).
Equation (11) is the error dynamics, which can be interpreted as a control problem where the system to be controlled is a linear system with control inputs v1(f)=v1(e1(f),e2(f)) and v2(f)=v2(e1(f),e2(f)). As long as these feedbacks stabilize the system, limf→∞ei(f)→0 for i=1,2. This simply implies that the two systems (7) and (8) evolving from different initial conditions are synchronized. As functions of e1(f) and e2(f), we choose v1(f) and v2(f) as follows:
(12)(v1(f)v2(f))=D(e1(t)e2(t)),
where D=(abcd) is a 2×2 constant feedback matrix to be determined. Hence the error system (11) can be written as
(13)(e1′(f)e2′(f))=C(e1(f)e2(f)),
where C=(a1+bcd) is the coefficient matrix.According to the Lyapunov stability theory and the Routh-Hurwitz criteria, if(14)a+d<0,c(1+b)-ad<0,
then the eigen values of the coefficient matrix of error system (11) must be real or complex with negative real parts and, hence, stable synchronized dynamics between systems (7) and (8) is guaranteed. Let
(15)a+d=c(1+b)-ad=-E,
where E>0 is a real number which is usually set equal to 1. There are several ways of choosing the constant elements a, b, c, d of matrix D in order to satisfy the Lyapunov stability theory and the Routh-Hurwitz criteria (14).


## 4. Numerical Simulation for Synchronization

For the constant elements of feedback matrix, choosinga=d=-0.5 and for the parameters involved in system under investigation, e=0.15, ω0=0.3 with the initial conditions  [x1(0),y1(0)]=[0,0], and [x2(0),y2(0)]=[1,1.5], we have simulated the system under consideration using mathematica. The results obtained show that the system under consideration achieved synchronization. Phase portraits, time series analysis, and error diagrams are the witness of achieving synchronization between master and slave system. Further, it also has been confirmed by the convergence of the synchronization quality defined by
(16)e(f)=e12(f)+e22(f).
Figures 1–6 confirms the convergence of the synchronization quality defined by (16).Figure 1
Phase plot of master system.Figure 2
Phase plot of slave system.Figure 3
Time series analysis of master system.Figure 4
Time series analysis of slave system.Figure 5
Time series analysis ofe1(f).Figure 6
Time series analysis ofe2(f).


## 5. Antisynchronization via Active Control

Antisynchronization of two coupled systemsẋ=f(x,y) (master system) and ẏ=g(x,y) (slave system) means limt→∞|x(t)+y(t)|→0. This phenomenon has been investigated both experimentally and theoretically in many physical systems [18, 19, 32–37]. A recent study of the antisynchronization phenomenon in nonequilibrium systems suggests that it could be used as a technique for particle separation in a mixture of interacting particles [19].In order to formulate the active controllers for Antisynchronization, we need to redefine the error functions ase1(f)=y1(f)+x1(f) and e2(f)=y2(f)+x2(f), where e1(f) and e2(f) are called the antisynchronization errors such that limf→∞ei(f)(i=1,2)→0. From (7) and (8), error dynamics can be written as
(17)e1′(f)=e2(f)+u1(f),e2′(f)=2ee1′(f)sin⁡f+ω022(1-ecos⁡⁡f)×[sin⁡2(f-y1)+sin⁡2(f-x1)]+u2(f).
In order to express (17) as only linear terms in e1(f) and e2(f), we redefine the control functions as follows:
(18)u1(f)=v1(f),u2(f)=-2ee1′(f)sin⁡f-ω022(1-ecos⁡⁡f)×[sin2(f-y1)+sin⁡2(f-x1)]+v2(f).
From (17) and (18), we have
(19)e1′(f)=e2(f)+v1(f),e2′(f)=v2(f).
Equation (19) is the error dynamics, which can be interpreted as a control problem where the system to be controlled is a linear system with control inputs v1(f)=v1(e1(f),e2(f)), and v2(f)=v2(e1(f),e2(f)). As long as these feedbacks stabilize the system, limf→∞ei(f)→0 for i=1,2. This simply implies that the two systems (7) and (8) evolving from different initial conditions are synchronized. As functions of e1(f) and e2(f), we choose v1(f) and v2(f) as follows:
(20)(v1(f)v2(f))=D(e1(f)e2(f)),
where D=(abcd) is a 2×2 constant feedback matrix to be determined. Hence the error system (19) can be written as
(21)(e1′(f)e2′(f))=C(e1(f)e2(f)),
where C=(a1+bcd), is the coefficient matrix.According to the Lyapunov stability theory and the Routh-Hurwitz criteria, if(22)a+d<0,c(1+b)-ad<0,
then the eigen values of the coefficient matrix of error system (19) must be real or complex with negative real parts, and, hence, stable synchronized dynamics between systems (7) and (8) is guaranteed. Let
(23)a+d=c(1+b)-ad=-E,
where E>0 is a real number which is usually set equal to 1. There are several ways of choosing the constant elements a, b, c, d of matrix D in order to satisfy the Lyapunov stability theory and the Routh-Hurwitz criteria (19).


## 6. Numerical Simulation for Antisynchronization

For the constant elements of feedback matrix, choosinga=d=-0.5 and for the parameters involved in system under investigation, e=0.15, ω0=0.3 with the initial conditions [x1(0),y1(0)]=[0,0] and [x2(0),y2(0)]=[1,1.5], we have simulated the system under consideration using mathematica. The results obtained show that the system under consideration achieved synchronization. Phase portraits, time series analysis, and error diagrams are the witness of achieving synchronization between master and slave system. Further, it also has been confirmed by the convergence of the synchronization quality defined by
(24)e(f)=e12(f)+e22(f).
Figures 7–12 confirms the convergence of the synchronization quality defined by (24).Figure 7
Phase plot of slave system.Figure 8
Time series analysis of slave system.Figure 9
Time series analysis ofe1(f).Figure 10
Time series analysis ofe2(f).Figure 11
Convergence of errors in synchronization.Figure 12
Convergence of errors in antisynchronization.


## 7. Conclusion

In this paper, we have investigated the synchronization and antisynchronization behaviour of the two identical planar oscillation of a satellite in elliptic orbit evolving from different initial conditions via the active control technique based on the Lyapunov stability theory and the Routh-Hurwitz criteria. The results were validated by numerical simulations usingmathematica. For the errors in synchronization and antisynchronization behavior of the system under study, we have observed that the rate of convergence of errors is faster in antisynchronization. Figures 1–12 are proof of it.



---

*Source: 816432-2011-08-22.xml*