兩類非線性Schr?dinger型方程的局部臨界周期分支

?

兩類非線性Schr?dinger型方程的局部臨界周期分支

引文格式： 楊劍，黃文韜.兩類非線性Schr?dinger型方程的局部臨界周期分支[J].桂林電子科技大學(xué)學(xué)報，2016,36(1)：71-75.

楊劍1，黃文韜1,2

(1.桂林電子科技大學(xué) 數(shù)學(xué)與計算科學(xué)學(xué)院，廣西 桂林 541004；

2.賀州學(xué)院 數(shù)學(xué)系，廣西 賀州 542800)

摘要：針對非線性Schr?dinger型方程的局部臨界周期分支問題，通過行波變換將兩類Schr?dinger型方程轉(zhuǎn)換為同一等價Hamiltonian系統(tǒng)，應(yīng)用Mathematica計算Hamiltonian系統(tǒng)的周期常數(shù)，得到原點(diǎn)為一階細(xì)中心的充要條件，并證明了該系統(tǒng)在原點(diǎn)鄰域存在一個局部臨界周期分支。分析結(jié)果表明，兩類非線性Schr?dinger型方程均恰有一個局部臨界周期分支，即在原點(diǎn)附近鄰域內(nèi)閉軌周期的單調(diào)性變換一次。

關(guān)鍵詞：周期常數(shù)；細(xì)中心；局部臨界周期分支

在微分方程定性理論中，極限環(huán)和局部臨界周期分支(簡稱臨界周期分支)問題是2個重要的研究方向，具有代表性的研究成果見文獻(xiàn)[1-14]。動力系統(tǒng)臨界周期分支描述了奇點(diǎn)附近鄰域內(nèi)閉軌周期的單調(diào)性變換的次數(shù)，單調(diào)周期函數(shù)已被運(yùn)用于機(jī)械振動和波理論等，如文獻(xiàn)[15-17]研究了閉軌周期單調(diào)性和臨界周期分支問題。

非線性Schr?dinger型方程

(1)

是奧地利物理學(xué)家Schr?dinger[18]于1925年提出的量子力學(xué)的一個基本方程，是將物質(zhì)波和波動方程相結(jié)合建立起來的非線性偏微分方程，被廣泛應(yīng)用于物理、化學(xué)等領(lǐng)域。

考慮兩類Schr?dinger型方程：

(2)

(3)

其中r2、s0、s2、c3、c5為實(shí)常數(shù)。

Guo等[19]用輔助方程法研究了方程(2)的精確行波解，Huang等[20]用雙曲輔助方程法研究了方程(3)的行波解，Zhu[21]用擴(kuò)展的雙曲輔助方程法研究了方程(3)的行波解。Geng等[22]通過行波變換將方程(2)和方程(3)轉(zhuǎn)換為同一等價Hamiltonian系統(tǒng)：

(4)

其中：

其對應(yīng)的Hamiltonian函數(shù)為：

(5)

Geng等[22]應(yīng)用動力系統(tǒng)方法和分支理論對系統(tǒng)(4)進(jìn)行了研究，給出了方程的參數(shù)空間分岔集和相應(yīng)的相圖，并得到了不同參數(shù)條件下所有可能的行波解。

目前，關(guān)于方程(4)奇點(diǎn)鄰域內(nèi)閉軌周期單調(diào)性變換次數(shù)問題的研究還很少。鑒于此，應(yīng)用復(fù)系統(tǒng)周期常數(shù)計算方法，研究系統(tǒng)(4)的中心情況和臨界周期分支問題，并得到了方程(2)和方程(3)的臨界周期分支情況，即周期波解單調(diào)性變換次數(shù)。

1預(yù)備知識

考慮系統(tǒng)原點(diǎn)為非退化中心且在原點(diǎn)鄰域解析的實(shí)系統(tǒng)：

(6)

其中參數(shù)λ∈Λ。

由文獻(xiàn)[23]可知，通過

變換，系統(tǒng)(6)變?yōu)橄铝袕?fù)系統(tǒng)：

(7)

其中：

顯然系統(tǒng)(7)的系數(shù)滿足共軛關(guān)系：

(8)

稱系統(tǒng)(6)與系統(tǒng)(7)互為伴隨。

通過極坐標(biāo)變換x=rcosθ,y=rsinθ，系統(tǒng)(6)變形為：

(9)

其中：

(10)

(11)

