# Robust Adaptive Exponential Synchronization of Two Different Stochastic Perturbed Chaotic Systems with Structural Perturbations

Authors: A. Soup Tewa Kammogne; H. B. Fotsin
Journal: International Journal of Engineering Mathematics
 (2015)

---

## Abstract

The robust adaptive exponential synchronization problem of stochastic chaotic systems with structural perturbations is investigated in mean square. The stochastic disturbances are assumed to be Brownian motions that act on the slave system and the norm-bounded uncertainties exist in all parameters after decoupling. The stochastic disturbances could reflect more realistic dynamical behaviors of the coupled chaotic system presented within a noisy environment. By using a combination of the Lyapunov functional method, the robust analysis tool, the stochastic analysis techniques, and adaptive control laws, we derive several sufficient conditions that ensure the coupled chaotic systems to be robustly exponentially synchronized in the mean square for all admissible parameter uncertainties. This approach cannot only make the outputs of both master and slave systems reachH
∞ synchronization with the passage of time between both systems but also attenuate the effects of the perturbation on the overall error system to a prescribed level. The main results are shown to be general enough to cover many existing ones reported in the literature.


---

## Body

## 1. Introduction

Synchronization is a fundamental phenomenon that enables coherent behavior in coupled dynamical systems. Since Pecora and Carroll [1] proposed a method to synchronize two identical chaotic systems with differential initial conditions, chaos synchronization has attracted much attention from various fields during the last decades [2, 3]. Many methods and experimental techniques have been presented to realize the synchronization of two identical or different chaotic systems [4], such as adaptive control [5], sliding mode control [6, 7], nonlinear feedback control [8, 9], and fuzzy system based control [10]. Among all these methods, adaptive control method has been used widely to treat the unknown parameters and uncertainties [11, 12]. For example, a strategy for adaptive synchronization of some electrical chaotic circuit based nonlinear control and robust adaptive generalized projective synchronization of Genesio-Tesi chaotic system with uncertain parameters are discussed in [13, 14], respectively. In [15], the stability of adaptive synchronization of chaotic systems was proposed by Sorrentino et al. A robust adaptive synchronization for a class of uncertain chaotic systems with application to Chua’s circuit was presented by Koofigar et al. [16]. Hu et al. developed a robust adaptive finite-time synchronization of two different chaotic systems with parameter uncertainties [17]. Bowong and Tewa investigated the unknown inputs’ adaptive observer for a class of chaotic systems with uncertainties [18]. However, we have noted that, in all of the above mentioned papers, the chaotic systems are deterministic differential equations without any random parameters. In addition, chaos synchronization is far from being straightforward because different aspects affect significantly the ability of systems to synchronize. A particularly important aspect is to consider the presence of many disturbances and uncertainties sources, which are unavoidable parts of any practical implementation. Therefore, robustness is a very desirable characteristic of a synchronization approach. In general, it is worth pointing out that the stochastic factors like uncertainty of the structure parameter, perturbation of external noise, and stochastic input are ubiquitous in nature, society, economy, and realistic engineering. Under normal circumstances, those stochastic factors just play a minor influence. However, when the designer of the system needs to make a choice, it will become a dominant factor which could affect the trend of deterministic system [19].Synchronization of chaotic systems affected by both structural and stochastic disturbances poses new challenges for the understanding of stability, sensitivity and robustness, bifurcations and chaos, and so forth. When analyzing the dynamical behaviors of chaotic systems, stochastic disturbances and modeling errors are probably two of the main sources that result in uncertainties. To overcome these difficulties, various adaptive synchronization schemes have been proposed and investigated (see an excellent text in [20]). An adaptive almost surely asymptotically synchronization for stochastic delayed neural networks with Markovian switching was discussed by Ding et al. [21] while Li and Fu proposed the synchronization of chaotic delayed neural networks with impulsive and stochastic perturbations [22]. M. Chen and W.-H. Chen proposed robust adaptive neural network synchronization controller design for a class of time delay uncertain chaotic systems in [23] and a robust stability and H
∞-control of uncertain systems with impulsive perturbations [24]. Recently, Fang et al. developed a robust adaptive exponential synchronization of stochastic perturbed chaotic delayed neural networks with parametric uncertainties [25]. It appears clearly from the above works that most did not consider structural perturbations. Some previous research such as [24, 25] took into account the structured uncertainties with the stochastic perturbations but supposed that the parametric uncertainties are time varying and norm-bounded and satisfy certain elementary conditions such that the uncertain matrices can be decomposed in product of Lebesque measurable function in order to develop the problem in terms of linear matrix inequality. However, this assumption seems bulky in engineering applications and also leads to the unpredictability of the LMI. According to the best of our knowledge, there are still few results about the synchronization of chaotic systems with both structural perturbations and stochastic disturbances. Therefore, robust adaptive synchronization analysis for uncertain stochastic systems has emerged as a challenging research issue.In this paper, we discuss the asymptotical synchronization and almost surely synchronization for two different chaotic systems with the consideration of both structural and stochastic perturbations. With the help from the Lyapunov functional method and adaption method, we employ the robust analysis tool and the stochastic analysis techniques to derive some relevant conditions under which the coupled systems is globally robustly synchronized in the mean square for all admissible parameter uncertainties. These conditions guarantee the robustness of the controller against the effect of exoteric perturbations. The rest of the paper is organized as follows: Section2 presents the Problem’s formulation and preliminaries. In Section 3 we present our main result that consists of a new controller that is robust enough against the effect of exoteric perturbation of modified Colpitts oscillator and Chua’s circuit. The computer simulation is given in Section 4. Finally, conclusions are presented in Section 5.


