# Synchronization of the Coupled Distributed Parameter System with Time Delay via Proportional-Spatial Derivative Control

Authors: Kun Yuan; Abdulaziz Alofi; Jinde Cao; Abdullah Al-Mazrooei; Ahmed Elaiw
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

By combining parabolic partial differential equation (PDE) theory with Lyapunov technique, the synchronization is studied for a class of coupled distributed parameter systems (DPS) described by PDEs. First, based on Kronecker product and Lyapunov functional, some easy-to-test sufficient condition is given to ensure the synchronization of coupled DPS with time delay. Secondly, in the case that the whole coupled system cannot synchronize by itself, the proportional-spatial derivative (P-sD) state feedback controller is designed and applied to force the network to synchronize. The sufficient condition on the existence of synchronization controller is given in terms of a set of linear matrix inequalities. Finally, the effectiveness of the proposed control design methodology is demonstrated in numerical simulations.


---

## Body

## 1. Introduction

Most practical systems are distributed in space and time, for example, systems related to heat flows, fluid flow, or flexible structure. The systems are called distributed parameter systems (DPSs), which are mathematically modelled by partial differential equations (PDEs) with boundary conditions. The key characteristic of DPSs is that their inputs, outputs, process states, and the relevant parameters may vary temporally as well as spatially. Due to such characteristic and infinite dimensional, the distributed parameter system is difficult to control. Recent research [1–12] show that the dynamical behavior of the parabolic PDEs can be described approximately in an array of low-order ordinary differential equations (ODEs). In [3], by integrating the Galerkin method with neural networks, the author proposes a simple and effective modeling method for DPS. In [13, 14], the K-L method is applied to model the distributed parameter system. In [4], by the Galerkin’s method and the geometric control, the authors stabilize DPSs. For nonlinear ODE system, fuzzy-model-based control technique is conceptually simple and effective. Recently, it has been successfully applied to fuzzy control design of nonlinear DPSs [15–20]. Despite these results, the control of distributed parameter systems is a broad area of research, and further results are demanded to design an effective control method, which motivates our study.The coupled system can be regarded as a large set of interconnected nodes, which can be expressed by the graph. Each node of the graph represents individual, and the edge of the graph represents the connections among them. The dynamical analysis of coupled system has become a focal topic of great interest, particularly the synchronization phenomena. Synchronization in an array of coupled dynamical systems was first investigated in [21]. Later, some results on synchronization in various coupled systems have been given in [22–26]. It is noted that the dynamical behavior of a coupled system is determined not only by the dynamical behaviors of the isolated nodes, but also by the edges, which depends on the coupled topology of the system. As a special case of coupled system, coupled distributed parameter systems have been an interest. If the coupled DPS is regarded as the special graph, the nodes of the network are linear DPS, the edges of the network are parabolic spatial differential operator with diffusive coupling. In other words, if every node is thermal process, the edges represent the thermal transmittance among the nodes. Synchronization of coupled system has been investigated by some researchers [27], and some synchronization criteria are given.In the case that the whole coupled system cannot synchronize by itself, some controllers should be designed and applied to force the network to synchronize. Due to the fact that the outputs, inputs, process states, and the relevant parameters of DPS may vary temporally as well as spatially, a special control strategy called proportional-spatial derivative (P-sD) control [17] is designed to achieve synchronization of coupled DPS, which is similar as the traditional PD (proportional-time derivative) control.Motivated by the above discussion, the aim of this paper is to synchronize coupled distributed parameter systems with time delay. First, based on Kronecker product and Lyapunov functional, some easy-to-test sufficient condition is given to ensure the synchronization of coupled DPS with time delay. Secondly, in the case that the whole coupled system cannot synchronize by itself, the proportional-spatial derivative (P-sD) state feedback controller is designed and applied to force the network to synchronize. The sufficient condition on the existence of synchronization controller is given in terms of a set of linear matrix inequalities. Finally, a numerical example is given to show effectiveness of the proposed method.Notation. Throughout this paper, for real symmetric matrices X and Y, the notation X
≥
Y (resp., X
>
Y) means that the matrix X
-
Y is positive semidefinite (resp., positive definite). The superscript “T” represents the transpose. Matrix dimensions, if not explicitly stated, are assumed to be compatible for algebraic operations.


