删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

多种多翼吸引子共存的新型三维分数阶混沌系统

本站小编 哈尔滨工业大学/2020-12-05

多种多翼吸引子共存的新型三维分数阶混沌系统

徐昌彪1,2,何颖辉1,吴霞2,莫运辉2

(1.重庆邮电大学 光电工程学院,重庆 400065; 2.重庆邮电大学 通信与信息工程学院,重庆 400065)



摘要:

自然界的物理现象大多以分数阶的形式存在,整数阶微分方程正好是分数阶微分方程的特例.与整数阶模型相比,分数阶模型更接近真实的世界,具有更诱人的发展前景.为使分数阶混沌系统中共存的多翼混沌吸引子类型更加丰富,提出了一个新型三维分数阶混沌系统,此系统最大的特点是具有多种多翼混沌吸引子共存,即双翼、三翼和四翼混沌吸引子共存.通过相图、Lyapunov指数谱、分岔图等数值仿真,分析了系统的动力学特性,给出了其存在混沌吸引子的必要条件,即q>0.822 4.固定参数,阶数q=0.98时,系统有双翼、双翼、四翼等混沌吸引子共存;q=0.83时,系统有双翼、三翼、四翼等混沌吸引子共存,表明了系统具有丰富的混沌特性.对系统进行了Multisim模拟电路仿真,仿真结果与数值分析相符,进一步验证了其混沌行为.采用分数阶Lyapunov稳定性理论以及定理1,设计了系统的自适应同步控制器,仿真表明响应系统与驱动系统在0.2 s内达到同步,在0.2 s内完成对未知参数的识别,因此,所设计的控制器是有效的.

关键词:  分数阶混沌系统  共存  多翼混沌吸引子  Multisim电路  自适应同步控制

DOI:10.11918/201905213

分类号:TP273

文献标识码:A

基金项目:教育部重点基金(02152)



A novel 3D fractional order chaotic system with coexistence of multiple types of multi-wing attractors

XU Changbiao1,2,HE Yinghui1,WU Xia2,MO Yunhui2

(1.School of Optoelectronic Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China; 2.School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China)

Abstract:

Most of the physical phenomena in the nature exist in the form of fractional order, and integer order equation is just a special case of fractional order equation. Compared with the integer order model, the fractional order model is closer to the real world and has a more attractive development prospect. In order to enrich the types of multi-wing chaotic attractors that coexist in fractional chaotic systems, a novel 3D fractional chaotic system is proposed in this paper. The most significant feature of the system is that there is a coexistence of multiple types of multi-wings chaotic attractors, namely the coexistence of two-wing, three-wing, and four-wing chaotic attractors. The dynamic characteristics of the system were analyzed via numerical simulations of phase diagram, Lyapunov exponent spectrum, and bifurcation diagram. A necessary condition for the existence of chaotic attractors was given, which is q>0.822 4. For fixed system parameters, when q=0.98, there is a coexistence of two-wing, two-wing, and four-wing chaotic attractors in the system; when q=0.83, there is a coexistence of two-wing, three-wing, and four-wing chaotic attractors, which indicates that the chaotic characteristics of the system are complex. The simulation results of an analog circuit of the system by the Multisim were consistent with those by the numerical analysis, which further verified the chaotic behaviors of the system. Based on the fractional Lyapunov stability theory and Theorem 1, an adaptive synchronous controller was designed. The simulation shows that the response system and the drive system reached the synchronization within the range of 0.2 s, and the identification of unknown parameters was completed within the range of 0.2 s, thus the controller is effective.

Key words:  fractional order chaotic system  coexistence  multi-wing chaotic attractors  Multisim circuit  adaptive synchronization control


徐昌彪, 何颖辉, 吴霞, 莫运辉. 多种多翼吸引子共存的新型三维分数阶混沌系统[J]. 哈尔滨工业大学学报, 2020, 52(5): 92-98. DOI: 10.11918/201905213.
XU Changbiao, HE Yinghui, WU Xia, MO Yunhui. A novel 3D fractional order chaotic system with coexistence of multiple types of multi-wing attractors[J]. Journal of Harbin Institute of Technology, 2020, 52(5): 92-98. DOI: 10.11918/201905213.
基金项目 教育部重点基金(02152) 作者简介 徐昌彪(1972-), 男, 教授, 硕士生导师 通信作者 何颖辉, 156089909@qq.com 文章历史 收稿日期: 2019-05-28



Abstract            Full text            Figures/Tables            PDF


多种多翼吸引子共存的新型三维分数阶混沌系统
徐昌彪1,2, 何颖辉1, 吴霞2, 莫运辉2    
1. 重庆邮电大学 光电工程学院, 重庆 400065;
2. 重庆邮电大学 通信与信息工程学院, 重庆 400065

收稿日期: 2019-05-28
基金项目: 教育部重点基金(02152)
作者简介: 徐昌彪(1972-), 男, 教授, 硕士生导师
通信作者: 何颖辉, 156089909@qq.com


