# State Observer Synchronization Used in the Three-Dimensional Duffing System

Authors: Jian-Qun Han; Yue-Chao Ma; Hong Sun
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Synchronization of chaotic systems has attracted extensive concern in the past few years. In this study, we investigate a new structure of Duffing system by the variable decomposition method. Then, we analyze the state observer synchronization based on the new Duffing system. It is proved theoretically that the designed observer can keep synchronization with Duffing chaotic system in transmitter. The design is presented reasonably with the conditional Lyapunov exponents, and its effectiveness is clearly shown in simulation results.


---

## Body

## 1. Introduction

Duffing equations are well-known nonlinear equations of applied science which are used as a powerful tool to discuss some important practical phenomena such as periodic orbit extraction, nonuniformity caused by an infinite domain, and nonlinear mechanical oscillators.Synchronization is a very common phenomenon of nonlinear dynamic system. It shows that two moving objects can make their states converge through the medium strong or weak coupling. Since Pecora and Carroll [1] put forward the principle of chaotic synchronization, synchronization of chaotic systems has become one of research hotspots in the field of nonlinear dynamics. Chaos synchronization has been widely applied to many fields such as biology, chemistry, medicine, electronics, information science, and secure communication. So far, there are various synchronous types [2] including complete synchronization, lag synchronization, anticipated synchronization, and generalized synchronization. Furthermore, many synchronous implementation methods have been studied which included drive-response synchronization, adaptive synchronization [3], impulsive synchronization [4], observer synchronization [5], and the sliding mode control [6]. The above methods are commonly used in data-driven fault diagnosis [7–11].For obtaining synchronization, both the structure and the parameters of target system are often supposed to be available, while the only unknown is the initial condition. As the model of a system is commonly imprecise, obviously, it is not necessarily a realistic case. Therefore, the synchronization or autosynchronization of chaotic systems with unknown parameters has been widely studied in recent years. The synchronization of Lorenz system, one of the most popular chaotic systems, is stressed, and the conditions of synchronization have been studied by various researchers [12–20]. The main strategy is based on the adaptive feedback control. It is illustrated that synchronization is achievable, while the global convergence is guaranteed by the Lyapunov convergence theorem [12–16]. However, the proof of global convergence poses a restrictive form of the system. Thus, the unknown parameters can only appear when the system has to be transformed into some particular forms, such as state-affine form, or in the dynamical equation of the observable state.Due to easy engineering implementation, the observer synchronization method has gained extensive attention. Observer synchronization has simple structure, fast synchronous speed, and high precision. The method can be widely used in communications and information processing.At present, the study of Duffing equation is in the form of two-dimensional state equation [21–24]. Three-dimensional Duffing state equation is presented in this paper. We combined the state observer with chaotic synchronization and proved the rationality of the nonlinear synchronous observer.This paper is organized as follows. In Section2, three-dimensional Duffing state equation is presented. The phase diagrams of new Duffing system are also given. In Section 3, the basic idea of the adaptive observer design is described in detail. With supportive simulation results, the conditional Lyapunov exponents of the system are then studied in Section 4. Finally, conclusions are drawn, respectively, in Section 5.


## 2. Three-Dimensional Duffing System