由文獻(xiàn)[24]可知，在原點(diǎn)充分小的鄰域內(nèi)，設(shè)P(h,λ)為系統(tǒng)(9)通過(h,0)的閉軌的最小周期函數(shù)，則P(h,λ)局部可積且泰勒展開式為：

(12)

定義1針對系統(tǒng)(6)原點(diǎn)的中心情況，若存在λ*，使得

(13)

則稱原點(diǎn)為k(k∈R)階細(xì)中心；當(dāng)k=0時，原點(diǎn)為粗中心；若存在任一正整數(shù)m，使得P2m=0，則原點(diǎn)為等時中心。

定義2對于參數(shù)λ*∈Λ，系統(tǒng)(6)的原點(diǎn)為細(xì)中心，當(dāng)ε0>0時，對任意ε∈(0,ε0)和任一λ*的充分小鄰域W，有λ1∈W，使得P′(h,λ)=0在U=(0,a)內(nèi)恰好有k個解，則稱λ*系統(tǒng)有k個局部臨界周期分支。

定義3考慮如下集合：

記

V(f1,f2,…,fl)={λ|fi(λ)=0(i=1,2,…,l),λ∈RN}，f∶RN→R,稱f1,f2,…,fl相對于f在λ*∈V(f1,f2,…,fl)上是無關(guān)的。

1)λ*的任意鄰域包含一個λ∈V(f1,f2,…,fl)，使得fl(λ)f(λ)<0。

2)對于變量V(f1,f2,…,fj),2≤j≤l-1，存在λ∈V(f1,f2,…,fj)，fj+1(λ)≠0，且λ的任一鄰域W都包含一個σ∈V(f1,f2,…,fj-1)，使得fj(σ)fj+1(λ)<0。

3)若λ∈V(f1)，f2(λ)≠0，則λ的任一鄰域W*均包含一個σ，使得f1(σ)f2(λ)<0。

顯然，若f1,f2,…,fl相對于f在λ*∈V(f1,f2,…,fl)是無關(guān)的，則任意λ∈V(f1,f2,…,fk-1)使得fk(λ)≠0時，對于每個k，f1,f2,…,fk-1相對于fk是無關(guān)的。

引理1[1]設(shè)系統(tǒng)參數(shù)為λ*時，系統(tǒng)(6)的原點(diǎn)為一個k階細(xì)中心，則至多有k個臨界周期從原點(diǎn)分支出來。若原點(diǎn)的周期常數(shù)P2,P4,…,P2k相對于P2k+2是無關(guān)的，則對于滿足m≤k的任意正整數(shù)m，恰好有m個臨界周期分支。

2細(xì)中心與局部臨界周期分支

由文獻(xiàn)[22]可知，當(dāng)

λ=(d0,d1,d2)∈Ci(i=1,2),

其中,

(14)

系統(tǒng)(4)的原點(diǎn)O(0,0)為中心。

通過

變換，系統(tǒng)(4)可變換為其伴隨復(fù)系統(tǒng)：

(15)

其中：

應(yīng)用文獻(xiàn)[23]中的遞推算法，利用Mathematica軟件分別計算在C1、C2條件下系統(tǒng)(15)原點(diǎn)的復(fù)周期常數(shù)，可得到原點(diǎn)細(xì)中心的階數(shù)及臨界周期分支的個數(shù)。

(16)

其中d0≠0。

由τ1=0可知，d1=0，τ2=τ3=…=0，則系統(tǒng)(4)可變形為：

(17)

顯然系統(tǒng)(17)是線性系統(tǒng)，此時原點(diǎn)為等時中心。

定理2當(dāng)λ=(d0,d1,d2)∈C1時，原點(diǎn)為系統(tǒng)(4)的等時中心的充分必要條件是λ=(d0,d1,d2)∈S1。其中：

(18)

(19)

其中d0≠0,τ1=…=τk-1=0,k=2。

從τ1=0可知，d1=0。由λ=(d0,d1,d2)∈C2可知，d1、d2不同時為零，則τ2≠0。系統(tǒng)(4)的原點(diǎn)至多是一階細(xì)中心，則最多有一個局部臨界周期分支。

定理4當(dāng)λ=(d0,d1,d2)∈C2時，系統(tǒng)(4)的原點(diǎn)為一階細(xì)中心，當(dāng)且僅當(dāng)λ=(d0,d1,d2)∈S2，其中，

(20)

