Asymptotic Stability of Solutions to a Class of Second-Order Delay Differential Equations

: We consider a class of second-order nonlinear delay differential equations with periodic coefﬁcients in linear terms. We obtain conditions under which the zero solution is asymptotically stable. Estimates for attraction sets and decay rates of solutions at inﬁnity are established. This class of equations includes the equation of vibrations of the inverted pendulum, the suspension point of which performs arbitrary periodic oscillations along the vertical line.

In the case γ = σ = 0 (1) coincides with the equation of the "inverted pendulum", (ϕ(t) = sin t, ψ(y) = sin y), with the oscillating suspension point along the vertical line according to the law x(t) = a sin(ωt) [1,2]. In the literature, this pendulum is often called the Kapitsa pendulum. (1) describes oscillations of a "controlled inverted pendulum", the suspension point of which performs arbitrary periodic oscillations.
It should be noted that a large number of works are devoted to the study of inverted pendulums with vibrating suspension points, pendulums on carts, systems consisting of several pendulums, with control, etc. (for example, refer to [3,4] and the bibliography therein). One of the main problems is the stability problem for these mechanical systems. Many papers are devoted to finding stability conditions, but very few papers contain estimates characterizing stabilization rates to equilibrium positions.
The present paper continues a cycle of authors' research [5][6][7][8][9][10][11] on the stability of solutions to ordinary differential equations and delay equations with periodic coefficients in linear terms. In [5], the following systems of linear differential equations were considered: where A(t) is an n × n matrix with continuous T-periodic entries. The following criterion for the asymptotic stability of the zero solution to (3) was established in terms of solvability of a special boundary value problem for the Lyapunov differential equation Theorem 1 ([5]). The zero solution to (3) is asymptotically stable if and only if there is a Hermitian solution to the boundary value problem: where the Hermitian matrix C(t) with continuous entries on [0, T] is positive definite.
Here and further, S > 0 (S < 0) denotes that S is a positive (negative) definite Hermitian matrix. We use the spectral matrix norm everywhere. Note that differential equations of the form (4) have been considered in many papers when studying qualitative behavior of solutions to differential equations with variable coefficients (for example, see [12][13][14][15][16][17]).
Using the solution H(t) to the boundary value problem (5), we can point out the radius of the disk containing all the multipliers of (3) and obtain estimates characterizing the decay rate of solutions to (3) at infinity [5]. Theorem 2 ([5]). If the zero solution to (3) is asymptotically stable, then the following estimate holds for its solutions: where H(t) is a T-periodic extension of the solution to (5) on the semi-axis {t ≥ 0}, c 1 (t) > 0 is the minimal eigenvalue of the matrix C(t).

Remark 2.
Theorems 1 and 2 allow us to study the asymptotic stability of the zero solution to (3) without finding multipliers.
Using the results mentioned above, the authors investigated the asymptotic stability of the zero solution to the following nonlinear systems [6]: where f (t, x) is a smooth real-valued vector-function, f (t, 0) = 0. Using the solution H(t) to the boundary value problem (5), estimates for attraction sets of the zero solution were obtained and estimates of exponential decrease of solutions to (7) at infinity were established. On the basis of these studies, introducing special Lyapunov-Krasovskii functionals, the authors investigated the stability of solutions to some classes of delay differential equations with periodic coefficients in linear terms (for example, see [7][8][9][10][11]). Conditions for the asymptotic stability of the zero solution and estimates characterizing decay rates of solutions as t → ∞ were obtained in these papers. These estimates are analogs of Krein's estimate.
We will use some assertions that follow from [7,8,11]. Consider the following system of delay differential equations: where A(t), B(t) are n × n matrices with continuous real-valued T-periodic entries, τ > 0 is a delay and F(t, u, v) is a continuous real-valued vector-function on [0, ∞) × R n × R n . We assume that F(t, u, v) satisfies the Lipschitz condition with respect to u on every compact set G ⊂ [0, ∞) × R n × R n and the following estimate: for some constants p ≥ 0, δ > 0. At first we formulate our results for the linear systems.
such that and the following differential inequality holds. Then the zero solution to (8) is asymptotically stable.
We introduce the following matrix: If the conditions of Theorem 3 hold, then P(t) > 0. Denote by p 1 (t) > 0 the minimal eigenvalue of the matrix P(t) and by k > 0 the maximal number such that the following inequality holds: Consider the initial value problem for (8): ) is a given vector-function.
Suppose that there are matrices H(t) and K(s) satisfying the conditions of Theorem 3. Then the following estimate holds for the solution to (12): where h 1 (t) > 0 is the minimal eigenvalue of the matrix H(t).
Inequality (13) is an analog of Krein's estimate. Consider nonlinear systems of the form (8). Supposing that the conditions of Theorem 3 hold, we define R as follows: Theorem 5 ([7,11]). Let the conditions of Theorem 3 hold. Then, the zero solution to (8) is asymptotically stable and the set of real-valued functions is an attraction set of the zero solution. Moreover, the solution to the initial value problem (12) with x 0 (t) ∈ E satisfies the following estimate: Using the results mentioned above, auxiliary transforms are performed and preliminary assertions are proved in Section 2. The main results are formulated and proved in Section 3.

