Loading [MathJax]/jax/output/HTML-CSS/jax.js
  数学物理学报  2015, Vol. 35 Issue (1): 83-96   PDF (376 KB)    
扩展功能    
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章    
尚德生
张耀明
一类非光滑对称三次系统的全局分支
尚德生, 张耀明    
山东理工大学理学院, 淄博 255049
摘要    :该文对一类对称三次 Hamilton 系统在非光滑对称摄动下产生的极限环数目进行研究.通过多参数摄动理论和定性分析方法,得到这类在非光滑摄动下的三次系统可以存在至少19个极限环
关键词摄动     非光滑三次系统     同宿轨     极限环    
Analysis of A Diffusive Predator-prey Model with Disease and Holling II Functional Response
SHANG De-Sheng, ZHANG Yao-Ming    
School of Science, Shandong University of Technology,Shandong Zibo 255049
Abstract    : The number of limit cycles of a class of symmetric cubic near-Hamiltonian system under symmetric non-smooth perturbations are investigated in this paper. At least 19 limit cycles are found in this class of perturbed non-smooth cubic system by using the method of multi-parameter perturbation theory and qualitative analysis.
Key words: Perturbation     Non-smooth cubic system     Homoclinic loop      Limit cycle    
1 引言和结果

Hilbert 第十六问题是关于平面n次向量场的极限环个数及其分布的. 若用 H(n) 表示一个给定的n次平面向量场存在极限环的个数, 则迄今为止,我们仅知道一个给定的多项式向量场的环性有限,且 H(2)4, H(3)13,H(4)20,并出现了许多关于这方面的研究, 见文献 [1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17] 等. 一般来说,对一般的平面向量场研究其极限环的个数是非常困难的, 为了规避困难许多论文转向对具有某些对称特性的系统进行研究, 如文献 [12, 15, 16]及其引用的文献等. 在文献[12]中,我们对一类Z2对称的五次系统讨论, 得到系统有 17 个极限环存在,并给出这些环的分布.

为了得到更多的环及这些环的相对位置分布,许多数学工作者提供了很 多简单有效的方法,比如李继斌等人给出了一个用 detection 函数法来摄动对称 Hamilton 系统的方法,见文献[11]等. 还有一个最先由韩茂安在文献[3]中给出的利用研究同宿轨稳定 性发现极限环的方法,后来又通过 文献 [4] 等推广到双同宿轨和异宿轨中,使该方法进一步完善. 运用这一方法的基本思路大体可以分为以下三步:

i)~ 计算决定同宿轨或双同宿轨稳定性的判定量;

ii)~ 通过改变参数的值来改变这些量的符号, 从而改变这些奇闭轨的稳定性,来产生极限环;

iii)~ 最后让同宿轨破裂再出现一个环.

