第16卷第6期 2011年12月 文章编号:1007-0249(201 1)06 0066-09 电路与系统学报 JOURNAL 0F CIRCUITS AND SYSTEMS Vo1.16 No.6 December, 201 1 一个新的四维超混沌系统的动力学分析及混沌反同步 黄苏海 , 田立新2 (1.淮海工学院理学院,江苏连云港222005;2.江苏大学非线性科学研究中心,江苏镇江212013) 摘要;提出了一个新的四维Chen—Qi-like混沌系统。通过计算该系统的时间序列的Lyapunov指数谱、Lyapunov 维数、分岔图、Poincar6截面图等分析了控制参数变化时,系统的非线性动力学特征。结果表明该新系统不但和Chen.Qi 系统族有类似的性质,而且又呈现不同的非线性特征。在微分方程不变性原理基础上,运用LMI技术和Riccati方程, 设计了一类新的非线性反馈控制器,实现了超混沌系统的反同步。仿真结果验证了该方法的可行性和有效性。 关键词t新的Chen-Qi-like系统;Lyapunov指数;Poincar6截面图;反同步;黎卡提方程 中图分类号t TP273;023I 文献标识码t A 1 引言 混沌是二十世纪重大发现之一。由于混沌行为是一种非常复杂的动力行为,它具有一些特别的动 力学特性,包括对初值条件极端地敏感、至少有一个正的李雅普诺夫指数、有连续功率谱、至少有一 个分数维等,使其在通信加密技术有很大的应用潜力,而且在电子学、电子系统、密码学、流体混合、 生物系统、图像数据加密、机械震动故障诊断、神经网络等越来越多的领域和方向都有快速的发展。 近年来混沌同步及其控制已成为非线性科学中的一个重要研究课题[1 ]。 在1963年,Lorenz[63在三维自治系统中发现了第一个混沌吸引子。近年来,随着对混沌的深入研究 和实际工程需要,各种非线性混沌系统也被相继提出,并得到了广泛研究[7-10]。 本 ̄Chen.Qi[11,12]混沌系统的基础上,构造了一个新的四维混沌系统。并利用理论推导、数值仿 真、Lyapunov指数谱、Lyapunov维数、分岔图、Poincar6截面图等分析了该系统的基本动力学特性。 最后,在微分方程不变性原理 3,14]基础上,利用Lyapunov稳定性理论 ̄HRiccati方程设计出一种新的非 线性反馈控制器来实现新混沌系统的反同步。 在本文中, 定义为 维Euclidean空间,R 定义为 ×m矩阵。X<0表示x是实对称负正定 矩阵。J定义为恰当维数的的实矩阵。 指的是Euclidean范数。 2 新类Chen—Qi系统的模型及基本动力学特性 根据Qi吸引子和Chen吸引子非线性部分的特征,构造了一个四维非线性动力学系统。系统的模 型如下: =a(y— 、+yz =( 一 )x一 + (1) =X 一bz—vZ—W =xyz一 一 w 其中 =( ,Y,z,w) ∈R 为系统的状态变量, ,b,C,d是正值参数。系统(1)中共含有6个非线性项, 具有很强的非线性特征。 2.1 基本性质 2.1.1对称性和不变性 收稿日期,2010-05—31 修订日期:2011—11-08 基金项目:国家自然科学基金资助项目(】0771088);淮海工学院自然科学基金资助项目(Z2009042) 第6期 黄苏海等:一个新的四维超混沌系统的动力学分析及混沌反同步 67 首先,注意到诸多类Lorenz系统在变换 :。x,Y,z) (一x,-y,Z)下对于所有的参数具有不变性。但是 对于本文提出并研究的系统(1)来说,并不满足上述条件,并且从图2和图3中的相图也能看出,新系 统(1)不再具有对称性,这与以往研究的类Lorenz系统有所不同。 2.1.2耗散性和吸引子的存在性 同时,考虑系统(1)的向量场散度(2),也就是系统的Jacobin矩阵(3)的迹(4) V = =dV dt iV =塞+ + + = (、 = 4 =:ZLES一(2) (3) 口 y 一 0 0 C一 一Z 一J= 2x —b—v 一1 v— —d z+Vz Tr(J):一a+c—b—Y—d (4) 又由于所有Lyapunov特征指数之和反映相空间体积元随时问演化的变化率。根据Liouville定理,变化 率反映为系统的Jacobin矩阵的迹,则有v(t)=V(O)e ,其中,4(i=1,2,3,4)为矩阵(3)的特征 根,LES为系统的4个Lyapunov特征指数。所以只要a+b+d—c+Y>0,则系统(1)始终是耗散的, 并以指数形式收敛。这就意味着,当t ∞时,包含系统轨线的每一个体积元都以指数的速率 一 +b+d—C+Y)收缩No。因此,系统的所有轨线最终都会被限制在一个体积为0的点集合上,并且它 的渐近动力学行为会被固定在一个吸引子上,这就说明了吸引子的存在性。 2.1.3微分同胚和拓扑等价 iO,Chen—Qi系统为: =f(x1 (5) (6) 新系统(1)为: =g(J,) 若系统(5)和系统(6)是微分同胚,则存在一个微分变换Y=h(x)满足: ( )=M~g( ( )) (7) 这里 ( )=_dh(x)是函数 ( )的Jacobian矩阵在点 的取值。令Xo和h(xo)分别是各个系统的平衡点。 A(x。)和s(y。)分别是对应的Jacobian矩阵,则对方程(7)两边求微得到: A(x0)=M ( ) (J,0) ( 0) (8) 因此,矩阵A(x。)和B(y。)的特征多项式应该相同。 不难验证,系统(1)与Chen.Qi系统有不同的特征多项式。因此该系统与Lorenz系统、Chen系统 , ̄lQi系统均不是微分同胚。可以通过严格的数学证明系统(1)与上J ̄Lorenz系统族中每一个系统都不 具有拓扑等价性,是一个完全新类动力学系统。 2.2平衡点及其稳定性分析 令 = =三= =0此时系统只有一个平衡点 (0,0,0,0)。 在 0处对系统(1)线性化,有Jaccobian矩阵: 口 0 c J= 0 0 0 (9) 0一b 一1 0 0 一 令det(J 一 )=0,当口=37,b=3,c=26,d=38,解得矩阵J的特征根: =一=一29.6919,22=18.6919, 3,24=-38,因为 是正实数, , 和 是三个负实数,所以平衡点是一个不稳定鞍点。 第6期 黄苏海等:一个新的四维超混沌系统的动力学分析及混沌反同步 69 控制已成为非线性科学中的一个重要研究内容。本文在微分方程不变性原理基础上,运用LMI技术和 Riccati方程,设计了一类新的非 线性反馈控制器来控制混沌系 统,并在参数未知的情况下实现 了超混沌系统的有限时间反同 步。最后,采用数值方法验证了 控制规律和参数变化规律的正 确性。 3.1 广义方法反同步 x 假定一般混沌系统的动态 方程表示为: =(a)在 y平面上的投影 (b)在 平面上的投影 f(x)+4厂+Or ( ) (11) 其中:X∈ ,f(x)∈ ”,0 ∈R 是未知参数向量, ( )∈R ”, Af∈R 为系统扰动项。(11)作 为驱动系统。 响应系统为: ji,= ( )+U (12) (c)在 —W平面上的投影 (d)在 _z平面上的投影 其中:Y∈R 是输出向量,U∈R 控制输入,h(y)∈R 。 图3系统(1)的超混沌吸引子在不同平面上的投影 定义系统(11)与系统(12) 的反同步状态误差为e=x+Y, 则反同步误差动力学方程表示 为: = ( )+Af+h(y)+0 ( )+ (13) 因此本文所研究的问题可 表述为:驱动系统和响应系统 在不同的初始条件下开始运 行,设计合适的非线性控制器, (a)Poincar6截面图在 平面(z=c) (b)Poincar6截面图在 _z平面 =0) 图4在不同平面上f ̄Poincar6截面图 使得在有限时间内驱动系统与响应系统的状态反同步。即存在时间 ,使m in (oll- ̄o。 一为J,设计合理的控制器,本文百先引八如卜1炭砹: 引理: (Schur线性矩阵不等式) 1l 】, ( ))、 ( )l >0, 等同于 (。 。 。。 )>0, Z(x)一Y(x)W ( )l, ( )>0,这里Z(x)=z ( ),W(x)=W ( )和 Y(x)都依赖于X。 假设1:假定存在一个变换矩阵厂使得,(P)=I I,并且在 厂作用下: ㈥)= l, 图5控关制于参Z轴的数C变分化岔时图 的 70 电路与系统学报 第l6卷 啪 ㈤,_[ : ;: , … ) )= ]∈R 和z= , +:,…, +。J∈ +g: ) Xs= , ,..., ]∈R 和 =Xis+1,Xis+2 ̄" ̄ ̄,Xis+q J 十g=,z) =ei, :,…,e‘(14) ( ) R“ , ( ) R ,其中O, , , 为非线性分解函数,初始值均为0。 假设2:存在矩阵C>0和正定矩阵Ql=Ql >0,A = T,R1=Rl >0,满足Riccati方程: case1: CA +A;C+CR1c+Q1=0,且 co(s,z, ,J') (cR1c)s Case2: CA +A;C+(万1+4)cc+J+Q1=0,有万1,4>0 且 c[o(s,Z,x,y)+ ( , )】≤ 万 +4 J,) I I+,) 根据Schur. ̄,角定理,不等式组等价于如下LMI组: 『CA + ;c+J C ] 、 l c 南 。 假设3:存在矩阵D=diag[d/]>0,-,=1,…,q和正定矩阵Q = >0,R2= >0,A = ;>0满足 Riccati方程 case1:DAz+ArzD+DR2。+Q z 。 ( z (DR2D)Z Case2:DAz+ ;D+( 2+ )DJD+I+Q2=0,有刃2, >0 且z D ( ,z, ,J,)+ ( ,J,)] 1 z 【(万:+ lI( ,J,) lDD+J)Jz 根据Schur ̄角定理,不等式组等价于如下LMI绢: DAz+ D+,D l 。 I 。 假设4:系统扰动项4厂满足: Il4厂II LlIxII 上 I l+三:l :I+…+三 I l 其中厶>0,f:1,2,…, 。 定理1:设定对角阵 =diag ]∈R 满足: (1)自适应律:专=一 , >0,其中 十l’eis+2, ̄..eis+q J,其它 ;0; (2)未知参数0的估计值 自适应律: 一『l )( )T I『I ]DZ I (3)三的估计值三自适应律: 一一 [ ,sTcD!]xXz1U] 其中: : I。 L J 针对系统(13)构造控制器 5 6 (17) (18) (19) 第6期 黄苏海等:一个新的四维超混沌系统的动力学分析及混沌反同步 71 Case1: =一/c 十万c ,一 c 一[ : ]一£『 f+ c2。 Case2: U=一/( )+ ( )一 ( )一ZIxI+ (21 误差系统(13)全局稳定到原点,即系统(1 1) 系统(12)全局反同步。 证明(Case1):令 = 一0,z=三一厶 =diagk+l, +2,…,鼠+ J。 将控制器(20)代入方程(13),得到: :4厂+ ( )+ ( )+ ( )一ZlXI+ ̄e (22) 在变换矩阵,作用下,得到: = + ( ,z, , )+4, 一/-slXsl+ ( ) (23、 =AzZ+ ( ,Z, ,J,)+4,Z.—Lzlx I+ ( ) + z 构建Lyapunov函数: = CS+ZTDZ+ + 三+杰. ( ) (24) J日 s+{ 对V求导: =SrC(A S+O(S,z,x,y)+Afs—Lslx {+ (x) )+( sS+O(S,z,x, )+4 一£ l l+ ( ) )c +z D(A z+ ( ,Z, ,y)+Afz一£zIx I+ ( ) +孑z_)+( zz+ ( ,Z, ,y)+AA一£zIxzI + ( ) + ̄Z)DZ+2 +2 三一2 . ( 一 +,)玑+,z =Sr(CA +A ̄C)S+2S co(s,z, ,y)+2S CAfs一2S cl£ l 1+2s c ( ) +z (DAz+A ̄D)Z+2Z DO(S,Z, ,y)+2Z D4厂Z一2Z zIxz1+2z D ( ) +2 0+2LT/.+2Z 啦。一2 J‘一 、 ( I一占) S+i z J =Sr(CA + ;c+CRlc)s一2S c Ix I+2S c ( ) +Z (DAz+ D+DR2D)Z +2ZrD ̄2(x)O' 姗 2 =Sr(一Q1) +z (一Q2+2De)Z 0(25) 其中 是很小的正数。 3.1.1 实例验证:考虑含未知参数的系统(1)作为驱动系统 I 1=01(x2一 1)+X2 3+4 l 2=( 一 ) 1一X1 3+ 2+4 l 3= 一02 3一 2 3一 4+ (26) I 4=XlX2X3一XlX3一O4x4+4‘ 和系统(27)作为响应系统 I 1=37(y2一Y1)+Y2Y3 』l 2一¨Y夕3=Y1Y2—3y3一Y21一YlY3十2 sgn(6Y2 y2)一Y4 + (27) Ij=,4=Y1Y2Y3一YlY3—38y4 以方程(1 1)和(12)形式, 改写系统(26)和(27) =f(x)+4厂+Ore(x) (28) 电路与系统学报 72 第l6卷 0 『 1 : x1+x2 (x) 0 0 l j 0 (29) j,= ( )+U 37(y2一Y1)+Y2Y3 一llyl—YlY3+26y2 ( ) y1y2-3y3一Y2 sgn(y2)一 4 一一 一 YlY2Y3一YlY3—38Y4 f0 1 0 0 I 选择变换矩阵:J1:= 一 : I:0。1。 I,定义系统(26)与系统(27)的反同步状态误差为: 1 0 0 0 1j r l 一 一令e,: +yl,e2: 2+y2, = + 3, : 4+ 4,注意到: ● 一0 0 O 0 O 37 一37el+e2Y3+x2e3—2x2Y3 l ㈤㈤ 一1le1+26e2一elY3一Xle3+2x ̄y3 l 一卜一1 一 .~ (30) ( )+而(x) Il 一3 一P4+ 1 2一 2 sgn(y2)+ 1 2一 2 sgn( 2)I 38 一elY3一xle3+ 1Y2Y3+x1X2X3+2XlY3 j O X 0 0 3 在J1作用下:S=eI,e,r O O 0 ]●●j●● = 一 .0 0+0 0 0 0 r 37e2—37e1+e2Y3+x2e3—2x2Y3 -]●● ● X z .● 一 1 1●J +Y24y 1,( (J,)十 ( ):--3e3--e一2sg n( Y 2) + x 1xY 2一一- 6e1le+2+2x 2 sgn( 2’I l -38e4一 3一 1P3+YlY2Y3+ l 2 3+2 lY3 j ( ) .I ,一v, sgn。L y_,2,x 十2y,3 -X2 Sgn ]( z)j,J’ J s:[l曩] j ,4 :1L4f ̄ 31, x,y :0 I 1.1457j 系统参数和状态初始值的选取: 数值仿真:以(26)系统作为驱动系统, (27)系统作为响应系统。 第6期 黄苏海等:一个新的四维超混沌系统的动力学分析及混沌反同步 73 ( I,X2,X3, 4)=(一0.1,0.2,一0.5,0.3), (Y1,Y2,Y3,Y4)=(0.7,一0.6,一0.2,0.8) 控制参数: (0)=(3o,0,0,0), 2(0)=0, (0):0, r/1=1,r/2=1 和 (三1,三2,上3,L4)=(0,0,0,0) = =0.5cos ̄t)llxll (a)XI ̄Yl的反同步曲线 (b)X2 ̄Dy2的反同步曲线 +0.3sin(2t), = =0 图6为系统(26)的状态 xi(j=1,2,3,4)和系统(27)的状态 Yi(f=1,2,3,4)之间的反同步仿真曲 线。由仿真结果可以看出两个混沌 系统的状态轨迹振幅相等运动方 向相反,即两个超混沌系统在控制 器(20)的作用下达到了反同步。 由图7可以看出,系统(26)和系 统(27)的相空间中轨迹在二维平 面的投影是反对称的,这亦说明此 时两者已达到了反同步。 结论:同以前不同系统的反同 (c)X3和y3的反同步曲线 (d)X4和Y 的反同步曲线 图6混沌系统(26)和系统(27)反同步过程的仿真曲线 步方法[ 相比,本文提出的控 制方法是基于严格数学推理的基 础上,具有广泛适用性,能够应用 于工程,具有一定的鲁棒性。 4 结论 本文构造并研究了一类新的 Chen.Qi.1ike系统。较为细致地研 究了该系统的一些非线性动力学 (a)xl-x2和yl-y2的反同步曲线 (b)xl-x3和 l 3的反同步曲线 行为,其中包括一些基本的动力学 特征、分岔、周期窗口和通向混沌 道路等。在微分方程不变性原理基 础上,根据LMI技术和Riccati方 程,设计了一类新的非线性反馈控 制器来控制混沌系统,实现了两个 系统在参数未知的情况下的有限 时间反同步。该控制器形式简单, 易于实现且收敛速度快,控制范围 (c)xl-x4. ̄JyI-y4的反同步曲线 (d) 2一 3和y2-y3的反同步曲线 宽。仿真结果验证了该方法的可行 性。 图7系统(26)和系统(27)的相轨迹在二维平面的投影 参考文献; [1] Jae Hun Kim,Chang Woo Park.Fuzzy adaptive synchronization ofuncertain chaotic systems【jr].Physics Letters A,2005,334:295・305. [2】 Ju H Park.synchr0nizati0n ofGenesio chaotic system via back stepping approach[J】_Chaos,Solitons and Fractals,2006,27:1369—1375 74 电路与系统学报 第16卷 [3】 Fa qiang Wang,Chong xin Liu.Synchronization ofunified chaotic system based on passive control[JJ.Physica D,2007,225:55.60. [4] Heng hui Chen.Global synchronization of chaotic systems via linear balanced feedback control[J]l Applied Mathematics and Computation。 2007,186:923-931. [5】 Wen wu Yu,Jin de Cao.Synchronization control of stochastic delayed neural networks[J】Physica A,2007,373:252.260. 【6】 Lorenz E N.Deterministic nonperiodic flow【J】l J Atmos Sci,1963,20:130.141. [7] Elhadj Z—Analysis of a new chaotic system with three quadratic nonlinearities[J].Dynamics of Continuous,Discrete and Impulsive Systems, Series B,2007,14(4):603—610. [8] Guang yi Wang,Hal lian He-A new Rosslor hyperchaotic system and its realization with systematic circuit parameter design[J】l Chinese Physics B,2008,17(11):4014—4021. [9] Guang yi Wang,Shui sheng Qiu,Hong wei Li.A new chaotic system and its circuit realization[J】.Chinese Physics B,2006,15(12): 2872—2877. 0】 王繁珍,齐国元,陈增强,等.一个新的三维混沌系统的分析、电路实现及同步[J].物理学报,2006,55(8):4005—4012. 1] G R Chen,Ucta T.Yet another chaotic attractor【J].International Journal ofBifurcation and Chaos,1999,9:1465.1466. 2】 G Y Qi,G R Chen,S Z Du,et a1.Analysis of a new chaotic system[J1.Physica A,2005,352(2-4):295—308. 3] J P Lasalle The extent of asymptotic stability[J].Proceedings of the National Academy of Sciences USA,1'960,46:363.365. C Wen,Y C Soh,J W Xiao.Adaptive control and synchronization for a class of nonlinear chaotic systems using partial system 4] Y W Wang,states[J].Physics Letters A,2006,35 1:79—84. [15] AI-Sawalha MM,Noorani MSM.Chaos anti。synchronization between two novel different hyperchaotic systems….Chin Phys Lett,2008, 25(8):2743—2747. [16] H Zhang,X K Ma.Synchronization of uncertain chaotic systems with parameters perturbation via active control[J].Chaos Solitons& Fractals,2004,21:39—47. [17] Codreanu S.Synchronization of spatiotemporal nonlinear dynamical systems by an active control[J].Chaos,Solitons&Fractals,2003,15: 507—5 l0. [18】 Vincent UE,Laoye JA.Synchronization,anti synchronization and current transports in non.identical chaotic ratchets[J].Physica A,2007384: ,230.240 作者简介;黄苏海(1973.),男,江苏省扬州市人,博士生, 讲师,主要从事非线性动力系统的研究;m-Or新(1963 ), 男,江苏省姜堰市人,江苏大学理学院教授,博士生导师, 从事无穷维动力系统方向研究。 、 Dynamical analysis and anti—synchrOnizati0n for a new four dimensional hyperchaotic system HUANG Su—hal ,TIAN Li.xin (1.School ofScience,HuaiHai Institute ofTechnology,Lianyungang 222005,China; 2.Nonlinear Scientific Research Center,Jiangsu University,Zhenjiang 212013,China) Abstract:A new Chen—Qi-like four dimensional chaotic system is presented.The nonlinear characteristic of the new svstem versus the control parameters is illustrated by Lyapunov exponent,Lyapunov dimensionbifurcation diagrams and Poincare ,section etc.The results indicate that this new system has some similar characteristics to the Chen—Qi family,and presents some distinct nonlinear properties.A new nonlinear controller is designed based on the invariance principle of differentia1 equations Riccati equations and LMI.With this method,the anti—synchronization of hyperchaotic systems can be achieved .Numerical simulation shows the feasibility and effectivenes of the proposed scbeme. Key words:new Chen—Qi—like system;Lyapunov exponent;Poincard section;anti—synchronization;Riccati equation