New Stability Results for Periodic Solutions of Generalized Van der Pol Oscillator via Second Bogolyubov’s Theorem

: A certain class of nonlinear differential equations representing a generalized Van der Pol oscillator is proposed in which we study the behavior of the existing solution. After using the appropriate variables, the ﬁrst Levinson’s change converts the equations into a system with two equations, and the second converts these systems into a Lipschitzian system. Our main result is obtained by applying the Second Bogolubov’s Theorem. We established some integrals, which are used to compute the average function of this system and arrive at a new general condition for the existence of an asymptotically stable unique periodic solution. One of the well-known results regarding asymptotic stability appears, owing to the Second Bogolubov’s Theorem, and the advantage of this method is that it can be applied not only in the periodic dynamical systems, but also in non-almost periodic dynamical systems.


Introduction
The mathematical model for the system of the Van der Pol equation is a well-known second-order ordinary differential equation with cubic non-linearity.The Van der Pol oscillator is a classical example of a self-oscillatory system and is now considered as a very useful mathematical model that can be used in much more complicated and modified systems.Examples of systems in which the generalized Van der Pol oscillator have recently been proposed are included in [1,2].In fact, from the point of view of practical implementation, these systems are simpler than the class of systems based on nonlinear differential equations.From a mathematical point of view, they are more complicated due to the presence of nonlinear means.Accurate mathematical analysis of the nature of attractors in such systems, including a rigorous substantiation of the non-linearity hypothesis, is a difficult problem that requires the development of new approaches.In [3], a system based on a non-autonomous Van der Pol oscillator with a delay is considered, which is alternately in the excitation and decay modes due to the periodic change in the parameter responsible for the bifurcation of the birth of a limit cycle.The excitation of oscillations at each new stage of activity is stimulated by the signal generated at the previous stage of activity, which enters through the delay line and undergoes a non-linear transformation using an auxiliary external signal at the main operating frequency, so that at each successive stage of activity, the phase variable is multiplied by a factor of 2. As a result, for the phases of successive trains of oscillations, a stretching-mapping of a circle (Bernoulli mapping) with chaotic dynamics is realized.Another scheme was proposed in [4], the advantage of which is that it operates in the absence of an auxiliary signal.The system contains two delayed feedback circuits with different delay times.In these circuits, the signal undergoes a nonlinear transformation, respectively, into the second and third harmonics.When they are mixed, a signal appears on a quadratic nonlinear element, which is used to stimulate the excitation of the oscillator at the next stage of its activity.The disadvantage of this scheme is the need to use a high degree of nonlinear transformations, which can be difficult in practice.We now proceed to discuss on the Second Bogolyubov's Theorem studied in [5][6][7][8] and also on the Van der Pol equations.For example, in [8], both are studied under the form u + (u 2 − 1)u + (1 + β )u = sin(s), where > 0, s ≥ 0 and β, are real constants.It is proven that the unique periodic solution is asymptotically stable by using the Second Bogolyubov's Theorem.In addition, the authors showed that the amplitude of ( 1) is |K| > √ 2 and the amplitude of ( 2) is |K| > π 2 .In our paper, we proved some results by using the Second Bogolyubov's Theorem to show the existence, uniqueness, and asymptotic stability for the periodic solution of the oscillatory Van Der Pol system where β, , c 1 , c 2 , c, α, c 3 ∈ R, ξ, ν ∈ N and 0 < < 1.This article introduces and analyzes new stability results for periodic solutions of the generalized Van der Pol oscillator via the Second Bogolyubov's Theorem.It consists of three sections.In the first one, we presented the general framework of our study and also introduced our main system.In the second section, we recalled most of the preliminary material and we presented the Second Bogolyubov's Theorem, which will be the key of our proofs.Finally, in Section 3, we stated and proved our results regarding the existence, uniqueness, and asymptotic stability of nonlinear second-order differential equations in the general form.Our study, of course, concluded in Section 4 by outlining the difficulties we faced, as well as novelties and perspectives.In addition, we stated the history of the Van der Pol equation.

Let
x = g(s, x, ), > 0, (5) where the function g periodic for t and locally Lipschitz with respect to x.We prepare to use the Second Bogolyubov's Theorem for (5) and give the next notation In the case that g is of class C 1 , we note the periodic case of the Second Bogolyubov's Theorem, which represents a part of the averaging principle: det(g 0 ) (v 0 ) = 0 and assures the existence and uniqueness, for > 0 small, of a T -periodic solution of system (5) in a neighborhood of v 0 , while the fact that all the eigenvalues of the Jacobian matrix (g 0 ) (v 0 ) have a negative real part and also provides its asymptotic stability.We will search for the periodic solution in a neighborhood of v 0 ∈ g −1 0 (0), with g 0 ∈ C 1 .Here, we consider The next theorem, Theorem 1, presents the periodic case of the Second Bogolyubov's Theorem, see [8], which is based on the Jacobian matrix of the average function in the vicinity of the value v 0 , where det(g 0 ) (v 0 ) > 0 and trace(g 0 ) (v 0 ) < 0 confirm the existence and uniqueness of the T-periodic solution of the system (5); > 0 is small enough, so it is asymptotically stable if det(g 0 ) (v 0 ) < 0.Then, the system (5) contains at least one non T-periodic solution that is asymptotically unstable.
Remark 1.We notice here that the first condition (1) of Theorem 1 ensures the existence of a periodic solution, while the second one (2) shows that the solution is unique.Condition (3) gives the asymptotic stability.

Main Results and Proofs
In this section, we prove the main results in the present paper, that both Equations (3) and (4) have a unique periodic and asymptotically stable solution.We present these main results in Theorems 2 and 3.
The first claim

2.
The second claim
We are going to prove Theorem 2, both for this end and to be ready to apply to the Second Bogolubov's Theorem; as in [9], we utilize Levinson's changes in order to rewrite (3) as a coupled system.
Proof.We have Using Proposition 1, we have ∀ξ, ν ∈ N Then It is not difficult to prove the other integrals with the same method.
By substituting these results in (13), we obtain the next Corollary.
We are now ready to prove Theorem 2.
Proof of Theorem 2. With the functions in (10), we have (i) and (ii) of claim (2) of Theorem 1 verified for Ω = R 2 .We will treat the question of existence of only one limit cycle (one periodic solution).We should now check (7); for this end, the Jacobian matrix g 0 (x 1 , x 2 ) is given by and Regarding det(B), trace(B), let and Then and We discuss the different cases or and Now, we check condition g 0 (v 0 ) = 0 in Theorem 1.
and also Then Example 1.As in [8], the Van der Pol equation is considered The existence and stability of the periodic solution was given.The results in our theorem, Theorem 2, show that (15) has exactly one cycle limit that is asymptotically stable if ) has a unique periodic asymptotically stable solution.
To prove Theorem 3, as in [9], we will use Levinson's changes to apply Bogolubov's Theorem.Let (z 1 , z 2 ) = (u, u ), so Equation (4) becomes Then , with the help of the Lipschitz property for x in the function g.The system (16) becomes where with By (6), we find the average function g 0 as The next Lemma will be useful to find the average function.
We are now in position to prove Theorem 3.
Proof of Theorem 3. Condition (i) for claim (2) of Theorem 1.In the case Ω = R 2 is verified, since the function in ( 17) is an absolute value function with respect to (x 1 , x 2 ).We will prove the existence of only one limit cycle (one periodic solution).We check conditions (7) of Theorem 1 and then demonstrate that g 0 (v 0 ) = 0 by using Corollary 2.
Let the Jacobian matrix g 0 (x 1 , x 2 ) given as and We will find det(B) and trace(B).Let u = −α

3.
We have overcome the difficulty of having a polynomial function and treated it with the Jacobian matrix, Lipschitz property, and some useful integrals.4.
We improved the early results in the literature.
This type of equation dates back to the famous Dutch scientist Balthazar Van der Pol, which was one of his best works with a major contribution to the development of some branches of modern mathematics and physics, as well as radio engineering to be more precise.His work was mainly related with the equation that bears his name and it has a surprisingly wide range of applications in daily life and the natural sciences.The equation discusses modeling processes in the human body in a good and basic way that includes most studies related to the internal organs of the human body and the nervous system, in addition to the vocal cords.The well-known result regarding the asymptotic stability of such models arises due to the Second Bogolobov's Theorem, and the advantage of this method is that it can be applied not only in periodic dynamical systems, but also in almost nonperiodic dynamical systems, see [11][12][13][14][15][16].
Open problem.It is interesting to determine where the methods in this paper can be extended and also be applicable for other finite time stability analysis like "Finite-Time Synchronization of Quantized Markovian-Jump Time-Varying Delayed Neural Networks via an Event-Triggered Control Scheme under Actuator Saturation".This open problem will be our next project.