Auxiliary Results
Consider the following linear equation: where the function ϕ(s) satisfies (2). In the next section we prove that the zero solution to (14) is asymptotically stable under certain conditions on the function ϕ(s) and parameters α, β, γ, σ, and ω. In order to prove this result, we need to perform a number of transforms. Using the following change of variables: we rewrite (14) as follows: We introduce the vector-function: This equation is equivalent to the system of linear delay differential equations: where We introduce the following functions: Then, (15) is reduced to the following system: where the entries of A(t, µ) have the following form: , .
We show below that there are T-periodic functions a(t), b(t), c(t) such that the vectorfunction u(t) = u 1 (t) u 2 (t) satisfies the following system: By analogy with [6], we take c(t) ≡ 1 and choose the functions a(t), b(t) such that the following equalities hold: Hence, we have (17), where Using the following representation: we rewrite A(t, µ) as follows: We show below that there are the functions a(t), b(t) satisfying (18), such that the spectrum of the matrix A 1 belongs to the left half-plane, the matrix A 2 (t, µ) is T-periodic, and the following equality holds: Taking into account (18), we seek the functions a(t), b(t) in the following form: Clearly, t 1 , t 2 ∈ [0, T] must be such that the following equalities hold: At first we construct the function b(t). Obviously, by the mean value theorem, there exists t 1 ∈ [0, T] such that the following equality holds: Taking into account the following equalities: Take the function b(t) in the form (23) with t 1 defined above. By (2), the function b(t) is T-periodic; moreover, from (26), we have We now construct the function a(t). At first, we show that there exists t 2 ∈ [0, T] such that the relation holds: Indeed, we can represent the right-hand side in the following form: Obviously, for the first summand, the equality is valid: By the mean value theorem, the second summand can be written as follows: Indeed, integrating by parts, we have Consequently, (27) is equivalent to (28). Therefore, the function a(t) can be taken in the form (24) with t 2 defined above. We now show that a(t) is T-periodic. Indeed, according to the definition, we have the following identities: On the other hand, by the properties of the function b(t), we have Hence, a(t + T) ≡ a(t). Thus, we have specified the functions a(t) and b(t) satisfying (18) such that the matrix A 2 (t, µ) is T-periodic and the following equality holds: Proof. From the explicit form (20) of the matrix A 1 , it follows that its spectrum belongs to the left half-plane if and only if α > 0 and We show below that this inequality is equivalent to (30). To perform this, we make a series of transforms of the integral in (31). Integrating by parts and taking into account (2), we have By (24), this equality can be written as follows: Using (25), we obtain Taking into account (2), we have The equivalence of (30) and (31) follows immediately from this formula. The lemma is proved.

Main Results
At first we consider (14). We introduce the following notation: Taking into account the transforms performed in the previous section, the matrix A 1 given in (20) has the following form: If α > 0, β < 0, then its spectrum belongs to the left half-plane. In this case the following Lyapunov equation: H 1 A 1 + A * 1 H 1 = −I has a unique solution H 1 = H * 1 > 0; moreover, it can be written as follows (for example, see [13]): We establish below some conditions for the coefficients of (14) under which it is possible to guarantee the asymptotic stability of the zero solution to (14).
Proof. We will use the transforms performed in the previous section. According to the notation, µ = 1 ω . We assume that 0 < µ < . By (24), this condition guarantees that 1 + µa(t) > 0. Taking into account these transforms, in order to prove the theorem, it is sufficient to establish the asymptotic stability of the zero solution to the linear system of delay differential equations: where the matrices A 2 (t, µ), A 3 (t, µ), and B(t, µ) are defined in (19), (21), and (22), respectively. According to Theorem 3, if there are matrices H(t, µ) and K(s) satisfying (10) and (11) and such that the following matrix inequality holds: then the zero solution to (32) is asymptotically stable. It is not difficult to calculate the spectral norm of the matrix B(t, µ) given in (19); moreover, Using (23) and (24), we have the following estimates: Hence, by the conditions of the theorem, we obtain We take k > 0 such that ρ = qe k τ/2 < 1. Consider the following matrix: where H 1 is the above solution to the Lyapunov equation, Obviously, H(0, µ) = 1 µ H 1 > 0 and, by (29), we have H(T, µ) = H(0, µ). Clearly, where By construction, there is , such that the following inequality holds: for 0 < µ ≤ µ 1 . Then there exists µ 0 , 0 < µ 0 ≤ µ 1 , such that the following estimate is valid: for µ ≤ µ 0 . Hence, the matrix H(t, µ) is a Hermitian solution to the problem of the form (5).
In this case, the condition β < 0 has the following form: Using Theorem 4, we have estimates characterizing decay rates of solutions to (32) at infinity. Making the inverse change of variables, we obtain asymptotic estimates for solutions to (15) and for solutions to (14).
Hence, (42) is included in the class of systems of the form (8) with the right-hand side satisfying (9). Then, by Theorem 5, we obtain an attraction set for the zero solution to (1) and the estimates of the exponential decrease of solutions to (1) with the initial data from this attraction set. and then the zero solution to (1) is asymptotically stable for sufficiently small µ > 0.

Remark 3.
As was mentioned above, in the case γ = σ = 0, (1) coincides with the equation of the "inverted pendulum" (ϕ(t) = sin t, ψ(y) = sin y) with the oscillating suspension point along the vertical line by according to the law x(t) = a sin(ωt) [1,2]. The conditions of the asymptotic stability of the form (43) coincide with known N.N. Bogolyubov's conditions obtained by the averaging method (for instance, refer to [18,19]).