## 2. Model Description and Preliminaries

Consider a dynamical network consisting ofN identical and diffusively coupled nodes, with each node being an n-dimensional delayed distributed parameter system. The state equations of the network are
(1)
∂
T
k
∂
t
=
D
∂
2
T
k
∂
z
2
+
A
T
k
(
t
,
z
)
+
B
T
k
(
t
-
τ
,
z
)
+
∑
j
=
1
,
j
≠
k
N
‍
G
k
j
∂
2
(
T
j
-
T
k
)
∂
z
2
,
h
h
z
∈
(
z
1
,
z
2
)
,
k
=
1,2
,
…
,
N
subject to the boundary conditions
(2)
∂
T
k
(
t
,
z
)
∂
t
|
z
=
z
1
=
0
,
∂
T
k
(
t
,
z
)
∂
t
|
z
=
z
2
=
0
and the initial condition
(3)
T
k
(
t
,
z
)
=
ϕ
k
(
t
,
z
)
,
t
∈
[
-
τ
,
0
]
,
where z
∈
[
z
1
,
z
2
] is the spatial coordinate, T
k
(
z
,
t
)
:
[
z
1
,
z
2
]
×
ℝ
→
ℝ
n denotes state variables, k
=
1,2
,
…
,
N is positive constant, and τ
>
0 denotes the state delay; the vector function f
(
·
) and g
(
·
) are nonlinear locally Lipschitz continuous functions. G
=
(
G
i
j
)
N
×
N is the coupling configuration matrix representing the topological structures of the networks which satisfies the diffusive coupling connection:
(4)
G
i
j
≥
0
,
i
≠
j
,
G
i
i
=
-
∑
j
=
1
,
j
≠
i
N
‍
G
i
j
,
which ensures the diffusion that ∑
j
=
1
N
‍
G
i
j
=
0.Equivalently, system (1) can be rewritten in a form as follows:
(5)
∂
T
k
∂
t
=
D
∂
2
T
k
∂
z
2
+
A
T
k
(
t
,
z
)
+
B
T
k
(
t
-
τ
,
z
)
+
∑
j
=
1
N
‍
G
k
j
∂
2
T
j
∂
z
2
,
hhhhh
z
∈
(
z
1
,
z
2
)
,
k
=
1,2
,
…
,
N
.LetT
(
t
,
z
) be a function to which all T
k
(
t
,
z
) are expected to synchronize and T
(
t
,
z
) satisfies the following equation:
(6)
∂
T
∂
t
=
D
∂
2
T
∂
z
2
+
A
T
(
t
,
z
)
+
B
T
(
t
-
τ
,
z
)
,
z
∈
(
z
1
,
z
2
)
subject to the boundary conditions
(7)
∂
T
(
t
,
z
)
∂
t
|
z
=
z
1
=
0
,
∂
T
(
t
,
z
)
∂
t
|
z
=
z
2
=
0
and the initial condition
(8)
T
(
t
,
z
)
=
ϕ
(
t
,
z
)
,
t
∈
[
-
τ
,
0
]
.Our objective is to investigate the synchronization of network (1) with respect to the isolated node T
(
t
,
z
). Denote e
i
(
t
,
z
)
=
T
i
(
t
,
z
)
-
T
(
t
,
z
); then the following error dynamical system is obtained:
(9)
∂
e
k
∂
t
=
D
∂
2
e
k
∂
z
2
+
A
e
k
(
t
,
z
)
+
B
e
k
(
t
-
τ
,
z
)
+
∑
j
=
1
N
‍
G
k
j
∂
2
e
j
∂
z
2
,
h
h
h
h
h
h
h
z
∈
(
0
,
l
)
,
k
=
1,2
,
…
,
N
.
The error system is subject to the boundary conditions and initial values as follows:
(10)
∂
e
k
(
t
,
z
)
∂
t
|
z
=
z
1
=
0
,
∂
e
k
(
t
,
z
)
∂
t
|
z
=
z
2
=
0
and the initial condition
(11)
e
k
(
t
,
z
)
=
ϕ
k
(
t
,
z
)
-
ϕ
(
t
,
z
)
,
t
∈
[
-
τ
,
0
]
.
Suppose e
(
t
,
z
)
=
[
e
1
T
(
t
,
z
)
,
e
2
T
(
t
,
z
)
,
…
,
e
N
T
(
t
,
z
)
]
T
∈
R
n
N; then the error system (9) can be written as
(12)
∂
e
∂
t
=
(
I
N
⊗
D
)
∂
2
e
∂
z
2
+
(
I
N
⊗
A
)
e
(
t
,
z
)
+
(
I
N
⊗
B
)
e
(
t
-
τ
,
z
)
+
(
G
⊗
I
n
)
∂
2
e
∂
z
2
,
z
∈
(
z
1
,
z
2
)
,
where ⊗ denotes the Kronecker product of the two matrices.Definition 1.
The coupled distributed parameter system in (1) is asymptotically synchronous if the synchronization error satisfies lim
⁡
t
→
∞
e
k
(
t
,
z
)
=
0 for all z
∈
[
z
1
,
z
2
] and k
=
1,2
,
…
,
N.
LetD
-
=
I
N
⊗
D
+
G
⊗
I
n, A
-
=
I
N
⊗
A and B
-
=
I
N
⊗
B; the system (12) can be rewritten as follows:
(13)
∂
e
∂
t
=
D
-
∂
2
e
∂
z
2
+
A
-
e
(
t
,
z
)
+
B
-
e
(
t
-
τ
,
z
)
,
A
-
∈
ℝ
n
N
×
n
N,B
-
∈
ℝ
n
N
×
n
N.Theorem 2.
For the coupled DPS (5), synchronize to the system (6) if there are matrices P
=
P
T
>
0
,
Q
=
Q
T
>
0, satisfying the following linear matrix inequalities:
(14)
Ω
=
[
P
A
-
+
A
-
T
P
+
Q
P
B
-
0
B
-
T
P
-
Q
0
0
0
-
2
P
D
-
]
<
0
.
Consider the following Lyapunov functional candidate for system
(15)
V
(
t
)
=
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
e
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
∫
t
-
τ
t
‍
e
T
(
θ
,
z
)
Q
e
(
θ
,
z
)
d
θ
d
z
,
d
V
d
t
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
e
t
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
[
e
T
(
t
,
z
)
Q
e
(
t
,
z
)
-
e
T
(
t
-
τ
,
z
)
h
h
h
h
h
h
×
Q
e
(
t
-
τ
,
z
)
e
T
]
d
z
+
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
[
D
-
e
z
z
(
t
,
z
)
d
z
+
A
-
e
(
t
,
z
)
hhhhhhhhhhhhhhhh
+
B
-
e
(
t
-
τ
,
z
)
]
d
z
+
∫
z
1
z
2
‍
[
e
T
(
t
,
z
)
Q
e
(
t
,
z
)
hhhhhhh
-
e
T
(
t
-
τ
,
z
)
Q
e
(
t
-
τ
,
z
)
]
d
z
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
D
-
e
z
z
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
e
T
(
t
,
z
)
[
2
P
A
-
+
Q
]
e
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
e
T
(
t
,
z
)
2
P
B
-
e
(
t
-
τ
,
z
)
d
z
-
∫
z
1
z
2
‍
e
T
(
t
-
τ
,
z
)
Q
e
(
t
-
τ
,
z
)
d
z
.
From integrating by parts and taking into account (10), we have
(16)
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
D
-
e
z
z
(
t
,
z
)
d
z
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
D
-
d
(
e
z
(
t
,
z
)
)
=
2
e
T
(
t
,
z
)
P
D
-
e
z
(
t
,
z
)
|
z
1
z
2
-
2
∫
z
1
z
2
‍
e
z
T
(
t
,
z
)
P
D
-
e
z
(
t
,
z
)
d
z
=
-
2
∫
z
1
z
2
‍
e
z
T
(
t
,
z
)
P
D
-
e
z
(
t
,
z
)
d
z
.
Using (26) and (16), we find that
(17)
d
V
d
t
=
∫
z
1
z
2
‍
ξ
T
(
t
,
z
)
Ω
ξ
(
t
,
z
)
d
z
,
where ξ
(
t
,
z
)
=
[
e
T
(
t
,
z
)
,
e
T
(
t
-
τ
,
z
)
,
e
z
T
(
t
,
z
)
]
T and
(18)
Ω
=
[
P
A
-
+
A
-
T
P
+
Q
P
B
-
0
B
-
T
P
-
Q
0
0
0
-
2
P
D
-
]


