余维三fold-fold-Hopf 分岔下簇发振荡及其分类1)

薛 淼 葛亚威 张正娣 毕勤胜

∗(江苏大学土木工程与力学学院,江苏镇江 212013)

†(江苏大学理学院,江苏镇江 212013)

引言

不同尺度耦合具有广泛的工程背景,涉及科学和工程技术的各个领域,如机翼的颤振[1]、绳系卫星中的张弛运动[2]、催化反应中的复合模态振荡[3]、神经元随机簇放电行为[4-8]、忆阻器电路系统混沌振荡[9-13],不同尺度耦合的机械系统复杂的响应特性[14-18]等.多尺度的耦合往往使系统表现出大幅振荡和微幅振荡交替出现,这种现象被称为簇发振荡.不同尺度的耦合,可以来自时间上的快慢[19],也可以来自空间上的大小[20],从而在所建立的相应数学模型中,状态变量或者其某种形式的组合可以分为若干不同的组,而各组之间随时间变化的速率存在着明显量级上的差异[21].

不同时间尺度耦合的工作最早可以追溯到van der Pol 电路方程中的张弛振动(relaxation oscillation)[22],直到诺贝尔奖获得者Hodgkin 和Huxley[23]建立了两时间尺度神经元(HH)模型,成功再现了实测中的神经元放电行为,不同尺度耦合问题逐渐引起重视.从那时起,如何处理不同尺度耦合系统的复杂动力学行为一直是困扰学者们的课题之一[24].围绕这一问题,国内外学者开展了大量的工作.由于缺乏有效的分析方法,早期的相关工作大都是围绕着近似求解开展的,如准静态法[25]、奇异摄动法[26],这些方法虽然能够近似描述系统的运动特性,但是无法解释不同尺度之间的相互作用[27].为此,基于不变流形等概念,以Fenichel 等[28]为代表,建立了几何奇异摄动理论,其核心思想是将系统分解为两个耦合的快慢子系统,通过快子系统的主导和慢子系统的调谐在几何结构上的交替,揭示其复杂运动的产生本质[29],但该方法仅能对相关的簇发振荡产生机制给出大致的结果.随后,以Rinzel[30]为代表的学者发展了这一理论,提出了快慢分析法,利用该方法,可以将簇发振荡按照两种方式分类,一是按照其几何结构的分类,如点-点簇发[31]、点-环簇发[32],二是根据快慢转换的分岔模式分类,如fold/fold[33],fold/Hopf 和Hopf/LPC[34]型簇发振荡.

近年来,基于分岔的分析方法,成为研究不同耦合系统中的簇发振荡的主要手段,如Razvan 和Yasaman[35]探讨了两神经元之间耦合强度对于簇发结构的影响规律,Vo 等[36]分析了反应扩散方程中的多模态吸引子的产生机理,Baldemir 等[37]给出了三维神经元模型中的伪平衡态在沉寂态和激发态转换时的作用,其他还有电路系统中的簇发振荡等.

需要指出的是,Rinzel[30]的快慢分析法仅对自治系统也即时域上的两尺度耦合系统有效,同时,由于高维和高余维系统的复杂性,相关工作大都局限在低维耦合系统存在余维一分岔情形,Izhikevich[38]给出了两快一慢耦合系统可能产生的各种簇发振荡及其机制.而对于非自治系统,如周期激励系统,当激励频率远小于系统的固有频率,也即存在频域两尺度耦合时,系统也可能会产生簇发振荡[39].为揭示这类系统的簇发振荡产生机制,参考几何奇异摄动理论,提出了频域两尺度快慢分析方法[40],并据此得到了诸如周期外激励和周期参数激励下的各种簇发振荡模式及其产生机理[41-43].

迄今为止,大部分结果都是针对具体系统存在余维一分岔时得到的,缺乏系统性.人们不禁要问,在不同分岔模式下,到底存在着哪几种簇发振荡形式,其产生机制又是什么.为此,本文从一类余维三foldfold-Hopf 分岔标准型出发,以慢变周期激励取得某一开折参数,分析该慢变激励值访问不同持久集时存在的各种簇发振荡行为及其机制.需要指出的是,在此仅考虑一开折参数,当对3 个不同开折参数加以分析,则可以得到开折参数空间中所有持久集中的行为对簇发振荡结构的影响,从而可以揭示该分岔点邻域内可能存在的所有簇发振荡模式.

