閆曉芳, 尚華輝, 楊紀(jì)華

(1. 永城職業(yè)學(xué)院 基礎(chǔ)部, 河南 永城 476600; 2. 寧夏師范學(xué)院 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 寧夏 固原 756000)

1 引言和主要結(jié)果

對(duì)于一個(gè)擾動(dòng)可積非Hamilton系統(tǒng)

(1)

其中,0<|ε|?1,Hy(x,y)μ(x,y)和Hx(x,y)μ(x,y)是關(guān)于x、y的m次多項(xiàng)式,f(x,y)和g(x,y)是關(guān)于x、y的n次多項(xiàng)式,H(x,y)稱(chēng)為系統(tǒng)(1)(ε=0)具有積分因子1/μ(x,y)的首次積分.假設(shè)未擾動(dòng)系統(tǒng)(1)(ε=0)至少有一個(gè)中心,確定系統(tǒng)(1)極限環(huán)的個(gè)數(shù)和分布稱(chēng)為Hilbert16問(wèn)題[1],通過(guò)擾動(dòng)一個(gè)系統(tǒng)的中心得到極限環(huán)是一種經(jīng)典的方式.一般來(lái)說(shuō),研究從一個(gè)系統(tǒng)的中心周期環(huán)域分支出極限環(huán)的個(gè)數(shù)問(wèn)題,共有4種方法:第一種是Poincare回歸映射法[2-3],第二種是Poincare-Pontryagin-Melnikov法(或稱(chēng)為Abelian積分法)[4-7],第三種是可逆積分因子法[8-9],第四種是平均法[10-13].在平面上,第二種和第四種方法等價(jià),見(jiàn)文獻(xiàn)[11].應(yīng)用平均法研究系統(tǒng)(1)的極限環(huán)分支問(wèn)題,由于需要做極坐標(biāo)變換,所以目前文獻(xiàn)要求H(x,y)=x2+y2,例如文獻(xiàn)[10]研究了μ(x,y)=x+1的情形,文獻(xiàn)[11]研究了μ(x,y)是平面二次圓錐曲線的情形,文獻(xiàn)[12]研究了μ(x,y)=(x+a)(y+b)的情形.

本文研究如下擾動(dòng)可積非Hamilton系統(tǒng)

(2)

其中

2 平均法和平均函數(shù)

首先介紹微分方程的平均法,更詳細(xì)的介紹參考文獻(xiàn)[6].

定理2考慮一個(gè)微分方程的初值問(wèn)題

(3)

其中,F0(x,t)和G0(x,t,ε)是關(guān)于t的T-周期函數(shù),x,x0∈U,T是不依賴(lài)于ε的常數(shù),U是R中的開(kāi)區(qū)間.定義平均函數(shù)

再考慮平均方程的初值問(wèn)題

(4)

注1由定理2可知,如果方程(3)滿(mǎn)足定理中的條件,則平均函數(shù)f0(y)的一個(gè)簡(jiǎn)單零點(diǎn)對(duì)應(yīng)方程(3)的一個(gè)極限環(huán).所以計(jì)算出平均函數(shù)f0(y)至關(guān)重要.

所以

ε2G0(r,θ,ε).

(5)

容易驗(yàn)證,方程(5)滿(mǎn)足定理2中的條件.根據(jù)定理2,與方程(5)相對(duì)應(yīng)的平均函數(shù)為

下面化簡(jiǎn)平均函數(shù)f0(r),計(jì)算可得

(7)

其中λi,j=ai-1,j+bi,j-1,這里假設(shè)λ0,0=a-1,j=bi,-1=0.因?yàn)閍i,j和bi,j是任意的,所以λi,j也是任意的.定義

(8)

引理1下列關(guān)系式成立：

(i)Ii,2j+1(r)=0,

(ii)r2iI2i,0(r)=(-1/a)iI0,0(r)-

(ii)由

(ar2cos2θ)i-(-1)i=(ar2cos2θ+1)×

得

把上式兩端同時(shí)關(guān)于θ從0到π積分即可得(ii)成立.引理1證畢.

(9)

其中

因?yàn)棣薸,j是任意的,所以Ai,j也是任意的.又λ0,0=0,所以A0,0=0.再由(9)式和引理1中(ii)可得

其中

i=1,2,…,n-1.

(10)

注意到,當(dāng)k是奇數(shù)時(shí),M(k)=0,即β2k+1=0,所以