定理5若系統(tǒng)(4)的原點(diǎn)為一階細(xì)中心，則至多有一個局部臨界周期從原點(diǎn)分支出來，且恰有一個局部臨界周期分支。

證明由文獻(xiàn)[24]知，對于系統(tǒng)(4)在原點(diǎn)充分小鄰域內(nèi)通過(h,0)的閉軌周期(最小正周期)函數(shù)可表示為：

(21)

其中：P2k為實(shí)系統(tǒng)(4)原點(diǎn)的第k個周期常數(shù)；τk為伴隨復(fù)系統(tǒng)(15)原點(diǎn)的第一個非零復(fù)周期常數(shù)；P2k與τk滿足P2k=-πτk。

P(h,y)的一階導(dǎo)數(shù)為：

(22)

其中G4(h)為關(guān)于h的二階多項(xiàng)式。

設(shè)d0=1,d1=0,d2=1,可知τ1=0,τ2=0.375,G4(h)無正零點(diǎn)，因此需要通過擾動得到。令ε1=-0.000 000 01,ε2=-0.000 1，有d0=1,d1=0+ε1=-0.000 000 01,d2=1+ε2=0.199 9。在ε1、ε2的擾動下，可知τ1=-1.5×10-8,τ2=0.374 813,經(jīng)計算可知，G4(h)有一個正零解h=0.000 141 457。

3結(jié)束語

通過行波變換將兩類Schr?dinger型方程(2)、方程(3)轉(zhuǎn)換為同一等價Hamiltonian系統(tǒng)(4)，應(yīng)用Mathematica計算該系統(tǒng)的周期常數(shù)，得到了原點(diǎn)為一階細(xì)中心的充要條件，并證明了該系統(tǒng)在原點(diǎn)鄰域存在一個局部臨界周期分支。分析結(jié)果表明，在式(14)成立時，系統(tǒng)(4)原點(diǎn)可通過擾動分支出一個局部臨界周期分支，方程(2)和方程(3)在連續(xù)周期波解u(x,t)≡0鄰域內(nèi)存在一個局部臨界周期分支，即在原點(diǎn)附近鄰域內(nèi)閉軌周期的單調(diào)性變換一次。其他非線性Schr?dinger型方程的中心情況及臨界周期分支問題仍需進(jìn)一步研究。

參考文獻(xiàn):

[1]CHICONEC,JACOBSM.Bifurcationofcriticalperiodsforplanevectorfields[J].TransactionsoftheAmericanMathematicalSociety,1989,312(2):433-486.

[2]SIBIRSKIKS.Onthenumberoflimitcyclesintheneighborhoodofasingularpoint.[J].Differentsial’nyeUravnenija,1965(1):53-66.

[3]LINYiping,LIJibin.Normalformandcriticalpointsoftheperiodofclosedorbitsforplanarautonomoussystems[J].ActaMathematicaSinica,1991,34:490-501.

[4]ROUSSEAUC,TONIB.Localbifurcationofcriticalperiodsinvectorfieldswithhomogeneousnonlinearitiesofthethirddegree[J].CanadianMathematicalBulletin,1993,36(4):473-484.

[5]ROMANOVSKIVG,HANMao’an.Criticalperiodbifurcationsofacubicsystem[J].JournalofPhysicsAGeneralPhysics,2003,36(18):5011-5022.

[6]ZHANGWeinian,HUOXiaorong,ZENGZhenbing.Weakcentersandbifurcationofcriticalperiodsinreversiblecubicsystems[J].Computers&MathematicswithApplications,2000,40(6):771-782.

[7]CHENXingwu,ZHANGWeinian.Decompositionofalgebraicsetsandapplicationstoweakcentersofcubicsystems[J].JournalofComputationalandAppliedMathematics,2009,232(2):565-581.

[8]YUPei,HANMaoan.Criticalperiodsofplanarrevertiblevectorfieldwiththird-degreepolynomialfunctions[J].InternationalJournalofBifurcationandChaos,2009,19(1):419-433.

[9]ROUSSEAUC,TONIB.LocalbifurcationsofcriticalperiodsinthereducedKuklessystem[J].CanadianJournalofMathematics,1997,49(2):338-359.

[10]MANOSASF,VILLADELPRATJ.Anoteonthecriticalperiodsofpotentialsystems[J].InternationalJournalofBifurcationandChaos,2006,16(3):765-774.