摘要: 自然界的物理现象大多以分数阶的形式存在,整数阶微分方程正好是分数阶微分方程的特例.与整数阶模型相比,分数阶模型更接近真实的世界,具有更诱人的发展前景.为使分数阶混沌系统中共存的多翼混沌吸引子类型更加丰富,提出了一个新型三维分数阶混沌系统,此系统最大的特点是具有多种多翼混沌吸引子共存,即双翼、三翼和四翼混沌吸引子共存.通过相图、Lyapunov指数谱、分岔图等数值仿真,分析了系统的动力学特性,给出了其存在混沌吸引子的必要条件,即q>0.822 4.固定参数,阶数q=0.98时,系统有双翼、双翼、四翼等混沌吸引子共存;q=0.83时,系统有双翼、三翼、四翼等混沌吸引子共存,表明了系统具有丰富的混沌特性.对系统进行了Multisim模拟电路仿真,仿真结果与数值分析相符,进一步验证了其混沌行为.采用分数阶Lyapunov稳定性理论以及定理1,设计了系统的自适应同步控制器,仿真表明响应系统与驱动系统在0.2 s内达到同步,在0.2 s内完成对未知参数的识别,因此,所设计的控制器是有效的.
关键词: 分数阶混沌系统    共存    多翼混沌吸引子    Multisim电路    自适应同步控制    
A novel 3D fractional order chaotic system with coexistence of multiple types of multi-wing attractors
XU Changbiao1,2, HE Yinghui1, WU Xia2, MO Yunhui2    
1. School of Optoelectronic Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China;
2. School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China



Abstract: Most of the physical phenomena in the nature exist in the form of fractional order, and integer order equation is just a special case of fractional order equation. Compared with the integer order model, the fractional order model is closer to the real world and has a more attractive development prospect. In order to enrich the types of multi-wing chaotic attractors that coexist in fractional chaotic systems, a novel 3D fractional chaotic system is proposed in this paper. The most significant feature of the system is that there is a coexistence of multiple types of multi-wings chaotic attractors, namely the coexistence of two-wing, three-wing, and four-wing chaotic attractors. The dynamic characteristics of the system were analyzed via numerical simulations of phase diagram, Lyapunov exponent spectrum, and bifurcation diagram. A necessary condition for the existence of chaotic attractors was given, which is q>0.822 4. For fixed system parameters, when q=0.98, there is a coexistence of two-wing, two-wing, and four-wing chaotic attractors in the system; when q=0.83, there is a coexistence of two-wing, three-wing, and four-wing chaotic attractors, which indicates that the chaotic characteristics of the system are complex. The simulation results of an analog circuit of the system by the Multisim were consistent with those by the numerical analysis, which further verified the chaotic behaviors of the system. Based on the fractional Lyapunov stability theory and Theorem 1, an adaptive synchronous controller was designed. The simulation shows that the response system and the drive system reached the synchronization within the range of 0.2 s, and the identification of unknown parameters was completed within the range of 0.2 s, thus the controller is effective.
Keywords: fractional order chaotic system    coexistence    multi-wing chaotic attractors    Multisim circuit    adaptive synchronization control    
分数阶微积分已有300多年的历史,自然界的物理现象大多以分数阶的形式存在,整数阶微分方程正好是分数阶微分方程的特例.与整数阶模型相比,分数阶模型更接近真实的世界,具有更诱人的发展前景,近年来已得到了越来越多的关注[1-2].值得注意的是,与只有固定翼混沌吸引子的分数阶混沌系统相比,具有多种多翼混沌吸引子共存的分数阶混沌系统显示出更复杂的动力学行为和更好的性能[3-4].在安全通信[5]和图像加密[6]中,此类分数阶混沌系统具有更高的序列复杂度以及更大的密钥空间,提高了系统安全性能.因此,发现和构造具有多种多翼混沌吸引子共存的低维分数阶混沌系统具有更大的价值.Zhou等人基于四翼整数阶忆阻混沌系统,构造了相应的分数阶忆阻系统,出现三翼与三翼、三翼与四翼混沌吸引子的共存[7];Borah和Roy设计了一种新的分数阶混沌系统,此系统具有三翼和四翼混沌吸引子共存[8];Xian等人构造了一个双翼与四翼混沌吸引子共存的分数阶混沌系统[9].目前,构造具有更多种多翼混沌吸引子共存的分数阶混沌系统仍然存在一定挑战.在已有文献中,能够产生双翼到四翼混沌吸引子共存的三维分数阶混沌系统还比较少见.

本文设计了一个新型三维分数阶混沌系统,对其进行了动力学分析.在q=0.98时,出现双翼、双翼、四翼等混沌吸引子的共存;q=0.83时,出现双翼、三翼、四翼等混沌吸引子的共存,表明系统具有丰富的混沌特性.与文献[10]和文献[11]相比,此系统不存在不光滑的非线性函数,更易于用硬件电路实现.对系统进行了Multisim模拟电路设计和仿真,仿真结果与数值分析相符,进一步验证了系统的混沌行为.基于分数阶Lyapunov稳定理论,设计了一个自适应同步控制器,仿真结果表明了所设计控制器的有效性.

1 系统模型及其特性 1.1 分数阶系统模型增广Lü系统[12]状态方程为