文献 [12, 16, 17]等运用此方法对系统 {˙x=y(1+x2+cy2)+εf(x,y),˙y=x(1ax2y2)+εg(x,y),(1.1)degree(f,g)=3,4,5的情形分别进行研究,得到系统 (1.1) 有11,15,17个极限环的结论.

本文运用这一方法讨论下面的非光滑三次系统 ˙x=y(1+x2+cy2),˙y=x(1ax2y2)+εg(x,y),(1.2) 其中 a<1<c<0,ac>1,而且 ε>0 充分小, g(x,y)=y(a0+a1|x|+a2|y|+a3x2+a4|xy|+a5y2).(1.3)

这里为了便于数值计算,取定 a=4,c=12,且把 (1.3)式 中的 ai 视为参数, 再记 δ=(a0,a1,,a5).

经过对系统(1.2)全面的定性分析得到下面的定理:

定理 1     当 a<1<c<0,ac>1 时, 非光滑三次摄动系统 (1.2) 至少存在 19 个极限环,其分布情况见图 4.

迄今为止,作者还没有看到三次非光滑系统的环性数有更好的结果出现. 论文后面的组织是这样安排的: 在第二节中先给出系统 (1.2) 中当ε=0时的未扰动系统(2.1)的结构,再给出一些全局分析所需的引理.

在第三节中利用第二节得到的结果利用定性分析的方法和多参数扰动理论对系统(1.2)的扰动Melnikov函数的系数符号分析完成定理1的证明.

     严格来说,系统 (1.2) 不是多项式系统扰动问题,所以不能说扰动三次系统可以有 {19} 个极限环存在,因为这里的 g(x,y) 仅是一个分段光滑的函数. 也就是说我们讨论的系统(1.2)是分段光滑系统, 但是这里的g(x,y)及其对变量 y 的偏导数gy(x,y) 都在XY -平面内连续,而且对 x 的偏导数gx(x,y) 也是在除去y轴的集合 R2{x=0}内连续, 这就需要证明对平面解析系统摄动的理论方法推广的非光滑系统 (1.2)上成立. 为此, 我们借鉴了韩茂安在文献[1] 中对同宿轨 Melnikov 函数展开式证明的方法,而得到定理 A,并把这一部分放在最后的附录中. 根据这一结论及其证明,我们知道分段光滑系统(1.2) Melnikov 函数展开式的前三个量的计算公式及相关引理和解析系统的情况相同. 因为对于同宿轨 L3 及其关于原点对称部分的情形可以根据附录中定理 A 的证明得到,而对于其他的扰动同宿轨问题, 我们根据 g(x,y) 在每一象限内解析,从而相应Hamilton系统的同宿轨解析扰动的公式也可以用到系统(1.2)中.

2 未扰系统的结构和必要的引理

ε=0 时,系统 (1.2) 成为 ˙x=y(1+x212y2),˙y=x(1+4x2y2).(2.1) 显然系统(2.1)有七个有限远奇点,分别是 O(0,0), S1(0,2), ˉS1(0,2), A1(12,3), A2(12,3), ˉA1(12,3), ˉA2(12,3),其中 O,S1,ˉS1 是鞍点,而 Ai,ˉAi,i=1,2 为系统(2.1)的中心. 系统(2.1)是以 H(x,y)=12(y2x2)+14(2x2y24x412y4) 为Hamilton函数的 Hamiltion 系统.

对于系统(2.1)的相图和奇点分布情况在文献 [12][16] 中都有描绘, 这里为了便于理解和后面分析, 再次给出其相图如图 1 所示.

图 1 未扰系统(2.1)的主要轨线分布示意图

根据文献 [12][16],得 L1:y=2(x2+1±x2x4),0x1; L2:y=2(x2+1±x2x4),1x0; L3:y=2(x2+1±1+x2x4),x(1+52,0)(0,1+52).

若记 P(x,y)=y(1+x212y2),Q(x,y)=x(1+4x2y2)+εg(x,y), 则系统 (1.2) 可以表示为 ˙x=P(x,y),˙y=Q(x,y).(2.2) 由于系统 (2.2) 关于原点对称,我们只需要对上半平面进行讨论即可,至于下半平面情形可以利用系统对称性得到相同的结构.

对于充分小的 ε>0,系统(2.2)在Lk,k=1,2,3附近有分界线 Lsk,Luk,使得 Ls1Ls2Lu1Lu2 是鞍点S1(ε)的稳定和不稳定流形,Ls3,Lu3 是鞍点Oε的稳定和不稳定流形, 其中 S1(ε)Oε 分别是未扰动系统(2.1)原来的鞍点S1(0,2)O(0,0) 扰动后在其附近产生的鞍点.

我们知道不稳定流形 Luk 和 稳定流形 Lsk 之间的有向距离可以由下面的式子计算 dk=εNkMk(δ)+O(ε2),(2.3) 其中 Nk>0 为常数,Mk 称为系统 (2.2) 的一阶 Melnikov 函数,且 Mk(δ)=Lkg(x,y)dx=6i=1Ai,kai1,(2.4) 这里 k=1,2,3,且 A1,k=Lkydx,A2,k=Lk|x|ydx,A3,k=Lky|y|dx,A4,k=Lkx2ydx,A5,k=Lk|xy|ydx,A6,k=Lky3dx.

利用数学软件Mathematica7.0和简单的计算,得到系数 Ai,k,i=1,2,,6,k=1,2,3 的数值值如下(为了统一起见,我们精确到20位有效数字) A1=A1,i=2laystyle10(1+x2x1x21+x2+x1x2)dx=0.40825717556690997860; A1,3=22laystyle1+520(1+x21+x2x41+x2+1+x2x4)dx=3.5255326669096916993; A2=A2,i=2laystyle10x(1+x2x1x21+x2+x1x2)dx=0.23188188802560381785; A2,3=22laystyle1+520x[1+x21+x2x41+x2+1+x2x4]dx=1.8476046323584734558; A3=B3,i=4laystyle10x1x2dx=43=1.33333333333333333333; A3,3=8laystyle1+5201+x2x4dx=10.0987596075361900248; A4=A4,i=2laystyle10x2(1+x2x1x21+x2+x1x2)dx=0.15345173010947017589; A4,3=22laystyle1+520x2(1+x21+x2x41+x2+1+x2x4)dx=1.36897205570547346050; A5=A5,i=4laystyle10x21x2dx=π4=0.78539816339744830962; A5,3=8laystyle1+5201+x2x4dx=6.0861098394892568386; A6=A6,i=22laystyle10[(1+x2x1x2)3(1+x2+x1x2)3]dx=3.3197985708311904876; A6,3=42laystyle1+520[(1+x21+x2x4)3(1+x2+1+x2x4)3]dx=25.050325702209475317, 这里 i=1,2.

根据上面的系数值,利用隐函数定理得到下面的引理1.

引理 1     存在函数 ϕ1(ε),ϕ4(ε) 使得 d1=d20(<0)  当且仅当  a1(>)ϕ1,(2.5) 并且 d1=d2=0,d30(<0)  当且仅当  a1=ϕ1,a4ϕ4(<ϕ4),(2.6) 其中 ϕ1(ε)=1.7606255453717508235a05.750053808368638102a20.6617667788375365471a33.3870612753969236686a414.316765311418486623a5+O(ε)=ϕ1+O(ε), ϕ4(ε)=1.5863148056671902a03.0555495894870569a2+0.85130624812062093a38.1552261557300819a5+O(ε).

     根据 (2.3),(2.4) 式 利用隐函数定理借助数学软件辅助可得结论. \hfill\rule{0.8mm}{3.5mm}

显然当dk=0,即 ai=ϕi(ε) 时, 在Lk附近存在同宿连接Lk,其中k=1,2,3,i=1,4.

下面对同宿轨 Lk,k=1,2的稳定性进行研究. 在系统(2.1)的鞍点S1附近扰动得到的鞍点用 S1(ε) 表示,并用 σ0,1 表示在扰动鞍点 S1(ε) 处的散度,则有 S1(ε)=(2(a0+2a2+2a5)ε+O(ε2),2+(a0+2a2+2a5)ε+O(ε2)),(2.7)σ0,1=div(S1(ε))=(Px+Qy)|S1(ε)=(a0+22a2+6a5)ε+O(ε2).(2.8)

这样根据 (2.8)式,直接运用隐函数定理得到

引理2     当ai=ϕi(ε),i=1,4 成立时, 存在函数 ϕ0(ε)=22a26a5+O(ε),使得 σ0,10(<0) 当且仅当 a0ϕ0(<ϕ0).

为了研究原点附近的鞍点Oε的同宿轨 L3 及其对称部分ˉL3的稳定性, 利用与前相同的方法可得Oε的坐标(O(ε2),O(ε2)),而 Oε 处的散度(记为σ0,2)为a0ε+O(ε2),故有

引理3     当 ai=ϕi(ε),i=0,1,4 成立时, 存在函数 ϕ2(ε)=322a5+O(ε),使得 σ0,20(<0) 当且仅当 a2ϕ2(>ϕ2).

另外,根据文献 [3, 4, 5],得到下面的引理(详见文献[4]).

引理 4     对于同宿轨Lk,k=1,2,3的稳定性问题,有下述结论成立.

(I)~ 如果 σ0,i>0(<0),i=1,2,则同宿轨 Lk,k=1,2,3 是内不稳定 (稳定)的,而双同宿轨L=L1L2L3ˉL3是外不稳定(稳定)的.

(II)~ 在 ai=ϕi,i=0,1,2,4 成立的条件下,积分 Lj(Px+Qy)dt=σ1,j(δ,ε) 收敛于 σ1,j(δ,0)=Lj(Px+Qy)dt, j=1,2,3.

(III)~ 若 σ1,j(δ,ε)<0(>0),j=1,2,3,则 Lj 是内稳定 (不稳定) 的; 若σ1,1(δ,ε)+σ1,2(δ,ε)>0(<0),则双同宿轨 L 是外不稳定 (稳定) 的.

如果 σ0,1=0 成立,根据 dt=1y(1+x212y2)dx,通过直接计算可得 σ1,j(δ,0)=Li(a0+a1|x|+2a2|y|+a3x2+2a4|xy|+3a5y2)dt=5k=1Bk,jak,  j=1,2,3,(2.9) 其中的系数 Bk,j 有如下计算结果 B1=B1,i=L1|x|y(1+x212y2)dx=22101+x2x1x2+1+x2+x1x21x21+x2+2x4dx=1.9002519479919369908; B1,3=L3|x|y(1+x212y2)dx=21+5201+x2+1+x2x4+1+x21+x2x41+2x21+x2x4dx=3.1110789472030908010; B2=B2,i=L12|y|22y(1+x212y2)dx=21021+x2+2x41+x2x1x21+x2+x1x2x1x21+x2+2x4dx=1.1192315758708453725; B2,3=L32|y|y(1+x212y2)dx=41+52011+x2x4dx=5.6487502839182275990; B3=B3,i=L1x2y(1+x212y2)dx=2210x(1+x2x1x2+1+x2+x1x2)1x21+x2+2x4dx=1.1380237675149813266; B3,3=L3x2y(1+x212y2)dx=21+520x(1+x2+1+x2x4+1+x21+x2x4)1+2x21+x2x4dx=2.2625258492532367229; B4=B4,i=L12|xy|y(1+x212y2)dt=41011x2dx=6.2831853071795864769; B4,3=L32|xy|y(1+x212y2)dx=41+520x1+x2x4dx=4.0688878715914054709; B5=B5,i=L13y26y(1+x212y2)dx=3210(x+1x2)1+x2x1x2+(x1x2)1+x2+x1x21x21+x2+2x4dx=5.8636143559464038513; B5,3=L33y2y(1+x212y2)dx=321+5201+x2+1+x2x4+1+x21+x2x41+x2x4dx=13.9703867746089301822, 这里 i=1,2.

\vskip 2mm

%\hskip 5mm \rule{1mm}{3mm}

%\qed

下面我们要对系统(2.1)的中心Ai,i=1,2经扰动后得到的奇点Ai(ε)(xi0,yi0)及其散度进行计算. 由于结果和方法都相同,所以这里只给出扰动 A1(12,3) 后得到的 A1(ε) 及其散度的计算结果.

(x10,y10)=(12+Aε+O(ε2),3+Bε+O(ε2)),代入系统 (2.2)得 A=14(23a06a16a23a332a463a5), B=14(22a02a126a22a323a462a5). 从而 div(A1(ε))=[Px(x,y)+Qy(x,y)]|(x10,y10)=(a0+22a1+23a2+12a3+6a4+9a5)ε+O(ε2).(2.10)

ai=ϕi,i=0,1,2,4 代入到 (2.9),(2.10)式中,再利用 Mathematica7.0计算可以得到

引理 5     在条件 ai=ϕi,i=0,1,2,4下,对分别同宿于鞍点 S1(ε),和 O(ε) 的同宿轨的散度积分以及 奇点 Ai(ε),i=1,2 处的散度有如下结果 σ1,i=0.2498216168581104a30.28120506343469a5+O(ε),σ1,3=5.302980427691662a3+6.21962317701866a5+O(ε),div(Ai(ε))=(0.07843576661900a3+0.06196848156750a5)ε+O(ε2),i=1,2.(2.11)

为了讨论包围所有有限远奇点的极限环的存在性及其个数问题,我们需要研究对充分大的 b0(>0),系统(2.2)从B(0,b0)出发的正半轨是否有界,也就是要证明正半轨线γ+By轴正半轴的交点 γ+B{x=0,y>0}=B(0,b) 是否满足 b>b0. 为此,对|h|充分大,且h<0时,取 B 为闭曲线 H(x,y)=hy轴正半轴的交点,即H(B)=h,则有 H(B)H(B)=εM4+O(ε2)=εM4+O(ε2),(2.12) 其中 M4=H=hQ(x,y)dx,且 M4=M4|ai=ϕi,i=0,1,2,4.

如果ε>0充分小,且 M4<0 则有不等式 b>b0 成立. 若用 L4 表示闭曲线 H(x,y)=h,为了方便运用数学软件 Mathematica7.0 计算,取 h=100. 把 Q(x,y) 代入到 M4,直接计算可得 M4=6i=1Kiai1,其中 K1=laystyleL4:H=hydx=42[58h+12116h12x2+1x2x4+(12h)dx58h+120x2+1+x2x4+(12h)dx]=84.00984790394565933; K2=laystyleL4:H=h|x|ydx=42[58h+12116h12xx2+1x2x4+(12h)dx58h+120xx2+1+x2x4+(12h)dx]=152.54983565169497551; K3=L4:H=hy|y|dx=8[58h+12116h12(x2+1x2x4+(12h))dx58h+120(x2+1+x2x4+(12h))dx]=531.79575392451528915; K4=laystyleL4:H=hx2ydx=42[58h+12116h12x2x2+1x2x4+(12h)dx58h+120x2x2+1+x2x4+(12h)dx]=364.7222181897085984; K5=L4:H=hy|xy|dx=8[58h+12116h12x(x2+1x2x4+(12h))dx58h+120x(x2+1+x2x4+(12h))dx]=1037.1160146886281834; K6=laystyleL4:H=hy3dx=42[58h+12116h12(x2+1x2x4+(12h))3dx58h+120(x2+1+x2x4+(12h))3dx]=84.00984790394565933.ak=ϕk,k=0,1,2,4 代入到 M4中,得 M4=706.8069206903958a3987.357032085631a5+O(ε).(2.13)

因为在 (2.11) 和 (2.13) 各式中都有两个参数 a3,a5, 我们限定 a3=a5,可得 σ1,i=0.03138344657658a5+O(ε);σ1,3=11.52260360471032a5+O(ε);div(Ai(ε))=0.01646728505150a5ε+O(ε2);M4=280.550111395235a5+O(ε),}(2.14) 其中 i=1,2.

