陸安山
(欽州學(xué)院物理與電子工程學(xué)院,廣西 欽州535000)
混沌是一種具有初值敏感性、遍歷性等極端復(fù)雜動力學(xué)行為的確定性系統(tǒng),因在生物學(xué)、醫(yī)學(xué)、電子信息等方面有巨大的應(yīng)用價值和潛力,吸引了各領(lǐng)域科學(xué)工作者廣泛研究混沌的建構(gòu)、非線性電路實現(xiàn)、保密通信、控制同步[1-3]等,并取得了豐碩的成果?;煦绲耐阶?990年Ott、Grobogi和Yorker提出OGY方法[4],成功控制混沌穩(wěn)定在周期軌道后,Pecora和Corroll首次提出混沌同步[5]的概念,并在實驗室用電路實現(xiàn)這一同步控制思想,引起人們的廣泛關(guān)注,迅速成為人們研究的熱點之一。研究者按照混沌不同的同步形式如延遲同步、完全同步、廣義同步、投影同步等,提出了如非線性控制、滑模控制、自適應(yīng)控制等多種混沌同步方法[6-8],混沌現(xiàn)象也開始從負(fù)面的作用走上正面利用的道路。文獻[9]提出一個包含H部件[10]、僅含一個非線性項的新三維非線性系統(tǒng)(Wu系統(tǒng)),并研究其記憶效應(yīng)和混沌特性。通過數(shù)值分析和設(shè)計電路實現(xiàn),驗證該系統(tǒng)是一個含廣義記憶元件、單渦卷吸引子,存在分岔、混沌和陣發(fā)混沌等復(fù)雜動力學(xué)行為的混沌系統(tǒng)。該文主要研究在非線性耦合函數(shù)作用下,該系統(tǒng)(Wu系統(tǒng))的同步以及同步控制器的尋獲問題。
武花干等人在文獻[9]提出一個含有H部件、僅含一個非線性項的新三維非線性動力學(xué)系統(tǒng),其模型如下。
理論分析了該系統(tǒng)的廣義記憶元件特性、混沌動力學(xué)行為,并用數(shù)值模擬和電路設(shè)計手段,相互驗證了該系統(tǒng)是一個含有H部件的簡單記憶混沌系統(tǒng)。對不同μ的參數(shù)值有不同的動力學(xué)行為,當(dāng)μ≥0.831時系統(tǒng)完全進入混沌狀態(tài)。系統(tǒng)(1)式狀態(tài)變量x隨參數(shù)μ變化的分岔圖如圖1所示。當(dāng)μ=0.95時,式(1)的部件H的廣義記憶元件的z(t)h(t)曲線如圖2所示,吸引子yz相圖如圖3所示,時域波形圖如4所示。
圖1 μ-x分岔圖
圖2 z(t)-h(t)曲線
圖3 y-z曲線
圖4 t-z曲線
研究還表明該簡單憶阻混沌系統(tǒng)僅存在一個確定的平衡點,其動力學(xué)行為依賴于可調(diào)參數(shù)μ。μ取不同值時,系統(tǒng)的行為可由倍周期分岔通向混沌道路,在μ=0.476、0.6、0.856、0.898處,存在陣發(fā)混沌現(xiàn)象。該簡單憶阻混沌系統(tǒng)(1)存在復(fù)雜的動力學(xué)行為。
該文用非線性函數(shù)耦合同步法[11-12]進行控制同步。設(shè)給定連續(xù)混沌系統(tǒng)為X˙(t)=F(X(t),t),適當(dāng)分離為:
其中X(t)∈Rn是系統(tǒng)的n維狀態(tài)矢量,G(X(t))為X˙(t)的線性部分,且有G(X(t))=AX(t),A為滿秩的常數(shù)矩陣,其特征值的實部全部為負(fù)。是的非線性部分。同理,可構(gòu)造一個新的連續(xù)混沌響應(yīng)系統(tǒng):
Y(t)∈Rn是系統(tǒng)的n維狀態(tài)矢量,則系統(tǒng)(3)和(4)的同步誤差定義為e(t)=Y(t)-X(t),可得=G(Y(t))-G(X(t))。
對于簡單憶阻混沌系統(tǒng)(1),驅(qū)動系統(tǒng)為:
響應(yīng)系統(tǒng)為:
設(shè)
可得驅(qū)動與響應(yīng)系統(tǒng)的誤差函數(shù)為
當(dāng)取非線性控制器為
則
其中滿秩的常數(shù)矩陣為A=則其三個特征根分別為λ1=-c、λ2=-d、λ3=-1。當(dāng)c>0、d>0時,滿秩矩陣A的三個特征根均小于零,根據(jù)Routh-Hurwitz判據(jù)知,誤差系統(tǒng)(9)漸近穩(wěn)定,則驅(qū)動系統(tǒng)(6)和響應(yīng)系統(tǒng)(7)達到同步。
該文用matlab數(shù)學(xué)軟件對簡單憶阻混沌系統(tǒng)(1)的誤差系統(tǒng)進行數(shù)值仿真研究。取驅(qū)動系統(tǒng)、響應(yīng)系統(tǒng)為和誤 差 系 統(tǒng) 初 值 分 別 為:x1(0)、y1(0)、z1(0)=-10,x2(0)、y2(0)、z2(0)=10,e1(0)、e2(0)、e3(0)=20,μ=0.95,c=2.5,d=2,則時間t與驅(qū)動和響應(yīng)系統(tǒng)對應(yīng)變量的變化圖如圖5所示,同步誤差效果如圖6所示??梢妕e1=3.426s、te2=2.815s和te3=3.836s時,誤差e1(t)、e2(t)和e3(t)已穩(wěn)定到零點,驅(qū)動和響應(yīng)系統(tǒng)達到同步。
該文利用非線性函數(shù)耦合混沌同步法,研究簡單記憶混沌系統(tǒng)的同步問題。通過誤差系統(tǒng)構(gòu)建特征值實部皆為負(fù)值的常滿秩矩陣,從而求解同步控制器,實現(xiàn)簡單記憶混沌系統(tǒng)的驅(qū)動與響應(yīng)系統(tǒng)達到同步目標(biāo)。該方法使得同步控制器更易于尋獲,理論和數(shù)值分析都驗證了該同步方法的正確性,為同步控制器易尋獲性提供了有益的思路。
圖5 時間與驅(qū)動、響應(yīng)系統(tǒng)對應(yīng)變量的變化圖
圖6 誤差系統(tǒng)同步響應(yīng)曲線
[1]Lorenz EN.Deterministic Nonperiodic Flow[J].Journalof the Atmospheric Sciences,1963:20(2):130–141.
[2]Buscarino,Arturo,F(xiàn)amoso et al..,Luigi.Chua's Circuit Dynamics Controlled by Using aWater Electrolytic Cell Component[J].INTERNATIONAL JOURNAL OF BIFURCATION AND CHAOS,2014,24(9):1910-1914.
[3]陸安山,陸益民.一種變形Liu混沌系統(tǒng)的分析及電路實現(xiàn)[J].江西師范大學(xué)學(xué)報:自然科學(xué)版,2014,38(5):32-35.
[4]Ott E,GrebogiC,Yorke JControlling chaos[J].Physical Review Letters,1990,64(11):1196-1199.
[5]Pecora L M,Carroll T L.synchronization in chaotic systems[J].PhysicalReview Letters,1990,64(8):821-824.
[6]Moreno A,Osorio M.Lyapunov approach to second-order sliding mode controllers and observers.In:Proceedings of the IEEE Conference on Dicision and Control.Cancun,Mexico:2008,2856-2861.
[7]Li Dekui,Lian Yuping.Hybrid function projective synchr onization and parameters identification of different dimension chaotica system[J].Journal of Qingdao Technological University,2013,34(3):95-101.
[8]Wang Jiankui,Chen Xiangqing,F(xiàn)u Jiankuan.Adaptive fi nite-time control of chaos in permanentmagnet synchronousmotor with uncertain parameters[J].NONLINEAR DYNAMICS,2014,78(2):1321-1328.
[9]武花干,包伯成,馮霏.簡單記憶混沌系統(tǒng)的動力學(xué)分析與電路實現(xiàn)[J].四川大學(xué)學(xué)報:工程科學(xué)版,2014(1):134-139.
[10]Di Ventra M,Pershin Y V,Chua L O.Circuit elements with memory:Memristors,memcapacitors,and meminductors[J].Proceedingsof the IEEE,2009,97(10):1717-1724.
[11]張旭東,俞建寧,郭蘭平,等.一個新混沌系統(tǒng)的非線性耦合同步的研究[J].云南民族大學(xué)學(xué)報:自然科學(xué)版,2011(02):149-151,156.
[12]NIEWen Jie,LANYue Heng,LIYong,etal..Generating large steady-state optomechanical entanglement by the action of Casimir force[J].Science China(Physics,Mechanics&Astronomy),2014,12:2276-2284.
[13]LiXiaojuan,Xu Zhenyuan,Xie Qingchun,etal..Generali zed synchronization of two different unidirection coupled Lorenz systems[J].Actaphysica sinica,2010,59(3):1532-1539.