$\left\{ {\begin{array}{*{20}{l}}{\dot x = - (bc/(b + c))x - yz + d,}\\{\dot y = by + xz,}\\{\dot z = cz + xy.}\end{array}} \right.$ (1)

将式(1)中第一个方程右边添加非线性项-kx-9yz-d,第3个方程右边添加非线性项-9xy+y2.再根据Caputo分数阶积分定义,构建了一个新型三维分数阶混沌系统,系统状态方程为

$\left\{ {\begin{array}{*{20}{l}}{{D^q}x = - ax - 10yz,}\\{{D^q}y = bx + xz,}\\{{D^q}z = cz - 8xy + {y^2}.}\end{array}} \right.$ (2)

式中:xyz为系统状态变量,ab为系统参数,q为系统阶数.

固定参数a=2,b=1,c=-0.78,取q=0.98,其最大Lyapunov指数(MAXLE)为1.065 96.利用预估-校正数值计算方法[13],结合matlab软件,得到系统的仿真相图.初始值为[2, 1, 1]、[5, 1, 1]和[6, 1, 1]时,分别得到两个孤立的双翼混沌吸引子和一个四翼混沌吸引子,其x-y面和y-z面如图 1所示.其中,图 1(a)和(b)中的蓝色部分为初始值[2, 1, 1]诱发的双翼混沌吸引子,红色部分为初始值[5, 1, 1]诱发的四翼混沌吸引子.

Fig. 1
图 1 q=0.98时系统的相 Fig. 1 Phase diagram of the system at q=0.98


q=0.83,其MAXLE为1.494 3.初始值为[2, 1, 1]、[5, 2, 1]和[11, 1, 1]时,分别得到双翼、三翼和四翼混沌吸引子,其x-y面和y-z面如图 2所示.

Fig. 2
图 2 q=0.83时系统的相 Fig. 2 Phase diagram of the system at q=0.83


由此可见,系统阶数q=0.98时,系统有双翼、双翼、四翼等混沌吸引子共存;q=0.83时,系统有双翼、三翼、四翼等混沌吸引子共存,表明了阶数q对系统的混沌特性具有较大影响.

1.2 平衡点及其稳定性参数不变,令系统(2)方程式的左边等于0,即

$\left\{ {\begin{array}{*{20}{l}}{0 = - 2x - 10yz,}\\{0 = y + xz,}\\{0 = - 0.78z - 8xy + {y^2}.}\end{array}} \right.$ (3)

求得平衡点,如表 1所示.

表 1
表 1 系统的平衡点 Tab. 1 Equilibrium points of the system 平衡点 状态变量

x y z

S0 0.303 87 -0.135 9 0.447 2

S1 0 0 0

S2 -0.303 87 0.135 9 0.447 2

S3 -0.321 36 -0.143 7 -0.447 2

S4 0.321 36 0.143 7 -0.447 2



表 1 系统的平衡点 Tab. 1 Equilibrium points of the system


在平衡点$\left( {\begin{array}{*{20}{l}}{\hat x}&, &{\hat y, }&{\hat z}\end{array}} \right)$处线性化系统(2),得Jacobian矩阵

$\mathit{\boldsymbol{J}} = \left[ {\begin{array}{*{20}{c}}{ - 2}&{ - 10\hat z}&{ - 10\hat y}\\{\hat z}&1&{\hat x}\\{ - 8\hat y}&{2\hat y - 8\hat x}&{ - 0.78}\end{array}} \right].$

令det(λEJ)=0,得到其特征多项式为

$f\left( \lambda \right) = {\lambda ^3} + {A_2}{\lambda ^2} + {A_1}\lambda + {A_0}.$ (4)

其中

${A_2} = 1.78,$

${A_1} = 8{{\hat x}^2} + 10{{\hat z}^2} - 80{{\hat y}^2} - 2\hat x\hat y - 1.22,$

${A_0} = 16{{\hat x}^2} + 80{{\hat y}^2} + 7.8{{\hat z}^2} + 20{{\hat y}^2}\hat z - 4\hat x\hat y - 160\hat x\hat y\hat z - 1.56.$

若平衡点稳定,其特征值满足Re[λ] < 0.依据Routh-Hurwitz判据,当且仅当A2>0,A0>0,A2A1-A0>0时,平衡点稳定[14-15].求得每个平衡点所对应的特征值及平衡点的稳定性如表 2所示.

表 2
表 2 系统的平衡点稳定性 Tab. 2 Stability of equilibrium points of systems 平衡点 特征根 稳定性

S0 λ1=-2.636为负实数;λ2=0.428-1.490i,
λ3=0.428+1.490i均有正实部 不稳定

S1 λ1=-2.0, λ2=-0.78为负实数;
λ3=1.0为正实数 不稳定

S2 λ1=-2.636为负实数;λ2=0.428-1.490i,
λ3=0.428+1.490i均有正实部 不稳定

S3 λ1=-2.692为负实数;λ2=0.456-1.452i
λ3=0.456+1.452i均有正实部 不稳定

S4 λ1=-2.692为负实数;λ2=0.456-1.452i
λ3=0.456+1.452i均有正实部 不稳定



表 2 系统的平衡点稳定性 Tab. 2 Stability of equilibrium points of systems


表 2表明,平衡点S1是第一类鞍点,S0S2S3S4是第二类鞍点.混沌吸引子的环围绕第二类鞍点产生,第一类鞍点起到连接环的作用[16-17].

λ是第二类鞍点的共轭复根,为使分数阶系统存在混沌吸引子,λ必须保持在不稳定区域.由文献[18]可得,需满足

$\left| {\arg \left( \lambda \right)} \right| < 0.5{\rm{ \mathit{ π} }}q.$ (5)

其中arg(λ)是特征值λ的辐角,所以有

$q > \frac{2}{{\rm{ \mathit{ π} }}}\left| {\arg \left( \lambda \right)} \right|.$ (6)

由此,a=2,b=1,c=-0.78时,可得分数阶系统(2)存在混沌吸引子的必要条件,即q>0.822 4.

1.3 时序图参数不变,取q=0.83, 初始值为[5, 2, 1],得到相应的时序图,如图 3(a)所示,从图中能直观地看出系统的混沌行为.

Fig. 3
图 3 系统的基本动力学分析 Fig. 3 Basic dynamics analysis of the system


1.4 最大Lyapunov指数谱和分岔图系统阶数的变化会改变平衡点的稳定性, 引起系统状态的变化.参数不变,选择阶数q为控制变量.采用Adams-Bashforth-Moulton数值算法的改进版本[19],初始值为[5, 2, 1]时,得到系统变量z随系统阶数q变化的最大Lyapunov指数谱和分岔图,如图 3(b)和图 3(c)所示.最大Lyapunov指数谱和分岔图可以比较直观地反映非线性动力学系统随阶数变化的动态特性.当最大Lyapunov指数>0的时候,系统处于混沌状态.由图 3(c)可知,当q∈[0.5,0.819]时,系统处于周期状态;当q∈(0.819,1]时,系统处于混沌状态.q从0.5向1增加的过程中,结合图 3(c)可知,系统通向混沌的道路为倍周期分岔道路,进一步验证了q对系统具有较大影响.

2 系统的电路仿真参数不变,取初始值为[5, 2, 1],根据文献[20],当q=0.83时,分数阶算子逼近

$H\left( s \right) = \frac{{3.1869\left( {s + 0.1981} \right)(s + 5.181)(s + 87.72)}}{{(s + 0.01324)(s + 0.299)(s + 6.792)(s + 153.8)}}.$ (7)

文献[19]中的分数阶等效串并联RC模块电路单元如图 4所示.

Fig. 4
图 4 分数阶算子的等效电路模块 Fig. 4 Equivalent circuit module of the fractional order operator


其中Ra=44.310 6 MΩ,Rb=1.307 8 MΩ,Rc=0.085 4 MΩ,Rd=0.005 9 MΩ,Ca=1.709 7 μF,Cb=2.215 7 μF,Cc=1.298 μF,Cd=0.724 5 μF.

采用线性电阻、电容、LM2924N型运算放大器、MULTIPLIER模拟乘法器,其中乘法器的输出增益为1,从而设计出了系统的模拟电子电路,如图 5所示.运算放大器U1A相当于加法器,其输出端电压为vU1A=-(R3/R2)x-(R3/R1)yz,U2A相当于积分器,利用分数阶算子单元模块,得到

Fig. 5
图 5 系统的电路原理 Fig. 5 Circuit schematic diagram of the system


${D^{0.83}}x = {v_{{\rm{U1A}}}} = - \left( {{R_3}/{R_2}} \right)x - \left( {{R_3}/{R_1}} \right)yz,$ (8)

同理

${D^{0.83}}x = {v_{{\rm{U}}4{\rm{A}}}} = \left( {{R_9}/{R_7}} \right)xz + \left( {{R_9}/{R_8}} \right)y,$ (9)

$\begin{array}{*{20}{c}}{{D^{0.83}}x = {v_{{\rm{U}}6{\rm{A}}}} = - \left( {{R_{14}}/{R_{13}}} \right)z - \left( {{R_{14}}/{R_{12}}} \right)xy + }\\{\left( {{R_{14}}/{R_{11}}} \right){y^2}.}\end{array}$ (10)

根据式(2)中的系数,a=2,b=1,c=-0.78,令R1=100 Ω,R2=500 Ω,R3=R7=R8=R9=1 kΩ,R5=R6=R16=R17=10 kΩ,R4=R10=R15=100 kΩ,R11=R14=1.56 kΩ,R12=195 Ω,R13=2 kΩ,Ra1=Ra2=Ra3=44.310 6 MΩ,Rb1=Rb2=Rb3=1.307 8 MΩ,Rc1=Rc2=Rc3=0.085 4 MΩ,Rd1=Rd2=Rd3=0.005 9 MΩ,Ca1=Ca2=Ca3=1.709 7 μF,Cb1=Cb2=Cb3=2.215 7 μF,Cc1=Cc2=Cc3=1.298 μF,Cd1=Cd2=Cd3=0.724 5 μF.

图 6为示波器上观察到的结果,可以看出实验结果与数值仿真结果完全相符,进一步验证了系统的混沌行为.

Fig. 6
图 6 电路实验结果 Fig. 6 Experimental results of the circuit


3 系统的自适应同步采用自适应控制方法[21-23]实现分数阶混沌系统的同步.令驱动系统为

$\left\{ {\begin{array}{*{20}{l}}{{D^{\rm{q}}}{x_1} = - a{x_1} - 10{x_2}{x_3},}\\{{D^{\rm{q}}}{x_2} = b{x_2} + {x_1}{x_3},}\\{{D^{\rm{q}}}{x_3} = c{x_3} - 8{x_1}{x_2} + x_2^2.}\end{array}} \right.$ (11)

其中x1x2x3为系统变量.响应系统为

$\left\{ {\begin{array}{*{20}{l}}{{D^{\rm{q}}}{y_1} = - a{y_1} - 10{y_2}{y_3} + {u_1},}\\{{D^{\rm{q}}}{y_2} = b{y_2} + {y_1}{y_3} + {u_2},}\\{{D^{\rm{q}}}{y_3} = c{y_3} - 8{y_1}{y_2} + y_2^2 + {u_3}.}\end{array}} \right.$ (12)

其中y1y2y3是状态变量,u1u2u3是自适应控制器.

令驱动系统(11)和响应系统(12)的同步误差为ei=yixi(i=1, 2, 3),则同步误差系统为

$\left\{ \begin{array}{l}{D^{\rm{q}}}{e_1} = - a{e_1} - 10{e_2}{e_3} - 10{e_2}{x_3} - \\\;\;\;\;\;\;\;\;\;10{e_3}{x_2} + {u_1},\\{D^{\rm{q}}}{e_2} = b{e_2} + {e_1}{e_3} + {e_1}{x_3} + \\\;\;\;\;\;\;\;\;\;{e_3}{x_1} + {u_2},\\{D^{\rm{q}}}{e_3} = c{e_3} - 8\left( {{e_1}{x_2} + {e_2}{x_1} + {e_1}{e_2}} \right) + \\\;\;\;\;\;\;\;\;\;\;e_2^2 + 2{e_2}{x_2} + {u_3}.\end{array} \right.$ (13)

参数估计误差定义为

$\left\{ {\begin{array}{*{20}{l}}{{e_{\rm{a}}}(t) = a - \tilde a(t),}\\{{e_{\rm{b}}}(t) = b - \tilde b(t),}\\{{e_{\rm{c}}}(t) = c - \tilde c(t).}\end{array}} \right.$ (14)

对参数估计误差(14)求q阶导,得

$\left\{ {\begin{array}{*{20}{l}}{{D^{\rm{q}}}{e_{\rm{a}}}(t) = - {D^{\rm{q}}}\tilde a(t),}\\{{D^{\rm{q}}}{e_{\rm{b}}}(t) = - {D^{\rm{q}}}\tilde b(t),}\\{{D^{\rm{q}}}{e_{\rm{c}}}(t) = - {D^{\rm{q}}}\tilde c(t).}\end{array}} \right.$ (15)

设计系统的自适应控制器为

$\left\{ \begin{array}{l}{u_1} = \tilde a\left( t \right){e_1} + 18{e_3}{x_2} + 17{e_2}{x_3} + \\\;\;\;\;\;\;9{e_2}{x_3} - {k_1}{e_1},\\{u_2} = - \tilde b\left( t \right){e_2} + 7{e_3}{x_1} - {e_2}{e_3} - \\\;\;\;\;\;\;2{e_3}{x_2} - {k_2}{e_2},\\{u_3} = - \tilde c\left( t \right){e_3} - {k_3}{e_3}.\end{array} \right.$ (16)

式中k1, k2, k3为正的增益常数.

定义参数更新定律为

$\left\{ {\begin{array}{*{20}{l}}{{D^{\rm{q}}}\tilde a\left( t \right) = {\lambda _1}\left( { - e_1^2 + {e_{\rm{a}}}} \right),}\\{{D^{\rm{q}}}\tilde b\left( t \right) = {\lambda _2}\left( {e_2^2 + {e_{\rm{b}}}} \right),}\\{{D^{\rm{q}}}\tilde c\left( t \right) = {\lambda _3}\left( {e_3^2 + {e_{\rm{c}}}} \right).}\end{array}} \right.$ (17)

式中λ1, λ2, λ3为正的增益常数.

对于分数阶混沌系统,由参考文献[24]给出的分数阶Lyapunov稳定性理论,得出引理1.

引理1??对于分数阶系统${D^{\rm{q}}}\mathit{\boldsymbol{X}}(t) = f(\mathit{\boldsymbol{X}})$,其中X=(x1, x2, …, xn)T为系统状态变量.当系统阶数q∈(0, 1)时,如果存在实对称正定矩阵P,使对任意状态变量矩阵X,恒有$J = \mathit{\boldsymbol{XP}}{D^9}\mathit{\boldsymbol{X}} \le 0$成立,则分数阶系统渐近稳定.

定理1??在自适应控制器(16)和参数更新定律(17)的作用下,同步误差系统(13)渐近稳定.

证明??构造J函数如下:

$J = \mathit{\boldsymbol{eP}}{D^{\rm{q}}}e = \left[ {e,{e_{\rm{a}}},{e_{\rm{b}}},{e_{\rm{c}}}} \right]\mathit{\boldsymbol{P}}{\left[ {{D^{\rm{q}}}e,{D^{\rm{q}}}{e_{\rm{a}}},{D^{\rm{q}}}{e_{\rm{b}}},{D^{\rm{q}}}{e_{\rm{c}}}} \right]^{\rm{T}}}$

其中,e为同步误差,ea, eb, ec为参数估计误差,P为实对称正定矩阵,选择正定矩阵P=diag(1, 1, 1, 1/λ1, 1/λ2, 1/λ3),λ1, λ2, λ3为特征值.

将同步误差系统(13)、参数估计误差(14)、自适应控制器(16)和参数更新定律(17)代入上式,简化可得

$\begin{array}{l}J = \left[ {e,{e_{\rm{a}}},{e_{\rm{b}}},{e_{\rm{c}}}} \right]\mathit{\boldsymbol{P}}{\left[ {{D^{\rm{q}}}e,{D^{\rm{q}}}{e_{\rm{a}}},{D^{\rm{q}}}{e_{\rm{b}}},{D^{\rm{q}}}{e_{\rm{c}}}} \right]^{\rm{T}}} = \\\;\;\;\;\;{e_1}{D^{\rm{q}}}{e_1} + {e_2}{D^{\rm{q}}}{e_2} + {e_3}{D^{\rm{q}}}{e_3} + \left( {1/{\lambda _1}} \right){e_{\rm{a}}}{D^{\rm{q}}}{e_{\rm{a}}} + \\\;\;\;\;\;\left( {1/{\lambda _2}} \right){e_{\rm{b}}}{D^{\rm{q}}}{e_{\rm{b}}} + \left( {1/{\lambda _3}} \right){e_{\rm{c}}}{D^{\rm{q}}}{e_{\rm{c}}} = \\\;\;\;\;\; - ae_1^2 - 10{e_1}{e_2}{e_3} - 10{e_1}{e_2}{x_3} - 10{e_1}{e_3}{x_2} + \tilde a\left( t \right)e_1^2 - \\\;\;\;\;\;{k_1}e_1^2 + 17{e_1}{e_2}{e_3} + 9{e_1}{e_2}{x_3} + 18{e_1}{e_3}{x_3} + be_2^2 + \\\;\;\;\;\;{e_1}{e_2}{e_3} + {e_1}{e_2}{x_3} + {e_2}{e_3}{x_1} - \tilde b\left( t \right)e_2^2 - {e_2}{e_3}{x_1} + \\\;\;\;\;\;8{e_2}{e_3}{x_1} - e_2^2{e_3} - 2{e_2}{e_3}{x_2} - {k_2}e_2^2 - ce_3^2 - \\\;\;\;\;\;8{e_1}{e_2}{e_3} - 8{e_1}{e_3}{x_2} - 8{e_2}{e_3}{x_1} + e_2^2{e_3} + \\\;\;\;\;\;2{e_2}{x_2}{e_3} - \tilde c\left( t \right)e_3^2 - {k_3}e_3^2 + \\\;\;\;\;\;\left( {1/{\lambda _1}} \right){e_{\rm{a}}}{D^{\rm{q}}}{e_{\rm{a}}} + \left( {1/{\lambda _2}} \right){e_{\rm{b}}}{D^{\rm{q}}}{e_{\rm{b}}} + \left( {1/{\lambda _3}} \right){e_{\rm{c}}}{D^{\rm{q}}}{e_{\rm{c}}} = \\\;\;\;\;\; - {e_{\rm{a}}}e_1^2 - {k_1}e_1^2 + {e_{\rm{b}}}e_2^2 - {k_2}e_2^2 + {e_{\rm{c}}}e_3^2 - {k_3}e_3^2 + \\\;\;\;\;\;\left( {1/{\lambda _1}} \right){e_{\rm{a}}}{D^{\rm{q}}}{e_{\rm{a}}} + \left( {1/{\lambda _2}} \right){e_{\rm{b}}}{D^{\rm{q}}}{e_{\rm{b}}} + \left( {1/{\lambda _3}} \right){e_{\rm{c}}}{D^{\rm{q}}}{e_{\rm{c}}} = \\\;\;\;\;\; - {k_1}e_1^2 - {k_2}e_2^2 - {k_3}e_3^2 - e_a^2 - e_b^2 - e_c^2.\end{array}$

显然,J≤0,据引理1误差系统渐近稳定.

证毕.

采用四阶Runge-Kutta算法,仿真步长为0.001.取q=0.83,a=2,b=1,c=-7.8,k1=k2=k3=20, λ1=λ2=λ3=20.x1(0)=1,x2(0)=2, x3(0)=3作为驱动系统(11)的初始条件;y1(0)=2, y2(0)=3, y3(0)=4作为响应系统(12)的初始条件;$\tilde a(0) = 1, \tilde b(0) = 2, \tilde c(0) = 3$作为参数估计的初始条件.

图 7描述了驱动系统(11)和响应系统(12)的同步误差e1, e2, e3和对系统未知参数的识别.结果表明响应系统与驱动系统在0.2 s内达到同步,在0.2 s内完成对未知参数$\tilde{a}, \tilde{b}, \tilde{c}$的识别.

Fig. 7
图 7 同步误差和未知参数收敛曲线 Fig. 7 Convergence curves of synchronization errors and unknown parameters


4 结论本文设计了一个新型三维分数阶混沌系统,分析了此系统的动力学特性,系统具有复杂的动态行为,能够出现双翼、三翼、四翼等不同类型多翼混沌吸引子共存.对系统进行了Multisim模拟电路设计与仿真,仿真结果与数值分析相符,进一步验证了系统的混沌行为.基于分数阶Lyapunov稳定理论及定理1,设计了一个自适应同步控制器,仿真结果表明了所设计控制器的有效性.


参考文献
[1] KHAN A, KHATTER D, AGRAWAL N. Synchronization of a new fractional order chaotic system[J]. International Journal of Dynamics and Control, 2018, 6(4): 1585. DOI:10.1016/j.chaos.2018.05.015


[2] 孙美美, 胡云安, 韦建明. 一种新型不确定分数阶混沌系统滑模同步控制方式[J]. 电子科技大学学报, 2017, 46(3): 555.
SUN Meimei, HU Yunan, WEI Jianming. A novel sliding mode synchronization method of uncertain fractional-order chaotic systems[J]. Journal of University of Electronic Science and Technology of China, 2017, 46(3): 555. DOI:10.3969/j.ssn.1001-0548.2017.03.013


[3] MUNOZ-PACHECO J M, ZAMBRANO-SERRANO E, VOLOS C, et al. A fractional order chaotic system with a 3D grid of variable attractors[J]. Chaos Solitons & Fractals, 2018, 113: 69. DOI:10.1016/j.chaos.2018.05.015


[4] LAI Qiang, NESTAR T, KENGNE J, et al. Coexisting attractors and circuit implementation of a new 4D chaotic system with two equilibria[J]. Chaos Solitons & Fractals, 2018, 107: 92. DOI:10.1016/j.chaos.2017.12.023


[5] LI Ruiguo, WU Haining. Secure communication on fractional-order chaotic systems via adaptive sliding mode control with teaching-learning-feedback-based optimization[J]. Nonlinear Dynamics, 2018, 1.


[6] SHUKLA M K, SIVA D, MAHAJAN A, et al. A novel scheme of image encryption based on synchronization of fractional order chaotic systems[C]//Proceedings of 2018 International Conference on Intelligent Circuits and Systems (ICICS). Piscataway: IEEE, 2018: 433. DOI: 10.1109/ICICS.2018.00094


[7] PING Z, KE M. An integer-order memristive system with two-to four-scroll chaotic attractors and its fractional-order version with a coexisting chaotic attractor[J]. Complexity, 2018, 2018: 1. DOI:10.1155/2018/4970152


[8] BORAH M, ROY B K. An enhanced multi-wing fractional-order chaotic system with coexisting attractors and switching hybrid synchronisation with its nonautonomous counterpart[J]. Chaos, Solitons & Fractals, 2017, 102: 372. DOI:10.1016/j.chaos.2017.03.055


[9] 鲜永菊, 夏诚, 徐昌彪. 新混沌系统及其分数阶系统的自适应同步控制[J]. 控制理论与应用, 2018, 35(6): 159.
XIAN Yongju, XIA Cheng, XU Changbiao. New chaotic system and the adaptive synchronization control of its fractional order system[J]. Control Theory & Applications, 2018, 35(6): 159. DOI:10.7641/CTA.2017.70567


[10] ZHANG Chaoxia, YU Simin. Generation of multi-wing chaotic attractor in fractional order system[J]. Chaos, Solitons & Fractals, 2011, 44(10): 845. DOI:10.1016/j.chaos.2011.06.017


[11] TANG Zilong, ZHANG Chaoxia, YU Simin. Design and circuit implementation of fractional-order multiwing chaotic attractors[J]. International Journal of Bifurcation and Chaos, 2012, 22(11): 1250269. DOI:10.1142/S0218127412502690


[12] LYU Jinhu, CHEN Guangrong, CHENG Daizhan. A new chaotic system and beyond:The generalized Lorenz-like system[J]. International Journal of Bifurcation and Chaos, 2014, 14(5): 1507. DOI:10.1142/S021812740401014X


[13] DIETHELM K, FORD N J, FREED A D. A predictor-corrector approach for the numerical solution of fractional differential equations[J]. Nonlinear Dynamics, 2002, 29(1-4): 3.


[14] BASHKIRTSEVA I A, ZAITSEVA S. The phenomenon of stochastic excitability in the enzymatic reaction model[J]. Institute of Natural Science & Mathematics, 2018, 3. DOI:10.20537/vm180101


[15] SAHA S, MAITI A, SAMATA G P. A Michaelis-Menten predator-prey model with strong Allee effect and disease in prey incorporating prey refuge[J]. International Journal of Bifurcation and Chaos, 2018, 28(06): 1850073. DOI:10.1142/S0218127418500736


[16] 周平, 危丽佳, 程雪峰. 一个分数阶混沌系统及其EWB仿真结果[J]. 西南大学学报(自然科学版), 2009, 31(1): 69.
ZHOU Ping, WEI Lijia, CHENG Xuefeng. A fractional-order chaotic system and its simulation result electronic work-bench (EW-B)[J]. Journal of Southwest University (Natural Science), 2009, 31(1): 69. DOI:10.13718/j.cnki.xdzk.2009.01.014


[17] OGATA K. Modern control engineering[M]. New Jersey: Prentice Hall, 2002.


[18] TAVAZOEI M S, HAERI M. A necessary condition for double scroll attractor existence in fractional-order systems[J]. Physics Letters A, 2007, 367(1-2): 102. DOI:10.1016/j.physleta.2007.05.081


[19] SOHAIL A, MAQBOOL K, ELLAHI R. Stability analysis for fractional-order partial differential equations by means of space spectral time Adams-Bashforth Moulton method[J]. Numerical Methods for Partial Differential Equations, 2018, 34(1): 19. DOI:10.1002/num.22171


[20] 邵书义.新型混沌系统的动力学分析、同步控制及应用[D].南京: 南京师范大学, 2014
SHAO Shuyi. Dynamic analysis, synchronous control and application of new chaotic system[D]. Nanjing: Nanjing Normal University, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10319-1014343213.htm


[21] LUO Shaohua, LI Shaobo, TAJADDODIANFAR F. Chaos and adaptive control of the fractional-order magnetic-field electromechanical transducer[J]. International Journal of Bifurcation and Chaos, 2017, 27(13): 1750203. DOI:10.1142/S0218127417502030


[22] LUO Shaohua, LI Shaobo, TAJADDODIANFAR F, et al. Adaptive synchronization of the fractional-order chaotic arch micro-electro-mechanical system via Chebyshev neural network[J]. IEEE Sensors Journal, 2018, 18(9): 3524. DOI:10.1109/JSEN.2018.2812859


[23] LI Ruiguo, WU Haining. Adaptive synchronization control based on QPSO algorithm with interval estimation for fractional-order chaotic systems and its application in secret communication[J]. Nonlinear Dynamics, 2018, 1.


[24] HUANG Yu, WANG Dongfeng, ZHANG Jinying, et al. Controlling and synchronizing a fractional-order chaotic system using stability theory of a time-varying fractional-order system[J]. PloS One, 2018, 13(3). DOI:10.1371/journal.pone.0194112



相关话题/系统 分数

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 液压轮毂马达辅助驱动系统控制策略实车验证
    液压轮毂马达辅助驱动系统控制策略实车验证曾小华,崔臣,张轩铭,宋大凤,李立鑫(汽车仿真与控制国家重点实验室(吉林大学),长春130025)摘要:为充分提升重型牵引车辆通过不良路面的能力,对国内某款重型牵引车在传统结构的基础上加装了前轴液压轮毂马达辅助驱动系统,并针对该混合动力系统,开发了工程化的控制 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 铁锰复合阴极MFC-EF耦合系统产电及降解RhB效能
    铁锰复合阴极MFC-EF耦合系统产电及降解RhB效能史珂1,2,赵庆良1,2,王维业2,王琨1,2(1.城市水资源与水环境国家重点实验室(哈尔滨工业大学),哈尔滨150090;2.哈尔滨工业大学环境学院,哈尔滨150090)摘要:为实现微生物燃料电池(MFC)微电的原位利用,结合电芬顿(EF)技术的 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 偏最小二乘法在系统故障诊断中的应用
    偏最小二乘法在系统故障诊断中的应用梁北辰,戴景民(哈尔滨工业大学仪器科学与工程学院,哈尔滨150001)摘要:为研究结合关键性能指标(KPI)的偏最小二乘法(PLS)在系统故障诊断中的作用和应用,通过活性污泥数学模型(BSM1),结合偏最小二乘算法,建立了可以模拟污水处理的仿真模型.通过收集系统在正 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 椭圆轨道绳系卫星系统释放的类反步法控制
    椭圆轨道绳系卫星系统释放的类反步法控制仲小清1,金雪松2,王敏1,李晓磊3,孙光辉3(1.中国空间技术研究院通信卫星事业部,北京100094;2.中国航天科技集团宇航部,北京100048;3.哈尔滨工业大学航天学院,哈尔滨150001)摘要:为解决椭圆轨道上绳系卫星系统的稳定和快速释放问题,在传统反 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 船用极区格网惯性导航系统综合校正方法
    船用极区格网惯性导航系统综合校正方法方涛,黄卫权,王宗义(哈尔滨工程大学自动化学院,哈尔滨150001)摘要:为抑制适用于极区的船用格网惯性导航系统随时间积累的导航误差,提出一种基于格网坐标系的综合校正方法对陀螺常值漂移进行估计和补偿.该方法首先基于格网坐标系推导了P方程,该方程建立起位置误差、格网 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 菱形翼布局无人机自适应分数阶滑模姿态控制
    菱形翼布局无人机自适应分数阶滑模姿态控制刘海波1,王和平1,2,孙俊磊1(1.西北工业大学航空学院,西安710072;2.西北工业大学深圳研究院,广东深圳518057)摘要:为研究存在复合干扰的非常规布局菱形翼长航时侦察无人机姿态控制问题,针对系统存在强耦合、非线性、多输入多输出等特点,结合滑模变结 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 星载固态存储系统自适应闪存转换层设计
    星载固态存储系统自适应闪存转换层设计张伟东1,2,董振兴1,朱岩1,安军社1(1.中国科学院复杂航天系统电子信息技术重点实验室(中国科学院国家空间科学中心),北京100190;2.中国科学院大学,北京100190)摘要:传统星载存储系统闪存转换层(FlashTranslationLayer,FTL) ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 航姿参考系统的改进杆臂效应补偿方法
    航姿参考系统的改进杆臂效应补偿方法谭强俊,程永生,唐彬,刘显学,周浩,李寅鑫(中国工程物理研究院电子工程研究所,四川绵阳621000)摘要:振动环境下叠加在航姿参考系统(AHRS)加速度计敏感轴上的杆臂效应会对加速度计测姿精度造成极大的影响.本文对杆臂效应处理算法中的低通滤波法和力学补偿法进行分析. ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 一种四进制Duffing混沌数字通信系统
    一种四进制Duffing混沌数字通信系统李雅楠,付永庆(哈尔滨工程大学信息与通信工程学院,哈尔滨150001)摘要:为提高多进制混沌通信系统的性能,解决Duffing混沌系统受混沌同步技术限制的问题,本文设计了基于Duffing振子的四进制混沌数字通信系统.利用二进制Duffing混沌调制方法与正交 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 可调磁通电机系统及其关键技术发展
    可调磁通电机系统及其关键技术发展郑萍,王明峤,乔光远,刘法亮,张书宽(哈尔滨工业大学电气工程及自动化学院,哈尔滨150080)摘要:针对传统汽车工业发展面临的能源危机和环境污染等问题,为满足新能源汽车驱动电机在宽调速范围内高效率运行的需求,对可调磁通电机(variable-fluxmachine,V ...
    本站小编 哈尔滨工业大学 2020-12-05