## 3. Synchronization Control

If the coupled distributed parameter system cannot synchronize to the isolate node (6), the controller will be designed. Consider
(19)
∂
T
k
∂
t
=
D
∂
2
T
k
∂
z
2
+
A
T
k
(
t
,
z
)
+
B
T
k
(
t
-
τ
,
z
)
+
∑
j
=
1
,
j
≠
k
N
‍
G
k
j
∂
2
(
T
j
-
T
k
)
∂
z
2
+
C
u
(
t
,
z
)
,
h
h
h
h
h
h
h
h
z
∈
(
z
1
,
z
2
)
,
k
=
1,2
,
…
,
N
,
(20)
∂
e
∂
t
=
D
-
∂
2
e
∂
z
2
+
A
-
e
(
t
,
z
)
+
B
-
e
(
t
-
τ
,
z
)
+
C
-
u
(
t
,
z
)
,
i
∈
S
≜
1,2
,
…
,
r
,
where C
-
=
I
n
⊗
C, and
(21)
u
(
t
,
z
)
=
K
e
(
t
,
z
)
+
L
e
z
(
t
,
z
)
,
where K and L are m
×
n
N gain matrices to be determined.With the control law (21), the overall close-loop system can be written as
(22)
e
t
(
t
,
z
)
=
D
-
e
z
z
(
t
,
z
)
+
(
A
-
+
C
-
K
)
e
(
t
,
z
)
+
B
-
e
(
t
-
τ
,
z
)
+
C
-
L
e
z
(
t
,
z
)
.Next, the asymptotical stability will be analyzed for the closed loop system (22) and the controller will be designed.Theorem 3.
For the coupled DPS (19), synchronize to the system (6) if there are matrices X
=
X
T
>
0, S
=
S
T
>
0, Y and R satisfying the following linear matrix inequalities:
(23)
Ω
-
=
[
A
-
X
+
X
A
-
T
+
C
-
Y
+
Y
T
C
-
+
S
B
-
X
C
-
R
*
-
S
0
*
0
-
2
D
-
X
]
<
0
,
where * represents blocks that are readily inferred by symmetry. Then the controller gains can be constructed as
(24)
K
=
Y
X
-
1
,
L
=
R
X
-
1
.Proof.
Consider the following Lyapunov functional candidate for system:(25)
V
(
t
)
=
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
e
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
∫
t
-
τ
t
‍
e
T
(
θ
,
z
)
Q
e
(
θ
,
z
)
d
θ
d
z
,
where P
=
P
T
>
0 and Q
=
Q
T
>
0. Consider
(26)
d
V
d
t
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
e
t
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
[
e
T
(
t
,
z
)
Q
e
(
t
,
z
)
-
e
T
(
t
-
τ
,
z
)
hhhhhhh
×
Q
e
(
t
-
τ
,
z
)
e
T
]
d
z
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
[
D
-
e
z
z
(
t
,
z
)
d
z
+
(
A
-
+
C
-
K
)
e
(
t
,
z
)
hhhhhhhhhhhhhh
+
B
-
e
(
t
-
τ
,
z
)
+
C
-
L
e
z
(
t
,
z
)
]
d
z
+
∫
z
1
z
2
[
e
T
(
t
,
z
)
Q
e
(
t
,
z
)
-
e
T
(
t
-
τ
,
z
)
‍
hhhhhhh
×
Q
e
(
t
-
τ
,
z
)
e
T
]
d
z
=
2
∫
z
1
z
2
‍
e
T
(
t
,
z
)
P
D
-
e
z
z
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
e
T
(
t
,
z
)
[
2
P
(
A
-
+
C
-
K
)
+
Q
]
e
(
t
,
z
)
d
z
+
∫
z
1
z
2
‍
e
T
(
t
,
z
)
2
P
B
-
e
(
t
-
τ
,
z
)
d
z
+
∫
z
1
z
2
‍
e
T
(
t
,
z
)
2
P
C
-
L
e
z
(
t
-
τ
,
z
)
d
z
-
∫
z
1
z
2
‍
e
T
(
t
-
τ
,
z
)
Q
e
(
t
-
τ
,
z
)
d
z
.
From the proof of Theorem 2, we can derive that
(27)
d
V
d
t
=
∫
z
1
z
2
‍
ξ
T
(
t
,
z
)
Ω
^
ξ
(
t
,
z
)
d
z
<
0
,
where ξ
(
t
,
z
)
=
[
e
T
(
t
,
z
)
,
e
T
(
t
-
τ
,
z
)
,
e
z
T
(
t
,
z
)
]
T
(28)
Ω
^
=
[
P
(
A
-
+
C
-
K
)
+
(
A
-
+
C
-
K
)
T
P
+
Q
P
B
-
P
C
-
L
B
-
T
P
-
Q
0
L
T
C
-
T
P
0
-
2
P
D
-
]
<
0
.
Next, we prove that (28) holds if LMIs (23) are satisfied. Pre- and post-Kronecker product P
-
1 to (28), we can derive that  P
-
1
⊗
Ω
^
⊗
P
-
1
<
0 and apply the change of variables X
=
P
-
1 and S
=
P
-
1
Q
P
-
1 such that
(29)
Ω
-
=
[
A
-
X
+
X
A
-
T
+
C
-
Y
+
Y
T
C
-
+
S
B
-
X
C
-
R
X
B
-
T
-
S
0
R
T
C
-
0
-
2
D
-
X
]
<
0
,
where Y
=
K
X and R
=
L
X.


