李 銀，呂顯瑞，趙 昕，姜 博

（1.韶關(guān)學(xué)院數(shù)學(xué)與信息科學(xué)學(xué)院，廣東韶關(guān) 512005；2.中山大學(xué)數(shù)學(xué)與計(jì)算機(jī)學(xué)院，廣州 510275；3.吉林大學(xué)數(shù)學(xué)學(xué)院，長(zhǎng)春 130012；4.吉林農(nóng)業(yè)大學(xué)信息技術(shù)學(xué)院，長(zhǎng)春 130118；5.新加坡國(guó)立大學(xué)理學(xué)院，新加坡119077，新加坡）

非線性Lorenz－Like系統(tǒng)的定性分析及控制

李 銀1，2，呂顯瑞3，趙 昕4，姜 博5

（1.韶關(guān)學(xué)院數(shù)學(xué)與信息科學(xué)學(xué)院，廣東韶關(guān) 512005；2.中山大學(xué)數(shù)學(xué)與計(jì)算機(jī)學(xué)院，廣州 510275；3.吉林大學(xué)數(shù)學(xué)學(xué)院，長(zhǎng)春 130012；4.吉林農(nóng)業(yè)大學(xué)信息技術(shù)學(xué)院，長(zhǎng)春 130118；5.新加坡國(guó)立大學(xué)理學(xué)院，新加坡119077，新加坡）

利用定性分析和優(yōu)化控制方法，實(shí)現(xiàn)了Lorenz－Like模型的Hopf分支和優(yōu)化控制，得到了模型的動(dòng)力學(xué)演化機(jī)理和分支條件及變量同步控制器和誤差系統(tǒng)的全局漸近穩(wěn)定條件.

Lorenz－Like系統(tǒng)；定性分析；Hopf分支；變量廣義控制

Lorenz－Like系統(tǒng)是Lorenz混沌系統(tǒng)的一種延拓，混沌系統(tǒng)在保密通信、激光物理、化學(xué)反應(yīng)和生物醫(yī)學(xué)等領(lǐng)域應(yīng)用廣泛［1－2］.本文針對(duì)Lorenz－Like模型［3］的動(dòng)力特性和狀態(tài)變量給出定性分析和優(yōu)化控制.首先，給出Lorenz－Like模型的吸引子、平衡點(diǎn)和耗散性等動(dòng)力學(xué)特性；其次，研究模型的Hopf分支，給出判斷分支的表達(dá)式和判斷依據(jù)；最后，利用優(yōu)化控制方法給出模型部分變量快速同步的控制器.數(shù)值模擬表明，其誤差系統(tǒng)是全局漸近穩(wěn)定的，證實(shí)了本文方法的有效性與可行性.

1 Lorenz－Like系統(tǒng)模型及定性分析

考慮如下Lorenz－Like系統(tǒng)模型［3－5］：其中（x，y，z）T∈R3為系統(tǒng)的狀態(tài)變量.當(dāng)系統(tǒng)參數(shù)p，q，k，h，s取正值時(shí)，文獻(xiàn)［3］證明了系統(tǒng)的最大Lyapunov指數(shù)是正的；當(dāng)p＝10，k＝1，h＝4，s＝2.5，q＝40時(shí)，系統(tǒng)吸引子如圖1（A）所示；當(dāng)p＝a，k＝a，h＝1，s＝c，q＝ab時(shí)，式（1）即為文獻(xiàn)［6－7］的系統(tǒng)；而當(dāng)a＝5，b＝4，c＝2時(shí)，系統(tǒng)吸引子如圖1（B）所示；當(dāng)p＝a，k＝a，h＝1，s＝b，q＝c－a時(shí)，系統(tǒng)（1）即為文獻(xiàn)［8］的系統(tǒng)；當(dāng)a＝2.1，b＝0.6，c＝30時(shí)，系統(tǒng)具有如圖1（C）所示的混沌吸引子；當(dāng)p＝a，k＝1，h＝1，s＝c，q＝b時(shí)；系統(tǒng)（1）即為文獻(xiàn)［9－10］的系統(tǒng)；當(dāng)a＝10，b＝16，c＝8／3時(shí)，系統(tǒng)具有如圖1（D）所示的混沌吸引子.

圖1 混沌系統(tǒng)的動(dòng)力學(xué)演化Fig.1 Evolutionary dynamics of the chaotic system

2 混沌系統(tǒng)（1）的Hopf分支

下面考慮混沌系統(tǒng)（1）的Hopf分支情況.設(shè)Cn是復(fù)數(shù)域C上的線性空間.數(shù)量積〈u，v〉中u，v∈Cn滿足如下條件：

3 部分變量廣義控制的應(yīng)用

基于文獻(xiàn)［13］的方法，考慮兩個(gè)混沌系統(tǒng)：一個(gè)狀態(tài)變量為x∈Rp，另一個(gè)狀態(tài)變量為y∈Rq，假設(shè)x m∈Rs和y s∈Rs為其狀態(tài)變量的子集.

由式（26）可知，Ω＜p有很多類型的解，系統(tǒng)（24）特征根有負(fù)實(shí)部，

由于參數(shù)和初值的選擇為任意，因此不失一般性，可選擇參數(shù)如下：p＝10，k＝1，h＝4，s＝2.5，q＝40，a＝2.1，b＝0.6，c＝30，θ（z2（t－τ））＝3z2（t－τ），初值x1（0）＝2，y1＝0，z1x2（0）＝0，y2＝0.經(jīng)過簡(jiǎn)單計(jì)算誤差可知，e1（0）＝－13，e2（0）＝－136.73，e3（0）＝－252.94.圖2（A）～（C）分別為部分同步τ＝0，0.3，－0.3的數(shù)值仿真誤差曲線.

