Skip to main content

Integro-differential inequality and stability of BAM FCNNs with time delays in the leakage terms and distributed delays

Abstract

In this paper, a class of impulsive bidirectional associative memory (BAM) fuzzy cellular neural networks (FCNNs) with time delays in the leakage terms and distributed delays is formulated and investigated. By establishing an integro-differential inequality with impulsive initial conditions and employing M-matrix theory, some sufficient conditions ensuring the existence, uniqueness and global exponential stability of equilibrium point for impulsive BAM FCNNs with time delays in the leakage terms and distributed delays are obtained. In particular, the estimate of the exponential convergence rate is also provided, which depends on the delay kernel functions and system parameters. It is believed that these results are significant and useful for the design and applications of BAM FCNNs. An example is given to show the effectiveness of the results obtained here.

1 Introduction

The bidirectional associative memory (BAM) neural network models were first introduced by Kosko [1]. It is a special class of recurrent neural networks that can store bipolar vector pairs. The BAM neural network is composed of neurons arranged in two layers, the X-layer and Y-layer. The neurons in one layer are fully interconnected to the neurons in the other layer. Through iterations of forward and backward information flows between the two layer, it performs a two-way associative search for stored bipolar vector pairs and generalize the single-layer autoassociative Hebbian correlation to a two-layer pattern-matched heteroassociative circuits. Therefore, this class of networks possesses good application prospects in some fields such as pattern recognition, signal and image process, and artificial intelligence [2]. In such applications, the stability of networks plays an important role; it is of significance and necessary to investigate the stability. It is well known, in both biological and artificial neural networks, the delays arise because of the processing of information. Time delays may lead to oscillation, divergence or instability which may be harmful to a system. Therefore, study of neural dynamics with consideration of the delayed problem becomes extremely important to manufacture high-quality neural networks. In recent years, there have been many analytical results for BAM neural networks with various axonal signal transmission delays, for example, see [311] and references therein. In addition, except various axonal signal transmission delays, time delay in the leakage term also has great impact on the dynamics of neural networks. As pointed out by Gopalsamy [12, 13], time delay in the stabilizing negative feedback term has a tendency to destabilize a system. Recently, some authors have paid attention to stability analysis of neural networks with time delays in the leakage (or "forgetting") terms [1218].

Since FCNNs were introduced by Yang et. al [19, 20], many researchers have done extensive works on this subject due to their extensive applications in classification of image processing and pattern recognition. Specially, in the past few years, the stability analysis on FCNNs with various delays and fuzzy BAM neural networks with transmission delays has been the highlight in the neural network field, for example, see [2127] and references therein. On the other hand, in respect of the complexity, besides delay effect, impulsive effect likewise exists in a wide variety of evolutionary processes in which states are changed abruptly at certain moments of time, involving such fields as medicine and biology, economics, mechanics, electronics and telecommunications. Many interesting results on impulsive effect have been gained, e.g., Refs. [2837]. As artificial electronic systems, neural networks such as CNNs, bidirectional neural networks and recurrent neural networks often are subject to impulsive perturbations, which can affect dynamical behaviors of the systems just as time delays. Therefore, it is necessary to consider both impulsive effect and delay effect on the stability of neural networks. To the best of our knowledge, few authors have considered impulsive BAM FCNNs with time delays in the leakage terms and distributed delays.

Motivated by the above discussions, the objective of this paper is to formulate and study impulsive BAM FCNNs with time delays in the leakage terms and distributed delays. Under quite general conditions, some sufficient conditions ensuring the existence, uniqueness and global exponential stability of equilibrium point are obtained by the topological degree theory, properties of M-matrix, the integro-differential inequality with impulsive initial conditions and analysis technique.

The paper is organized as follows. In Section 2, the new neural network model is formulated, and the necessary knowledge is provided. The existence and uniqueness of equilibrium point are presented in Section 3. In Section 4, we give some sufficient conditions of exponential stability of the impulsive BAM FCNNs with time delays in the leakage terms and distributed delays. An example is given to show the effectiveness of the results obtained here in Section 5. Finally, in Section 6, we give the conclusion.

2 Model description and preliminaries

In this section, we will consider the model of impulsive BAM FCNNs with time delays in the leakage terms and distributed delays, it is described by the following functional differential equation:

i ( t ) = - a i x i ( t - δ i ) + j = 1 m a i j g j ( y j ( t ) ) + j = 1 m ã i j v j + I i + j = 1 m α i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 m α ̃ i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 m T i j v j + j = 1 m H i j v j , t t k x i ( t + ) = x i ( t - ) + P i k ( x i ( t - ) ) , t = t k , k N = { 1 , 2 , } , j ( t ) = - b j y j ( t - θ j ) + i = 1 n b j i f i ( x i ( t ) ) + i = 1 n b ̃ j i u i + J j + i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x i ( t - s ) ) d s + i = 1 n β ̃ j i 0 + K ̄ i j ( s ) f i ( x i ( t - s ) ) d s + i = 1 n T ̄ j i u i + i = 1 n H ̄ j i u i , t t k y j ( t + ) = y j ( t - ) + Q j k ( y j ( t - ) ) , t = t k , k N = { 1 , 2 , } ,
(1)

for i = 1, 2, ..., n, j = 1, 2,..., m, t > 0, where x i (t) and y j (t) are the states of the i th neuron and the j th neuron at time t, respectively; δ i ≥ 0 and θ j ≥ 0 denote the leakage delays, respectively; f i and g j denote the signal functions of the i th neuron and the j th neuron at time t, respectively; u i , v j and I i , J j denote inputs and bias of the i th neuron and the j th neuron, respectively; a i > 0, b j > 0, a i j , ã i j , α i j , α ̃ i j , b j i , b ̃ j i , β j i , β ̃ j i are constants, a i and b j represent the rate with which the i th neuron and the j th neuron will reset their potential to the resting state in isolation when disconnected from the networks and external inputs, respectively; a ij , b ji and ã i j , b ̃ j i denote connection weights of feedback template and feedforward template, respectively; α ij , β ji and α ̃ i j , β ̃ j i denote connection weights of the distributed fuzzy feedback MIN template and the distributed fuzzy feedback MAX template, respectively; T i j , T ̄ j i and H i j , H ̃ j i are elements of fuzzy feedforward MIN template and fuzzy feedforward MAX template, respectively; and denote the fuzzy AND and fuzzy OR operations, respectively; K ij (s) and K ̄ j i ( s ) correspond to the delay kernel functions, respectively. t k is called impulsive moment and satisfies 0 < t1 < t2 < , lim k + t k =+; x i ( t k - ) and x i ( t k + ) denote the left-hand and right-hand limits at t k , respectively; P ik and Q jk show impulsive perturbations of the i th neuron and j th neuron at time t k , respectively.

We always assume x i ( t k + ) = x i ( t k ) and y j ( t k + ) = y j ( t k ) , k N . The initial conditions are given by

x i ( t ) = ϕ i ( t ) , - t 0 , y j ( t ) = φ j ( t ) , - t 0 ,

where ϕ i (t), φ j (t) (i = 1, 2, ..., n; j = 1, 2, ..., m) are bounded and continuous on (-∞, 0], respectively.

If the impulsive operators P ik (x i ) = 0, Q jk (y j ) = 0, i = 1, 2, ..., n, j = 1, 2, ..., m, k N, then system (1) may reduce to the following model:

i ( t ) = - a i x i ( t - δ i ) + j = 1 m a i j g j ( y j ( t ) ) + j = 1 m ã i j v j + I i + j = 1 m α i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 m α ̃ i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 m T i j v j + j = 1 m H i j v j , j ( t ) = - b j y j ( t - θ j ) + i = 1 n b j i f i ( x i ( t ) ) + i = 1 n b ̃ j i u i + J j + i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x i ( t - s ) ) d s + i = 1 n β ̃ j i 0 + K ̄ i j ( s ) f i ( x i ( t - s ) ) d s + i = 1 n T ̄ j i u i + i = 1 n H ̄ j i u i .
(2)

System (2) is called the continuous system of model (1).

Throughout this paper, we make the following assumptions:

(H1) For neuron activation functions f i and g j (i = 1, 2, ..., n; j = 1, 2, ..., m), there exist two positive diagonal matrices F = diag(F1, F2, ..., F n ) and G = diag(G1, G2, ..., G m ) such that

F i = sup x y f i ( x ) - f i ( y ) x - y , G j = sup x y g j ( x ) - g j ( y ) x - y

for all x, y R (xy).

(H2) The delay kernels K ij : [0, +∞) → R and K ̄ j i : [ 0 , + ) R are real-valued piecewise continuous, and there exists δ > 0 such that

k i j ( λ ) = 0 + e λ s | K i j ( s ) | d s , k ̄ j i ( λ ) = 0 + e λ s | K ̄ j i ( s ) | d s

Are continuous for λ [0,δ), i = 1,2, ..., n, j = 1,2, ..., m.

(H3) Let P ̄ k ( x ) =x+ P k ( x ) and Q ̄ k ( y ) =y+ Q k ( y ) be Lipschitz continuous in Rn and Rm , respectively, that is, there exist nonnegative diagnose matrices Γ k = diag(γ1k, γ2k, ..., γ nk ) and Γ ¯ k = diag ( γ ¯ 1 k , γ ¯ 2 k , , γ ¯ m k ) such that

| P ̄ k ( x ) - P ̄ k ( y ) | Γ k | x - y | , f o r a l l x , y R n , k N , | Q ̄ k ( u ) - Q ̄ k ( v ) | Γ ̄ k | u - v | , f o r a l l u , v R m , k N ,

where

P ̄ k ( x ) = ( P ̄ 1 k ( x 1 ) , P ̄ 2 k ( x 2 ) , , P ̄ n k ( x n ) ) T , Q ̄ k ( x ) = ( Q ̄ 1 k ( y 1 ) , Q ̄ 2 k ( y 2 ) , , Q ̄ m k ( y m ) ) T , P k ( x ) = ( P 1 k ( x 1 ) , P 2 k ( x 2 ) , , P n k ( x n ) ) T , Q k ( y ) = ( Q 1 k ( y 1 ) , Q 2 k ( y 2 ) , , Q m k ( y m ) ) T .

To begin with, we introduce some notation and recall some basic definitions.

PC[J, Rl ] = {z(t): JRl |z(t) is continuous at tt k , z ( t k + ) =z ( t k ) , and z ( t k - ) exists for t, t k J, k N}, where J R is an interval, l N.

PC = {ψ: (-∞, 0] → Rl | ψ(s) is bounded, and ψ(s+) = ψ(s) for s (-∞, 0), ψ(s-) exists for s (-∞, 0], ϕ(s-) = ϕ(s) for all but at most a finite number of points s (-∞, 0]}.

For an m × n matrix A, |A| denotes the absolute value matrix given by |A| = (|a ij |)m ×n. For A = (a ij )m × n, B = (b ij )m × n Rm × n, AB (A > B) means that each pair of corresponding elements of A and B such that the inequality a ij b ij (a ij > b ij ).

Definition 1 A function (x, y) T : (-∞, +∞) → Rn+mis said to be the special solution of system (1) with initial conditions

x ( s ) = ϕ ( s ) , y ( s ) = φ ( s ) s ( - , 0 ] ,

if the following two conditions are satisfied

  1. (i)

    (x, y) T is piecewise continuous with first kind discontinuity at the points t k , k K. Moreover, (x, y) T is right continuous at each discontinuity point.

  2. (ii)

    (x, y) T satisfies model (1) for t ≥ 0, and x(s) = ϕ(s), y(s) = φ(s) for s (-∞, 0].

Especially, a point (x*, y*) T Rn+mis called an equilibrium point of model (1), if (x(t), y(t)) T = (x*, y*) T is a solution of (1).

Throughout this paper, we always assume that the impulsive jumps P k and Q k satisfy (referring to [2837])

P k ( x * ) = 0 a n d Q k ( y * ) = 0 , k N ,

i.e.,

P ̄ k ( x * ) = x * a n d Q ̄ ( y * ) = y * , k N ,
(3)

where (x*, y*) T is the equilibrium point of continuous systems (2). That is, if (x*, y*) T is an equilibrium point of continuous system (2), then (x*, y*) T is also the equilibrium of impulsive system (1).

Definition 2 The equilibrium point (x*, y*) T of model (1) is said to be globally exponentially stable, if there exist constants λ > 0 and M ≥ 1 such that

| | x ( t ) - x * | | + | | y ( t ) - y * | | M ( | | ϕ - x * | | + | | φ - y * | | ) e - λ t

for all t ≥ 0, where (x(t), y(t)) T is any solution of system (1) with initial value (ϕ(s), φ(s)) T and

| | x ( t ) - x * | | = i = 1 n | x i ( t ) - x i * | , | | y ( t ) - y * | | = j = 1 m | y j ( t ) - y j * | , | | ϕ - x * | | = sup - < s 0 i = 1 n | ϕ i ( s ) - x i * | , | | φ - y * | | = sup - < s 0 j = 1 m | φ j ( s ) - y j * | .

Definition 3 A real matrix D = (d ij )n × nis said to be a nonsingular M-matrix if d ij ≤ 0, i, j = 1, 2, ..., n, ij, and all successive principal minors of D are positive.

Lemma 1[38]Let D = (d ij )n × nwith d ij ≤ 0 (ij), then the following statements are true:

  1. (i)

    D is a nonsingular M-matrix if and only if D is inverse-positive, that is, D -1 exists and D -1 is a nonnegative matrix.

  2. (ii)

    D is a nonsingular M-matrix if and only if there exists a positive vector ξ = (ξ 1, ξ 2, ..., ξ n ) T such that Dξ > 0.

Lemma 2[20]For any positive integer n, let h j : RR be a function (j = 1, 2, ..., n), then we have

| j = 1 n α j h j ( u j ) - j = 1 n α j h j ( v j ) | j = 1 n | α j | | h j ( u j ) - h j ( v j ) | , | j = 1 n α j h j ( u j ) - j = 1 n α j h j ( v j ) | j = 1 n | α j | | h j ( u j ) - h j ( v j ) |

for all α = (α1, α2, ..., α n ) T , u = (u1, u2, ..., u n ) T , v = (v1, v2, ..., v n ) T Rn .

3 Existence and uniqueness of equilibrium point

In this section, we will proof the existence and uniqueness of equilibrium point of model (1). For the sake of simplification, let

Ĩ i = j = 1 m ã i j v j + I i + j = 1 m T i j v j + j = 1 m H i j v j , i = 1 , 2 , , n , J ̃ j = i = 1 n b ̃ j i u i + J j + i = 1 n T ̄ j i u i + i = 1 n H ̄ j i u i , j = 1 , 2 , , m ,