1 频域两尺度快慢分析方法

对于自治快慢耦合系统

式中,x∈RM,y∈RN,µ ∈RK,0 <ε ≪1.传统的快慢分析法,也即奇异几何摄动分析方法,是将慢状态变量y视为分岔参数,从而由隐函数定理可得到快子系统随慢变量变化的平衡曲线,记为x=φ(y,µ),以及其分岔集,记为B(x,y,µ)=0.快子系统的平衡曲线及其吸引子决定耦合系统的主要行为,包含沉寂态(quiescent state,QS)和激发态(spiking state,SP),以及沉寂态和激发态之间的分岔连接,而慢子系统则决定其相应的局部调谐行为[44].通过耦合系统真实轨迹和平衡曲线及其分岔的叠加,得到其相应的分岔机制[45].而对于非自治系统,如周期激励系统

其中x∈RM,µ ∈RK,当周期激励频率远小于系统的固有频率ωN,Ω ≪ωN,也即存在频域上的两尺度耦合时,同样会产生快慢效应[46].由于不存在明显的快慢子系统,传统的分析方法不再适用.为此提出了两种方法,一是直接法[47],引入慢变参数θ=Ωt,也即用代数方程取代微分方程构成慢子系统,从而得到快子系统为

称为广义自治系统.另一为间接法[48],引入慢变参数w=Asin(Ωt),构成慢子系统,从而有广义自治系统

构成快子系统.在具体应用过程中,直接法需要引入一系列放大图,而间接法则相对简便一些,因此,本文采用间接法.由广义自治系统可以得到随w变化的平衡曲线及其分岔集,分别记为

为揭示平衡曲线及其分岔集对系统轨迹的影响机制,引入真实轨迹的转换相图概念[49],定义为

也即将慢变量w视为广义状态变量.通过平衡曲线及其分岔集和耦合系统轨迹的转换相图的叠加,揭示各种快慢效应的产生机制.

2 余维三fold-fold-Hopf 分岔

分岔的余维数由其算子的补空间的维数决定,根据向量场标准型理论,Yu 等[50]推导了平衡点Jacobi 矩阵存在双零及一对纯虚根特征值时相应分岔的标准型及其普适开折,其三阶截断的一般形式可以表示为

其中µi∈R,i=1,2,3 为普适开折参数.当所有非线性项不存在退化情形时,该分岔则可以称为余维三fold-fold-Hopf 分岔.

对于双零特征值的Bogdanov-Takens(BT) 分岔,双开折参数平面可以划分为4 个不同的区域,而对于一对纯虚根的Hopf 分岔,单开折参数可以分为两个区间,不同区域或区间分别对应于不同的动力学行为.而对于含三开折参数的余维三fold-fold-Hopf分岔而言,在三开折参数空间上系统的动力学特性要复杂许多[51].从系统(7) 中可以看到,其平衡点E(x0,y0,u0,v0)满足

上,则可能产生超临界Hopf 分岔,导致周期振荡,其振荡频率分别约为而当参数位于两余维一分岔集的相切点上,则可能产生诸如fold-fold,fold-Hopf,Hopf-Hopf 分岔,同样,也会存在着更高阶的退化,导致余维三分岔甚至余维四分岔.取定参数

图1 给出了广义自治系统在(µ1,µ2)参数平衡上的分岔集.

图1 广义自治系统在(µ1,µ2)参数平面上的分岔集Fig.1 Bifurcation sets of the generalized autonomous system on the(µ1,µ2)plane

分岔集将参数平面划分为3 个区域,对应于不同的动力学行为.在区域(I)中系统在(x,y)平面上表现为稳定的平衡点,且在稳定的平衡点局部存在不稳定的极限环,而在(u,v)平面上表现为稳定的极限环.当参数穿越分岔集HB1进入区域(II)时,稳定的平衡点发生了亚临界Hopf 分岔失稳,且同时不稳定的极限环消失,系统在(x,y)平面上表现为稳定的极限环,而在(u,v)平面上表现为概周期振荡,也即2-D环面.当参数穿越分岔集FB1进入区域(III)时,系统表现出两个稳定的平衡点和一个不稳定的平衡点共存.而当参数从区域(II)穿越分岔集HB2进入区域(IV)时,系统在(x,y)平面上表现为稳定的平衡点,而在(u,v)平面上表现为稳定的极限环.

