# Almost Periodic Dynamics for Memristor-Based Shunting Inhibitory Cellular Neural Networks with Leakage Delays

Authors: Lin Lu; Chaoling Li
Journal: Computational Intelligence and Neuroscience
 (2016)

---

## Abstract

We investigate a class of memristor-based shunting inhibitory cellular neural networks with leakage delays. By applying a new Lyapunov function method, we prove that the neural network which has a unique almost periodic solution is globally exponentially stable. Moreover, the theoretical findings of this paper on the almost periodic solution are applied to prove the existence and stability of periodic solution for memristor-based shunting inhibitory cellular neural networks with leakage delays and periodic coefficients. An example is given to illustrate the effectiveness of the theoretical results. The results obtained in this paper are completely new and complement the previously known studies of Wu (2011) and Chen and Cao (2002).


---

## Body

## 1. Introduction

It is common knowledge that shunting inhibitory cellular neural networks (SICNNs) have a wide application in many fields such as image processing, signal processing, pattern recognition, psychophysics, speech, perception, robotics, and vision [1–3]. Thus, the theoretical analysis and applied research on SICNNs have attracted worldwide attention. During the past decades, memristor which is a new circuit element has received much attention due to its wide range of applications in computer, physics, electronic engineering, and so on [4, 5]. In particular, memristor has memory function and nanometer dimensions. The former can help us to deal with nanocomputing and the latter can provide a very high density and is less power hungry. The memristor can exhibit features as what the neurons in the human brain possess [4].In practical implementation, the time delays often occur in neural networks due to the finite switching speed of the neuron amplifiers and the finite signal transmission velocity. Here, we would like to point out that a typical time delay called leakage (or forgetting) delay may occur in the negative feedback term of the neural networks and plays an important role in characterizing the dynamical behavior of neural networks [6–11]. For example, time delay in the stabilizing negative feedback term may destabilize a system [12]. Balasubramaniam et al. [13] argued that the existence and uniqueness of the equilibrium point have nothing to do with time delays and initial conditions. Thus, it is important to study the leakage delays’ effect on the dynamical behavior of memristor-based neural networks. In recent years, there is some work on this topic. We refer the readers to [14–16].As is known to us, periodic oscillation of neural networks plays an important role in the daily life of human beings. Periodic oscillation of neural networks has been widely applied in many biological and cognitive activities. For example, periodic oscillatory or chaotic phenomena often occur in the human brain. Thus, some authors investigate the periodic oscillatory dynamical behavior of neural networks for grasping the mechanism of the human brain. We refer the readers to [17–19]. However, in many cases, the periodic parameters of neural networks may experience certain perturbations and then they may be not periodic. Thus, it is more reasonable to characterize the reality of neural networks with almost periodic parameters. In recent years, many authors consider the almost periodic oscillation of neural networks with or without delay and numerous good results have been available. For example, Liu et al. [20] focused on the almost periodic solution of impulsive Hopfield neural networks with finite distributed delays by applying fixed point theorems, Lyapunov functional, and some inequality techniques. Li et al. [21] investigated the existence and global exponential stability of almost periodic solution for high-order BAM neural networks with delays on time scales. By using a fixed point theorem and by constructing a suitable Lyapunov functional, authors established some sufficient conditions to ensure the existence and global exponential stability of almost periodic solution for high-order bidirectional associative memory neural networks with delays on time scales. Huang [22] presented some sufficient conditions for the existence and exponential stability of almost periodic solutions for fuzzy cellular neural networks with time-varying delays. Li et al. [23] established some sufficient conditions to ensure the existence and stability of pseudo almost periodic solution for neutral type high-order Hopfield neural networks with delays in leakage terms on time scales by means of fixed point theorem and the theory of calculus on time scales. For more results on this aspect, we refer the readers to [24–34]. To the best of our knowledge, there are no results on the existence and stability of almost periodic solution of memristor-based shunting inhibitory cellular neural networks with leakage delays.Inspired by the discussions above, in this article, we considered the following memristor-based shunting inhibitory cellular neural networks with leakage delays:(1)
x
˙
i
j
t
=
-
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
+
L
i
j
t
,where i
=
1,2
,
…
,
n,  j
=
1,2
,
…
,
m,  C
i
j represents the cell at the (
i
,
j
) position of the lattice, the r-neighborhood N
r
(
i
,
j
) of C
i
j is(2)
N
r
i
,
j
=
C
i
j
k
l
:
max
⁡
k
-
l
,
l
-
j
≤
ι
,
1
≤
k
≤
m
,
1
≤
l
≤
n
,
x
i
j is the activity of the cell C
i
j
,
L
i
j
(
t
) is the external input to C
i
j, the constant a
i
j
>
0 represents the passive decay rate of the cell activity, the activation function f
(
x
k
l
) is a positive continuous function representing the output or firing rate of the cell C
i
j, σ
i
j
(
t
)
≥
0 and τ
i
j
(
t
)
≥
0 denote the leakage delay and transmission delay at time t, C
i
j
k
l
t
,
x
i
j
t
≥
0 is memristive synaptic weights (which means the connection or coupling strength of postsynaptic activity of the cell transmitted to the cell C
i
j), which is defined as follows:(3)
C
i
j
k
l
t
,
ξ
=
C
i
j
k
l
∗
t
,
h
i
j
ξ
>
T
i
j
,
C
i
j
k
l
∗
∗
t
,
h
i
j
ξ
<
T
i
j
,for t
∈
R, a
i
j
(
t
,
x
)
=
C
i
j
k
l
∗
(
t
) or C
i
j
k
l
∗
∗
(
t
) when h
i
j
(
ξ
)
=
T
i
j, when h
i
j
:
R
n
+
m
→
R
(
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
) are threshold level functions, T
i
j
∈
R
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m are threshold level, and C
i
j
k
l
∗ and C
i
j
k
l
∗
∗ are all continuous functions. The initial states associated with (3) are given by(4)
x
i
j
s
=
φ
i
j
s
,
s
∈
-
δ
,
0
,
φ
i
j
s
∈
C
-
δ
,
0
,
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
, where δ
=
m
a
x
t
∈
R
{
τ
i
j
(
t
)
,
η
i
j
(
t
)
}
.The main purpose of this article is to investigate the existence and exponential stability of the almost periodic solutions for system (1). With the aid of new Lyapunov function techniques, we establish some new sufficient criteria which guarantee the existence, uniqueness, and exponential stability of the almost periodic solution of system (1). Also, the derived results on the almost periodic solution are applied to prove the existence and stability of periodic solution for memristor-based shunting inhibitory cellular neural networks with leakage delays and periodic coefficients. The obtained results of this article are new and complement previously known publications.The remainder of the paper is organized as follows. In Section2, some necessary definitions and lemmas are stated. In Section 3, a set of sufficient criteria which guarantee the global existence and boundedness of any solutions and the existence and exponential stability of an almost periodic solution of neural networks (1) are established. The global exponential periodicity and stability of system (1) are analyzed in Section 4. An example is given to show the correctness of the theoretical predictions in Section 5. A brief conclusion is drawn in Section 6.


## 2. Preliminaries

