鐘 溢， 陳鳳娟

(浙江師范大學(xué) 數(shù)理與信息工程學(xué)院,浙江 金華 321004)

0 引 言

同宿軌在動力系統(tǒng)的研究中有著悠久的歷史,可以追溯到Poincaré在研究三體問題時發(fā)現(xiàn)的同宿纏結(jié)現(xiàn)象[1-3].該現(xiàn)象由穩(wěn)定流形與不穩(wěn)定流形橫截相交,從而產(chǎn)生馬蹄意義和似Hénon吸引子意義下的混沌[4-5].令d為穩(wěn)定流形與不穩(wěn)定流形之間的距離,經(jīng)典的Melnikov函數(shù)[6-7]作為d的主項(xiàng),在一定程度上反映了d的特性.文獻(xiàn)[8]推導(dǎo)了高階Melnikov函數(shù),對距離d有了更精細(xì)的刻畫,其中一階Melnikov函數(shù)吻合于經(jīng)典的Melnikov函數(shù).本文通過高階Melnikov函數(shù)分析距離d,研究了同宿軌向量場的一個新的擾動現(xiàn)象.首先給出如下定義:

定義1設(shè)X為光滑流形,F為X上的向量場,且F存在同宿軌.如果存在函數(shù)G,使得向量場F+εG(ε充分小)也存在同宿軌,那么稱G為F的特征擾動函數(shù).

定義2稱所有特征擾動函數(shù)組成的集合為F的特征擾動空間,記為E(F),即

E(F)={G:G為F的特征擾動函數(shù)}.

下面考慮如下二維自治微分方程:

式(1)中:f(x,y),g(x,y),P(x,y),Q(x,y)為有界閉球V?R2上關(guān)于(x,y)的解析函數(shù),且P(x,y),Q(x,y)不顯含時間t.假設(shè)當(dāng)ε=0時,方程(1)存在鞍點(diǎn)O(x0,y0)和同宿到O的同宿軌道(t)={(a(t),b(t)):t∈R}.

式(2)中,

為Melnikov積分.下面將M和d(ε)視為擾動函數(shù)(P(x,y),Q(x,y))的泛函,記為M(P,Q)和d(P,Q),即

由文獻(xiàn)[8]知,式(2)可進(jìn)一步表示為

d(P,Q)=M0(P,Q)ε+M1(P,Q)ε2+M2(P,Q)ε3+…+Mk(P,Q)εk+1+….

(4)

式(4)中:Mk(P,Q)為高階Melnikov函數(shù);k=0,1,2,….由式(4)知,若對任意的非負(fù)整數(shù)k,Mk(P,Q)=0,則d(P,Q)=0.從而(f,g)的特征擾動空間為

E(f,g)={(P,Q):Mk(P,Q)=0,k=0,1,2,…}.

1 主要結(jié)論

下面給出本文的主要結(jié)果.

定理1對于方程(1),設(shè)a(t)是t的奇函數(shù),b(t)是t的偶函數(shù),且f(x,y)是x的偶函數(shù),g(x,y)是x的奇函數(shù),則(f,g)的特征擾動空間為

E(f,g)={(P,Q):P(x,y)是x的偶函數(shù),Q(x,y)是x的奇函數(shù)}.

在證明定理1之前,先給出E(f,g)的性質(zhì).定義(P,Q)的范數(shù)為

‖(P,Q)‖=max{‖P‖0,‖Q‖0}.

定理2(E(f,g),‖5‖)是Banach空間.

證明 分兩步證明.

1)?α,β∈R,?(P1,Q1),(P2,Q2)∈E(f,g),因?yàn)棣罰1+βP2是x的偶函數(shù),αQ1+βQ2是x的奇函數(shù),所以α(P1,Q1)+β(P2,Q2)∈E(f,g).因此,E(f,g)是線性子空間.

2)證明完備性.設(shè){(Pn,Qn)}?E(f,g)為Cauchy序列.由于E(f,g)?C(V)×C(V),且C(V)×C(V)按范數(shù)‖5‖是Banach空間,因此,存在(P,Q)∈C(V)×C(V),使得(Pn,Qn)依范數(shù)收斂于(P,Q),即‖(Pn,Qn)-(P,Q)‖=‖(Pn-P,Qn-Q)‖→ 0,n→∞.于是

‖Pn-P‖0→0, ‖Qn-Q‖0→0,n→∞.

對任意(x,y)∈V,

|(P(-x,y),Q(-x,y))-(P(x,y),-Q(x,y))|≤

|(P(-x,y),Q(-x,y))-(Pn(-x,y),Qn(-x,y))|+

|(Pn(-x,y),Qn(-x,y))-(P(x,y),-Q(x,y))|=

|(Pn(-x,y),Qn(-x,y))-(P(-x,y),Q(-x,y))|+

|(Pn(x,y),-Qn(x,y))-(P(x,y),-Q(x,y))|=

因此,

(P(-x,y),Q(-x,y))=(P(x,y),-Q(x,y)),

即P(x,y)關(guān)于x是偶函數(shù),Q(x,y)關(guān)于x是奇函數(shù).從而(P,Q)∈E(f,g).定理2證畢.

2 Mk(P,Q)的推導(dǎo)

