Skip to main content

Global exponential stability and existence of periodic solutions for delayed reaction-diffusion BAM neural networks with Dirichlet boundary conditions

Abstract

In this paper, both global exponential stability and periodic solutions are investigated for a class of delayed reaction-diffusion BAM neural networks with Dirichlet boundary conditions. By employing suitable Lyapunov functionals, sufficient conditions of the global exponential stability and the existence of periodic solutions are established for reaction-diffusion BAM neural networks with mixed time delays and Dirichlet boundary conditions. The derived criteria extend and improve previous results in the literature. A numerical example is given to show the effectiveness of the obtained results.

1 Introduction

Neural networks (NNs) have been extensively studied in the past few years and have found many applications in different areas such as pattern recognition, associative memory, combinatorial optimization, etc. Delayed versions of NNs were also proved to be important for solving certain classes of motion-related optimization problems. Various results concerning the dynamical behavior of NNs with delays have been reported during the last decade (see, e.g., [17]). Recently, the authors in [1] and [2] considered the problem of exponential passivity analysis for uncertain NNs with time-varying delays and passivity-based controller design for Hopfield NNs, respectively.

Since NNs related to bidirectional associative memory (BAM) were proposed by Kosko [8], the BAM NNs have been one of the most interesting research topics and have attracted the attention of researchers. In the design and applications of networks, the stability of the designed BAM NNs is one of the most important issues (see, e.g., [912]). Many important results concerning mainly the existence and stability of equilibrium of BAM NNs have been obtained (see, e.g., [915]).

However, strictly speaking, diffusion effects cannot be avoided in the NNs when electrons are moving in asymmetric electromagnetic fields. So, we must consider that the activations vary in space as well as in time. In [1634], the authors considered the stability of NNs with diffusion terms which were expressed by partial differential equations. In particular, the existence and attractivity of periodic solutions for non-autonomous reaction-diffusion Cohen-Grossberg NNs with discrete time delays were investigated in [20]. The authors derived sufficient conditions on the stability and periodic solutions of delayed reaction-diffusion NNs (RDNNs) with Neumann boundary conditions in [2125]. In these works, due to the divergence theorem employed, a negative integral term with gradient was removed in their deduction. Therefore, the stability criteria acquired by them do not contain diffusion terms; that is to say, the diffusion terms do not have any effect on their deduction and results. Meanwhile, some conditions dependent on the diffusion coefficients were given in [30, 3234] to ensure the global exponential stability and periodicity of RDNNs with Dirichlet boundary conditions based on 2-norm.

To the best of our knowledge, there are few reports about global exponential stability and periodicity of RDNNs with mixed time delays and Dirichlet boundary conditions, which are very important in theories and applications and also are a very challenging problem. In this paper, by employing suitable Lyapunov functionals, we shall apply inequality techniques to establish global exponential stability criteria of the equilibrium and periodic solutions for RDNNs with mixed time delays and Dirichlet boundary conditions. The derived criteria extend and improve previous results in the literature [22, 29].

Throughout this paper, we need the following notations. R n denotes the n-dimensional Euclidean space. We denote

u ( t , x ) u = Ω i = 1 m | u i u i | r d x , φ u ( s , x ) u = sup s 0 [ Ω i = 1 m | φ u i ( s , x ) u i | r d x ]

and

v ( t , x ) v = Ω j = 1 n | v j v j | r d x , φ v ( s , x ) v = sup s 0 [ Ω j = 1 n | φ v j ( s , x ) v j | r d x ] , r 2 .

Let u i = u i (t,x), v j = v j (t,x).

The remainder of this paper is organized as follows. In Section 2, the basic notations, model description and assumptions are introduced. In Sections 3 and 4, criteria are proposed to determine global exponential stability, and periodic solutions are considered for reaction-diffusion recurrent neural networks with mixed time delays, respectively. An illustrative example is given to illustrate the effectiveness of the obtained results in Section 5. We also conclude this paper in Section 6.

2 Model description and preliminaries

In this paper, the RDNNs with mixed time delays are described as follows:

(1)

The RDNNs model given in (1) can be regarded as RDNNs with two layers, where m is the number of neurons in the first layer and n is the number of neurons in the second layer. x= ( x 1 , x 2 , , x l ) T Ω R l , Ω is a compact set with smooth boundary Ω and mesΩ>0 in the space R l ; u= ( u 1 , u 2 , , u m ) T R m , v= ( v 1 , v 2 , , v n ) T R n . u i (t,x) and v j (t,x) represent the state of the i th neuron in the first layer and the j th neuron in the second layer at time t and in the space x, respectively. b j i , b ˜ j i , b ¯ j i , d i j , d ¯ i j and d ˜ i j are known constants denoting the synaptic connection strengths between the neurons in the two layers, respectively; f j , f ˜ j , f ¯ j , g i , g ˜ i and g ¯ i denote the activation functions of the neurons and the signal propagation functions, respectively. I i and J j denote the external inputs on the i th neuron and j th neuron, respectively; p i and q j are differentiable real functions with positive derivatives defining the neuron charging time; τ i j (t) and θ j i (t) represent continuous time-varying delay and discrete delay, respectively; D i k 0 and D j k 0, i=1,2,,m, k=1,2,,l and j=1,2,,n, stand for the transmission diffusion coefficient along the i th neuron and j th neuron, respectively.

System (1) is supplemented with the following boundary conditions and initial values:

u i (t,x)=0, v j (t,x)=0,t0,xΩ,
(2)
u i (s,x)= φ u i (s,x), v j (s,x)= φ v j (s,x),(s,x)(,0]×Ω
(3)

for any i=1,2,,m and j=1,2,,n, where n ¯ is the outer normal vector of Ω, φ= ( φ u φ v ) = ( φ u 1 , , φ u m , φ v 1 , , φ v n ) T C are bounded and continuous, where C= { φ | φ = ( φ u φ v ) , φ : ( ( , 0 ] × R m ( , 0 ] × R n ) R m + n } . It is the Banach space of continuous functions which maps ( ( , 0 ] ( , 0 ] ) into R m + n with the topology of uniform convergence for the norm

φ= ( φ u φ v ) = sup s 0 [ Ω i = 1 m | φ u i | r d x ] + sup s 0 [ Ω j = 1 n | φ v j | r d x ] .

Remark 1 Some famous NN models became a special case of system (1). For example, when D i k =0 and D j k =0 (i=1,2,,m, k=1,2,,l), the special case of model (1) is the model which has been studied in [1315]. When b ˜ j i =0 and d ˜ i j =0, i=1,2,,m, j=1,2,,n, system (1) became NNs with distributed delays and reaction-diffusion terms [18, 22, 29].

Throughout this paper, we assume that the following conditions are made.

(A1) The functions τ i j (t), θ j i (t) are piecewise-continuous of class C 1 on the closure of each continuity subinterval and satisfy

0 τ i j ( t ) τ i j , 0 θ j i ( t ) θ j i , τ ˙ i j ( t ) μ τ < 1 , θ ˙ j i ( t ) μ θ < 1 , τ = max 1 i m , 1 j n { τ i j } , θ = max 1 i m , 1 j n { θ j i }