需要指出的是,在FB1,2上,会产生fold 分岔,导致轨迹趋向(x,y)平面上的不同的稳定平衡点,引起轨迹在不同平衡点之间的跳跃.另外,在fold 分岔和Hopf 分岔曲线的相切点处,会产生余维二Bogdanov-Takens 分岔.

为揭示这些不同区域中的行为对不同尺度耦合系统动力特性的影响,不失一般性,让开折参数µ1随时间作缓慢周期变化,令µ1=Asin(Ωt),其中A和Ω分别对应于激励的振幅和频率,其中令Ω 远小于系统的固有频率.显然,此时系统存在着频域上的两尺度耦合.

显然,µ1∈[-A,A].随着激励幅值A的逐渐增大,系统轨迹访问(µ1,µ2)平面上不同区域的行为仅存在两种可能,也即

情形一:(I)→(II)→(IV);

情形二:(III)→(II)→(IV).

因此,上述两种情形下的簇发振荡即为该分岔向量场标准型截断到三次非线性项时可能的簇发振荡模式.

当然,也可以取开折参数µ2随时间作缓慢周期变化,随A的变化,轨迹访问的区域也会发生变化.在此仅考虑µ1作周期变化.

由频域两尺度分析法,将µ1=Asin(Ωt) 视为慢变参数,构成慢子系统.而基于上述分岔分析,可以得到相应广义自治系统随慢变量变化的平衡曲线和分岔行为.

分别考虑当µ2=-0.4(情形一)和µ2=0.1(情形二)两种典型情形下的尺度效应.图2 分别给出了其相应广义自治系统随慢变量µ1=Asin(Ωt)变化的分岔图.

从图2 可以看出,两种情形下的分岔图几乎接近,但其中依然存在着细微的区别.当µ2=-0.4 时,随着慢变参数µ1的增大,状态变量x沿稳定平衡曲线EB1运动,而在µ1=0 附近处,也即HB1点,穿越亚临界Hopf 分岔系统进入周期振荡,并在HB2点分岔值附近,稳定的极限环消失而稳定于平衡曲线EB3.在平衡曲线上不存在fold 分岔(见图2(a)).需要指出的是,由本文后续的分析,稳定极限环LC分别在两个亚临界Hopf 分岔HB1和HB2的分岔值附近消失的原因都是由LPC 分岔导致稳定的极限环和不稳定的极限环发生碰撞消失.

图2 随µ1 变化的分岔图Fig.2 Bifurcation diagram with the variation ofµ1

当µ2=0.1 时,平衡曲线上存在着两fold 分岔点FB1,2(见图2(b)),当慢变参数µ1增大到µ1=0 时,产生fold 分岔,导致周期振荡,在HB2点分岔值附近,发生LPC 分岔导致稳定的极限环和不稳定的极限环碰撞消失而稳定于平衡曲线EB5.需要指出的是,由本文后续的分析,稳定的极限环LC在fold 分岔值附近消失的原因是发生了鞍点同宿分岔.

进一步从(µ1,µ2) 参数平面上的分岔集看,当µ2=-0.4 时,慢变参数µ1在-A和+A之间作周期变化,系统轨迹可能会进入图1 中的区域(I),(II)和(IV),而µ2=0.1 时,系统轨迹可能会进入图1 中的区域(III),(II)和(IV),从而导致尺度效应的不同.

3 簇发振荡演化及其分岔机制

为分析图1 所示分岔集划分的不同区域中的动力学行为对系统尺度效应的影响,在此分别考虑上述两种情形下随激励幅值增大时系统的动力学行为演化过程.

3.1 情形一:µ2=−0.4

当µ2=-0.4 且A较小时,随着慢变参数µ1的变化,轨迹将周期访问图1 中的区域(I)和(II),使得系统的动力学行为受到两区域中广义自治系统的特性的影响.图3 给出了A=0.10 时两尺度耦合系统在不同平面上的相图及其相应的时间历程.

从图3 可以看出,系统表现为周期簇发振荡,每一周期振荡的轨迹可以分为两部分,分别对应于沉寂态和激发态(见图3(b)).为分析激发态本质特性,定义

