Skip to main content

Stability and estimate of solution to uncertain neutral delay systems

Abstract

The coefficients and delays in models describing various processes are usually obtained as a results of measurements and can be obtained only approximately. We deal with the question of how to estimate the influence of ‘mistakes’ in coefficients and delays on solutions’ behavior of the delay differential neutral system x i (t) q i (t) x i (t θ i (t))+ j = 1 n ( p i j (t)Δ p i j (t)) x j (t τ i j (t)Δ τ i j (t))= f i (t), i=1,,n, t[0,). This topic is known in the literature as uncertain systems or systems with interval defined coefficients. The goal of this paper is to obtain stability of uncertain systems and to estimate the difference between solutions of a ‘real’ system with uncertain coefficients and/or delays and corresponding ‘model’ system. We develop the so-called Azbelev W-transform, which is a sort of the right regularization allowing researchers to reduce analysis of boundary value problems to study of systems of functional equations in the space of measurable essentially bounded functions. In corresponding cases estimates of norms of auxiliary linear operators (obtained as a result of W-transform) lead researchers to conclusions about existence, uniqueness, positivity and stability of solutions of given boundary value problems. This method works efficiently in the case when a ‘model’ used in W-transform is ‘close’ to a given ‘real’ system. In this paper we choose, as the ‘models’, systems for which we know estimates of the resolvent Cauchy operators. We demonstrate that systems with positive Cauchy matrices present a class of convenient ‘models’. We use the W-transform and other methods of the general theory of functional differential equations. Positivity of the Cauchy operators is studied and then used in the analysis of stability and estimates of solutions.

Results: We propose results about exponential stability of the given system and obtain estimates of difference between the solution of this uncertain system and the ‘model’ system x i (t) q i (t) x i (t θ i (t))+ j = 1 n p i j (t) x j (t τ i j (t))= f i (t), i=1,,n, t[0,). New tests of stability and in the future of existence and uniqueness of boundary value problems for neutral delay systems can be obtained on the basis of this technique.

1 Introduction

In systems of differential equations, modeling real processes, the coefficients and delays are usually known only approximately. For example, models in economics are obtained as a result of corresponding long time observation and regression procedures (see, for example, [1]). The coefficients and delays in model describing various processes in physics are usually obtained as a results of measurements and also can be obtained only approximately [2]. The same situation can be noted in models of population dynamics [3]. A natural question is to estimate the influence of ‘mistakes’ in coefficients and delays on solutions’ behavior of delay differential systems. This topic is known in the literature as uncertain systems or systems with interval defined coefficients.

Important applications present the main reason of attention of many authors to the topic. Systems with uncertain coefficients were considered in the papers [422].

Stability of systems with uncertain coefficients was studied in the papers [68, 1114, 1618, 23]. Stabilization of uncertain systems with unknown delay was studied in [4, 19, 22]. Stability of systems with uncertain delay is studied in [12, 15, 20, 21]. Stability of neutral uncertain systems was considered in [10, 16, 17]. The basic technique of these publications is the use of Lyapunov-Krasovskii functionals.

Consider the system of the neutral functional differential equations

x i ( t ) ( S i x i ) ( t ) + j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t τ i j ( t ) Δ τ i j ( t ) ) = f i ( t ) , i = 1 , , n , t [ 0 , ) ,
(1)

with uncertain coefficients p i j (t) and delays τ i j (t) and the model system

x i (t) ( S i x i ) (t)+ j = 1 n p i j (t) x j ( t τ i j ( t ) ) = f i (t),i=1,,n,t[0,),
(2)

where

x i (ξ)=0,i=1,,n,ξ<0.

We assume here that p i j ,Δ p i j , f i L , where L is the space of essentially bounded functions on [0,+], τ i j and Δ τ i j :[0,+][0,+] are measurable essentially bounded functions, i,j=1,,n. Concerning the operators S i we assume that S i : L L are linear continuous Volterra operators and the spectral radius ρ( S i ) of the operator S i is less than 1 for i=1,,n. Our main purpose is to obtain conditions on the smallness of Δ p i j (t) and Δ τ i j (t) such that the exponential stability of (2) is inherited by (1). Our second aim is to estimate the modulus of the difference of solutions to systems (1) and (2).

The studies of neutral functional differential equations have their own history. Equations in the form

( x ( t ) q ( t ) x ( g ( t ) ) ) + i = 1 m b i (t)x ( h i ( t ) ) =f(t),t[0,),
(3)

