<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.37115</article-id><article-id pub-id-type="publisher-id">AM-19865</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Limit Cycle Bifurcations in a Class of Cubic System near a Nilpotent Center
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>iao</surname><given-names>Jiang</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Mathematics, Shanghai Maritime University, 1550 Haigang Avenue in New Harbor City, Shanghai 201306, PR China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jiaojiang08@yahoo.cn</email></corresp></author-notes><pub-date pub-type="epub"><day>30</day><month>07</month><year>2012</year></pub-date><volume>03</volume><issue>07</issue><fpage>772</fpage><lpage>777</lpage><history><date date-type="received"><day>April</day>	<month>25,</month>	<year>2012</year></date><date date-type="rev-recd"><day>June</day>	<month>4,</month>	<year>2012</year>	</date><date date-type="accepted"><day>June</day>	<month>11,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper we deal with a cubic near-Hamiltonian system whose unperturbed system is a simple cubic Hamiltonian system having a nilpotent center. We prove that the system can have 5 limit cycles by using bifurcation theory.
 
</p></abstract><kwd-group><kwd>Near-Hamiltonian System; Nilpotent Center; Hopf Bifurcation; Limit Cycle</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In the International Congress of Mathematics held in Paris in 1900, Hilbert made a list of 23 problems. The second part of Hilbert’s 16th problem is still an open and difficult question: to find a upper bound of the number of limit cycles and their relative locations in polynomial vector fields of order n.</p><p>If the singular point of the system is a non-saddle, nor nilpotent, the related Hopf bifurcations are elementary, see [1-3] and their references. Hopf bifurcations from the elementary focus type of singularities have found broad and important applications in biology, chemistry and physics and engineering, see [4-7] for examples. Yet for the bifurcation of limit cycles from a non-elementary center in a more general planar vector field, its intrinsic dynamics is still far away from understanding due to the complexity and technical difficulties in dealing with such bifurcations.</p><p>Then it was natural to restrict the study of the nilpotent center by assuming the system is a perturbation of a Hamiltonian system. Consider the following system</p><disp-formula id="scirp.19865-formula38073"><label>(1.1)</label><graphic position="anchor" xlink:href="15-31777\cd14862d-3847-4af4-8f9d-02122a371cf3.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="15-31777\4e3e3bcb-8a73-46bc-ab01-3be063e56c33.jpg" />, <img src="15-31777\ea638df8-2539-426f-9b7a-d0a19c113cd9.jpg" />and <img src="15-31777\f7057bb3-207d-410a-a64c-2141561b64af.jpg" /> are <img src="15-31777\f33f5e99-b41f-41fa-a70a-61cbb24300d7.jpg" /> functions, <img src="15-31777\9bea8729-ab6a-438c-810b-7639e4d46143.jpg" />is small and <img src="15-31777\d7e57722-b304-46c4-a9bf-25349abee11a.jpg" /> with D a compact set.</p><p>When<img src="15-31777\168aef4b-6261-4472-bfc1-338a53529380.jpg" />, system (1.1) becomes</p><disp-formula id="scirp.19865-formula38074"><label>(1.2)</label><graphic position="anchor" xlink:href="15-31777\73f4cb0b-844c-41b6-b55a-d55bf67946b6.jpg"  xlink:type="simple"/></disp-formula><p>which is Hamiltonian system. Now suppose that the Hamiltonian system (1.2) has a nilpotent center at the origin, namely the function H satisfies the following conditions:</p><p>(H1) <img src="15-31777\836b9a29-619f-4561-a238-def1fc2806a9.jpg" />is a <img src="15-31777\a96100c8-7b1d-4c7e-ab3a-b61f2a2205f1.jpg" /> function, satisfying</p><p><img src="15-31777\15ac29c0-2380-427d-843a-5d29599bf5ab.jpg" />;</p><p>(H2)<img src="15-31777\2c089ba3-9006-47c7-93b5-68985c6dfb20.jpg" />, the equation <img src="15-31777\44032259-85b6-4e41-9c7e-25763e331166.jpg" /> defines a closed curve L<sub>h</sub> surrounding the origin and L<sub>h</sub> approaches the origin as h goes to zero;</p><p>(H3)<img src="15-31777\e0aba154-3600-4f9e-b405-8c37c9390a67.jpg" /> ,<img src="15-31777\180602d1-56a9-4bb4-8950-983583d7a631.jpg" />.</p><p>It follows that the expansion of H at the origin has the form</p><p><img src="15-31777\74119c01-4920-48d7-894d-7dbedacce3cc.jpg" /></p><p>Assume that the equation <img src="15-31777\8eb0cfff-0ee6-47c4-8e6c-a8ab4127a352.jpg" /> intersects the positive x-axis at<img src="15-31777\ba9ae5d5-7f2f-469e-a4cf-c2dd260ae69a.jpg" />. Let <img src="15-31777\5eecfbc8-ec51-4191-9e41-ebff8a460b26.jpg" /> denote the first intersection point of the positive orbit of (1.1) starting at <img src="15-31777\e4883948-d636-4e61-a2b3-7862c5509480.jpg" /> with the positive x-axis. Then, we have</p><disp-formula id="scirp.19865-formula38075"><label>(1.3)</label><graphic position="anchor" xlink:href="15-31777\fff30607-5aac-4f98-bcb3-a0272503ddce.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.19865-formula38076"><label>(1.4)</label><graphic position="anchor" xlink:href="15-31777\c138f7dd-913a-42a5-97da-82a4f38d9578.jpg"  xlink:type="simple"/></disp-formula><p>The Abelian integral M above is called the first order Melnikov function of system (1.1). From Han [<xref ref-type="bibr" rid="scirp.19865-ref8">8</xref>], we have a general theorem as follows.</p><p>Theorem 1.1. Suppose that the origin is nilpotent singular point <img src="15-31777\ba80a2fe-6178-49f9-b2be-fa6353ddfc13.jpg" /> and that <img src="15-31777\7a551fde-c244-4e9e-ac64-2de8e86d15a8.jpg" /> approaches the origin as h goes to zero. If there exist an integer <img src="15-31777\8e1e298a-f5df-4316-b4fd-aeb952dc6044.jpg" /> and <img src="15-31777\e03fefea-0f80-47e2-b723-87c30660d8a9.jpg" /> such that</p><p><img src="15-31777\146242d5-8884-43f2-bf5a-b6d89dc8fbb2.jpg" /></p><p>and</p><p><img src="15-31777\3042bebf-7f67-4b19-bf0f-5338c6dfc90c.jpg" /></p><p>then we have 1) <img src="15-31777\cb12b755-1149-4af3-a3bc-bee780c94048.jpg" />has at most k zeros near <img src="15-31777\d1e7255b-615d-418e-9dbd-3db6d095189d.jpg" /> for <img src="15-31777\e6574012-d8d5-4487-8733-2b5e97d7029b.jpg" /> and all <img src="15-31777\a0fab2eb-ae7a-4a09-939f-dbc5bc6b9121.jpg" /> near<img src="15-31777\668744c9-c3ea-401f-9a16-42adaecaf5b3.jpg" />, and k zeros can appear for some <img src="15-31777\06ed629b-958f-4e15-afdd-b8592ebed37a.jpg" /> near<img src="15-31777\57b22341-e2c6-4b02-8290-c121c1a0e206.jpg" />.</p><p>2) System (1.1) has at least k limit cycles near the origin for some <img src="15-31777\466e9334-c51d-498e-a934-0cf07d292810.jpg" /> near<img src="15-31777\df3027bf-6e1b-4ba2-86af-c0a4c6de7fbd.jpg" />.</p></sec><sec id="s2"><title>2. Main Results and Proof</title><p>Consider the following near-Hamiltonian system:</p><disp-formula id="scirp.19865-formula38077"><label>(2.1)</label><graphic position="anchor" xlink:href="15-31777\11822ee4-d43a-48f6-bba6-0778d65136d9.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="15-31777\330bb26e-7650-42f0-9af6-8b1cfbe45981.jpg" /> and p and q are cubic polynomials. We can write</p><disp-formula id="scirp.19865-formula38078"><label>(2.2)</label><graphic position="anchor" xlink:href="15-31777\c439834d-9122-4cb0-91d2-381bebb3846b.jpg"  xlink:type="simple"/></disp-formula><p>Then unperturbed system <img src="15-31777\8670ea65-732c-44dc-9403-a6ac5039f9ba.jpg" /> is a Hamiltonian system with Hamiltonian</p><disp-formula id="scirp.19865-formula38079"><label>(2.3)</label><graphic position="anchor" xlink:href="15-31777\be086551-0f87-4373-9f54-53d82d75d3e5.jpg"  xlink:type="simple"/></disp-formula><p>system <img src="15-31777\4e7971bc-02d2-4e29-ad46-4476061d04a9.jpg" /> has a nilpotent center at the origin. Let <img src="15-31777\7db9d2d2-d572-475e-ae80-e96dc2d9eb34.jpg" /> be the closed curve defined by<img src="15-31777\42559581-0e3a-4a29-84c1-a392cb9890b8.jpg" />. Then it can be presented as</p><disp-formula id="scirp.19865-formula38080"><label>(2.4)</label><graphic position="anchor" xlink:href="15-31777\51aa51ab-fbd4-4dd5-9db4-38247d74a97a.jpg"  xlink:type="simple"/></disp-formula><p>Assume that the positive solution of the above equation in y is</p><disp-formula id="scirp.19865-formula38081"><label>(2.5)</label><graphic position="anchor" xlink:href="15-31777\03f0ad31-5267-4f70-be7e-f8ac797d747c.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="15-31777\a99a9a23-ee5f-41ca-924b-f2e753403d08.jpg" /> and<img src="15-31777\a7c9b81e-63f2-44b4-b86a-1356740088b5.jpg" />. Then by (2.4) and (2.5) we obtain</p><p><img src="15-31777\beaf9980-9b01-4ae6-a40b-8d25b851f110.jpg" /></p><p><img src="15-31777\d1cabc66-0bce-4878-a159-61393df2a5e5.jpg" /></p><p><img src="15-31777\99ae392a-acf1-4e7b-8162-9ec580b08246.jpg" /></p><p><img src="15-31777\cf48317a-c978-4e27-8bba-a1bfe1121baf.jpg" /></p><p><img src="15-31777\206ab8ff-d8f3-48f3-9dac-52313c1be7ba.jpg" /></p><p>By [<xref ref-type="bibr" rid="scirp.19865-ref8">8</xref>] the negative solution of (2.4) in y satisfies<img src="15-31777\0d7e0caf-d918-48c2-8ec1-91e45dcaf1a9.jpg" />. Thus, two solutions of (2.4) are</p><disp-formula id="scirp.19865-formula38082"><label>(2.6)</label><graphic position="anchor" xlink:href="15-31777\03840109-9131-4e0d-94a9-c8e8204f03e7.jpg"  xlink:type="simple"/></disp-formula><p>On the other hand, the intersection points of L<sub>h</sub> and xaxis have the x-coordinates <img src="15-31777\d65e2543-1cfd-4e23-9b49-9d83bd348178.jpg" /> and<img src="15-31777\b1beab74-57c0-4cd6-b81b-fc8bfa487407.jpg" />. Then by (2.2) we can write</p><disp-formula id="scirp.19865-formula38083"><label>(2.7)</label><graphic position="anchor" xlink:href="15-31777\0f6ffae7-309c-4696-8ab9-939602d2d1ee.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.19865-formula38084"><label>(2.8)</label><graphic position="anchor" xlink:href="15-31777\542f5140-ea79-4aa7-87c7-cea15787d2a2.jpg"  xlink:type="simple"/></disp-formula><p>Here,</p><p><img src="15-31777\12692aa3-04f6-43e8-8ba7-256c04f3c0eb.jpg" /></p><p>Introduce</p><disp-formula id="scirp.19865-formula38085"><label>(2.9)</label><graphic position="anchor" xlink:href="15-31777\e6ba5801-37bd-4344-a3b4-5017130a1674.jpg"  xlink:type="simple"/></disp-formula><p>Then, similar to the method of Han [<xref ref-type="bibr" rid="scirp.19865-ref8">8</xref>] we have</p><disp-formula id="scirp.19865-formula38086"><label>(2.10)</label><graphic position="anchor" xlink:href="15-31777\281a236f-f89b-460b-85bc-ab3f83d80533.jpg"  xlink:type="simple"/></disp-formula><p>Therefore, in turn by (2.6)-(2.10) we have</p><p><img src="15-31777\61782ce8-4929-4238-bb25-ff760228e72b.jpg" /></p><p><img src="15-31777\3848b566-344b-4ee0-af60-c13bc7f628cc.jpg" /></p><p><img src="15-31777\eafac5d7-2fb5-433e-8ffc-d3fe035796c7.jpg" /></p><p>Noting that<img src="15-31777\2dfe1300-275a-4a6f-8c6c-e6913e5c61d9.jpg" />, then similarly we have</p><p><img src="15-31777\c7bd9364-e33d-4cbc-9844-8d496f857dc5.jpg" /></p><p><img src="15-31777\298fb5f4-27ad-4b48-a466-e8df348a592d.jpg" /></p><p>In the same way, using<img src="15-31777\dd2c32b6-a52e-4ce8-bd52-d50288a5f86f.jpg" />, we have</p><p><img src="15-31777\c40bb88a-3921-4351-ad61-4e3b175e10eb.jpg" /></p><p>Hence, we have</p><p><img src="15-31777\ede644e4-3365-4c62-b298-b00bc91b5acb.jpg" /></p><p>where <img src="15-31777\336d7430-92a0-4d7a-bbe6-581867b3019f.jpg" /> And</p><p><img src="15-31777\77183b3b-d2af-4a4b-a37d-05daded9fc85.jpg" /></p><p><img src="15-31777\92ae910c-459c-493c-ba9f-27b83a4bb82e.jpg" /></p><p><img src="15-31777\031af091-e553-410e-8839-20dddc0d6c9d.jpg" /></p><p><img src="15-31777\83d9c693-4cd2-482e-a613-4c775c53332b.jpg" /></p><p><img src="15-31777\807878ec-b4f1-49b8-a4ea-e7a1700a8640.jpg" /></p><p><img src="15-31777\da4fda6b-063c-43bf-9633-3f3dd9fd14fc.jpg" /></p><p>Now it is direct that</p><p><img src="15-31777\513c5485-ad0c-40b0-a1d4-3bde73588639.jpg" /></p><p><img src="15-31777\7367e600-b338-4e25-a032-078d7361be26.jpg" /></p><p>Here, if let<img src="15-31777\bfdbef28-6838-428d-897b-fae1be853958.jpg" />, <img src="15-31777\9e86d685-c480-4802-bdb8-21aeed98f5c9.jpg" />, <img src="15-31777\b3e18131-e091-4fda-855d-a4ec7d25a732.jpg" />, then for some cubic system (2.1) we can obtain the above determinant is not zero, then the function M can have 5 simple zeros in h &gt; 0 near h = 0 for some <img src="15-31777\3d038903-2c1d-4620-8dc6-efca3f0db4c3.jpg" /> near<img src="15-31777\7e4257b2-0115-476e-b82c-0512b975f992.jpg" />. For example, let</p><p><img src="15-31777\733b2680-795a-4f2d-a651-6d8d95594b09.jpg" />, <img src="15-31777\826ac176-da3b-4ae0-8641-d4a0bad81d35.jpg" />, <img src="15-31777\5d86fa98-e259-4b13-a9cf-6d0927b40eb4.jpg" />, we obtain from the above formula</p><p><img src="15-31777\b4b92a95-72e8-48a6-b18b-d2015c636404.jpg" /></p><p>Here,</p><p><img src="15-31777\30c1dbe8-90c4-425d-8782-6ba89d0e7fb4.jpg" /><img src="15-31777\7b2ed31e-364c-4a4f-8002-e3b9772ac3ce.jpg" /><img src="15-31777\03aa4fe4-d9c5-4a90-b439-96597fcf79d5.jpg" /><img src="15-31777\d3572220-fb34-4bf4-b1e6-8d9856e54927.jpg" /><img src="15-31777\aca7f428-4940-4b6e-9f34-9ad2d6398abf.jpg" /><img src="15-31777\263a5e8d-de01-4bd4-a318-1fdca70fc4b4.jpg" /><img src="15-31777\b836bd7d-384c-4568-b206-f896f70d0fcd.jpg" /><img src="15-31777\015c5f2b-3ccf-4cb9-a1da-8a9d654c5b34.jpg" /><img src="15-31777\105bce3e-bf7c-4a8a-87d5-d0e33e26fb69.jpg" /><img src="15-31777\39e5ff0e-ada0-4d88-b3a4-4bac9ea90f7d.jpg" /><img src="15-31777\c9e21f58-f50e-4e56-bf13-79ef96d249de.jpg" /><img src="15-31777\7db21a32-57e1-4ecb-a389-d81656e5c0a4.jpg" /></p><p>then we can obtain</p><disp-formula id="scirp.19865-formula38087"><label>(2.11)</label><graphic position="anchor" xlink:href="15-31777\290f6a23-a449-4575-abb0-1a9e483b71cb.jpg"  xlink:type="simple"/></disp-formula><p>By Theorem 1.1 we have:</p><p>Theorem 2.1. The function <img src="15-31777\239132d6-4ace-445a-959f-bcfa333c36cf.jpg" /> has at most 5 zeros in <img src="15-31777\b0846c10-370d-4ab3-9466-d1fec491cfe0.jpg" /> near<img src="15-31777\cb5c4076-32c7-4863-8779-2a68a8cbf703.jpg" />, and for <img src="15-31777\783aa142-d967-4a60-95d6-1d09d1e431a8.jpg" /> small the cubic system (2.1) can have 5 limit cycles near the origin.</p></sec><sec id="s3"><title>REFERENCES</title></sec><sec id="s4"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.19865-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">A. Andronov, E. Leontovich, I. Gordon and A. Maier, “Theory of Bifurcations of Dynamical Systems on a Plane,” Israel Program for Scientific Translations, Jerusalem, 1971.</mixed-citation></ref><ref id="scirp.19865-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">M. Han, “On Hopf Cyclicity of Planar Systems,” Journal of Mathematical Analysis and Applications, Vol. 245, No. 2, 2000, pp. 404-422. doi:10.1006/jmaa.2000.6758</mixed-citation></ref><ref id="scirp.19865-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Y. A. Kuznetsov, “Elements of Applied Bifurcation Theory,” Springer-Verlag, New York, 1995.</mixed-citation></ref><ref id="scirp.19865-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">T. Carmon, R. Uzdin, C. Pigier, Z. Musslimani, M. Segev and A. Nepomnyashchy, “Rotating Propeller Solitons,” Physical Review Letters, Vol. 87, No. 14, 2001, p. 143901. 
doi:10.1103/PhysRevLett.87.143901</mixed-citation></ref><ref id="scirp.19865-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">J. Guckenheimer and P. Holmes, “Non-Linear Oscillations, Dynamical Systems and Bifurcation of Vector Fields,” Springer-Verlag, New York, 1983.</mixed-citation></ref><ref id="scirp.19865-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">B. D. Hassard, N. D. Kazarinoff and Y. H. Wan, “Theory and Applications of Hopf Bifurcation,” Cambridge University Press, Cambridge, 1981.</mixed-citation></ref><ref id="scirp.19865-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">S. Wiggins, “Global Bifurcations and Chaos: Analytical Methods,” Springer-Verlag, New York, 1988.</mixed-citation></ref><ref id="scirp.19865-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">M. Han J. Jiang and H. Zhu, “Limit Cycle Bifurcations in Near-Hamiltonian Systems by Perturbing a Nilpotent Center,” International Journal of Bifurcation and Chaos, Vol. 18, No. 10, 2008, pp. 3013-3027.  
doi:10.1142/S0218127408022226</mixed-citation></ref></ref-list></back></article>