图4 给出了轨迹在(u,v) 平面上的Poincar´e映射.从图4 中可以发现,激发态表现为概周期振荡,进一步的数值计算表明,其相应的两振荡频率分别近似于区域(I)中(u,v)平面上稳定极限环的频率和区域(II)中µ1=0.0 附近(x,y)平面上稳定的极限环相应频率.

为揭示簇发振荡的产生机制,图5 给出了(µ1,x)平面上转换相图和平衡曲线的叠加图.假设轨迹从P1点出发,对应于µ1取最小值-0.10,轨迹几乎严格沿着稳定平衡曲线EB1运动,表现为沉寂态.当轨迹运动到亚临界Hopf 分岔点HB1时,由于存在慢变过程的影响,轨迹将穿越该分岔点,直到P2点,亚临界Hopf 分岔的效应才会产生,使得轨迹快速跳向相应的稳定极限环LC,进入激发态.需要指出的是在亚临界Hopf 分岔点附近,稳定的平衡点与不稳定的极限环共存,且不稳定的极限环在亚临界Hopf 分岔点HB1处消失.

图3 A=0.10 时的簇发振荡Fig.3 Bursting oscillations for A=0.10

图4 (u,v)平面上的Poincar´e映射Fig.4 Poincar´e map on the(u,v)plane

图5 A=0.1 时转换相图和平衡曲线(µ1,x)平面上的叠加图Fig.5 Overlap of the transformed phase portrait and equilibrium branches on the(µ1,x)plane for A=0.1

轨迹很快几乎严格按照稳定极限环振荡.当慢变参数µ1增大到最大值0.10 时,轨迹转向.极限环的LPC 分岔导致稳定的极限环和不稳定的极限环碰撞消失使得轨迹振荡幅值逐渐减小,最终趋于稳定平衡曲线EB1.

当轨迹回到P1点时,完成一个周期的簇发振荡.从分岔机制上,该振荡可以称为周期Hopf/LPC 簇发,从几何结构上,则可以称为周期2-D torus-point 簇发.

随着激励幅值的增大,其转换相图将向左右两边扩展,而其簇发振荡的模式保持不变.当激励幅值增加到A=4.881 0 时,在HB2点也会产生Hopf 分岔,导致簇发振荡吸引子的定性改变(见图2(a)).图6 给出了A=8.0 时的簇发振荡吸引子.

从图6 可以看出,系统轨迹同样由Hopf 分岔导致在沉寂态和激发态之间转换,而激发态依然表现为概周期振荡模式(见图7).

图6 A=8.0 时的簇发振荡Fig.6 Bursting oscillations for A=8.0

图6 A=8.0 时的簇发振荡(续)Fig.6 Bursting oscillations for A=8.0(continued)

进一步分析可以发现,在外激励幅值较小时,沉寂态和激发态之间的转换均由µ1=0.0 处亚临界Hopf 分岔及其附近的LPC 分岔引起,而当激励幅值较大时,相应转换由HB1和HB2两点的亚临界Hopf分岔及其附近的两个LPC 分岔引起,这也可以从其转换相图和平衡曲线的叠加图中得到证实(见图8).

图8 A=8.0 时转换相图和平衡曲线(µ1,x)平面上的叠加图Fig.8 Overlap of the transformed phase portrait and equilibrium branches on the(µ1,x)plane for A=8.0

从P1点出发的轨迹几乎严格沿EB1运动,表现为沉寂态.亚临界Hopf 分岔导致轨迹进入激发态,而LPC 分岔则使得轨迹回到沉寂态,使得轨迹几乎严格地沿稳定平衡曲线EB3运动.

簇发振荡的轨迹大致可以分为3 部分,也即分别沿EB1和EB3的运动和按照稳定的极限环LC的振荡,而由Hopf 分岔和LPC 分岔相互转换,因此,从分岔机制上,这种振荡可以称为周期Hopf/LPC/Hopf/LPC 簇发振荡,而从几何结构上则可以称为周期point-torus-point 簇发.

由于广义自治系统在区域(I) 中为稳定的平衡点,在(II)中为2-D 环面运动,而在区域(IV)中则又为稳定的平衡点,当激励幅值较小时,轨迹会连接区域(I)和(II)中的行为,而当激励幅值较大时,轨迹则会表现出3 个区域中的不同特性,从而导致不同的簇发振荡模式,因此,当轨迹访问的区域越多,簇发振荡模式可能会越复杂,因为其受到更多区域中广义自治系统行为的影响.