with some constants τ i j 0, θ j i 0, τ>0, θ>0 for all t0.

(A2) The functions p i () and q j () are piecewise-continuous of class C 1 on the closure of each continuity subinterval and satisfy

a i = inf ζ R p i ( ζ ) > 0 , p i ( 0 ) = 0 , c j = inf ζ R q j ( ζ ) > 0 , q j ( 0 ) = 0 .

(A3) The activation functions and the signal propagation functions are bounded and Lipschitz continuous, i.e., there exist positive constants L j f , L j f ˜ , L j f ¯ , L i g , L i g ˜ and L i g ¯ such that for all η 1 , η 2 R,

| f j ( η 1 ) f j ( η 2 ) | L j f | η 1 η 2 | , | f ˜ j ( η 1 ) f ˜ j ( η 2 ) | L j f ˜ | η 1 η 2 | , | f ¯ j ( η 1 ) f ¯ j ( η 2 ) | L j f ¯ | η 1 η 2 | , | g i ( η 1 ) g i ( η 2 ) | L i g | η 1 η 2 | , | g ˜ i ( η 1 ) g ˜ i ( η 2 ) | L i g ˜ | η 1 η 2 | , | g ¯ i ( η 1 ) g ¯ i ( η 2 ) | L i g ¯ | η 1 η 2 | .

(A4) The delay kernels K j i (s), K ¯ i j (s):[0,)[0,) (i=1,2,,m, j=1,2,,n) are real-valued non-negative continuous functions that satisfy the following conditions:

  1. (i)

    0 + K j i (s)ds=1, 0 + K ¯ j i (s)ds=1;

  2. (ii)

    0 + s K j i (s)ds<, 0 + s K ¯ i j (s)ds<;

  3. (iii)

    There exist a positive μ such that

    0 + s e μ s K j i (s)ds<, 0 + s e μ s K ¯ i j (s)ds<.

Let ( u , v )=( u 1 , u 2 ,, u n , v 1 , v 2 ,, v n ) be the equilibrium point of system (1).

Definition 1 The equilibrium point of system (1) is said to be globally exponentially stable if we can find r2 such that there exist constants α>0 and β1 such that

u ( t , x ) u + v ( t , x ) v β e 2 α t ( φ u ( s , x ) u + φ v ( s , x ) v )
(4)

for all t0.

Remark 2 It is well known that bounded activation functions always guarantee the existence of an equilibrium point for system (1).

Lemma 1 [33]

Let Ω be a cube | x l |< d l (l=1,,m), and let h(x) be a real-valued function belonging to C 1 (Ω) which vanishes on the boundary Ω of Ω, i.e., h(x) | Ω = 0 . Then

Ω h 2 (x)dx d l 2 Ω | h x l | 2 dx.
(5)

3 Global exponential stability

Now we are in a position to investigate the global exponential stability of system (1). By constructing a suitable Lyapunov functional, we arrive at the following conclusion.

Theorem 1 Let (A1)-(A4) be in force. If there exist w i >0 (i=1,2,,n+m), r2, γ i j >0, β j i >0 such that

w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + ( r 1 ) j = 1 n a i r β j i r r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ ) r ) < 0

and

w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r ) < 0 ,
(6)

in which i=1,2,,m, j=1,2,,n, L j f , L j f ˜ , L j f ¯ , L i g , L i g ˜ and L i g ¯ are Lipschitz constants, D i = min 1 k l { D i k / d k 2 }, D j = min 1 k l { D j k / d k 2 }, then the equilibrium point ( u , v ) of system (1) is unique and globally exponentially stable.

Proof If (6) holds, we can always choose a positive number δ>0 (may be very small) such that

w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + ( r 1 ) j = 1 n a i r β j i r r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ i ) r ) + δ < 0

and

w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r ) + δ < 0 ,
(7)

where i=1,2,,m, j=1,2,,n.

Let us consider the functions

F i ( x i ) = w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + ( r 1 ) j = 1 n a i r β j i r r 1 0 + k j i ( s ) d s + 2 x i n a i r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ i ) r 0 + e 2 x i s k ¯ i j ( s ) d s )

and

G j ( y j ) = w m + j [ r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 0 + k ¯ i j ( s ) d s + 2 y j m c j r 1 ] + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r 0 + e 2 y j s k j i ( s ) d s ) ,
(8)

where x i , y j [0,+), i=1,2,,m, j=1,2,,n.

From (8) and (A4), we derive F i (0)<δ<0, G j (0)<δ<0; F i ( x i ) and G j ( y j ) are continuous for x i , y j [0,+). Moreover, F i ( x i )+ as x i + and G j ( y j )+ as y j +. Thus there exist constants ε i , σ j [0,+) such that

F i ( ε i ) = w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r F i ( ε i ) = + ( r 1 ) j = 1 n a i r β j i r r 1 0 + k j i ( s ) d s + 2 ε i n a i r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r F i ( ε i ) = + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ i ) r 0 + e 2 ε i s k ¯ i j ( s ) d s ) = 0

and

G j ( σ j ) = w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 0 + k ¯ i j ( s ) d s + 2 σ j m c j r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r 0 + e 2 σ j s k j i ( s ) d s ) = 0 ,
(9)

where i=1,2,,m, j=1,2,,n.

By using α= min 1 i m , 1 j n { ε i , σ j }, obviously, we get

F i ( α ) = w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r F i ( α ) = + ( r 1 ) j = 1 n a i r β j i r r 1 0 + k j i ( s ) d s + 2 α n a i r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r F i ( α ) = + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ i ) r 0 + e 2 α s k ¯ i j ( s ) d s ) 0

and

G j ( α ) = w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 0 + k ¯ i j ( s ) d s + 2 α m c j r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r 0 + e 2 α s k j i ( s ) d s ) 0 ,
(10)

where i=1,2,,m, j=1,2,,n.

Suppose (u,v)= ( u 1 , u 2 , , u n , v 1 , v 2 , , v n ) T is any solution of model (1). Rewrite model (1) as

(11)
(12)

Multiplying (11) by u i u i and integrating over Ω yield

1 2 d d t Ω ( u i u i ) 2 d x = Ω k = 1 l ( u i u i ) x k ( D i k ( u i u i ) x k ) d x p i ( ξ i ) Ω ( u i u i ) 2 d x + Ω j = 1 n ( b j i ( u i u i ) ( f j ( v j ) f j ( v j ) ) ) d x + j = 1 n Ω ( b ˜ j i ( u i u i ) ( f ˜ j ( v j ( t θ j i ( t ) , x ) ) f ˜ j ( v j ) ) ) d x + j = 1 n Ω [ ( b ¯ j i ( u i u i ) t k j i ( t s ) ( f ¯ j ( v j ( s , x ) ) f ¯ j ( v j ) ) ) d s ] d x .
(13)

According to Green’s formula and the Dirichlet boundary condition, we get