3 系统 (1.2) 的定性分析和定理 1 的证明

为了完成对系统(1.2)的定性分析和定理 1 的证明,我们假设 ε>0 充分小,a3=a5a5>0 固定.

根据 (2.12) 及 (2.14)式得 M4<0,从而系统(1.2)的所有轨道有界, 而且根据(2.14)式得 σ1,i<,σ1,3>0,以及 div(Ai(ε))<0,其中 i=1,2, 故焦点 Ai(ε) 是稳定的,且同宿轨 L1L2 都是内稳定的,而双同宿轨 L 为外稳定的,但 L3ˉL3 都是内不稳定的, 双同宿轨 Γ=L3ˉL3 是外不稳定的. 这样根据环域定理得到在L1L2 内分别存在一个不稳定的极限环,记为γ1,i,i=1,2 (见图 2).

图 2 系统(1.2)在 ε>0 充分小,a5=a3>0,且 ai=ϕi,i=0,1,2,4 条件下的极限环分布示意图

对固定的 a5(=a3)>0 ,如果 0<a2ϕ2a5 成立, 根据引理 3和引理 4 得 σ0,2<0, 这样在同宿轨L3内部以及双同宿轨 L3ˉL3 外部都改变稳定性,即由不稳定变成稳定,从而在 L3内部(但在双同宿轨 L1L2 外部) 有一个不稳定的极限环 γ2,1 出现,而且在双同宿轨 L3ˉL3外部也出现一个大的不稳定极限环 Γ1 (见图 3(a)).