were considered in the well-known books [3, 24, 25] (see also the bibliography therein), where existence and uniqueness of solutions and especially stability and oscillation results for these equations were obtained. There exist problems in applications whose models can be written in the form [2]

x (t)q(t) x ( g ( t ) ) + i = 1 m b i (t)x ( h i ( t ) ) =f(t),t[0,).
(4)

This equation is a particular case of (1).

Let us note here that the operator S i : L L in (1) can be, for example, of the following forms:

( S i y ) ( t ) = q i ( t ) y ( g i ( t ) ) , where  g i ( t ) t , t [ 0 , ) , y ( ξ ) = 0 , for  ξ < 0
(5)

or

( S i y)(t)= 0 t k i (t,s)y(s)ds,t[0,+),
(6)

where q j (t) are essentially bounded measurable functions, τ j (t) are measurable functions for j=1,,m, and k i (t,s) are summable with respect to s and measurable essentially bounded with respect to t for i=1,,n. All linear combinations of the operators (5) and (6) and their superpositions are also allowed.

The study of the neutral functional differential equations is essentially based on the questions of the action and estimates of the spectral radii of the operators in the spaces of discontinuous functions, for example, in the spaces of summable or essentially bounded functions. The operator (5) is a key object in this topic. Properties of this operator were studied, for example, in [26, 27]. In order to achieve the action of the operator (5) in the space of essentially bounded functions L , we have to assume that mes{t: g i (t)=c}=0 for every constant c. Let us suppose everywhere below that this condition is fulfilled. It is known that the spectral radius of the integral operator (6), considered on every finite interval t[0,ω], is equal to 0 (see, for example, [1]) and the inequality

sup t 0 0 t | k i (t,s)|ds<1,t[0,),i=1,,n,

implies that the spectral radii of the integral operators S i defined on the semiaxis are less than 1. Concerning the operator (5), we can note that the inequality

esssup t 0 | q i (t)|<1,t[0,),i=1,,n,

is a sufficient condition that the spectral radii ρ( S i ) of the operators S i is less than 1. Below we assume that this inequality is fulfilled.

Various results on existence and uniqueness of solutions to boundary value problems for (4) and its stability were obtained in [1], where also the basic results about the representation of solutions were presented. Note also in this connection the papers [2831], where results on nonoscillation and positivity of Green’s functions for neutral functional differential equations were obtained.

The goal of this paper is to obtain stability of uncertain systems and to estimate the difference between solutions of a ‘real’ system with uncertain coefficients and/or delays and corresponding ‘model’ system with fixed coefficients and delays. Instead of the traditional Lyapunov functionals technique, we propose an approach based on the idea of Azbelev’s W-transform presented in the book [1]. This transform is a corresponding right regularization allowing researchers to reduce the analysis of boundary value problems to the study of systems of functional equations in the space of measurable essentially bounded functions. In corresponding cases estimates of norms of auxiliary linear operators (obtained as a result of W-transform) lead researchers to conclusions about existence, uniqueness, positivity and stability of solutions of given boundary value problems. This method works efficiently in the case when a ‘model’ used in W-transform is ‘close’ to a given ‘real’ uncertain system. In this paper we choose as a ‘model’, system for which we know solutions or estimates of the resolvent (Green’s) operators. In this paper we demonstrate that systems with positive Cauchy operators [32, 33] represent a class of very convenient ‘models’.

It should be noted that, although results are below formulated for the systems with neutral terms in the form (5), we can, using the approach of the papers [30, 31], obtain corresponding results with linear combinations and superpositions of the operators defined by (5) and (6).

2 Stability and estimation of solution of neutral systems with uncertain coefficients

Consider two systems

x i (t) q i (t) x i ( g i ( t ) ) + j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t τ i j ( t ) ) = f i (t),t[0,),
(7)
x i (t) q i (t) x i ( g i ( t ) ) + j = 1 n p i j (t) x j ( t τ i j ( t ) ) = z i (t),t[0,),
(8)

where

x i (ξ)= x i (ξ)=0,i=1,,n,ξ<0.

Here p i j ,Δ p i j , f i L , and L is the space of essentially bounded functions on [0,]. τ i j :[0,][0,] are measurable essentially bounded functions for i,j=1,,n.

It is known [1] that the general solution of (8) has the representation

x i (t)= 0 t j = 1 n W i j (t,s) z j (s)ds+ j = 1 n W i j 0 (t,0) x j (0),i=1,,n,
(9)

where W(t,s)= { W i j ( t , s ) } i , j = 1 n is the Cauchy matrix of system (8) and W 0 (t,0)={ W i j 0 (t,0)} is the fundamental matrix of system (8) satisfying the condition W 0 (0,0)=I.

