计算机科学2007V01.34N0.12
分数阶Lii系统中的混沌及其控制*)
武相军1王兴元2
(河南大学计算中心
开封475004)1
(大连理工大学电信学院计算机系大连116024)2
摘要研究了分数阶Ln系统的混沌动力学行为、数值模拟证明分数阶Lil系统存在混沌,并且得出分数阶Ln系统能产生混沌吸引子的最低阶数为2.5阶。利用线性反馈控制法研究了分数阶Lil混沌系统的混沌控制问题,得出了受控分数阶LIi混沌系统的混沌轨道达到不稳定平衡点时的条件,数值模拟进一步验证了该方法的有效性.关键词分数阶Ln系统,线性反馈控制法,混沌控制
ChaosintheFractionalOrderLfiSystemanditsControl
WUXiang-JunIWANGXing-Yuan2
(ComputingCenter,HenanUniversity,Kaifeng475001)1
(Schoolof日ectronic&InformationEngineering,DalianUniversityofTechnology。Dali珊ai16024)z
Abstract
ThechaoticbehaviorsinthefractionalorderLUsystemisstudiedinthispaper.Chaosexistsinthefractional
to
orderLUsystemwithorderles¥than3verifiedbynumericalsimuhtions.Thelowe.storderwefoundhavechaos
in
thissystemis2.5.TheproblemofcontrolchaoticbehaviorofthefractionalorderLnchaoticsystemusingthelinearfeedbackcontrolmethodisaddressed.Theconditionssuppressingchaostounstableequilibriumpointmericalsimulationsshowtheeffectivenessofthelinearfeedbackmethod.
Keywords
are
derived.Nu.
FractionalorderLnsystem。Linearfeedbackcontrolmethod,Chaoscontr01
分数阶微积分的理论研究已有300多年的历史,但是分数阶微积分理论由于长期没有实际应用背景而发展缓慢。近年来,国内外学者研究发现,在黏弹性系统、电介质偏振现象、电极一电解偏振现象、电磁波中存在分数阶动力系统,从而使得分数阶微积分理论应用到物理和工程领域成为一个热点研究课题E],23。最近,分数阶混沌系统引起人们广泛的兴趣和深入的研究。在Lorenz混沌系统,Duffing混沌系统、Chua’s混沌电路、ROssler混沌和超混沌系统以及临界混沌系统[刀中,通过数值仿真,发现当系统的阶数为分数时仍出现混沌状态E3-q,且更能反映系统呈现的工程物理现象,从而促进了分数阶混沌的研究以及分数阶微积分理论的发展。近来,人们
它的分数阶微分就不一定具有这样的意义了,因此函数,(£)的初始值常常难以测度和给定。
下面我们看另一种分数阶微分定义——caputo定义的分数阶微分‘2川],其数学表达式如下,
业2:,一业.
de
。
舻’
(3)
式中,玎一1≤<起,口为分数,挖为正整数,J,是卢(胗0)阶硒一
eamrm-Liouville积分算子,定义如下
珊(£)一南J
值容易测度和确定。
1
6尚dr
∥一、
(4)
尝试研究分数阶混沌系统的控制与同步问题。文E83研究了
分数阶Chen系统的混沌特性和分数阶Chen系统的混沌控
使用Caputo定义的分数阶微分,对于初始值,(o),/(o)’..・,,”D(o),由于它们具有各自的含义,因此这些初始
如果使用Riemann-Liouville定义的分数阶微分,需要预先给定初始条件,并且要求这些初始条件具有相同的意义或性质;而Caputo定义的分数阶微分允许非同质的初始条件,在初始条件完全相同的情况下,RiemanrrLiouviUe定义的分
制;文[93讨论了分数阶混沌系统的混沌同步问题。
分数阶微分有多种定义[1’2’1引,在以上文献的分数阶微分研究中,应用的都是Riemann-Liouville定义的分数阶微分,其数学表达式如下:
生世:!丑业f£.丛互2。drr(起一口)de'’o(--r)。一。+l’
d,
数阶微分和Caputo定义的分数阶微分是等价的[1.z].基于
f1、
、17
以上的分析,作者使用Caputo定义的分数阶微分研究了分数阶Lti系统的混沌动力学行为。给出了数值仿真结果,并且利用线性反馈控制法研究了分数阶Ln系统的混沌控制,数值
式中,r(・)是伽马函数,玎一1≤口<他,口为分数,露为正整数.
如果函数八£)的初始值均为零,则(1)式的拉普拉斯变换表
达式可表示为
模拟进一步验证了该控制方法的有效性。
(2)
L{掣}=川彻)
1
分数阶微分及其近似计算
2002年,吕金虎等人发现了介于Lorenz系统m1和Chert
使用Riemann-Liouville定义的分数阶微分.必须精确知道分数阶微分的未知函数在初始时间t=0时的初始值m】。在实际应用中,函数,(£)可能有某种特定的物理意义,但是
系统‘1钉之间的一个新系统,即Lu系统‘・‘],它由如下的三维自治方程组来描述
*)本文得到国家自然科学基金(编号:60573172)和河南省自然科学基金(编号。0511011400和0611054700)资助.・204・
万方数据
fz&=n(y—z)其中黛、Y、=为系统的状恣变量,a、b和f为系统的控制参数。
.{y&=--xz+cy,
(5)
当参数a----36、b=3和c=20时,系统(5)进入混沌状态,图1
【z&=xy一搬
所承为如系统的混沌吸弓l子。
Z
J
),
匿1糙系统酶混淹吸萼|子
f鲁=口(y--x)
斟t毛+—‰口(豫t一珠・)十{辫一鞘◆锣.氛磊h.=elF两妻口l幽一laC。yj一而)国
跚2殉÷志e一粕翻+嘲汁㈣・茅石而岩∥1幽一l一而)
悟--'ry
O’Z
龈t=翱+志(臻,蠕+t一62‰t)+
氚瓦he+z2)二∥2.m+1(一而z,+cM)
根据文[1s~17],我识霹浚褥戮系统<6>鹊溪溺一修正算
法,该方法是Adams-Bashforth-Mouhon方法的一般化。微亍哥差等丽查舔一种l(乃弱一酝i)赢瓦硒暑∥3一种l
t乃弱一%’分方程
式巾,
。等=“∽),o≤≤T
.
f磁+,=动+亍巧1再差狰小抖t口(∞一而)
z疆’(O);森酚,k=O,1,2。…,n--1
{难t2骛+i孓瑟j羞》J。一t(一妨≈+cyj),
等价于如下的Volterra积分方程㈣
z=勖)l!}+忐Js磐抟如
㈣
【锑+t=知+志盏p西一t(而"一慨)
dti,j。#+l。
瑕设矗2:y炎,气=nh(n=0,1,2,…,N),弱式(7)霹用蘩
f嘏十‘~(携一壤)≤露+1)赣,
j--O下式子表示;
_{(n--j+2)8i+1+(n--j)8i+1--2(n--j+1)e,+1.
1≤0≤孵,。
魏(“-)一磊毋,锊+穴轰丽““t.霸(“一))
11.
j=n+1
编∑鸯,一lf(ti,磊<务))
届小时l=警((靠一I『+1)dI一妇一』)吨)。
O≤j≤露,i=1,2,3.
式中.
2分数阶L矗混沌系统的数值模拟
aj.一12
fn"+1--(n--a)(n+1)。+1。
j=O
在以下的数值模拟中,分数阶系统(6)的参数仍然取口一{(撵一歹+2)f}1+(n--j)。州--2(n--j+1).+1,l≤,《洳,
36、b=3翔c=20。佟者分别利用第2节分数阶系统的离{敦
【1。j=n+l
纯方法耱Wolf求系统最大Lyapunov擐羧瞻方法疆《,褒究了
当初始值分别取z(O)-----1、,(o)=z和z(O)=3,受(o<绕≤l,i霸(“)一窟t--O舻智+志磐州f(ti州枷,
;l,2,3)取不同值时,分数阶Ln系统的幼力学特性。
岛,抖z一警((抖~歹+1).一<蔸一J)。)。
(1)01+岛+岛>/-2.5
銎巍=o。9、岛=o.94和岛=o。96嚣寸,分数阶系统(g)的
误差德诗麦£一Maxlz(tj)一蕊(务)|=0(妒)0一O,1,…,耀黧魏烫2瑟示,其最大跏辨nov揍数必筑5468,霉凳分数
^D,其中p=Min(2,】+口).
阶系统(6)确实处予混沌状态。当蘸一O.9、受一o.8和秘一使用上述方法,可将分数阶系统(6)离散化为如下表达0.8时,分数阶系统(6)的相图如图3所承,其最大Lyapunov式:
指数为0.3954,可见分数阶系统(6)是混沌系统。当融一O.9、瑰---0.7和岛=o.9时,分数阶系统(6)的榴匿如蓬4所零,
万
方数据・205・
其最大Lyapunov指数为0.3216,可见系统(6)是混沌系统。同理,当0<岛≤1(i=1,2,3)且a+岛+岛≥2.5时,系统同
样处于混沌状态。
Z
),
图2岛=O.9,&=O.94和岛=O.96时分数阶Lil系统的相图
),
图3岛=o.9,晓=O.8和岛=0.8时分数阶Lll系统的相图
),
图4岛----09,岛=O.7和岛----0.9时分数阶Lil系统的相图
),
图5毋=O.9,如----0.7和侥=0.7时分数阶Lil系统的相图
当兜=0.9、&=0.7和岛=O.7时,系统(6)的相图如图5所示,系统的运动轨道趋向于一个极限环(不稳定周期轨道),
・206.
可见系统(6)不是混沌的。同理,当O<晚≤1(i=1,2,3)且岛
(下转第210页)
万方数据
VM更为显著的识别效果。
参考文献
l
themargindistribu—
tion.In:Proceedingsofthe12thConferenceonComputationalTaylor
on
4
HandwritingRecog・
nition。2004.167~172
GoughJ.ChothiaCSUPERR~MILY:HMMSrepresentingall
er
identificationandverification.Frontiersin
JS,CristianiniN.Furtherresults
proteinsofknownstructure.900Psequencesearches,alignmentsandgenomeassignments.NueleicAcidsResearch,2002,30(1)l
268~272
56
BaumLE.SellGR.Growthfunctionsfor
LearningTheory.1999
23
transformationz
Off
J&Optimizingclassifiersforimbalanced
sets.NIPS,1999,11:253~259training
SchiapbachA.BunkeH.UsingHMMbasedRecognizersforwrit’
KarakoulasGJ.Taylor
manifolds.PacJMath。1968,(27)2:21l~227
BakerJkThedragonsystem-Anoverview.IEEEAcoustSpeechSignalProcessing,1975,AASSP-23(1):24~29
(上接第206页)
+如+03<2.5时,系统同样也不是混沌的。
综合(1)、(2),可知分数阶Ln系统能产生混沌吸引子的最低阶数为2.5阶。
(6)的混沌运动轨道镇定到不稳定平衡点E上。设计线性反馈控制器为
r让15五l(z一;)
<砌=乜(y--y),
(10)
3分数阶Ln混沌系统的控制
在本节中,我们研究将分数阶Ln混沌系统镇定到它的不稳定平衡点。取研一岛=岛=O.9,显然,分数阶Ln系统处于混沌状态,并且E(0,O,O)是它的一个不稳定平衡点。下面我们利用线性反馈控制法将分数阶Ln混沌系统的运动轨道镇定到不稳定平衡点。
受控分数阶Ln系统为
f掸z
,
、
【均=五,(2一;)
式(10)中(至,;,三)代表系统(6)的不稳定平衡点E,愚1、岛
和岛是反馈增益。
受控系统(9)在平衡点处的Jacobi矩阵为
r一口一愚1
a
00
产1。0
其特征方程为
c一是2
O
一6一岛
l面2口(Y一工)一“l
(A+口+量1)(A—c+忌2)(A+b+k3)=O
显然,特征方程的特征根为
(9)
A1=一n一愚l’A2=c一愚z,A3=一6一岛
{参一xz+cx飞。
I水z
.
当特征方程的实特征根均为负值时,受控系统(9)渐近稳定地趋于平衡点E,可求得岛>一36、量2>20和五3>--3。
【面。xy—bz一均
这里0=29,t/./(i=1,2,3),是反馈外部输入控制,可使系统
t/s
图6利用线性反馈法镇定分数阶Ln系统到平衡点时工(£)、y(f)和2(£)的变化曲线
图6给出了当选取受控系统(9)的初始点为:z(0)=2、Y(O)一1和z(O)=一1,取kl—O、量2----25和愚3=O时,受控系统(9)镇定到平衡点E(0,0,o)上的结果。由图6可见:当t分别接近5.8s、7.8s和13.8s时,z(f)、y(£)和z(t)分别稳定到了零点,即受控系统(9)被镇定到平衡点E(0,0,o)上。
结论本文研究了分数阶Ln系统的混沌动力学特性,数值模拟证明分数阶Ln系统确实存在混沌,并且得出分数阶Ln系统能产生混沌吸引子的最低阶数为2.5阶。作者利
67891345
GrigorenkoI,GrigorenkoEPhys.Rev.Lett.9134101,2003GaoX,YuJBChaosSolitonsFract..241097,2005HartlyTT。LorenzoC485,1995
LiCG.ChertGRPhysicaA.34155,2004
F,QamlnerHk
IE既Trans.CAS
I,42
王发强,刘崇新.物理学报.553922,2006
LiCG.ChenGRChaosSolitonzFract.。22549,2004
LiCG。LiaoXF。Yu
JBPhys.Rev.,E6867203,2003
SamkoSG.KlihasAA,MafichevOLFractionalintegralsand
derivativesltheoryBreach),1993
andapplications(Amsterdam:Gordanand
用线性反馈控制器成功地将分数阶n混沌系统镇定到平衡
点。
1112
Caputo乩ILGeophys.J.IL
Astron.S[k,13529,1967
LorenzEN.J.Atmos.Sci.20130,1993
参考文献
1
PodlubnyLFractionaldemiePress),19992
HiferILApplicationsof
fractional
calculus
in
physics(New
differential
131415
16
ChenGR.UetaT-Int.J.Bifur.Chaos,91465,1993
Lti
J,ChenGILInt.J.Bifur.Chaos,12659,2002
equations(NewYork:Aca—
Diethelml(.Elec.Trans.Numer.AnaL,5l,1997DiethelmK,FordNDiethelmK,F0rdN
J.J.Math.AnaL
1718
J.Nonlinear脚L,29
AppL,265229・2002
3,2002
Jersey:WbrldScientific),2001
WolfA,SwinneyJB.SwinneyHLPhysicaD.16285,1985
・210.
万方数据
分数阶Lü系统中的混沌及其控制
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
武相军, 王兴元, WU Xiang-Jun, WANG Xing-Yuan
武相军,WU Xiang-Jun(河南大学计算中心,开封,475004), 王兴元,WANG Xing-Yuan(大连理工大学电信学院计算机系,大连,116024)计算机科学
COMPUTER SCIENCE2007,34(12)4次
参考文献(18条)
1. Podlubny I Fractional differential equations 1999
2. Hifer R Applications of fractional calculus in physics 20013. Grigorenko I;Grigorenko E 查看详情 20034. Gao X;Yu J B 查看详情 2005
5. Hartly T T;Lorenzo C F;Qammer H K 查看详情 19956. Li C G;Chen G R 查看详情 2004
7. 王发强;刘崇新 分数阶临界混沌系统及电路实验的研究[期刊论文]-物理学报 2006(8)8. Li C G;Chen G R 查看详情 20049. Li C G;Liao X F;Yu J B 查看详情 2003
10. Samko S G;Klibas A A;Mariehev O I Fractional integrals and derivatives:theory and applications1993
11. Caputo M.II;Geophys.J.R 查看详情 196712. Lorenz E N 查看详情 1993
13. Chen G R;Ueta T 查看详情[外文期刊] 199314. Lu J;Chen G R 查看详情 200215. Diethelm K 查看详情 1997
16. Diethelm K;Ford N J 查看详情 200217. Diethelm K;Ford N J 查看详情 2002
18. Wolf A;Swinney J B;Swinney H L 查看详情 1985
本文读者也读过(10条)
1. 张济仕. 史蕊. ZHANG Ji-shi. SHI Rui 分数阶Rikitake系统中的混沌及其控制[期刊论文]-河南大学学报(自然科学版)2009,39(1)
2. 乔宗敏. 金永容. QIAO Zong-min. JIN Yong-rong 分数阶Lorenz系统的混沌控制方法[期刊论文]-安徽大学学报(自然科学版)2010,34(6)
3. 周平. 张年英. ZHOU Ping. ZHANG Nian-ying 利用线性反馈控制分数阶Rǒssler混沌系统到达任意稳定状态[期刊论文]-重庆邮电大学学报(自然科学版)2007,19(6)
4. 孙克辉. 任健. 丘水生 分数阶Chen混沌系统的PC同步控制研究[会议论文]-2007
5. 周亚非. 王中华. ZHOU Ya-fei. WANG Zhong-hua 分数阶混沌激光器系统的同步[期刊论文]-半导体光电2008,29(5)
6. 陈向荣. 刘崇新. 王发强. 李永勋. Chen Xiang-Rong. Liu Chong-Xin. Wang Fa-Qiang. Li Yong-Xun 分数阶Liu混沌系统及其电路实验的研究与控制[期刊论文]-物理学报2008,57(3)
7. 陈保颖. 张家军. 苑占江. Chen Baoying. Zhang Jiajun. Yuan Zhanjiang 分数阶Rucklidge混沌系统的同步研究[期刊论文]-动力学与控制学报2010,08(3)
8. 张俊. ZHANG Jun 带时延的分数阶神经网络系统的混沌与同步研究[期刊论文]-电脑知识与技术2010,06(36)9. 武向军 承载技术的革命——PTN[期刊论文]-通信世界B2009(40)
10. 吴耿. 丰建文. 张维强. Wu Geng. Feng Jianwen. Zhang Weiqiang 分数阶Genesio系统的混沌同步[期刊论文]-三峡大学学报(自然科学版)2008,30(4)
引证文献(4条)
1. 张济仕. 史蕊 分数阶Rikitake系统中的混沌及其控制[期刊论文]-河南大学学报(自然科学版) 2009(1)2. 渠慎明. 侯松鹂 R(O)ssler系统中的混沌控制[期刊论文]-光盘技术 2009(4)
3. 冯伟 基于FPGA的超混沌PN序列产生与性能分析[期刊论文]-杭州电子科技大学学报 2011(1)
4. 董俊. 张广军. 姚宏. 王相波 分数阶动力学模型的混沌特性及投影同步控制[期刊论文]-计算机仿真 2012(12)
引用本文格式:武相军. 王兴元. WU Xiang-Jun. WANG Xing-Yuan 分数阶Lü系统中的混沌及其控制[期刊论文]-计算机科学 2007(12)