then model (2) is reduced to

i ( t ) = - a i x i ( t - δ i ) + j = 1 m a i j g j ( y j ( t ) ) + j = 1 m α i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 m α ̃ i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + Ĩ i , j ( t ) = - b j y j ( t - θ j ) + i = 1 n b j i f i ( u i ( t ) ) + i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x i ( t - s ) ) d s + i = 1 n β ̃ j i 0 + K ̄ j i ( s ) f i ( u i ( t - s ) ) d s + J ̃ j .
(4)

It is evident that the dynamical characteristics of model (2) are as same as of model (4).

Theorem 1 Under assumptions (H1) and (H2), system (1) has one unique equilibrium point, if the following condition holds,

(C1) there exist vectors ξ = (ξ1, ξ2, ..., ξ n ) T > 0, η = (η1, η2, ..., η m ) T > 0 and positive number λ > 0 such that

( λ - a i e λ δ i ) ξ i + j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( λ ) G j η j < 0 , i = 1 , 2 , , n , ( λ - b j e λ θ j ) η j + i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( λ ) F i ξ i < 0 . j = 1 , 2 , , m .

Proof. Let h ( x 1 , , x n , y 1 , , y m ) = ( h 1 , , h n , h ¯ 1 , , h ¯ m ) T , where

h i = a i x i - j = 1 m a i j g j ( y j ) - j = 1 m α i j k i j ( 0 ) g j ( y j ) - j = 1 m α ̃ i j k i j ( 0 ) g j ( y j ) - Ĩ i , h ¯ j = b j y j - i = 1 n b j i f i ( x i ) - i = 1 n β j i k ̄ j i ( 0 ) f i ( x i ) - i = 1 n β ̃ j i k ̄ j i ( 0 ) f i ( x i ) - J ̃ j

for i = 1, 2, ..., n; j = 1, 2, ..., m. Obviously, from assumption (H2), the equilibrium points of model (4) are the solutions of system of equations:

h i = 0 , i = 1 , 2 , , n , h ¯ j = 0 , j = 1 , 2 , , m .
(5)

Define the following homotopic mapping:

H(x1, ..., x n , y1, ..., y m ) = θh(x1, ..., x n , y1, ..., y m ) + (1 - θ)(x1, ..., x n , y1, ..., y m ) T , where θ [0, 1]. Let H k (k = 1, 2, ..., n + m) denote the k th component of H(x1, ..., x n , y1, ..., y m ), then from assumption (H1) and Lemma 2, we have

| H i | [ 1 + θ ( a i - 1 ) ] | x i | - θ j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) G j | y j | - θ j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) | g j ( 0 ) | - θ | Ĩ i | , | H n + j | [ 1 + θ ( b j - 1 ) ] | y j | - θ i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) F i | x i | - θ i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) | f i ( 0 ) | - θ | J ̃ j |
(6)

for i = 1, 2, ..., n, j = 1, 2, ..., m. Denote

H ̄ = ( | H 1 | , | H 2 | , , | H n + m | ) T , z = ( | x 1 | , , | x n | , | y 1 | , , | y m | ) T , C = diag ( a 1 , , a n , b 1 , , b m ) , L = diag ( F 1 , , F n , G 1 , , G m ) , P = ( | Ĩ 1 | , , | Ĩ n , | , | J ̃ 1 | , , | J ̃ m | ) T Q = ( | f 1 ( 0 ) | , , | f n ( 0 ) | , | g 1 ( 0 ) | , , | g m ( 0 ) | ) T , A = | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) n × m , B = | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) m × n , T = 0 A B 0 , ω = ( ξ 1 , , ξ n , η 1 , , η m ) T > 0 .

Then, the matrix form of (6) is

H ̄ [ E + θ ( C - E ) ] z - θ T L z - θ ( P + T Q ) = ( 1 - θ ) z + θ [ ( C - T L ) z - ( P + T Q ) ] .

Since condition (C1) holds, and k ij (λ), k ̄ j i ( λ ) are continuous on [0, δ ), when λ = 0 in (C1), we obtain

- a i ξ i + j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) G j η j < 0 , i = 1 , 2 , , n , - b j η j + i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) F i ξ i < 0 . j = 1 , 2 , , m .

or in matrix form,

( - C + T L ) ω < 0 .
(7)

From Lemma 1, we know that C - TL is a nonsingular M-matrix, so (C - TL)-1 is a nonnegative matrix. Let

Γ = { z = ( x 1 , , x n , y 1 , , y m ) T | z ω + ( C T L ) 1 ( P + T Q ) } ,

then Γ is nonempty, and from (6), for any z = (x1,..., x n , y1,..., y m ) T ∂Γ, we have

H ̄ ( 1 - θ ) z + θ ( C - T L ) [ z - ( C - T L ) - 1 ( P + T Q ) ] (1) = ( 1 - θ ) [ ω + ( C - T L ) - 1 ( P + T Q ) ] + θ ( C - T L ) ω > 0 , θ [ 0 , 1 ] . (2) (3)

Therefore, for any (x1, ..., x n , y1, ..., y m ) T ∂Γ and θ [0, 1], we have H ≠ 0. From homotopy invariance theorem [39], we get

deg ( h , Γ , 0 ) = deg ( H , Γ , 0 ) = 1 ,

by topological degree theory, we know that (5) has at least one solution in Γ. That is, model (4) has at least an equilibrium point.

Now, we show that the solution of the system of Equations (5) is unique. Assume that ( x 1 * , , x n * , y 1 * , , y m * ) T and ( x ^ 1 , , x ^ n , ŷ 1 , , ŷ m ) T are two solutions of the system of Equations (5), then

a i ( x i * - x ^ i ) = j = 1 m a i j [ g j ( y j * ) + g j ( ŷ j ) ] + j = 1 m α i j k i j ( 0 ) g j ( y j * ) - j = 1 m α i j k i j ( 0 ) g j ( ŷ j ) + j = 1 m α ̃ i j k i j ( 0 ) g j ( y j * ) - j = 1 m α ̃ i j k i j ( 0 ) g j ( ŷ j ) , b j ( y j * - ŷ j ) = i = 1 n b j i [ f i ( x i * ) - f i ( x ^ i ) ] + i = 1 n β j i k ̄ j i ( 0 ) f i ( x i * ) - i = 1 n β j i k ̄ j i ( 0 ) f i ( x ^ i ) + i = 1 n β ̃ j i k ̄ j i ( 0 ) f i ( x i * ) - i = 1 n β ̃ j i k ̄ j i ( 0 ) f i ( x ^ i ) ,

it follows that

a i | x i * - x ^ i | j = 1 m | a i j | | g j ( y j * ) + g j ( ŷ j ) | + | j = 1 m α i j k i j ( 0 ) g j ( y j * ) - j = 1 m α i j k i j ( 0 ) g j ( ŷ j ) | + | j = 1 m α ̃ i j k i j ( 0 ) g j ( y j * ) - j = 1 m α ̃ i j k i j ( 0 ) g j ( ŷ j ) | , b j | y j * - ŷ j | i = 1 n | b j i | | f i ( x i * ) - f i ( x ^ i ) | + | i = 1 n β j i k ̄ j i ( 0 ) f i ( x i * ) - i = 1 n β j i k ̄ j i ( 0 ) f i ( x ^ i ) | + | i = 1 n β ̃ j i k ̄ j i ( 0 ) f i ( x i * ) - i = 1 n β ̃ j i k ̄ j i ( 0 ) f i ( x ^ i ) | .

