分数阶非线性系统中的共存吸引子_夏鸿鸣.pdf

第4 7卷/第3期/2 0 2 3年5月河北师范大学学报/自然科学版/J O U R N A LO FH E B E IN O R M A LU N I V E R S I T Y(N a t u r a lS c i e n c e)V o l.4 7N o.3M a y.2 0 2 3文章编号:1 0 0 0-5 8 5 4(2 0 2 3)0 3-0 2 3 9-0 6收稿日期:2 0 2 2-0 9-2 2;修回日期:2 0 2 2-1 0-2 8基金项目:国家自然科学基金(1 1 7 0 2 1 9 5);甘肃省自然科学基金(2 2 J R I I R E 1 9 3)作者简介:夏鸿鸣(1 9 6 8),男,甘肃天水人,教授,硕士,研究方向为非线性动力学.分数阶非线性系统中的共存吸引子夏鸿鸣(天水师范学院 数学与统计学院,甘肃 天水 7 4 1 0 0 1)摘要:分数阶混沌系统具有非常有趣和复杂的动力学行为.首先提出了一个新的分数阶非线性系统,并对该系统的一些基本动力学特性进行了数值模拟和理论分析,借助平衡点、相图、分岔图进一步研究了分数阶非线性系统的复杂动力学行为,通过改变初值条件和系统参数,研究发现新的分数阶非线性系统产生不同的共存吸引子.关键词:分数阶;共存吸引子;分岔;混合同步中图分类号:O3 2 2;F4 0 7.6 7 文献标志码:A d o i:1 0.1 3 7 6 3/j.c n k i.j h e b n u.n s e.2 0 2 3 0 2 0 0 6S t u d i e so nC o e x i s t i n gA t t r a c t o r s i nF r a c t i o n a lN o n l i n e a rS y s t e m sX I A H o n g m i n g(S c h o o l o fM a t h e m a t i c sa n dS t a t i s t i c s,T i a n s h u iN o r m a lU n i v e r s i t y,G a n s uT i a n s h u i 7 4 1 0 0 1,C h i n a)A b s t r a c t:F r a c t i o n a l o r d e rs y s t e mh a sm a n y i n t e r e s t i n ga n dc o m p l e xd y n a m i c a lb e h a v i o r s.An e wf r a c-t i o n a l-o r d e rn o n l i n e a r s y s t e mi sp r o p o s e d,a n ds o m eb a s i cd y n a m i c c h a r a c t e r i s t i c so f t h e s y s t e ma r en u m e r-i c a l l ys i m u l a t e da n dt h e o r e t i c a l l ya n a y z e d.F u r t h e r m o r e,v a r i o u sa t t r a c t o r sc a nb eo b t a i n e df r o mt h i sf r a c-t i o n a lo r d e rs y s t e m.T h ec o m p l e xd y n a m i c sb e h a v i o r sa r es t u d i e db ys t e a d yp o i n t,p h a s ep o r t r a i t s,a n db i-f u r c a t i o nd i a g r a m sn u m e r i c a l l y.I na d d i t i o n,t h e c o e x i s t i n ga t t r a c t o r s a r eo b s e r v e db ya d j u s t i n g i n i t i a l v a l u e sa n ds y s t e mp a r a m e t e r s.K e yw o r d s:f r a c t i o n a l o r d e r;c o e x i s t i n ga t t r a c t o r;b i f u r c a t i o n;h y b r i ds y n c h r o n i z a t i o n分数阶微积分是整数阶微积分的推广,已有3 0 0多年的历史,早期因为分数阶微积分缺乏应用背景和计算困难等原因,分数阶微积分理论及应用没有得到实质性进展1.近2 0a来,由于计算机的快速发展,学者发现分数阶微积分有很多的优点,可以更加准确地描述与历史相关的物理变化过程,实际系统具有这种动态特性的现象很多,因此,很多学者在不同的领域采用分数阶系统描述,并得到了很好的结果2-5.另一方面,分数阶混沌系统与对应的整数阶混沌系统相比有着更加丰富的动力学行为,目前,学者研究了很多的分数阶混沌系统并取得了大量的优秀成果6-8.最近几年,很多新的混沌系统不断被发现,如多翼混沌系统、级联混沌系统等,研究表明,许多系统具有多重稳定性,共存吸引子是非线性系统中存在的一种新的动力学行为,通过调整初始状态,系统运行轨线可能渐近趋向点、混沌、周期、准周期等不同的稳定状态,存在共存吸引子的混沌系统在保密通信领域有更好的应用价值9.基于此,本文中,笔者提出了一个新的分数阶混沌系统,不断调整系统参数和初值条件,得到了系统的各种共存吸引子,并进一步研究了此分数阶混沌系统的混合同步控制问题.1 分数阶微积分定义定义1 分数阶微积分是对非整数阶积分 微分算子的推广,如下定义2:aDqt=dqdtq,R(q)0,1,R(q)=0,ta(dt)-q,R(q)0.(1)其中,q表示阶数,R(q)表示q的实部,a和t是积分的上下限.由于C a p u t o分数阶微分定义几何和物理意义明确,在工程实际中应用广泛,因此采用C a p u t o定义.定义2 函数f(x)的C a u p t o分数阶导数定义为 CaDqtf(t)=1(n-q)ta(t-)n-q-1f(n)()d,n-1qn.(2)为了简单起见,关于分数阶非线性系统的表达式中积分下限统一用*代替.2 分数阶系统描述文献1 0 提出了一个整数阶的四维非线性系统,方程如下描述:x=-a x+y z,y=x z-y3,z=-b x y+c z+w,w=y-d z,(3)其中x,y,z,w是状态变量,a,b,c,d是正的参数.结果表明系统(1)具有丰富的动力学行为.考虑系统的阶数为分数阶,则相应的动力学方程表示如下:Dq1*x=-a x+y z,Dq2*y=x z-y3,Dq3*z=-b x y+c z+w,Dq4*w=y-d z.(4)首先,计算系统的平衡点 -a x+y z=0,x z-y3=0,-b x y+c z+w=0,y-d z=0.(5)通过计算可得,系统的平衡点为s=(0,0,0,0),进一步计算雅克比矩阵 JJ=-azy0z-3yx0-b y-b xc101-d0.(6)令d e t(I-JJ)=0,可得下列多项式:f()=4+P33+P22+P1+P0,(7)其中 P3=3y2+a-c,P2=b x2+(3a+b+3c)y2-z2-a c+d,P1=3b y4+a b x2+3(a c+d)y2+c z2+b xyz-x+a d,P0=3a d y2+z2-yz-a x.(8)042当系统参数选取a=6,b=4,c=8,d=2时,计算可得在平衡点的特征值1=0,2=-6,3=-(2 1+4 3 3)/2,4=(2 1-4 3 3)/2.根据稳定性理论,此平衡点是不稳定的.2.1 分数阶系统的混沌和分岔分析当系统参数选取a=6,b=4,c=8,d=2时,调整系统的阶数qi(i=1,2,3,4),这里考虑等阶的情形q1=q2=q3=q4=q.图1 a,b给出了阶数不同时系统的相图.当系统的阶数为q=0.6 4时,系统呈现出周期态,随着阶数的增加,系统产生混沌吸引子.a.q=0.6 4;b.q=0.9 6.图1 阶数不同时,分数阶系统(4)的相图F i g.1 2-DP r o j e c t i o n so fN e wS y s t e m(4)i nt h ex O zP l a n e f o rD i f f e r e n tO r d e r s 为了进一步分析系统的动力学行为,固定系统的阶数,调整系统的参数,利用分岔图研究系统(4)的分岔过程.系统参数选取b=4,c=8,d=2,q=0.9 6时,改变参数2.0a3.2,系统随着参数a变化的分岔图见图2.可以看出,当参数值a3.0 53.2 时,系统呈现混沌吸引子,当参数值a=2.6时,系统发生一系列的倍周期分岔.图2 当a(2.0,3.2)时,分数阶系统(4)的分岔图F i g.2 B i f u r c a t i o nD i a g r a mo f t h eF r a c t i o n a lO r d e rS y s t e mS h o w i n gP e r i o d-d o u b l i n gR o u t e t oC h a o s f o ra(2.0,3.2)系统(4)对应分岔图的二维相图如图3 a,b所示.142qi=0.9 6(i=1,2,3,4);a.a=2.9;b.a=3.2.图3 当a取不同值时,分数阶系统(4)的不同吸引子F i g.3 P h a s eP o r t r a i t so fV a r i o u sA t t r a c t o r s t h eF r a c t i o n a l-o r d e rS y s t e m(4)f o rD e r i v a t i v eO r d e r 接下来,系统参数选取a=6,c=8,d=2,不断调整参数1 1.0b1 5.5,系统(4)的分岔图见图4.图4 分数阶系统(4)随着参数b变化的分岔图F i g.4 B i f u r c a t i o nD i a g r a mo f t h eF r a c t i o n a lO r d e rS y s t e m(4)f o rb(1 1.0,1 5.5)由图4可见,当b=1 2时,系统发生倍周期分岔,其对应的相图见图5 a,b.a.b=1 4;b.b=1 5.5.图5 当b取不同值时,分数阶系统(4)的不同吸引子F i g.5 A t t r a c t o r so f t h eF r a c t i o n a lO r d e rS y s t e mw i t hD i f f e r e n tP a r a m e t e r s242由图5可见,当参数取一定值时,系统(4)从混沌态跳跃到周期态.2.2 分数阶系统的共存吸引子主要考虑改变初始条件和系统参数,分析系统(4)的共存吸引子.为简单起见,固定系统的阶数为q=0.9 6.对于此系统的共存吸引子,分以下几种取值形式.1)系统参数取值为a=6,b=4,c=1.2 5,d=2,初值条件分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=-2,z(0)=3,w(0)=4时,共存吸引子如图6 a所示.可以看出,系统呈现1周期.2)当参数取值为a=6,b=4,c=2,d=2,初值条件分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=-2,z(0)=3,w(0)=4时,系统(4)的共存吸引子如图6 b所示.可以发现,系统呈现4周期.3)当参数取值为a=6,b=1 0.1,c=8,d=2,初值条件分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=2,z(0)=-3,w(0)=4时,共存吸引子如图7 a所示.4)当参数取值为a=8.9,b=4,c=8,d=2,初值条件分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=2,z(0)=-3,w(0)=4时,共存吸引子如图7 b所示.红色和蓝色吸引子对应的初值分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=-2,z(0)=3,w(0)=4;参数值分别为a.a=6,b=4,c=1.2 5,d=2;b.a=6,b=4,c=2,d=2.图6 分数阶系统(4)随着系统参数和初值条件变化的共存吸引子F i g.6 C o e x i s t i n gA t t r a c t o r so fF r a c t i o n a lO r d e rS y s t e m(4)w i t hV a r i e dS y s t e mP a r a m e t e r s红色和蓝色吸引子对应的初值分别为x(0)=1,y(0)=2,z(0)=3,w(0)=4和x(0)=1,y(0)=2,z(0)=-3,w(0)=4;参数值分别为a.a=6,b=1 0.1,c=8,d=2;b.a=8.9,b=4,c=8,d=2.图7 分数阶系统(4)随参数值和初值条件变化的各类共存吸引子F i g.7 V a r i o u sC o e x i s t i n gA t t r a c t o r so fF r a c t i o n a lO r d e rS y s t e m(4)w i t hV a r i e dS y s t e mP a r a m e t e r s342从上述数值模拟中可以看出,当改变系统的初值和参数值时,新的分数阶系统(4)呈现出各类不同的共存吸引子,包括周期、倍周期、混沌吸引子.3 结 论提出了一个新的分数阶非线性系统,并研究了此系统的丰富的动力学行为.主要包括稳定点分析以及不同参数下的相图、系统随着不同的参数变化的分岔图;进一步调整系统的初始条件和参数取值,数值仿真得到了新系统不同状态的共存吸引子.参考文献:1 L ICP,D E NG W H,XUD.C h a o sS y n c h r o n i z a t i o no f t h eC h u aS y s t e m w i t haF r a c t i o n a lO r d e rJ.P h y s i c aA:S t a t i s t i c a lM e c h a n i c sa n dI t sA p p l i c a t i o n s,2 0 0 6,3 6 0(2):1 7 1-1 8 5.d o i:1 0.1 0 1 6/j.p h y s a.2 0 0 5.0 6.0 7 82 P O D L U B N YI.F r a c t i o n a lD i f f e r e n t i a lE q u a t i o n so fM a t h e m a t i c s i nS c i e n c e a n dE n g i n e e r i n gM.C a l i f:A c a d e m i cP r e s s,1 9 9 9:1 9 8.3 YAN GXY,L I U H,L ISG.S y n c h r o n i z a t i o no fF r a c t i o n a l-o r d e r a n dI n t e g e r-o r d e rC h a o t i c(H y p e r-c h a o t i c)S y s t e m sw i t hD i f f e r e n tD i m e n s i o n sJ.A d v a n c e s i nD i f f e r e n c eE q u a t i o n s,2 0 1 7(1):3 4 4.d o i:1 0.1 1 8 6/s 1 3 6 6 2-0 1 7-1 3 9 9-44 史传宝,王光义,臧寿池.一个新的混沌系统及其共存吸引子的研究J.杭州电子科技大学学报(自然科学版),2 0 1 7,3 7(4):1-5.d o i:1 0.1 3 9 5 4/j.c n k i.h d u.2 0 1 7.0 4.0 0 1S H IC h u a n b a o,W A N G G u a n g y i,Z A N GS h o u c h i.A N e wC h a o t i cS y s t e ma n dI t sC o e x i s t i n gA t t r a c t o rJ.J o u r n a lo fH a n g z h o uD i a n z iU n i v e r s i t y(N a t u r a lS c i e n c e),2 0 1 7,3 7(4):1-5.5 S U NKH,W A N G X.B i f u r c a t i o n sa n dC h a o si nF r a c t i o n a l-o r d e rS i m p l i f i e dL o r e n zS y s t e mJ.I n tJB i f u r c a t i o na n dC h a o s,2 0 1 0,2 0:1 2 0 9-1 2 1 9.d o i:1 0.1 1 4 2/s 0 2 1 8 1 2 7 4 1 0 0 2 6 4 1 16 MA G I NRL.F r a c t i o n a lC a l c u l u s i nB i o e n g i n e e r i n gM.N e wY o r k:B e g e l lH o u s e,2 0 0 4:3 2 1-7 7.7 闫少辉,王尔童,孙溪,等.一个吸引子共存的混沌系统及其同步电路实现J.深圳大学学报理工版,2 0 2 1,3 8(6):6 4 9-6 5 7.d o i:1 0.3 7 2 4/s p.j.1 2 4 9.2 0 2 1.0 6 6 4 9YANS h a o h u i,WAN G E r t o n g,S UN X i,e ta l.A C h a o t i cS y s t e m w i t h A t t r a c t o rC o e x i s t e n c ea n dI t sS y n c h r o n i z a t i o nC i r c u i t I m p l e m e n t a t i o nJ.J o u r n a l o fS h e n z h e nU n i v e r s i t y(N a t u r a lS c i e n c e),2 0 2 1,3 8(6):6 4 9-6 5 7.8 鲜永菊,莫运辉,徐昌彪,等.具有多种吸引子共存类型的新型四维混沌系统J.华南理工大学学报(自然科学版),2 0 2 0,4 8(3):3 2-4 3.X I A NY o n g j u,MOY u n h u i,X UC h a n g b i a o,e ta l.N e wF o u r-d i m e n s i o n a lC h a o t i cS y s t e m w i t hM u l t i p l eT y p e so fC o e x i s t e n c eo fA t t r a c t o rJ.J o u r n a l o fS o u t hC h i n aU n i v e r s i t yo fT e c h n o l o g y(N a t u r a lS c i e n c e),2 0 2 0,4 8(3):3 2-4 3.9 CHE N M,WANGC,WU H G,e t.a l.AN o n-a u t o n o m o u sC o n s e r v a t i v eS y s t e ma n dI t sR e c o n s t i t u t i o ni nI n t e g r a lD o m a i nJ.N o n l i n e a rD y n a m i c s,2 0 2 1,1 0 3:6 4 3-6 5 5.d o i:1 0.1 0 0 7/s 1 1 0 7 1-0 2 0-0 6 1 1 5-21 0 L E UT CHOGD,K E N G N EJ,KAMD J E UKL.D y n a m i c a lA n a l y s i so f aN o v e lA u t o n o m o u s4-DH y p e r j e r kC i r c u i tw i t hH y p e r b o l i cS i n eN o n l i n e a r i t yC h a o s,A n t i m o n o t o n i c i t ya n daP l e t h o r ao fC o e x i s t i n g A t t r a c t o r sJ.C h a o sS o l i t o n s&F r a c t a l s,2 0 1 7,1 0 7:6 7-8 7.d o i:1 0.1 0 1 6/j.c h a o s.2 0 1 7.1 2.0 0 8(责任编辑 刘新喜)442