3.2 情形二:µ2=0.1

增大参数µ2到0.1,随着外激励幅值的变化,系统轨迹将可能访问区域(III)、(II) 和(IV),从而导致簇发振荡模式发生定性改变.图9 给出了A=0.10 时系统的动力学行为.

图9 A=0.1 时的簇发振荡Fig.9 Bursting oscillations for A=0.1

簇发振荡的轨迹大致可以分为两个部分,分别对应沉寂态和激发态.定义截面Σ ≡{(x,y,u,v)|x=1.0,∀t∈R},图10 给出了激发态在截面上的映射,可以看出激发态表现为概周期振荡模式.进一步的数值分析可以发现,两激发态的振荡频率分别近似对应于(x,y)平面上由稳定的极限环的频率和(u,v)平面上极限环的振荡频率.

图10 (u,v)平面上的poincar´e映射Fig.10 Poincar´e map on the plane(u,v)

同样,图11 给出了(µ1,x)平面上的转换相图与平衡曲线的叠加图.从P1点出发的轨迹几乎严格沿稳定平衡曲线EB1运动,表现为沉寂态,当轨迹运动到分岔点时,产生Fold 分岔,平衡点失稳导致轨迹跳向稳定极限环LC,产生大幅振荡,表现为激发态.

图11 A=0.1,µ2=0.1 时转换相图和平衡曲线在(µ1,x)平面上的叠加图Fig.11 Overlap of the transformed phase portrait and equilibrium branches on(µ1,x)plane forA=0.1,µ2=0.1

轨线几乎严格按照稳定的极限环LC振荡,直到慢变参数µ1达到其最大值0.10.轨迹转而向左,按照极限环LC振荡,并最终由鞍点同宿分岔导致稳定的极限环LC破坏消失,跳向稳定平衡曲线EB1.

类似地,由分岔机制,该振荡可以称为周期fold/Homoclinic 簇发振荡.

比较A=0.1,而µ2分别取µ2=-0.4 和µ2=0.1时的两种簇发振荡,虽然两吸引子相似,但本质上是不同的.在µ2=-0.4 时,随着慢变参数µ1的周期变化,簇发振荡(见图3) 的轨迹访问图1 中的区域(I)和区域(II),当(µ1,µ2)位于区域(I)时,由于此时广义自治系统表现为(x,y)平面上的平衡点和(u,v)平面上的极限环,此时轨迹在(x,y)平面上将沿稳定平衡曲线EB1运动,而在(u,v)平面上则按照相应的极限环振荡,这也可以从图12 中x和u的时间历程得到证实.随着µ1的变化在平衡曲线EB1产生Hopf 分岔,导致激发态和沉寂态之间的转换.激发态表现为概周期振荡是因为系统在(x,y)平面上表现出稳定的极限环和(u,v)平面上的极限环的相互作用,对应于区域(II)中的2-D 环面运动.

在µ2=0.1 时,随着慢变参数µ1的周期变化,簇发振荡(见图3)的轨迹访问图1 中的区域(III)和区域(II),当(µ1,µ2)位于区域(III)时,由于此时广义自治系统表现为(x,y)和(u,v)两子平面上均趋于稳定的平衡点,此时轨迹在两子平面上均沿稳定平衡曲线运动,而当(µ1,µ2)位于区域(II)时,轨迹则趋于该域内的稳定环面解,表现为概周期振荡模式.

在µ1=0.0 附近的fold 分岔导致从区域(III)到区域(II)中的二维环面振荡,使得轨迹在沉寂态(平衡点) 和激发态(二维环面) 之间转换,这也可以从图12(b)中x和u的时间历程得到证实.

图12 x 和u 的时间历程Fig.12 Time history of x and u

随着外激励幅值的增大,µ1的变化范围也会增大,一旦穿越图1 中的分岔集HB2时,簇发振荡模式会发生定性改变.图13 给出了A=5.0 时系统的振动特性.

随着A的继续增大,系统的轨迹访问的区域不再发生变化,而其运动模式却会发生改变.图14 给出了A=8.0 时的簇发振荡.

从图13 可以看出,其簇发振荡分别会出现两种形式的沉寂态和激发态,记为QS1,2和S P1,2(见图13(c)).而随着激励幅值的增大,激发态S P2消失,如A=8.0(见图14(c)).

