W-transform for exponential stability of second order delay differential equations without damping terms

In this paper a method for studying stability of the equation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x^{\prime \prime }(t)+\sum_{i=1}^{m}a_{i}(t)x(t- \tau_{i}(t))=0$\end{document}x″(t)+∑i=1mai(t)x(t−τi(t))=0 not including explicitly the first derivative is proposed. We demonstrate that although the corresponding ordinary differential equation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$x^{\prime \prime}(t)+\sum_{i=1}^{m}a_{i}(t)x(t)=0$\end{document}x″(t)+∑i=1mai(t)x(t)=0 is not exponentially stable, the delay equation can be exponentially stable.


Introduction
In this paper we apply the W-method to obtain explicit conditions for exponential stability of linear second order delay differential equations. The method consists in a transformation of a given differential equation to an operator equation by the substitution where W (t, s) is the Cauchy function for some known exponentially stable equation.
The main object of this paper is the second order delay differential equation with a corresponding initial function defining what should be put in the equation instead of x(tτ i (t)) when tτ i (t) <  or x(tθ i (t)) when tθ i (t) < . For simplicity and without loss of generality we can consider the zero initial function Various applications of equation (.) and its generalizations can be found, for example, in the theory of self-excited oscillations, in oscillation processes in a vacuum tube, in dynamics of an auto-generator, in description of processes of in-feed grinding and cutting (see []); on position control in mechanical engineering, on electromechanical systems, and on combustion engines []. The problem of stabilizing the rolling of a ship by the activated tanks method in which ballast water is pumped from one position to another was reduced in [] to analysis of stability of the second order delay equation.
Asymptotic properties of second order equations without damping term were studied in ([], Chapter III, Section , pp.-), where instability of the equation for every pair of positive constants b and τ was obtained. Conditions of unboundedness of solutions to the equation with variable coefficients and delays were obtained in []. The condition ∞  τ (t) dt < ∞ is necessary and sufficient for the boundedness of all solutions to the equation . Results as regards the boundedness of solutions for vanishing delays (τ i (t) →  for t → ∞) and as regards asymptotic representations of solutions were obtained in [, ], see also [], Chapter III, Section . Boundedness of solutions for equations with advanced arguments (τ i (t) ≤ ) was studied in []. First results on the exponential stability for the equation x (t) + ax(t)bx(tτ ) =  with constant coefficients and delay were obtained in [-]. First results on the exponential stability of the second order equation (.) without damping term and with variable coefficients and delays were obtained recently in []. All results of [] concern only the case of a non-oscillation equation, and cannot be efficiently applied to the oscillation case. Our present paper fills this gap.
Our method reduces the study of stability of equations without damping term to the study of stability of corresponding equations with damping term. The stability of autonomous delay differential equations of the second order with damping terms was studied in [, , ], which in the case of delay differential equations apply quasi-polynomials and not polynomials as in the case of ordinary differential equations. That is why a special approach for the analysis of the characteristic equations is used []. The technique of the where x  (t), x  (t) are two solutions of homogeneous equation (.), (.) satisfying the conditions the kernel C(t, s) in this representation is called the Cauchy function (fundamental function in other terminology) of equation (.). Consider equation (.) with the following initial conditions: Let us formulate several definitions concerning stability.
Definition  Equation (.) is uniformly exponentially stable if there exist N >  and α > , such that the solution of (.), (.) satisfies the estimates where N and α do not depend on t  and initial function ϕ.
Definition  We say that the Cauchy function C(t, s) of equation (.) satisfies the exponential estimate if there exist positive N  and α  such that It is well known that for equation (.) with bounded delays these two definitions are equivalent [].
The paper is built as follows. In the first section we describe known results on asymptotic properties of second order delay equations. In Section , we formulate the main results of the paper and compare them with known results. Auxiliary assertions can be found in Section . The proofs of the assertions, formulated in Section , can be found in Section .

Formulation of main results
Let us consider the following ordinary differential equation: with constant positive coefficients A and B. We demonstrate that the exponential stability of this equation, under corresponding conditions on coefficients and delays, implies the exponential stability of delay differential equation (.). If x is a solution of equation (.) satisfying (.), we can write the equality

