Abstract
In this note we consider the classical Massera theorem, which proves the uniqueness of a periodic solution for the Liénard equation
and investigate the problem of the existence of such a periodic solution when f is monotone increasing for and monotone decreasing for but with a single zero, because in this case the existence is not granted. Sufficient conditions for the existence of a periodic solution and also a necessary condition, which proves that with this assumptions actually it is possible to have no periodic solutions, are presented.
MSC: 34C05, 34C25.
Keywords:
Liénard equation; limit cycles; Massera theorem1 Preliminaries and wellknown results
The problem of the existence and uniqueness of a periodic solution for the Liénard equation
has been widely investigated in the literature. Among the uniqueness results, a special place is, without any doubt, taken by the classical Massera theorem. This is due to the simple geometrical ideas and the fact that this result, despite several efforts, is in most cases no more valid for the generalized Liénard equation
For related results concerning the more general equation
we refer to [1].
Equation (1.1) is equivalent to the phaseplane system
and to the Liénard system
We present the classical Massera theorem which is a milestone among the results of limit cycles uniqueness for system (1.2).
Theorem 1.1 (Massera [2])
System (1.2) has at most one limit cycle which is stable, and hence equation (1.1) has at most one nontrivial periodic solution which is stable, provided thatfis continuous and
1. is monotone decreasing for,
2. is monotone increasing for.
The theorem of Massera improved a previous result due to Sansone [3] in which there was the additional assumption . This assumption comes from the fact that Sansone was using the polar coordinates. Such a strong restriction on f is clearly not satisfied in the polynomial case, and hence Massera’s result is much more powerful. We recall the recent papers [4] and [5] in which a discussion concerning these two results, as well as related results, may be found.
We observe that in his paper Massera was proving the uniqueness of limit cycles, regardless of the existence, because only the monotonicity properties and the continuity were required. It is well known that in order to fulfill the necessary conditions for the existence of limit cycles, the only two cases to be considered are the following:
A has two zeros . In this case the existence of limit cycles is granted (see, for instance, [6]). Moreover, has three zeros at , 0 and , is negative for () and positive outside this interval, and F is monotone increasing for and .
B remains negative for (or for ), while it intersects the xaxis once in (or for ).
This is because if we consider the distance from the origin, namely , and evaluate in systems (1.2) and (1.3), we get
Hence the distance from the origin increases if f is negative and decreases if f is positive, and a necessary condition for the existence of limit cycles is that changes sign.
As for [A], in a very recent paper [5], the monotonicity assumptions, which are required in the whole line, were relaxed to the interval , where α and δ can be easily determined. In particular, the following result holds.
Theorem 1.2 ([5])
The Liénard system (1.3) has exactly one limit cycle, which is stable, provided thatfis continuous and
1. if, is monotone decreasing for, is monotone increasing for;
2. if, is monotone decreasing for, is monotone increasing forwith
wherea, bare the zeros ofandα, βare the nontrivial zeros of. (If, the uniqueness result comes from classical results due to LevinsonSmith[7]and Sansone[3]and no monotonicity assumptions are required.)
Let us consider case [B]. Now the existence of the limit cycle is no more granted, and, as far as we know, this case was not explicitly investigated in the literature.
The aim of this note is to discuss such a situation. In the next section, we present some sufficient conditions for the existence, and hence the uniqueness, of the limit cycle as well as a necessary condition. Such a condition shows that actually there are situations in which there is no limit cycle despite the monotonicity properties and the fact that changes sign.
2 Existence of the limit cycle
We consider the case
B remains negative for , while it intersects the xaxis once in .
The case ‘ is negative for , while it intersects the xaxis once in ’ is treated in exactly the same way.
At first we work in the phaseplane
and we give a constructive result based on the methods used in [8].
Theorem 2.1The phaseplane system (1.2) has exactly one limit cycle, which is stable, provided thatfis continuous and
1. is monotone decreasing for,
2. is monotone increasing forand has a single zero at,
3.
4.
Proof We just notice that the continuity of f guarantees the property of uniqueness for the solutions to the Cauchy problem associated to system (1.2), and therefore the trajectories of such a system cannot intersect. Consider the 0isocline of system (1.2), that is, the points in the phaseplane were and .
For these points, we have . Let Δ be the graph of the function . Clearly, we get
and that is bounded as . Consider the slope of the trajectories of system (1.2) given by
In , is negative and the slope at the points below Δ is negative. Hence we can select a point , with and arbitrarily small, such that is smaller than all the other values of for . Let be a negative semitrajectory which starts at P. Inspection of the slopes shows that such a semitrajectory does not intersect Δ, and therefore is bounded away from the xaxis for . We now consider the positive semitrajectory which starts at P. At first we observe that being monotone, has a limit, , for . Hence, for , there exists such that for every . It is trivial to see that in the strip the slope is bounded for , and therefore intersects the line at a point with , or intersects the xaxis in . Clearly, we must consider the first situation being our goal to reach the xaxis. For , we get
This means that in , Δ lies between the lines and and hence the 0isocline is bounded away from the line . Such a line has the slope . On the other hand, in the slope of is
Therefore, if , eventually will intersect the line at some , and hence also Δ at some . A straightforward calculation yields that such an inequality holds if , and being ϵ arbitrary, this gives , that is, the value of the limit for in the assumptions. Now we recall that in the slope at the points above Δ is negative, and it is easy to see that eventually will intersect the xaxis in . In fact, any trajectory intersecting the xaxis in intersects the yaxis in . This is because the slope being bounded for , no vertical asymptote is possible. Then will intersect the xaxis in . This proves that the trajectory is winding (see Figure 1). Being , the origin is a source and the existence of a limit cycle is granted by the PoincaréBendixson theorem, while the uniqueness is given by monotonicity properties in virtue of the Massera theorem. □
Figure 1. Trajectoryin Theorem 2.1.
We observe that if we keep all the sign assumptions and the limits at infinity, but not the monotonicity, clearly the uniqueness is no more granted, but the existence result still holds, and one can relax the assumption
with
as in [8].
Hence we have the following corollary.
Corollary 2.2The phaseplane system (1.2) has at least a stable limit cycle provided thatfis continuous and, with, and
1.
2.
As already mentioned, the first result is based on the direct construction of a winding trajectory.
A more powerful approach can be adopted in the Liénard plane, and it is based on the results presented in [9], where the problem of intersection with the vertical isocline is investigated in full for the more general Liénard system
where
The following result holds.
Theorem 2.3The Liénard system (1.3) has exactly one limit cycle, which is stable, provided thatfis continuous and
1. is monotone decreasing for,
2. is monotone increasing forand has a single zero at,
3.
4.
Proof As already mentioned, the continuity of f guarantees the property of uniqueness for the solutions to the Cauchy problem associated to system (1.3), and therefore the trajectories of such a system cannot intersect.
As the proof is based on different applications of both necessary and sufficient conditions for the intersection of a semitrajectory with the vertical isocline (2.14) given in [9], we are going to introduce the same notation as in [9] and to recall the theorems we are going to use here.
As well as in the phaseplane, for any point , we denote by and the positive, respectively negative, semitrajectory of the Liénard system (2.13) passing through P. With respect to the intersection of with the vertical isocline (2.14), we observe that by a standard inspection of the direction of the field, it suffices to consider the cases
and
With respect to the intersection of with the vertical isocline (2.14), it suffices to consider the cases
and
As in [9], we define
and, finally, we introduce the function
For this particular case, we have the following result.
Theorem 2.4 ([[9], Theorem 2.5.])
If, for everyverifying (2.19), intersects the curveatwith, then
Theorem 2.5 ([[9], Theorem 2.3.])
If there is a constantsuch that
then, for everyverifying (2.18), intersects the curveatwith.
Let us see that hypothesis (2.15) implies that there exists such that (2.22) is verified. First of all, we observe that in the case we are dealing with . By hypothesis (2.15), there exists such that, for , we have . It is not restrictive to suppose . This is because we are studying the case , being the result already proved for , in virtue of the previous theorem. Then, for , we have
In our assumptions, for , . Then
A straightforward calculation yields
therefore, choosing , we can find such that for ,
At last
Then limit (2.22) is verified provided that
ϵ and η being arbitrary, the existence of such that the last inequality holds is ensured by . With a very similar calculation (with , where α is the nonnull zero of F), hypothesis (2.16) yields to prove that
As the necessary condition for the intersection of the negative semitrajectory with the vertical isocline is denied, there exists a point such that is not intersecting the curve defined in (2.14). By a standard inspection of the direction of the field, crosses the negative ysemiaxis and reaches a point verifying (2.18). By the existence of such that the sufficient condition (2.22) holds, the positive semitrajectory intersects at first the yaxis because again there are no vertical asymptotes and then the curve because the sufficient condition for intersection is fulfilled. By the uniqueness for the solutions to the Cauchy problem associated to system (1.3), trajectories cannot intersect and hence is a winding trajectory (see Figure 2). Being , the origin is a source and the existence of a limit cycle is granted by the PoincaréBendixson theorem, while the uniqueness is given by monotonicity properties in virtue of the Massera theorem. □
Figure 2. Trajectoryin Theorem 2.3.
If we relax the monotonicity assumptions, we lose uniqueness, but existence still holds. Hence we have the following corollary.
Corollary 2.6The Liénard system (1.3) has at least a stable limit cycle provided thatfis continuous andwith, and
1.
2.
Now we consider the case of nonexistence of limit cycles. We consider again the results of [9] and, in particular, the case in which the semitrajectory , starting at a point verifying (2.19), intersects the vertical isocline in the positive xhalf plane, but does not intersect the same curve in the negative xhalf plane, and we get the following necessary condition.
Theorem 2.7The Liénard system (1.3) has no limit cycles provided thatfis continuous and
1. is monotone decreasing for,
2. is monotone increasing forand has a single zero at,
3.
4.
Proof At first we need the following results.
Theorem 2.8 ([[9], Theorem 2.6.])
If there is a constantsuch that
then, for everyverifying (2.19), intersects the curveatwith.
Theorem 2.9 ([[9], Theorem 2.2])
If, for everyverifying (2.18), intersects the curve, then
As in the previous theorem, the assumptions (2.28) and (2.29) become respectively
where α is the nonnull zero of F.
With the same calculation, which we omit for the sake of simplicity, we get that if assumptions hold, we can provide a trajectory which is unwinding and goes to infinity without intersecting the curve for x negative (see Figure 3).
Figure 3. Trajectoryin Theorem 2.7.
Now assume, by contradiction, that there exists a limit cycle. The monotonicity property of the Massera theorem guarantees that such a limit cycle is unique and stable. Therefore there are trajectories surrounding the limit cycle which are winding and going to such a cycle. But with the same calculation already applied, which we omit for the sake of simplicity, if assumptions (2.26) and (2.27) hold, we can provide a large trajectory which is unwinding and goes to infinity without intersecting the curve for x negative. We can apply again the PoincaréBendixson theorem and get the existence of a second limit cycle which is unstable.
But as already mentioned, the monotonicity properties of the Massera theorem give the uniqueness, and hence in this situation we can conclude that there is no limit cycle and all trajectories are unwinding. □
We observe that there are cases which still are not investigated. For instance, if we consider such that
we can produce a trajectory which lies entirely below the curve in the Liénard plane (or below the xaxis in the phaseplane) and therefore is not winding or unwinding.
In order to describe in full the possible configurations, we consider the following example.
Example 2.10 Let f be continuous, monotone decreasing for and monotone increasing for with a single zero at .
Define and consider the system
If or , system (2.30) has no limit cycles, while if , system (2.30) has exactly one limit cycle which is stable.
This is because if , we can apply Theorem 2.7, while if , then is negative, and if , then is positive (for , is positive for ). And as already mentioned, if does not change sign, there is no limit cycle. On the other hand, if , then we can apply Theorem 2.3, and in this case there is a unique limit cycle which is stable.
We can describe the movement of such a cycle. At , changes sign and we have a standard Hopf bifurcation, which gives a small amplitude limit cycle bifurcating from the origin. When ϵ decreases, such a limit cycle enlarges and eventually collapses to a separatrix at . For a complete description of the bifurcation from a separatrix, we recall [10]. For the bifurcation of limit cycles we recall [11].
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
Both authors contributed equally in writing this paper.
Acknowledgements
We thank the anonymous referees for their valuable suggestions and Dr. Francesco Della Santa for the figures.
References

