# Global Asymptotic Stability of Switched Neural Networks with Delays

Authors: Zhenyu Lu; Kai Li; Yan Li
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

This paper investigates the global asymptotic stability of a class of switched neural networks with delays. Several new criteria ensuring global asymptotic stability in terms of linear matrix inequalities (LMIs) are obtained via Lyapunov-Krasovskii functional. And here, we adopt the quadratic convex approach, which is different from the linear and reciprocal convex combinations that are extensively used in recent literature. In addition, the proposed results here are very easy to be verified and complemented. Finally, a numerical example is provided to illustrate the effectiveness of the results.


---

## Body

## 1. Introduction

In the past thirty years, neural networks have found extensive applications in associative memory, pattern recognition, and image processing [1–3]. It is true that most applications of neural networks are heavily dependent on the dynamic behaviors of neural networks, especially on global asymptotic stability of neural networks. On the other hand, time delays are inevitably encountered in the hardware implementation due to the finite switching speed of amplifier, which may destroy the system performance and become a source of oscillation or instability in neural networks. Therefore, stability of neural networks with delays has attracted increasing attention and lots of stability criteria have been reported in the literature [4, 5].As a special class of hybrid systems, switched systems are organized by a switching rule that orchestrates the switching. In reality, neural networks sometimes have finite modes that switch from one to another at different times according to a switching law. In [6, 7], the authors studied the stability problem of different kinds of switched neural networks with time delays. Different from the model in these works, in this paper, we consider a class of neural networks with state-dependent switchings. Our switched neural networks model is general and it generalizes the conventional neural networks.Recently, convex analysis has been significantly employed in the stability analysis of time-delay systems [8–17]. According to the feature of different convex functions, different convex combination approaches are adopted in the literature, such as the linear convex combination [8–10], reciprocal convex combination [11–13], and quadratic convex combination [15–17]. In [8, 9, 12, 14, 16, 17], convex combination technology was successfully used to derive some stability criteria for neural networks with time delays. It should be pointed out that the lower bound of the time delay in [8, 9, 16] is zero, which means the information on the lower bound of the time delay cannot be sufficiently used. Namely, the conditions obtained in [9, 14, 16] fail to take effect on the stability of neural networks when the lower bound of the time delay is strictly greater than zero.In this paper, some delay-dependent stability criteria in terms of LMIs are derived. The advantages are as follows. Firstly, differential inclusions and set-valued maps are employed to deal with the switched neural networks with discontinuous right-hand sides. Secondly, our results employ the quadratic convex approach, which is different from the linear and reciprocal convex combinations that are extensively used in recent literature on stability. Thirdly, the lower boundτ
1 of the time-varying delays is not zero and its information is adequately used to construct the Lyapunov-Krasovskii functional. Fourthly, we resort to neither Jensen’s inequality with delay-dividing approach nor the free-weighting matrix method compared with previous results.The organization of this paper is as follows. Some preliminaries are introduced in Section2. In Section 3, based on the quadratic convex approach, delay-dependent stability criteria in terms of LMIs are established for switched neural networks with time-varying delays. Then, an example is given to demonstrate the effectiveness of the obtained results in Section 4. Finally, conclusions are given in Section 5.Notations. Throughout this paper, R
n denotes the n-dimensional Euclidean space. A
T and A
-
1 denote the transpose and the inverse of the matrix A, respectively. A
>
0
(
A
≥
0
) means that the matrix A is symmetric and positive definite (semipositive definite). ∗ represents the elements below the main diagonal of a symmetric matrix. The identity and zero matrices of appropriate dimensions are denoted by I and 0, respectively. S
Y
M
(
A
) is defined as S
Y
M
(
A
)
=
A
+
A
T. d
i
a
g
{
⋯
} denotes a block-diagonal matrix.


## 2. System Description and Preliminaries