图13 A=5.0 时的簇发振荡Fig.13 Bursting oscillations for A=5.0

图14 A=8.0 时的簇发振荡Fig.14 Bursting oscillations for A=8.0

图14 A=8.0 时的簇发振荡(续)Fig.14 Bursting oscillations for A=8.0(continued)

由于轨迹在图1 中的区域(III),(II)和(IV)来回运动,因此,三区域中广义自治系统的行为将影响簇发振荡的特性.不难想象,在(x,y)子平面上,在(III)和(IV) 中,轨迹趋向其中相应的平衡曲线,而在(II)中,轨迹将趋向其中的二维环面解;在(u,v) 子平面上,在(III)中,轨迹将趋向其中稳定平衡曲线,在(II)中,轨迹将趋向其中在(u,v)的二维环面解,而在(IV)中,轨迹将趋向其中在(u,v)子平面上的稳定极限环.

为分析激发态的振动特性,同样引入截面Σ ≡{(x,y,u,v)|x=1.0,∀t∈R},图15 给出了A=5.0 和A=8.0 时Poincar´e映射在(u,v)平面上的投影,两种情形下Poincar´e映射的投影点均分布于固定圆上,说明簇发振荡的激发态均表现为概周期特性,但两种激发态的具体振动特性存在着一定的差别.

图15 (u,v)平面上的Poincar´e映射Fig.15 Poincar´e map on the(u,v)plane

为分析这些差别的产生原因,转向其转换相图与平衡曲线的叠加图,见图16.

在A=5.0 时,当轨迹从P1点出发,由于轨迹处在图1 中的区域(III) 中,轨迹几乎严格沿其中的稳定平衡曲线EB1运动,表现为沉寂态(见图16(a)).当慢变参数µ1增大到µ1=0.0 时,轨迹进入区域(II),由fold 分岔导致轨迹跳向2-D 环面,进入激发态(见图16(b)).当慢变参数µ1增大到HB2时,轨迹进入区域(IV),LPC 分岔使得稳定的极限环和不稳定的极限环碰撞消失导致激发态的振荡幅值逐渐减小,最终轨迹稳定到该域内的平衡曲线.

当轨迹抵达P2点后,随着时间的继续延长,µ1缓慢减小,轨迹几乎严格沿稳定平衡曲线EB5,当轨迹抵达集HB2时,由于慢变效应,亚临界Hopf 分岔的影响并没有立即出现,轨迹穿越分岔点,转而几乎严格沿不稳定平衡曲线EB4运动,直到抵达P3点,Hopf 分岔的影响才开始出现,使得轨迹跳向区域(II)中的2-D 环面,产生激发态(见图16(b)).

图16 (a)A=5.0,µ2=0.1 时转换相图和平衡曲线(µ1,x)平面上的叠加图;(b)(a)图的局部放大图;(c)A=8.0,µ2=0.1 时转换相图和平衡曲线(µ1,x)平面上的叠加图;(d)(c)图的局部放大图Fig.16 (a)Overlap of the transformed phase portrait and equilibrium branches on the(µ1,x)plane for A=5.0,µ2=0.1;(b)locally enlarged parts in(a);(c)overlap of the transformed phase portrait and equilibrium branches on the(µ1,x)plane for A=8.0,µ2=0.1;(d)locally enlarged parts in(c)

当慢变参数µ1减少到µ1=-0.05 左右时,轨迹经由鞍点同宿分岔导致极限环消失,逐渐稳定于该区域内的稳定平衡曲线EB1,转回沉寂态.

当轨迹回到出发点P1时,一周期的簇发振荡过程结束.从分岔机制上,该运动模式可以称为周期fold/LPC/Hopf/Homoclinic 簇发振荡.

当A=8.0 时,轨迹随当慢变参数µ1从-8.0 增大到8.0 这半个周期内,轨迹运动的模式与A=5.0 时相应段机理完全一样.而在µ1从8.0 减小到-8.0 这半个周期内,轨迹运动的模式表现出不一样的特性.

从P2出发的轨迹,几乎严格沿稳定平衡曲线EB5运动,穿越分岔点HB2,转而几乎严格沿不平衡曲线EB4运动,再穿越鞍点同宿的分岔值,沿稳定平衡曲线EB3运动.这一现象的产生,主要是因为随着A的增大,轨迹沿EB5运动的惯性增大,从而导致两分岔点HB2的亚临界Hopf 分岔效应还来不及出现,轨迹已经穿越区域(II).

