Skip to main content

Nonoscillation, maximum principles, and exponential stability of second order delay differential equations without damping term

Abstract

Delays, arising in nonoscillatory and stable ordinary differential equations, can induce oscillation and instability of their solutions. That is why the traditional direction in the study of nonoscillation and stability of delay equations is to establish a smallness of delay, allowing delay differential equations to preserve these convenient properties of ordinary differential equations with the same coefficients. In this paper, we find cases in which delays, arising in oscillatory and asymptotically unstable ordinary differential equations, induce nonoscillation and stability of delay equations. We demonstrate that, although the ordinary differential equation x (t)+c(t)x(t)=0 can be oscillating and asymptotically unstable, the delay equation x (t)+a(t)x(th(t))b(t)x(tg(t))=0, where c(t)=a(t)b(t), can be nonoscillating and exponentially stable. Results on nonoscillation and exponential stability of delay differential equations are obtained. On the basis of these results on nonoscillation and stability, the new possibilities of non-invasive (non-evasive) control, which allow us to stabilize a motion of single mass point, are proposed. Stabilization of this sort, according to common belief, requires a damping term in the second order differential equation. Results obtained in this paper refute this delusion.

MSC:34K20.

1 Introduction

Let us start with one of motivations of this study considering a simplified model for motion of a single mass point

X (t)=f(t),t[0,+),
(1.1)

where X(t)=col{ x 1 (t), x 2 (t), x 3 (t)}, f(t)=col{ f 1 (t), f 2 (t), f 3 (t)}. Let a solution Y(t)=col{ y 1 (t), y 2 (t), y 3 (t)}, which defines the trajectory we want to ‘hold’, be known and our purpose is to hold our object ‘close’ to this trajectory. It is impossible to ‘achieve’ this proximity on the semiaxis, because of instability of system (1.1). As a result, we have to make corrections of the motion changing the right-hand side f. We try to construct a control which makes this correction automatically. A standard approach is to construct a feedback control u(t) which depends on the difference X(t)Y(t) or more exactly on the difference X(tτ(t))Y(tτ(t)) since in real systems delay in receiving signal and in reaction on this signal arises. The control terms like

u(t)= i = 1 m P i (t) { X ( t τ i ( t ) ) Y ( t τ i ( t ) ) } ,t[0,+)
(1.2)

(here P i (t) are 3×3 matrices, i=1,,m) are usually called non-invasive (non-evasive by other authors). Note that in laser physics the term ‘Pyragas control’ is widely used. Adding this control u(t) into the equation of motion, we get

X (t)+ i = 1 m P i (t)X ( t τ i ( t ) ) =g(t),t[0,+),
(1.3)

where

g(t)=f(t)+ i = 1 m P i (t)Y ( t τ i ( t ) ) ,t[0,+),
(1.4)

is the known right-hand side. Using the fact that Y(t) satisfies (1.1), i.e. Y (t)=f(t), and subtracting this equality from (1.3), we get the following system:

Z (t)+ i = 1 m P i (t)Z ( t τ i ( t ) ) =0,t[0,+)
(1.5)

for the difference Z(t)=X(t)Y(t).

The problem of stabilization has been reduced to the exponential stability of this system. In the simplest case of diagonal matrices P i , we get three separate equations. In the case of absence of delays ( τ i (t)0, t[0,+), i=1,,m) and constant diagonal entries, we get existence of positive roots of characteristic equation in the case of negative entries and pure imaginary roots in the case of positive ones. Thus, even in the case of the positive entries, each of these equations is Lyapunov stable, but not exponentially stable and, as a result, there is no stability with respect to the right-hand side (in the sense that small changes in the right-hand side could imply very essential changes in the solutions). In this paper we limit ourselves to diagonal matrices P i (t) (i=1,,m) in order to be concentrated on the scalar second order equation

x (t)+px(tτ)=0,t[0,+),
(1.6)

and its natural generalizations with variable coefficients and delays

x (t)+p(t)x ( t τ ( t ) ) =0,t[0,+)
(1.7)

and