Ω k = 1 l ( u i u i ) x k ( D i k ( u i u i ) x k ) dx= k = 1 l Ω D i k ( ( u i u i ) x k ) 2 dx.
(14)

Moreover from Lemma 1, we have

k = 1 l Ω D i k ( ( u i u i ) x k ) 2 dx Ω k = 1 l D i k d k 2 ( u i u i ) 2 dx D i l u i u i 2 2 .
(15)

From (11)-(15), (A2), (A3) and the Holder integral inequality, we obtain that

d d t Ω | u i u i | 2 d x 2 D i l Ω | u i u i | 2 d x 2 a i Ω | u i u i | 2 d x + 2 Ω j = 1 n ( | b j i | | u i u i | L j f | v j v j | ) d x + 2 j = 1 n Ω ( | b ˜ j i | | u i u i | | f ˜ ( v j ( t θ j i ( t ) , x ) ) f ˜ ( v j ) | ) d x + 2 j = 1 n Ω [ ( | b ¯ j i | t k j i ( t s ) | u i u i | | f ¯ j ( v j ( s , x ) ) f ¯ j ( v j ) | ) d s ] d x .
(16)

Multiplying both sides of (12) by v j v j , similarly, we also have

d d t Ω | v j v j | 2 d x 2 D j l Ω | v j v j | 2 d x 2 c j Ω | v j v j | 2 d x + 2 Ω i = 1 m ( | d i j | L i g | u i u i | | v j v j | ) d x + 2 i = 1 m Ω ( | d ˜ i j | | g ˜ ( u i ( t τ i j ( t ) , x ) ) g ˜ ( u i ) | | v j v j | ) d x + 2 i = 1 m Ω [ | d ¯ i j | t k ¯ i j ( t s ) | g ¯ i ( u i ( s , x ) ) g ¯ i ( u i ) | | v j v j | d s ] d x .
(17)

Choose a Lyapunov functional as follows:

V ( t ) = Ω i = 1 m w i [ n a i r 1 | u i u i | r e 2 α t + j = 1 n | b ˜ j i | r n r e θ 1 μ θ t θ j i ( t ) t e 2 α ξ | f ˜ j ( v j ( ξ , x ) ) f ˜ j ( v j ) | r d ξ + j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) t s t e 2 α ( s + ξ ) | f ¯ j ( v j ( ξ , x ) ) f ¯ j ( v j ) | r d ξ d s ] d x + Ω j = 1 n w m + j [ m c j r 1 | v j v j | r e 2 α t + i = 1 m | d ˜ i j | r m r e τ 1 μ τ t τ i j ( t ) t e 2 α ξ | g ˜ i ( u i ( ξ , x ) ) g ˜ i ( u i ) | r d ξ + i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) t s t e 2 α ( s + ξ ) | g ¯ i ( u i ( ξ , x ) ) g ¯ i ( u i ) | r d ξ d s ] d x .

Its upper Dini-derivative along the solution to system (1) can be calculated as follows:

D + V ( t ) Ω i = 1 m w i [ r n a i r 1 | u i u i | r 1 | u i u i | t e 2 α t + 2 α e 2 α t n a i r 1 | u i u i | r + e 2 α t j = 1 n | b ˜ j i | r n r e θ 1 μ θ | f ˜ j ( v j ( t , x ) ) f ˜ j ( v j ) | r j = 1 n | b ˜ j i | r n r e θ 1 μ θ ( 1 θ ˙ j i ( t ) ) e 2 α ( t θ j i ( t ) ) | f ˜ j ( v j ( t θ j i ( t ) , x ) ) f ˜ j ( v j ) | r + e 2 α t j = 1 n | b ¯ j i | r n r β j i r 0 + e 2 α s k j i ( s ) | f ¯ j ( v j ( t , x ) ) f ¯ j ( v j ) | r d s e 2 α t j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) | f ¯ j ( v j ( t s , x ) ) f ¯ j ( v j ) | r d s ] d x + Ω j = 1 n w m + j [ r m c j r 1 | v j v j | r 1 | v j v j | t e 2 α t + 2 α e 2 α t m c j r 1 | v j v j | r + e 2 α t i = 1 m | d ˜ i j | r m r e τ 1 μ τ | g ˜ i ( u i ( t , x ) ) g ˜ i ( u i ) | r i = 1 m | d ˜ i j | r m r e τ 1 μ τ e 2 α ( t τ i j ( t ) ) ( 1 τ ˙ i j ) | g ˜ i ( u i ( t τ i j ( t ) , x ) ) g ˜ i ( u i ) | r + e 2 α t i = 1 m | d ¯ i j | r m r γ i j r 0 + e 2 α s k ¯ i j ( s ) | g ¯ i ( u i ( t , x ) ) g ¯ i ( u i ) | r d s e 2 α t i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) | g ¯ i ( u i ( t s , x ) ) g ¯ i ( u i ) | r d s ] d x Ω i = 1 m w i [ r n a i r 1 | u i u i | r 2 e 2 α t ( D i l | u i u i | 2 a i | u i u i | 2 + j = 1 n ( | b j i | | u i u i | L j f | v j v j | ) + j = 1 n ( | b ˜ j i | | u i u i | | f ˜ ( v j ( t θ j i ( t ) , x ) ) f ˜ ( v j ) | ) + j = 1 n ( | b ¯ j i | t k j i ( t s ) | u i u i | | f ¯ j ( v j ( s , x ) ) f ¯ j ( v j ) | ) d s ) + 2 α e 2 α t n a i r 1 | u i u i | r + e 2 α t j = 1 n | b ˜ j i | r n r e θ 1 μ θ | f ˜ j ( v j ( t , x ) ) f ˜ j ( v j ) | r e 2 α t j = 1 n | b ˜ j i | r n r | f ˜ j ( v j ( t θ j i , x ) ) f ˜ j ( v j ) | r + e 2 α t j = 1 n | b ¯ j i | r n r β j i r 0 + e 2 α s k j i ( s ) | f ¯ j ( v j ( t , x ) ) f ¯ j ( v j ) | r d s e 2 α t j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) | f ¯ j ( v j ( t s , x ) ) f ¯ j ( v j ) | r d s ] d x + Ω j = 1 n w m + j [ r m c j r 1 | v j v j | r 2 e 2 α t ( D j l | v j v j | 2 c j | v j v j | 2 + i = 1 m ( | d i j | L i g | u i u i | | v j v j | ) + i = 1 m ( | d ˜ i j | | g ˜ ( u i ( t τ i j , x ) ) g ˜ ( u i ) | | v j v j | ) + i = 1 m ( | d ¯ i j | t k ¯ i j ( t s ) | g ¯ i ( u i ( s , x ) ) g ¯ i ( u i ) | | v j v j | ) d s ) + 2 α e 2 α t m c j r 1 | v j v j | r + e 2 α t i = 1 m | d ˜ i j | r m r e τ 1 μ τ | g ˜ i ( u i ( t , x ) ) g ˜ i ( u i ) | r e 2 α t i = 1 m | d ˜ i j | r m r | g ˜ i ( u i ( t τ i j , x ) ) g ˜ i ( u i ) | r + e 2 α t i = 1 m | d ¯ i j | r m r γ i j r 0 + e 2 α s k ¯ j i ( s ) | g ¯ i ( u i ( t , x ) ) g ¯ i ( u i ) | r d s e 2 α t i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) | g ¯ i ( u i ( t s , x ) ) g ¯ i ( u i ) | r d s ] d x .
(18)