In this section, we list several definitions and notations. SupposeE
∈
R
n
+
m; then x
→
F
(
x
) is called a set-valued map from E to R
n
+
m, if, for each point x
∈
E, there exists a nonempty set F
(
x
)
⊂
R
n
+
m
. A set-valued map F with nonempty values is said to be upper semicontinuous at x
0
∈
E, if, for any open set N containing F
(
x
0
), there exists a neighborhood M of x
0 such that F
(
M
)
⊆
N
. The map F
(
x
) is said to have a closed (convex, compact) image if, for each x
∈
E
,
F
(
x
) is closed (convex, compact). For φ
∈
C
(
-
δ
,
0
], let φ
=
s
u
p
s
∈
-
δ
,
0
φ
s
. Given the function V
:
R
n
+
m
→
R,  ∇
V denotes the gradient of V and ∂
V denotes Clarke’s generalized gradient of V
.In (1), since C
i
j
k
l
(
t
,
ξ
) is discontinuous, the classical definition of the solution for differential equations cannot apply here. To handle this problem, Filippov developed a solution concept for the differential equation with a discontinuous right-hand side. Based on this definition, a differential equation with a discontinuous right-hand side has the same solution set as a certain differential inclusion. In what follows, we use this definition to discuss dynamical behavior of (1). Let the set-valued maps be as follows:(5)
c
o
¯
C
i
j
k
l
t
,
ξ
=
C
i
j
k
l
∗
t
,
h
i
j
ξ
>
T
i
j
,
c
o
¯
C
i
j
k
l
∗
t
,
C
i
j
k
l
∗
∗
t
,
h
i
j
ξ
=
T
i
j
,
C
i
j
k
l
∗
∗
t
,
h
i
j
ξ
<
T
i
j
,for t
∈
R
,
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, where c
o
¯ denotes the convex closure of a set. Obviously, c
o
¯
(
a
i
j
(
t
,
ξ
)
) is all closed, convex, and compact in ξ for each t
∈
R
,
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
. We define the Filippov solution of (1) as the following.Definition 1.
A functionx
(
t
) is said to be a solution of (1) on [
0
,
T
) with initial condition (4), if x
(
t
) is absolutely continuous on any compact interval of [
0
,
T
) and satisfies differential inclusions(6)
d
x
i
j
t
d
t
∈
-
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
+
L
i
j
t
, or, equivalently, there exist D
i
j
k
l
(
t
)
∈
c
o
¯
(
C
i
j
k
l
(
t
,
x
i
j
(
t
)
)
) satisfying(7)
d
x
i
j
t
d
t
=
-
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
∑
D
k
l
∈
N
ι
i
,
j
D
i
j
k
l
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
+
L
i
j
t
,for a.e. t
∈
[
0
,
T
], i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
.Definition 2.
A continuous functionx
(
t
)
:
R
→
R
n
+
m is said to be almost periodic on R if, for any ε
>
0, it is possible to find a real number l
=
l
(
ε
)
>
0
. For any interval with length l, there exists a number ω
=
ω
(
ε
) in this interval, such that x
t
+
ω
-
x
t
<
ε for all t
∈
R
.Definition 3.
A continuous functionx
(
t
)
:
R
→
R
n
+
m is said to be asymptotically almost periodic on R if, for any ε
>
0, there exist T
>
0,  l
=
l
(
ε
), and ω
=
ω
(
ε
) in any interval with the length of l, such that x
t
+
ω
-
x
t
<
ε for all t
≥
T
.Definition 4.
The neural networks model is said to be globally exponentially almost periodic if the statex
(
t
,
φ
) of the neural networks model is globally exponentially convergent to an almost periodic state x
∗
t
,
φ; that is, there are constants M
≥
1 and μ
>
0 such that, for any φ
∈
C
(
-
δ
,
0
], x
(
t
,
φ
)
-
x
∗
(
t
,
φ
)
≤
M
φ
-
ψ
e
-
μ
t
,
t
≥
0
. In addition, if x
∗
(
t
,
φ
) is a periodic solution (equilibrium), then the neural networks model is said to be globally exponentially periodic (stable).Definition 5 (see [35]).
V
(
x
)
:
R
n
+
m
→
R is said to be regular, if, for each x
∈
R
n
+
m and ν
∈
R
n
+
m,(i)
there exists the usual right or left directional derivative(8)
D
+
V
x
,
ν
=
l
i
m
h
→
0
+
V
x
+
h
ν
-
V
x
h
;
(ii)
the generalized directional derivative ofV at x in the direction ν
∈
R
n
+
m is defined as (9)
D
+
+
V
x
,
ν
=
l
i
m
y
→
x
,
h
→
0
+
V
y
+
h
ν
-
V
y
h
;then D
+
V
(
x
,
ν
)
=
D
+
+
V
(
x
,
ν
)
.Definition 6.
For a locally Lipschitz functionV
:
R
n
+
m
→
R, one can define Clarke’s generalized gradient of V at point x, as follows: (10)
∂
V
x
=
c
o
¯
lim
k
→
∞
⁡
∇
V
x
k
:
x
k
⟶
x
,
x
k
∉
N
,
x
k
∉
Ω
,where Ω
⊂
R
n
+
m is the set of points where V is not differentiable and N
⊂
R
n
+
m is an arbitrary set with measure zero.Lemma 7 (see [36]).
IfV
(
x
)
:
R
n
+
m
→
R is Clarke’s regular and x
(
t
)
:
[
0
,
+
∞
)
→
R
n
+
m is absolutely continuous on any compact interval of [
0
,
+
∞
), then x
(
t
) and V
(
x
(
t
)
)
:
[
0
,
+
∞
)
∈
R
n
+
m are differential for a.a. t
∈
[
0
,
+
∞
), and one will have d
v
(
t
)
/
d
t
=
γ
(
t
)
T
x
˙
(
t
)
,
∀
γ
(
t
)
∈
∂
V
x
t, where ∂
V
(
x
(
t
)
) is Clarke’s generalized gradient.Lemma 8 (see [37]).
Let matrixM
=
(
m
i
j
)
n
×
n have nonpositive off-diagonal elements. Then, M is a nonsingular M-matrix if and only if one of the following conditions holds:(1)
There existn positive constants α
1
,
α
2
,
…
,
α
n such that(11)
α
i
m
i
i
-
∑
j
=
1
,
j
≠
i
n
α
j
m
j
i
>
0
,
i
=
1,2
,
…
,
n
.
(2)
There existn positive constants β
1
,
β
2
,
…
,
β
n such that (12)
β
i
m
i
i
-
∑
j
=
1
,
j
≠
i
n
β
j
m
i
j
>
0
,
i
=
1,2
,
…
,
n
.Denotea
+
=
s
u
p
t
∈
R
a
t
,
a
-
=
i
n
f
t
∈
R
a
t, where a
:
R
→
R is a bounded continuous function.Throughout this paper, we assume that the following conditions are satisfied:(H1)
Fori
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, C
i
j
k
l
∗
(
t
)
,
C
i
j
k
l
∗
∗
(
t
)
,
L
i
j
(
t
)
,
a
i
j
(
t
)
,
σ
i
j
(
t
), and τ
i
j
(
t
) are continuous functions and are almost periodic; that is, for any ε
>
0, it is possible to find a real number l
=
l
ε
>
0; for any interval with length l, there exists a number ω
=
ω
(
ε
) in this interval, such that C
i
j
k
l
∗
t
+
ω
-
C
i
j
k
l
∗
t
<
ε,  C
i
j
k
l
∗
∗
t
+
ω
-
C
i
j
k
l
∗
∗
t
<
ε,  L
i
j
t
+
ω
-
L
i
j
t
<
ε,  a
i
j
t
+
ω
-
a
i
j
t
<
ε,  τ
i
j
t
+
ω
-
τ
i
j
t
<
ε, and  σ
i
j
t
+
ω
-
σ
i
j
t
<
ε for all t
∈
R
.
(H2)
Fori
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, σ
i
j
(
t
) is bounded above and below by positive constants and σ
i
j
′
(
t
) is a bounded continuous function and σ
i
j
+
a
i
j
+
<
1
.
(H3)
There exists constantL
>
0 such that f
u
-
f
v
≤
L
u
-
v for u
,
v
∈
R
.
(H4)
There exist positive constantsδ
11
,
δ
12
,
…
,
δ
1
m
,
δ
21
,
δ
22
,
…
,
δ
n
m and ς such that (13)
-
a
i
j
t
1
-
2
a
i
j
+
σ
i
j
+
-
a
i
j
t
-
1
-
σ
i
j
′
t
a
i
j
t
-
σ
i
j
t
δ
k
l
1
-
a
k
l
+
σ
k
l
+
+
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
L
δ
i
j
ϱ
/
ς
1
-
a
i
j
+
σ
i
j
+
2
<
-
ς
,
for allt
>
0 and i
,
j
=
1,2
,
…
,
n, j
=
1,2
,
…
,
m, where (
C
i
j
k
l
)
u
(
t
)
=
max
⁡
C
i
j
k
l
∗
t
,
C
i
j
k
l
∗
∗
t and ϱ
=
m
a
x
1
≤
i
≤
n
,
1
≤
j
≤
m
δ
i
j
(
C
i
j
k
l
)
+
u
f
0
+
L
i
j
+
.
(H5)
There exists a nonempty subsetΛ
i
j
⊂
R
(
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
) satisfying the following property: if h
i
j
(
x
)
≤
T
i
j
≤
h
i
j
(
y
), then there exists λ
i
j
∈
Λ
i
j such that x
i
j
≤
λ
i
j
≤
y
i
j or y
i
j
≤
λ
i
j
≤
x
i
j
.
(H6)
Fori
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, f
(
λ
i
j
)
=
0 for any λ
i
j
∈
Λ
i
j
.