By using of Lemma 2 and hypothesis (H1), we have

a i | x i * - x ^ i | - j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) G j | y j * - ŷ j | 0 , b j | y j * - ŷ j | - i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) F i | x i * - x ^ i | 0 .
(8)

Let Z = d i a g ( | x 1 * x ^ 1 | , , | x n * x ^ n | , | y 1 * y ^ 1 | , , | y m * y ^ m | ) , then the matrix form of (8) is (C -TL)Z ≤ 0. Since C - TL is a nonsingular M-matrix, (C - TL) -1 ≥ 0, thus Z ≤ 0, accordingly, Z = 0, i.e., x i * = x ^ i , y j * = ŷ j ( i = 1 , 2 , , n , j = 1 , 2 , , m ) . This shows that model (4) has one unique equilibrium point. According to (3), this implies that system (1) has one unique equilibrium point. The proof is completed.

Corollary 1 Under assumptions (H1) and (H2), system (1) has one unique equilibrium point if C - TL is a nonsingular M-matrix.

Proof. Since that C - TL is a nonsingular M-matrix, from Lemma 1, there exists a vector ω = (ξ1, ... ξ n , η1, ..., η m ) T > 0 such that (C TL) ω > 0, or (-C + TL) ω < 0. It follows that

- a i ξ i + j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( 0 ) G j η j < 0 , i = 1 , 2 , , n , - b j η j + i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( 0 ) F i ξ i < 0 , j = 1 , 2 , , m .

From the continuity of k ij (λ) and k ̄ j i ( λ ) , it is easy to know that there exists λ > 0 such that

( λ - a i e λ δ i ) ξ i + j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( λ ) G j η j < 0 , i = 1 , 2 , , n , ( λ - b j e λ θ j ) η j + i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( λ ) F i ξ i < 0 , j = 1 , 2 , , m .

That is, condition (C1) holds. This completes the proof.

4 Exponential stability and exponential convergence rate

In this section, we will discuss the global exponential stability of system (1) and give an estimation of exponential convergence rate.

Lemma 3 Let a < b ≤ +∞, and u(t) = (u1(t), ..., u n (t)) T PC[[a, b), Rn ] and v(t) = (v1(t), ..., v m (t)) T PC[[a, b), Rm ] satisfy the following integro-differential inequalities with the initial conditions u(s) PC[(-∞, 0], Rn ] and v(s) PC[(-∞, 0], Rm ]:

D + u i ( t ) - r i u i ( t - δ i ) + j = 1 m p i j v j ( t ) + j = 1 m q i j 0 + | K i j ( s ) | v j ( t - s ) d s , D + v j ( t ) - r ̄ j v j ( t - θ j ) + i = 1 n p ̄ j i u i ( t ) + i = 1 n q ̄ j i 0 + | K ̄ j i ( s ) | u i ( t - s ) d s
(9)

for i = 1, 2, ..., n, j = 1, 2, ..., m, where r i > 0, p ij > 0, q ij > 0, r ̄ j >0, p ̄ j i >0, q ̄ j i >0, i = 1, 2,...,n, j = 1, 2, ..., m. If the initial conditions satisfy

u ( s ) κ ξ e - λ ( s - a ) , s ( - , a ] , v ( s ) κ η e - λ ( s - a ) , s ( - , a ] ,
(10)

in which λ > 0, ξ = (ξ1, ξ2, ..., ξ n ) T > 0 and η = (η1, η2, ..., η m ) T > 0 satisfy

( λ - r i e λ δ i ) ξ i + j = 1 m ( p i j + q i j k i j ( λ ) ) η j < 0 , i = 1 , 2 , , n , ( λ - r ̄ j e λ θ j ) η j + i = 1 n ( p ̄ j i + q ̄ j i k ̄ j i ( λ ) ) ξ i < 0 , j = 1 , 2 , , m .
(11)

Then

u ( t ) κ ξ e - λ ( t - a ) , t [ a , b ) , v ( t ) κ η e - λ ( t - a ) , t [ a , b ) .

Proof. For i {1, 2, ..., n}, j {1, 2, ..., m} and arbitrary ε > 0, set z i (t) = (κ + ε) ξ i e-λ(t - a), z ̄ j ( t ) = ( κ + ε ) η j e - λ ( t - a ) , we prove that

u i ( t ) z i ( t ) = ( κ + ε ) ξ i e - λ ( t - a ) , t [ a , b ) , i = 1 , 2 , , n , v j ( t ) z ̄ j ( t ) = ( κ + ε ) η j e - λ ( t - a ) , t [ a , b ) , j = 1 , 2 , , m .
(12)

If this is not true, no loss of generality, suppose that there exist i0 and t* [a, b) such that

u i 0 ( t * ) = z i 0 ( t * ) , D + u i 0 ( t * ) ż i 0 ( t * ) , u i ( t ) z i ( t ) , v j ( t ) z ̄ j ( t )
(13)

for t [a, t*], i = 1, 2,..., n, j = 1, 2,..., m.

However, from (9) and (12), we get

D + u i 0 ( t * ) - r i 0 u i 0 ( t * - δ i 0 ) + j = 1 m p i 0 j v j ( t * ) + j = 1 m q i 0 j 0 + | K i 0 j ( s ) | v j ( t * - s ) d s - r i 0 ( κ + ε ) ξ i 0 e - λ ( t * - δ i 0 - a ) + j = 1 m p i 0 j η j ( κ + ε ) η j e - λ ( t * - a ) + j = 1 m q i 0 j ( κ + ε ) η j e - λ ( t * - a ) 0 + e λ s | K i 0 j ( s ) | d s = [ - r i 0 ξ i 0 e λ δ i 0 + j = 1 m ( p i 0 j + q i 0 j k i 0 j ( λ ) ) η j ] ( κ + ε ) e - λ ( t * - a ) .

Since (11) holds, it follows that - r i 0 ξ i 0 e λ δ i 0 + j = 1 m ( p i 0 j + q i 0 j k i 0 j ( λ ) ) η j <-λ ξ i 0 <0. Therefore, we have

D + u i 0 ( t * ) < - λ ξ i 0 ( κ + ε ) e - λ ( t * - a ) = ż i 0 ( t * ) ,

which contradicts the inequality D + u i 0 ( t * ) ż i 0 ( t * ) in (13). Thus (12) holds for all t [a, b). Letting ε → 0, we have

u i ( t ) κ ξ i e - λ ( t - a ) , t [ a , b ) , i = 1 , 2 , , n , v j ( t ) κ η j e - λ ( t - a ) , t [ a , b ) , j = 1 , 2 , , m .

The proof is completed.

Remark 1. Lemma 3 is a generalization of the famous Halanay inequality.

Theorem 2 Under assumptions (H1)-(H3), if the following conditions hold,

(C1) there exist vectors ξ = (ξ1, ξ2, ..., ξ n ) T > 0, η = (η1, η2, ..., η m ) T > 0 and positive number λ > 0 such that

( λ - a i e λ δ i ) ξ i + j = 1 m | a i j | + ( | α i j | + | α ̃ i j | ) k i j ( λ ) G j η j < 0 , i = 1 , 2 , , n , ( λ - b j e λ θ j ) η j + i = 1 n | b j i | + ( | β j i | + | β ̃ j i | ) k ̄ j i ( λ ) F i ξ i < 0 , j = 1 , 2 , , m ;