图 3 在参数 a0,a1,a2,a4 相应变化下, 系统 (2.1) 出现极限环的局部结构示意图

对固定的 ak,k=2,3,5,如果 0<a0ϕ0a2ϕ2a5,根据引理 2 和引理 4,得 σ0,1>0,即同宿轨Li内部以及双同宿轨 L1L2 外部的稳定性都由稳定变成不稳定的. 从而分别有稳定极限环 γ3,i,i=1,2 出现在 γ1,iLi,i=1,2之间, 而且还有一个稳定的极限环γ4,1出现在L3γ2,1之间(见图 3(b)).

对固定的 ak,k=0,2,3,5,如果 0<ϕ4a4a0ϕ0a2ϕ2a5,根据引理 1,得 d3<0,从而大同宿轨 L3 破裂,在L3γ4,1 之间出现一个稳定的极限环 γ5,1. 最后,对固定的 ak,k=0,2,,5,如果 0<ϕ1a1ϕ4a4a0ϕ0a2ϕ2a5, 根据引理 1得 d1=d2>0,这样同宿轨 Li,i=1,2 都破裂, 在Liγ3,i 之间又产生一个不稳定的极限环γ6,i (分别见图 3(c),(d)).

这样综合以上分析,并根据系统 (1.2) 的对称性即可完成对系统的分析, 最终得到系统 (1.2) 有19个极限环存在(见图 4).

