脉冲作用下的生物力学系统非光滑分岔分析.docxVIP

脉冲作用下的生物力学系统非光滑分岔分析 1 基于周期脉冲的分岔控制 1999年，陈关荣等。 现实世界中的许多现象是呈状态突变的,如动物体中的心脏跳动、生物神经网络和病理学中的并发节拍等.这些突变的持续时间与整个运动过程的持续时间相比是非常短暂的,可以认为是瞬间发生的,人们提出了脉冲系统以描述这些特征.另一方面,基于脉冲系统理论的脉冲控制技术具有易于实施、控制成本相对较低等优点,而且有时能取得良好的控制性能,所以脉冲控制技术得到了广泛的应用,如应用于卫星运行轨道的改变、生态系统的管理与控制、混沌系统的脉冲控制与同步等 目前,周期解的稳定性及分岔受到了许多学者的关注 Chen系统的方程为 其中a,b,c为Chen系统的参数. 在Chen系统的第一个方程上加一个周期脉冲,系统可以表示为 其中δ(.)为Diracδ函数,即δ(t)=0,t≠0且 当t=kτ,k∈N时,系统(2)可表示为 并假设 即系统的解(x(t),y(t),z(t))在kτ时刻右连续. 参数a,b,c分别固定为a=35,b=3,c=19.此时未加脉冲作用的Chen系统,即系统(1)具有三个平衡点A(0,0,0),B(3,3,3),C(-3,-3,3),其中A为鞍点,B与C为稳定焦点;且未加脉冲作用的Chen系统具有自然的对称性,即它在变换(x,y,z)→(-x,-y,z)下保持不变.系统初值均取为(-3.5,-3.5,1).未加脉冲作用时Chen系统的时间序列及相图如图1所示. 2 闪映射属性与原系统的关系 首先固定周期τ=1.57,取h作为分岔参数,利用频闪映射研究幅值h对系统的影响.针对h的某个变化范围中的每一个值,求出系统的数值解,然后取出系统状态的频闪映射的点,得到h对系统的影响.对每一个h的值,求出长度为1000个脉冲周期的数值解.然后对每一个h的值,画出该状态的最后200个频闪点.系统(2)随h变化的分岔图如图2所示. 频闪映射是Poincaré映射的一种特殊情形.利用Poincaré映射的性质及其与相应系统的关系,可得频闪映射与原系统具有如下关系 从图2可以发现,当h≤3.55时,系统(2)有1T周期解.当h=0.01时,由于C为稳定焦点,系统经过周期τ后接近焦点,然后受脉冲作用后稍离开焦点,最终形成1T周期解,如图3所示.当h3.55时,由于A为鞍点,B和C为稳定焦点,系统(2)具有3T周期解,如图4所示;然后经6T周期解,12T周期解等倍周期分岔到达混沌,如图5和图6所示.在混沌区域中有周期窗口.当h=4.6718时系统开始由混沌状态突变为2T周期解,如图7所示.当h=9.3511时系统开始由2T周期解突变为混沌状态.在混沌区域中也有周期窗口.当h=10.7767时系统开始由混沌状态突变为1T周期解,如图8所示.当h=16.4755时系统开始由1T周期解经倍周期分岔变为2T周期解.当h=21.815时系统经倒倍周期分岔变为1T周期解. 总之,脉冲作用下的Chen系统向混沌的演化过程可能是倍周期分岔到达混沌,也可能是由周期解直接到达混沌,中间伴随阵发性混沌出现.在混沌区域中,出现不同的周期窗口. 3 脉冲作用下的系统分岔机理 由于脉冲作用下的Chen系统为非光滑系统,故不能用连续系统的分岔分析方法进行分析.下面利用文献[9—11]所提出的非光滑系统的分岔分析方法对脉冲作用下的Chen系统的分岔行为进行分析. 设系统(2)的解为 其中(x 定义映射P P 则整个系统的Poincaré映射P可定义为 即 设Poincaré映射P有平衡点u(x,y,z),即u=T(u).此时系统具有周期τ的解,即1T周期解.相应地,如果Poincaré映射P有k个周期点u(x,y,z),即u=T 映射P (i)当所有特征值的模都小于1时,系统(2)的周期解是渐近稳定的. (ii)当一个特征值通过-1穿过单位圆,而其他特征值的模都小于1时,系统(2)的周期解产生倍周期分岔. (iii)当一个特征值通过1穿过单位圆,而其他特征值的模都小于1时,系统(2)的周期解产生鞍结分岔. (iv)当一对共轭复特征值穿过单位圆,而其他特征值的模都小于1时,系统(2)的周期解产生Hopf分岔. 下面利用Floquet理论揭示脉冲作用下的Chen系统的分岔机理. 3.1 系统周期运动情况 当0h3.5480时,系统所有的Floquet特征乘子都在单位圆内,而在h=3.5480处有一Floquet特征乘子通过1穿出单位圆,如表1所列.根据Floquet理论可判断,当0h3.5480时,系统做稳定的1T周期解运动,在h=3.5480处系统发生鞍结分岔,系统由稳定的1T周期解变为稳定的3T周期解,此结果与图2(b)一致. 当4.6718h≤9.2874时,系统所有的Floquet特征乘子都在单位圆