Substituting representation (9) into (7), we get

z i ( t ) j = 1 n Δ p i j ( t ) σ ( t τ i j ( t ) , 0 ) 0 t τ i j ( t ) k = 1 n W j k ( t τ i j ( t ) , s ) z k ( s ) d s = j = 1 n Δ p i j ( t ) σ ( t τ i j ( t ) , 0 ) k = 1 n W j k 0 ( t τ i j ( t ) , 0 ) x k ( 0 ) + f i ( t ) ,
(10)

where σ(t,s)={ 1 , t s , 0 , t < s and W j k (ξ,s)=0 if ξ<s.

We can write (10) in the vector form

z(t)=Ωz(t)+ f (t),
(11)

where the operator Ω: L n L n ( L n is the space of n vector-functions with measurable essentially bounded components) is defined as follows:

Ωz(t)= { j = 1 n Δ p i j ( t ) σ ( t τ i j ( t ) , 0 ) 0 t τ i j ( t ) k = 1 n W j k ( t τ i j ( t ) , s ) z j ( s ) d s } i = 1 n ,
(12)

and

f (t)= { f i ( t ) + j = 1 n Δ p i j ( t ) σ ( t τ i j ( t ) , 0 ) k = 1 n W j k 0 ( t τ i j ( t ) , 0 ) x k ( 0 ) } i = 1 n .

Estimating the norm of the operator Ω: L n L n , we see that

Ω max 1 i n esssup t 0 j = 1 n |Δ p i j (t)| esssup t 0 0 t k = 1 n | W j k (t,s)|ds.
(13)

Estimates of esssup t 0 0 t | W i j (t,s)|ds is the key problem in our approach. Now describe the cases in which the Cauchy matrix W(t,s) and the fundamental matrix W 0 (t,0) of system (8) can be estimated.

Let us prove several assertions for the neutral system

x i ( t ) q i ( t ) x ( g i ( t ) ) + j = 1 n p i j ( t ) x j ( t τ i j ( t ) ) = f i ( t ) , t [ 0 , ) , i = 1 , , n , x i ( ξ ) = x i ( ξ ) = 0 , i = 1 , , n , ξ < 0 .
(14)

Theorem 1 Let p i i 0, q i 0,

h i ( t ) t p i i (s)ds 1 e ln e 1 + γ i ( t ) q i ( t ) ,t[0,),
(15)

where

h i (t)=min { t τ i i ( t ) , g i ( t ) } ,
(16)
γ i (t)= { p i i ( g i ( t ) ) p i i ( t ) if  g i ( t ) 0 , 0 if  g i ( t ) < 0 ,
(17)

p i j 0, ij, i,j=1,,n, then all elements W i j (t,s) of the Cauchy matrix and all elements W i j 0 (t,0) of the fundamental matrix of system (14) are nonnegative.

Proof According to the result of [31], the conditions of Theorem 1 imply the positivity of the Cauchy functions and fundamental solutions of the diagonal equations

x i (t) q i (t) x i ( g i ( t ) ) + p i i (t) x i ( t τ i i ( t ) ) = f i (t),t[0,),i=1,,n.
(18)

Let us consider the system of diagonal equations denoting by X(t,s)= { X i j ( t , s ) } i , j = 1 n its Cauchy and X 0 (t,0)= { X i j 0 ( t , 0 ) } i , j = 1 n fundamental matrices, respectively. These matrices are the diagonal ones, i.e. X i j (t,s)0 and X i j 0 (t,0)0 for ji, i,j=1,,n, and the elements X i i (t,s) and X i i 0 (t,0) coincide with the Cauchy functions and the fundamental solutions of the scalar diagonal equations, respectively. We can write the system

x i ( t ) = 0 t X i i ( t , s ) j = 1 , j i n p i j ( s ) x j ( s τ i j ( s ) ) d s + 0 t X i i ( t , s ) f i ( s ) d s + X i i 0 ( t , 0 ) x i ( 0 ) , t [ 0 , ) , i = 1 , , n , x i ( ξ ) = 0 , i = 1 , , n , ξ < 0 ,
(19)

which is equivalent to system (14). On every finite interval [0,ω] the spectral radius of the integral operator

(Kx)(t)= { 0 t X i i ( t , s ) j = 1 , j i n p i j ( s ) x j ( s τ i j ( s ) ) d s } i , j = 1 n ,t[0,ω],
(20)

where

x i (ξ)=0,i=1,,n,ξ<0,

is 0. Denoting

ϕ(t)= 0 t X i i (t,s) f i (s)ds+ X i i 0 (t,0) x i (0),t[0,),i=1,,n,
(21)

we can write (19) as x(t)=( ( I K ) 1 ϕ)(t)=ϕ(t)+(Kϕ)(t)+( K 2 ϕ)(t)+ and consequently every component of solution vector x(t)=col{ x 1 (t),, x n (t)} are nonnegative for nonnegative vectors f=col{ f 1 ,, f n } and x(0)=col{ x 1 (0),, x n (0)}. As ω tends to infinity we complete the proof of Theorem 1. □

Definition 1 The Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,s) satisfies the exponential estimates, if there exist real numbers N,α>0 such that | W i j (t,s)|N e α ( t s ) and | W i j 0 (t,0)|N e α t for 0st<, i,j=1,,n.

Theorem 2 Let the conditions of Theorem  1 be fulfilled for system (14) and there exist y i j 0 for i,j=1,,n, such that j = 1 n y i j >0, and the system of inequalities

k = 1 n p i k (t) y k j δ i j ,t[0,),i=1,,n,
(22)

where δ i j ={ 1 , i = j , 0 , i j , be satisfied.

Then the Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates and the integral estimate

0 t W i j (t,s)ds y i j for i,j=1,,n.
(23)

Proof Let us extend the coefficients and delays in system (14) on the interval [H,0), where H= max i , j = 1 , , n esssup t 0 τ i j (t) as follows: p i i (t)= esssup t 0 p i i (t), p i j (t)0 for ij, τ i j (t)0, q i (t)0 for i,j=1,,n, t[H,0), and consider system (14) on [H,), calling it the extent system. The constant vector-function Y j (t) { y i j } i = 1 n satisfies the inequality x i (t) q i (t) x ( g i (t))+ j = 1 n p i j (t) x j (t τ i j (t)) δ i j , t[H,), i=1,,n. All entries of the Cauchy matrix of the extent system are nonnegative according to Theorem 1. The Cauchy matrices of system (14) and the extent system coincide in 0st<. Using nonnegativity of all entries of the Cauchy matrix W(t,s) of system (14) we get the inequalities H t W i j (t,s)ds y i j for i,j=1,,n, t[H,) that implies inequalities (23). For every bounded right hand side, the solution of system (14) is bounded. From this the exponential estimates of entries of the Cauchy matrix W(t,s) and the fundamental matrix W 0 (t,s) follow (see [1], Paragraph 5.3). □

Analogously we can obtain the following assertion.

Theorem 3 Let the conditions of Theorem  1 be fulfilled for system (14), let there exist a vector γ=col{ γ 1 ,, γ n } such that γ i >0 for i=1,,n, and

k = 1 n p i j (t) γ j 1,t[0,),i=1,,n.
(24)

Then the Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates and the integral estimate

0 t j = 1 n W i j (t,s)ds γ i for i=1,,n.
(25)

Theorem 4 Let all coefficients p i j be constants and the conditions of Theorem  1 be fulfilled for system (14). Then the Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates if and only if all components of the solution vector γ=col{ γ 1 ,, γ n } of the algebraic system

k = 1 n p i j γ j =1,i=1,,n,
(26)

are positive.

Proof Sufficiency follows from Theorem 3. To prove necessity we extend the coefficients and delays in system (14) on the interval [H,0), where H= max i , j = 1 , , n esssup t 0 τ i j (t) such that τ i j (t)0, q i (t)0, p i j (t)= p i j , f i (t)1 for t[H,0), i,j=1,,n. The constant vector-function γ=col{ γ 1 ,, γ n } is now the solution of the system

x i (t) q i (t) x i ( g i ( t ) ) + j = 1 n p i j x j ( t τ i j ( t ) ) =1,t[H,),
(27)

satisfying the condition x i (H)= γ i , i=1,,n. Using the formula of the solution’s representation (see (9)) on [H,), we obtain

γ i = H t j = 1 n W i j (t,s)ds+ j = 1 n W i j 0 (t,H) γ j ,i=1,,n.
(28)

All elements W i j (t,s) of the Cauchy matrix of system (27) are nonnegative. From the fact that W i j 0 (t,0) satisfies the exponential estimate we obtain γ i >0 for i=1,,n. □

Remark 1 Under the conditions of Theorem 4, the equalities

lim t 0 t j = 1 n W i j (t,s)ds= γ i for i=1,,n,
(29)
lim t 0 t W i j (t,s)ds= y i j for i,j=1,,n,
(30)

are true.

Theorem 5 Let q i 0, the inequality (15) be fulfilled and let there exist a vector γ=col{ γ 1 ,, γ n } with all positive components such that

p i i (t) γ i j = 1 , i j n | p i j (t)| γ j 1,t[0,),i=1,,n.
(31)

Then the Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates and the integral estimates

0 t j = 1 n | W i j (t,s)|ds γ i for i=1,,n,
(32)
0 t | W i j (t,s)|ds y i j for i,j=1,,n,
(33)

where the constants γ i and y i j are defined by inequalities (31) and (22) (where there are | p i k (t)| instead of p i k (t)), respectively.

Proof Consider the system

x i ( t ) q i ( t ) x i ( g i ( t ) ) + p i i ( t ) x i ( t τ i i ( t ) ) j = 1 , i j n | p i j ( t ) | x j ( t τ i j ( t ) ) = f i ( t ) , t [ 0 , ) ,
(34)

where

x i (ξ)= x i (ξ)=0,i=1,,n,ξ<0.

We can write the equivalent integral system

x i ( t ) = 0 t X i i ( t , s ) j = 1 , j i n | p i j ( s ) | ( s ) x j ( s τ i j ( s ) ) d s + 0 t X i i ( t , s ) f i ( s ) d s + X i i 0 ( t , 0 ) x i ( 0 ) , t [ 0 , ) , i = 1 , , n .
(35)

The operator

( | K | x ) (t)= { 0 t X i i ( t , s ) j = 1 , j i n | p i j ( s ) | x j ( s τ i j ( s ) ) d s } i , j = 1 n ,t[0,),
(36)

where

x i (ξ)=0,i=1,,n,ξ<0,

is positive. Condition (31) implies according to Theorem 5.6 [34] that its spectral radius ρ(|K|)<1. According to Theorem 5.3 [34] we obtain the inequality ρ(K)<ρ(|K|) for the spectral radius of the operator K defined by equality (20) for t[0,). The exponential estimates of X i i (t,s) and X i i 0 (t,0) imply that

φ i (t) 0 t X i i (t,s) f i (s)ds+ X i i 0 (t,0) x i (0),t[0,),i=1,,n,
(37)

is bounded for every bounded right hand sides f(t)=col{ f 1 (t),, f n (t)}. Comparing the operators K and |K| we see that ( I | K | ) 1 φ = (I+|K|+| K 2 |+)φ ( I K ) 1 φ = (I+K+ K 2 +)φ(I+|K|+| K 2 |+)φ = ( I | K | ) 1 φ and consequently the solution vector x(t)=col{ x 1 (t),, x n (t)} is bounded for every bounded φ=col{ φ 1 (t),, φ n (t)}. These inequalities prove the estimates (32) and (33). According to Theorem 3.5 [1], the Cauchy and fundamental matrices satisfy the exponential estimates. □

Theorem 6 Let the conditions of Theorem  5 be fulfilled and the norm of the operator Ω, defined by equality (12), satisfy the inequality Ω<1, then the Cauchy matrix W(t,s) of system (14) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates and the difference Δx of solution x of system (7) and x ¯ of system (14) satisfies the inequality

Δx max 1 i n f 1 Ω γ i ,i=1,,n,
(38)

for sufficiently large t.

In order to prove Theorem 6 we can note that in the case Ω<1 we get z 1 1 Ω f. A reference to (9) completes the proof.

3 Stability and estimation of solution of neutral systems with uncertain delays

Consider the system

x i (t) q i (t) x i ( g i ( t ) ) + j = 1 n p i j (t) x j ( t θ i j ( t ) ) = f i (t).
(39)

Rewrite (39) in the form

x i ( t ) q i ( t ) x i ( g i ( t ) ) + j = 1 n p i j ( t ) x j ( t τ i j ( t ) ) + j = 1 n p i j ( t ) ( x j ( t θ i j ( t ) ) x j ( t τ i j ( t ) ) ) = f i ( t ) .
(40)

We have

z i (t)+ j = 1 n p i j (t) ( x j ( t θ i j ( t ) ) x j ( t τ i j ( t ) ) ) = f i (t),
(41)

where

x i (t) q i (t) x i ( g i ( t ) ) + j = 1 n p i j (t) x j ( t τ i j ( t ) ) = z i (t).
(42)

It is clear that the solution of (42) has the representation

x i (t)= 0 t j = 1 n W i j (t,s) z j (s)ds+ j = 1 n W i j 0 (t,0) x j (0),i=1,,n,
(43)

where W(t,s)={ W i j (t,s)}, i,j=1,,n is the Cauchy matrix of system (42) and W i j 0 (t,0) is its fundamental matrix. Let us assume for simplicity that t τ i j (t)0 and t θ i j (t)0 for t0. If this is not fulfilled we can extend coefficients on the interval [H,0) and consider the system on [H,) as in the Section 2.

It is clear that x j (t Θ i j (t)) x j (t τ i j (t))= t Θ i j ( t ) t τ i j ( t ) x j (s)ds and

z i (t)+ j = 1 n p i j (t) t Θ i j ( t ) t τ i j ( t ) x j (s)ds= f i (t).
(44)

From (42) we get x i (t)= z i (t) j = 1 n p i j (t) x j (t τ i j (t)). Substituting this into (43), we obtain

z i (t)= j = 1 n p i j (t) t Θ i j ( t ) t τ i j ( t ) ( I S j ) 1 ( k = 1 n p j k ( s ) x j ( s τ j k ( s ) ) z j ( s ) ) ds+ f i (t),
(45)

where the operator S j is defined by (5), and after substituting x from (43) we get

z i ( t ) = j = 1 n p i j ( t ) t Θ i j ( t ) t τ i j ( t ) ( I S j ) 1 ( k = 1 n p j k ( s ) ( 0 s τ j k ( s ) m = 1 n W k m ( s , ξ ) z m ( ξ ) d ξ + m = 1 n W k m 0 ( s τ k m ( s ) , 0 ) x m ( 0 ) ) z j ( s ) ) d s + f i ( t ) .
(46)

Here the operator Ω: L n L n , is defined as Ω=col{ Ω 1 ,, Ω n }, where every component Ω i is defined as follows:

( Ω i z i ) ( t ) = j = 1 n p i j ( t ) t Θ i j ( t ) t τ i j ( t ) ( I S j ) 1 × ( k = 1 n p j k ( s ) ( 0 s τ i k ( s ) m = 1 n W k m ( s , ξ ) z m ( ξ ) d ξ ) d s z j ( s ) ) d s .
(47)

We can write (46) in the vector form

z(t)=Ωz(t)+f(t).
(48)

It is clear that

Ω max 1 1 n esssup t 0 j = 1 n | p i j ( t ) | 1 1 q esssup t 0 { k = 1 n p j k ( t ) γ k + 1 } × esssup t 0 | θ i j ( t ) τ i j ( t ) | ,
(49)

where q= max 1 i n esssup t 0 q i (t) and the components γ i of the vector γ=col{ γ 1 ,, γ n } satisfies the system of inequalities (31).

Theorem 7 Let the conditions of Theorem  5 be fulfilled and Ω<1, then the Cauchy matrix W(t,s) of system (39) and its fundamental matrix W 0 (t,0) satisfy the exponential estimates and the difference Δx of solution x of system (39) and x ¯ of system (14) satisfies the estimate

Δx max 1 i n f 1 Ω γ i ,i=1,,n,
(50)

for sufficiently large t.

4 Stability and estimation of solution of neutral systems with uncertain coefficients and delays

Consider the system

x i (t) q i (t) x ( g i ( t ) ) + j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t θ i j ( t ) ) = f i (t)
(51)

and the auxiliary system

x i (t) q i (t) x ( g i ( t ) ) + j = 1 n p i j (t) x j ( t τ i j ( t ) ) = z i (t).
(52)

The solution of (52) has the representation

x i (t)= 0 t j = 1 n W i j (t,s) z j (s)ds+ j = 1 n W i j 0 (t,0) x j (0),i=1,,n.
(53)

Rewrite (51) in the form

x i ( t ) q i ( t ) x ( g i ( t ) ) + j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t τ i j ( t ) ) j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t τ i j ( t ) ) + j = 1 n ( p i j ( t ) Δ p i j ( t ) ) x j ( t θ i j ( t ) ) = f i ( t ) .
(54)