沿稳定平衡曲线EB3运动在FB2点由fold 分岔跳向平衡曲线EB1,随后抵达P1点,完成一周期的簇发振荡.

同样,按照轨迹单周期内的分岔行为,该运动则可以称为周期fold/LPC 簇发.

再回到簇发振荡过程中x和u的时间历程上(见图17)的叠加图,从图17 中可以看出,当t∈T1时,x和u均沿平衡曲线运动,当t∈T2,此时轨迹位于区域(IV),x沿平衡曲线运动,而u则按照其中的极限环振荡.

图17 x 和u 的时间历程的叠加图Fig.17 Overlap of the time history of x and u

对于A=5.0,当t∈T3时,轨迹位于区域(II),x和u均按照其中的2-D 环面运动,当t∈T4,轨迹位于区域(III),x和u均沿其中的稳定平衡曲线运动,t∈T5,轨迹位于区域(II),x按照其中的2-D 环面振荡,而u则沿其中的稳定平衡曲线运动,t∈T1时,x和u均沿其中的稳定平衡曲线运动,而当t∈T2时,x沿稳定平衡曲线运动,u按照区域(IV)中的极限环振荡.

而在A=8.0 时,当t∈S1时,轨迹位于区域(IV),x和u均按照稳定平衡曲线运动,t∈S2时,x沿稳定平衡曲线运动,u按照区域(IV) 中的极限环振荡.t∈S3时,x和u均沿区域(III)中稳定平衡曲线运动,t∈S4时,x和u均按照(II)中稳定环面运动振荡.

4 结论与讨论

对于不同的余维分岔,开折参数可以由分岔集划分为不同的区域,分别对应于不同的动力学行为.通过引入周期激励项取代某一开折参数,随着激励幅值的不同,该周期变化的开折参数会穿越不同的区域,从而导致相应区域中的行为交替出现.当激励频率远小于系统的固有频率时,通常会产生簇发振荡.定义整个激励项为慢变量,对于原系统,也就是某一开折参数会在正负激励幅值之间周期缓慢变化,从而导致系统轨迹会访问开折参数空间中的不同区域.这些区域中的各种平衡态及区域边界上的分岔行为会直接影响到簇发振荡的动力特性,一般而言,轨迹访问的区域越多,簇发振荡的行为越复杂.

由于存在着慢过效应,在快慢耦合系统轨迹上的真实分岔,往往与理论上的分岔点之间存在着滞后效应,这种滞后效应的滞后时间也会随之激励幅值的增大而延长,因为激励幅值的增大,会导致轨迹沿相应平衡态运动的惯性增大,特别是,当激励幅值增加到一定值,会导致轨迹沿某平衡态运动并穿越该区域,也即相关分岔效应来不及出现从而导致了不同的簇发振荡类型.本文采用Izhikevich EM 对簇发振荡类型的命名方式,即采用导致沉寂态和激发态交替转迁的分岔来命名簇发振荡的类型.然而,对于沉寂态或激发态过程中发生不会导致转迁行为的某些分岔对簇发振荡模式的影响,尤其是对实际的物理系统的影响,需要进一步研究探讨.

本文仅考虑了双开折参数平面上不同区域内运动对系统快慢效应的影响,对于余维三局部分岔问题,更一般的情形,应该考虑三开折参数空间不同子空间中行为的影响,同时,本文仅考虑了非线性项的三次截断,且非线性项不存在退化行为,当非线项存在退化,会导致分岔余维数的增大,而更高次的非线性项的截断会导致平衡态数目及其分岔类型的增加,使得相应的簇发振荡模式更加丰富.

对于局部分岔下的快慢效应,通过向量场标准型的开折参数的周期扰动,在一定程度可以对该分岔所导致的所有可能的各种簇发振荡进行归类,相关工作尚需进一步深入探讨.

最后,需要指出的是,本文仅考虑了余维三foldfold-Hopf 分岔的标准型及其普适开折的理论上的一般形式,而对于实际系统,当其平衡态产生该类分岔并存在各种扰动时,可以将向量场约化为这一标准型及其普适开折,标准型中的系数和开折参数经过系列推导,可为表示为实际系统各系数的表达式,一般情况下,这些具体表达式可能会非常繁杂.