There are many problems in engineering, such as packaging systems based on displacement excitation of nonlinear vibrations, pressure sensors nonlinear vibration. They can be simplified into a forced Duffing equation with cubic non-linearity, governed by(1)
x
¨
(
t
)
+
μ
x
˙
(
t
)
+
ω
0
2
x
(
t
)
+
ε
x
(
t
)
3
=
F
cos
⁡
Ω
t
,
where x
(
t
) is an unknown real function; μ, ω
0, ε, F, Ω are known physical parameters; and the dot represents differentiation with respect to t.Definition 1.
The particular form of Duffing system related to system (1) is described as
(2)
x
˙
1
=
x
2
,
x
˙
2
=
x
1
-
x
1
3
-
k
x
2
+
r
cos
⁡
ω
t
,
where r and ω are the amplitude and frequency of driving force, respectively, k is damping ratio, and x
1
-
x
1
3 is the nonlinear restoring force.
When given the system initial statex
1
(
0
)
=
0, x
2
(
0
)
=
1 in (2), two-dimensional Duffing system under sinusoidal signal drive produces the chaotic phase diagram as shown in Figure 1.Figure 1
Two-dimensional Duffing chaotic system phase diagram withx
1
(
0
)
=
0 and x
2
(
0
)
=
1.Definition 2.
Given a constantb, system (2) can be converted into state equation
(3)
x
˙
1
=
x
2
,
x
˙
2
=
-
k
x
2
-
x
3
+
b
+
r
cos
⁡
⁡
ω
t
,
x
˙
3
=
3
x
1
2
x
2
-
x
2
.Proof.
Choose the following differential equation:(4)
x
˙
3
=
3
x
1
2
x
2
-
x
2
.
Intergrating x
˙
3 taken in (4) for system (2) yields
(5)
x
3
=
x
1
3
-
x
1
+
c
,
where c is undetermined constant, and the value can be determined by the system initial state.Remark 3.
From the above analysis, it follows thatc and b are unbounded, but b is used to counteract the effects of c. When given x
1
(
0
)
=
0 and x
3
(
0
)
=
0 in the system (3), c
=
0 and b
=
0 hold. The three-dimensional Duffing chaotic dynamic system (3) is the same as system (2). Under the general conditions, system (2) can be modified to system (3) by adjusting variable b to offset the impact of variable c.
When givenx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
0, and b
=
0 in (3), three-dimensional Duffing system produces the chaotic phase diagram shown in Figure 2. Its projections on the xoy, yoz, and xoz coordinate plane are shown in Figures 3, 4, and 5. By comparing Figures 1 and 3, we can obtain the result that two-dimensional Duffing system and three-dimensional Duffing system are consistent. Figures 4 and 5 can only produce in three-dimensional Duffing system.
The simulation results of equivalent condition are the premise of system initial statesx
1
(
0
)
=
0 and x
3
(
0
)
=
0; otherwise the system status may not be chaotic. When given x
1
(
0
)
=
0, x
2
(
0
)
=
1, and x
3
(
0
)
=
2, simulation result is shown in Figure 6. The system state is not chaotic.
When it is assumed thatb
=
2 in (3), the system phase diagram as shown in Figure 7 is chaotic. From the above analysis, results counteracting the effects of initial condition on system state and keeping the system chaotic can be obtained by adjusting the variable b.Figure 2
Three-dimensional Duffing chaotic system phase diagram withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
0, and b
=
0.Figure 3
Three-dimensional Duffing chaotic system phase diagram in the xoy plane projection withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
0, and b
=
0.Figure 4
Three-dimensional Duffing chaotic system phase diagram in the yoz plane projection withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
0, and b
=
0.Figure 5
Three-dimensional Duffing chaotic system phase diagram in the xoz plane projection withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
0, and b
=
0.Figure 6
Three-dimensional Duffing chaotic system phase diagram in the xoy plane projection withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
2, and b
=
0.Figure 7
Three-dimensional Duffing chaotic system phase diagram in the xoy plane projection withx
1
(
0
)
=
0, x
2
(
0
)
=
1, x
3
(
0
)
=
2, and b
=
2.


## 3. State Observer Synchronization

