曹 明
(陜西學(xué)前師范學(xué)院 數(shù)學(xué)系, 陜西 西安 710100)
?
一類(n+1)次多項(xiàng)式系統(tǒng)極限環(huán)的存在性與唯一性
曹明
(陜西學(xué)前師范學(xué)院 數(shù)學(xué)系, 陜西 西安710100)
摘要:基于Liapunov級(jí)數(shù)理論對(duì)一類(n+1)次多項(xiàng)式系統(tǒng)在原點(diǎn)的焦點(diǎn)量和極限環(huán)問題進(jìn)行了研究,得到了原點(diǎn)是該系統(tǒng)的焦點(diǎn)或中心的一個(gè)充分條件和該系統(tǒng)依賴于參數(shù)δ的Hopf分支問題,且分析了當(dāng)參數(shù)滿足δma=0時(shí)系統(tǒng)極限環(huán)的存在性與唯一性.
關(guān)鍵詞:多項(xiàng)式系統(tǒng); Liapunov級(jí)數(shù); 分支; 極限環(huán)
0引言
在文獻(xiàn)[1-3]中,作者提出了相伴系統(tǒng)的概念,并詳細(xì)研究了二次系統(tǒng)的相伴系統(tǒng)的極限環(huán)問題.基于文獻(xiàn)[4,5]的啟發(fā),本文討論的是另一類(n+1)次多項(xiàng)式系統(tǒng)
(1)
(方程均為實(shí)系數(shù)多項(xiàng)式,且bn≠0).
系統(tǒng)(1)是在三次系統(tǒng)(2)
(2)
基礎(chǔ)上添加一些實(shí)或虛的不變直線y=yi得到的,其中yi是方程φ(y)=1+b1y+b2y2+…+bnyn=0(bn≠0)的根,所以稱系統(tǒng)(1)和(2)為一對(duì)相伴系統(tǒng).令方程φ(y)=0的最大負(fù)實(shí)根為A(若A不存在,則記A=-∞),最小正實(shí)根為B(若B不存在,則記B=-∞).
(3)
取新的變換x=x,y=H-1(u),系統(tǒng)(3)進(jìn)一步可化為廣義Lienard系統(tǒng)(4)
(4)
其中H(A)
1中心和焦點(diǎn)問題
為了詳細(xì)地研究系統(tǒng)(1)的拓?fù)浣Y(jié)構(gòu),首先考察系統(tǒng)(1)對(duì)應(yīng)線性系統(tǒng)
(5)
易知系統(tǒng)(5)有唯一的奇點(diǎn)O(0,0).若給系統(tǒng)(5)附加非線性項(xiàng)lx2+mx2y+ax3和b1xy+b2xy2+…+bnxyn,(bn≠0),就得到系統(tǒng)(1).應(yīng)用文獻(xiàn)[5]的定理4.7,我們有如下引理:
引理1對(duì)于自治微分系統(tǒng)(1),我們有(Ⅰ)若|δ|≥2,則O(0,0)為結(jié)點(diǎn);(Ⅱ)若0<|δ|<2,則O(0,0)為(粗)焦點(diǎn).在系統(tǒng)(1)中,當(dāng)|δ|≥2時(shí)O(0,0)為結(jié)點(diǎn),當(dāng)0<|δ|<2時(shí)屬于中心焦點(diǎn)判定問題.
引理2對(duì)于系統(tǒng)(1),O(0,0)的焦點(diǎn)量是:ω0=δ;若ω0=0,則ω1=a;若ω0=ω1=0,則ω2=ma;若ω0=ω1=ω2=0,則O(0,0)是中心.若ωi>0(<0),ω0=ω1=ω2=…=ωi-1=0,那么O(0,0)是不穩(wěn)定(穩(wěn)定)細(xì)焦點(diǎn)(其中ωi是O的第i階焦點(diǎn)量).
證明:為了獲得O的焦點(diǎn)量,我們采用文獻(xiàn)[2]和[6]的方法和記號(hào).在系統(tǒng)(4)中令δ=0,F(x),g(x)按冪級(jí)數(shù)展開如下列形式:
F(x)=lx2+ax3+lmx4+amx5+lm2x6+
ax7+lm3x8+……
f(x)=F′(x)=2lx1+3ax2+4lmx3+5amx4+…
≡a1x+a2x2+a3x3+…+
g(x)=x+mx3+m2x5+m3x7+…+mnx2n+1+…
≡C0x+C1x2+C2x3+…+
g3=2β1,g5=2(-4c1β2+β3),
通過計(jì)算可得,g3=2a,g5=2am.由文獻(xiàn)[6]定理5知,如果ω0=δ,g3=2β1=2a>0(<0)則O(0,0)是一階不穩(wěn)定(穩(wěn)定)的細(xì)焦點(diǎn);δ=g3=a=0,g5=0,則O(0,0)是中心.記O的第i階焦點(diǎn)量為ωi,由于系統(tǒng)(1)化為系統(tǒng)(4)時(shí)O的穩(wěn)定性改變了,因此有ω0=δ;若ω0=0,則ω1=a;若ω0=ω1=0,則ω2=ma;若ω0=ω1=ω2=0,則O(0,0)是中心.
如果ωi>0(<0),ω0=ω1=ω2=…ωi-1=0,
那么O(0,0)是不穩(wěn)定(穩(wěn)定)細(xì)焦點(diǎn)(其中ωi是O的第i階焦點(diǎn)量).
若ω0=ω1=ω2=0,
即δ=a=ma=0時(shí),系統(tǒng)(1)可化為
(6)
因?yàn)閄(-x,y)=X(x,y),Y(-x,y)=-Y(x,y),系統(tǒng)(6)的向量場(chǎng)是關(guān)于y軸對(duì)稱的,由平面系統(tǒng)的對(duì)稱原理可得O(0,0)是中心.
2系統(tǒng)的極限環(huán)
不失一般性,我們不妨假設(shè)l≥0(否則可通過變換(x,y,t)→(x,-y,-t)來改變符號(hào)),a≥0(否則可通過變換(x,y,t)→(-x,y,-t)來改變符號(hào)).同時(shí)在系統(tǒng)(4)中我們?nèi)砸?x,y,t)記(x,μ,τ).
考慮l=0和l≠0兩種情形下系統(tǒng)(1)的極限環(huán).
定理1當(dāng)l=0時(shí),考慮一類自治微分系統(tǒng)(1)有:
(Ⅰ)若a>0,則原點(diǎn)O為一階穩(wěn)定細(xì)焦點(diǎn),則當(dāng)0<δ?1,在原點(diǎn)O的外圍至少存在(1)的一個(gè)不穩(wěn)定極限環(huán).
(Ⅱ)a<0,則原點(diǎn)O為一階不穩(wěn)定細(xì)焦點(diǎn),則當(dāng)-1<δ?0,在原點(diǎn)O的外圍至少存在(1)的一個(gè)穩(wěn)定極限環(huán).
證明: 首先基于Liapunov形式級(jí)數(shù)法理論,得到了O是系統(tǒng)(1)的焦點(diǎn)或中心的一個(gè)充分條件.
假設(shè)形式級(jí)數(shù)
其中Fk(x,y)是x,y的k次其次式,k=3,4,5…,于是
-2mcos3θsinθ-2acos4θ+2(b12-b2)cosθsin3θ
自然地,可以得到如下結(jié)論:
引入極坐標(biāo)變換x=rcosθ,y=rsinθ后,經(jīng)計(jì)算得結(jié)果
不失一般性,令α=-2b13+4b1b2-2b3,
引入極坐標(biāo)變換x=rcosθ,y=rsinθ后,經(jīng)計(jì)算得結(jié)果
F6(x,y)=
基于Liapunov形式級(jí)數(shù)法理論自然得到當(dāng)a=0時(shí),原點(diǎn)O是系統(tǒng)(1)的中心.
注:當(dāng)δ=0且l=0時(shí),一類多項(xiàng)式系統(tǒng)(1)在有限遠(yuǎn)奇點(diǎn)O(0,0)處的拓?fù)浣Y(jié)構(gòu)完全取決于參數(shù)α的符號(hào).
繼續(xù)證明定理1的后半部分,只證明情形(Ⅰ),情形(Ⅱ)類似可得.
顯然O(0,0)是系統(tǒng)(1)|-2<δ<0穩(wěn)定(粗)焦點(diǎn),當(dāng)l=0,a<0,且參數(shù)δ從-1增大到0時(shí),O(0,0)是系統(tǒng)(1)|δ=l=0的不穩(wěn)定細(xì)焦點(diǎn).根據(jù)Hopf分支理論,在原點(diǎn)O(0,0)的外圍至少存在(1)的一個(gè)不穩(wěn)定極限環(huán).
定理2如果δma=0,則系統(tǒng)(1)或(4)至多有一個(gè)極限環(huán).具體內(nèi)容是:
(1)如果a=0,則當(dāng)系統(tǒng)(1)或(4)滿足條件δm≠0時(shí)沒有極限環(huán).
(2)如果m=0,則當(dāng)δ≥0時(shí),系統(tǒng)(1)或系統(tǒng)(4)沒有極限環(huán);當(dāng)δ<0時(shí),系統(tǒng)至多有一個(gè)極限環(huán),若存在,必穩(wěn)定.
(3)如果δ=0時(shí),則當(dāng)系統(tǒng)滿足ml≠0時(shí)系統(tǒng)(1)或(4)沒有極限環(huán).
此定理可分為三個(gè)部分來證明:
第一部分,可視為命題1:如果a=0,則當(dāng)系統(tǒng)(1)或(4)滿足條件δm≠0時(shí)沒有極限環(huán).
證明:當(dāng)a=0時(shí),則在系統(tǒng)(4)中
我們定義曲線L和H如下
(u,v)∈D={(u,v:v<0
mu2<1)}
由F(u)=F(v)得:
(u-v)[δ+δmuv+l(u+v)]=0
即 [δ+δmuv+l(u+v)]=0
(7)
(u-v)[-δ+mδ(u+v)2+2lm(u+v)+
m2δu2v2]=0
即
[-δ+mδ(u+v)2+2lmuv(u+v)+
m2δu2v2]=0
(8)
假設(shè)曲線L和H在區(qū)域D內(nèi)相交,交點(diǎn)(u,v),即聯(lián)立(7)、(8)可得:
δ(1+2m2u2v2-mu2-mv2)=0
(9)
當(dāng)δ≠0且m≠0時(shí),(9)等價(jià)于下面兩式
m(u2+v2)=1+2m2u2v2
(10)
如果L和H相交,則由(10)可知mu2和mv2一定滿足一元二次方程
z2-(1+2m2u2v2)z+m2u2v2=0
(11)
可是由(10)、(11)知:
1+m2u2v2
當(dāng)δ≠0且m≠0時(shí),
顯然mu2v2≠1+m2u2v2,故L和H在區(qū)域內(nèi)D沒有交點(diǎn),從而由文獻(xiàn)[8]知系統(tǒng)(1)α=0或(4)α=0沒有極限環(huán).
第二部分可視為命題2:如果m=0,則當(dāng)δ≥0時(shí),系統(tǒng)(1)或(4)沒有極限環(huán);當(dāng)δ<0時(shí),系統(tǒng)至多有一個(gè)極限環(huán),若存在,必穩(wěn)定.
證明:記系統(tǒng)(1)的右邊分別P(x,y)為Q(x,y)和則有
所以當(dāng)A
如果δ<0,則在系統(tǒng)(4)m=0中
F(x)=δx+lx2+ax3
f(x)=F′(x)=δ+2lx+3ax2
g(x)=x
易于驗(yàn)證滿足下列條件
(a)φ(0)=0,φ′(y)>0;
(b)xg(x)>0,x≠0,G(±∞)=+∞;
由文獻(xiàn)[8]的定理5.9可知系統(tǒng)(4)m=0至多存在一個(gè)極限環(huán);若存在,必穩(wěn)定.在下面討論系統(tǒng)(1)或(4)的極限環(huán)時(shí),不妨令a>0,m≠0.
如果δ=0則系統(tǒng)(4)中
同理定義曲線L和H如下
(u,v)∈D={(u,v:v<0
mv2<1,mu2<1)}
由F(u)=F(v),得:
(u-v)[l(u+v)+a(u2+uv+v2)-
amu2v2]=0
[l(u+v)+a(u2+uv+v2)-amu2v2]=0
(12)
(u-v)[3a-am(u2+uv+v2)+
2lm(u+v)+3amuv+am2u2v2]=0,
即
[3a-am(u2+uv+v2)+2lm(u+v)+
3amuv+am2u2v2]=0
(13)
假設(shè)曲線L和H在區(qū)域D內(nèi)相交,交點(diǎn)(u,v),即聯(lián)立(12)、(13)可得:
a+ml(u+v)+amuv=0
(14)
第三部分可視為命題3:如果δ=0時(shí),則當(dāng)系統(tǒng)滿足ml≠0時(shí),系統(tǒng)(1)或(4)沒有極限環(huán).
證明:因?yàn)楫?dāng)δ=0且m≠0時(shí),(14)等價(jià)于下面兩式
(15)
如果L和H相交,則由(10)可知mu和mv一定滿足一元二次方程
(16)
可是由(15)、(16)知:
muv=-1
(17)
綜上三個(gè)命題的證明,定理2得證.
參考文獻(xiàn)
[1] 謝向東,張劍峰.平面多項(xiàng)式系統(tǒng)及其相伴系統(tǒng)[J].數(shù)學(xué)研究,2004,37(2):161-166.
[2] Xie Xiangdong,Chen Fengde.The uniqueness of limit cycle and the structure of critical point at infinity for a class of cubic system[J].Ann of Diff Eqs,2005,21(3):474-479.
[3] 謝向東,陳鳳德.一類具有兩虛不變直線的三次系統(tǒng)的極限環(huán)與分支[J].數(shù)學(xué)物理學(xué)報(bào),2005,25(4):538-545.
[4] 金山,魯世平.一類多項(xiàng)式系統(tǒng)極限環(huán)的唯一性與分支[J].數(shù)學(xué)物理學(xué)報(bào),2011,31(6):1 669-1 673.
[5] Du Jia,Xao Jian.Existence of limit cycles for a class of (n+1)-th polynomial systems[J].Chinese Journal of Engineering Mathematics,2014,31(2):274-285.
[6] 杜佳,肖箭.一類多項(xiàng)式系統(tǒng)中心焦點(diǎn)及Hopf分支問題[J].合肥師范學(xué)院學(xué)報(bào),2012,30(6):14-16.
[7] 蔡健林,馬暉.廣義方程的奇點(diǎn)的中心焦點(diǎn)判定問題[J].浙江大學(xué)學(xué)報(bào),1991,25(5):562-569.
[8] 張芷芬,丁同仁.微分方程定性理論[M].北京:科學(xué)出版社,1985.
【責(zé)任編輯:陳佳】
The (n+ 1) polynomial systems in the existence and uniqueness of limit cycles
CAO Ming
(Department of Mathematics, Shaanxi Xueqian Normal University, Xi′an 710100, China)
Abstract:The current essay tries to investigate the focus at the origin and limit cycles of the kind of (n+1)-th polynomial system based on the Liapunov series theory.We can conclude the origin is a sufficient condition of the focus or center in the system and the Hopf branch problem of the system relies on the parameter δ,and analysis when the parameters meet δma=0 the system the existence and uniqueness of limit cycles.
Key words:polynomial system; Liapunov series; branch; limit cycles
中圖分類號(hào):O175.14
文獻(xiàn)標(biāo)志碼:A
文章編號(hào):1000-5811(2016)02-0169-05
作者簡(jiǎn)介:曹明(1986-),女,陜西咸陽人,講師,研究方向:微分方程與動(dòng)力系統(tǒng)
基金項(xiàng)目:陜西省教育科學(xué)“十二五”規(guī)劃項(xiàng)目(SGH140708); 陜西學(xué)前師范學(xué)院科研基金重點(diǎn)項(xiàng)目(2014ZDKJ010)
收稿日期:2016-01-04