From (18) and the Young inequality, we can conclude

D + V ( t ) Ω e 2 α t i = 1 m w i [ ( r n a i r 1 D i l | u i u i | r r n a i r | u i u i | r + ( r 1 ) j = 1 n a i r | u i u i | r + j = 1 n ( n r | b j i | r ( L j f ) r | v j v j | r ) + ( r 1 ) j = 1 n a i r | u i u i | r + j = 1 n | b ˜ j i | r n r | f ˜ j ( v j ( t θ j i ( t ) , x ) ) f ˜ j ( v j ) | r + ( r 1 ) j = 1 n ( a i r β j i r r 1 t k j i ( t s ) | u i u i | r d s ) + j = 1 n ( | b ¯ j i | r n r β j i r t k j i ( t s ) | f ¯ j ( v j ( s , x ) ) f ¯ j ( v j ) | r d s ) ) + 2 α n a i r 1 | u i u i | r + j = 1 n | b ˜ j i | r n r e θ 1 μ θ | f ˜ j ( v j ( t , x ) ) f ˜ j ( v j ) | r j = 1 n | b ˜ j i | r n r | f ˜ j ( v j ( t θ j i ( t ) , x ) ) f ˜ j ( v j ) | r + j = 1 n | b ¯ j i | r n r β j i r 0 + e 2 α s k j i ( s ) | f ¯ j ( v j ( t , x ) ) f ¯ j ( v j ) | r d s j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) | f ¯ j ( v j ( t s , x ) ) f ¯ j ( v j ) | r d s ] d x + Ω e 2 α t j = 1 n w m + j [ ( r m c j r 1 D j l | v j v j | r r m c j r | v j v j | r + ( r 1 ) i = 1 m c j r | v j v j | r + i = 1 m ( | d i j | r m r ( L i g ) r | u i u i | r ) + ( r 1 ) i = 1 m ( c j r | v j v j | r ) + i = 1 m | d ˜ i j | r m r | g ˜ i ( u i ( t τ i j ( t ) , x ) ) g ˜ i ( u i ) | r + ( r 1 ) i = 1 m ( c j r γ i j r r 1 t k ¯ i j ( t s ) | v j v j | r d s ) + i = 1 m ( | d ¯ i j | r m r γ i j r t k ¯ i j ( t s ) | g ¯ i ( u i ( s , x ) ) g ¯ i ( u i ) | r d s ) ) + 2 α m c j r 1 | v j v j | r + i = 1 m | d ˜ i j | r m r e τ 1 μ τ | g ˜ i ( u i ( t , x ) ) g ˜ i ( u i ) | r i = 1 m | d ˜ i j | r m r | g ˜ i ( u i ( t τ i j ( t ) , x ) ) g ˜ i ( u i ) | r + i = 1 m | d ¯ i j | r m r γ i j r 0 + e 2 α s k ¯ i j ( s ) | g ¯ i ( u i ( t , x ) ) g ¯ i ( u i ) | r d s i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) | g ¯ i ( u i ( t s , x ) ) g ¯ i ( u i ) | r d s ] d x Ω e 2 α t i = 1 m [ w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + 2 α n a i r 1 + ( r 1 ) j = 1 n ( a i r β j i r r 1 t k j i ( t s ) d s ) ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r 0 + e 2 α s k ¯ i j ( s ) ( L i g ¯ i ) r d s ) ] | u i u i | r d x + Ω e 2 α t j = 1 n [ w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m ( c j r γ i j r r 1 t k ¯ i j ( t s ) d s ) + 2 α m c j r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r 0 + e 2 α s k j i ( s ) d s ) ] | v j v j | r d x .
(19)

From (6), we can conclude

D + V(t)0and soV(t)V(0),t0.
(20)

Since

V ( 0 ) = Ω i = 1 m w i [ n a i r 1 | u i ( 0 , x ) u i | r + j = 1 n | b ˜ j i | r n r e θ 1 μ θ θ j i ( t ) 0 | f ˜ j ( v j ( ξ , x ) ) f ˜ j ( v j ) | r d ξ + j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) s 0 e 2 α ( s + ξ ) | f ¯ j ( v j ( ξ , x ) ) f ¯ j ( v j ) | r d ξ d s ] d x + Ω j = 1 n w m + j [ m c j r 1 | v j ( 0 , x ) v j | r + i = 1 m | d ˜ i j | r m r e τ 1 μ τ τ i j ( t ) 0 | g ˜ i ( u i ( ξ , x ) ) g ˜ i ( u i ) | r d ξ + i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) s 0 e 2 α ( s + ξ ) | g ¯ i ( u i ( ξ , x ) ) g ¯ i ( u i ) | r d ξ d s ] d x Ω i = 1 m max 1 i m { w i } [ n a i r 1 | u i ( 0 , x ) u i | r + j = 1 n | b ˜ j i | r n r e θ 1 μ θ ( L j f ˜ ) r θ j i 0 | v j ( ξ , x ) v j | r d ξ + j = 1 n | b ¯ j i | r n r ( L j f ¯ ) r β j i r 0 + k j i ( s ) s 0 e 2 α ( s + ξ ) | v j ( ξ , x ) v j | r d ξ d s ] d x + Ω j = 1 n max 1 j n { w m + j } [ m c j r 1 | v j ( 0 , x ) v j | r + i = 1 m | d ˜ i j | r ( L i g ˜ ) r m r e τ 1 μ τ τ i j 0 | u i ( ξ , x ) u i | r d ξ + j = 1 n | d ¯ i j | r m r ( L i g ¯ ) r γ i j r 0 + k ¯ j i ( s ) s 0 e 2 α ( s + ξ ) | u i ( ξ , x ) u i | r d ξ d s ] d x { max 1 i m { w i } + max 1 j n { w m + j } max 1 j n [ i = 1 m | d ¯ i j | r m r ( L i g ¯ ) r γ i j r 0 + k ¯ i j ( s ) s e 2 α s d s ] + max 1 j n { w m + j } max 1 j n [ i = 1 m | d ˜ i j | r ( L i g ˜ ) r m r e τ τ 1 μ τ ] } φ u ( s , x ) u r + { max 1 j n { w m + j } + max 1 i m { w i } max 1 i m [ j = 1 n | b ¯ j i | r n r ( L j f ¯ ) r β j i r 0 + s e 2 α s k j i ( s ) d s ] + max 1 i m { w i } max 1 i m [ j = 1 n | b ˜ j i | r n r ( L j g ˜ ) r e θ θ 1 μ θ ] } φ v ( s , x ) v r .
(21)

Noting that