根據(jù)式(4),要證明定理1,只需證對任意的k≥0,Mk(P,Q)=0.令D={(s,z)∈R2:s∈(-∞,+∞),|z|0是一個與ε無關(guān)的小常量.作坐標(biāo)變換

(x,y)=

式(6)中:(s,z)∈D;A(s,z),U(s,z),B(s,z),V(s,z)為

把A(s,z),U(s,z),B(s,z),V(s,z)在(s,z)=(t,0)處作Taylor展開,得

引理11)當(dāng)n是奇數(shù)時,An,m(t),Un,m(t)是t的偶函數(shù);當(dāng)n是偶數(shù)時,An,m(t),Un,m(t)是t的奇函數(shù).

2)當(dāng)n是奇數(shù)時,Bn,m(t),Vn,m(t)是t的奇函數(shù);當(dāng)n是偶數(shù)時,Bn,m(t),Vn,m(t)是t的偶函數(shù).

證明 不失一般性,只證明An,m(t)的奇偶性,其他函數(shù)的奇偶性可類似證明.

根據(jù)定理1的假設(shè),可得上式分子部分是s的奇函數(shù),分母部分是s的偶函數(shù).因此,A(s,z)是s的奇函數(shù).即對于任意的(s,z)∈D,

A(-s,z)=-A(s,z).

等式兩邊關(guān)于s求導(dǎo),根據(jù)鏈?zhǔn)椒▌t得

A1,0(-s,z)=A1,0(s,z);A2,0(-s,z)=-A2,0(s,z).

因?yàn)樽兞縵與s相互獨(dú)立,所以A1,0(s,z)關(guān)于z求導(dǎo)不改變其對s變量的奇偶性.這就意味著對于所有的整數(shù)m≥0,

A1,m(-s,z)=A1,m(s,z);A2,m(-s,z)=-A2,m(s,z).

歸納得

令(s,z)=(t,0),得引理1中的結(jié)論.引理1證畢.

下面求解方程(8).將式(9)代入式(8)得

由常數(shù)變易法,解得方程(11)在穩(wěn)定流形上的解為

類似地,解得方程(11)在不穩(wěn)定流形上的解為

由式(12)和式(13)知,Z+(t),S+(t),Z-(t),S-(t)可寫成如下級數(shù)形式:

其中,K為常數(shù),視具體得到的系數(shù)項(xiàng)而定.

于是,穩(wěn)定流形與不穩(wěn)定流形之間的距離為Z+(t)及Z-(t)在t=0處的值之差,即

引理2對于t≥0,有

證明 運(yùn)用數(shù)學(xué)歸納法.首先,當(dāng)k=0時,

在式(7)中令(s,z)=(t,0)可知,U(t,0)是t的奇函數(shù),B(t,0),V(t,0)是t的偶函數(shù),E(t,τ)關(guān)于t,τ都是偶函數(shù).由積分變量代換得

假設(shè)引理2的結(jié)論對正整數(shù)k

而ρ

于是

結(jié)合引理1得到

因此

引理2證畢.

由式(4)得d(P,Q)=0.于是方程(1)存在同宿軌.根據(jù)定義1 知,(P,Q)為向量場(f,g)的特征擾動函數(shù).從而證得E(f,g)是向量場(f,g)的特征擾動空間.定理1證畢.

3 具體例子

令S1=R/{2nπ},n∈Z.考慮如下系統(tǒng):

式(20)中,(θ,v)∈S1×R.當(dāng)ε=0 時,方程(20)有1個鞍點(diǎn)O(±π,0)和2條同宿到O的同宿軌道±,表達(dá)式為

±=(±2arctan(sinht),±2secht),t∈(-∞,+∞).

比較式(1),此時a(t)=±2arctan(sinht)為t的奇函數(shù),b(t)=±2secht為t的偶函數(shù).f(θ,v)=v,P(θ,v)=0為θ的偶函數(shù),g(θ,v)=-sinθ,Q(θ,v)=vsinθ為θ的奇函數(shù).從而由定理1知,(0,vsinθ)為同宿軌向量場(v,-sinθ)的一個特征擾動函數(shù).圖1為方程(20)當(dāng)ε=0時在柱面S1×R上的相圖,圖2為圖1在yoz平面上的投影.

圖1 方程(20)在ε=0時的同宿軌 圖2 圖1在yoz平面上的投影

當(dāng)ε=0.8時,圖3和圖4分別給出了方程(20)的同宿軌及其在yoz平面上的投影.

圖3 方程(20)在ε=0.8時的同宿軌 圖4 圖3在yoz平面上的投影

當(dāng)ε=-0.8時,圖5和圖6分別給出了方程(20)的同宿軌及其在yoz平面上的投影.

圖5 方程(20)在ε=-0.8時的同宿軌 圖6 圖5在yoz平面上的投影

4 結(jié) 語

同宿軌的不穩(wěn)定性導(dǎo)致了復(fù)雜的同宿纏結(jié)和混沌現(xiàn)象.本文通過高階Melnikov函數(shù)研究了一類向量場自治擾動下的同宿軌持續(xù)性問題,得到了它的特征擾動空間,并證明了它是Banach空間.特征擾動空間對結(jié)構(gòu)不穩(wěn)定向量場的理解具有積極意義,相關(guān)結(jié)果可推廣至異宿軌向量場.