Carletti, T, Rosati, L, Villari, G: Qualitative analysis of the phase portrait for a class of planar vector fields via the comparison method. Nonlinear Anal.. 67(1), 39–51 (2007). Publisher Full Text

Massera, JL: Sur un Théoreme de G. Sansone sur l’equation de Liénard. Boll. Unione Mat. Ital.. 9(3), 367–369 (1954)

Sansone, G: Soluzioni periodiche dell’equazione di Liénard. Calcolo del periodo. Rend. Semin. Mat. (Torino). 10, 155–171 (1951)

Sabatini, M, Villari, G: On the uniqueness of limit cycles for Liénard equation: the legacy of G. Sansone. Matematiche. LXV(II), 201–214 (2010)

Villari, G: An improvement of Massera’s theorem for the existence and uniqueness of a periodic solution for the Liénard equation. Rend. Ist. Mat. Univ. Trieste. 44, 187–195 (2012)

Villari, G: On the existence of periodic solutions for Liénard’s equation. Nonlinear Anal. TMA. 7, 71–78 (1983). Publisher Full Text

Levinson, N, Smith, O: A general equation for relaxation oscillations. Duke Math. J.. 9, 382–403 (1942). Publisher Full Text

Villari, G: Periodic solutions of Lénard’s equation. J. Math. Anal. Appl.. 86, 379–386 (1982). Publisher Full Text

Villari, G, Zanolin, F: On a dynamical system in the Liénard plane. Necessary and sufficient conditions for the intersection with the vertical isocline and applications. Funkc. Ekvacioj. 33, 19–38 (1990)

Villari, G, Villarini, M: Limit cycles and bifurcation from a separatrix for a polynomial Liénard system in the plane. Differ. Equ. Dyn. Syst.. 5, 423–437 (1997)

Perko, LM: Bifurcation of limit cycles: geometric theory. Proc. Am. Math. Soc.. 114, 225–236 (1992). Publisher Full Text