e 2 α t ( min 1 i m + n w i ) ( u ( t , x ) u + v ( t , x ) v ) V(t),t0.
(22)

Let

β = max { max 1 i m { w i } + max 1 j n { w m + j } max 1 j n [ i = 1 m | d ¯ i j | r m r ( L i g ¯ ) r γ i j r 0 + k ¯ i j ( s ) s e 2 α s d s ] β = + max 1 j n { w m + j } max 1 j n [ i = 1 m | d ˜ i j | r ( L i g ˜ ) r m r e τ τ 1 μ τ ] , max 1 j n { w m + j } + max 1 i m { w i } max 1 i m [ j = 1 n | b ¯ j i | r n r ( L j f ¯ ) r β j i r 0 + s e 2 α s k j i ( s ) d s ] + max 1 i m { w i } max 1 i m [ j = 1 n | b ˜ j i | r n r ( L j g ˜ ) r e θ θ 1 μ θ ] } / min 1 i m + n { w i } .

Clearly, β1.

It follows that

u ( t , x ) u + v ( t , x ) v β e 2 α t ( φ u ( s , x ) u + φ v ( s , x ) v ) ,

for any t0, where β1 is a constant. This implies that the solution of (1) is globally exponentially stable. This completes the proof of Theorem 1. □

Remark 3 In this paper, the derived sufficient condition includes diffusion terms. Unfortunately, in the proof in the previous papers [2124], a negative integral term with gradient is left out in their deduction. This leads to the fact that those criteria are irrelevant to the diffusion term. Obviously, Lyapunov functional to construct is more general and our results expand the model in [22, 29].

When b ˜ j i =0 and d ˜ i j =0 (i=1,2,,m, j=1,2,,n), system (1) becomes the following BAM NNs with distributed delays and reaction-diffusion terms:

u i t = k = 1 l x k ( D i k u i x k ) p i ( u i ( t , x ) ) + j = 1 n ( b j i f j ( v j ( t , x ) ) ) u i t = + j = 1 n b ¯ j i t k j i ( t s ) f ¯ j ( v j ( s , x ) ) d s + I i ( t ) , v j t = k = 1 l x k ( D j k v j x k ) q j ( v j ( t , x ) ) + i = 1 m ( d i j g i ( u i ( t , x ) ) ) v j t = + i = 1 m d ¯ i j t k ¯ i j ( t s ) g ¯ i ( u i ( s , x ) ) d s + J j ( t ) .
(23)

For (23), we get the following result.

Corollary 1 Let (A1)-(A4) be in force. If there exist w i >0 (i=1,2,,n+m), r2, γ i j >0, β j i >0 such that

w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + ( r 1 ) j = 1 n a i r β j i r r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ¯ i j | r γ i j r ( L i g ¯ ) r ) < 0

and

w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r ) < 0 ,
(24)

where i=1,2,,m, j=1,2,,n, L j f , L j f ˜ , L j f ¯ , L i g , L i g ˜ and L i g ¯ are Lipschitz constants. Then the equilibrium point ( u , v ) of system (1) is unique and globally exponentially stable.

4 Periodic solutions

In this section, we consider the stability criterion for periodic oscillatory solutions of system (1), in which external input I i : R + R, i=1,2,,m, and J j : R + R, j=1,2,,n, are continuously periodic functions with period ω, that is,

I i (t+ω)= I i (t), J j (t+ω)= J j (t),i=1,2,,m,j=1,2,,n.

By constructing a Poincaré mapping, the existence of a unique ω-periodic solution and its stability are readily established.

Theorem 2 Let (A1)-(A4) be in force. There exists only one ω-periodic solution of system (1), and all other solutions converge exponentially to it as t+ if there exist constants w i >0 (i=1,2,,n+m), r2, γ i j >0, β j i >0 (i=1,2,,m, j=1,2,,n) such that

w i ( r n a i r 1 D i l r n a i r + 2 ( r 1 ) j = 1 n a i r + ( r 1 ) j = 1 n a i r β j i r r 1 ) + j = 1 n w m + j m r ( | d i j | r ( L i g ) r + | d ˜ i j | r e τ 1 μ τ ( L i g ˜ ) r + | d ¯ i j | r γ i j r ( L i g ¯ ) r ) < 0

and

w m + j ( r m c j r 1 D j l r m c j r + 2 ( r 1 ) i = 1 m c j r + ( r 1 ) i = 1 m c j r γ i j r r 1 ) + i = 1 m w i n r ( | b j i | r ( L j f ) r + | b ˜ j i | r e θ 1 μ θ ( L j f ˜ ) r + | b ¯ j i | r β j i r ( L j f ¯ ) r ) < 0 ,
(25)

where i=1,2,,m and j=1,2,,n, L j f , L j f ˜ , L j f ¯ , L i g , L i g ˜ and L i g ¯ are Lipschitz constants in (A3).

Proof For any ( φ u φ v ) , ( ψ u ψ v ) C, we denote the solutions of system (1) through ( ( 0 0 ) ) , ( ( φ u φ v ) ) and ( ( 0 0 ) ) , ( ( ψ u ψ v ) ) as

u(t, φ u ,x)= ( u 1 ( t , φ u , x ) , , u m ( t , φ u , x ) ) T ,v(t, φ v ,x)= ( v 1 ( t , φ v , x ) , , v n ( t , φ v , x ) ) T

and

u(t, ψ u ,x)= ( u 1 ( t , ψ u , x ) , , u m ( t , ψ u , x ) ) T ,v(t, ψ v ,x)= ( v 1 ( t , ψ v , x ) , , v n ( t , ψ v , x ) ) T ,

respectively. Define

u t ( φ u , x ) = u ( t + θ , φ u , x ) , θ ( , 0 ] , t 0 , v t ( φ v , x ) = v ( t + θ , φ v , x ) , θ ( , 0 ] , t 0 .

Clearly, for any t0, ( u t ( φ u ) v t ( φ v ) ) C. Now, we define

y i = u i (t, φ u ,x) u i (t, ψ u ,x), z j = v j (t, φ v ,x) v j (t, ψ v ,x).

Thus, we can obtain from system (1) that

y i t = k = 1 l x k ( D i k y i x k ) ( p i ( u i ( t , φ u , x ) ) p i ( u i ( t , ψ u , x ) ) ) y i t = + j = 1 n ( b j i ( f j ( v j ( t , φ v , x ) ) f j ( v j ( t , ψ v , x ) ) ) ) y i t = + j = 1 n ( b ˜ j i ( f ˜ j ( v j ( t θ j i ( t ) , φ v , x ) ) f ˜ j ( v j ( t θ j i ( t ) , ψ v , x ) ) ) ) y i t = + j = 1 n ( b ¯ j i t k j i ( t s ) ( f ¯ j ( v j ( s , φ v , x ) ) f ¯ j ( v j ( s , ψ v , x ) ) ) ) d s , z j t = k = 1 l x k ( D j k z j x k ) ( q j ( v j ( t , φ v , x ) ) q j ( v j ( t , ψ v , x ) ) ) z j t = + i = 1 m ( d i j ( g i ( u i ( t , φ u , x ) ) g i ( u i ( t , ψ u , x ) ) ) ) z j t = + i = 1 m ( d ˜ i j ( g ˜ i ( u i ( t τ i j ( t ) , φ u , x ) ) g ˜ i ( u i ( t τ i j ( t ) , ψ u , x ) ) ) ) z j t = + i = 1 m ( d ¯ i j t k ¯ i j ( t s ) ( g ¯ i ( u i ( s , φ u , x ) ) g ¯ i ( u i ( s , ψ u , x ) ) ) d s ) .