Definition 4.
Suppose that there exist(6)
A
=
[
0
1
0
0
-
k
-
1
1
1
-
1
]
,
B
=
[
1
1
1
]
,
E
=
[
0
1
0
]
such that
(7)
x
˙
=
[
0
1
0
0
-
k
-
1
1
1
-
1
]
x
+
[
1
1
1
]
×
[
0
b
3
x
1
2
x
2
-
x
1
-
2
x
2
+
x
3
]
+
[
0
1
0
]
r
cos
⁡
⁡
ω
t
.Proof.
For (7), written in general form as follows:
(8)
x
˙
=
A
x
+
B
u
+
E
w
,
y
=
C
T
x
,
u
=
g
(
x
)
,
where x, u, and w are, respectively, the system state variables, the nonlinear state feedback input, and system driver, x
∈
R
n, u
∈
R
k, and w
∈
R
m; y
∈
R
p is the output of sending system; coefficients are n
=
3, k
=
3, m
=
1, and p
=
1.
There existrank
⁡
E
=
m, rank
⁡
C
T
=
p, and rank
⁡
C
T
E
=
m in three-dimensional Duffing system. Notice that p
≥
m here.
At the receiving system, when the nonlinear input feedback can be measured, system (8) for state observer can be defined in (9). Consider
(9)
z
˙
=
G
^
z
+
F
^
u
+
J
^
y
,
x
^
=
z
-
H
^
y
,
where z
∈
R
n and x
^
∈
R
n; G
^, F
^, J
^, and H
^ are the unknown constant matrix.
Now definee
=
x
^
-
x
=
z
-
(
I
n
+
H
^
C
T
)
x and N
^
=
I
n
+
H
^
C
T. Then, we obtain the dynamic equation of the error system that
(10)
e
˙
=
G
^
e
+
(
J
^
C
T
-
N
^
A
+
G
^
N
^
)
x
+
(
F
^
-
N
^
B
)
u
-
N
^
E
w
.
If system (10) has
(11)
(
I
n
+
H
^
C
T
)
E
=
0
,
F
^
-
N
^
B
=
0
,
J
^
C
T
-
N
^
A
+
G
^
N
^
=
0
,
then it can be rewritten as
(12)
e
˙
=
G
^
e
.Remark 5.
If all eigenvalues ofG
^ are less than zero, then x
^ will exponentially converge at x.
ByJ
^
C
T
-
N
^
A
+
G
^
N
^
=
0, we can obtain
(13)
G
^
=
N
^
A
-
(
J
^
+
G
^
H
^
)
C
T
.
Let M
^
=
J
^
+
G
^
H
^. Then, we can obtain
(14)
J
^
=
M
^
+
M
^
C
T
H
^
-
N
^
A
H
^
,
which leads to
(15)
z
˙
=
(
N
^
A
-
M
^
C
T
)
z
+
F
^
u
+
J
^
y
.
So the state observer design problem turns into finding the rightH
^ and M
^. To make it meet (
I
n
+
H
^
C
T
)
E
=
0 and eig
(
G
^
)
<
0.
From (15), we can obtain
(16)
G
^
=
N
^
A
-
M
^
C
T
.Remark 6.
Alleig
(
G
^
)
<
0, if ∃ for all M
^
∈
R
n
×
p and B
=
{
x
∣
x
=
eig
(
N
^
A
-
M
^
C
T
)
,
for
all
x
∈
R
}. Namely, the eigenvalues of N
^
A
-
M
^
C
T can be set arbitrarily if (
N
^
A
,
C
T
) can be completely observable.Proof.
By (11), we can get H
^
C
T
E
=
-
E and rank
⁡
C
T
E
=
m. So we can obtain
(17)
H
^
=
-
E
(
C
T
E
)
+
+
Y
(
I
p
+
(
C
T
E
)
(
C
T
E
)
+
)
,
where Y
∈
R
n
×
p, n, p
∈
R; Y is the right dimension of arbitrary matrix.
Consider the expressionN
^
=
I
n
+
H
^
C
T; we can obtain
(18)
N
^
=
(
I
n
+
Y
C
T
)
(
I
n
-
E
(
C
T
E
)
+
C
T
)
.Theorem 7.
Assumerank
⁡
N
^
=
n
-
m if |
I
n
+
Y
C
T
|
≠
0 and |
I
p
+
C
T
Y
|
≠
0. One can obtain
(19)
rank
⁡
[
C
T
C
T
N
^
A
⋮
C
T
(
N
^
A
)
n
-
1
]
=
n
which yields that (
N
^
A
,
C
T
) can be completely observable and the eigenvalues of G
^ can be arbitrarily set.By Theorem7, we know that |
I
p
+
C
T
Y
|
≠
0 if rank
⁡
C
=
p and Y
=
0. Then, we have
(20)
N
^
=
I
n
-
E
(
C
T
E
)
+
C
T
.Remark 8.
WhenN
^
=
I
n
-
E
(
C
T
E
)
+
C
T, we can always find a suitable M
^, making eig
(
G
^
)
<
0. Namely, x
^
→
x for ∀
t
→
∞.
Now let(21)
A
=
[
0
1
0
0
-
0.5
-
1
1
1
-
1
]
,
B
=
[
1
1
1
]
,
E
=
[
0
1
0
]
,
C
=
[
1
1
0
]
,
u
(
t
)
=
[
0
b
3
x
1
2
x
2
-
x
1
-
2
x
2
+
x
3
]
,
w
(
t
)
=
r
cos
⁡
⁡
ω
t
in
(
7
)
.
From the above analysis, we can obtain
(22)
H
^
=
-
E
(
C
T
E
)
+
=
[
0
-
1
0
]
,
N
^
=
I
n
+
H
^
C
T
=
[
1
0
0
-
1
0
0
0
0
1
]
if Y
=
0 in (17).
LetM
^
=
[
1
1
1
]. Then we obtain G
^
=
[
-
1
0
0
-
1
-
2
0
0
0
-
1
] and eig
(
G
^
)
=
(
-
2
,
-
1
,
-
1
). All eigenvalues of G
^ are negative. Then x
^
→
x as t
→
∞. So we can obtain the receiving system:
(23)
z
˙
=
[
-
1
0
0
-
1
-
2
0
0
0
-
1
]
z
+
[
1
0
0
-
1
0
0
0
0
1
]
u
(
t
)
+
[
1
-
1
1
]
y
,
x
^
=
z
-
[
0
-
1
0
]
y
.