## 2. Problem’s Formulation and Preliminaries

Let two classes of nonlinear systems be given in the following form:(1)
d
x
τ
=
A
x
τ
+
D
φ
x
τ
+
B
d
τ
,
y
τ
=
C
x
τ
,
(2)
d
x
^
τ
=
F
x
^
τ
+
G
f
x
^
τ
d
τ
,
y
~
τ
=
C
x
^
τ
,where x
∈
R
n and x
^
∈
R
n are the systems states and y
∈
R
n
×
1 and y
~
∈
R
n
×
1 are the output vectors of systems (1) and (2), respectively. A
,
F
,
D
,
G
,
C
∈
R
n
×
n; B
∈
R
n
×
1 are known matrices such that (
A
,
C
) and (
F
,
C
) are observable. (
φ
,
f
)
:
R
n
→
R
n are nonlinear continuously differentiable vector function.Considering the effect of structural and exoteric perturbation on system’s parameters, the master system is given by(3)
d
x
τ
=
A
+
Δ
A
τ
x
τ
+
D
φ
x
τ
+
B
+
Δ
B
τ
d
τ
,
y
τ
=
C
x
τ
.Based on the observer method, the slave is constructed with nonidentical parametric uncertainties and stochastic perturbation given by(4)
d
x
^
τ
=
F
+
Δ
F
τ
x
^
τ
+
G
f
x
^
τ
+
u
τ
d
τ
+
σ
t
,
x
τ
-
x
^
τ
d
ω
~
τ
,where Δ
A
(
τ
) and Δ
F
(
τ
) are bounded structural variations of the system which satisfy the condition Δ
(
•
)
(
τ
)
≤
δ (δ being a positive constant) and u
(
τ
) is the nonlinear controller. y is measurable output. ω
~
(
τ
)
=
(
ω
~
1
(
τ
)
,
…
,
ω
~
n
(
τ
)
)
T is n-dimensional Brownian motion defined on complete probability space ℧
,
F
,
P with a natural filtration F
τ
τ
≥
0 generated by ω
~
(
s
)
:
0
≤
s
≤
τ, where ℧ is associated with the canonical space generated by ω
~
(
τ
) with the probability measure P. σ
:
R
+
×
R
n
×
R
n
→
R
n
×
n is the noise intensity matrix.Define the synchronization error vector as(5)
e
τ
=
x
τ
-
x
^
τ
,
y
e
τ
=
y
-
y
~
τ
.So the error dynamical system is(6)
d
e
τ
=
A
+
F
e
τ
-
Δ
F
τ
x
^
τ
+
Δ
A
τ
x
τ
+
g
x
τ
,
x
^
τ
+
Δ
B
τ
-
u
τ
d
τ
+
σ
t
,
e
τ
d
ω
~
τ
,where (7)
g
x
τ
,
x
^
τ
=
D
φ
x
τ
-
G
f
x
^
τ
-
F
x
τ
+
A
x
^
τ
+
B
. It is clear that the synchronization problem can be transformed to be the equivalent problem of stabilizing the error system (6).Definition 1.
Systems (3) and (4) are said to be globally stable in mean squares if for any given initial condition the following equality is verified: (8)
lim
τ
→
∞
⁡
E
e
τ
2
=
0
,where E
· is the mathematical expectation.Assumption 2.
Let us defineh
(
τ
)
=
Δ
A
(
τ
)
-
Δ
F
(
τ
) and assume h
(
τ
)
∈
L
2; that is, ∫
0
∞
h
T
(
τ
)
h
(
τ
)
d
τ
≤
∞.
Given a scalarγ
>
0, the error system is said to be robustly stochastically stable with noise attenuation γ if it is robustly stochastically stable and under zero initial conditions: (9)
y
e
τ
E
2
<
γ
n
τ
,
h
2
.For nonzero n
(
τ
)
∈
L
2
0
,
∞, (10)
y
e
τ
E
2
=
E
∫
0
∞
y
e
τ
2
d
τ
1
/
2
.In this paper, the following performance index for a prescribed γ
>
0 is considered: (11)
J
e
τ
=
E
∫
0
∞
y
e
T
τ
y
e
τ
-
γ
2
n
T
τ
n
τ
d
τ
.The problem to be tackled is formulated as follows. For the stochastic systems (3) and (4) and a given scalar γ
>
0, we like to design an adaptive feedback controller K such that the resulting master-slave trajectories robustly stochastically synchronize with noise attenuation level γ which is the disturbance attenuation level [26] defined by (12)
γ
=
E
∫
0
∞
y
e
T
τ
y
e
τ
d
τ
∫
0
∞
n
T
τ
n
τ
d
τ
.The control requirement is that the response of error to the disturbances with the stabilization controller should be brought to less than γ which is the suboptimal H
∞ limits. For the standard H
∞ problem, a controller exists if and only if a unique positive definite solution to the two Riccati equations exists [27]. State space formulas can be derived for all controllers such that the H
∞-norm of the closed loop transfer function is less than γ.Assumption 3.
Assume that the noise intensity function matrixσ
:
R
+
×
R
n
→
R
n
×
n is uniformly Lipschitz continuous in terms of the norm induced by the trace inner product on the matrices as follows:(13)
t
r
a
c
e
σ
T
t
,
v
σ
t
,
v
≤
v
T
J
T
J
vfor all v
∈
R
n. J is a real positive n
×
n matrix.Definition 4.
There exists continuous functionV
∈
C
1,2
R
+
×
R
n
;
R
+ such that an operator I
V from R
+
×
R
n along the trajectory of the error system (6) is defined as (14)
I
V
τ
,
e
=
V
τ
τ
,
e
+
V
e
τ
,
e
f
τ
,
e
+
1
2
t
r
a
c
e
σ
T
τ
,
e
V
e
e
σ
τ
,
e
,where(15)
V
τ
τ
,
e
=
∂
V
τ
,
e
τ
∂
τ
,
V
e
τ
,
e
=
∂
V
τ
,
e
τ
∂
e
1
,
…
,
∂
V
τ
,
e
τ
∂
e
n
,
V
e
e
τ
,
e
=
∂
2
V
τ
,
e
τ
∂
e
i
∂
e
j
n
×
n
.Definition 5.
Systems (3) and (4) are said to be exponentially synchronized in mean square if there exist constants ℏ
>
0 and ƛ
>
0, such that(16)
E
e
τ
2
≤
ƛ
exp
⁡
-
ℏ
τwith τ
>
0 and ℏ is called the decay rate of exponential synchronization.Lemma 6.
Suppose two continuous functionsf
Q
∈
(
R
+
×
R
→
R
) and π
(
τ
,
x
(
τ
)
) a nonzero decreasing function such that f
Q
(
τ
,
x
(
τ
)
)
=
-
x
T
(
τ
)
Q
x
(
τ
)
+
x
T
π
(
τ
,
x
(
τ
)
) and then there exists a suitable positive constant η
′ satisfying the following inequality:(17)
f
Q
+
η
′
I
τ
,
x
τ
≤
-
η
′
x
T
τ
x
τ
.