We consider the following Lyapunov functional:

V ( t ) = Ω i = 1 m w i [ n a i r 1 | y i | r e 2 α t + j = 1 n | b ˜ j i | r n r ( 1 μ θ ) t θ j i ( t ) t | f ˜ j ( v j ( ξ , φ v , x ) ) f ˜ j ( v j ( ξ , ψ v , x ) ) | r d ξ + j = 1 n | b ¯ j i | r n r β j i r 0 + k j i ( s ) t s t e 2 α ( s + ξ ) × | f ¯ j ( v j ( ξ , φ v , x ) ) f ¯ j ( v j ( ξ , ψ v , x ) ) | r d ξ d s ] d x + Ω j = 1 n w m + j [ m c j r 1 | z j | r e 2 α t + i = 1 m | d ˜ i j | r m r ( 1 μ τ ) t τ i j ( t ) t | g ˜ i ( u i ( ξ , φ u , x ) ) g ˜ i ( u i ( ξ , ψ u , x ) ) | r d ξ d s + i = 1 m | d ¯ i j | r m r γ i j r 0 + k ¯ i j ( s ) t s t e 2 α ( s + ξ ) × | g ¯ i ( u i ( ξ , φ u , x ) ) g ¯ i ( u i ( ξ , ψ u , x ) ) | r d ξ d s ] d x .

By a minor modification of the proof of Theorem 1, we can easily get

u ( t , φ u , x ) u ( t , ψ u , x ) + v ( t , φ v , x ) v ( t , ψ v , x ) β e 2 α t ( φ u ψ u + φ v ψ v )
(26)

for t0,in which β1 is a constant. Now, we can choose a positive integer N such that

β e α N ω 1 4 ,β e α N ω 1 4 .
(27)

Defining a Poincaré mapping P:CC by

P ( φ u φ v ) = ( u ω ( φ u ) v ω ( φ v ) ) ,
(28)

due to the periodicity of system, we have

P N ( φ u φ v ) = ( u N ω ( φ u ) v N ω ( φ v ) ) .
(29)

Let t=Nω, then from (26)-(29) we can derive that

| P N ( φ u φ v ) P N ( ψ u ψ v ) | 1 2 | ( φ u φ v ) ( ψ u ψ v ) |,

which shows that P N is a contraction mapping. Therefore, there exists a unique fixed point ( φ u φ v ) C, namely, P N ( φ u φ v ) = ( φ u φ v ) .

Since P N ( P ( φ u φ v ) ) =P ( P N ( φ u φ v ) ) =P ( φ u φ v ) , then P ( φ u φ v ) is also a fixed point of P N . Because of the uniqueness of a fixed point of P N , then P ( φ u φ v ) = ( φ u φ v ) .

Let (u(t, φ u ,x),v(t, φ v ,x)) be the solution of system (1) through ( ( 0 0 ) , ( φ u φ v ) ) , then (u(t+ω, φ u ,x),v(t+ω, φ v ,x)) is also a solution of system (1). Clearly,

( u t + ω ( φ u , x ) v t + ω ( φ v , x ) ) = ( u t ( u ω ( φ u ) ) v t ( v ω ( φ v ) ) ) = ( u t ( φ u , x ) v t ( φ v , x ) )

for t0. Hence ( u ( t + ω , φ u , x ) , v ( t + ω , φ v , x ) ) T = ( u ( t , φ u , x ) , v ( t , φ v , x ) ) T for t0.

This shows that ( u ( t , φ u , x ) , v ( t , φ v , x ) ) T is exactly one ω-periodic solution of system (1), and it is easy to see that all other solutions of system (1) converge exponentially to it as t+. The proof is completed. □

5 Illustration example

In this section, a numerical example is given to illustrate the effectiveness of the obtained results.

Example 1 Consider the following system on Ω={ ( x 1 , x 2 ) T |0< x k < 0.2 π,k=1,2} R 2 :

(30)

where k j i (t)= k ¯ i j (t)=t e t , i,j,l=1,2. f 1 (η)= f 2 (η)= f ˜ 1 (η)= f ˜ 2 (η)= f ¯ 1 (η)= f ¯ 2 (η)= g 1 (η)= g 2 (η)= g ˜ 1 (η)= g ˜ 2 (η)= g ¯ 1 (η)= g ¯ 2 (η)=tanh(η), n=m=l=2, λ 1 =2.5, θ j i (t)= τ i j (t)=0.020.01sin(2πt), L j f = L j f ˜ = L j f ¯ = L i g = L i g ˜ = L i g ¯ =1, i,j=1,2. p i ( u i (t,x))= u i (t,x), q j ( v j (t,x))=2 v j (t,x), D 1 = D 2 =1, D 1 = D 2 =2, a 1 = a 2 =1, c 1 = c 2 =2, r=2, μ τ = μ θ =0.2, d 11 =0.5, d 12 =1, d 21 =0.5, d 22 =0.2, d ˜ 11 =0.1, d ˜ 12 =0.2, d ˜ 21 =0.3, d ˜ 22 =0.5, d ¯ 11 =0.2, d ¯ 12 =0.6, d ¯ 21 =0.5, d ¯ 22 =0.8, b 11 =0.3, b 12 =0.6, b 21 =0.5, b 22 =0.8, b ˜ 11 =1, b ˜ 12 =0.5, b ˜ 21 =0.3, b ˜ 22 =0.3, b ¯ 11 =0.1, b ¯ 12 =0.2, b ¯ 21 =0.5, b ¯ 22 =0.4. By simple calculation with w 1 = w 2 = w 3 = w 4 =1, β 11 = β 12 = β 21 = β 22 =1 and γ 11 = γ 12 = γ 21 = γ 22 =1, we have

r n a 1 r 1 D 1 λ 1 r n a 1 r + 2 ( r 1 ) j = 1 2 a 1 r + ( r 1 ) j = 1 2 a 1 r β j 1 r r 1 + j = 1 2 m r ( | d 1 j | r ( L 1 g ) r + | d ˜ 1 j | r 1 1 μ τ ( L 1 g ˜ ) r + | d ¯ 1 j | r γ 1 j r ( L 1 g ¯ ) r ) = 1.15 < 0 ,
(31)
r n a 2 r 1 D 2 λ 1 r n a 2 r + 2 ( r 1 ) j = 1 2 a 2 r + ( r 1 ) j = 1 2 a 2 r β j 2 r r 1 + j = 1 2 m r ( | d 2 j | r ( L 2 g ) r + | d ˜ 2 j | r 1 1 μ τ ( L 2 g ˜ ) r + | d ¯ 2 j | r γ 2 j r ( L 2 g ¯ ) r ) = 2.38 < 0 ,
(32)
r m c 1 r 1 D 1 λ 1 r m c 1 r + 2 ( r 1 ) i = 1 2 c 1 r + ( r 1 ) i = 1 2 c 1 r γ i 1 r r 1 + i = 1 2 n r ( | b 1 i | r ( L 1 f ) r + | b ˜ 1 i | r 1 1 μ θ ( L 1 f ˜ ) r + | b ¯ 1 i | r β 1 i r ( L 1 f ¯ ) r ) = 19.15 < 0
(33)

