DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS
DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS
DIFFERENtIAl & DIFFERENCE EqUAtIONS ANd APPlICAtIONS
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
LIMIT CYCLES OF LIÉNARD SYSTEMS<br />
M. AMAR AND Y. BOUATTIA<br />
We calculate the amplitude and period of the limit cycle of the following classes of Liénard<br />
equations by using the method of Lopez and Lopez-Ruiz: ẍ + ε(x 2 − 1)ẋ + x 2n−1 = 0; ẍ +<br />
ε(x 2m − 1)ẋ + x 2n−1 = 0; ẍ + ε(|x| n − 1)ẋ + sign(x) ·|x| m = 0; ẍ + ε(x 2m − 1)ẋ +<br />
x 1/(2n−1) = 0, where m,n ∈ N.Wegivenumericalresults.<br />
Copyright © 2006 M. Amar and Y. Bouattia. This is an open access article distributed under<br />
the Creative Commons Attribution License, which permits unrestricted use, distribution,<br />
and reproduction in any medium, provided the original work is properly cited.<br />
1. Introduction<br />
Limit-cycle behavior is observed in many physical and biological systems. The problem<br />
of determining when a nonlinear dynamical system exhibits limit cycle has been of great<br />
interest for more than a century. Limit cycles cannot occur in linear systems, conservative<br />
systems, and gradient systems. The limit cycles are caused by nonlinearities. It was found<br />
that many oscillatory circuits can be modeled by the Liénard equation<br />
where<br />
ẍ + f (x)ẋ + g(x) = 0, (1.1)<br />
·= d dt . (1.2)<br />
It can be interpreted mechanically as the equation of motion for a unit mass subject<br />
to a nonlinear damping force − f (x)ẋ and a nonlinear restoring force −g(x). Applications<br />
of Liénard’s equation can be found in many important examples including chemical<br />
reactions, growth of a single species, predator-prey systems, and vibration analysis.<br />
In Section 2, we give the theorem of existence and uniqueness of the limit cycle for the<br />
Liénard equation. In Section 3, we give the method of calculation of the amplitude of the<br />
limit cycle of the perturbed centers (see [2]),<br />
ẋ =−y 2l−1 + εP(x, y),<br />
ẏ = x 2k−1 + εQ(x, y),<br />
(1.3)<br />
Hindawi Publishing Corporation<br />
Proceedings of the Conference on Differential & Difference Equations and Applications, pp. 743–755