## 3. Main Results

Assumption 7.
We can choose the controller in this following form:(18)
u
τ
=
1
2
K
L
y
-
C
x
^
τ
,where L is the matrix with appropriate dimension and K is the estimated feedback gain which is updated according to the following adaptation algorithm: (19)
K
˙
=
ρ
y
-
C
x
^
τ
2
,where ρ is a positive constant to be specified by the designer.
The control problem with a feedback (18) is to find an admissible internally stabilizing control K which would be attenuating disturbances such that the norm of the stable closed loop system from the disturbances to the controlled outputs is less than γ (γ is equal to 1 for optimal and slightly greater than 1 for suboptimal control). The control objective is stated in the mathematical form in (9). One advantage of this type of controller is that it can be easily constructed through time varying resistors, capacitors, or operational amplifiers and their combinations or using a digital signal processor together with the appropriate converters.Assumption 8.
There exist two known positive functionsρ
x
(
τ
)
,
τ and ρ
x
^
(
τ
)
,
τ such that g
x
(
τ
)
,
x
^
(
τ
)
=
M
ρ
(
x
,
τ
)
-
ρ
(
x
^
,
τ
). We assume the following conditions: (1)
For boundedx and x
^, ρ
(
x
,
τ
)
-
ρ
(
x
^
,
τ
) is also bounded.
(2)
The functiong
(
x
,
x
^
) enters the closed ball Λ and remains there with (20)
Λ
=
x
,
x
^
∈
R
n
×
n
and
k
0
∈
R
/
g
x
,
x
^
=
M
ρ
x
,
τ
-
ρ
x
^
,
τ
≤
M
k
0
x
-
x
^
.
M and k
0 are positive constants.Assumption 9.
The “robust performant” gain scheduled system is assumed to be bounded in response to certain added inputs. Moreover, the performance level must not degrade discontinuously for arbitrarily small errors in controller (18).Lemma 10 (Gronwall inequality (see [28])).
Let us consider three continuous positive functionsφ, ϕ, and y on [
0
,
T
′
] satisfying the inequality ∀
τ
∈
0
,
T
′, y
τ
≤
φ
τ
+
∫
0
τ
ϕ
τ
y
s
d
s, and then ∀
t
∈
0
,
T
′, y
(
τ
)
≤
φ
(
τ
)
+
∫
0
τ
φ
s
ϕ
s
exp
⁡
∫
s
τ
ϕ
(
u
)
d
u
d
s.Lemma 11.
For two integrable functionsV
∈
C
1,2
R
+
×
R
n
;
R
+ and e
(
τ
)
∈
R
3, there exist two positive constants μ
1 and μ
2 such that (21)
μ
1
E
e
τ
2
≤
E
V
τ
≤
μ
2
E
e
τ
2
.Lemma 12.
For any matricesQ
1
∈
R
n
×
m, Q
2
∈
R
n
×
m, Θ
=
Θ
T
>
0, and Θ
∈
R
n
×
n we have the following inequalities: (22)
Q
1
T
Q
2
+
Q
2
T
Q
1
≤
Q
1
T
Θ
Q
1
+
Q
2
T
Θ
-
1
Q
2
.Assumption 13.
We can choose a constant matrixL
∈
R
n
×
n and two positive definite matrices P
=
P
T
>
0 and Q
=
Q
T such that(23)
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
+
2
P
L
C
+
Ω
2
P
2
+
σ
I
=
-
Q
,
(24)
L
T
P
=
C
,where (25)
Ω
2
=
η
+
γ
-
2
d
2
+
δ
1
+
α
1
M
,
σ
=
η
-
1
δ
2
+
δ
1
-
1
δ
2
+
α
1
-
1
k
0
2
M
.Theorem 14.
System (3) exponentially synchronizes system (4) under controller (18) in mean square.Proof.
The closed loop error system with control (18) can be written in the following form: (26)
d
e
τ
=
A
+
F
e
τ
-
Δ
F
τ
x
^
τ
+
Δ
A
τ
x
τ
+
g
x
τ
,
x
^
τ
+
Δ
B
τ
-
1
2
K
L
y
~
-
C
x
^
τ
d
τ
+
σ
t
,
e
τ
d
ω
~
τ
.Consider the Lyapunov function: (27)
V
τ
=
e
T
τ
P
e
τ
+
1
2
ρ
K
τ
-
K
0
2
,where K
0 is a constant to be defined later.
Relation (25) yields by differentiating V with respect to time τ
(28)
I
V
τ
=
e
˙
T
τ
P
e
τ
+
e
T
τ
P
e
˙
τ
+
1
ρ
K
τ
-
K
0
K
˙
=
A
+
F
-
L
C
e
τ
+
Δ
A
τ
e
τ
-
Δ
F
x
^
τ
+
Δ
B
τ
T
P
e
τ
+
e
T
τ
P
A
+
F
-
L
C
e
τ
+
Δ
A
τ
x
τ
-
Δ
F
x
^
τ
+
Δ
B
τ
T
+
2
x
-
T
P
L
C
e
-
K
e
T
P
L
C
x
-
+
2
e
T
P
g
x
,
x
^
+
1
ρ
K
τ
-
K
0
K
˙
+
t
r
a
c
e
σ
T
τ
,
e
τ
P
σ
τ
,
e
τ
=
e
T
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
e
τ
+
2
e
T
τ
P
Δ
A
τ
x
τ
-
Δ
F
τ
x
^
τ
+
2
e
T
τ
P
Δ
B
τ
+
2
e
T
P
L
C
e
-
K
e
T
P
L
C
e
+
2
e
T
P
g
x
,
x
^
+
1
ρ
K
τ
-
K
0
K
˙
+
t
r
a
c
e
σ
T
τ
,
e
τ
P
σ
τ
,
e
τ
=
e
T
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
e
τ
+
2
e
T
τ
P
Δ
F
e
τ
+
2
e
T
τ
P
Δ
A
τ
-
Δ
F
τ
e
τ
+
2
e
T
τ
P
Δ
B
τ
+
2
e
T
P
L
C
e
-
K
e
T
P
L
C
e
+
2
e
T
P
g
x
,
x
^
+
1
ρ
K
-
K
0
K
˙
+
t
r
a
c
e
σ
T
τ
,
e
τ
P
σ
τ
,
e
τ
.It is clear from Lemma 12 that (29)
2
e
T
P
Δ
B
τ
≤
δ
1
e
τ
T
P
2
e
τ
+
δ
1
-
1
δ
2
e
τ
T
e
τ
,
2
e
T
τ
P
Δ
F
x
-
τ
≤
η
e
T
τ
P
2
e
τ
+
η
-
1
δ
2
e
T
τ
e
τ
,
2
e
T
τ
P
Δ
A
-
Δ
F
e
τ
≤
γ
-
2
d
2
e
T
τ
P
2
e
τ
+
2
γ
2
e
T
τ
P
h
T
h
,
2
e
T
P
g
x
,
x
-
≤
α
1
M
e
T
τ
P
2
e
τ
+
α
1
-
1
k
0
2
M
e
T
τ
e
τ
.For any time τ
≥
0, the system state is bounded (as the system is chaotic); hence we may write x
(
τ
)
≤
d. Here d, η, δ
1, α
1, and γ are positive constants.
Relation (28) leads to(30)
I
V
τ
≤
e
T
τ
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
e
τ
+
η
e
T
τ
P
2
e
τ
+
η
-
1
δ
2
e
T
τ
e
τ
+
γ
-
2
d
2
e
T
τ
P
2
e
τ
+
2
γ
2
e
T
τ
P
h
T
τ
h
τ
+
δ
1
-
1
δ
2
e
T
τ
e
τ
+
δ
1
e
T
τ
P
2
e
τ
+
2
x
-
T
P
L
C
e
-
J
T
J
+
K
e
T
P
D
C
e
+
α
1
M
e
T
τ
P
2
e
τ
+
α
1
-
1
k
0
2
M
e
T
τ
e
τ
+
1
ρ
K
-
K
0
K
˙
≤
e
T
τ
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
+
η
P
2
+
η
-
1
δ
2
I
+
γ
-
2
d
2
P
2
+
δ
1
-
1
δ
2
I
+
δ
1
P
2
e
τ
+
2
γ
2
e
T
τ
P
h
T
τ
h
τ
+
2
e
T
P
L
C
e
-
K
e
T
P
L
C
e
+
α
1
M
e
T
τ
P
2
e
τ
+
α
1
-
1
k
0
2
M
e
T
τ
e
τ
+
e
T
τ
J
T
J
e
τ
+
1
ρ
K
-
K
0
K
˙
≤
e
T
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
+
2
P
L
C
+
η
+
γ
-
2
d
2
+
δ
1
+
α
1
M
P
2
+
η
-
1
δ
2
+
δ
1
-
1
δ
2
+
α
1
-
1
k
0
2
M
I
+
J
T
J
e
+
2
γ
2
e
T
τ
P
h
T
τ
h
τ
-
K
L
T
P
e
2
+
1
ρ
K
-
K
0
K
˙
.Let one recall relations (23): (31)
Ω
2
=
η
+
γ
-
2
d
2
+
δ
1
+
α
1
M
,
σ
=
η
-
1
δ
2
+
δ
1
-
1
δ
2
+
α
1
-
1
k
0
2
M
and  let
2
γ
2
=
K
0
.
Remark
1. There exist a lot of papers that study H
-
∞ synchronization problems of master-slave chaotic systems, in which the choice of the performance index which satisfies condition (9) remains a hard task in practice and leads sometimes to the poor performance of the control. Our method does not need the arbitrary choice of this quantity. Besides, the control objective is obtained by the selection of the initial value of the feedback gain which enables approaching the synchronization. This partly shows the less conservativeness of our control strategy.
Using inequality (31), we obtain (32)
I
V
τ
≤
e
T
τ
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
+
2
P
L
C
+
Ω
2
P
2
+
J
T
J
+
σ
I
e
τ
+
2
γ
2
e
T
τ
P
h
T
τ
h
τ
-
K
˙
ρ
.Finally, by letting (33)
ϖ
T
τ
ϖ
τ
=
γ
2
h
T
τ
h
τ
-
K
˙
ρ
,
A
+
F
-
L
C
T
P
+
P
A
+
F
-
L
C
+
2
P
L
C
+
Ω
2
P
2
+
J
T
J
+
σ
I
=
-
Qone obtains(34)
I
V
τ
≤
-
e
T
Q
e
+
2
e
T
τ
P
ϖ
T
τ
ϖ
τ
.Then, if the inequality -
e
T
Q
e
+
2
e
T
(
τ
)
P
ϖ
T
(
τ
)
ϖ
(
τ
)
<
0 holds, we have(35)
I
V
τ
≤
e
τ
ϖ
τ
T
-
Q
P
P
0
e
τ
ϖ
τ
≤
0
.For all nonzero ϖ
(
τ
)
∈
L
2
0
,
∞, under zero initial conditions, we have E
V
(
0
)
=
V
0
a
n
d
E
V
(
τ
)
≥
0. Integrating both sides of (35) from 0 to τ
>
0 and taking the expectation, we have(36)
E
V
τ
=
E
∫
0
τ
I
V
t
d
t
.We derive the cost function as(37)
J
e
=
E
∫
0
τ
y
e
T
t
y
e
t
-
γ
2
ϖ
T
t
ϖ
t
+
I
V
t
d
t
-
E
V
t
≤
E
∫
0
τ
y
e
T
t
y
e
t
-
γ
2
ϖ
T
t
ϖ
t
+
I
V
t
d
t
.Substituting (35) into (37) yields(38)
J
e
e
τ
,
ϖ
τ
≤
E
e
τ
ϖ
τ
T
Ψ
e
τ
ϖ
τ
,where(39)
Ψ
=
-
Q
+
C
-
P
P
-
γ
2
I
,
C
-
=
C
T
C
0
0
0
.For a positive constant ƛ, if the following condition holds:(40)
Ψ
^
=
Ψ
+
ƛ
P
0
0
0
=
-
Q
+
C
+
ƛ
P
P
P
-
γ
2
I
<
0
,then we have (41)
J
e
e
τ
,
ϖ
τ
≤
E
e
τ
ϖ
τ
T
Ψ
e
τ
ϖ
τ
=
E
e
τ
ϖ
τ
T
Ψ
^
-
ƛ
P
0
0
0
e
τ
ϖ
τ
<
-
E
e
τ
ϖ
τ
T
ƛ
P
0
0
0
e
τ
ϖ
τwhich leads to the following conclusion: J
e
(
e
(
τ
)
,
ϖ
(
τ
)
)
<
0. This implies that the error system (20) via control inputs u
(
τ
) with the update law (19) is asymptotically stable. Then, (9) follows immediately from (37) and the proof of Definition 1 is completed.
From Lemma6, we get(42)
I
V
τ
≤
-
η
′
e
T
e
.It follows from Ito’s formula that (43)
E
V
τ
,
e
τ
-
E
V
0
,
e
0
=
∫
0
τ
E
I
V
≤
-
η
′
∫
0
τ
E
e
s
T
e
s
d
s
,
τ
>
0
.From Lemma 11, (43) becomes (44)
λ
m
i
n
P
E
e
τ
2
≤
E
V
τ
,
e
τ
≤
E
V
0
,
e
0
-
η
′
∫
0
τ
E
e
s
2
d
s
.Using Gronwall’s inequality (Lemma 10), we have (45)
E
e
τ
2
≤
ε
~
exp
⁡
-
ε
^
τ
,where ε
~
=
E
V
(
0
,
e
(
0
)
)
/
λ
m
i
n
(
P
) and ε
^
=
η
′
/
λ
m
i
n
(
P
)
.
Relation (45) guarantees the exponential synchronization. This completes the proof.Remark
2. It is important to show that the trajectories of the slave system are bounded under the feedback in the presence of perturbations: (46)
sup
τ
≥
0
⁡
x
^
1,3
≤
π
.
(R-1)
Since the systems are assumed to operate in chaotic mode without feedback, their trajectories converge to compact invariant set. Let one assume thatξ is the spectral radius of the perturbation on the system dynamics. We set r
+
ξ
>
0, and let the closed ball B
-
r
+
ξ strictly contains such a compact set (since the chaotic trajectories are bounded, we assume that they are contained in B
-
r
+
ξ); let ∞
>
T
≥
0 be the smallest number such that e
(
τ
)
∈
B
-
r
+
ξ
∀
τ
≥
T. The previous development shows that V is positive definite Lyapunov function with a negative definite derivative for any values of the states contained in a compact set. Hence, there always exist control gains such that the system under feedback is forward complete.
(R-2)
It may be shown that the system under control is forward complete; that is, if there exists a set of initial conditions and gains such that, together, they generate solutions that tend to infinity, these solutions may unboundedly grow only in infinite time. From this, it follows that for eachθ
>
0 and T
>
0, there exists M
(
θ
,
T
) such that e
(
0
)
≤
θ
+
ξ
⇒
e
(
τ
)
≤
M
(
θ
+
ξ
,
T
)
,
∀
τ
∈
0
,
T, where M is in general a nondecreasing function of its arguments. We assume that T is the synchronization time since, by assumption, the system operates in open loop for all τ
∈
0
,
T for any T
>
0
,
e
(
τ
)
≤
max
⁡
M
(
0
+
ξ
,
T
)
,
r
+
ξ for all τ
∈
0
,
T for any T
>
0. That is, the solutions are bounded. Note that the bound max
⁡
M
(
θ
+
ξ
,
T
)
,
r
+
ξ is independent of the gains, and for T
≥
T
∗, we can safely assume that max
⁡
M
(
0
+
ξ
,
T
)
,
r
+
ξ
=
r
+
ξ and π depends only on r; hence, (R-2) of the proposition follows. If T and θ
ξ
θ
ξ
=
m
a
x
(
θ
+
ξ
) are given and max
⁡
M
(
θ
+
ξ
,
T
)
,
r
+
ξ
=
M, then π depends on T and θ
ξ and so does K; hence, (R-1) of the proposition follows. In either case, K does not depend on the initial conditions e
(
0
), nor is it required that e
(
0
)
≤
r
+
ξ.