and

r m c 2 r 1 D 2 λ 1 r m c 2 r + 2 ( r 1 ) i = 1 2 c 2 r + ( r 1 ) i = 1 2 c 2 r γ i 2 r r 1 + i = 1 2 n r ( | b 2 i | r ( L 2 f ) r + | b ˜ 2 i | r 1 1 μ θ ( L 2 f ˜ ) r + | b ¯ 2 i | r β 2 i r ( L 2 f ¯ ) r ) = 20.98 < 0 ,
(34)

that is, (6) holds.

The simulation results are shown in Figures 1-8. When x 2 =0.5018, the states surfaces of u( x 1 ,0.5018,t) are shown in Figures 1-2, while x 1 =0.5018, the states surfaces of u(0.5018, x 2 ,t) are shown in Figures 3-4. When x 2 =0.5018, the states surfaces of v( x 1 ,0.5018,t) are shown in Figures 5-6, while x 1 =0.5018, the states surfaces of v(0.5018, x 2 ,t) are shown in Figures 7-8, which illustrates that the system states in (30) converge to equilibrium solution. Therefore, it follows from Theorem 1 and the simulation study that (30) has one unique equilibrium solution which is globally exponentially stable.

Figure 1
figure 1

The surface of u 1 ( x 1 ,0.5018,t) when x 2 =0.5018 .

Figure 2
figure 2

The surface of u 2 ( x 1 ,0.5018,t) when x 2 =0.5018 .

Figure 3
figure 3

The surface of u 1 (0.5018, x 2 ,t) when x 1 =0.5018 .

Figure 4
figure 4

The surface of u 2 (0.5018, x 2 ,t) when x 1 =0.5018 .

Figure 5
figure 5

The surface of v 1 ( x 1 ,0.5018,t) when x 2 =0.5018 .

Figure 6
figure 6

The surface of v 2 ( x 1 ,0.5018,t) when x 2 =0.5018 .

Figure 7
figure 7

The surface of v 1 (0.5018, x 2 ,t) when x 1 =0.5018 .

Figure 8
figure 8

The surface of v 2 (0.5018, x 2 ,t) when x 1 =0.5018 .

Remark 4 Since a 1 + 1 2 j = 1 2 ( L j f ) 2 (| b j 1 | 2 +| b ¯ j 1 | 2 )=0.25>0, the conditions of Corollary 3.2 in [22] and r a 1 +(r1) j = 1 2 L j f (| b j 1 |+| b ¯ j 1 |)+ L 1 g j = 1 2 (| d 1 j |+| d ¯ 1 j |)=3.3>0, under the conditions of Example 1, the conditions of Theorem 1 in [29] are not satisfied. However, by (31)-(34) and Theorem 1, we can derive that (30) has one unique equilibrium solution which is globally exponentially stable.

6 Conclusions

In this paper, by employing suitable Lyapunov functionals, Young’s inequality and Hölder’s inequality techniques, global exponential stability criteria of the equilibrium point and periodic solutions for RDNNs with mixed time delays and Dirichlet boundary conditions have been derived, respectively. The derived criteria contain and extend some previous NNs in the literature. Hence, our results have an important significance in design as well as in applications of periodic oscillatory NNs with mixed time delays. An example has been given to show the effectiveness of the obtained results.

