In the present paper, we have investigated the non Liénard system. We have shown that limit cycles may bifurcate at the origin. Bendixons theorem has been used in our study to prove non-existence of limit cycles. We have also proved that the system has unique limit cycle through change of the parameters.
In the present investigation, we revisit the problem of bifurcation of limit cycles. The problem of limit cycle was studied intensively. For Liénard, we can read [
We give criterion for the non Liénard system to have or not to have limit cycles with some parameters. We also demonstrate that the system exhibits a Hopf-bifurcation. Now we consider the following Liénard equation
x ¨ + f ( x ) x ˙ + g ( x ) = 0. (1)
The above equation may be written in two dimensional autonomous dynamical system
x ˙ = y , y ˙ = − g ( x ) − f ( x ) y . (2)
Therefore, the above equations can be written in the Liénard plane as
x ˙ = y − F ( x ) , y ˙ = − g ( x ) , (3)
where F ( x ) = ∫ 0 x f ( t ) d t .
Theorem 1.1 [
Theorem 1.2 [
In Section 2, the main system equations results have been presented, the section has been divided in two cases.
The case I considered the conditions that the system has a limit cycle when O ( 0,0 ) is an anti saddle.
Finally, the case of saddle point with limit cycle is presented in theorems and lemmas in Section 4 along with the concluding remakes.
The main part of this paper is devoted to explain the existence and uniqueness of limit cycles of the following differential equations system
x ˙ = − x + a y , y ˙ = b x − a y − x 2 y , (4)
the singular points of the system are x = 0 and x = ± a − a b .
The Jacobian matrix
A = [ − 1 a b − a ]
has the determinant | A | = a − a b for a − a b ≥ 0 the origin O is anti saddle and for a − a b < 0 the origin O is saddle for more details (see [
The system (2.2) needs to change to the Liénard system (1.1).
Let z = − x + a y so z ˙ = − x ˙ + a y ˙ after simplify and substitute a y = x + z so that we have
x ˙ = z , z ˙ = − ( 1 + a + x 2 ) z − ( ( a − a b ) x + x 3 ) . (5)
After change z to y we can get system (1.1) as follows
x ˙ = y − ( ( 1 + a ) x + 1 3 x 3 ) , y ˙ = − ( ( a − a b ) x + x 3 ) . (6)
The system (2.4) is considered in two cases.
Case I: The Origin is an anti Saddle
The case under consideration is a − a b ≥ 0 , in this case and as above the system (2.4) has unique equilibrium point O ( 0,0 ) which is an anti saddle.
Lemma 2.1
For a = 0 system (2.4) has no limit cycle.
Proof
Let a = 0 , then F ( x ) = x + 1 3 x 3 and g ( x ) = x 3 . Thus g ( x ) F ( x ) = x 4 ( 1 + 1 3 x 2 ) ≥ 0 and by using theorem (1.2) there is no limit cycle so we just look for a ≠ 0 .
In the case of a < 0 O becomes saddle also for a 2 − 4 > 0 or ( a > 2 ) O is node in two cases no limit cycles surround O. Thus, in the sequel, we only need to consider 0 < a < 2 .
Consider the polynomial Liénard system of degree n
x ˙ = y − ( a 1 x + a 2 x 2 + ⋯ + a m x n ) , y ˙ = − ( b 1 x + b 2 x 2 + ⋯ + b k x k ) . (7)
Lemma 2.2 [
For system (2.5) with b 1 = 1 , the first three focal values at O ( 0,0 ) are
η 2 = − a 1 , η 4 = 1 8 ( 2 a 2 b 2 − 3 a 3 ) ,
η 6 = c 0 ( 6 a 2 a 4 + 20 a 4 b 2 − 15 a 3 b 3 − 15 a 5 ) .
where c 0 is positive constant.
By scaling x → x a − a b and t → t a − a b [where new x = x a − a b , x ˙ = x ˙ and y ˙ = a − a b y ˙ ], then system (2.4) becomes
x ˙ = y − ( a − 1 a − a b x + 1 3 ( a − a b ) 3 a − a b x 3 ) , y ˙ = − x ( 1 + 1 ( a − a b ) 2 x 2 ) . (8)
Therefore the three focal values of O ( 0,0 ) and by using Lemma 2.2 namely are
η 2 = − 1 + a a − a b , η 4 = − 1 ( a − a b ) a − a b ,
η 6 = − 5 c 0 ( a − a b ) 2 a − a b .
If a ≠ − 1 then O, is strong focus which is unstable for a < − 1 and stable if a > − 1 , and for a = 1 , then O is weak focus of order one which is stable.
By using Hopf-bifurcation (by changing of stability), for a > − 1 no limit cycle because no change of stability if a = − 1 , then O is weak focus of order one which is stable. Thus as a decreasing from −1 O becomes unstable and one stable limit cycle appears from Hopf-bifurcation.
Theorem 2.3
For 1 < a < 2 the system (2.4) has a unique stable limit cycle.
Proof:
Now we apply theorem (1.1) consider g ( x ) = ( a − a b ) x + x 3 since a − a b > 0 So g ( x ) has only one root which is x = 0 . For
F ( x ) = 1 3 x 3 + ( a + 1 ) x the roots are a 1 = − − 3 ( a + 1 ) < 0 < a 2 = − 3 ( a + 1 ) . The roots of f ( x ) are − − ( a + 1 ) < 0 < − ( a + 1 ) and f ( x ) has minimum at ( 0, a + 1 ) .
Since we have − − 3 ( a + 1 ) < − − ( a + 1 ) < 0 < − ( a + 1 ) < − 3 ( a + 1 ) , then we deduce that f ( x ) ≥ 0 for x ∉ ( a 1 , a 2 )
g ( x ) F ( x ) = x 2 ( a − a b + x 2 ) ( a + 1 + 1 3 x 3 ) since a − a b > 0 then the term a − a b + x 2 > 0 and the value of ( a + 1 + 1 3 x 3 ) < 0 in the interval ( − − 3 ( a + 1 ) , − 3 ( a + 1 ) ) so g ( x ) F ( x ) ≤ 0 for x ∈ ( a 1 , a 2 ) . Finally since a 1 = − a 2 so we have G ( a 1 ) = G ( a 2 ) .
Case II: The Origin is a saddle
In this case, we discuss system (2.4) when a − a b < 0 and as above the system has three equilibrium points O ( 0,0 ) and ± α where α = a b − a trance the ( α ,0 ) to the Origin by the relation x → ( x − α )
x ˙ = y − ( ( 1 + a b ) x − ( a b − a ) x 2 + 1 3 x 3 ) , y ˙ = − x ( 2 ( a b − a ) − 3 a b − a x + x 2 ) . (9)
Let t = − τ , y → − y , (2.4) is converted into
x ˙ = y − ( − ( 1 + a b ) x + ( a b − a ) x 2 − 1 3 x 3 ) , y ˙ = − x ( 2 ( a b − a ) − 3 a b − a x + x 2 ) . (10)
η 2 = a b + 1 2 ( a b − a ) , η 4 = − 1 4 ( a b − a ) a − a b ,
η 6 = − 5 8 ( a b − a ) 2 ( a b − a ) .
By using Hopf-bifurcation, for a b + 1 < 0 no limit cycle because no change of stability if a b + 1 = 0 , then O becomes weak focus of order one which is stable.
Thus for fixe b a * = − 1 b is bifurcate value so as a * increasing, O becomes unstable and one stable limit cycle appear from Hopf-bifurcation.
Lemma 3.1
a b + a + 2 > 0 equivalent to a b + 1 > 0 .
Proof
Assume that a b + 1 < 0 since a + 1 < 0 , then we have a b + a + 2 = a b + 1 + a + 1 < 0 contradiction. Thus for a b + a + 2 > 0 also we get a b + 1 > 0 .
Lemma 3.2 [
If there exists a constant m ≥ 0 such that F ′ ( x ) G ( x ) − m F ( x ) g ( x ) ≥ 0 for x ( ≠ 0 ) , System (2.8) has at most one limit cycle.
We have by putting c = a b + 1 , α = a b − a and after simplify we have
ϕ ( x , m ) = ( m 3 − 1 4 ) x 4 + ( 3 2 α − 2 m α ) x 3 + ( ( c + 11 3 α 2 ) m − c 4 − 3 α 2 ) x 2 ( c α + 2 α 3 − ( 3 c α + 2 α 3 ) m ) x + 1 2 c α 2 .
Let m = 3 4 so we have
ϕ ( x , 3 4 ) = ( 1 2 c − 1 4 α 2 ) x 2 + ( 1 2 α 3 − 1 4 c α ) x + 1 2 c α 2
Δ = 1 4 ( a b − a ) 3 + 1 4 ( a b + 1 ) ( a b − a ) 2 − 15 16 ( a b + 1 ) 2 ( a b − a ) .
Since ( a b − a ) > 0 , we can delete from upper equation and for suitable a as small enough we have
Δ = 1 4 ( a b − a ) 2 + 1 4 ( a b + 1 ) ( a b − a ) − 15 16 ( a b + 1 ) 2 < 0 .
A non-Liénard system is studied and analyzed by adapting Hopf-bifurcation theory. It has been proved that the system has unique limit cycle under some change of parameters under two cases. Bendixons theorem is used to prove non-existence of limit cycles.
I would like to express my thanks to Prof. V.P. Sing, Depr. of Mathematics Albaha University for his voluble suggestion for improving the paper.
The authors declare no conflicts of interest regarding the publication of this paper.
Saeed, A.E.M. and Hamid, A.B. (2019) On Non Liénard Type with One Limit Cycle. Journal of Applied Mathematics and Physics, 7, 3031-3036. https://doi.org/10.4236/jamp.2019.712213