In this paper, we consider a class of switched neural networks with delays as follows:(1)
x
˙
i
t
=
-
d
i
x
i
t
x
i
t
+
∑
j
=
1
n
a
i
j
x
i
t
f
j
x
j
t
+
∑
j
=
1
n
b
i
j
x
i
t
f
j
x
j
t
-
τ
j
t
,
t
≥
0
,
i
=
1,2
,
…
,
n
,where (2)
d
i
x
i
t
=
d
i
∗
,
x
i
t
≤
0
,
d
i
∗
∗
,
x
i
t
>
0
,
a
i
j
x
i
t
=
a
i
j
∗
,
x
i
t
≤
0
,
a
i
j
∗
∗
,
x
i
t
>
0
,
b
i
j
x
i
t
=
b
i
j
∗
,
x
i
t
≤
0
,
b
i
j
∗
∗
,
x
i
t
>
0
,
x
i
(
t
) is the state variable of the ith neuron, and a
i
j
(
x
i
(
t
)
) and b
i
j
(
x
i
(
t
)
) denote the feedback connection weight and the delayed feedback connection weight, respectively. f
j
:
R
→
R is bounded continuous function; τ
j
(
t
) corresponds to the transmission delay and satisfies i
,
j
=
1,2
,
…
,
n. d
i
∗
>
0, d
i
∗
∗
>
0, a
i
j
∗
,
a
i
j
∗
∗
,
b
i
j
∗, and b
i
j
∗
∗
,
i
,
j
=
1,2
,
…
,
n, are all constant numbers. The initial condition of system (1) is x
(
s
)
=
ϕ
(
s
)
=
(
ϕ
1
(
s
)
,
ϕ
2
(
s
)
,
…
,
ϕ
n
(
s
)
)
T
∈
C
(
[
-
τ
2
,
0
]
,
R
n
).The following assumptions are given for system (1):(H1)
Forj
∈
1,2
,
…
,
n
,
f
j is bounded and there exist constants h
j
-
,
h
j
+ such that(3)
h
j
-
≤
f
j
s
1
-
f
j
s
2
s
1
-
s
2
≤
h
j
+
,
f
j
0
=
0
,
for alls
1
,
s
2
∈
R, s
1
≠
s
2.
(H2)
The transmission delayτ
j
(
t
) is a differential function and there exist constants 0
≤
τ
1
<
τ
2, μ such that (4)
0
≤
τ
1
≤
τ
j
t
≤
τ
2
,
τ
˙
j
t
≤
μ
,
for allt
≥
0, j
=
1,2
,
…
,
n.Obviously, system (1) is a discontinuous system; then its solution is different from the classic solution and cannot be defined in the conventional sense. In order to obtain the solution of system (1), some definitions and lemmas are given.Definition 1.
For a system with discontinuous right-hand sides,(5)
d
x
d
t
=
F
t
,
x
,
x
0
=
x
0
,
x
∈
R
n
,
t
≥
0
.A set-valued map is defined as (6)
Φ
t
,
x
=
⋂
δ
>
0
⋂
μ
N
=
0
c
o
¯
F
t
,
B
x
,
δ
∖
N
,where c
o
¯
[
E
] is the closure of the convex hull of set E, E
⊂
R
n, B
(
x
,
δ
)
=
{
y
:
y
-
x
<
δ
,
x
,
y
∈
R
n
,
δ
∈
R
+
}, and N
⊂
R
n, μ
(
N
) is Lebesgue measure of set N.
A solution in Filippov’s sense of system (5) with initial condition x
(
0
)
=
x
0
∈
R
n is an absolutely continuous function x
(
t
)
,
t
∈
[
0
,
T
]
,
T
>
0, which satisfy x
(
0
)
=
x
0 and differential inclusion: (7)
d
x
d
t
∈
Φ
t
,
x
,
for
a
.
a
.
t
∈
0
,
T
.
IfF
(
t
,
x
) is bounded, then the set-valued function Φ
(
t
,
x
) is upper semicontinuous with nonempty, convex, and compact values [18]. Then, the solution x
(
t
) of system (5) with initial condition exists and it can be extended to the interval [
0
,
+
∞
) in the sense of Filippov.
By applying the theories of set-valued maps and differential inclusions [18–20], system (1) can be rewritten as the following differential inclusion: (8)
x
˙
i
t
∈
-
d
_
i
,
d
¯
i
x
i
t
+
∑
j
=
1
n
a
_
i
j
,
a
¯
i
j
f
j
x
j
t
+
∑
j
=
1
n
b
_
i
j
,
b
¯
i
j
f
j
x
j
t
-
τ
j
t
,
for
a
.
a
.
t
≥
0
,
i
=
1,2
,
…
,
n
,where [
ξ
_
i
,
ξ
¯
i
] is the convex hull of [
ξ
_
i
,
ξ
¯
i
], ξ
_
i
,
ξ
¯
i
∈
R. d
_
i
=
min
⁡
{
d
i
∗
,
d
i
∗
∗
}, d
¯
i
=
max
⁡
{
d
i
∗
,
d
i
∗
∗
}, a
_
i
j
=
min
⁡
{
a
i
j
∗
,
a
i
j
∗
∗
}, a
¯
i
j
=
max
⁡
{
a
i
j
∗
,
a
i
j
∗
∗
}, b
_
i
j
=
min
⁡
{
b
i
j
∗
,
b
i
j
∗
∗
}, and b
¯
i
j
=
max
⁡
{
b
i
j
∗
,
b
i
j
∗
∗
}. The other parameters are the same as in system (1).Definition 2.
A constant vectorx
∗
=
(
x
1
∗
,
x
2
∗
,
…
,
x
n
∗
)
T is called an equilibrium point of system (1), if, for i
=
1,2
,
…
,
n, (9)
0
∈
-
d
_
i
,
d
¯
i
x
i
∗
+
∑
j
=
1
n
a
_
i
j
,
a
¯
i
j
f
j
x
j
∗
+
∑
j
=
1
n
b
_
i
j
,
b
¯
i
j
f
j
x
j
∗
.
It is easy to find that the origin(
0,0
,
…
,
0
)
T is an equilibrium point of system (1).Definition 3 (see [18]).
A functionx
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
,
…
,
x
n
(
t
)
)
T is a solution of (1), with the initial condition x
(
s
)
=
ϕ
(
s
)
=
(
ϕ
1
(
s
)
,
ϕ
2
(
s
)
,
…
,
ϕ
n
(
s
)
)
T
∈
C
(
[
-
τ
2
,
0
]
,
R
n
), if x
(
t
) is an absolutely continuous function and satisfies differential inclusion (8).Lemma 4 (see [18]).
Suppose that assumption (H1) is satisfied; then solutionx
(
t
) with initial condition ϕ
(
s
)
=
(
ϕ
1
(
s
)
,
ϕ
2
(
s
)
,
…
,
ϕ
n
(
s
)
)
T
∈
C
(
[
-
τ
2
,
0
]
,
R
n
) of (1) exists and it can be extended to the interval [
0
,
+
∞
).Before giving our main results, we present the following important lemmas that will be used in the proof to derive the stability conditions of the switched neural networks.Lemma 5 (see [21]).
Given constant matricesΣ
1
,
Σ
2, and Σ
3, where Σ
1
T
=
Σ
1, Σ
2
T
=
Σ
2, (10)
Σ
1
Σ
3
Σ
3
T
-
Σ
2
<
0is equivalent to the following conditions: (11)
Σ
2
>
0
,
Σ
1
+
Σ
3
Σ
2
-
1
Σ
3
T
<
0
.Lemma 6 (see [17]).
For real symmetric matricesW
1
,
W
2
,
W
3
∈
R
m
×
m and a scalar continuous function τ satisfy τ
1
≤
τ
≤
τ
2, where τ
1 and τ
2 are constants satisfying 0
≤
τ
1
≤
τ
2. If W
1
≥
0, then (12)
τ
2
W
1
+
τ
W
2
+
W
3
<
0
,
∀
τ
∈
τ
1
,
τ
2
⟺
τ
i
2
W
1
+
τ
i
W
2
+
W
3
<
0
,
i
=
1,2
.Lemma 7 (see [16]).
LetW
>
0, and let y
(
s
) be an appropriate dimensional vector. Then, one has the following facts for any scalar function η
(
s
)
≥
0, ∀
s
∈
[
t
1
,
t
2
]: (i)
(13)
-
∫
t
1
t
2
y
T
s
W
y
s
d
s
≤
t
2
-
t
1
ξ
t
T
F
1
T
W
-
1
F
1
ξ
t
+
2
ξ
t
T
F
1
T
∫
t
1
t
2
y
s
d
s
,
(ii)
(14)
-
∫
t
1
t
2
η
s
y
T
s
W
y
s
d
s
≤
∫
t
1
t
2
η
s
d
s
ξ
t
T
F
2
T
W
-
1
F
2
ξ
t
+
2
ξ
t
T
F
2
T
∫
t
1
t
2
η
s
y
s
d
s
,
(iii)
(15)
-
∫
t
1
t
2
η
2
s
y
T
s
W
y
s
d
s
≤
t
2
-
t
1
ξ
t
T
F
3
T
W
-
1
F
3
ξ
t
+
2
ξ
t
T
F
3
T
∫
t
1
t
2
η
s
y
s
d
s
,where matrices F
i  (i
=
1,2
,
3) and vector ξ
t independent of the integral variable are appropriate dimensional arbitrary ones.