(C2) μ= sup k N { ln μ k t k - t k - 1 } <λ, where μ k = max 1 i n , 1 j m { 1 , γ i k , γ ̄ j k } , k N,

then system (1) has exactly one globally exponentially stable equilibrium point, and its exponential convergence rate equals λ - μ.

Proof. Since (C1) holds, from Theorem 1, we know that system (1) has one unique equilibrium point ( x 1 * , , x n * , y 1 * , , y m * ) T . Now, we assume that (x1(t), ..., x n (t), y1(t), ..., y m (t)) T is any solution of system (1), let x ̄ i ( t ) = x i ( t ) - x i * , i = 1, 2, ..., n, ȳ j ( t ) = y j ( t ) - y j * , j = 1, 2, ..., m. It is easy to see that system (1) can be transformed into the following system

x ̄ ° i ( t ) = - a i x ̄ i ( t - δ i ) + j = 1 m a i j g j ( ȳ j ( t ) + y j * ) - g j ( y j * ) + j = 1 m α i j 0 + K i j ( s ) g j ( ȳ j ( t - s ) + y j * ) d s - j = 1 m α i j 0 + K j i ( s ) g j ( y j * ) d s + j = 1 m α ̃ i j 0 + K i j ( s ) g j ( ȳ j ( t - s ) + y j * ) d s - j = 1 m α ̃ i j 0 + K i j ( s ) g j ( y j * ) d s , t t k , x ̄ i ( t k + ) = P ̃ i k ( x ̄ i ( t k - ) ) , k N ȳ ° j ( t ) = - b j ȳ j ( t - θ j ) + i = 1 n b j i ( f i ( x ̄ i ( t ) + x i * ) - f i ( x i * ) ) + i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x ̄ i ( t - s ) + x i * ) d s - i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x i * ) d s + i = 1 n β ̃ j i 0 + K ̄ i j ( s ) f i ( x ̄ i ( t - s ) + x i * ) d s - i = 1 n β ̃ j i 0 + K ̄ i j ( s ) f i ( x i * ) d s , t t k , ȳ j ( t k + ) = Q ̃ j k ( y j ( t k - ) ) , k N ,
(14)

where P ̃ i k ( x ̄ i ( t ) ) = P ̄ i k ( x ̄ i ( t ) + x i * ) - P ̄ i k ( x i * ) , Q ̃ j k ( ȳ j ( t ) ) = Q ̄ j k ( ȳ j ( t ) + y j * ) - Q ̄ j k ( y j * ) , and the initial conditions of (14) are

ϕ ̃ ( s ) = x ( s ) - x * = ϕ ( s ) - x * , s ( - , 0 ] , φ ̃ ( s ) = y ( s ) - y * = φ ( s ) - y * , s ( - , 0 ] .

From (H1) and Lemma 2, we calculate the upper right derivative along the solutions of first equation and third equation of (14), we can obtain

D + | x ̄ i ( t ) | - a i | x ̄ i ( t - δ i ) | + j = 1 m | a i j | G j | ȳ j ( t ) | + j = 1 m ( | α i j | + | α ̃ i j | ) G j 0 + | K i j ( s ) | | ȳ j ( t - s ) | d s , D + | ȳ j ( t ) | - b j | ȳ j ( t - θ j ) | + i = 1 n | b j i | F i | x ̄ i ( t ) | + i = 1 n ( | β j i | + | β ̃ j i | ) F i 0 + | K ̄ j i ( s ) | | x ̄ i ( t - s ) | d s

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

Let u i ( t ) =| x ̄ i ( t ) |, v j ( t ) =| ȳ j ( t ) |, r i = a i , p ij = |a ij |G j , q i j = ( | α i j | + | α ̃ i j | ) G j , r ̄ j = b j , p ̄ j i =| b j i | F i , q ̄ j i = ( | β j i | + | β ̃ j i | ) F i ( i = 1 , 2 , , n ; j = 1 , 2 , , m ) , then we have

D + u i ( t ) - r i u i ( t - δ i ) + j = 1 m p i j v j ( t ) + j = 1 m q i j 0 + | K i j ( s ) | v j ( t - s ) d s , D + v j ( t ) - r ̄ j v j ( t - θ j ) + i = 1 n p ̄ j i u i ( t ) + i = 1 n q ̄ j i 0 + | K ̄ j i ( s ) | u i ( t - s ) d s
(15)

for i = 1, 2, ..., n, j = 1, 2, ..., m, and from (C1), there exist vectors ξ = (ξ1, ξ2, ..., ξ n ) T > 0, η = (η1, η2, ..., η m ) T > 0 and positive number λ > 0 such that

( λ - r i e λ δ i ) ξ i + j = 1 m p i j + q i j k i j ( λ ) G j η j < 0 , i = 1 , 2 , , n , ( λ - r ̄ j e λ θ j ) η j + i = 1 n p ̄ j i + q ̄ j i k ̄ j i ( λ ) F i ξ i < 0 , j = 1 , 2 , , m .
(16)

Taking | | ϕ ̃ | | + | | φ ̃ | | min 1 i n , 1 j m { ξ i , η j } , it is easy to prove that

u ( t ) κ ξ e - λ t , - t 0 = t 0 , v ( t ) κ η e - λ t , - t 0 = t 0 .
(17)

From Lemma 3, we obtain that

u ( t ) κ ξ e - λ t , t 0 t < t 1 , v ( t ) κ η e - λ t , t 0 t < t 1 .
(18)

Suppose that for lk, the inequalities

u ( t ) κ μ 0 μ 1 μ l - 1 ξ e - λ t , t l - 1 t < t l , v ( t ) κ μ 0 μ 1 μ l - 1 η e - λ t , t l - 1 t < t l .
(19)

hold, where μ0 = 1. When l = k + 1, we note that

u ( t k ) = | P ̃ k ( u ( t k - ) ) | Γ k u ( t k - ) κ μ 0 μ 1 μ k - 1 Γ k ξ lim t t k - e - λ t (1) κ μ 0 μ 1 μ k - 1 μ k ξ e - λ t k , (2) (3)
(20)

and

v ( t k ) = | Q ̃ k ( v ( t k - ) ) | Γ ̄ k v ( t k - ) κ μ 0 μ 1 μ k - 1 Γ ̄ k η lim t t k - e - λ t (1) κ μ 0 μ 1 μ k - 1 μ k η e - λ t k . (2) (3)
(21)

From (20), (21) and μ k ≥ 1, we have

u ( t ) κ μ 0 μ 1 μ k - 1 μ k ξ e - λ t , - t t k , v ( t ) κ μ 0 μ 1 μ k - 1 μ k η e - λ t , - t t k .
(22)

Combining (15),(16),(22) and Lemma 3, we obtain that

u ( t ) κ μ 0 μ 1 μ k ξ e - λ t , t k t < t k + 1 , v ( t ) κ μ 0 μ 1 μ k η e - λ t , t k t < t k + 1 .
(23)

Applying the mathematical induction, we can obtain the following inequalities

u ( t ) κ μ 0 μ 1 μ k ξ e - λ t , t [ t k , t k + 1 ) , k N , v ( t ) κ μ 0 μ 1 μ k η e - λ t , t [ t k , t k + 1 ) , k N .
(24)

According to (C2), we have μ k e μ ( t k - t k - 1 ) < e λ ( t k - t k - 1 ) , so we have