x (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =0,t[0,+),
(1.8)

where

x(ξ)=0for ξ<0,
(1.9)

p i and τ i are measurable essentially bounded functions (below we write this as p i L ) (i=1,,m).

We understand the solution of (1.7) in the traditional sense like in the book [1], i.e. as a function x:[0,+)(,+) with absolutely continuous derivative x on every finite interval such that x L , satisfying (1.7).

Note that the delay equation (1.7) usually inherits the oscillation properties of the corresponding ODE

x (t)+p(t)x(t)=0.
(1.10)

For example, it was proven by Brands [2] that for every nonnegative p(t) and bounded delay τ(t), (1.7) is oscillatory if and only if corresponding ordinary differential equation (1.10) is oscillatory. The asymptotic behavior of ODE (1.10) is not inherited by (1.7): Hale in his well-known book [[3], Chapter 5, Section 5.4], discussed the characteristic equation

λ 2 +p e τ λ =0
(1.11)

for (1.6) and demonstrated that there exist at least two roots with positive real parts in the case pτ> π 2 . Myshkis, analyzing the roots of the characteristic equation (1.11), proved instability of equation (1.6) for each couple of positive constants p and τ (see [[4], Chapter III, Section 16, pp.105-106]). Note that the degree of instability of (1.6) (the maximum among positive real parts of the roots to (1.11)) tends to zero when τ0 and τ+ [4]. The results on existence of unbounded solutions in the case of variable coefficients and delays were obtained on the basis of a growth of the Wronskians in [5] and developed in [6, 7]. One of the results can be formulated as follows [6]: if there exists a positive constant ε such that τ(t)>ε and p(t)>ε, then there exist unbounded solutions to (1.7). In the paper [8] results on the instability of the system

X (t)+P(t)X(tω)=0,t(,+),
(1.12)

where ω>0 and P(t) is an ω-periodic symmetric matrix function, were obtained.

Results on the boundedness of solutions to the delay equation (1.7) in the case of nondecreasing coefficient p(t) and 0 p α (t)τ(t)dt< for a corresponding α were obtained in the paper [9]. The asymptotic formulas of solutions to the second order equation (1.6) are presented in [[4], Chapter III, Section 16] and to (1.7) with a delay τ(t) summable on the semiaxis in [10, 11]. In [6] it was found that all solutions of (1.7) with positive nondecreasing and bounded coefficient p(t) and nondecreasing h(t)tτ(t) are bounded if and only if

0 τ(t)dt<.
(1.13)

Note that an approach for studying asymptotic properties of equations with linear transformations of arguments

x (t)= j = l , j 0 l a j x ( q j t ) +λx(t),t(,+),
(1.14)

where a j , q and λ are constants, was proposed in [12] and developed in [13], where the equations with a combination of delayed and advanced arguments are considered. A study of advanced equations (see, for example, (1.7), where τ(t)0) can be found in the paper [14] in which results on boundedness, stability, and asymptotic representations of solutions are obtained.

Summarizing, we can conclude that it appears to be impossible to obtain results as regards the exponential stability for one term of (1.7). It means that the control of the form

u(t)=P(t) { X ( t τ ( t ) ) Y ( t τ ( t ) ) } ,
(1.15)

without damping term cannot help us in the stabilization of system (1.1). Stability and instability of the second order delay differential equation

x (t)+ q 1 x (t)+ q 2 x (tτ)+ p 1 x(t)+ p 2 x(tτ)=0
(1.16)

with constant coefficients was studied in [15, 16], where it was noted that this equation in the case p 1 p 2 <0 is of interest in machine tool analysis, in biology in explaining self-balancing of the human body and in robotics in constructing biped robots (see the bibliography therein). Results on the stability of the equation

x (t)+q(t) x (t)+p(t)x ( t τ ( t ) ) =0,q>0,p>0,
(1.17)

were obtained in [17, 18]. First results on the exponential stability of the equation

x (t)+q(t) x ( t θ ( t ) ) +p(t)x ( t τ ( t ) ) =0,q>0,p>0,
(1.18)

without the assumption θ(t)0, as far as we know, were obtained in [19] and developed in [20, 21]. The asymptotic stability of the equation

x (t)+ p 1 x(t)+ p 2 x(tτ)=0,
(1.19)

with constant delay τ and coefficients p 1 and p 2 and without damping term were obtained in the paper [15]. These results were based on Pontryagin’s technique for analysis of the roots of quasi-polynomials [22] and could not be used in the case of (1.8) with variable coefficients and/or delays. Note also that the assumption as regards the absence of the delay in the first term does not allow us to use this result of [15] in the stabilization in the model of the motion described in (1.1) by the delay control (1.2). Let us try to imagine situations in which variable delays and coefficients arising in the delayed feedback control may be important: (1) the case of control for missiles, where the delay depends on their distance from the controller and is variable; (2) spending of fuel implies the change of the mass of the missiles, which leads to variable coefficients in the delay system. As a conclusion, it should be stressed that there are no results for the exponential stability of the second order equation (1.8) without damping term in the case of variable coefficients and/or delays.

It will be demonstrated that the control of the form (1.2), where m2, can stabilize system (1.1). We obtain results as regards the exponential stability of (1.8), which are based on the maximum principles for the second order delay differential equation (1.8). Denoting (Mx)(t)= x (t)+ i = 1 m p i (t)x(t τ i (t)), we can formulate the maximum inequalities principle [23] in the form: it follows from the inequalities (My)(t)(Mx)(t), y(0)x(0), y (0) x (0) that y(t)x(t) for all t[0,). Using the solutions’ representation (see (2.5) below), we can see that this principle is reduced to positivity of the Cauchy (fundamental) function of (1.8) and positivity of its nontrivial solution with nonnegative initial conditions. Nonoscillation of solutions in many cases leads to this property of the solution of second order equation (1.8). The positivity of the Cauchy function will open a way to the analysis of the asymptotic stability of nonlinear delay differential equations on the basis of the known schemes of upper and lower functions. The maximum boundedness principle [23] can be described in the form: there exists a constant Ω such that the solution of the initial value problem (Mx)(t)=f(t), x(0)=α, x (0)=β satisfies the inequality |x(t)|Ω(|f(t)|+|α|+|β|) for t[0,). A generalization of the Perron theorem (see, for example, Theorem 3.5 in [1]) claims the exponential stability in this case.

The paper is organized as follows. After formulation of stabilization problem in Section 1, we discuss formulas of solution representations in Section 2. Simple results on stability of delay differential equations in the case m=2, which will be proven on the basis of general assertions obtained in Section 5, are formulated in Section 3. Auxiliary results used in the proof of our main assertions can be found in Section 4. Results on the positivity of the Cauchy function C(t,s), nonoscillation and exponential stability are obtained in Sections 5 and 6. Discussion and open problems are presented in Section 7.

2 About representation of solutions for second order delay differential equations

The main object of this paper is the second order delay differential equation

x (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =f(t),t[0,+),
(2.1)
x(ξ)=φ(ξ)for ξ<0,
(2.2)

with measurable essentially bounded f, p i , φ, τ i (i=0,1,,n1), and τ i (t)0 for t0. In the traditional approach (see, for example, [4, 24, 25]), the homogeneous equation is considered as the equation

x (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =0,t[0,+),
(2.3)

with the initial function (2.2). If (2.3), (2.2) are considered for all possible continuous functions φ, the space of its solutions is infinite-dimensional. In [26] a homogeneous equation was defined as (2.3) with the zero initial function

x(ξ)=0for ξ<0.
(2.4)

This definition of the homogeneous equation allows researchers to build an analog of the classical general theory of ordinary differential equations (ODEs) for functional differential equations using the key notions of the classical ODEs theory, and it is based on the fact that the space of the solutions of the second order equation (2.3), (2.4) becomes two-dimensional. Note that (2.3), (2.2) can be written as (2.1), (2.4), where f(t)= i = 1 m p i (t)φ(t τ i (t)){1χ(t τ i (t))}, χ(t)=1 for t0, and χ(t)=0 for t<0. Thus (2.1), (2.4) is a nonhomogeneous one according to the approach of the book [1].

The general solution of (2.1), (2.4) (we understand the notion of the general solution exactly like in the theory of ODEs) can be represented in the form [26]

x(t)= 0 t C(t,s)f(s)ds+ x 1 (t)x(0)+ x 2 (t) x (0),
(2.5)

where x 1 (t), x 2 (t) are two solutions of the homogeneous equation (2.3), (2.4) satisfying the conditions

x 1 (0)=1, x 1 (0)=0, x 2 (0)=0, x 2 (0)=1,
(2.6)

and the kernel in this representation is called the Cauchy function C(t,s) of (2.1). Comparing this representation with the variation-of-constant formula of Hale [[3], Chapter 6, Section 6.1], we see that they almost coincide; the kernel in the integral term is called the fundamental function by his numerous followers. For every fixed s the function C(t,s) as a function of the variable t satisfies the equation [1, 26]

x (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =0,t[s,+),
(2.7)
x(ξ)=0for ξ<s,
(2.8)

and the initial conditions C(s,s)=0, C t (s,s)=1. Note that we will use this definition of the Cauchy function for the construction of C(t,s) below. The behavior of the fundamental system x 1 , x 2 of solutions of (2.3), (2.4) determines the existence and uniqueness of solutions of the boundary value problems for this equation. Positivity of the Cauchy function C(t,s) of (2.1), (2.4) is a basis of various theorems about differential inequalities (under corresponding conditions, the solution of an inequality is greater than the solution of an equation). These theorems help essentially in the study of stability for delay differential equations and their natural generalization - functional differential equations (FDEs) [1].

Let us formulate several definitions concerning stability.

Definition 2.1 Equation (2.3) is uniformly exponentially stable if there exist N>0 and α>0, such that the solution of (2.3), (2.9), where

x(ξ)=φ(ξ),ξ< t 0 ,x( t 0 )= x 0 , x ( t 0 )= x 0 ,
(2.9)

satisfies the estimate

|x(t)|N e α ( t t 0 ) , t 0 t<+,
(2.10)

where N and α do not depend on t 0 .

Definition 2.2 We say that the Cauchy function C(t,s) of (2.1) satisfies the exponential estimate if there exist positive N and α such that

|C(t,s)|N e α ( t s ) ,0st<+.
(2.11)

It is well known that for (2.1) with bounded delays these two definitions are equivalent [1].

3 Tests of positivity of the Cauchy functions and stability

Let us try to understand in which direction nonoscillation of (1.8) can be established. Taking into account the well-known oscillation results [27], we come to the conclusion: only such smallness of the nonnegative coefficient p(t) that its integral on the semiaxis is bounded may be suitable for nonoscillation of (1.7). In the case of (2.1) the analogous smallness of the sum of the nonnegative coefficients p i (t) should be required for nonoscillation. Regarding stability, we know that the smallness of the delay (1.13) is necessary and sufficient for boundedness of all solutions to (1.7) with constant coefficient p(t)p>0. It seems strange even to hope that assertions on nonoscillation and exponential stability can be achieved in the cases of non-tending to zero coefficients and delays. We will demonstrate that, although the ordinary differential equation

x (t)+ { i = 1 m p i ( t ) } x(t)=0,t[0,+),
(3.1)

can be oscillating and asymptotically unstable, the delay equation (2.3), (2.4) under corresponding conditions on the coefficients p i (t) and delays τ i (t) is nonoscillating and exponentially stable. The basic idea of our approach is to avoid the condition on the nonnegativity of the coefficients p i (t) for all i=1,,m, and to allow terms with positive and terms with negative coefficients p i (t) to compensate each other.

In order to formulate several simple corollaries of the main results proven below in Section 5, let us consider the equation

x (t)+a(t)x ( t τ ( t ) ) b(t)x ( t θ ( t ) ) =0,t[0,+),
(3.2)
x(ξ)=0for ξ<0,
(3.3)

where a(t), b(t), τ(t), and θ(t) are measurable essentially bounded nonnegative functions. We denote

q = essinf t 0 q(t), q = esssup t 0 q(t).
(3.4)

Theorem 3.1 Assume that 0τ(t)θ(t), 0b(t)a(t),

4 { a ( t ) b ( t ) } [ b ( θ τ ) ] 2 ,t[0,+),
(3.5)
0< [ b ( θ τ ) ] θ 1 e .
(3.6)

Then

  1. (1)

    the Cauchy function C(t,s) of (3.2) is nonnegative for 0s<t<+;

  2. (2)

    if there exists such positive ε such that

    a(t)b(t)ε,
    (3.7)

then the Cauchy function C(t,s) of (3.2) satisfies the exponential estimate (2.11) and the integral estimate

sup t 0 0 t C(t,s)ds 1 ε ;
(3.8)
  1. (3)

    if there exists lim t {a(t)b(t)}=k, with k>0, then

    lim t 0 t C(t,s)ds= 1 k .
    (3.9)

Corollary 3.1 Assume that the delays τ(t)τ, θ(t)θ are constants and

0<ε4 { a ( t ) b ( t ) } b 2 ( θ τ ) 2 ,t[0,+),
(3.10)
0< b (θτ)θ 1 e .
(3.11)

Then the assertions of Theorem  3.1 are true.

Consider the equation

x (t)+a(t)x(tτ)=f(t),t[0,+),a(t) a >0.
(3.12)

There exist unbounded solutions of corresponding homogeneous equation x (t)+a(t)x(tτ)=0, t[0,+) (see, for example, [[4], Chapter III, Section 16, p.106] in the case of constant a, and noted above in Section 1 sufficient conditions for existence of unbounded solutions from the paper [6]). Together with the fact of the oscillation [27] of all solutions this leads to chaos in behavior of solutions in the sense that small errors in a measurement of x( t 0 ) and x ( t 0 ) can imply unpredictable changes in x( t 0 +ω) and x ( t 0 +ω). Thus (3.12) is unstable. To stabilize its solution x(t) to the given ‘trajectory’ y(t) satisfying this equation, we choose the control in the form

u(t)=b(t) [ x ( t θ ) y ( t θ ) ] .
(3.13)

Subtracting the equality

y (t)+a(t)y(tτ)=f(t),t[0,+),

from the equation

x (t)+a(t)x(tτ)=b(t) [ x ( t θ ) y ( t θ ) ] +f(t),t[0,+),

we get the equation

z (t)+a(t)z(tτ)b(t)z(tθ)=0,t[0,+)

for the difference z(t)=x(t)y(t). The stabilization has been reduced to the stability of this equation. A possible algorithm to construct this stabilizing control is clear now: first we choose the delay θ close to τ such that condition (3.11) is fulfilled, then we choose b(t) close to a(t) such that condition (3.10) is fulfilled.

Example 3.1 Stabilizing (3.12), where a(t)a, let us choose the control in the form (3.13) with constant coefficient b(t)b. We come to the study of the exponential stability of the equation

x (t)+ax(tτ)bx(tθ)=g(t),t[0,+),
(3.14)

with constant coefficient and delays and g(t)=f(t)+by(tθ). We can choose θτ= 1 e b θ from (3.11); then from (3.10) we get the following condition of the exponential stability:

0<4{ab} 1 e 2 θ 2 .
(3.15)

Example 3.2 The equation

x (t)+a(t)x(tτ)=0,a(t)+,t[0,+),τ=const,
(3.16)

where a(t) a >0, possesses oscillating solutions with amplitudes tending to infinity [6] that leads to the chaos in behavior of its solutions. This equation can also be stabilized by the control in the form (3.13). Consider, for example, the equation

x (t)+tx(tτ)=0,t[1,+),τ=const.
(3.17)

If we choose b(t)=tΔ, θ(t)=τ+ γ t , then the stabilization can be achieved by the control (3.13) with the parameters satisfying the inequalities

0<2 Δ <γ< 1 τ e .
(3.18)

In the following assertion we assume the smallness of the difference of the delays θτ and do not assume the smallness of the delay θ.

Consider for simplicity the equation

x (t)+a(t)x(tτ)b(t)x(tθ)=0,t[0,+),
(3.19)
x(ξ)=0for ξ<0,
(3.20)

with constant delays τ and θ and positive coefficients a(t) and b(t).

Theorem 3.2 Assume that 0<τ<θ, there exists a positive ε such that

ε { a ( t ) b ( t ) } 1 4 [ b ( θ τ ) ] 2 ,t[0,+)
(3.21)

and

1 a exp { b ( θ τ ) 2 4 } arctg b ( θ τ ) 2 a exp { b ( θ τ ) 2 4 } >θτ
(3.22)

is fulfilled.

Then

  1. (1)

    the Cauchy function C(t,s) of (3.19) is nonnegative for 0s<t<+;

  2. (2)

    the solutions x 1 (t), x 2 (t) of (3.19), satisfying initial conditions (2.6), are nonegative for 0<t<+;

  3. (3)

    the Cauchy function C(t,s) of (3.19) satisfies the exponential estimate (2.11) and the integral estimate (3.8);

  4. (4)

    if there exists lim t {a(t)b(t)}=k, then equality (3.9) is fulfilled.

Example 3.3 Consider the equation

x (t)+x(t)bx(tθ)=0,t[0,+),
(3.23)
x(ξ)=0,ξ<0.
(3.24)

It is clear from the definition of the Cauchy function (see (2.7) and (2.8)) that for t<θ, this equation is equivalent to the ordinary differential equation

x (t)+x(t)=0,t[0,θ],

which Cauchy function is C(t,s)=sin(ts). Let us demonstrate that condition (3.6) is essential for positivity of the Cauchy function C(t,s) in Theorem 3.1, assuming that the numbers b and θ are chosen such that all other conditions of Theorem 3.1 are fulfilled. If π<θ, then in the triangle 0st<θ, its Cauchy function C(t,s)=sin(ts) changes its sign. The same example demonstrates that condition (3.22) is essential in Theorem 3.2 for the positivity of C(t,s). Note that in the case π 2 <θ, the solution x 1 (t)=cost satisfying conditions (2.6) changes its sign on the interval [0,θ]. It is clear that inequality (3.22) implies π 2 >θ.

Remark 3.1 Note the closely related results for the first order delay differential equation

x ( t ) + a ( t ) x ( t τ ( t ) ) b ( t ) x ( t θ ( t ) ) = 0 , t [ 0 , + ) , x ( ξ ) = 0 , ξ < 0 ,
(3.25)

which were obtained in [23]; see also Chapter 14 in the book [28]. Although the solutions of the equation

x (t)+a(t)x ( t τ ( t ) ) =0,t[0,+),
(3.26)

can oscillate and the amplitudes of their solutions tend to infinity, the solutions of (3.25) cannot oscillate and tend to zero exponentially. For example, consider the equation

x (t)+ax(tτ)=0,aτ> π 2 ,t[0,+),
(3.27)

with constant coefficient and delay, which has oscillating solutions with amplitudes tending to infinity, but the equation

x (t)+ax(tτ)bx(tθ)=0,t[0,+),
(3.28)

is nonoscillating and exponentially stable if b>0, 0<(ab)θ 1 e , 0<a(θτ) 1 e .

4 Auxiliary results

Let us formulate several results on systems of functional differential equations with Volterra operators. We mean the Volterra operators according to the following Tikhonov definition. Assume that B: C [ 0 , ) L [ 0 , ) is a linear bounded operator, where C [ 0 , ) and L [ 0 , ) are the spaces of continuous and essentially bounded functions x:[0,+)(,+) respectively.

Definition 4.1 We say that an operator B: C [ 0 , ) L [ 0 , ) is a Volterra one if for every two functions y 1 and y 2 and every ω, the equality y 1 (t)= y 2 (t) for t[0,ω] implies the equality (B y 1 )(t)=(B y 2 )(t) for t[0,ω].

Consider the system

( M i x)(t) x i (t)+ j = 1 2 ( B i j x j )(t)= f i (t),t[0,+),i=1,2,
(4.1)

where B i j : C [ 0 , ) L [ 0 , ) are linear bounded Volterra operators for i,j=1,2. The general solution of this system can be represented in the form [1]

x(t)= 0 t K(t,s)f(s)ds+K(t,0)x(0),
(4.2)

where K(t,s)= { K i j ( t , s ) } i , j = 1 , 2 is called the Cauchy matrix of system (4.1), x=col{ x 1 , x 2 }, f=col{ f 1 , f 2 }.

Lemma 4.1 Let the Cauchy function k(t,s) of the scalar equation

x 1 (t)+( B 11 x 1 )(t)= f 1 (t),t[0,+),
(4.3)

be positive in 0st<+, B 12 and ( B 21 ) be positive operators. Then the following assertions are equivalent:

  1. (1)

    K 21 (t,s)0 and K 22 (t,s)>0 for 0st<+;

  2. (2)

    there exists a vector function v=col{ v 1 , v 2 } with absolutely continuous components and essentially bounded derivatives v =col{ v 1 , v 2 } such that

    v 2 (t)>0,( M i v)(t)0,i=1,2,t[0,+), v 1 (0)0.
    (4.4)

This lemma follows from Theorem 16.6 for the case n=2 (see [[28], p.409]).

In order to formulate the conditions of positivity of the Cauchy function k(t,s) of the scalar first order equation (4.3), let us introduce the function H:[0,+)[0,+) which describes the ‘size’ of the memory of the operator B 11 .

Definition 4.2 Let us define the function H(t):[0,+)[0,+) as the maximal possible value on [0,t] for which the equality y 1 (s)= y 2 (s) for s[H(t),+) for every two continuous functions y 1 and y 2 :[0,+)(,+) implies the equality ( B 11 y 1 )(s)=( B 11 y 2 )(s) for almost every s[t,+).

Lemma 4.2 Let B 11 : C [ 0 , ) L [ 0 , ) be a positive linear Volterra operator and

H ( t ) t ( B 11 1)(s)ds 1 e ,t(0,+),
(4.5)

then nontrivial solutions of the homogeneous equation

x 1 (t)+( B 11 x 1 )(t)=0,t[0,+),
(4.6)

have no zeros for t[0,+) and its Cauchy function C(t,s)>0 for 0st<+.

This lemma follows from Theorem 15.7 (see [[28], p.358]).

Denote H = esssup t 0 {tH(t)}.

Denoting B 11 the norm of the operator B 11 : C [ 0 , ) L [ 0 , ) , as a corollary of this assertion we get the following.

Lemma 4.3 If B 11 H 1 e , then the Cauchy function C(t,s) of (4.3) is positive for 0st<+.

Remark 4.1 Let us consider the second order scalar functional differential equation

x (t)+ ( B 11 x ) (t)+( B 12 x)(t)= f 1 (t),t[0,ω],
(4.7)

and the corresponding system

x 1 (t)+ j = 1 2 ( B 1 j x j )(t)= f 1 (t), x 2 (t) x 1 (t)=0,t[0,ω].
(4.8)

The Cauchy function C(t,s) of (4.7) coincides with the entry K 21 (t,s) of the Cauchy matrix K(t,s) of system (4.8) and C t (t,s) coincides with the entry K 11 (t,s). If the function y(t) is the solution of the initial value problem

x (t)+ ( B 11 x ) (t)+( B 12 x)(t)=0,t[0,ω],
(4.9)
x(0)=1, x (0)=0,
(4.10)

then K 22 (t,0)=y(t) and K 12 (t,0)= y (t). Note that if z(t) is the solution of the initial value problem consisting of equation (4.9) and the conditions

x(0)=0, x (0)=1,
(4.11)

then K 21 (t,0)=z(t), K 11 (t,0)= z (t).

5 Main results

Let us consider the equation

(Mx)(t) x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) =f(t),t[0,+),
(5.1)

and the corresponding homogeneous equation

x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) =0,t[0,+),
(5.2)

where

x(ξ)=0for ξ<0.

Denote

χ(t,s)= { 1 , t s , 0 , t < s .

Theorem 5.1 Let ( 1 ) i + 1 p i (t)0 for i=1,,2m, p 2 i 1 (t)+ p 2 i (t)0 for i=1,,m, t[0,+), and the Cauchy function of the first order equation

y (t)+ i = 1 m | p 2 i (t)|χ ( t τ 2 i ( t ) , 0 ) t τ 2 i ( t ) t τ 2 i 1 ( t ) y(s)ds=0,t[0,+),
(5.3)

be positive for 0st<+. Then the following assertions are equivalent:

  1. (1)

    there exists a function v with absolutely continuous and essentially bounded derivative v and essentially bounded derivative v such that

    v(t)>0, v (t)0,(Mv)(t)0,t[0,+);
    (5.4)
  2. (2)

    there exists a bounded absolutely continuous function u with essentially bounded derivative u such that

    u ( t ) 0 , u 2 ( t ) u ( t ) + i = 1 2 m p i ( t ) χ ( t τ i ( t ) , 0 ) exp { t τ i ( t ) t u ( s ) d s } 0 , t [ 0 , + ) ;
    (5.5)
  3. (3)

    the Cauchy function C(t,s) of (5.1) is nonnegative for 0s<t<+, and solutions x 1 (t), x 2 (t) of (5.2), satisfying initial conditions (2.6), are such that x 1 (t)>0, x 2 (t)0 for 0<t<+.

The proof of Theorem 5.1 is based on several auxiliary results.

Lemma 5.1 Assertions (1) and (2) of Theorem  5.1 are equivalent.

Proof (2) (1). In order to prove this implication, let us choose the function v(t) in the condition (1) in the form

v(t)=exp { 0 t u ( s ) d s }
(5.6)

and demonstrate that it satisfies the condition (1). Differentiating, we get

v (t)=u(t)exp { 0 t u ( s ) d s }
(5.7)

and

v (t)= [ u 2 ( t ) u ( t ) ] exp { 0 t u ( s ) d s } .
(5.8)

We see that v(t)>0, v (t)0. In the calculation of (Mv)(t), one gets, carrying the exponent (5.6) out of the brackets,

( M v ) ( t ) = exp { 0 t u ( s ) d s } × { u 2 ( t ) u ( t ) + i = 1 2 m p i ( t ) χ ( t τ i ( t ) , 0 ) exp [ t τ i ( t ) t u ( s ) d s ] } 0 .
(5.9)

It is clear now that inequalities (5.5) imply that (Mv)(t)0. The implication (2) (1) has been proven.

(1) (2). Denote u(t)= v ( t ) v ( t ) , where the function v(t) from the assertion (1), and demonstrate that this function satisfies the assertion (2). Obviously the function u(t) is nonnegative. The function v(t) satisfies the equation v (t)+u(t)v(t)=0 and consequently it can be represented in the form (5.6). The derivatives of the function v(t) are defined by (5.7) and (5.8). From the inequality (Mv)(t)0 we get inequality (5.9) and consequently condition (5.5) is fulfilled.

The implication (1) (2) has been proven. □

Proof of Theorem 5.1 (1) (3). Let us rewrite (5.1) in the form

( M x ) ( t ) x ( t ) + i = 1 m | p 2 i ( t ) | χ ( t τ 2 i ( t ) , 0 ) { x ( t τ 2 i 1 ( t ) ) x ( t τ 2 i ( t ) ) } + i = 1 m | p 2 i ( t ) | x ( t τ 2 i 1 ( t ) ) χ ( t τ 2 i 1 ( t ) , 0 ) ( 1 χ ( t τ 2 i ( t ) , 0 ) ) + i = 1 m [ p 2 i 1 ( t ) + p 2 i ( t ) ] χ ( t τ 2 i 1 ( t ) , 0 ) x ( t τ 2 i 1 ( t ) ) = f ( t ) , t [ 0 , + ) .
(5.10)

Consider the system

x 1 ( t ) + i = 1 m | p 2 i ( t ) | χ ( t τ 2 i ( t ) , 0 ) t τ 2 i ( t ) t τ 2 i 1 ( t ) x 1 ( s ) d s + i = 1 m | p 2 i ( t ) | x 2 ( t τ 2 i 1 ( t ) ) χ ( t τ 2 i 1 ( t ) , 0 ) ( 1 χ ( t τ 2 i ( t ) , 0 ) ) + i = 1 m [ p 2 i 1 ( t ) + p 2 i ( t ) ] x 2 ( t τ 2 i 1 ( t ) ) = f ( t ) , x 2 ( t ) x 1 ( t ) = 0 , t [ 0 , + ) ,
(5.11)

where

x 1 (ξ)= x 2 (ξ)=0for ξ<0.
(5.12)

To compare this system with system (4.1), we observe that the operators B i j : C [ 0 , ) L [ 0 , ) are of the following forms:

( B 11 x 1 )(t) i = 1 m | p 2 i (t)|χ ( t τ 2 i ( t ) , 0 ) t τ 2 i ( t ) t τ 2 i 1 ( t ) x 1 (s)ds,
(5.13)
( B 12 x 2 ) ( t ) = i = 1 m [ p 2 i 1 ( t ) + p 2 i ( t ) ] χ ( t τ 2 i 1 ( t ) , 0 ) x 2 ( t τ 2 i 1 ( t ) ) ( B 12 x 2 ) ( t ) = + i = 1 m | p 2 i ( t ) | x 2 ( t τ 2 i 1 ( t ) ) χ ( t τ 2 i 1 ( t ) , 0 ) ( 1 χ ( t τ 2 i ( t ) , 0 ) ) ,
(5.14)
( B 21 x 1 )(t)= x 1 (t),( B 22 x 2 )(t)=0.
(5.15)

Under the conditions ( 1 ) i + 1 p i (t)0, p 2 i 1 (t)+ p 2 i (t)0, τ 2 i 1 (t) τ 2 i (t) of Theorem 5.1, the operators B 11 , B 12 , and ( B 21 ) are positive. Note that, according to the condition of Theorem 5.1, the Cauchy function of the first order equation (5.3) is positive for 0st+. The vector function col{ v 1 (t), v 2 (t)}, where v 1 (t)= v (t), v 2 (t)=v(t), satisfies all conditions of assertion (2) of Lemma 4.1. According to Lemma 4.1, we get the positivity of the element K 22 (t,s) and the nonnegativity of the element K 21 (t,s) of the Cauchy matrix of system (5.11) for 0st+. Remark 4.1 completes now the proof of the implication (1) (3).

To prove the implication (3) (1) we choose v(t)= y 1 (t), where y 1 (t) is the solution of the initial problem

(My)(t)=0,y(0)=1, y (0)=0.
(5.16)

Noting that the equivalence (1) (2) was obtained in Lemma 5.1, we can conclude that Theorem 5.1 has been proven. □

Denote H = esssup t 0 max 1 i 2 m τ i (t). For bounded τ i (t) for i=1,,m, we have H <+.

Remark 5.1 The inequality

t H t i = 1 m | p 2 i (s)| [ τ 2 i ( s ) τ 2 i 1 ( s ) ] ds 1 e for t0,
(5.17)

where we set p 2 i (s)=0 for s<0, implies the positivity of the Cauchy function of the first order equation (5.3) (see Theorem 15.7, p.358 in [28]).

Remark 5.2 There is a corresponding inconvenience in choosing the functions v(t) because of the equality v(t τ i (t))=0 for t τ i (t)<0. Trying to avoid this inconvenience, we can make the following trick.

Let us define the auxiliary operator B: C [ H , + ) L [ H , + ) , where C [ H , + ) and L [ H , + ) are the spaces of continuous and of essentially bounded functions x:[ H ,+)(,+), respectively, by the formula

(By)(t)= { β y ( t ) , H t 0 , i = 1 m p 2 i ( t ) t τ 2 i ( t ) t τ 2 i 1 ( t ) y ( s ) d s , t 0 ,
(5.18)

where the parameter β will be defined below in the formulation of Theorem 5.2. Define the norm of this operator as B=max{β, esssup t 0 i = 1 m | p 2 i (t)|[ τ 2 i (t) τ 2 i 1 (t)]}.

Let us define p 2 i 1 (t)0, p 2 i (t)0 for i=1,,m and τ j (t)0 for j=1,,2m, on the interval t[ H ,0), and consider the equation

(Mx)(t) x (t)+q(t) x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) =f(t),t[ H ,+),
(5.19)

where f:[ H ,+)(,+) is a measurable essentially bounded function and

q(t)= { β , H t 0 , 0 , t 0 ,

with the same β as in (5.18), which will be defined below.

It is clear from the definition of the Cauchy function by (2.7), (2.8) that the Cauchy functions of (5.1) and (5.19) coincide in 0st<+.

Theorem 5.2 Assume that ( 1 ) i + 1 p i (t)0 for i=1,,2m, p 2 i 1 (t)+ p 2 i (t)0, τ 2 i 1 (t) τ 2 i (t) for i=1,,m and t[0,+) and there exists a real positive number α such that

  1. (a)

    B H 1 e , where β=α in (5.18) defining the operator B;

  2. (b)

    the inequality

    α 2 + i = 1 2 m p i (t)exp { α τ i ( t ) } 0,t[0,+)
    (5.20)

is fulfilled.

Then

  1. (1)

    the Cauchy function C(t,s) of (5.19) is nonnegative for 0s<t<+;

  2. (2)

    if in addition there exists a positive ε such that

    i = 1 2 m p i (t)ε,
    (5.21)

then the Cauchy function C(t,s) of (5.1), its derivative C t (t,s)= K 11 (t,s) and solutions of the homogeneous equation Mx=0 satisfy the exponential estimate and the estimate

sup t 0 0 t C(t,s)ds 1 ε
(5.22)

is true.

Proof of Theorem 5.2 Let us start with the proof of the assertion (1). We set v(t)=exp{α(t H )} for t[ H ,+) in the assertion (1) of the following reformulation of the assertion (1) (3) of Theorem 5.1 for (5.19).

Lemma 5.2 Let ( 1 ) i + 1 p i (t)0 for 1=1,,2m, p 2 i 1 (t)+ p 2 i (t)0, τ 2 i 1 (t) τ 2 i (t) for i=1,,m, t[0,+) and the Cauchy function of the first order equation

y (t)+(By)(t)=f(t),t[ H ,+),
(5.23)

where β=α, be positive for H st<+. Then the assertion (2) follows from the assertion (1), where

  1. (1)

    there exists a function v with absolutely continuous and essentially bounded derivative v and essentially bounded derivative v such that

    v(t)>0, v (t)0,(Mv)(t)0,t[ H ,+);
  2. (2)

    the Cauchy function C(t,s) of (5.19) is nonnegative for 0s<t<+, and the solutions x 1 (t), x 2 (t) of the homogeneous equation

    (Mx)(t) x (t)+q(t) x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) =0,t[ H ,+),

satisfying the initial conditions (2.6), are such that x 1 (t)>0, x 2 (t)0 for H <t<+.

Substituting this function v(t)=exp{α(t H )} into the differential operation M defined by (5.19), we obtain the following two inequalities on [ H ,0] and [0,+). The first one is

exp ( α ( t H ) ) { α 2 β α } 0,t[ H ,0),
(5.24)

and it is satisfied for β=α. This explains how to choose β in the definition of the operator B by (5.18) and the coefficient q(t) in (5.19). The second inequality is

α 2 exp(αt)+ i = 1 2 m p i (t)exp { α ( t τ i ( t ) ) } 0,t[0,+),
(5.25)

and, after carrying the exponent out of the brackets, we get inequality (5.20). Thus inequalities (5.24) and (5.25) are results of a choice of the corresponding function v(t) in assertion (1) of Lemma 5.2.

Equation (5.19) can be rewritten in the form

(Mx)(t) x (t)+ ( B x ) (t)+ i = 1 2 m { p 2 i 1 ( t ) + p 2 i ( t ) } x ( t τ i ( t ) ) =f(t),t[ H ,+),

where the operator B is defined by (5.18). The condition (a) implies, according to Lemma 4.3, the positivity of the Cauchy function of the first order equation (5.23).

Now we see that all conditions of Lemma 5.2 and assertion (1) are fulfilled. According to Lemma 5.2, we obtain assertion (1) of Theorem 5.2.

Let us prove assertion (2) of Theorem 5.2. The general solution of (5.19) can be represented in the form

x(t)= H t C(t,s)f(s)ds+ x 1 (t)x ( H ) + x 2 (t) x ( H ) ,
(5.26)

where solutions x 1 and x 2 of the homogeneous equation (Mx)(t)=0, t[ H ,+) satisfy the conditions

x 1 ( H ) =1, x 1 ( H ) =0, x 2 ( H ) =0, x 2 ( H ) =1.
(5.27)

The function y(t) 1 ε satisfies the initial conditions

x ( H ) = 1 ε , x ( H ) =0,
(5.28)

and the equality

( M 1 ε ) (t)=0,t [ H , 0 ] .
(5.29)

Using condition (5.21), we obtain

( M 1 ε ) (t)1,t[0,+).
(5.30)

From representation (5.26), using (5.27)-(5.30), we get

1 ε = 0 t C(t,s) ( M 1 ε ) (s)ds+ x 1 (t) 1 ε ,t[0,+),
(5.31)

and, taking into account the nonnegativity of C(t,s), we get

0 0 t C(t,s)ds+ x 1 (t) 1 ε 1 ε ,t[0,+).
(5.32)

Positivity of x 1 (t) and nonnegativity of C(t,s) imply that

0 t C(t,s)ds 1 ε ,t[0,+);
(5.33)

this proves the integral estimate (5.22). It is clear that

0 t C(t,s)|f(s)|ds 1 ε ,t[0,+)
(5.34)

for every measurable essentially bounded f such that |f(t)|1 for t[0,+). It means that the function

y(t) 0 t C(t,s)|f(s)|ds,t[0,+),
(5.35)

is bounded for every function f measurable and bounded on the semiaxis t[0,+). Each of the solutions of (5.19) can be presented on [0,+) as a solution of the equation

(Mx)(t)=ψ(t),t[0,+),
(5.36)

where ψ(t) is a corresponding measurable bounded function such that ψ(t)=f(t) for t H . Then the solution x of (5.36) is bounded and consequently the solution of (5.19) is bounded for every measurable bounded function f(t), t[ H ,+) and

H t C(t,s)|f(s)|ds,t[0,+),
(5.37)

is bounded for every measurable bounded function f.

Let us prove that the function x 2 (t)=C(t, H ) is bounded on the semiaxis t[ H ,+). The function y(t)=1exp{t+ H } satisfies the Cauchy problem

(Mx)(t)=f(t),t[ H ,+),
(5.38)
x ( H ) =0, x ( H ) =1,
(5.39)

where the differential operation M defined by (5.19) and f(t)=(My)(t). It is clear that this function f is bounded on the semiaxis t[ H ,+). The solution y of (5.38), (5.39) can be written in the form

1exp { t + H } = H t C(t,s)f(s)ds+ x 2 (t),t[ H ,+).
(5.40)

It follows now from the boundedness of the integral (5.37) for every bounded f that x 2 (t) is bounded for t[ H ,+).

Consider the system

{ ( M 1 y ) ( t ) y 1 ( t ) + ( B y 1 ) ( t ) + i = 1 m { p 2 i 1 ( t ) + p 2 i ( t ) } y 2 ( t τ 2 i 1 ( t ) ) = f 1 ( t ) , ( M 2 y ) ( t ) y 2 ( t ) y 1 ( t ) = f 2 ( t ) , t [ H , + ) .
(5.41)

If α satisfies inequality (5.20), the functions v 1 (t)=αexp{α(t+ H )}, v 2 (t)=exp{α(t+ H )} satisfy the assertion (2) of Lemma 4.1. The condition (a) implies that all other conditions of Lemma 4.1 are also fulfilled. Now, according to Lemma 4.1, the entries in the second row of the Cauchy matrix K(t,s)= { K i j ( t , s ) } i , j = 1 2 of system (5.41) satisfy the inequalities K 21 (t,s)0, K 22 (t,s)>0 for 0st<+. Let us substitute the constant vector z=col{ z 1 , z 2 }, where z 1 =1 and z 2 = 1 + B ε instead of y 1 (t) and y 2 (t), respectively, into system (5.41). Condition (5.21) implies that ( M 1 y)(t)1, ( M 2 y)(t)=1. Consider now system (5.41), where f 1 (t)( M 1 z)(t), f 2 (t)1 for t[ H ,+). From representation (4.2) of the solution we get for the second component of the solution vector

z 2 = H t K 21 (t,s) f 1 (s)ds+ H t K 22 (t,s)ds+ K 22 (t,H) 1 + B ε K 21 ( t , H ) .
(5.42)

In Remark 4.1 we explained that K 21 (t,s)=C(t,s), K 21 (t, H )= x 2 (t), K 22 (t, H )= x 1 (t), where C(t,s) is the Cauchy function and x 1 , x 2 are solutions of the equation

(Mx)(t) x (t)+q(t) x (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =0,t[ H ,+),
(5.43)

satisfying conditions (5.27). Now it follows from the boundedness of the integral H t K 21 (t,s)ds and solutions x 1 , x 2 that the integral H t K 22 (t,s)ds is bounded, and, using the positivity of K 22 (t,s), we can claim that the integral H t K 22 (t,s)f(s)ds is bounded for every measurable bounded f(t), t[ H ,+). We see that the component y 2 of the solution vector to system (5.41) is bounded for every measurable bounded f 1 and f 2 . The positivity of the Cauchy function of the first order equation (5.25) together with condition (5.21) imply the boundedness of the component y 1 of solution vector to system (5.41). Now a generalization of the Perron theorem (see, for example, Theorem 3.5 in [1]) claims that all entries of the Cauchy matrix K(t,s)= { K i j ( t , s ) } i , j = 1 2 satisfy the exponential estimate (2.11). Each of the solutions of (5.1) coincides with a corresponding solution of (5.19) with the measurable essentially bounded right-hand side ψ(t) such that ψ(t)=f(t) for t H . Thus for every bounded right-hand side, the solution of (5.1) is bounded, and the exponential estimate (2.11) of the Cauchy function of (5.1), its derivative C t (t,s)= K 11 (t,s), and the solutions of the homogeneous equation Mx=0 satisfy the exponential estimate (2.11).  □

Corollary 5.1 If under the conditions of Theorem  5.2 there exists a positive limit

lim t i = 1 2 m p i (t)=k,
(5.44)

then

lim t 0 t C(t,s)ds= 1 k .
(5.45)

The proof of Corollary 5.1 can be obtained from equality (5.31) and the fact of the exponential estimate of all elements of the Cauchy matrix of system (5.41).

Remark 5.3 The positivity of ε in condition (5.21) is essential, as the following example demonstrates. Consider the equation

x (t)+x ( t | sin 1 2 t | ) x ( t 2 | sin 1 2 t | ) =0,t[0,+).
(5.46)

Here ε=0 and one of solutions is a constant and does not tend to zero when t+.

Consider the equation

(Mx)(t) x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) + j = 1 n q j (t)x ( t θ j ( t ) ) =f(t),t[0,+),
(5.47)

and the corresponding homogeneous equation

x (t)+ i = 1 2 m p i (t)x ( t τ i ( t ) ) + j = 1 n q j (t)x ( t θ j ( t ) ) =0,t[0,+),
(5.48)

where

x(ξ)=0for ξ<0;
(5.49)

q j , θ j , p i and τ i are measurable essentially bounded functions for j=1,,n, i=1,,m.

Theorem 5.3 Let all assumptions of Theorem  5.2 be fulfilled.

Then

  1. (1)

    if q j (t)0 for t[0,+), then the Cauchy function C(t,s) of (5.47) is positive for 0s<t<+;

  2. (2)

    if there exist positive ε 0 and ε such that

    i = 1 2 m p i (t)ε,ε ε 0 j = 1 n | q j (t)|,t[0,+),
    (5.50)

then the Cauchy function C(t,s) of (5.47) satisfies the exponential estimate (2.11) and the integral estimate

sup t 0 0 t |C(t,s)|ds 1 ε 0 .
(5.51)

Proof Using the formula of the solutions’ representation (2.5), we can conclude that the solution of the initial problem (5.47), (5.49), (5.52), where

x(0)=0, x (0)=0,
(5.52)

satisfies the integral equation

x(t)=(Kx)(t)+ψ(t),
(5.53)

and the operator K: C [ 0 , ) C [ 0 , ) is defined by equalities (5.54), (5.49), where

(Kx)(t)= 0 t W(t,s) j = 1 n q j (s)x ( s θ i ( s ) ) ds,
(5.54)

where W(t,s) is the Cauchy function of (5.1), and

ψ(t)= 0 t W(t,s)f(s)ds.
(5.55)

According to Theorem 5.2, we obtain the positivity of W(t,s). The operator K in the case of nonpositive q j (t) is positive. If we consider this operator K on every finite interval [0,ω] its spectral radius ρ(K) is zero, and one can write

x(t)= ( I K ) 1 ψ(t)=ψ(t)+(Kψ)(t)+ ( K 2 ψ ) (t)+.
(5.56)

If f0, then ψ0 and consequently x0. Reference to (2.5) completes the proof of the assertion (1).

In order to prove the assertion (2), we have to demonstrate that the solution x is bounded for every bounded right-hand side f and then the extension of the Bohl-Perron implies the exponential estimate (2.11) of the Cauchy matrix (see Theorem 3.5 in the book [1]). According to Theorem 5.2, the Cauchy function W(t,s) of (5.1) satisfies the exponential estimate (2.11). This implies the boundedness of ψ for every bounded f. If the norm of the operator K: C [ 0 , ) C [ 0 , ) is less than one, then the bounded operator ( I K ) 1 exists. We have only to prove that K<1.

All conditions of Theorem 5.2 are fulfilled. According to Theorem 5.2, we have

W sup t 0 0 t W(t,s)ds 1 ε
(5.57)

and

KW esssup t 0 j = 1 n | q j (t)| 1 ε esssup t 0 j = 1 n | q j (t)| ε ε 0 ε <1.
(5.58)

This, according to our observations above, completes the proof of exponential estimate (2.11) of the Cauchy function of (5.47).

Now let us prove estimate (5.51). For the solution x of the initial problem (5.47), (5.49), (5.52), we have two representations

x(t)= 0 t C(t,s)f(s)ds,
(5.59)

and (5.56). Estimating for |f(t)|1, we obtain

|x(t)|=| 0 t C(t,s)f(s)ds| ( I K ) 1 W 1 1 ε ε 0 ε 1 ε = 1 ε 0 .
(5.60)

This completes the proof of the assertion (2) of Theorem 5.3. □

6 Proofs of Theorems 3.1 and 3.2

Proof of Theorem 3.1 For (3.2), inequality (5.20) is of the following form:

α 2 +a(t)exp { α τ ( t ) } b(t)exp { α θ ( t ) } 0,t[0,+),
(6.1)

which can be written as

α 2 + b ( t ) { exp [ α τ ( t ) ] exp [ α θ ( t ) ] } + ( a ( t ) b ( t ) ) exp { α τ ( t ) } 0 , t [ 0 , + ) ,
(6.2)

or, using the Lagrange theorem,

α 2 α b ( t ) ( θ ( t ) τ ( t ) ) exp [ α η ( t ) ] + ( a ( t ) b ( t ) ) exp { α τ ( t ) } 0 , t [ 0 , + ) ,
(6.3)

where τ(t)η(t)θ(t). Setting α= 1 2 { b ( t ) ( θ ( t ) τ ( t ) ) } we get

1 4 { b ( t ) ( θ ( t ) τ ( t ) ) } 2 1 2 { b ( t ) ( θ ( t ) τ ( t ) ) } b ( t ) ( θ ( t ) τ ( t ) ) × exp [ 1 2 { b ( t ) ( θ ( t ) τ ( t ) ) } η ( t ) ] + ( a ( t ) b ( t ) ) exp { 1 2 [ b ( t ) ( θ ( t ) τ ( t ) ) ] τ ( t ) } 0 , t [ 0 , + ) ,
(6.4)
( a ( t ) b ( t ) ) exp { 1 2 [ b ( t ) ( θ ( t ) τ ( t ) ) ] τ ( t ) } 1 4 { b ( t ) ( θ ( t ) τ ( t ) ) } b ( t ) ( θ ( t ) τ ( t ) ) × exp { 1 2 [ b ( t ) ( θ ( t ) τ ( t ) ) ] η ( t ) } , t [ 0 , + ) ,
(6.5)

which is fulfilled if

4 ( a ( t ) b ( t ) ) { b ( t ) ( θ ( t ) τ ( t ) ) } 2 ,t[0,+).
(6.6)

The inequality in the condition (a) of Theorem 5.2 is fulfilled if inequality (3.6) is satisfied.

Now Theorem 3.1 follows from Theorem 5.2 and Corollary 5.1. □

Remark 6.1 We have not yet used the assertion (2) of Theorem 5.1 to obtain sufficient conditions for the positivity of the Cauchy function C(t,s) of (3.2). Let us build the function u(t), proving the following assertion.

Proof of Theorem 3.2 To construct the function u(t) on the interval [0,τ] we solve the equation u (t)= u 2 (t). Its solution is

u(t)= 1 c 1 t ,
(6.7)

where c 1 is a constant. We will explain below how to choose c 1 . On the interval [τ,θ] we solve the inequality

u (t)= u 2 (t)+ a exp { b ( θ τ ) 2 2 } .
(6.8)

Separating variables and integrating, one gets

arctg u a exp { b ( θ τ ) 2 4 } = a exp { b ( θ τ ) 2 4 } (t+ c 2 ),
(6.9)

where c 2 is a constant, and

u(t)= a exp { b ( θ τ ) 2 4 } tg [ a exp { b ( θ τ ) 2 4 } ( t + c 2 ) ] .
(6.10)

Choosing c 2 such that

u(θ)= b ( θ τ ) 2 ,
(6.11)

we obtain

c 2 = 1 a exp { b ( θ τ ) 2 4 } arctg b ( θ τ ) 2 a exp { b ( θ τ ) 2 4 } θ.
(6.12)

For positivity of the function u(t) for t[τ,θ], we have to assume that c 2 >τ, which follows from inequality (3.22).

Let us now choose the constant c 1 in (6.7) such that

u(τ)= 1 c 1 τ ,
(6.13)

where u(τ) is defined by (6.10), (6.12).

Setting

u(t) 1 2 b (θτ)for tθ,
(6.14)

we have built the positive absolutely continuous function u(t) satisfying the condition (2) of Theorem 5.1. Now reference to Theorem 5.1 completes the proof of the assertions (1) and (2) of Theorem 3.2. To prove assertions (3) and (4) we actually repeat the proof of the assertion (2) of Theorem 5.2 (in the case H =0).

The proof of Theorem 3.2 have been completed. □

7 Open problems

Let us formulate a list of problems to be solved in the future.

  1. 1.

    To obtain results on the exponential stability of the equation

    x (t)+ j = 1 n q j (t) x ( t θ j ( t ) ) + i = 1 m p i (t)x ( t τ i ( t ) ) =f(t),t[0,+),
    (7.1)

where q j (t)0 for t[0,+). Results of this type were considered impossible. It was assumed in the previous works [1721] that q j (t)>0, and in the presented paper q j (t)0 for t[0,+), j=1,,m.

  1. 2.

    To obtain results as regards stabilization of the equation x ( n ) (t)=f(t), where n>2, to the trajectory y(t) by the control of the form

    u(t)= i = 1 m p i (t) { x ( t τ i ( t ) ) y ( t τ i ( t ) ) } ,t[0,+),
    (7.2)

without derivatives, i.e. to obtain results about the exponential stability of the equation

x ( n ) (t)+ i = 1 m p i (t)x ( t τ i ( t ) ) =f(t),t[0,+).
(7.3)

Results of this type were considered impossible.

  1. 3.

    To obtain results on oscillation/nonoscillation, existence of solutions tending to zero or tending to infinity for the second order equation (2.3) without the assumption of the nonnegativity of the coefficients.

  2. 4.

    To obtain results on stabilization of system (1.1) to the trajectory Y(t) by the control of the form

    u(t)= i = 1 m P i (t) { X ( t τ i ( t ) ) Y ( t τ i ( t ) ) } ,t[0,+),
    (7.4)

in the case of more general than diagonal matrices P i (t), i.e. to obtain results on the exponential stability of the system

X (t)+ i = 1 m P i (t)X ( t τ i ( t ) ) =g(t),t[0,+).
(7.5)
  1. 5.

    To obtain results as regards the exponential stability of systems of functional differential equations of different orders.

References

  1. Azbelev NV, Maksimov VP, Rakhmatullina LF: Introduction to Theory of Functional-Differential Equations. Nauka, Moscow; 1991.

    MATH  Google Scholar 

  2. Brands JJAM: Oscillation theorems for second-order functional-differential equations. J. Math. Anal. Appl. 1978, 63: 54-64. 10.1016/0022-247X(78)90104-X

    Article  MathSciNet  MATH  Google Scholar 

  3. Hale JK: Theory of Functional Differential Equations. Springer, New York; 1977.

    Book  MATH  Google Scholar 

  4. Myshkis AD: Linear Differential Equations with Delayed Argument. Nauka, Moscow; 1972.

    MATH  Google Scholar 

  5. Domoshnitsky A: Properties of Wronskian. In Kraevye Zadachi. Politekh. Inst. Perm, Perm; 1983: (in Russian)57-60. (in Russian)

    Google Scholar 

  6. Domoshnitsky A: Unboundedness of solutions and instability of differential equations of the second order with delayed argument. Differ. Integral Equ. 2001,14(5):559-576.

    MathSciNet  MATH  Google Scholar 

  7. Paatashvili DV: About unbounded solutions of linear differential equations of second order with delay argument. Differ. Uravn. 1989, 25: 774-780.

    MathSciNet  Google Scholar 

  8. Dolgii YF, Nikolaev SG: Instability of a periodic delay system. Differ. Uravn. 1998, 34: 463-468.

    MathSciNet  Google Scholar 

  9. Izjumova DV: About boundedness and stability of solutions of nonlinear functional-differential equations of the second order. Proc. Georgian Acad. Sci. 1980 (in Russian),100(2):285-288. (in Russian)

    MathSciNet  Google Scholar 

  10. Pinto M: Nonlinear delay-differential equations with small lag. Int. J. Math. Math. Sci. 1997, 20: 137-146. 10.1155/S0161171297000203

    Article  MathSciNet  MATH  Google Scholar 

  11. Pinto M: Asymptotic solutions for second order delay differential equations. Nonlinear Anal. 1999, 28: 1729-1740.

    Article  MathSciNet  MATH  Google Scholar 

  12. Romanenko EY, Sharkovskii AN: Asymptotic Behavior of Solutions of Linear Functional-Differential Equations. Akad. Nauk of Ukraina, Inst. Math., Kiev; 1978.

    MATH  Google Scholar 

  13. Derfel GA, Molchanov SA: On T. Kato’s problem regarding bounded solutions of differential-functional equations. Funkc. Anal. Prilozh. 1990, 24: 67-69. 10.1007/BF01077925

    Article  MathSciNet  MATH  Google Scholar 

  14. Dosla Z, Kiguradze IT: On boundedness and stability of solutions of second order linear differential equations with advanced arguments. Adv. Math. Sci. Appl. 1999,9(1):1-24.

    MathSciNet  MATH  Google Scholar 

  15. Cahlon B, Schmidt D: Stability criteria for second-order delay differential equations with mixed coefficients. J. Comput. Appl. Math. 2004, 170: 79-102. 10.1016/j.cam.2003.12.043

    Article  MathSciNet  MATH  Google Scholar 

  16. Cahlon B, Schmidt D: Stability criteria for certain second-order delay differential equations. Dyn. Contin. Discrete Impuls. Syst. 2003, 10: 593-621.

    MathSciNet  MATH  Google Scholar 

  17. Burton TA Mathematics in Science and Engineering 178. In Stability and Periodic Solutions of Ordinary Differential and Functional Differential Equations. Academic Press, San Diego; 1985.

    Google Scholar 

  18. Burton TA: Stability by Fixed Point Theory for Functional Differential Equations. Dover, New York; 2006.

    MATH  Google Scholar 

  19. Domoshnitsky A: About applicability of Chaplygin’s theorem to one component of the solution vector. Differ. Uravn. 1990, 26: 1699-1705.

    MathSciNet  Google Scholar 

  20. Bainov D, Domoshnitsky A: Stability of a second-order differential equations with retarded argument. Dyn. Stab. Syst. 1994,9(2):145-151. 10.1080/02681119408806173

    Article  MathSciNet  MATH  Google Scholar 

  21. Berezansky L, Braverman E, Domoshnitsky A: Stability of the second order delay differential equations with a damping term. Differ. Equ. Dyn. Syst. 2008,16(3):185-205. 10.1007/s12591-008-0012-4

    Article  MathSciNet  MATH  Google Scholar 

  22. Pontryagin LS: On zeros of some transcendental functions. Izv. Akad. Nauk SSSR, Ser. Mat. 1942 (the English translation in Am. Math. Soc. Transl., Ser. 2 1, 95-110 (1955)), 6: 115-134. (the English translation in Am. Math. Soc. Transl., Ser. 2 1, 95-110 (1955))

    MathSciNet  MATH  Google Scholar 

  23. Domoshnitsky A: Maximum principles and nonoscillation intervals for first order Volterra functional differential equations. Dyn. Contin. Discrete Impuls. Syst., Ser. A Math. Anal. 2008, 15: 769-814.

    MathSciNet  MATH  Google Scholar 

  24. Ladde GS, Lakshmikantham V, Zhang BG: Oscillation Theory of Differential Equations with Deviating Argument. Dekker, New York; 1987.

    MATH  Google Scholar 

  25. Norkin SB Transl. Math. Monographs 31. In Differential Equations of the Second Order with Retarded Argument. Am. Math. Soc., Providence; 1972.

    Google Scholar 

  26. Azbelev NV: About zeros of solutions of linear differential equations of the second order with delayed argument. Differ. Uravn. 1971,7(7):1147-1157.

    MathSciNet  MATH  Google Scholar 

  27. Koplatadze RG: On oscillatory properties of solutions of functional differential equations. Mem. Differ. Equ. Math. Phys. 1994, 3: 1-179.

    MathSciNet  MATH  Google Scholar 

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

    Book  MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander Domoshnitsky.

Additional information

Competing interests

The author declares that he has no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Domoshnitsky, A. Nonoscillation, maximum principles, and exponential stability of second order delay differential equations without damping term. J Inequal Appl 2014, 361 (2014). https://doi.org/10.1186/1029-242X-2014-361

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2014-361

Keywords