## 3. Main Results

For presentation convenience, in the following we denoteτ
21
=
τ
2
-
τ
1, H
-
=
d
i
a
g
{
h
1
-
,
h
2
-
,
…
,
h
n
-
}, H
+
=
d
i
a
g
{
h
1
+
,
h
2
+
,
…
,
h
n
+
}, E
i
=
0
n
×
i
-
1
n
,
I
n
,
0
n
×
10
-
i
n
T
,
i
=
1,2
,
…
,
10. A
c
=
-
D
,
0,0
,
0,0
,
0,0
,
A
,
B
,
0
T, D
=
d
i
a
g
(
D
i
)
n
×
n, D
i
=
m
i
n
{
d
_
i
,
d
¯
i
}, A
=
(
A
i
j
)
n
×
n, A
i
j
=
m
a
x
{
a
_
i
j
,
a
¯
i
j
}, B
=
(
B
i
j
)
n
×
n, B
i
j
=
m
a
x
{
b
_
i
j
,
b
¯
i
j
}.Theorem 8.
Suppose assumptions (H1) and (H2) hold; the origin of system (1) is globally asymptotically stable if there exist matrices P
>
0, S
>
0, Q
i
>
0  (i
=
1,2
,
3,4
,
5), R
j
>
0  (j
=
1,2
,
3,4), and Γ
=
d
i
a
g
(
γ
1
,
γ
2
,
…
,
γ
n
)
>
0, K
1
=
d
i
a
g
(
k
11
,
k
12
,
…
,
k
1
n
)
>
0, K
2
=
d
i
a
g
(
k
21
,
k
22
,
…
,
k
2
n
)
>
0, and F
i  (i
=
1,2
,
…
,
12), such that the following two LMIs hold:(16)
Ξ
=
Ξ
11
τ
1
F
4
T
3
τ
1
F
6
T
τ
1
F
5
T
τ
21
F
1
T
3
τ
21
F
3
T
τ
21
F
2
T
τ
21
F
7
T
3
τ
21
F
9
T
τ
21
F
8
T
∗
-
τ
1
Q
3
0
0
0
0
0
0
0
0
∗
∗
-
τ
1
R
2
0
0
0
0
0
0
0
∗
∗
∗
-
R
1
0
0
0
0
0
0
∗
∗
∗
∗
-
τ
21
Q
3
0
0
0
0
0
∗
∗
∗
∗
∗
-
τ
21
R
2
0
0
0
0
∗
∗
∗
∗
∗
∗
-
R
1
0
0
0
∗
∗
∗
∗
∗
∗
∗
-
τ
21
Q
5
0
0
∗
∗
∗
∗
∗
∗
∗
∗
-
τ
21
R
4
0
∗
∗
∗
∗
∗
∗
∗
∗
∗
-
R
3
<
0
,
Π
=
Π
11
τ
2
F
4
T
3
τ
2
F
6
T
τ
2
F
5
T
τ
21
F
10
T
3
τ
21
F
12
T
τ
21
F
11
T
∗
-
τ
2
Q
3
0
0
0
0
0
∗
∗
-
τ
2
R
2
0
0
0
0
∗
∗
∗
-
R
1
0
0
0
∗
∗
∗
∗
-
τ
21
Q
5
0
0
∗
∗
∗
∗
∗
-
τ
21
R
4
0
∗
∗
∗
∗
∗
∗
-
R
3
<
0
,where Ξ
11
=
Δ
0
+
S
Y
M
(
(
τ
21
E
2
-
E
7
)
(
2
F
2
+
3
F
3
+
2
F
8
+
3
F
9
)
)
+
S
Y
M
(
(
τ
1
E
1
-
(
E
1
+
E
5
)
)
(
2
F
5
+
3
F
6
)
), Π
11 = Δ
0
+
S
Y
M
(
[
A
c
,
0
]
Q
1
τ
21
E
1
,
E
6
T
)
-
S
Y
M
(
E
7
(
2
F
2
+
3
F
3
+
2
F
8
+
3
F
9
)
)
+
S
Y
M
(
(
τ
2
E
1
-
(
E
5
+
E
6
)
)
(
2
F
5
+
3
F
6
)
) with (17)
Δ
0
=
S
Y
M
E
1
,
E
5
,
E
6
+
E
7
P
A
c
,
E
1
-
E
3
,
E
3
-
E
4
T
+
E
8
Γ
A
c
T
+
E
8
S
E
8
T
-
1
-
μ
E
9
S
E
9
T
+
E
1
,
E
3
Q
1
+
Q
4
E
1
,
E
3
T
-
1
-
μ
E
1
,
E
2
Q
1
E
1
,
E
2
T
+
S
Y
M
A
c
,
0
Q
2
τ
2
E
1
,
E
5
+
E
6
+
E
7
T
+
A
c
,
0
Q
4
τ
21
E
1
,
E
6
+
E
7
T
+
E
1
,
E
1
Q
2
E
1
,
E
1
T
-
E
1
,
E
4
Q
2
+
Q
4
E
1
,
E
4
T
+
τ
2
E
1
,
A
c
Q
3
E
1
,
A
c
T
+
τ
21
E
3
,
E
10
Q
5
E
3
,
E
10
T
+
A
c
τ
2
2
R
1
+
τ
2
3
R
2
A
c
T
+
E
10
τ
21
2
R
3
+
τ
21
3
R
4
E
10
T
+
S
Y
M
E
5
+
E
6
,
E
1
-
E
2
F
4
+
E
7
,
E
2
-
E
4
F
7
+
E
6
,
E
3
-
E
2
F
10
+
S
Y
M
E
8
K
1
H
+
E
1
T
+
E
1
H
-
K
1
E
8
T
+
E
9
K
2
H
+
E
2
T
+
E
2
H
-
K
2
E
9
T
-
S
Y
M
E
8
K
1
E
8
T
+
E
1
H
-
K
1
H
+
E
1
T
+
E
9
K
2
E
9
T
+
E
2
H
-
K
2
H
+
E
2
T
.Proof.
Define a vectorξ
t
∈
R
10
n as (18)
ξ
t
T
=
x
T
t
,
x
T
t
-
τ
t
,
x
T
t
-
τ
1
,
x
T
t
-
τ
2
,
∫
t
-
τ
1
t
x
T
s
d
s
,
∫
t
-
τ
t
t
-
τ
1
x
T
s
d
s
,
∫
t
-
τ
2
t
-
τ
t
x
T
s
d
s
,
f
T
x
t
,
f
T
x
t
-
τ
t
,
x
˙
T
t
-
τ
1
.
Consider the following Lyapunov-Krasovskii functional candidate as follows:(19)
V
t
=
V
1
t
+
V
2
t
+
V
3
t
+
V
4
t
,where (20)
V
1
t
=
ζ
T
t
P
ζ
t
,
V
2
t
=
2
∑
i
=
1
n
γ
i
∫
0
x
i
t
f
i
s
d
s
,
V
3
t
=
∫
t
-
τ
t
t
f
T
x
s
S
f
x
s
d
s
+
∫
t
-
τ
t
t
-
τ
1
x
T
t
x
T
s
Q
1
x
T
t
x
T
s
T
d
s
,
V
4
t
=
∫
t
-
τ
2
t
x
T
t
x
T
s
Q
2
x
T
t
x
T
s
T
+
τ
2
-
t
+
s
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
τ
2
-
t
+
s
2
x
˙
T
s
R
1
x
˙
s
+
τ
2
-
t
+
s
3
x
˙
T
s
R
2
x
˙
s
d
s
+
∫
t
-
τ
2
t
-
τ
1
x
T
t
x
T
s
Q
4
x
T
t
x
T
s
T
+
τ
2
-
t
+
s
x
T
s
x
˙
T
s
Q
5
x
T
s
x
˙
T
s
T
+
τ
2
-
t
+
s
2
x
˙
T
s
R
3
x
˙
s
+
τ
2
-
t
+
s
3
x
˙
T
s
R
4
x
˙
s
d
s
,where ζ
T
(
t
)
=
x
T
t
∫
t
-
τ
1
t
x
T
s
d
s
∫
t
-
τ
2
t
-
τ
1
x
T
s
d
s.
Calculating the time derivatives ofV
i
(
t
)
(
i
=
1,2
,
3,4
) along the trajectories of system (8), we obtain (21)
V
˙
1
t
=
2
ζ
T
t
P
ζ
˙
t
≤
2
x
T
t
P
-
D
x
t
+
A
f
x
t
+
B
f
x
t
-
τ
t
x
t
-
x
t
-
τ
1
x
t
-
τ
1
-
x
t
-
τ
2
≤
2
ξ
t
T
E
1
,
E
5
,
E
6
+
E
7
P
A
c
,
E
1
-
E
3
,
E
3
-
E
4
T
ξ
t
=
ξ
t
T
S
Y
M
E
1
,
E
5
,
E
6
+
E
7
P
A
c
,
E
1
-
E
3
,
E
3
-
E
4
T
ξ
t
,
V
˙
2
t
=
2
∑
i
=
1
n
γ
i
f
i
x
i
t
x
˙
t
=
2
f
T
x
t
Γ
x
˙
t
≤
2
ξ
t
T
E
8
Γ
A
c
T
ξ
t
,
V
˙
3
t
=
f
T
x
t
S
f
x
t
-
1
-
τ
˙
t
f
T
x
t
-
τ
t
S
f
x
t
-
τ
t
+
x
T
t
x
T
t
-
τ
1
Q
1
x
T
t
x
T
t
-
τ
1
T
-
1
-
τ
˙
t
x
T
t
x
T
t
-
τ
t
Q
1
x
T
t
x
T
t
-
τ
t
T
+
2
∫
t
-
τ
t
t
-
τ
1
∂
x
T
t
x
T
s
∂
t
Q
1
x
T
t
x
T
s
T
d
s
≤
f
T
x
t
S
f
x
t
-
1
-
μ
f
T
x
t
-
τ
t
S
f
x
t
-
τ
t
+
x
T
t
x
T
t
-
τ
1
Q
1
x
T
t
x
T
t
-
τ
1
T
-
1
-
μ
x
T
t
x
T
t
-
τ
t
Q
1
x
T
t
x
T
t
-
τ
t
T
+
2
x
˙
T
t
,
0
Q
1
∫
t
-
τ
t
t
-
τ
1
x
T
t
x
T
s
T
d
s
≤
ξ
t
T
E
8
S
E
8
T
-
1
-
μ
E
9
S
E
9
T
+
E
1
,
E
3
Q
1
E
1
,
E
3
T
-
1
-
μ
E
1
,
E
2
Q
1
E
1
,
E
2
T
+
2
A
c
,
0
Q
1
τ
t
-
τ
1
E
1
,
E
6
T
ξ
t
,
V
˙
4
t
=
x
T
t
x
T
t
Q
2
x
T
t
x
T
t
T
+
x
T
t
x
T
t
-
τ
1
Q
4
x
T
t
x
T
t
-
τ
1
T
-
x
T
t
x
T
t
-
τ
2
Q
2
x
T
t
x
T
t
-
τ
2
T
-
x
T
t
x
T
t
-
τ
2
Q
4
x
T
t
x
T
t
-
τ
2
T
+
2
x
˙
T
t
,
0
Q
2
∫
t
-
τ
2
t
x
T
t
x
T
s
T
d
s
+
2
x
˙
T
t
,
0
Q
4
∫
t
-
τ
2
t
-
τ
1
x
T
t
x
T
s
T
d
s
+
τ
2
x
T
t
x
˙
T
t
Q
3
x
T
t
x
˙
T
t
T
+
τ
21
x
T
t
-
τ
1
x
˙
T
t
-
τ
1
Q
5
x
T
t
-
τ
1
x
˙
T
t
-
τ
1
T
+
x
˙
T
t
τ
2
2
R
1
+
τ
2
3
R
2
x
˙
t
+
x
˙
T
t
-
τ
1
τ
21
2
R
3
+
τ
21
3
R
4
x
˙
t
-
τ
1
+
V
x
+
V
y
=
ξ
t
T
E
1
,
E
1
Q
2
E
1
,
E
1
T
+
E
1
,
E
3
Q
4
E
1
,
E
3
T
-
E
1
,
E
4
Q
2
+
q
4
E
1
,
E
4
T
+
2
A
c
,
0
Q
2
τ
2
E
1
,
E
5
+
E
6
+
E
7
T
+
2
A
c
,
0
Q
4
τ
2
-
τ
1
E
1
,
E
6
+
E
7
T
+
τ
2
E
1
,
A
c
Q
3
E
1
,
A
c
T
+
τ
21
E
3
,
E
10
Q
5
E
3
,
E
10
T
+
A
c
τ
2
2
R
1
+
τ
2
3
R
2
A
c
T
+
E
10
τ
21
2
R
3
+
τ
21
3
R
4
E
10
T
ξ
t
+
V
x
+
V
y
,where(22)
V
x
=
-
∫
t
-
τ
2
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
,
V
y
=
-
∫
t
-
τ
2
t
-
τ
1
x
T
s
x
˙
T
s
Q
5
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
3
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
4
x
˙
s
d
s
.
It is easy to obtain the following identities:(23)
τ
2
-
t
+
s
=
τ
t
-
t
+
s
+
τ
2
-
τ
t
,
τ
2
-
t
+
s
2
=
τ
t
-
t
+
s
2
+
τ
2
2
-
τ
2
t
+
2
τ
2
-
τ
t
s
-
t
.
Thus, we have(24)
V
x
=
-
∫
t
-
τ
2
t
-
τ
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
∫
t
-
τ
t
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
=
-
∫
t
-
τ
2
t
-
τ
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
∫
t
-
τ
t
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
t
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
t
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
∫
t
-
τ
t
t
x
˙
T
s
2
τ
2
-
τ
t
R
1
+
3
τ
2
2
-
τ
2
t
R
2
+
6
τ
2
-
τ
t
s
-
t
R
2
x
˙
s
d
s
≤
-
∫
t
-
τ
2
t
-
τ
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
∫
t
-
τ
t
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
t
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
t
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
2
∫
t
-
τ
t
t
x
˙
T
s
τ
2
-
τ
t
R
1
x
˙
s
d
s
-
3
∫
t
-
τ
t
t
x
˙
T
s
τ
2
-
τ
t
2
R
2
x
˙
s
d
s
≤
-
∫
t
-
τ
2
t
-
τ
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
-
∫
t
-
τ
t
t
x
T
s
x
˙
T
s
Q
3
x
T
s
x
˙
T
s
T
+
2
τ
t
-
t
+
s
x
˙
T
s
R
1
x
˙
s
+
3
τ
t
-
t
+
s
2
x
˙
T
s
R
2
x
˙
s
d
s
.
Similarly, we have(25)
V
y
≤
-
∫
t
-
τ
2
t
-
τ
t
x
T
s
x
˙
T
s
Q
5
x
T
s
x
˙
T
s
T
+
2
τ
2
-
t
+
s
x
˙
T
s
R
3
x
˙
s
+
3
τ
2
-
t
+
s
2
x
˙
T
s
R
4
x
˙
s
d
s
-
∫
t
-
τ
t
t
-
τ
1
x
T
s
x
˙
T
s
Q
5
x
T
s
x
˙
T
s
T
+
2
τ
t
-
t
+
s
x
˙
T
s
R
3
x
˙
s
+
3
τ
t
-
t
+
s
2
x
˙
T
s
R
4
x
˙
s
d
s
.
Applying Lemma7 to V
x and V
y, we get (26)
V
x
≤
ξ
t
T
τ
2
-
τ
t
F
1
T
Q
3
-
1
F
1
+
2
F
1
T
E
7
,
E
2
-
E
4
T
+
τ
2
-
τ
t
2
F
2
T
R
1
-
1
F
2
+
4
F
2
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
3
T
R
2
-
1
F
3
+
6
F
3
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
F
4
T
Q
3
-
1
F
4
+
2
F
4
T
E
5
+
E
6
,
E
1
-
E
2
T
+
τ
2
t
F
5
T
R
1
-
1
F
5
+
4
F
5
T
τ
t
E
1
-
E
5
+
E
6
T
+
3
τ
t
F
6
T
R
2
-
1
F
6
+
6
F
6
T
τ
t
E
1
-
E
5
+
E
6
T
ξ
t
,
V
y
≤
ξ
t
T
τ
2
-
τ
t
F
7
T
Q
5
-
1
F
7
+
2
F
7
T
E
7
,
E
2
-
E
4
T
+
τ
2
-
τ
t
2
F
8
T
R
3
-
1
F
8
+
4
F
8
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
9
T
R
4
-
1
F
9
+
6
F
9
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
-
τ
1
F
10
T
Q
5
-
1
F
10
+
2
F
10
T
E
6
,
E
3
-
E
2
T
+
τ
t
-
τ
1
2
F
11
T
R
3
-
1
F
11
+
4
F
11
T
τ
t
-
τ
1
E
3
-
E
6
T
+
3
τ
t
-
τ
1
F
12
T
R
4
-
1
F
12
+
6
F
12
T
τ
t
-
τ
1
E
3
-
E
6
T
ξ
t
.
On the other hand, from assumption (H1), we have(27)
-
2
f
i
x
i
t
-
h
i
-
x
i
t
K
1
f
i
x
i
t
-
h
i
+
x
i
t
≥
0
,
-
2
f
i
x
i
t
-
τ
t
-
h
i
-
x
i
t
-
τ
t
K
1
f
i
x
i
t
-
τ
t
-
h
i
+
x
i
t
-
τ
t
≥
0
,where K
1
=
d
i
a
g
(
k
11
,
k
12
,
…
,
k
1
n
)
>
0
,
K
2
=
d
i
a
g
(
k
21
,
k
22
,
…
,
k
2
n
)
>
0
,
i
=
1,2
,
…
,
n.
From (27), it is easy to get that (28)
2
ξ
t
T
-
E
8
K
1
E
8
T
+
E
8
K
1
H
+
E
1
T
+
E
1
H
-
K
1
E
8
T
-
E
1
H
-
K
1
H
+
E
1
T
-
E
9
K
2
E
9
T
+
E
9
K
2
H
+
E
2
T
+
E
2
H
-
K
2
E
9
T
-
E
2
H
-
K
2
H
+
E
2
T
≥
0
.
It follows from (19)–(26) and (28) that (29)
V
˙
t
≤
ξ
t
T
2
E
1
,
E
5
,
E
6
+
E
7
P
A
c
,
E
1
-
E
3
,
E
3
-
E
4
T
+
2
E
8
Γ
A
c
T
+
E
8
S
E
8
T
-
1
-
μ
E
9
S
E
9
T
+
E
1
,
E
3
Q
1
E
1
,
E
3
T
-
1
-
μ
E
1
,
E
2
Q
1
E
1
,
E
2
T
+
2
A
c
,
0
Q
1
τ
t
-
τ
1
E
1
,
E
6
T
+
E
1
,
E
1
Q
2
E
1
,
E
1
T
+
E
1
,
E
3
Q
4
E
1
,
E
3
T
-
E
1
,
E
4
Q
2
+
q
4
E
1
,
E
4
T
+
2
A
c
,
0
Q
2
τ
2
E
1
,
E
5
+
E
6
+
E
7
T
+
2
A
c
,
0
Q
4
τ
21
E
1
,
E
6
+
E
7
T
+
τ
2
E
1
,
A
c
Q
3
E
1
,
A
c
T
+
τ
21
E
3
,
E
10
Q
5
E
3
,
E
10
T
+
A
c
τ
2
2
R
1
+
τ
2
3
R
2
A
c
T
+
E
10
τ
21
2
R
3
+
τ
21
3
R
4
E
10
T
+
τ
2
-
τ
t
F
1
T
Q
3
-
1
F
1
+
2
F
1
T
E
7
,
E
2
-
E
4
T
+
τ
2
-
τ
t
2
F
2
T
R
1
-
1
F
2
+
4
F
2
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
3
T
R
2
-
1
F
3
+
6
F
3
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
F
4
T
Q
3
-
1
F
4
+
2
F
4
T
E
5
+
E
6
,
E
1
-
E
2
T
+
τ
2
t
F
5
T
R
1
-
1
F
5
+
4
F
5
T
τ
t
E
1
-
E
5
+
E
6
T
+
3
τ
t
F
6
T
R
2
-
1
F
6
+
6
F
6
T
τ
t
E
1
-
E
5
+
E
6
T
+
τ
2
-
τ
t
F
7
T
Q
5
-
1
F
7
+
2
F
7
T
E
7
,
E
2
-
E
4
T
+
τ
2
-
τ
t
2
F
8
T
R
3
-
1
F
8
+
4
F
8
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
9
T
R
4
-
1
F
9
+
6
F
9
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
-
τ
1
F
10
T
Q
5
-
1
F
10
+
2
F
10
T
E
6
,
E
3
-
E
2
T
+
τ
t
-
τ
1
2
F
11
T
R
3
-
1
F
11
+
4
F
11
T
τ
t
-
τ
1
E
3
-
E
6
T
+
3
τ
t
-
τ
1
F
12
T
R
4
-
1
F
12
+
6
F
12
T
τ
t
-
τ
1
E
3
-
E
6
T
-
2
E
8
K
1
E
8
T
+
2
E
8
K
1
H
+
E
1
T
+
2
E
1
H
-
K
1
E
8
T
-
2
E
1
H
-
K
1
H
+
E
1
T
-
2
E
9
K
2
E
9
T
+
2
E
9
K
2
H
+
E
2
T
+
2
E
2
H
-
K
2
E
9
T
-
2
E
2
H
-
K
2
H
+
E
2
T
ξ
t
=
ξ
t
T
Δ
0
+
Δ
1
ξ
t
,where Δ
0 is defined in the theorem context, and (30)
Δ
1
=
2
A
c
,
0
Q
1
τ
t
-
τ
1
E
1
,
E
6
T
+
τ
2
-
τ
t
F
1
T
Q
3
-
1
F
1
+
τ
2
-
τ
t
2
F
2
T
R
1
-
1
F
2
+
4
F
2
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
3
T
R
2
-
1
F
3
+
6
F
3
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
F
4
T
Q
3
-
1
F
4
+
τ
2
t
F
5
T
R
1
-
1
F
5
+
4
F
5
T
τ
t
E
1
-
E
5
+
E
6
T
+
3
τ
t
F
6
T
R
2
-
1
F
6
+
6
F
6
T
τ
t
E
1
-
E
5
+
E
6
T
+
τ
2
-
τ
t
F
7
T
Q
5
-
1
F
7
+
τ
2
-
τ
t
2
F
8
T
R
3
-
1
F
8
+
4
F
8
T
τ
2
-
τ
t
E
2
-
E
7
T
+
3
τ
2
-
τ
t
F
9
T
R
4
-
1
F
9
+
6
F
9
T
τ
2
-
τ
t
E
2
-
E
7
T
+
τ
t
-
τ
1
F
10
T
Q
5
-
1
F
10
+
τ
t
-
τ
1
2
F
11
T
R
3
-
1
F
11
+
4
F
11
T
τ
t
-
τ
1
E
3
-
E
6
T
+
3
τ
t
-
τ
1
F
12
T
R
4
-
1
F
12
+
6
F
12
T
τ
t
-
τ
1
E
3
-
E
6
T
.
It is easy to see thatΔ
0
+
Δ
1 is a quadratic convex combination of matrices on τ
(
t
)
∈
[
τ
1
,
τ
2
].
Applying Lemma5 to (16), we have (31)
Δ
0
+
Δ
1
τ
t
=
τ
1
=
Δ
0
+
Δ
1
τ
t
=
τ
1
=
Δ
0
+
2
τ
21
E
2
-
E
7
2
F
2
+
3
F
3
+
2
F
8
+
3
F
9
+
2
τ
1
E
1
-
E
1
+
E
5
2
F
5
+
3
F
6
+
τ
1
F
4
T
Q
3
-
1
F
4
+
3
F
6
T
R
2
-
1
F
6
+
τ
1
2
F
5
T
R
1
-
1
F
5
+
τ
21
F
1
T
Q
3
-
1
F
1
+
F
7
T
Q
5
-
1
F
7
+
3
F
3
T
R
2
-
1
F
3
+
3
F
9
T
R
4
-
1
F
9
+
τ
21
2
F
2
T
R
1
-
1
F
2
+
F
8
T
R
3
-
1
F
8
<
0
,
Δ
0
+
Δ
1
τ
t
=
τ
2
=
Δ
0
+
Δ
1
τ
t
=
τ
2
=
Δ
0
+
2
A
c
,
0
Q
1
τ
21
E
1
,
E
6
T
-
2
E
7
2
F
2
+
3
F
3
+
2
F
8
+
3
F
9
+
2
τ
2
E
1
-
E
5
+
E
6
2
F
5
+
3
F
6
+
τ
2
F
4
T
Q
3
-
1
F
4
+
3
F
6
T
R
2
-
1
F
6
+
τ
2
2
F
5
T
R
1
-
1
F
5
<
0
.Since ξ
t
T
{
F
2
T
R
1
-
1
F
2
+
F
5
T
R
1
-
1
F
5
+
F
8
T
R
3
-
1
F
8
+
F
11
T
R
3
-
1
F
11
}
ξ
t
≥
0, from Lemma 6, if LMIs (16) are true, then Δ
0
+
Δ
1
<
0
,
∀
τ
(
t
)
∈
[
τ
1
,
τ
2
]. Then, we can see that the origin of system (1) is asymptotically stable.
The proof is completed.Remark 9.
In [16], stability analysis for neural networks with time-varying delay is studied by using quadratic convex combination. Our results have two advantages compared with the results of that paper. On the one hand, the information on the lower bound τ
1 of the time-varying delays is considered. On the other hand, the augmented vector ζ
(
t
) includes the distributed delay terms.Remark 10.
We use three inequalities in Lemma7 combined with the quadratic convex combination implied by Lemma 6, rather than Jensen’s inequality and the linear convex combination. In addition, our theoretical proof is not concerned with free-weighting matrix method.Remark 11.
To use the quadratic convex approach, we construct the Lyapunov-Krasovskii functional with the following term:∫
t
-
τ
2
t
∑
j
=
1
3
(
τ
2
+
t
-
s
)
j
g
(
s
)
d
s. The degree increase of τ
2
-
t
+
s by 1 means the number increase of the integral by 1 due to the fact that ∫
τ
2
t
∫
t
+
θ
t
g
s
d
s
d
θ
=
∫
t
-
τ
2
t
(
τ
2
-
t
+
s
)
g
(
s
)
d
s.
In the caseτ
1
=
0, we have the following result.Corollary 12.
Suppose assumptions (H1) and (H2) hold withτ
1
=
0; the origin of system (1) is globally asymptotically stable if there exist matrices P
>
0, S
>
0, Q
i
>
0  (i
=
1,2
,
3), R
j
>
0  (j
=
1,2), and Γ
=
d
i
a
g
(
γ
1
,
γ
2
,
…
,
γ
n
)
>
0, K
1
=
d
i
a
g
(
k
11
,
k
12
,
…
,
k
1
n
)
>
0, K
2
=
d
i
a
g
(
k
21
,
k
22
,
…
,
k
2
n
)
>
0, and F
i  (i
=
1,2
,
…
,
6), such that the following two LMIs hold: (32)
Ξ
=
Δ
0
τ
2
F
1
T
3
τ
2
F
3
T
τ
2
F
2
T
∗
-
τ
2
Q
3
0
0
∗
∗
-
τ
2
R
2
0
∗
∗
∗
-
R
1
<
0
,
Π
=
Π
11
τ
2
F
4
T
3
τ
2
F
6
T
τ
2
F
5
T
∗
-
τ
2
Q
3
0
0
∗
∗
-
τ
2
R
2
0
∗
∗
∗
-
R
1
<
0
,where Π
11
=
Δ
0
+
τ
2
S
Y
M
(
[
A
c
,
0
]
Q
1
E
1
,
0
T
)
-
τ
2
S
Y
M
(
E
2
2
F
2
+
3
F
3
)
+
τ
2
S
Y
M
(
E
1
(
2
F
5
+
3
F
6
)
) with (33)
Δ
0
=
S
Y
M
E
1
,
E
4
+
E
5
P
A
c
,
E
1
-
E
3
T
+
E
6
Γ
A
c
T
+
E
6
S
E
6
T
-
1
-
μ
E
7
S
E
7
T
+
E
1
,
E
1
Q
1
+
Q
2
E
1
,
E
1
T
-
1
-
μ
E
1
,
E
2
Q
1
E
1
,
E
2
T
-
E
1
,
E
3
Q
2
E
1
,
E
3
T
+
S
Y
M
A
c
,
0
Q
1
0
,
E
5
T
+
A
c
,
0
Q
2
τ
2
E
1
,
E
4
+
E
5
T
+
τ
2
E
1
,
A
c
Q
3
E
1
,
A
c
T
+
A
c
τ
2
2
R
1
+
τ
2
3
R
2
A
c
T
+
S
Y
M
E
5
,
E
1
-
E
2
F
1
+
τ
2
E
2
-
E
5
2
F
2
+
3
F
3
+
E
4
,
E
1
-
E
2
F
4
-
S
Y
M
E
4
2
F
5
+
3
F
6
+
S
Y
M
E
8
K
1
H
+
E
1
T
+
E
1
H
-
K
1
E
8
T
+
E
9
K
2
H
+
E
2
T
+
E
2
H
-
K
2
E
9
T
-
S
Y
M
E
8
K
1
E
8
T
+
E
1
H
-
K
1
H
+
E
1
T
+
E
9
K
2
E
9
T
+
E
2
H
-
K
2
H
+
E
2
T
.In addition, when the information of the time derivative of delays is unknown or the derivative of the delays does not exist, then we have the following result.Corollary 13.
Suppose assumption (H1) holds withτ
1
=
0; the origin of system (1) is globally asymptotically stable if there exist matrices P
>
0, Q
i
>
0  (i
=
2,3), R
j
>
0  (j
=
1,2), and Γ
=
d
i
a
g
(
γ
1
,
γ
2
,
…
,
γ
n
)
>
0, K
1
=
d
i
a
g
(
k
11
,
k
12
,
…
,
k
1
n
)
>
0, K
2
=
d
i
a
g
(
k
21
,
k
22
,
…
,
k
2
n
)
>
0, and F
i  (i
=
1,2
,
…
,
6), such that the following two LMIs hold:(34)
Ξ
=
Δ
0
τ
2
F
1
T
3
τ
2
F
3
T
τ
2
F
2
T
∗
-
τ
2
Q
3
0
0
∗
∗
-
τ
2
R
2
0
∗
∗
∗
-
R
1
<
0
,
Π
=
Π
11
τ
2
F
4
T
3
τ
2
F
6
T
τ
2
F
5
T
∗
-
τ
2
Q
3
0
0
∗
∗
-
τ
2
R
2
0
∗
∗
∗
-
R
1
<
0
,where Π
11
=
Δ
0
+
τ
2
S
Y
M
(
[
A
c
,
0
]
Q
1
E
1
,
0
T
)
-
τ
2
S
Y
M
(
E
2
2
F
2
+
3
F
3
)
+
τ
2
S
Y
M
E
1
2
F
5
+
3
F
6 with (35)
Δ
0
=
S
Y
M
E
1
,
E
4
+
E
5
P
A
c
,
E
1
-
E
3
T
+
E
6
Γ
A
c
T
+
E
1
,
E
1
Q
2
E
1
,
E
1
T
-
E
1
,
E
3
Q
2
E
1
,
E
3
T
+
S
Y
M
A
c
,
0
Q
1
0
,
E
5
T
+
A
c
,
0
Q
2
τ
2
E
1
,
E
4
+
E
5
T
+
τ
2
E
1
,
A
c
Q
3
E
1
,
A
c
T
+
A
c
τ
2
2
R
1
+
τ
2
3
R
2
A
c
T
+
S
Y
M
E
5
,
E
1
-
E
2
F
1
+
τ
2
E
2
-
E
5
2
F
2
+
3
F
3
+
E
4
,
E
1
-
E
2
F
4
-
S
Y
M
E
4
2
F
5
+
3
F
6
+
S
Y
M
E
8
K
1
H
+
E
1
T
+
E
1
H
-
K
1
E
8
T
+
E
9
K
2
H
+
E
2
T
+
E
2
H
-
K
2
E
9
T
-
S
Y
M
E
8
K
1
E
8
T
+
E
1
H
-
K
1
H
+
E
1
T
+
E
9
K
2
E
9
T
+
E
2
H
-
K
2
H
+
E
2
T
.Remark 14.
It is worth noting that when we consider system (1) without switching, that is, d
i
∗
=
d
i
∗
∗, a
i
j
∗
=
a
i
j
∗
∗, and b
i
j
∗
=
b
i
j
∗
∗, then Corollary 12 is the main Theorem 1 of [16].Remark 15.
Compared with the results on stability of neural networks with continuous right-hand side [6], our results on stability of neural networks are with discontinuous right-hand sides. So the results of this paper are less conservative and more general.