## 4.  Simulation Investigation

In this section, we will provide simulation results for a system where the master oscillator is the perturbed modified Colpitts oscillators and the observer is a perturbed and stochastic Chua system.From (1), the modified Colpitts systems whose dimensionless equations are given in [29] is represented by the values as follows: (47)
x
=
x
1
x
2
x
3
,
A
=
0
0
1
0
-
b
0
1
-
1
-
1
-
c
11
,
B
=
0
-
b
0
1
,
D
=
-
1
0
0
0
0
0
0
0
0
,
φ
x
=
φ
x
2
,
x
3
0
0
.The slave based on the observer method is the Chua’s circuit and can be generalized in form (2) [30], where (48)
F
=
-
α
α
0
0
-
1
1
-
1
-
β
1
-
β
2
,
G
=
-
α
0
0
0
1
0
0
0
1
,
f
x
^
=
f
x
^
1
0
0
.
φ
(
x
2
,
x
3
) and f
(
x
^
1
) are continuous-time functions defined, respectively, by φ
(
x
2
,
x
3
)
=
a
2
exp
⁡
(
-
a
x
3
-
b
x
2
) and f
(
x
^
1
)
=
b
′
x
^
1
+
0.5
a
′
-
b
′
(
x
^
1
+
1
-
x
^
1
-
1
) with a
′
<
b
′
<
0.If the parameters are taken asα
=
10, β
1
=
15, β
2
=
0.0385, a
′
=
-
1.28, b
′
=
-
0.69, a
=
2.251362, b
=
192.3, b
0
=
0.1064814815, c
11
=
0.934, and a
2
=
8.518518
×
1
0
-
11, the systems specified by (46) and (47) exhibit chaotic dynamics. Figures 1 and 2 give the time history and phase portrait of both systems.Figure 1
Time history and phase portraits of modified Colpitts oscillator.Figure 2
Time history and phase portraits of Chua’s circuit.Remark
3. Systems (1) and (2) are quite different to make our illustration general. Most previous studies in the literature have predominantly concentrated on standard systems such as the Lorentz, the classic Colpitts oscillator, the Chua system, the Chen system, the Lu system, or the Rössler system either in the studies of their stability analysis and periodic oscillations or in the studies of their synchronization (the general form of these dynamical systems is provided in (1)). It has been shown that the modified Colpitts oscillator presents different topology and can exhibit complicated dynamics with reference to circuits previously mentioned. If B
≠
0, the problem becomes complex and requires some supplementary conditions by comparison to those employed in the literature; this aspect makes the present result more general with respect to those encountered in the literature.For convenience, we select the following matrices:(49)
C
=
Λ
-
1
0
-
1
0
0
0
0
0
0
,
L
=
1
0
0
0
1
0
0
0
1
,where Λ
≠
0.If we letΛ
=
-
1
0
-
2, then the estimated feedback gain is updated according to the following adaptation algorithm: (50)
K
˙
=
ρ
x
1
-
x
^
1
2
+
x
3
-
x
^
3
2
.Since structural variations may take any forms in theory depending on the sources of disturbances, we illustrate the proposed control scheme with a more global example where many perturbation sources are taken into consideration [29–31] under the assumption that Δ
(
•
)
(
τ
)
≤
δ. Hence we choose an arbitrary structural variation as follows:(51)
Δ
A
=
ϕ
τ
0.1
0
-
0.3
0
-
0.5
0
0.4
0
0.8
;
Δ
F
0
=
ϕ
τ
0.8
0
-
2.4
0
-
4
0
3.2
0
6.4
,where ϕ
(
τ
)
∈
-
1
,
1.Let us select the perturbed matrix simply in the following form:(52)
Δ
B
τ
=
ε
1
τ
ε
2
τ
0
.We adopt ε
1
(
τ
) as an additive Gaussian noise such that ε
1
(
τ
)
≤
1.75 and ε
2
(
τ
) is estimated as ε
2
(
τ
)
=
0.8
sin
⁡
8
τ
+
0.3
sin
⁡
(
12
τ
). The noise perturbation can be taken as σ
(
τ
,
e
(
τ
)
)
=
diag
⁡
0.002
e
x
(
τ
)
,
0.065
e
y
(
τ
)
,
0.042
e
z
(
τ
) and ω
(
τ
) is a 3-dimensional Brownian motion satisfying E
d
ω
(
τ
)
=
0 and E
d
ω
(
τ
)
2
=
d
τ. The initial values of master system and the slave system are adopted as (
x
10
,
x
20
,
x
30
)
=
(
2
×
1
0
-
5
,
2
×
1
0
-
5
,
2
×
1
0
-
6
) and x
^
10
,
x
^
20
,
x
^
30
=
8
×
1
0
-
3
,
6
,
5, respectively. In addition, we select ρ
=
0.5 and K
(
0
)
=
1.28 which lead to γ
=
0.80.Figure3 shows the time response of the synchronization errors without the adaptive robust feedback controller u
(
τ
) designed in (18). From these figures, one can see that the synchronization errors diverge. Figure 4 depicts the output errors with the robust adaptive feedback controller u
(
τ
), from which we can see that the master system (3) and slave system (4) can be exponentially synchronized in mean square. From the theoretical analysis, it appears that the complexity of perturbations cannot give any information about the performance of the system. It is hard to obtain an optimal value of feedback control gain. In order to explore the synchronization behavior of the scheme, K is taken as a variable of the output error (we recall that the initial value of the system variables does not affect the evolution of the feedback gain) to see how the solution of the output error evolves with the variation of the parameter. From Figure 5, it can be observed that y
e
x and y
e
z converge to zero at about K
=
1.827 which means that the synchronization can be achieved when K
>
1.827.Figure 3
The time history of the output error with disturbances when no control input is applied.Figure 4
The time history of the output errors when the control input is applied in presence of structural perturbation.Figure 5
Time evolution of the estimated feedback gainK in presence of structural perturbation.


## 5. Conclusion

Because the adaptive control has some interesting features such as low sensitivity to external disturbances, robustness to the plant uncertainties, and easy realization, in this paper, we use this method to realize exponential synchronization of two different uncertain chaotic systems in which the slave system is noise perturbed. By employing the Lyapunov functional method and adaptive control, several sufficient conditions have been obtained which ensure the coupled chaotic systems to be exponentially robustly synchronized in the mean square. Furthermore, the theoretical analysis is easily verified by using the standard numerical software. We have selected two perturbed systems consisting of modified Colpitts oscillator and Chua’s oscillator. It was found that the controller maintains robust stable synchronization in the presence of exoteric perturbations, structural uncertainties, and noises.



---

*Source: 535317-2015-09-29.xml*