安徽理工大學(xué)電氣與信息工程學(xué)院 徐甜甜

本文研究動力學(xué)特性更為復(fù)雜的新三維混沌系統(tǒng)。首先利用數(shù)值建模分析了三維混沌系統(tǒng)的基本動力學(xué)特性，然后搭建新混沌系統(tǒng)硬件電路，通過Multisim軟件進(jìn)行硬件電路仿真模擬，最后驗證了系統(tǒng)的物理可行性，結(jié)果表明仿真實驗與理論分析結(jié)論吻合。

1963年MIT(Massachusetts Institute of Technology)氣象學(xué)家Loren發(fā)現(xiàn)已確定的三階微分方程具有不規(guī)則的解，提出了“蝴蝶效應(yīng)”理論，開啟了研究混沌現(xiàn)象的序幕?；煦缱鳛榉蔷€性動力學(xué)的一個分支，在很多領(lǐng)域具有廣泛應(yīng)用。復(fù)雜混沌系統(tǒng)的產(chǎn)生、分析和控制近年來引起了國內(nèi)外同行的廣泛關(guān)注。經(jīng)典的混沌系統(tǒng)諸如：R?ssler系統(tǒng)、Chen系統(tǒng)及Lü系統(tǒng)等被提出，一些新的混沌系統(tǒng)被發(fā)現(xiàn)，它們具有更大的Lyapunov指數(shù)和更強(qiáng)的混沌特性。本文基于文獻(xiàn)中Lorenz-Like系統(tǒng)，搭建了新三維混沌系統(tǒng)，發(fā)現(xiàn)此系統(tǒng)的混沌特性比原系統(tǒng)復(fù)雜，在不同參數(shù)值下不僅折疊吸引子的渦卷數(shù)增加；并且發(fā)現(xiàn)在4.28＜b＜10.5時，系統(tǒng)產(chǎn)生新的兩翼折疊混沌吸引子，其最大Lyapunov指數(shù)高達(dá)6.7872，比上述文獻(xiàn)中混沌系統(tǒng)的Lyapunov指數(shù)值均大。

1 混沌系統(tǒng)模型及特性分析

1.1 混沌系統(tǒng)模型

本文基于Lorenz-Like系統(tǒng)構(gòu)建了一個新三維自治混沌系統(tǒng)，該系統(tǒng)的數(shù)學(xué)模型可描述為：

式中，x，y，z為狀態(tài)變量。當(dāng)初值為(10，10，60)，參數(shù)a=10、b=3、c=50、h=-1時，系統(tǒng)存在一個典型混沌吸引子如圖1所示。

圖1 系統(tǒng)(1)相圖

1.2 三維系統(tǒng)參數(shù)的影響

系統(tǒng)動力學(xué)特性隨參數(shù)的變化而變化，系統(tǒng)的運(yùn)行狀態(tài)可以直觀的由Lyapunov指數(shù)譜及分岔圖反映。當(dāng)固定參數(shù)a、c、h，參數(shù)b變化。圖2(a)反映在0＜b≤2.256及b＞12.39區(qū)域Lyapunov指數(shù)譜符號為(-，-，-)或(0，-，-)，系統(tǒng)處于周期運(yùn)動狀態(tài)；在2.257＜b≤12.35區(qū)域Lyapunov指數(shù)譜符號為(+，0，-)，系統(tǒng)處于混沌運(yùn)動狀態(tài)。圖2(b)與圖2(a)的動力學(xué)特性相對應(yīng)，反映在4.28＜b＜10.5區(qū)域，系統(tǒng)表現(xiàn)出不同的動力學(xué)特性。取b=8時，系統(tǒng)(1)的相圖如圖3所示，其吸引子軌跡在特定的區(qū)域內(nèi)具有遍歷性，與系統(tǒng)(1)典型混沌吸引子及文獻(xiàn)中所描述吸引子均不相同；且在此區(qū)間內(nèi)計算系統(tǒng)(1)的最大Lyapunov指數(shù)，得知在b=7.5時的最大Lyapunov指數(shù)為6.7872，比原Lorenz-Like系統(tǒng)及典型系統(tǒng)參數(shù)值下Lyapunov指數(shù)值均大，系統(tǒng)表現(xiàn)出更復(fù)雜的動力學(xué)狀態(tài)。

圖2 三維系統(tǒng)隨b變化的Lyapunov指數(shù)譜及分岔圖

圖3 b=8三維系統(tǒng)相圖

2 三維混沌系統(tǒng)電路的設(shè)計與實現(xiàn)

混沌系統(tǒng)的物理可行性可通過模擬電子電路實現(xiàn)，采用模塊化的設(shè)計方法驗證系統(tǒng)的動態(tài)行為。系統(tǒng)由運(yùn)算放大器、模擬乘法器等電子元件構(gòu)成。由于電子元件工作電壓幅值受限，將系統(tǒng)輸出電壓幅值降低為原來的1/10，并進(jìn)行時間尺度變換。令，當(dāng)b=3時，得：

由圖4得其狀態(tài)方程為：

圖4中C1=C2=C3=100nF，比較(2)、(3)式得電阻R1=R2=R4=R7=R9=R10=R12=R13=R14=R15=10kΩ，R3=5kΩ、R5=20kΩ、R6=R8=100kΩ、R11=33.3kΩ。此時Multisim仿真相圖如圖5(a)、(b)所示。當(dāng)b=8時，由圖5系統(tǒng)的相圖可知，為避免飽和現(xiàn)象將系統(tǒng)狀態(tài)變量z調(diào)整為原來的1/20，此時硬件系統(tǒng)中R3=1kΩ、R7=0.5kΩ、R8=50kΩ、R10=2kΩ、R11=12.5kΩ，其余參數(shù)不變，仿真相圖如圖5(c)、(d)所示。

結(jié)論：本文通過改進(jìn)Lorenz-Like系統(tǒng)方程產(chǎn)生一個新的三維自治混沌系統(tǒng)，研究了系統(tǒng)的復(fù)雜動力學(xué)特性。從系統(tǒng)相圖發(fā)現(xiàn)新混沌系統(tǒng)混沌吸引子比Lorenz-Like系統(tǒng)復(fù)雜，折疊吸引子的渦卷數(shù)增加。從系統(tǒng)參數(shù)的Lyapunov指數(shù)圖譜和分岔圖可以看出，參數(shù)變化系統(tǒng)所經(jīng)歷的周期與混沌運(yùn)動狀態(tài)。為驗證新混沌系統(tǒng)的物理可行性，根據(jù)系統(tǒng)數(shù)學(xué)模型搭建系統(tǒng)硬件電路，電路仿真結(jié)果與Matlab仿真結(jié)果相同，表明新建混沌系統(tǒng)的正確性。

圖4 電路原理圖

圖5 系統(tǒng)(1)電路實現(xiàn)仿真相圖