## 3. Boundedness and Almost Periodicity

In this section, we will prove the existence of bounded solution and the exponential stability of almost periodic solution for (1).Theorem 9.
Assume that assumptions (H1)–(H4) hold. Letx
(
t
) be the solution of (1) with initial condition(14)
x
i
j
t
=
φ
i
j
t
,
φ
i
j
t
-
∫
t
-
σ
i
j
t
t
a
i
j
s
φ
i
j
s
d
s
<
δ
i
j
ϱ
ς
,
t
∈
-
δ
,
0
,where i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, and (15)
ϱ
=
m
a
x
1
≤
i
≤
n
,
1
≤
j
≤
m
δ
i
j
C
i
j
k
l
+
u
f
0
+
L
i
j
+
.Then,(16)
x
i
j
t
-
∫
t
-
σ
i
j
t
t
a
i
j
s
x
i
j
s
d
s
<
δ
i
j
ϱ
ς
,
(17)
x
i
j
t
≤
δ
i
j
ϱ
/
ς
1
-
a
i
j
+
σ
i
j
+
,where t is in the interval of existence and i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
.Proof.
Letx
(
t
) be a solution of (1) with initial condition (14). For i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, there exists D
i
j
k
l
(
t
)
∈
c
o
¯
C
i
j
k
l
t
,
x
i
j
t satisfying(18)
x
˙
i
j
t
=
-
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
∑
D
k
l
∈
N
ι
i
,
j
D
i
j
k
l
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
+
L
i
j
t
,for a.a. t
∈
[
0
,
T
]
. By (5) and (H3), we have(19)
D
i
j
k
l
t
≤
m
a
x
C
i
j
k
l
∗
t
,
C
i
j
k
l
∗
∗
t
≤
C
i
j
k
l
u
,
f
x
i
j
t
≤
f
x
i
j
t
-
f
0
+
f
0
≤
L
x
i
j
t
+
f
0
.For t in the interval of existence and i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, denote(20)
u
i
j
t
=
x
i
j
t
-
∫
t
-
σ
i
j
t
t
a
i
j
s
x
i
j
s
d
s
.Suppose (16) holds; then, for a given t
∗
>
0 in the interval of existence and i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, we get(21)
x
i
j
t
∗
≤
u
i
j
t
∗
+
∫
t
∗
-
σ
i
j
t
∗
t
∗
a
i
j
s
x
i
j
s
d
s
≤
δ
i
j
ϱ
ς
+
a
i
j
+
σ
i
j
+
sup
s
∈
-
δ
,
t
∗
⁡
x
i
j
s for any t
∗
∈
(
-
δ
,
t
]. Then,(22)
s
u
p
s
∈
-
δ
,
t
∗
x
i
j
s
<
δ
i
j
ϱ
ς
+
a
i
j
+
σ
i
j
+
s
u
p
s
∈
-
δ
,
t
∗
x
i
j
s for any t
∗
∈
(
-
δ
,
t
]
. It follows from (22) that(23)
x
i
j
t
≤
δ
i
j
ϱ
/
ς
1
-
a
i
j
+
σ
i
j
+
.Therefore, (17) holds. Thus, it suffices to prove (16). Assume that (16) does not hold. Then, there exist k
∈
{
1,2
,
…
,
n
},  l
∈
{
1,2
,
…
,
m
}, and t
0
>
0 such that u
k
l
t
0
=
δ
k
l
ϱ
/
ς and (16) holds for all t
∈
[
-
δ
,
t
0
] and i
=
1,2
,
…
,
n,  j
=
1,2
,
…
,
m and hence D
-
u
k
l
t
0
≥
0. From system (1), we have(24)
d
u
k
l
t
d
t
=
x
k
l
′
t
-
a
k
l
t
x
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
x
k
l
t
-
σ
k
l
t
=
-
a
k
l
t
x
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
x
k
l
t
-
σ
k
l
t
-
a
k
l
t
x
k
l
t
-
σ
k
l
t
-
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
f
x
k
l
t
-
τ
k
l
t
x
k
l
t
+
L
k
l
t
=
-
a
k
l
t
x
k
l
t
-
a
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
x
k
l
t
-
σ
k
l
t
-
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
f
x
k
l
t
-
τ
k
l
t
x
k
l
t
+
L
k
l
t
=
-
a
k
l
t
u
k
l
t
-
a
k
l
t
∫
t
-
σ
k
l
t
t
a
k
l
s
x
k
l
s
d
s
-
a
k
l
t
-
1
-
η
k
l
′
t
a
k
l
t
-
σ
k
l
t
x
k
l
t
-
σ
k
l
t
-
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
f
x
k
l
t
-
τ
k
l
t
x
k
l
t
+
L
k
l
t
. It follows from (H3) that(25)
D
-
u
k
l
t
0
≤
-
a
k
l
t
0
u
k
l
t
0
+
a
k
l
t
0
∫
t
-
σ
k
l
t
0
t
0
a
k
l
s
x
k
l
s
d
s
+
a
k
l
t
0
-
1
-
σ
k
l
′
t
0
a
k
l
t
0
-
σ
k
l
t
0
x
k
l
t
0
-
σ
k
l
t
0
+
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
0
f
x
k
l
t
0
-
τ
k
l
t
0
x
k
l
t
0
+
L
k
l
t
0
≤
-
a
k
l
t
0
u
k
l
t
0
+
a
k
l
t
0
∫
t
-
σ
k
l
t
0
t
0
a
k
l
s
x
k
l
s
d
s
+
a
k
l
t
0
-
1
-
σ
k
l
′
t
0
a
k
l
t
0
-
σ
k
l
t
0
x
k
l
t
0
-
σ
k
l
t
0
+
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
0
L
x
k
l
t
0
-
τ
k
l
t
0
+
f
0
x
k
l
t
0
+
L
k
l
t
0
≤
-
a
k
l
t
0
δ
k
l
ϱ
ς
+
a
k
l
t
0
a
k
l
+
η
k
l
+
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
+
a
k
l
t
0
-
1
-
σ
k
l
′
t
0
a
k
l
t
0
-
σ
k
l
t
0
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
+
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
t
0
L
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
+
f
0
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
+
L
k
l
+
≤
-
a
k
l
t
0
1
-
2
a
k
l
+
σ
k
l
+
-
a
k
l
t
0
-
1
-
σ
k
l
′
t
0
a
k
l
t
0
-
σ
k
l
t
0
δ
k
l
1
-
a
k
l
+
σ
k
l
+
+
∑
D
k
l
∈
N
r
i
,
j
D
i
j
k
l
u
t
0
L
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
2
ϱ
ς
+
ϱ
<
0
,which is a contradiction and shows that (16) holds. The proof of Theorem 9 is complete.Lemma 10.
Assume that (H3), (H5), and (H6) hold; then, for anyx
=
(
x
11
,
x
12
,
…
,
x
n
m
)
T,  y
=
(
y
11
,
y
12
,
…
,
y
n
m
)
T
∈
R
n
m, one has(26)
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
-
c
o
¯
C
i
j
k
l
t
,
y
i
j
t
f
y
i
j
≤
C
i
j
k
l
u
t
L
x
i
j
-
y
i
j
,where i
=
1,2
,
…
,
n,  j
=
1,2
,
…
,
m
.Proof.
For any giveni
∈
{
1,2
,
…
,
n
},  j
∈
{
1,2
,
…
,
m
}, and x
i
j
,
y
i
j
∈
R
n
+
m, we consider three cases.
Ifh
i
j
(
x
)
,
h
i
j
(
y
)
<
T
i
j, then (27)
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
-
c
o
¯
C
i
j
k
l
t
,
y
i
j
t
f
y
i
j
=
C
i
j
k
l
∗
t
f
x
i
j
-
C
i
j
k
l
∗
t
f
y
i
j
≤
C
i
j
k
l
u
t
L
x
i
j
-
y
i
j
.If h
i
j
(
x
)
,
h
i
j
(
y
)
>
T
i
j, then (28)
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
-
c
o
¯
C
i
j
k
l
t
,
y
i
j
t
f
y
i
j
=
C
i
j
k
l
∗
∗
t
f
x
i
j
-
C
i
j
k
l
∗
∗
t
f
y
i
j
≤
C
i
j
k
l
u
t
L
x
i
j
-
y
i
j
.If h
i
j
(
x
)
≤
T
i
j
≤
h
i
j
(
y
) or h
i
j
(
y
)
≤
T
i
j
≤
h
i
j
(
x
), then it follows from (H4) that there exists λ
j
∈
Λ
j such that x
i
j
≤
λ
i
j
≤
y
i
j or y
i
j
≤
λ
i
j
≤
x
i
j
. Let x
i
j
≤
λ
i
j
≤
y
i
j. In this case, from (H5), we get (29)
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
-
c
o
¯
C
i
j
k
l
t
,
y
i
j
t
f
y
i
j
=
C
i
j
k
l
∗
∗
t
f
x
i
j
-
C
i
j
k
l
∗
t
f
y
i
j
≤
C
i
j
k
l
∗
∗
t
f
x
i
j
-
f
λ
j
+
C
i
j
k
l
∗
t
f
λ
j
-
f
y
i
j
≤
C
i
j
k
l
u
t
L
x
i
j
-
y
i
j
.Based on all the cases above, we can conclude that (26) holds. The proof of Lemma 10 is complete.Now, we state our main result.Theorem 11.
If (H1)–(H6) hold, then there exists a unique almost periodic solutionx
∗
(
t
,
ψ
) for system (1) which is globally exponentially stable; that is, for any other solution x
(
t
,
ψ
) of system (1), there exist constants M
,
μ
>
0 such that x
t
,
φ
-
x
∗
t
,
ψ
≤
M
φ
-
ψ
e
-
μ
t for all t
>
0.Proof.
First, we prove that any solution of (1) is asymptotically almost periodic; that is, for any ε
>
0, there exist T
>
0,  l
=
l
(
ε
), and ω
=
ω
(
ε
) in any interval with the length of l, such that x
t
+
ω
-
x
t
≤
ε for all t
≥
T
.
For anyε
>
0, let ω
=
ω
(
ε
) and y
i
j
(
t
)
=
x
i
j
(
t
+
ω
)
-
x
i
j
(
t
),  i
=
1,2
,
…
,
n,  j
=
1,2
,
…
,
m, and then we get(30)
d
y
i
j
t
d
t
∈
-
a
i
j
t
x
i
j
t
+
ω
-
σ
i
j
t
+
ω
+
L
i
j
t
+
ω
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
+
ω
,
x
i
j
t
+
ω
f
x
i
j
t
+
ω
-
τ
i
j
t
+
ω
x
i
j
t
+
ω
+
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
L
i
j
t
+
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
=
-
a
i
j
t
y
i
j
t
+
A
i
j
t
,
x
t
,
y
t
+
Θ
i
j
t
,
ω
,where A
i
j
(
t
,
x
(
t
)
,
y
(
t
)
) and Θ
i
j
(
t
,
ω
) are defined as follows:(31)
A
i
j
t
,
x
t
,
y
t
=
c
o
¯
C
i
j
k
l
t
,
u
+
v
f
u
+
v
-
c
o
¯
C
i
j
k
l
t
,
u
f
u
,for all u
,
v
∈
R
n
+
m and(32)
Θ
i
j
t
,
ω
=
-
a
i
j
t
+
ω
-
a
i
j
t
+
ω
-
σ
i
j
t
+
L
i
j
t
+
ω
-
L
i
j
t
+
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
+
ω
,
x
i
j
t
+
ω
f
x
i
j
t
+
ω
-
τ
i
j
t
+
ω
x
i
j
t
+
ω
.In view of (H1) and the boundedness of x
(
t
), we can conclude that, for any ε
>
0, there exist l
=
l
(
ε
)
>
0, and ω
=
ω
(
ε
) in any interval with the length of l, such that, for any Θ
∗
(
t
,
ω
)
∈
c
o
¯
(
Θ
(
t
,
ω
)
), Θ
∗
t
,
ω
<
N
ε
/
2 for all t
≥
0, where N
>
0 is a constant. Let(33)
d
y
i
j
t
d
t
=
-
a
i
j
t
y
i
j
t
+
A
i
j
∗
t
,
x
t
,
y
t
+
Θ
i
j
∗
t
,
ω
,where A
i
j
∗
(
t
,
x
(
t
)
,
y
(
t
)
)
∈
c
o
¯
(
A
i
j
(
t
,
x
(
t
)
,
y
(
t
)
)
) and Θ
∗
(
t
,
ω
)
∈
c
o
¯
(
Θ
(
t
,
ω
)
)
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
.
In view of (H3), we can chooser
>
0 and ς
>
0 such that(34)
-
a
k
l
t
-
r
1
-
2
a
k
l
+
σ
k
l
+
-
a
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
δ
k
l
1
-
a
k
l
+
σ
k
l
+
+
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
e
r
τ
i
j
t
L
δ
k
l
ϱ
/
ς
1
-
a
k
l
+
σ
k
l
+
2
<
-
ς
. Let T
0
≥
m
a
x
{
0
,
ω
}
. For i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m and t
∈
R, denote(35)
U
i
j
t
=
e
r
t
y
i
j
t
-
∫
t
-
σ
i
j
t
t
a
i
j
s
e
r
s
y
i
j
s
d
s
.Then, for t
≥
T
0 and i
=
1,2
,
…
,
n,  j
=
1,2
,
…
,
m, we have(36)
d
U
i
j
t
d
t
=
r
e
r
t
y
i
j
t
+
e
r
t
y
i
j
′
t
-
a
i
j
t
e
r
t
y
i
j
t
+
1
-
σ
i
j
′
t
a
i
j
t
-
σ
i
j
t
e
r
t
-
σ
i
j
t
y
i
j
t
-
σ
i
j
t
=
-
a
i
j
t
-
r
e
r
t
y
i
j
t
+
e
r
t
-
a
i
j
t
-
1
-
σ
i
j
′
t
a
i
j
t
-
σ
i
j
t
e
-
r
σ
i
j
t
×
y
i
j
t
-
σ
i
j
t
+
A
i
j
∗
t
,
x
t
,
y
t
+
Θ
i
j
∗
t
,
ω
=
-
a
i
j
t
-
r
U
i
j
t
-
a
i
j
t
-
r
∫
t
-
σ
i
j
t
t
a
i
j
s
e
r
s
y
i
j
s
d
s
+
e
r
t
-
a
i
j
t
-
1
-
σ
i
j
′
t
a
i
j
t
-
σ
i
j
t
e
-
r
σ
i
j
t
×
y
i
j
t
-
σ
i
j
t
+
A
i
j
∗
t
,
x
t
,
y
t
+
Θ
i
j
∗
t
,
ω
.Now, we define a candidate Lyapunov function as follows:(37)
W
1
t
=
m
a
x
U
i
j
t
δ
i
j
,
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
,
W
1
∗
t
=
s
u
p
s
≤
t
W
1
s
. Obviously, W
1
∗
(
t
) is nondecreasing. It follows that(38)
e
r
ρ
y
i
j
ρ
δ
i
j
-
1
≤
e
r
ρ
y
i
j
ρ
-
∫
ρ
-
σ
i
j
ρ
ρ
a
i
j
s
e
r
s
y
i
j
s
d
s
δ
i
j
-
1
+
∫
ρ
-
σ
i
j
ρ
ρ
a
i
j
s
e
r
s
y
i
j
s
d
s
δ
i
j
-
1
≤
W
1
∗
t
+
a
i
j
+
σ
i
j
+
s
u
p
s
≤
t
e
r
s
y
i
j
s
δ
i
j
-
1 for all ρ
≤
t, where i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
. By (H3), we have(39)
e
r
ρ
y
i
j
ρ
δ
i
j
-
1
≤
W
1
∗
t
1
-
a
i
j
+
σ
i
j
+for all t
≥
0
,
ρ
≤
t, where i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
. For any given t
≥
0, there exist k
∈
{
1,2
,
…
,
n
},  l
∈
{
1,2
,
…
,
m
} such that(40)
W
1
t
=
U
k
l
t
δ
k
l
.Calculating the derivative d
W
1
(
t
)
/
d
t along the positive half trajectory of (1) yields(41)
d
W
1
t
d
t
≤
-
a
k
l
t
-
r
W
1
t
+
1
δ
k
l
-
a
k
l
t
-
r
∫
t
-
σ
k
l
t
t
a
k
l
s
e
r
s
y
k
l
s
d
s
+
e
r
t
-
a
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
e
-
r
σ
k
l
t
×
y
k
l
t
-
σ
k
l
t
+
A
k
l
∗
t
,
x
t
,
y
t
+
Θ
k
l
∗
t
,
ω
≤
-
a
k
l
t
-
r
W
1
t
+
a
k
l
t
-
r
σ
k
l
+
a
k
l
+
W
1
∗
t
1
-
σ
k
l
+
a
k
l
+
+
a
k
l
t
e
r
σ
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
W
1
∗
t
1
-
σ
k
l
+
a
k
l
+
+
1
δ
k
l
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
L
e
r
τ
k
l
t
W
1
∗
t
1
-
a
k
l
+
σ
k
l
+
2
+
e
r
t
Θ
k
l
∗
t
,
ω
≤
e
r
t
δ
k
l
Θ
k
l
∗
t
,
ω
-
a
k
l
t
-
r
1
-
2
σ
k
l
+
a
k
l
+
-
a
k
l
e
r
σ
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
δ
k
l
1
-
σ
k
l
+
a
k
l
+
-
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
L
e
r
τ
k
l
t
δ
k
l
1
-
a
k
l
+
σ
k
l
+
2
W
1
t
δ
k
l
≤
e
r
t
δ
k
l
Θ
k
l
∗
t
,
ω
≤
N
ε
2
δ
m
i
n
e
r
twhen W
1
∗
(
t
)
≤
W
1
(
t
)
. Thus,(42)
d
W
1
∗
t
d
t
≤
N
ε
2
δ
m
i
n
e
r
t
,
∀
t
∈
0
,
T
. Then, for all i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, we have(43)
y
i
j
t
≤
δ
i
j
W
1
∗
t
e
r
t
1
-
σ
i
j
+
a
i
j
+
≤
δ
i
j
W
1
∗
0
e
r
t
m
i
n
1
≤
i
≤
n
,
1
≤
j
≤
m
1
-
σ
i
j
+
a
i
j
+
+
δ
i
j
N
ε
2
δ
m
i
n
.Thus, there exists a constant T
>
0 such that, for any t
>
T,(44)
y
t
≤
N
ε
r
δ
m
i
n
∑
i
=
1
n
∑
j
=
1
m
δ
i
j
.Taking 0
<
N
<
r
δ
m
i
n
/
∑
i
=
1
n
‍
δ
i, we have y
(
t
)
<
ε for any t
>
T
. Namely, for any T
>
0, there exist T
>
0
,
l
=
l
(
ε
)
>
0, and ω
=
ω
(
ε
) in any interval with the length of l, such that x
(
t
+
ω
)
-
x
(
t
)
<
ε for all t
≥
T
. Therefore, any solution x
(
t
) of (1) with initial condition (14) is asymptotically almost periodic.
Next, we prove that there exists at least one almost periodic solution of (1).
Letx
(
t
) be any solution of (1) with initial conditions (4) and (14). It is easy to see that, for any sequence {
t
k
}
k
∈
N satisfying l
i
m
t
→
t
k
=
+
∞, the sequence {
x
(
t
+
t
k
)
}
k
∈
N is equicontinuous and uniformly bounded. In view of Arzela-Ascoli theorem and diagonal selection principle, we can select a subsequence of {
t
k
} (still denoted by {
t
k
}), such that x
(
t
+
t
k
) uniformly converges to a continuous function x
∗
(
t
) on any compact set of R
. We next prove that x
∗
(
t
) is a solution of (1).
Letz
i
j
(
t
,
t
k
)
=
x
i
j
(
t
+
t
k
)
-
x
i
j
∗
(
t
) and(45)
Θ
¯
i
j
t
,
t
k
=
Θ
i
j
t
,
t
k
-
a
i
j
t
z
i
j
t
,
t
k
+
c
o
¯
A
i
j
t
,
x
i
j
∗
t
,
z
i
j
t
,
t
k
;then(46)
c
o
¯
Θ
¯
i
j
t
,
t
k
≤
c
o
¯
Θ
i
j
t
,
t
k
+
a
i
j
+
z
i
j
t
,
t
k
+
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
+
u
z
i
j
t
,
t
k
.In view of the boundedness of x
(
t
), we can select the sequence {
t
k
} satisfying Θ
i
j
∗
t
,
t
k
≤
1
/
k for any Θ
i
j
∗
(
t
,
t
k
)
∈
c
o
¯
(
Θ
i
j
) and all t
≥
0
. From this and (48), we have l
i
m
k
→
+
∞
c
o
¯
(
Θ
¯
i
j
(
t
,
t
k
)
)
=
0 for all t
≥
0 and i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
. Applying Lebesgue’s dominated convergence theorem, we have(47)
x
∗
t
+
f
-
x
∗
t
=
l
i
m
k
→
∞
x
t
k
+
f
-
x
∗
t
k
=
l
i
m
k
→
∞
∫
t
t
+
f
x
˙
s
+
t
k
d
s
∈
l
i
m
k
→
∞
∫
t
t
+
f
-
a
i
j
s
+
t
k
x
i
j
s
+
t
k
-
σ
i
j
s
+
t
k
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
s
+
t
k
,
x
i
j
s
+
t
k
f
x
i
j
s
+
t
k
-
τ
i
j
s
+
t
k
×
x
i
j
s
+
t
k
+
L
i
j
s
+
t
k
d
s
∈
∫
t
t
+
f
l
i
m
k
→
∞
c
o
¯
Θ
i
j
∗
t
,
t
k
-
a
i
j
s
x
i
j
∗
t
-
σ
i
j
s
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
s
,
x
i
j
∗
s
f
x
i
j
∗
s
τ
i
j
s
+
L
i
j
s
d
s
=
∫
t
t
+
f
-
a
i
j
s
x
i
j
∗
t
-
σ
i
j
s
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
s
,
x
i
j
∗
s
f
x
i
j
∗
s
τ
i
j
s
+
L
i
j
s
d
sfor all t
∈
R and f
∈
R. Thus, x
∗
(
t
) is a solution of (1).
Here, we will prove thatx
∗
(
t
) is the almost periodic solution of (1). By the proof of the above step, for any ε
>
0, there exist T
>
0,  l
=
l
(
ε
), and ω
=
ω
(
ε
) in any interval with the length of l, such that x
t
+
ω
-
x
t
≤
ε for all t
≥
T
. Hence, there exists sufficiently large constant K
>
0 such that x
t
+
t
k
+
ω
-
x
t
+
t
k
≤
ε for all t
∈
R and k
>
K
. Let k
→
+
∞, and then it is easy to obtain that x
∗
t
+
ω
-
x
∗
t
≤
ε for all t
∈
R; that is, x
∗
(
t
) is the almost periodic solution of (1). Finally, we mainly prove that the almost periodic solution of (1) is unique and globally exponentially stable.
Letx
(
t
) be any solution of (1) with initial conditions (4) and (14), and let x
∗
(
t
) be an almost periodic solution of (1); that is,(48)
d
x
i
j
t
d
t
∈
-
a
i
j
t
x
i
j
t
-
σ
i
j
t
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
t
x
i
j
t
+
L
i
j
t
,
(49)
d
x
i
j
∗
t
d
t
∈
-
a
i
j
t
x
i
j
∗
t
-
σ
i
j
t
-
∑
C
k
l
∈
N
ι
i
,
j
c
o
¯
C
i
j
k
l
t
,
x
i
j
∗
t
f
x
i
j
∗
t
-
τ
i
j
t
x
i
j
∗
t
+
L
i
j
t
.Let w
i
j
(
t
)
=
x
i
j
(
t
)
-
x
i
j
∗
(
t
), and then(50)
d
u
i
j
t
d
t
∈
-
a
i
j
t
w
i
j
∗
t
-
σ
i
j
t
-
c
o
¯
A
i
j
t
,
x
i
j
∗
t
-
τ
i
j
t
,
w
i
j
t
-
τ
i
j
t
,where A
i
j
(
t
,
x
j
∗
(
t
)
,
u
j
(
t
)
) is defined by (31). Let(51)
d
u
i
j
t
d
t
=
-
a
i
j
t
w
i
j
∗
t
-
σ
i
j
t
-
c
o
¯
A
i
j
t
,
x
i
j
∗
t
-
τ
i
j
t
,
w
i
j
t
-
τ
i
j
t
,where A
i
j
t
,
x
i
j
t
-
τ
i
j
t
,
u
i
j
t
-
τ
i
j
t
∈
c
o
¯
A
i
j
t
,
x
i
j
t
-
τ
i
j
t
,
u
i
j
t
-
τ
i
j
t  i
,
j
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m.
Fori
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m and t
∈
R, denote(52)
W
i
j
t
=
e
r
t
w
i
j
t
-
∫
t
-
σ
i
j
t
t
a
i
j
s
e
r
s
w
i
j
s
d
s
.Then, for t
≥
T
0 and i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m, we have(53)
d
W
i
j
t
d
t
=
-
a
i
j
t
-
r
W
i
j
t
-
a
i
j
t
-
r
∫
t
-
σ
i
j
t
t
a
i
j
s
e
r
s
w
i
j
s
d
s
+
e
r
t
-
a
i
j
t
-
1
-
σ
i
j
′
t
a
i
j
t
-
σ
i
j
t
e
-
r
σ
i
j
t
×
w
i
j
t
-
σ
i
j
t
+
e
r
t
A
i
j
t
,
x
∗
t
,
w
t
.We define another candidate Lyapunov function as follows:(54)
V
2
t
=
max
⁡
W
i
j
t
δ
i
j
,
i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
,
V
2
∗
t
=
s
u
p
s
≤
t
V
2
s
.Obviously, V
2
∗
(
t
) is nondecreasing. It follows that(55)
e
r
ρ
w
i
j
ρ
δ
i
j
-
1
≤
e
r
ρ
w
i
j
ρ
-
∫
ρ
-
σ
i
j
ρ
ρ
a
i
j
s
e
r
s
w
i
j
s
d
s
δ
i
j
-
1
+
∫
ρ
-
σ
i
j
ρ
ρ
a
i
j
s
e
r
s
w
i
j
s
d
s
δ
i
j
-
1
≤
W
1
∗
t
+
a
i
j
+
σ
i
j
+
s
u
p
s
≤
t
e
r
s
w
i
j
s
δ
i
j
-
1for all ρ
≤
t, where i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m
. By (H3), we have(56)
e
r
ρ
w
i
j
ρ
δ
i
j
-
1
≤
W
2
∗
t
1
-
a
i
j
+
σ
i
j
+
.For any given t
≥
0, there exist k
∈
{
1,2
,
…
,
n
} and l
∈
{
1,2
,
…
,
m
} such that(57)
V
2
t
=
W
k
l
t
δ
k
l
.Calculating the derivative d
V
2
(
t
)
/
d
t along the positive half trajectory of (1) yields(58)
d
V
2
t
d
t
≤
-
a
k
l
t
-
r
V
2
t
+
1
δ
k
l
∫
t
-
σ
k
l
t
t
a
k
l
s
e
r
s
w
k
l
s
d
s
+
e
r
t
-
a
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
e
-
r
σ
k
l
t
×
w
k
l
t
-
σ
k
l
t
+
A
k
l
t
,
x
∗
t
,
w
t
≤
-
a
k
l
t
-
r
V
2
t
+
a
k
l
t
-
r
σ
k
l
+
a
k
l
+
V
2
∗
t
1
-
σ
k
l
+
a
k
l
+
+
a
k
l
t
e
r
σ
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
V
2
∗
t
1
-
σ
k
l
+
a
k
l
+
+
1
δ
k
l
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
L
e
r
τ
k
l
t
V
2
∗
t
1
-
a
k
l
+
σ
k
l
+
2
≤
-
a
k
l
t
-
r
1
-
2
σ
k
l
+
a
k
l
+
-
a
k
l
e
r
σ
k
l
t
-
1
-
σ
k
l
′
t
a
k
l
t
-
σ
k
l
t
δ
k
l
1
-
σ
k
l
+
a
k
l
+
-
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
t
L
e
r
τ
k
l
t
δ
k
l
1
-
a
k
l
+
σ
k
l
+
2
V
2
t
δ
k
l
≤
ς
V
2
t
δ
k
l
<
0when W
2
∗
(
t
)
≤
W
2
(
t
)
. Thus,(59)
d
V
2
∗
t
d
t
≤
0
,
∀
t
∈
R
.Then, for all i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m and for any t
>
0, we have(60)
w
i
j
t
≤
δ
i
j
V
2
∗
t
e
-
r
t
1
-
σ
i
j
+
a
i
j
+
≤
δ
i
j
V
2
∗
0
e
-
r
t
1
-
σ
i
j
+
a
i
j
+
.Thus, for any t
>
0, we have(61)
x
t
-
x
∗
t
≤
M
φ
-
ψ
e
-
r
t
,where (62)
M
=
∑
i
=
1
n
∑
j
=
1
m
δ
i
j
1
-
σ
i
j
+
a
i
j
+
δ
m
i
nand x
(
t
)
=
x
(
t
,
φ
) and x
∗
(
t
)
=
x
∗
(
t
,
ψ
)
. The proof of Theorem 11 is complete.