## 4. Simulation and Analysis

In this section, we evaluate the performance of state observer synchronization used in three-dimensional Duffing system. Under the condition of interference and no interference, we simulated transceiver system synchronization.To discuss the synchronization performance of system, we assumed coefficient of transmitting systemk
=
0.5, b
=
2, r
=
0.7, and ω
=
1 in (7). The synchronous state observer is designed as a receiving system. The simulation system based on MATLAB is shown in Figure 8.Figure 8
State observer synchronization simulation system used in three-dimensional Duffing system.Under the condition of no interference, synchronization error is shown in Figure9. There exist e
1
=
x
1
-
x
^
1, e
2
=
x
2
-
x
^
2, and e
3
=
x
3
-
x
^
3. Error curves tend to zero. Simulation results show that transceiver system can keep the good synchronization performance.The error curves driven by the outside periodic force0.7
cos
⁡
⁡
t and the damping ratio fixed to 0.5. The parameter b can be modified. There is no interference.
(a)
The curvee
1 for x
1
-
x
^
1
(b)
The curvee
2 for x
2
-
x
^
2
(c)
The curvee
3 for x
3
-
x
^
3When there is interference, we assume adding Gaussian white noise tou
(
t
) or w
(
t
) in (8). Assume noise power spectrum density n
0
=
0.1. Error curves tend to zero too. The simulation results are similar to the Figure 9 without interference.


## 5. Conclusions

The present work has studied a new three-dimensional chaotic synchronization of Duffing system. The chaotic system is different from the previous Duffing system. It contains three equations. The third equation describes nonlinear effect on the new Duffing system. This system redefined state variables, nonlinear feedback input, and driving force. Because the synchronization method based on observer does not need the initial state in the same basin of attraction, we prove observer synchronization in three-dimensional Duffing system. By adopting the method of simple linear observer, receiving system can successfully keep synchronization with transmitting chaotic system.



---

*Source: 986359-2014-04-03.xml*