Using the technique of two previous sections, we obtain the system of the following form:

z(t)=Ωz(t)+ f (t),

where f(t) f (t)0 when t, and estimating the norm of the corresponding operator Ω, we get

Ω max 1 i n { esssup t 0 j = 1 n | Δ p i j ( t ) | γ j } + max 1 i n { esssup t 0 j = 1 n ( | p i j ( t ) | + | Δ p i j ( t ) | ) } 1 1 q esssup t 0 | θ i j ( t ) τ i j ( t ) | × esssup t 0 ( k = 1 n | p j k ( t ) | γ k + 1 ) ,
(55)

where γ i , i=1,,n, are defined by inequalities (31).

Theorem 8 Let the conditions of Theorem  5 be fulfilled and Ω<1, then the Cauchy matrix C(t,s) of system (51) and its fundamental matrix C 0 (t,0) satisfy the exponential estimates and the difference Δx of solution x of system (51) and x ¯ of system (14) satisfies the inequality

Δx max 1 i n f 1 Ω γ i ,i=1,,n,

for sufficiently large t.

References

  1. Azbelev NV, Maksimov VP, Rakhmatullina LF Advanced Series in Mathematical Science and Engineering 3. In Introduction to the Theory of Linear Functional-Differential Equations. World Federation Publishers Company, Atlanta; 1995.

    Google Scholar 

  2. Kolmanovskii V, Myshkis AD: Introduction to Theory and Applications of Functional Differential Equations. Kluwer Academic, Dordrecht; 1999.

    Book  Google Scholar 

  3. Gopalsamy K: Stability and Oscillation in Delay Differential Equations of Population Dynamics. Kluwer Academic, Dordrecht; 1992.

    Book  Google Scholar 

  4. Cao Y, Sun Y, Cheng C: Delay-dependent robust stabilization of uncertain systems with multiple state delays. IEEE Trans. Autom. Control 1998, 43: 1608-1612. 10.1109/9.728880

    Article  MathSciNet  Google Scholar 

  5. Chen XW, Liu B: Existence and uniqueness theorem for uncertain differential equations. Fuzzy Optim. Decis. Mak. 2010, 9: 69-81. 10.1007/s10700-010-9073-2

    Article  MathSciNet  Google Scholar 

  6. Dong Y, Liu J: Exponential stabilization of uncertain nonlinear time-delay systems. Adv. Differ. Equ. 2012., 2012: Article ID 180 10.1186/1687-1847-2012-180

    Google Scholar 

  7. Dong Y, Mei S, Wang X: Novel stability criteria of nonlinear uncertain systems with time-varying delay. Abstr. Appl. Anal. 2011., 2011: Article ID 969674

    Google Scholar 

  8. Dzhalladova IA, Bastinec J, Diblik J, Khusainov DY: Estimates of exponential stability for solutions of stochastic control systems with delay. Abstr. Appl. Anal. 2011., 2011: Article ID 920412

    Google Scholar 

  9. Yao K, Gao JW, Gao Y: Some stability theorems of uncertain differential equation. Fuzzy Optim. Decis. Mak. 2013, 12: 3-13. 10.1007/s10700-012-9139-4

    Article  MathSciNet  Google Scholar 

  10. Khusainov D, Diblik J, Bavstinec J, Ryvolova J: Interval stability of linear systems neutral type. Vyčisl. Prikl. Mat. (Kiev) 2011, 3(106):148-160. (in Russian)

    Google Scholar 

  11. Kwon OM, Park JH: Exponential stability of uncertain dynamic systems including state delay. Appl. Math. Lett. 2006, 19: 901-907. 10.1016/j.aml.2005.10.017

    Article  MathSciNet  Google Scholar 

  12. Kwon OM, Park JH: Exponential stability analysis for uncertain neural networks with interval time-varying delays. Appl. Math. Comput. 2009, 212: 530-541. 10.1016/j.amc.2009.02.043

    Article  MathSciNet  Google Scholar 

  13. Phat VN, Nam PT: Exponential stability and stabilization of uncertain linear time-varying systems using parameter dependent Lyapunov function. Int. J. Control 2007, 80: 1333-1341. 10.1080/00207170701338867

    Article  MathSciNet  Google Scholar 

  14. Phat VN, Niamsup P: Stability analysis for a class of functional differential equations and applications. Nonlinear Anal., Theory Methods Appl. 2009, 71: 6265-6275. 10.1016/j.na.2009.06.028

    Article  MathSciNet  Google Scholar 

  15. Shao H-Y: New delay-dependent stability criteria for systems with interval delay. Automatica 2009, 45: 744-749. 10.1016/j.automatica.2008.09.010

    Article  Google Scholar 

  16. Shatyrko AV, Khusainov DY, Diblik J, Bastinec J, Rivolova A: Estimates of perturbation of nonlinear indirect interval control system of neutral type. J. Autom. Inf. Sci. 2011, 43: 13-28. 10.1615/JAutomatInfScien.v43.i1.20

    Article  Google Scholar 

  17. Shatyrko AV, Khusainov DY: Absolute interval stability of indirect regulating systems of neutral type. J. Autom. Inf. Sci. 2010, 42: 43-54.

    Article  Google Scholar 

  18. Shatyrko AV, Khusainov D: Interval stability investigations of different regulator systems with time delay by Lyapunov-Krasovskii functional. Bull. Taras Shevchenko Natl. Univ. Kyiv. Ser., Phys. Math. 2009, 3: 212-221. (in Ukrainian)

    Google Scholar 

  19. Sun YJ: Global stabilizability of uncertain systems with time-varying delays via dynamic observer-based output feedback. Linear Algebra Appl. 2002, 353: 91-105. 10.1016/S0024-3795(02)00292-6

    Article  MathSciNet  Google Scholar 

  20. Tian J, Zhou X: Improved asymptotic stability criteria for neural networks with interval time-varying delay. Expert Syst. Appl. 2010, 37: 7521-7525. 10.1016/j.eswa.2010.04.093

    Article  Google Scholar 

  21. Verriest EI, Florchinger P: Stability of stochastic systems with uncertain time delays. Syst. Control Lett. 1995, 24: 41-47. 10.1016/0167-6911(94)00030-Y

    Article  MathSciNet  Google Scholar 

  22. Yue D: Robust stabilization of uncertain systems with unknown input delay. Automatica 2004, 40: 331-336. 10.1016/j.automatica.2003.10.005

    Article  Google Scholar 

  23. Wang TY, Xie LH, De Souz CE: Robust control of a class of uncertain nonlinear systems. Syst. Control Lett. 1992, 19: 139-149. 10.1016/0167-6911(92)90097-C

    Article  Google Scholar 

  24. Erbe LN, Kong Q, Zhang BG: Oscillation Theory for Functional Differential Equations. Dekker, New York; 1995.

    Google Scholar 

  25. Gyori I, Ladas G: Oscillation Theory of Delay Differential Equations. Clarendon, Oxford; 1991.

    Google Scholar 

  26. Drakhlin ME: Operator of the internal superposition in the space of summable functions. Izv. Vysš. Učebn. Zaved., Mat. 1986, 5: 18-23.

    MathSciNet  Google Scholar 

  27. Drakhlin ME, Plyshevskaya TK: To the theory of functional differential equations. Differ. Uravn. 1978, 14: 1347-1361.

    MathSciNet  Google Scholar 

  28. Berezansky L, Braverman E: Non-oscillation properties of a linear neutral differential equations. Funct. Differ. Equ. 2002, 9(3-4):275-288.

    MathSciNet  Google Scholar 

  29. Berezansky L, Braverman E: Oscillation criteria for linear neutral differential equations. J. Math. Anal. Appl. 2003, 286: 601-617. 10.1016/S0022-247X(03)00502-X

    Article  MathSciNet  Google Scholar 

  30. Domoshnitsky A: Extension of Sturm’s theorem to equations with time-lag. Differ. Uravn. 1983, 19: 1475-1482.

    MathSciNet  Google Scholar 

  31. Gusarenko SA, Domoshnitskii AI: Asymptotic and oscillation properties of the first order scalar functional-differential equations. Differ. Uravn. 1989, 25: 1480-1491.

    MathSciNet  Google Scholar 

  32. Agarwal RP, Berezansky L, Braverman E, Domoshnitsky A: Nonoscillation Theory of Functional Differential Equations with Applications. Springer, New York; 2012.

    Book  Google Scholar 

  33. Domoshnitsky A, Sheina MV: Nonnegativity of Cauchy matrix and stability of systems with delay. Differ. Uravn. 1989, 25: 201-208.

    MathSciNet  Google Scholar 

  34. Krasnosel’skii MA, Vainikko GM, Zabreiko PP, Rutitskii JB, Stezenko VJ: Approximate Methods for Solving Operator Equations. Nauka, Moscow; 1969. (in Russian)

    Google Scholar 

Download references

Acknowledgements

This paper was prepared as an open auxiliary theoretical result in the frame of the project ‘Exploitation of a synergetic model for development of business of innovation type’ on the topic No. 2013/276-C ‘Development of a model of operation of innovation business as dynamic model with memory effect’, supported by Perm National Research Polytechnic University with financial support of Ministry of Science and Education of Russian Federation (agreement No. 02.G25.31.0068 of 23.05.2013).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander Domoshnitsky.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors worked on this paper together and have the same rights on the content of the paper. All authors read and approved the final manuscript.

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

Domoshnitsky, A., Gitman, M. & Shklyar, R. Stability and estimate of solution to uncertain neutral delay systems. Bound Value Probl 2014, 55 (2014). https://doi.org/10.1186/1687-2770-2014-55

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-55

Keywords