[11]YUPei,HANMao’an,ZHANGJizhou.Criticalperiodsofthird-orderplanarhamiltoniansystems[J].InternationalJournalofBifurcationandChaos,2010,20(7):2213-2224.

[12]ZOULan,CHENXingwu,ZHANGWeinian.LocalbifurcationsofcriticalperiodsforcubicLiénardequationswithcubicdamping[J].JournalofComputationalandAppliedMathematics,2008,222(2):404-410.

[13]XUQiujin,HUANGWentao.ThecenterconditionsandlocalbifurcationofcriticalperiodsforaLiénardsystem[J].AppliedMathematicsandComputation,2011,217(15):6637-6643.

[14]CHENTing,HUANGWentao,RENDacheng.WeakcentersandlocalcriticalperiodsforaZ2-equivariantcubicsystem[J].NonlinearDynamics,2014,78(4):2319-2329.

[15]DUZhengdong.OnthecriticalperiodsofLiénardsystemswithcubicrestoringforces[J].InternationalJournalofMathematicsandMathematicalSciences, 2004(61):3259-3274.

[16]GASULLA,ZHAOY.Bifurcationofcriticalperiodsfromtherigidquadraticisochronousvectorfield[J].BulletinDesSciencesMathematiques,2008,132(4):292-312.

[17]CIMAA,GASULLA,DASILVAPR.Onthenumberofcriticalperiodsforplanarpolynomialsystems[J].NonlinearAnalysis:Theory,Methods&Applications,2008,69(7):1889-1903.

[18]SCHR?DINGERE.Anundulatorytheoryofthemechanicsofatomsandmolecules[J].PhysicalReview,1926,28(6):1049-1070.

[19]GUOBolong,WUYaping.Orbitalstabilityofsolitarywavesforthenonlinearderivativeschr?dingerequation[J].JournalofDifferentialEquations,1995,123(1):35-55.

[20]HUANGDingjiang,LIDesheng,ZHANGHongqing.ExplicitandexacttravellingwavesolutionsforthegeneralizedderivativeSchr?dingerequation[J].Chaos,Solitons&Fractals,2007,31(3):586-593.

[21]ZHUShundong.Exactsolutionsforthehigh-orderdispersivecubic-quinticnonlinearSchr?dingerequationbytheextendedhyperbolicauxiliaryequationmethod[J].Chaos,Solitons&Fractals,2007,34(5):1608-1612.

[22]GENGYixiang,LIJibin,ZHANGLixiang.ExactexplicittravelingwavesolutionsfortwononlinearSchr?dingertypeequations[J].AppliedMathematicsandComputation,2010,217(4):1509-1521.

[23]黃文韜.微分自治系統(tǒng)的幾類極限環(huán)分支與等時中心問題[D].長沙：中南大學(xué)，2004.

[24]葉彥謙.多項(xiàng)式微分系統(tǒng)定性理論[M].上海科學(xué)技術(shù)出版社，1995.

編輯：張所濱

Local bifurcation of critical periods for two nonlinear Schr?dinger type equations

YANG Jian1, HUANG Wentao1,2

(1.School of Mathematics and Computational Science, Guilin University of Electronic Technology, Guilin 541004, China；

2.Department of Mathematics, Hezhou University, Hezhou 542800, China)

Abstract：To solve the local bifurcation of critical periods of nonlinear Schr?dinger type equation，the generalized nonlinear derivative Schr?dinger equation and the high-order dispersive nonlinear Schr?dinger equation are converted to a equivalence Hamiltonian system by using traveling wave transformation. The period constants of Hamiltonian system are calculated and the necessary and sufficient condition of the origin with one order weak center is obtained by using Mathematica software. It is proved that there is one local bifurcation of critical periods in the origin of the system．The result shows that the two nonlinear Schr?dinger type equations have only one local bifurcation of critical periods in the neighborhood of near origin, which means that the monotonicity of the periodic for two nonlinear Schr?dinger type equations changes once.

Key words：period constant; weak center; local bifurcation of critical periods

中圖分類號：O175.12

文獻(xiàn)標(biāo)志碼：A

文章編號：1673-808X(2016)01-0071-05

通信作者：黃文韜(1966-)，男，廣西永福人，教授，博士，研究方向?yàn)槲⒎址匠潭ㄐ岳碚摗-mail:huangwentao@163.com

基金項(xiàng)目：國家自然科學(xué)基金(11261013)；廣西自然科學(xué)基金(2012GXNSFAA053003)

收稿日期：2015-06-12