## 4. Simulation

Consider the following distributed parameter system:(30)
∂
T
k
∂
t
=
∂
2
T
k
∂
z
2
+
α
T
k
(
t
)
+
β
T
k
(
t
-
τ
)
)
+
∑
j
=
1
N
‍
G
k
j
∂
2
T
j
∂
z
2
,
z
∈
(
0
,
π
)
,
∂
T
∂
z
|
z
=
0
=
0
,
∂
T
∂
z
|
z
=
π
=
0
,
T
(
t
,
z
)
=
cos
⁡
z
+
cos
⁡
2
z
,
t
∈
[
-
τ
,
0
]
,
where time delay τ
=
1, i
=
1,2
,
3 are the state variable of the ith neural network. Choosing the coupling matrix and coupling matrices as follows:
(31)
G
=
0.01
×
[
-
2
1
1
1
-
2
1
1
1
-
2
]
.
From Figures 1, 2, and 3, we can derive that T
k
(
k
=
1,2
,
3
) cannot synchronize. Next the controller will be designed to synchronize three distributed parameter systems.Figure 1
State of the first distributed parameter system.Figure 2
State of the second distributed parameter system.Figure 3
State of the third distributed parameter system.Letα
=
0.01, β
=
-
0.02. Figures 1–3 give the open-loop state T
k
(
k
=
1,2
,
3
) responses. As shown in the simulation, we can find that the coupled distributed parameter system can not synchronize partially under u
=
0. From Theorem 3, the observer and controller gains are constructed via solving the matrix inequalities (23), we obtain the controller gain matrices as follows:
(32)
K
=
[
-
1.9706
-
0.0201
-
0.0366
-
0.0202
-
1.9705
-
0.0310
-
0.0033
-
0.0088
-
1.9701
]
and L is zero matrix.Under the designed controller, Figures4, 5, and 6 show the coupled distributed parameter system synchronization.Figure 4
State of the first distributed parameter system under the controller.Figure 5
State of the second distributed parameter system under the controller.Figure 6
State of the third distributed parameter system under the controller.


## 5. Conclusion

In this paper, we have addressed the problem of synchronization problem for a class of coupled DPS. First, some sufficient conditions are given to ensure the synchronization of the coupled DPS. Secondly, in the case that the coupled DPS cannot synchronize by itself, the P-sD state feedback control has been designed and applied to force the coupled DPS to synchronize. The controller has been developed in terms of LMIs based directly on the error system. The coupled DPS can synchronize to the isolated node under the controller. Finally, the developed design method is applied to the simulation example, and the achieved simulation results show the effectiveness and benefit of the proposed controller. It is worth pointing out that the proposed design method is easily developed for the coupled DPS with Lipschtiz nonlinear function.



---

*Source: 418258-2014-03-16.xml*