## 4. Numerical Example

In this section, an example is provided to verify the effectiveness of the results obtained in the previous section.Example 1.
Consider two-dimensional switched neural networks with time-varying delays(36)
x
˙
i
t
=
-
d
i
x
i
t
x
i
t
+
∑
j
=
1
2
a
i
j
x
i
t
f
j
x
j
t
+
∑
j
=
1
2
b
i
j
x
i
t
f
j
x
j
t
-
τ
j
t
,
t
≥
0
,
i
=
1,2
,where (37)
d
1
x
1
t
=
2
,
x
1
t
≤
0
,
2.1
,
x
1
t
>
0
,
d
2
x
2
t
=
2.1
,
x
2
t
≤
0
,
2
,
x
2
t
>
0
,
a
11
x
1
t
=
1
,
x
1
t
≤
0
,
0.9
,
x
1
t
>
0
,
a
12
x
1
t
=
1
,
x
1
t
≤
0
,
1.2
,
x
1
t
>
0
,
a
21
x
2
t
=
-
0.9
,
x
2
t
≤
0
,
-
1
,
x
2
t
>
0
,
a
22
x
2
t
=
-
1.2
,
x
2
t
≤
0
,
-
1
,
x
2
t
>
0
,
b
11
x
1
t
=
-
0.7
,
x
1
t
≤
0
,
-
0.8
,
x
1
t
>
0
,
b
12
x
1
t
=
1
,
x
1
t
≤
0
,
1.2
,
x
1
t
>
0
,
b
21
x
2
t
=
1
,
x
2
t
≤
0
,
0.9
,
x
2
t
>
0
,
b
22
x
2
t
=
-
1
,
x
2
t
≤
0
,
-
1.2
,
x
2
t
>
0
,with τ
j
(
t
)
=
0.5
cos
⁡
(
1.6
t
)
+
3
,
j
=
1,2, and take the activation function as f
(
x
)
=
sin
⁡
x
1
,
sin
⁡
x
2
T. We can obtain that τ
1
=
2.5, τ
2
=
3.5, τ
21
=
1, μ
=
0.8, (38)
D
=
2
0
0
2
,
A
=
1
1.2
1
1.2
,
B
=
0.8
1.2
1
1.2
H
-
=
-
1
0
0
-
1
,
H
+
=
1
0
0
1
.
For written simplification, we use the Matlab LMI Control Toolbox and then a solution to LMIs (32) is obtained (
F
i
,
i
=
1
,
…
,
6, is omitted since the dimension is too big) as follows: (39)
P
=
28.5171
-
20.9241
0.5484
-
0.4280
-
20.9241
27.7701
-
0.4064
0.4984
0.5484
-
0.4064
10.8451
-
3.9493
-
0.4280
0.4984
-
3.9493
10.6288
,
Q
1
=
7.3594
-
1.0832
0.0588
-
0.0060
-
1.0832
7.4606
-
0.0074
0.0723
0.0588
-
0.0074
3.5541
-
0.4045
-
0.0060
0.0723
-
0.4045
3.6751
,
Q
2
=
67.4363
-
64.0296
3.2893
-
3.0414
-
64.0296
65.2276
-
3.0862
0.0723
3.2893
-
3.0862
35.1371
-
5.2829
-
3.0414
0.0723
-
5.2829
34.8395
,
Q
3
=
15.7736
-
2.2877
0.8863
-
0.8539
-
2.2877
15.6801
-
0.8085
0.4984
0.8863
-
0.8085
7.6796
-
6.7268
-
0.8539
0.4984
-
6.7268
7.3682
,
S
=
2.4075
1.0722
1.0722
2.5766
,
G
=
3.5681
0
0
3.7071
,
K
1
=
57.4100
0
0
65.2979
,
K
2
=
26.5401
0
0
29.9701
.
Therefore, according to Corollary12, we can see that the origin of system (36) is globally asymptotically stable. The state trajectories of variables x
1
(
t
) and x
2
(
t
) are shown in Figure 1.Figure 1
The state curves of system (36).Remark 16.
Because the parameters of system (1) are discontinuous, the results obtained in [6] about neural networks with continuous right-hand sides cannot be used here. In addition, the lower bounds of the delays of system (36) are not zero, so the results obtained in [8, 9, 16] cannot be used here.


## 5. Conclusions

In this paper, the delay-dependent stability for a class of switched neural networks with time-varying delays has been studied by using the quadratic convex combination. Some delay-dependent criteria in terms of LMIs have been obtained. The lower boundτ
1 of the time-varying delays is considered to be nonzero so that the information of τ
1 can be used adequately. It is worth noting that we resort to neither Jensen’s inequality with delay-dividing approach nor the free-weighting matrix method compared with previous results.



---

*Source: 717513-2015-12-29.xml*