F(r)=α0[I0,0(r)+α/α0]+

(11)

注2由(9)式可知

所以

3 定理1的證明

首先把函數(shù)F(r)表示成若干個(gè)線性無(wú)關(guān)函數(shù)的線性組合.

I0,0(r)+α/α0,r2I0,0(r),r4I0,0(r),…,

(12)

證明由(11)式知F(r)可以由(12)式中函數(shù)線性表出.所以只需證明它們是線性無(wú)關(guān)的即可.對(duì)于它們的任意線性組合

δ0(I0,0(r)+α/α0)+(δ1r2+δ2r4+…+

(δ0+δ1r2+δ2r4+…+

(13)

δ0+δ1(-1/a)+δ2a-2+…+

所以(13)式變?yōu)?/p>

為了證明本文中定理1,還需要下面引理,見(jiàn)文獻(xiàn)[14].

引理3考慮n+1個(gè)線性無(wú)關(guān)的解析函數(shù)hi:U→R,i=1,2,…,n+1,其中U?R是開(kāi)區(qū)間.假設(shè)存在j∈{1,…,n+1}使得hj不變號(hào),則存在n+1個(gè)常數(shù)ci,i=1,…,n+1,使得

在U中至少有n個(gè)孤立零點(diǎn).

I0,0(r)+α/α0,r2I0,0(r),r4I0,0(r),…,

4 結(jié)論和展望

[1]ARNOLDVI.Tenproblemsin:theoryofsingularitiesanditsapplications[J].AdvSovietMath,1990,1:1-8.

[2]CHICONEC,JACOBSM.Bifurcationoflimitcyclesfromquadraticisochrones[J].JDiffEqns,1991,91(2):268-326.

[3]LICZ,LLIBREJ,ZHANGZF.Weakfocus,limitcyclesandbifurcationsforboundedquadraticsystems[J].JDiffEqns,1995,115(1):193-223.

[4]ATABAIGIA,NYAMORADIN,ZANGENEHH.Thenumberoflimitcyclesofaquinticpolynomialsystemwithacenter[J].NonlinearAnalysis,2009,71(7):3008-3017.

[5]YAOHY,HANMA.Thenumberoflimitcyclesofaclassofpolynomialdifferentialsystems[J].NonlinearAnalysis,2012,75(1):341-357.

[6]GASULLA,LICZ,TORREGROSAJ.Limitcyclesappearingfromtheperturbationofasystemwithamultiplelineofcriticalpoints[J].NonlinearAnalysis,2012,75(1):278-285.

[7]XIONGYQ.Thenumberoflimitcyclesinperturbationsofpolynomialsystemswithmultiplecirclesofcriticalpoints[J].JMathAnalAppl,2016,440(1):220-239.

[8]GIACOMINIH,LLIBREJ,VIANOM.Ontheshapeoflimitcyclesthatbifurcatefromnon-Hamiltoniancenters[J].NonlinearAnalysis,2001,43(7):837-59.

[9]VIANOM,LLIBREJ,GIACOMINIH.ArbitraryorderbifurcationsforperturbedHamiltonianplanarsystemsviathereciprocalofanintegratingfactor[J].NonlinearAnalysis,2002,48(1):117-36.

[10]LLIBREJ,RIOJ,RODRIGUEZJ.Averaginganalysisofaperturbatedquadraticcenter[J].NonlinearAnalysis,2001,46(1):45-51.

[11]GINEJ,LLIBREJ.Limitcyclesofcubicpolynomialvectorfieldsviatheaveragingtheory[J].NonlinearAnalysis,2007,66(8):1707-1721.

[12]BUICAA,LLIBREJ.Limitcyclesofaperturbedcubicpolynomialdifferentialcenter[J].Chaos,Solitons&Fractals,2007,32(3):1059-1069.

[13]COLLB,LLIBREJ,PROHENSR.Limitcyclesbifurcatingfromaperturbedquarticcenter[J].Chaos,Solitons&Fractals,2011,44(4):317-334.

[14]COLLB,GASULLA,PROHENSR.Bifurcationoflimitcyclesfromtwofamiliesofcenters[J].DynContinDiscreteImpulsSyst:MathAnal,2005,12(2):275-287.

[15]VERHULSTF.NonlinearDifferentialEquationsandDynamicalSystems[M].Berlin:Springer-Verlag,1996.