圖2 誤差同步演化曲線Fig.2 Error synchronizing evolution curves

綜上，本文利用控制理論和優(yōu)化控制方法，探討并定性分析了Lorenz－Like模型的動(dòng)力學(xué)特性及Hopf分支與部分變量的廣義控制，實(shí)現(xiàn)了誤差變量的快速同步控制.數(shù)值模擬結(jié)果表明：誤差系統(tǒng)是全局漸近穩(wěn)定的；本文方法簡(jiǎn)單易學(xué)，且實(shí)現(xiàn)同步控制的時(shí)間短，是實(shí)現(xiàn)系統(tǒng)模型控制的有效方法.

［1］ 劉式適，劉式達(dá).物理學(xué)中的非線性方程［M］.北京：北京大學(xué)出版社，2012：341－345.（LIU Shishi，LIU Shida.Nonlinear Equations in Physics［M］.Beijing：Peking University Press，2012：341－345.）

［2］ 王興元.混沌系統(tǒng)的同步及在保密通信中的應(yīng)用［M］.北京：科學(xué)出版社，2012：356－385.（WANG Xingyuan.Synchronization and Application of Chaos the Systems in Secure Communication［M］.Beijing：Science Press，2012：356－385.）

［3］ MU Chunlai，ZHANG Fuchen，SHU Yonglu，et al.On the Boundedness of Solutions to the Lorenz－Like Family of Chaotic Systems［J］.Nonlinear Dynamics，2012，67（2）：987－996.

［4］ Nusse H E，Yorke J A.Dynamics：Numerical Explorations［M］.New York：Springer－Verlag，1998.

［5］ LIU Yongjian，YANG Qigui.Dynamics of a New Lorenz－Like Chaotic System［J］.Nonlinear Analysis：Real World Applications，2010，11（4）：2563－2572.

［6］ JIANG Bo，HAN Xiujing，BI Qinsheng.Hopf Bifurcation Analysis in theTSystem［J］.Nonlinear Analysis：Real World Applications，2010，11（1）：522－527.

［7］ Tigan G，Opris D.Analysis of a 3D Chaotic System［J］.Chaos，Solitons ＆Fractals，2008，36（5）：1315－1319.

［8］ Tigan G，Constantinescu D.Heteroclinic Orbits in theTand the LüSystems［J］.Chaos，Solitons ＆Fractals，2009，42（1）：20－23.

［9］ LI Xianfeng，CHU Yandong，ZHANG Jiangang，et al.Nonlinear Dynamics and Circuit Implementation for a New Lorenz－Like Attractor［J］.Chaos，Solitons ＆Fractals，2009，41（5）：2360－2370.

［10］ LI Yin，CHEN Yong，LI Biao.Adaptive Control and Function Projective Synchronization in 2D Discrete－Time Chaotic Systems［J］.Communications in Theoretical Physics，2009，51（2）：270－278.

［11］ YAN Zhenya.Hopf Bifurcation in the Lorenz－Type Chaotic System［J］.Chaos，Solitons ＆Fractals，2007，31（5）：1135－1142.

［12］ Yuri A K.Elements of Applied Bifurcation Theory［M］.New York：Springer－Verlag，1995.

［13］ LI Yin，LI Biao，CHEN Yong.Anticipated Function Synchronization with Unknown Parameters of Discrete－Time Chaotic Systems［J］.International Journal of Modern Physics C，2009，20（4）：597－608.

Qualitative Analysis and Optimal Control of a Nonlinear Lorenz－Like Model

LI Yin1，2，LüXianrui3，ZHAO Xin4，JIANG Bo5
（1.SchoolofMathematicsandInformationScience，ShaoguanUniversity，Shaoguan512005，GuangdongProvince，China；2.SchoolofMathematicsandComputationalScience，SunYat－SenUniversity，Guangzhou510275，China；3.CollegeofMathematics，JilinUniversity，Changchun130012，China；4.CollegeofInformationTechnology，JilinAgriculturalUniversity，Changchun130118，China；5.FacultyofScience，NationalUniversityofSingapore，Singapore119077，Singapore）

We realized the Hopf bifurcation and optimal control of Lorenz－Like models via the qualitative analysis and optimal control methods，and obtained the dynamical evolution mechanism，the conditions of bifurcation；variable synchronous controller and the stable conditions of the error system.

Lorenz－Like model；qualitative analysis；Hopf bifurcation；variables generalized control

O175.12

A

1671－5489（2014）03－0429－06

10.13413／j.cnki.jdxblxb.2014.03.05

2013－09－12.

李 銀（1980—），男，漢族，博士研究生，講師，從事非線性動(dòng)力學(xué)的研究，E－mail：liyin2009521＠163.com.通信作者：趙 昕（1974—），男，漢族，博士，副教授，從事微分方程與動(dòng)力系統(tǒng)的研究，E－mail：jlndzx＠sina.com.

國(guó)家自然科學(xué)基金（批準(zhǔn)號(hào)：10871214）、吉林省自然科學(xué)基金（批準(zhǔn)號(hào)：201215184）和韶關(guān)市科技計(jì)劃項(xiàng)目（批準(zhǔn)號(hào)：313140546）.

趙立芹）