## 4. Periodicity and Stability

In this section, we will analyze the global exponential periodicity and stability of (1). By Theorem 11, we have the following results.Theorem 12.
In addition to (H1)–(H6), ifC
i
j
k
l
∗
(
t
)
,
C
i
j
k
l
∗
∗
(
t
)
,
L
i
j
(
t
)
,
a
i
j
(
t
)
,
τ
i
j
(
t
), and σ
i
j
(
t
) are ω-periodic functions, then there exists a unique ω-periodic solution x
∗
(
t
,
ψ
) for (1) which is globally exponentially stable.Next, we consider the following special form of (1):(63)
x
˙
i
j
t
=
-
a
i
j
x
i
j
t
-
σ
i
j
-
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
t
,
x
i
j
t
f
x
i
j
t
-
τ
i
j
x
i
j
t
+
L
i
j
,where i
=
1,2
,
…
,
n
,
j
=
1,2
,
…
,
m and C
i
j
k
l
∗
,
C
i
j
k
l
∗
∗
,
L
i
j
,
a
i
j, and τ
i
j are all constants and(64)
C
i
j
k
l
t
,
ξ
=
C
i
j
k
l
∗
t
,
h
i
j
ξ
>
T
i
j
,
C
i
j
k
l
∗
∗
t
,
h
i
j
ξ
<
T
i
j
,for t
∈
R,  a
i
j
(
t
,
x
)
=
C
i
j
k
l
∗
(
t
) or C
i
j
k
l
∗
∗
(
t
) when h
i
j
(
ξ
)
=
T
i
j
. Let E
∗
=
(
e
i
j
)
n
×
n, where (65)
e
i
j
=
1
1
-
a
i
j
+
σ
i
j
+
a
i
j
1
-
2
a
i
j
σ
i
j
-
∑
C
k
l
∈
N
ι
i
,
j
C
i
j
k
l
u
L
for
i
=
j
,
-
1
1
-
a
i
j
+
σ
i
j
+
∑
C
k
l
∈
N
r
i
,
j
C
i
j
k
l
u
L
for
i
≠
j
.By Theorem 11 and Lemma 8, we have the following result.Theorem 13.
In addition to (H1)–(H6), further assume that the matrixE
∗ is a nonsingular M-matrix. Then, there exists a unique equilibrium x
∗ for (63) which is globally exponentially stable; namely, for any other solution x
(
t
,
φ
) of (63), there exist constants M
,
r
>
0 such that x
(
t
,
φ
)
-
x
∗
≤
M
φ
-
ψ
e
-
r
t for all t
>
0.Remark 14.
In [1], Wu investigated the pseudo almost periodic solution of shunting inhibitory cellular neural networks with time-varying delay (without leakage delay). In [2], Chen and Cao studied the almost periodic solution of shunting inhibitory cellular neural networks with constant delay. In this paper, we study the almost periodic solutions for memristor-based shunting inhibitory cellular neural networks with leakage delays which is different from the work of [1, 2]. All the obtained results in [1, 2] cannot be applied to system (1) to obtain the exponential stability of almost periodic solutions for memristor-based shunting inhibitory cellular neural networks (1), which implies that our results obtained in this paper are completely new and complement the previous studies to some extent.