图 4 系统(1.2)的极限环个数及其分布示意图
附录

在本附录中,考虑 Hamilton 系统的非光滑摄动 {˙x=Hy(x,y)+εp(x,y),˙y=Hx(x,y)+εq(x,y),(A.1) 其中当 x0 时,(p(x,y)q(x,y))=(p+(x,y)q+(x,y)) ; 当x0时,(p(x,y)q(x,y))=(p(x,y)q(x,y))(x,y) 的非光滑函数,但 (p+(x,y)q+(x,y))(p(x,y)q(x,y)) 在其相应的左、右半平面内都解析.

ε=0时,系统 (A.1) 是一个 Hamilton 系统. 假设 Hamilton 函数对 H(x,y)=h0 相应的有同宿于初等鞍点 (设为原点O(0,0)) 的同宿轨 L:z=z0(t),且该同宿轨与y轴相交.

假设 h(z,ε)=(Hy+εp,Hx+εq), 而 A0(0,a0) 是同宿轨Ly轴的交点 (见图 5).

图 5 Hamilton系统的同宿轨在非光 滑扰动下的稳定流形和不稳定流形示意图

ε0 且充分小时,在原点O(0,0)附近有唯一一个双曲鞍点 S(ε,δ) (这里为了方便我们假设原点 O(0,0) 恒为系统的鞍点,即 S(ε,δ) 仍是原点), 在L附近有鞍点S(ε,δ)的稳定和不稳定的分界线Lu(ε,δ), Ls(ε,δ),并用参数表示为 Ls:z=zs(t,ε,δ)=(xs(t,ε,δ),ys(t,ε,δ)),t0,Lu:z=zu(t,ε,δ)=(xu(t,ε,δ),yu(t,ε,δ)),t0.(A.2)t0的情形,相应于左半平面内,可设为 zs(t,ε,δ)=z0(t)+εzs1(t,δ)+O(ε),t0;(A.3)t0的情形,相应的在右半平面内,可设为 zu(t,ε,δ)=z0(t)+εzu1(t,δ)+O(ε),t0.(A.4)

分别将式 (A.3) 和 (A.4) 代入式 (A.1),得 zu,s1 满足 ˙zs1(t)=hz(z0,0)zs1(t)+hε(z0,0),x0;˙zu1(t)=hz(z0,0)zu1(t)+hε(z0,0),x0.Qu,s(t)=r(t),zu,s1(t), 其中 r(t)=(Hx,Hy)|z=z0(t). 则 Qu,s(t) 满足 ˙Qu,s(t)=˙r(t),zu,s1(t)+r(t),hz(z0,0)zu,s1(t)+r(t),hε(z0,0).(A.5) 由于 ˙r(t)=(DHx,DHy)(Hy,Hx)=HxxHyHxyHx,HxyHyHyyHx,而 (r(t),hz(z0,0)zu,s1)=(HxyHxHxxHy,HyyHxHxyHy), 则根据式 (A.5) 得 ˙Qu,s(t)=r(t),hε(z0,0).

由此得 Qu(t)=Qu(0)+t0(qHy+pHx)dt,  Qs(t)=Qs(0)+t0(q+Hy+p+Hx)dt.(A.6)

利用 limtQu(t)=0,limt+Qs(t)=0. 可得如下的位置函数 d(ε,δ,A)=r(0)|r(0)|,zu(0,ε,δ)zs(0,ε,δ)=ε|r(0)|r(0),zu1(0,δ)zs1(0,δ)+o(ε). 从而位置函数满足 d(ε,δ,A)=εM(δ)|(Hy(A)),Hx(A)|+o(ε),(A.7) 其中 M(δ)=Qu(0)Qs(0)=0(qHy+pHx)dt++0(q+Hy+p+Hx)dt=Llqdxpdy+Lrq+dxp+dy,这里 Ll,Lr 分别指同宿轨 L 的左半部分和右半部分.      这里关于 M(δ) 的结果可以由文献 [10] 直接得到,但是为了整体性和后面的讨论, 我们也给出了简要证明.

A(0,a) 是正半轨 Ls 与y轴的交点附近的一点 (在这里根据图 5 中讨论的情况取 a<as). 设 As=A0+as(ε,δ)n0,Au=A0+au(ε,δ)n0,记 C 为过A的正半轨与y轴的第一个交点,再记 B 为过C的正半轨与y轴的第一个交点,则 A=A0+an0. 若再定义 B=A0+P(a,ε,δ)n0, 则当 aas 时,有 au(ε,δ)=P(a,ε,δ).

这样同宿轨扰动后的位置函数d满足 d(ε,δ,A)=au(ε,δ)as(ε,δ)=n0,zu(0,ε,δ)zs(0,ε,δ)=P(a,ε,δ)a=F(a,ε,δ).

完全类似于文献 [1] 中的证明,可得 Fa=1+11+μ(F,a,ε,δ)exp[ACε(p+x+q+y)dt+CBε(px+qy)dt].(A.8) 如果 σ0=div(S(ε,δ))0, 则在我们研究的系统 (A.1) 中有 AC(p+x+q+y)dt=Lr1(p+x+q+y)dt+Lr2(p+x+q+y)dt.(A.9) 但是当aas时,Lr1(p+x+q+y)dt,且Lr2(p+x+q+y)dt 有界.

类似地,有 CB(px+qy)dt=Ll1(px+qy)dt+Ll2(px+qy)dt.(A.10) 且在aas时,Ll1(px+qy)dt, 且 Ll2(px+qy)dt 有界. 这样在aas时,有 AC(p+x+q+y)dt+CB(px+qy)dt. 这样若σ0<0,在aas时, 有Fa1. 所以 F(a,0,δ) 有一个且至多有一个根, 从而同宿轨 L 是稳定的.

说明     这里的 Lr1Ll1 是包含在原点邻域内的稳定和不稳定流形Lu,s的右半部分和左半部分, 而Lr2Ll2是原点外的相应的部分(见图 5). Ll,Lr 是指同宿轨L左、右半部分.

对于σ0>0的情形下的证明方法完全相似,故从略.

下一步我们讨论在 M(δ)=0,及σ0=0 的情形下, 同宿轨L的稳定性.

根据y轴左、右两边函数 (p±q±)(HyHx) 都解析的事实,再由式 (A.8),可得 Fa=1+11+O(F)exp(1τ1+2τ2+ACG1dt+CBG2dt),(A.11) 这里 1=ε(p+x+q+y)|x=y=0, 2=ε(px+qy)|x=y=0, τ1=ACdt,τ2=CBdt,G1=ε(p+x+q+y)1,G2=ε(px+qy)2. 则在条件 M(δ)=0σ0=0 下,有 Fa=1+exp(ACG1dt+CBG2dt+O(F)). 所以 Fa=0 当且仅当 F1=ACG1dt+CBG2dt+O(F)=0. 但是在 (a,ε)(0,0) 时, ACG1dt+CBG2dtσ10. 而在 σ1<0 时,有 F1<0,所以 Fa<0,从而 L 是稳定同宿轨.

根据上述讨论,可得下面的定理 A.

定理 A     对于系统 (A.1),有下面的结论

(1)~ 在L附近存在同宿轨 L 当且仅当位置函数 d(ε,δ,A)=0, 相应的一阶 Melnikov 函数为 M(δ)=Llqdxpdy+Lrq+dxp+dy, 其中 Ll,Lr 分别是同宿轨L的左半部分和右半部分.

(2)~ 如果位置函数 d(ε,δ,A)=0 成立,而 σ0=div(S(ε,δ))0,则在 σ0<0 (相应地,σ0>0)时, 同宿轨 L 是稳定的 (相应地,是不稳定的).

(3)~ 如果位置函数 d(ε,δ,A)=0σ0=0 成立,则 σ1=Lr(p+x+q+y)dt+Ll(px+qy)dt 收敛,且在 σ1<0 (相应地,σ1>0),同宿轨 L 是稳定的 (相应地,是不稳定的).

说明     对于双同宿轨的非光滑摄动情形的证明和 同宿轨情形的证明完全相似,这里就不再重复了.

参考文献
[1] 韩茂安. 动力系统的周期解和分支理论. 北京:科学出版社, 2002
[2] Chen L, Wang M. On relative locations and the number of limit cycles for quadratic systems. Acta Math Sinica, 1979, 22: 751–758
[3] Han M. Cyclicity of planar homoclinic loops and quadratic integrable systems. Sci China (Ser A), 1997,40(12): 1247–1258
[4] Han M, Chen J. On the number of limit cycles in double homoclinic bifurcations. Sci China (Ser A), 2000,43(9): 914–928
[5] Han M, Hu S, Liu X. On the stability of double homoclinic and heteroclinic cycles. Nonlinear Analysis,2003, 53: 701–713
[6] Han M, Shang D, Wang Z, Yu P. Bifurcation of limit cycles in a fourth-order near-Hamiltonian system. International Journal of Bifurcation and Chaos, 2007, 17(11): 4117–4144
[7] Ilyashenko Yu. Finiteness Theorem for Limit Cycles. Providence, RI: American Mathematical Society,1991
[8] Li C Z, Liu C J, Yang J Z. A cubic system with thirteen limit cycles. J Diff Eqs, 2009, 246: 3609–3619
[9] Li J. Hilbert’s 16th problem and bifurcations of planar polynomial vector fields. Int J Bifurcation and Chaos, 2003, 13(1): 47–106
[10] Liang F, Han M, Zhang X. Bifurcation of limit cycles from generalized homoclinic loops in planar piecewise smooth systems. Journal of Differential Equations, 2013, 255(11): 4403–4436
[11] Liu Y, Huang W. A cubic system with twelve small amplitude limit cycles. Bull Sci Math, 2005, 129:83–98
[12] Shang D, Zhang T. Bifurcations of a cubic system perturbed by degree five. Acta Mathematica Scientia,2009, 29B(1): 11–24
[13] Shi S. Aconcrete example of a quadratic system of the existence of four limit cycles for plane quadratic systems. Sci Sinica A, 1980, 23: 153–158
[14] Smale S. Mathematical problems for the next century. The Math Intell, 1998, 20(2): 7–15
[15] Xu W, Han M. On the number of limit cycles of a Z4-equivariant quintic polynomial system. 2010, 216:3022–3034
[16] Zhang T, Zang H, Han M. Bifurcations of limit cycles in a cubic system. Chaos, Solitons and Fractals,2003, 20(3): 629–638
[17] Zang H, Chen W, Zhang T. Perturbation from a cubic Hamiltonian with three eight-loops. Chaos, Solitons and Fractals, 2004, 22: 61–74
一类非光滑对称三次系统的全局分支
尚德生, 张耀明