References

  1. Kwon OM, Park JH, Lee SM, Cha EJ: A new augmented Lyapunov-Krasovskii functional approach to exponential passivity for neural networks with time-varying delays. Appl. Math. Comput. 2011, 217(24):10231-10238. 10.1016/j.amc.2011.05.021

    Article  MathSciNet  Google Scholar 

  2. Ji DH, Koo JH, Won SC, Lee SM, Park JH: Passivity-based control for Hopfield neural networks using convex representation. Appl. Math. Comput. 2011, 217(13):6168-6175. 10.1016/j.amc.2010.12.100

    Article  MathSciNet  Google Scholar 

  3. Lee SM, Kwon OM, Park JH: A novel delay-dependent criterion for delayed neural networks of neutral type. Phys. Lett. A 2010, 374(17-18):1843-1848. 10.1016/j.physleta.2010.02.043

    Article  MathSciNet  Google Scholar 

  4. Cao J, Wang J: Global exponential stability and periodicity of recurrent neural networks with time delays. IEEE Trans. Circuits Syst. I, Regul. Pap. 2005, 52(5):920-931.

    Article  MathSciNet  Google Scholar 

  5. Huang C, Cao J: Convergence dynamics of stochastic Cohen-Grossberg neural networks with unbounded distributed delays. IEEE Trans. Neural Netw. 2011, 22(4):561-572.

    Article  Google Scholar 

  6. Ensari T, Arik S: Global stability of a class of neural networks with time varying delay. IEEE Trans. Circuits Syst. II, Express Briefs 2005, 52(3):126-130.

    Article  Google Scholar 

  7. Rakkiyappan R, Balasubramaniam P: Delay-dependent asymptotic stability for stochastic delayed recurrent neural networks with time varying delays. Appl. Math. Comput. 2008, 198(2):526-533. 10.1016/j.amc.2007.08.053

    Article  MathSciNet  Google Scholar 

  8. Kosko B: Bi-directional associative memories. IEEE Trans. Syst. Man Cybern. 1988, 18(1):49-60. 10.1109/21.87054

    Article  MathSciNet  Google Scholar 

  9. Park JH, Kwon OM: Delay-dependent stability criterion for bidirectional associative memory neural networks with interval time-varying delays. Mod. Phys. Lett. B 2009, 23(1):35-46. 10.1142/S0217984909017807

    Article  Google Scholar 

  10. Park JH, Park CH, Kwon OM, Lee SM: New stability criterion for bidirectional associative memory neural networks of neutral-type. Appl. Math. Comput. 2008, 199(2):716-722. 10.1016/j.amc.2007.10.032

    Article  MathSciNet  Google Scholar 

  11. Park JH, Kwon OM: On improved delay-dependent criterion for global stability of bidirectional associative memory neural networks with time-varying delays. Appl. Math. Comput. 2008, 199(2):435-446. 10.1016/j.amc.2007.10.001

    Article  MathSciNet  Google Scholar 

  12. Cao J, Wang L: Exponential stability and periodic oscillatory solution in BAM networks with delays. IEEE Trans. Neural Netw. 2002, 13(2):457-463. 10.1109/72.991431

    Article  Google Scholar 

  13. Park J, Lee SM, Kwon OM: On exponential stability of bidirectional associative memory neural networks with time-varying delays. Chaos Solitons Fractals 2009, 39(3):1083-1091. 10.1016/j.chaos.2007.05.003

    Article  MathSciNet  Google Scholar 

  14. Wu R: Exponential convergence of BAM neural networks with time-varying coefficients and distributed delays. Nonlinear Anal., Real World Appl. 2010, 11(1):562-573. 10.1016/j.nonrwa.2009.02.003

    Article  MathSciNet  Google Scholar 

  15. Liu X, Martin R, Wu M: Global exponential stability of bidirectional associative memory neural networks with time delays. IEEE Trans. Neural Netw. 2008, 19(2):397-407.

    Google Scholar 

  16. Zhang W, Li J: Global exponential synchronization of delayed BAM neural networks with reaction-diffusion terms and the Neumann boundary conditions. Bound. Value Probl. 2012., 2012: Article ID 2. doi:10.1186/1687-2770-2012-2

    Google Scholar 

  17. Zhang W, Li J, Shi N: Stability analysis for stochastic Markovian jump reaction-diffusion neural networks with partially known transition probabilities and mixed time delays. Discrete Dyn. Nat. Soc. 2012., 2012: Article ID 524187. doi:10.1155/2012/524187

    Google Scholar 

  18. Song Q, Zhao Z, Li YM: Global exponential stability of BAM neural networks with distributed delays and reaction-diffusion terms. Phys. Lett. A 2005, 335(2-3):213-225. 10.1016/j.physleta.2004.12.007

    Article  Google Scholar 

  19. Zhang W, Li J: Global exponential stability of reaction-diffusion neural networks with discrete and distributed time-varying delays. Chin. Phys. B 2011., 20(3): Article ID 030701

    Google Scholar 

  20. Pan J, Zhan Y: On periodic solutions to a class of non-autonomously delayed reaction-diffusion neural networks. Commun. Nonlinear Sci. Numer. Simul. 2011, 16(1):414-422. 10.1016/j.cnsns.2010.02.022

    Article  MathSciNet  Google Scholar 

  21. Song Q, Cao J: Global exponential stability and existence of periodic colutions in BAM with delays and reaction-diffusion terms. Chaos Solitons Fractals 2005, 23(2):421-430. 10.1016/j.chaos.2004.04.011

    Article  MathSciNet  Google Scholar 

  22. Cui B, Lou X: Global asymptotic stability of BAM neural networks with distributed delays and reaction-diffusion terms. Chaos Solitons Fractals 2006, 27(5):1347-1354. 10.1016/j.chaos.2005.04.112

    Article  MathSciNet  Google Scholar 

  23. Zhao H, Wang G: Existence of periodic oscillatory solution of reaction-diffusion neural networks with delays. Phys. Lett. A 2005, 343(5):372-383. 10.1016/j.physleta.2005.05.098

    Article  MathSciNet  Google Scholar 

  24. Song Q, Cao J, Zhao Z: Periodic solutions and its exponential stability of reaction-diffusion recurrent neural networks with continuously distributed delays. Nonlinear Anal., Real World Appl. 2006, 7(1):65-80. 10.1016/j.nonrwa.2005.01.004

    Article  MathSciNet  Google Scholar 

  25. Wang Z, Zhang H: Global asymptotic stability of reaction-diffusion Cohen-Grossberg neural network with continuously distributed delays. IEEE Trans. Neural Netw. 2010, 21(1):39-49.

    Article  Google Scholar 

  26. Zhang X, Wu S, Li K: Delay-dependent exponential stability for impulsive Cohen-Grossberg neural networks with time-varying delays and reaction-diffusion terms. Commun. Nonlinear Sci. Numer. Simul. 2011, 16(3):1524-1532. 10.1016/j.cnsns.2010.06.023

    Article  MathSciNet  Google Scholar 

  27. Wang L, Zhang R, Wang Y: Global exponential stability of reaction-diffusion cellular neural networks with S -type distributed time delays. Nonlinear Anal., Real World Appl. 2009, 10(2):1101-1113. 10.1016/j.nonrwa.2007.12.002

    Article  MathSciNet  Google Scholar 

  28. Zhu Q, Li X, Yang X: Exponential stability for stochastic reaction-diffusion BAM neural networks with time-varying and distributed delays. Appl. Math. Comput. 2011, 217(13):6078-6091. 10.1016/j.amc.2010.12.077

    Article  MathSciNet  Google Scholar 

  29. Lou X, Cui B, Wu W: On global exponential stability and existence of periodic solutions for BAM neural networks with distributed delays and reaction-diffusion terms. Chaos Solitons Fractals 2008, 36(4):1044-1054. 10.1016/j.chaos.2006.08.005

    Article  MathSciNet  Google Scholar 

  30. Zhang W, Li J, Chen M: Dynamical behaviors of impulsive stochastic reaction-diffusion neural networks with mixed time delays. Abstr. Appl. Anal. 2012., 2012: Article ID 236562. doi:10.1155/2012/236562

    Google Scholar 

  31. Wang Z, Zhang H, Li P: An LMI approach to stability analysis of reaction-diffusion Cohen-Grossberg neural networks concerning Dirichlet boundary conditions and distributed delays. IEEE Trans. Syst. Man Cybern., Part B, Cybern. 2010, 40(6):1596-1606.

    Article  Google Scholar 

  32. Lu J: Robust global exponential stability for interval reaction-diffusion Hopfield neural networks with distributed delays. IEEE Trans. Circuits Syst. II, Express Briefs 2007, 54(12):1115-1119.

    Article  Google Scholar 

  33. Lu J, Lu L: Global exponential stability and periodicity of reaction-diffusion recurrent neural networks with distributed delays and Dirichlet boundary conditions. Chaos Solitons Fractals 2009, 39(4):1538-1549. 10.1016/j.chaos.2007.06.040

    Article  MathSciNet  Google Scholar 

  34. Wang J, Lu J: Global exponential stability of fuzzy cellular neural networks with delays and reaction-diffusion terms. Chaos Solitons Fractals 2008, 38(3):878-885. 10.1016/j.chaos.2007.01.032

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the referees for valuable comments and suggestions in improving this paper. This work is partially supported by the National Natural Science Foundation of China under Grant No. 60974139, the Special Research Project in Shaanxi Province Department of Education (2013JK0578) and Doctor Introduced project of Xianyang Normal University under Grant No. 12XSYK008.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Weiyuan Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

WZ designed and performed all the steps of proof in this research and also wrote the paper. JL and MC participated in the design of the study and suggested many good ideas that made this paper possible and helped to draft the first manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhang, W., Li, J. & Chen, M. Global exponential stability and existence of periodic solutions for delayed reaction-diffusion BAM neural networks with Dirichlet boundary conditions. Bound Value Probl 2013, 105 (2013). https://doi.org/10.1186/1687-2770-2013-105

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-105

Keywords