## 5. Examples

In this section, we present an example to verify the analytical predictions obtained in the previous section. Consider the following memristor-based shunting inhibitory cellular neural networks:(66)
x
˙
11
t
=
-
a
11
t
x
11
t
-
σ
11
t
-
∑
C
k
l
∈
N
1
1,1
C
11
k
l
t
,
x
11
t
f
x
11
t
-
τ
11
t
x
11
t
+
L
11
t
,
x
˙
12
t
=
-
a
12
t
x
12
t
-
σ
12
t
-
∑
C
k
l
∈
N
1
1,2
C
12
k
l
t
,
x
12
t
f
x
12
t
-
τ
12
t
x
12
t
+
L
12
t
,
x
˙
21
t
=
-
a
21
t
x
21
t
-
σ
21
t
-
∑
C
k
l
∈
N
1
2,1
C
21
k
l
t
,
x
21
t
f
x
21
t
-
τ
21
t
x
21
t
+
L
21
t
,
x
˙
22
t
=
-
a
22
t
x
22
t
-
σ
22
t
-
∑
C
k
l
∈
N
1
2,2
C
22
k
l
t
,
x
22
t
f
x
22
t
-
τ
22
t
x
22
t
+
L
22
t
,where g
f
(
x
)
=
tanh
⁡
x
-
1,  a
11
(
t
)
=
1
/
10,  a
12
(
t
)
=
1
/
12,  a
21
(
t
)
=
1
/
15,  a
22
(
t
)
=
1
/
20,  L
11
(
t
)
=
0.2
cos
⁡
2
t
+
sin
⁡
5
t,  L
12
(
t
)
=
0.1
cos
⁡
3
t
+
sin
⁡
7
t,  L
21
(
t
)
=
0.6
cos
⁡
11
t
+
sin
⁡
3
t,  L
22
(
t
)
=
0.4
cos
⁡
11
t
+
sin
⁡
13
t
,
τ
11
(
t
)
=
cos
2
⁡
t
/
18,  τ
12
(
t
)
=
s
i
n
2
t
/
15,  τ
21
(
t
)
=
s
i
n
2
t
/
17,  σ
22
(
t
)
=
s
i
n
2
t
/
15,  σ
11
(
t
)
=
s
i
n
2
t
/
7,  σ
12
(
t
)
=
c
o
s
2
t
/
9, σ
21
(
t
)
=
c
o
s
2
t
/
23,  σ
22
(
t
)
=
c
o
s
2
t
/
19, and(67)
C
11
k
l
t
,
ξ
=
0.4
cos
⁡
t
,
ξ
>
1
,
unchanged
,
ξ
=
1
,
0.3
sin
⁡
t
,
ξ
<
1
,
C
12
k
l
t
,
ξ
=
cos
⁡
5
t
,
ξ
>
1
,
unchanged
,
ξ
=
1
,
sin
⁡
3
t
,
ξ
<
1
,
C
21
k
l
t
,
ξ
=
0.1
cos
⁡
3
t
,
ξ
>
1
,
unchanged
,
ξ
=
1
,
0.4
sin
⁡
5
t
,
ξ
<
1
,
C
22
k
l
t
,
ξ
=
cos
⁡
7
t
,
ξ
>
1
,
unchanged
,
ξ
=
1
,
0.2
sin
⁡
t
,
ξ
<
1
.Then, (68)
Ł
=
1
,
m
a
x
1
≤
i
≤
n
,
1
≤
j
≤
m
a
i
j
+
σ
i
j
+
≤
1
70
<
1
. Let δ
11
=
0.2,  δ
12
=
0.3,  δ
21
=
0.32,  δ
22
=
0.41, and  ς
=
0.3
. Then, (69)
-
a
11
t
1
-
2
a
11
+
σ
11
+
-
a
11
t
-
1
-
σ
11
′
t
a
11
t
-
σ
11
t
δ
11
1
-
a
11
+
σ
11
+
+
∑
C
11
∈
N
1
1,1
C
11
11
u
t
L
δ
11
ϱ
/
ς
1
-
a
11
+
σ
11
+
2
=
-
0.3089
<
-
0.3
=
-
ς
,
-
a
12
t
1
-
2
a
12
+
σ
12
+
-
a
12
t
-
1
-
σ
12
′
t
a
12
t
-
σ
12
t
δ
12
1
-
a
12
+
σ
12
+
+
∑
C
12
∈
N
1
1,2
C
12
12
u
t
L
δ
12
ϱ
/
ς
1
-
a
12
+
σ
12
+
2
=
-
0.4125
<
-
0.3
=
-
ς
,
-
a
21
t
1
-
2
a
21
+
σ
21
+
-
a
21
t
-
1
-
σ
21
′
t
a
21
t
-
σ
21
t
δ
21
1
-
a
21
+
σ
21
+
+
∑
C
21
∈
N
1
2,1
C
21
21
u
t
L
δ
21
ϱ
/
ς
1
-
a
21
+
σ
21
+
2
=
-
0.3907
<
-
0.3
=
-
ς
,
-
a
22
t
1
-
2
a
22
+
σ
22
+
-
a
22
t
-
1
-
σ
22
′
t
a
22
t
-
σ
22
t
δ
22
1
-
a
22
+
σ
22
+
+
∑
C
22
∈
N
1
2,2
C
22
22
u
t
L
δ
22
ϱ
/
ς
1
-
a
22
+
σ
22
+
2
=
-
0.5128
<
-
0.3
=
-
ς
.Thus, all the conditions in Theorem 11 are satisfied. Then, we can conclude that system (66) has a unique almost periodic solution x
∗
(
t
), which is globally exponentially stable. These results are shown in Figures 1–4.Figure 1
The state trajectory(
t
,
x
11
(
t
)
) of system (66) with the initial values 0.1,0.2,0.5,0.3, 0.6,0.2,0.2,0.9, and 0.3,0.8,0.9,0.7.Figure 2
The state trajectory(
t
,
x
12
(
t
)
) of system (66) with the initial values 0.1,0.2,0.5,0.3, 0.6,0.2,0.2,0.9, and 0.3,0.8,0.9,0.7.Figure 3
The state trajectory(
t
,
x
21
(
t
)
) of system (66) with the initial values 0.1,0.2,0.5,0.3, 0.6,0.2,0.2,0.9, and 0.3,0.8,0.9,0.7.Figure 4
The state trajectory(
t
,
x
22
(
t
)
) of system (66) with the initial values 0.1,0.2,0.5,0.3, 0.6,0.2,0.2,0.9, and 0.3,0.8,0.9,0.7.


## 6. Conclusions

In the article, we have investigated a class of memristor-based shunting inhibitory cellular neural networks with leakage delays. Using the concept of the Filippov solution and differential inclusion, we study the dynamical nature of the memristor-based shunting inhibitory cellular neural networks with leakage delays. Applying a new Lyapunov function technique, a set of sufficient criteria which ensure the existence, uniqueness, and global exponential stability of almost periodic solution of the neural networks are established. The obtained results on the almost periodic solution are applied to prove the existence and stability of periodic solution for this neural network with periodic coefficients and leakage delays. An example is presented to illustrate the effectiveness of the theoretical findings. The almost periodic fluctuations can help us process visual information and predict pathological brain states. At present, pseudo almost periodic solutions of neural networks have also been paid more attention by many authors. However, very few results on pseudo almost periodic solutions of memristor-based neural networks with leakage delays have been reported, which might be our future research topic.



---

*Source: 3587271-2016-10-19.xml*