u ( t ) κ e μ t 1 e μ ( t 2 - t 1 ) e μ ( t k - 1 - t k - 2 ) ξ e - λ t (1) = κ ξ e μ t k - 1 e - λ t κ ξ e - ( λ - μ ) t , t [ t k - 1 , t k ) , k N , (2) (3)

and

v ( t ) κ e μ t 1 e μ ( t 2 - t 1 ) e μ ( t k - 1 - t k - 2 ) η e - λ t (1) = κ η e μ t k - 1 e - λ t κ η e - ( λ - μ ) t , t [ t k - 1 , t k ) , k N . (2) (3)

That is

u ( t ) κ ξ e - ( λ - μ ) t , t ( - , t k ) , k N , v ( t ) κ η e - ( λ - μ ) t , t ( - , t k ) , k N .
(25)

It follows that

i = 1 n | x i ( t ) - x i * | + j = 1 m | y j ( t ) - y j * | = i = 1 n u i ( t ) + j = 1 m v j ( t ) (1) i = 1 n κ ξ i e - ( λ - μ ) t + j = 1 m κ η j e - ( λ - μ ) t (2) = i = 1 n ξ i + j = 1 m η j min 1 i n , 1 j m { ξ i , η j } ( | | ϕ ̃ | | + | | φ ̃ | | ) e - ( λ - μ ) t (3) = M ( | | ϕ - x * | | + | | φ - y * | | ) e - ( λ - μ ) t , (4) (5) 

where M= i = 1 n ξ i + j = 1 m η j min 1 i n , 1 j m { ξ i , η j } , then we have

| | x ( t ) - x * | | + | | y ( t ) - y * | | M | | ϕ - x * | | + | | φ - y * | | e - ( λ - μ ) t .

The proof is completed.

Remark 2. In Theorem 2, the parameters μ k and μ depend on the impulsive disturbance of system (1), and λ is actually an estimate of exponential convergence rate of continuous system (2), which depends on the delay kernel functions and system parameters. In order to obtain more precise estimate of the exponential convergence rate of system (1) (or system (2)), we suggest the following optimization problem:

( OP ) { max λ , s . t .  (C 1 ) h o l d s .

Obviously, for continuous system (2), we can immediately obtain the following corollaries.

Corollary 2 Under assumptions (H1) and (H2), if condition (C1) holds, then system (2) has exactly one globally exponentially stable equilibrium point, and its exponential convergence rate equals λ.

Corollary 3 Under assumptions (H1) and (H2), system (2) has exactly one globally exponentially stable equilibrium point if C - TL is a nonsingular M-matrix.

Remark 3. Note that Lemma 2 transforms the fuzzy AND () and the fuzzy OR () operation into the SUM operation (∑). So above results can be applied to the following classical impulsive BAM neural networks with time delays in the leakage terms and distributed delays:

i ( t ) = - a i x i ( t - δ i ) + j = 1 m a i j g j ( y j ( t ) ) + j = 1 m α i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + I i , t t k x i ( t + ) = x i ( t - ) + P i k ( x i ( t - ) ) , t = t k , k N , j ( t ) = - b j y j ( t - θ j ) + i = 1 n b j i f i ( x i ( t ) ) + i = 1 n β j i 0 + K ̄ j i ( s ) f i ( x i ( t - s ) ) d s + J j , t t k y j ( t + ) = y j ( t - ) + Q j k ( y j ( t - ) ) , t = t k , k N
(26)

for i = 1, 2,..., n; j = 1, 2,..., m.

For model (26), it is easy to obtain the following result:

Theorem 3 Under assumptions (H1)-(H3), if the following conditions hold,

(C1') there exist vectors ξ = (ξ1, ξ2, ..., ξ n ) T > 0, η = (η1, η2, ..., η m ) T > 0 and positive number λ > 0 such that

( λ - a i e λ δ i ) ξ i + j = 1 m | a i j | + | α i j | k i j ( λ ) G j η j < 0 , i = 1 , 2 , , n , ( λ - b j e λ θ j ) η j + i = 1 n | b j i | + | β j i | k ̄ j i ( λ ) F i ξ i < 0 , j = 1 , 2 , , m ;

(C2) μ = sup k N ln μ k t k - t k - 1 < λ , where μ k = max 1 i n , 1 j m { 1 , γ i k , γ ̄ j k } , k N ,

then system (26) has exactly one globally exponentially stable equilibrium point, and its exponential convergence rate equals λ - μ.

5 An illustrative example

In order to illustrate the feasibility of our above-established criteria in the preceding sections, we provide a concrete example. Although the selection of the coefficients and functions in the example is somewhat artificial, the possible application of our theoretical theory is clearly expressed.

Example. Consider the following impulsive BAM FCNNs with time delays in the leakage terms and distributed delays:

i ( t ) = - a i x i ( t - δ i ) + j = 1 2 a i j g j ( y j ( t ) ) + j = 1 2 ã i j v j + I i + j = 1 2 α i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 2 α ̃ i j 0 + K i j ( s ) g j ( y j ( t - s ) ) d s + j = 1 2 T i j v j + j = 1 2 H i j v j , t t k x i ( t + ) = x i ( t - ) + P i k ( x i ( t - ) ) = x i ( t - ) - ( 1 + e 0 . 1 2 5 k ) ( x i ( t - ) - 1 ) , t = t k , j ( t ) = - b j y j ( t - θ j ) + i = 1 2 b j i f i ( x i ( t ) ) + i = 1 2 b ̃ j i u i + J j + i = 1 2 β j i 0 + K ̄ j i ( s ) f i ( x i ( t - s ) ) d s + i = 1 2 β ̃ j i 0 + K ̄ i j ( s ) f i ( x i ( t - s ) ) d s + i = 1 2 T ̄ j i u i + i = 1 2 H ̄ j i u i , t t k y j ( t + ) = y j ( t - ) + Q j k ( y j ( t - ) ) = y j ( t - ) - ( 1 + e 0 . 1 2 5 k ) ( y j ( t - ) - 1 ) , t = t k
(27)

for k N, i = 1, 2, j = 1, 2, t > 0, t0 = 0, t k = tk-1+ 0.5k, k N, where

a 1 = 4 . 5 , a 2 = 4 . 5 , δ 1 = 0 . 2 , δ 2 = 0 . 3 , a 1 1 = 4 3 , a 1 2 = - 1 2 , a 2 1 = 1 2 , a 2 2 = 2 3 , ã 1 1 = 1 , ã 1 2 = - 2 , ã 2 1 = - 2 , ã 2 2 = 1 , I 1 = 6 7 1 2 , I 2 = 5 1 2 , α 1 1 = 1 3 , α 1 2 = - 1 4 , α 2 1 = 1 4 , α 2 2 = 2 3 , α ̃ 1 1 = 1 3 , α ̃ 1 2 = 1 4 , α ̃ 2 1 = - 1 4 , α ̃ 2 2 = 2 3 , T 1 1 = 1 , T 1 2 = 0 , T 2 1 = 0 , T 2 2 = 1 , H 1 1 = 1 , H 1 2 = 0 , H 2 1 = 0 , H 2 2 = 1 , v 1 = 1 , v 2 = 2 ; b 1 = 4 . 5 , b 2 = 4 . 5 , θ 1 = 0 . 2 , θ 2 = 0 . 1 , b 1 1 = 1 3 , b 1 2 = - 2 3 , b 2 1 = 4 3 , b 2 2 = 1 3 , b ̃ 1 1 = 1 , b ̃ 1 2 = 3 , b ̃ 2 1 = 2 , b ̃ 2 2 = - 2 , J 1 = - 1 2 , J 2 = 7 6 , β 1 1 = 1 3 , β 1 2 = - 1 6 , β 2 1 = 1 3 , β 2 2 = 1 3 , β ̃ 1 1 = 1 3 , β ̃ 1 2 = 1 6 , β ̃ 2 1 = 1 3 β ̃ 2 2 = 1 3 , T ̃ 1 1 = 1 , T ̃ 1 2 = 0 , T ̃ 2 1 = 0 , T ̃ 2 2 = 1 , H ̃ 1 1 = 1 , H ̃ 1 2 = 0 , H ̃ 2 1 = 0 , H ̃ 2 2 = 1 , u 1 = 1 , u 2 = 1 ; K i j ( s ) = K ̄ i j ( s ) = e - s , f i ( s ) = g j ( s ) = | s + 1 | - | s - 1 | 2 , i , j = 1 , 2 .

From above parameters, we have F1 = F2 = 1, G1 = G2 = 1, and ( k i j ( λ ) ) 2 × 2 = ( k ¯ j i ( λ ) ) 2 × 2 = ( 1 1 λ 1 1 λ 1 1 λ 1 1 λ ) , Γ k = Γ ¯ k = ( e 0.125 k e 0.125 k ) .

Solving the following optimization problem

{ max λ 0 > ( λ a 1 e λ δ 1 ) ξ 1 + ( | a 11 | + ( | α 11 | + | α ˜ 11 | ) k 11 ( λ ) ) G 1 η 1 + ( | a 12 | + ( | α 12 | + | α ˜ 12 | ) k 12 ( λ ) ) G 2 η 2 , 0 > ( λ a 2 e λ δ 2 ) ξ 1 + ( | a 21 | + ( | α 21 | + | α ˜ 21 | ) k 21 ( λ ) ) G 1 η 1 + ( | a 22 | + ( | α 22 | + | α ˜ 22 | ) k 22 ( λ ) ) G 2 η 2 0 > ( λ b 1 e λ θ 1 ) η 1 + ( | b 11 | + ( | β 11 | + | β ˜ 11 | ) k ¯ 11 ( λ ) ) F 1 ξ 1 + ( | b 12 | + ( | β 12 | + | β ˜ 12 | ) k ¯ 12 ( λ ) ) F 2 ξ 2 0 > ( λ b 2 e λ θ 2 ) η 2 + ( | b 21 | + ( | β 21 | + | β ˜ 21 | ) k ¯ 21 ( λ ) ) F 1 ξ 1 + ( | b 22 | + ( | β 22 | + | β ˜ 22 | ) k ¯ 22 ( λ ) ) F 2 ξ 2 , λ > 0 , ξ = ( ξ 1 , ξ 2 ) T > 0 , η = ( η 1 , η 2 ) T > 0 .

We obtain that λ ≈ 0.3868 > 0, ξ = (1082041, 1327618) T > 0 and η = (716212, 1050021) T > 0, so (C1) holds. From Theorem 1, we know system (27) has a unique equilibrium point, this equilibrium point is (1, 1, 1, 1) T . Also,

μ k = max 1 i 2 , 1 j 2 { 1 , γ i k , γ ̄ j k } = e 0 . 1 2 5 k , (1) μ = sup k N ln μ k t k - t k - 1 = 0 . 1 2 5 k 0 . 5 k = 0 . 2 5 < 0 . 3 8 6 8 = λ . (2) (3)

That is, (C2) holds. From Theorem 2, the unique equilibrium point (1, 1, 1, 1) T of system (27) is globally exponentially stable, and its exponential convergence rate is about 0.1368. The numerical simulation is shown in Figure 1 and 2.

Figure 1
figure 1

Behavior of the state variable x ( t ) with time impulses.

Figure 2
figure 2

Behavior of the state variable y ( t ) with time impulses.

6 Conclusions

In this paper, a class of impulsive BAM FCNNs with time delays in the leakage terms and distributed delays has been formulated and investigated. Some new criteria on the existence, uniqueness and global exponential stability of equilibrium point for the networks have been derived by using M-matrix theory and the impulsive delay integro-differential inequality. Our stability criteria are delay-dependent and impulse-dependent. The neuronal output activation functions and the impulsive operators only need to are Lipschitz continuous, but need not to be bounded and monotonically increasing. Some restrictions of delay kernel functions are also removed. It is worthwhile to mention that our technical methods are practical, in the sense that all new stability conditions are stated in simple algebraic forms and provided a more precise estimate of the exponential convergence rate, so their verification and applications are straightforward and convenient. The effectiveness of our results has been demonstrated by the convenient numerical example.

References

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

    Article  MathSciNet  Google Scholar 

  2. Mathai G, Upadhyaya BR: Performance analysis and application of the bidirectional associative memory to industrial spectral signatures. Proceedings of the IJCNN-1989 1989, 1: 33–37.

    Google Scholar 

  3. Li C, Liao X, Zhang R: Delay-dependent exponential stability analysis of bi-directional associative memory neural networks with time delay: an LMI approach. Chaos Solitons Fractals 2005, 24: 1119–1134. 10.1016/j.chaos.2004.09.052

    Article  MATH  MathSciNet  Google Scholar 

  4. Arik S: Global asymptotic stability analysis of bidirectional associative memory neural networks with time delays. IEEE Trans Neural Netw 2005, 16: 580–586. 10.1109/TNN.2005.844910

    Article  Google Scholar 

  5. Gopalsamy K, He X: Delay-independent stability in bi-directional associative memory with axonal signal transmission delays. IEEE Trans Neural Netw 1994, 5: 998–1002. 10.1109/72.329700

    Article  Google Scholar 

  6. Mohamad S: Global exponential stability in continuous-time and discrete-time delay bidirectional neural networks. Phys D 2001, 159: 233–251. 10.1016/S0167-2789(01)00344-X

    Article  MATH  MathSciNet  Google Scholar 

  7. Zhang J, Yang Y: Global stability analysis of bi-directional associative memory neural networks with time delay. Int J Circ Theory Appl 2001, 29: 185–196. 10.1002/cta.144

    Article  MATH  Google Scholar 

  8. Liang J, Cao J: Exponential stability of continuous-time and discrete-time bidirectional associative memory networks with delays. Chaos Solitons Fractals 2004, 22: 773–785. 10.1016/j.chaos.2004.03.004

    Article  MATH  MathSciNet  Google Scholar 

  9. Feng C, Plamondon R: Stability analysis of bidirectional associative memory networks with time delays. IEEE Trans Neural Netw 2003, 14: 1560–1565. 10.1109/TNN.2003.820829

    Article  Google Scholar 

  10. Rao VSH, Phaneendra BhRM: Global dynamics of bidirectional associative memory neural networks involving transmission delays and dead zones. Neural Netw 1999, 12: 455–465. 10.1016/S0893-6080(98)00134-8

    Article  Google Scholar 

  11. Song Q, Cao J: Global exponential stability of bidirectional associative memory neural networks with distributed delays. J Comput Appl Math 2007, 202: 266–279. 10.1016/j.cam.2006.02.031

    Article  MATH  MathSciNet  Google Scholar 

  12. Gopalsamy K: Stability and Oscillations in Delay Differential Equations of Population Dynamics. Kluwer, Dordrecht; 1992.

    Chapter  Google Scholar 

  13. Gopalsamy K: Leakage delays in BAM. J Math Anal Appl 2007, 325: 1117–1132. 10.1016/j.jmaa.2006.02.039

    Article  MATH  MathSciNet  Google Scholar 

  14. Peng S: Global attractive periodic solutions of BAM neural networks with continuously distributed delays in the leakage terms. Nonlinear Anal Real World Appl 2010, 11: 2141–2151. 10.1016/j.nonrwa.2009.06.004

    Article  MATH  MathSciNet  Google Scholar 

  15. Li X, Fu X, Balasubramaniam P, Rakkiyappan R: Existence, uniqueness and stability analysis of recurrent neural networks with time delay in the leakage term under impulsive perturbations. Nonlinear Anal Real World Appl 2010, 11: 4092–4108. 10.1016/j.nonrwa.2010.03.014

    Article  MATH  MathSciNet  Google Scholar 

  16. Li X, Cao J: Delay-dependent stability of neural networks of neutral type with time delay in the leakage term. Nonlinearity 2010, 23: 1709–1726. 10.1088/0951-7715/23/7/010

    Article  MATH  MathSciNet  Google Scholar 

  17. Balasubramaniam P, Kalpana M, Rakkiyappan R: Global asymptotic stability of BAM fuzzy cellular neural networks with time delay in the leakage term, discrete and unbounded distributed delays. Math Comput Model 2011, 53: 839–853. 10.1016/j.mcm.2010.10.021

    Article  MATH  MathSciNet  Google Scholar 

  18. Li X, Rakkiyappan R, Balasubramaniam P: Existence and global stability analysis of equilibrium of fuzzy cellular neural networks with time delay in the leakage term under impulsive perturbations. J Franklin Institute 2011, 348: 135–155. 10.1016/j.jfranklin.2010.10.009

    Article  MATH  MathSciNet  Google Scholar 

  19. Yang T, Yang L, Wu C, Chua L: Fuzzy cellular neural networks: theory. Proceedings of IEEE International Workshop on Cellular Neural Networks and Applications 1996, 181–186.

    Google Scholar 

  20. Yang T, Yang L: The global stability of fuzzy cellular neural networks. IEEE Trans Cric Syst I 1996, 43: 880–883. 10.1109/81.538999

    Google Scholar 

  21. Liu Y, Tang W: Exponential stability of fuzzy cellular neural networks with constant and time-varying delays. Phys Lett A 2004, 323: 224–233. 10.1016/j.physleta.2004.01.064

    Article  MATH  MathSciNet  Google Scholar 

  22. Yuan K, Cao J, Deng J: Exponential stability and periodic solutions of fuzzy cellular neural networks with time-varying delays. Neurocomputing 2006, 69: 1619–1627. 10.1016/j.neucom.2005.05.011

    Article  Google Scholar 

  23. Huang T: Exponential stability of fuzzy cellular neural networks with distributed delay. Phys Lett A 2006, 351: 48–52. 10.1016/j.physleta.2005.10.060

    Article  MATH  Google Scholar 

  24. Li K: Impulsive effect on global exponential stability of BAM fuzzy cellular neural networks with time-varying delays. Int J Syst Sci 2010, 41: 131–142. 10.1080/00207720903042921

    Article  MATH  Google Scholar 

  25. Song Q, Cao J: Dynamical behaviors of discrete-time fuzzy cellular neural networks with variable delays and impulses. J Franklin Institute 2007, 345: 39–59.

    Article  MathSciNet  Google Scholar 

  26. Wang S, Chung KFL, Fu D: Applying the improved fuzzy cellular neural network IFCNN to white blood cell detection. Neurocomputing 2007, 70: 1348–1359. 10.1016/j.neucom.2006.07.012

    Article  Google Scholar 

  27. Li K, Zhang L: Global exponential stability of BAM fuzzy cellular neural networks with distributed delays and impulses. J Appl Math Inform 2011, 29: 211–225.

    MATH  MathSciNet  Google Scholar 

  28. Akca H, Alassar R, Covachev V, Covacheva Z, Al-Zahrani E: Continuous-time additive Hopfield-type neural networks with impulses. J Math Anal Appl 2004, 290: 436–451. 10.1016/j.jmaa.2003.10.005

    Article  MATH  MathSciNet  Google Scholar 

  29. Xu D, Yang Z: Impulsive delay differential inequality and stability of neural networks. J Math Anal Appl 2005, 305: 107–120. 10.1016/j.jmaa.2004.10.040

    Article  MATH  MathSciNet  Google Scholar 

  30. Yang Z, Xu D: Impulsive effects on stability of Cohen-Grossberg neural networks with variable delays. Appl Math Comput 2006, 177: 63–78. 10.1016/j.amc.2005.10.032

    Article  MATH  MathSciNet  Google Scholar 

  31. Gopalsamy K: Stability of artificial neural networks with impulses. Appl Math Comput 2004, 154: 783–813. 10.1016/S0096-3003(03)00750-1

    Article  MATH  MathSciNet  Google Scholar 

  32. Mohamad S, Gopalsamy K, Akca H: Exponential stability of artificial neural networks with distributed delays and large impulses. Nonlinear Anal Real World Appl 2008, 9: 872–888. 10.1016/j.nonrwa.2007.01.011

    Article  MATH  MathSciNet  Google Scholar 

  33. Li Y: Global exponential stability of BAM neural networks with delays and impulses. Chaos, Solitions Fractals 2005, 24: 279–285.

    Article  MATH  Google Scholar 

  34. Li Y, Yang C: Global exponential stability analysis on impulsive BAM neural networks with distributed delays. J Math Anal Appl 2006, 324: 1125–1139. 10.1016/j.jmaa.2006.01.016

    Article  MATH  MathSciNet  Google Scholar 

  35. Liu B, Huang L: Global exponential stability of BAM neural networks with recent-history distributed delays and impulses. Neurocomputing 2006, 69: 2090–2096. 10.1016/j.neucom.2005.09.014

    Article  Google Scholar 

  36. Huang Z, Luo X, Yang Q: Global asymptotic stability analysis of bidirectional associative memory neural networks with distributed delays and impulse. Chaos, Solitons Fractals 2007, 34: 878–885. 10.1016/j.chaos.2006.03.112

    Article  MATH  MathSciNet  Google Scholar 

  37. Zhou Q: Global exponential stability of BAM neural networks with distributed delays and impulses. Nonlinear Anal Real World Appl 2009, 10: 144–153. 10.1016/j.nonrwa.2007.08.019

    Article  MATH  MathSciNet  Google Scholar 

  38. Horn RA, Johnson CR: Topics in Matrix Analysis. Cambridge University Press, Cambridge; 1991.

    Chapter  Google Scholar 

  39. Cronin J: Fixed Points and Topological Degree in Nonlinear Analysis. American Mathematical Society, Providence 1964.

    Google Scholar 

Download references

Acknowledgements

This work is supported by the Scientific Research Fund of Sichuan Provincial Education Department under Grant 09ZC057.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kelin Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ZX designed and performed all the steps of proof in this research and also wrote the paper. KL participated in the design of the study and helped to draft and revise manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

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, X., Li, K. Integro-differential inequality and stability of BAM FCNNs with time delays in the leakage terms and distributed delays. J Inequal Appl 2011, 43 (2011). https://doi.org/10.1186/1029-242X-2011-43

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2011-43

Keywords