Let us denote by
To connect (.) and (.) we suppose that the constants A and B are such that . . , m, and the following inequality be fulfilled: then the Cauchy function C(t, s) of equation (.) and the fundamental system x  (t), x  (t) of equation (.), (.) satisfy an exponential estimate.
. . , m, and the following inequality be fulfilled: then the Cauchy function C(t, s) of equation (.) and the fundamental system . . , m, and the following inequality be fulfilled: Let us formulate corollaries for the equation f is a measurable essentially bounded function.
then the Cauchy function C(t, s) of equation (.) and the fundamental system x  (t), x  (t) of equation (.), (.) satisfy an exponential estimate.

Estimates of integrals of the Cauchy functions for auxiliary equations
Let us consider the following auxiliary equation: The proofs of these lemmas can be found in [], Lemmas .-..

Proofs of main theorems
The following result, which we formulate here for equation (.), is known as the Bohl-Perron theorem.

Lemma  Suppose the solution of equation (.) with the zero initial conditions,
is bounded on [, ∞) for any essentially bounded on the semi-axes right-hand side f . Then equation (.) is uniformly exponentially stable.
Let us consider the ordinary differential equation with constant positive coefficients A and B. It is clear that the solution of equation (.) which satisfies the initial conditions can be written in the form where W (t, s) is the Cauchy function of equation (.). Its derivatives are the following: Theorem  Let A > , B > , τ i (t) ≤ θ i (t) and the following inequality be fulfilled: Proof Consider the equation Let us rewrite equation (.) in the following forms: and We have to prove that, for every essentially bounded function f (t), the solution x(t) is also bounded on the semiaxis t ∈ [, +∞). To prove exponential stability of (.) we assume the existence of an unbounded solution x(t) and demonstrate that this is impossible.
There exist measurable delay functions η i , τ i (t) ≤ η i (t) ≤ θ i (t) such that equality (.) can be written as and then in the form Let us rewrite it in the form where the operator K : L ∞ → L ∞ is defined by the equality and The inequality (.) implies that the norm K of the operator K : L ∞ → L ∞ is less than one. In this case, there exists the bounded operator (I -K) - . For every bounded righthand side the solution z of equation (.) is bounded.
In the case A >  and B > , the Cauchy function W (t, s) and its derivative W t (t, s) satisfy exponential estimates. The boundedness of the solution x of equation (.) and its derivative x follow now from the boundedness of z.
We have got a contradiction with our assumption that the solution x(t) is unbounded on the semiaxis.
Thus solutions x(t) are bounded on the semiaxis for a bounded right-hand side f (t). Then by Lemma  the Cauchy function C(t, s) of equation (.) and the solutions x  and x  satisfy the exponential estimates.
This completes the proof of Theorem .
To prove Theorems - we set the norms of W , W t and W tt obtained in Lemmas - into Theorem .
The proofs of Corollaries - are results of substitution of W , W t and W tt into Theorems -, when we take into account that

Conclusion
In this paper to obtain exponential stability conditions we use the substitution where Y (t, s) is the fundamental function of exponentially stable autonomous ordinary differential equation of the second order and then analyze the operator equation in some functional Banach spaces on semi-axes. This method is usually called the Wmethod and used in many problems for FDE such as stability, oscillation and nonoscillation, and boundary value problems. To apply this method in our research we need integral estimations of a fundamental function and its first and second derivative. Such estimates were obtained recently and allow us to obtain here new stability results for a wide class of delay differential equations of the second order. In particular we obtain new stability results for delay equations with positive and negative coefficients. Such equations without delay are unstable so the results obtained here one can use to stabilize mathematical models described by ordinary differential equations of the second order.

Discussion and some topics for future research
The idea of applications of W-transform for second order delay differential equations first appeared in paper []. Lemmas  and  were taken from []. All other results in the paper are new and have not been published before.
Explicit integral estimates of the fundamental function and its derivatives we obtain here only for the simplest equation: ordinary differential equation with constant coefficients. It is interesting to obtain such estimates for the ordinary differential equation with variable coefficients or for delay differential equation with constant coefficients. It will allow one to improve the results obtained in this paper.
We did not study here nonlinear equations. It is interesting to consider the Wtransformation method for nonlinear equations such thaẗ x(t) + f t,ẋ r(t) + G t, x h(t) = .
The next problem is to apply the W-transformation method for the instability of delay differential equations of the second order.
We suppose that the W-transformation method can also be applied for a vector delay differential equation

X(t) + A(t)Ẋ g(t) + B(t)X h(t) = ,
where A and B are n × n matrix-functions. To apply this method it is necessary to obtain explicit integral estimations for the fundamental matrix and its derivative of the following vector ordinary differential equation with constant coefficients: