宋 爽，湯 瓊，瞿民凱，趙思遠

（湖南工業(yè)大學(xué) 理學(xué)院，湖南 株洲 412007）

1 背景知識

1979年，學(xué)者O.E.R?ssler首次發(fā)現(xiàn)了一個可以同時產(chǎn)生兩個方向上雙曲不穩(wěn)定的混沌奇異吸引子的四維動力學(xué)系統(tǒng)[1]，并將之命名為超混沌現(xiàn)象。相較于常見的三維混沌系統(tǒng)，超混沌系統(tǒng)在實際應(yīng)用中更能滿足通信保密、圖像加密、人工神經(jīng)網(wǎng)絡(luò)、非線性電路等領(lǐng)域的需求。超混沌系統(tǒng)可以通過引入新的狀態(tài)變量反饋或基于憶阻器元件實現(xiàn)，得到四維及以上的高維動力學(xué)系統(tǒng)[2-5]，并通過改變非線性項取值范圍或憶阻器的控制參數(shù)，獲得倍周期分叉、單渦卷吸引子、雙渦卷吸引子、單方向多渦卷吸引子，甚至網(wǎng)格多渦卷吸引子等混沌現(xiàn)象，現(xiàn)已成為混沌系統(tǒng)中極具研究價值的一個方向。

相較于低維混沌系統(tǒng)較為簡單的動力學(xué)特性和不夠多的密鑰參數(shù)，超混沌系統(tǒng)具有更復(fù)雜的動力學(xué)特性和更豐富的密鑰參數(shù)。判定一個動力學(xué)系統(tǒng)是否為超混沌系統(tǒng)，最直觀的判斷標準就是具有兩個或兩個以上的正Lyapunov指數(shù)。

近幾十年來，國內(nèi)外學(xué)者一直在研究生成超混沌系統(tǒng)的方法。受學(xué)者Li Q.D.[6]發(fā)現(xiàn)超混沌Chen系統(tǒng)的啟發(fā)，科研工作者們多通過引入新狀態(tài)變量反饋法生成超混沌系統(tǒng)[7-14]。由于引入狀態(tài)變量的反饋控制器多為非線性項，它可進一步豐富動力學(xué)系統(tǒng)的混沌特征，增加系統(tǒng)的不可預(yù)測性。該系統(tǒng)的電路實現(xiàn)取決于反饋量非線性項的復(fù)雜程度，如何找到簡單且易實現(xiàn)的超混沌系統(tǒng)一直是混沌電路的研究目標。

雖然關(guān)于超混沌系統(tǒng)的構(gòu)造已經(jīng)有了大量的研究，但是由于混沌吸引子較為復(fù)雜的動力學(xué)特性，現(xiàn)有關(guān)于多渦卷混沌系統(tǒng)電路實現(xiàn)的報道仍然比較少。因此，對不同非線性控制函數(shù)構(gòu)造的超混沌系統(tǒng)進行詳細的動力學(xué)分析，并搭建出對應(yīng)的電路模型，是十分必要的。因此，本研究在三維Lü系統(tǒng)[7]的基礎(chǔ)上加以修改，構(gòu)建了一個基于狀態(tài)變量反饋的四維動力學(xué)系統(tǒng)，經(jīng)計算Lyapunov指數(shù)后，證實其為超混沌系統(tǒng)。然后，對其進行了基本的平衡點穩(wěn)定性分析、耗散性分析、對稱性分析等動力學(xué)分析，并通過計算Lyapunov指數(shù)及維數(shù)、繪制Poincaré截面圖等方法對該超混沌系統(tǒng)進行了數(shù)值仿真分析。最后，根據(jù)反相器、積分器等基本的運算放大電路搭建出其相應(yīng)的電路模型。本研究基于憶阻器反饋設(shè)計新的超混沌系統(tǒng)，該方法在增加系統(tǒng)維度的同時豐富了系統(tǒng)的混沌特征，且對應(yīng)的混沌電路設(shè)計更易實現(xiàn)，本研究可為混沌系統(tǒng)的電路實現(xiàn)提供參考。

2 超混沌系統(tǒng)的數(shù)學(xué)模型

在三維Lü系統(tǒng)基礎(chǔ)上進行修改，構(gòu)建了一個基于狀態(tài)變量反饋的四維動力學(xué)系統(tǒng)，新的超混沌系統(tǒng)的狀態(tài)方程如下：

式中：a、b、c、d、e、g均為實數(shù)，且a、b、c、d、e均為正實數(shù)，e為控制參數(shù)。

由式（1）可知，該混沌系統(tǒng)是基于狀態(tài)變量反饋控制實現(xiàn)的超混沌系統(tǒng)。

當(dāng)系統(tǒng)（1）取a=35、b=35、c=1、d=3、e=3、g=-27時，得到該超混沌系統(tǒng)的相軌跡圖如圖1所示。

圖1 超混沌系統(tǒng)（1）的相軌跡圖Fig.1 Phase diagram of the hyperchaotic system (1)

通過觀察圖1所示超混沌系統(tǒng)的相軌跡圖可以得知，該超混沌系統(tǒng)的渦卷數(shù)目為2個，有明顯的混沌現(xiàn)象。

3 超混沌系統(tǒng)的動力學(xué)分析及數(shù)值仿真分析

3.1 超混沌系統(tǒng)的動力學(xué)分析

由于非線性混沌系統(tǒng)的解析解通常無法輕易獲得，所以需要通過平衡點穩(wěn)定性分析、耗散性分析等定性分析的方法獲得系統(tǒng)的近似形態(tài)。

3.1.1 平衡點穩(wěn)定性分析

對于一個四階非線性自治系統(tǒng)，經(jīng)平衡點線性化后得到其雅可比矩陣，根據(jù)解得的特征值性質(zhì)的不同，會得到不同性質(zhì)的平衡點。下面針對每一個平衡點逐一展開討論。

令系統(tǒng)（1）的右邊全部等于0，則可以得到如下方程組：

解得系統(tǒng)（1）的平衡點有兩種情況：

1）d=0時，系統(tǒng)無平衡點；

2）d≠0時，系統(tǒng)存在平衡點P1（0, 0,, 0）。

3.1.1.1 系統(tǒng)（1）無平衡點時的超混沌吸引子

取a=8、b=0.1、c=1、d=0、e=2、g=-8時，初值為[0 1 1 0]時，系統(tǒng)（1）沒有平衡點，且此時該新超混沌系統(tǒng)的Lyapunov指數(shù)為：

可見，系統(tǒng)（1）滿足存在兩個正的Lyapunov指數(shù)的條件，故該系統(tǒng)存在隱藏的超混沌吸引子，其相圖如圖2所示。

圖2 系統(tǒng)（1）隱藏超混沌吸引子的相軌跡圖Fig.2 Phase diagram of the hidden hyperchaotic attractor of system (1)

通過觀察圖2所示的相軌跡圖可知，該系統(tǒng)有明顯的混沌現(xiàn)象。

3.1.1.2 系統(tǒng)（1）有一個平衡點時的超混沌吸引子

當(dāng)d≠0時，系統(tǒng)存在平衡點P1(0, 0,, 0)。其中，g可取任意實數(shù)，則該系統(tǒng)有無限平衡點集。在平衡點處做線性化處理，可得到系統(tǒng)（1）線性化后的雅可比矩陣為

將式（3）代入平衡點坐標后，令|λI-J|=0，則系統(tǒng)（1）在平衡點處線性化后得到的特征方程為

特征方程的第一個特征值為λ1=-d，設(shè)

則根據(jù)Routh-Hurwitz判據(jù)，可得：

則有

若系統(tǒng)（1）同時滿足

當(dāng)系統(tǒng)（1）取a=35、b=35、c=1、d=3、e=3、g=-27時，在平衡點P1處對系統(tǒng)（1）進行線性化處理后，求得對應(yīng)特征方程（4）的解為

此時，特征方程（4）的解為兩正兩負的實特征值，且所有特征值之和滿足λ1+λ2+λ3+λ4<0，故此時系統(tǒng)（1）的平衡點P1是具有二維穩(wěn)定流形和二維不穩(wěn)定流形的鞍焦平衡點。將此時系統(tǒng)（1）的各參數(shù)取值代入后，式子仍然成立，故P1是不穩(wěn)定的鞍焦點，系統(tǒng)呈現(xiàn)超混沌現(xiàn)象，此時系統(tǒng)（1）的相圖如圖3所示，可見結(jié)果符合前文論述結(jié)論。

圖3 系統(tǒng)（1）有一個平衡點時的超混沌吸引子相軌跡圖Fig.3 Phase diagram of the Hyperchaotic system (1)with an equilibrium point

3.1.2 對稱性分析

對于超混沌系統(tǒng)（1）做(x,y,z,w)～ (-x, -y, -z,-w)變換后，發(fā)現(xiàn)該系統(tǒng)的狀態(tài)方程仍保持不變。故該系統(tǒng)關(guān)于z軸對稱，且這種對稱關(guān)系對于系統(tǒng)的所有參數(shù)都成立。

3.1.3 耗散性分析

系統(tǒng)（1）的相空間容積率為

當(dāng)a-c-d<0時，系統(tǒng)是耗散的，相空間的體積元將以指數(shù)形式收斂：

解得V=V0e-(a+c+d)t。

即當(dāng)時間趨近于無窮大時，系統(tǒng)的每一條軌線都將以-(a+c+d)的速度指數(shù)收斂至0。最終，系統(tǒng)的所有軌線收斂至體積為0的集合，系統(tǒng)的所有軌線都固定于一個吸引子上。

當(dāng)系統(tǒng)（1）取a=35、b=35、c=1、d=3、e=3、g=-27時，=-(35+1+3)=-39<0，即此時的系統(tǒng)（1）是耗散的。

3.2 超混沌系統(tǒng)的數(shù)值仿真分析

由于混沌系統(tǒng)的解對初值具有敏感依賴性，從而很難證明所獲得的近似解是否會收斂于系統(tǒng)的真實解。因此，除了對系統(tǒng)進行定性分析之外，還需要通過數(shù)值分析方法，如Lyapunov指數(shù)譜、Poincaré圖等，對混沌系統(tǒng)進行進一步的定量分析，以彌補定性分析的某些不足，從而更全面地研究動力學(xué)系統(tǒng)的混沌特征。

3.2.1 Lyapunov指數(shù)譜及維數(shù)

當(dāng)取a=35、b=35、c=1、d=3、e=3、g=-27時，可以計算得：

此時，系統(tǒng)（1）滿足存在兩個正的Lyapunov指數(shù)，處于超混沌狀態(tài)，其Lyapunov指數(shù)圖見圖4。

圖4 系統(tǒng)（1）的Lyapunov指數(shù)圖Fig.4 Lyapunov exponential diagram of system (1)

Lyapunov維數(shù)計算公式為

其中，k滿足以下條件：

將Lyapunov指數(shù)代入后算得，

故可知該超混沌系統(tǒng)為分數(shù)維，具有豐富的動力學(xué)特性。

3.2.2 Poincaré截面圖

當(dāng)選取參數(shù)a=35、b=35、c=1、d=3、e=3、g=-27時，繪制出該系統(tǒng)的Poincaré截面圖，如圖5所示。

圖5 系統(tǒng)（1）的Poincaré截面圖Fig.5 Poincaré diagram of system (1)

觀察圖5所示的系統(tǒng)（1）Poincaré截面圖可以發(fā)現(xiàn)，不論是yz相平面還是xw相平面的Poincaré截面圖，上面都有著一片密密麻麻雜亂分布的點，結(jié)合該系統(tǒng)有兩個正Lyapunov指數(shù)的特征，可以得知此時系統(tǒng)（1）處于超混沌狀態(tài)。

4 超混沌系統(tǒng)的電路實現(xiàn)

本文采取模塊化設(shè)計思路，運用簡單的電感、電容、運算放大器等元件實現(xiàn)電路的搭建。

由于要運用到反相器等運算放大電路，需先將方程組（1）改寫為

為了使信號在時域上擴展而幅度不發(fā)生變化，需要先對各變量取值范圍進行核查，若各變量取值超過元件額定值，則需要對系統(tǒng)的狀態(tài)方程進行時間尺度變換：

令τ=τ0，其中τ0是時間尺度變換因子，取τ0=100 μs，則方程組可改寫為

根據(jù)式（7），可搭建出其對應(yīng)的輸出為x、y、z、w的電路圖，如圖6所示。

圖6 系統(tǒng)（1）的Simulink電路圖Fig.6 Simulink simulation circuit of system (1)

圖6中，當(dāng)運算放大器工作于線性狀態(tài)時，其兩個輸入端可視為等電位。由于運算放大器的差模輸入電阻較大，故流入運算放大器輸入端的電流可以小到忽略不計，即可以將運算放大器的兩個輸入端視作開路，得到：

圖中的電容是為了避免因參數(shù)選取過大導(dǎo)致進入混沌前時間過長，特選取C1=C2=C3=C4=100 nF，比較式（7）和（8），可得：

解得：R1=R2=R3=2.86 kΩ，

R4=R6=R7=100 kΩ，

R5=R9=0.1 kΩ，

R8=R10=33.33 kΩ。

故可得系統(tǒng)（1）的電路圖如圖7所示。

圖7 系統(tǒng)（1）的電路圖Fig.7 Circuit realization diagram of system (1)

5 超混沌系統(tǒng)電路的仿真驗證

為了證明設(shè)計電路的正確性，在Simulink中搭建了對應(yīng)的模擬電路。本文采用了模塊化設(shè)計的思路，運用運算放大器、電容、電感等基本元件組成反相器，與積分器一起搭建出模擬電路。其中，給定運算放大器的門限電壓設(shè)置為±15 V，考慮到乘法器AD633JN的乘法因子為0.1/V，R5、R9的取值變更為1 kΩ，并按照計算出的參數(shù)設(shè)置元件的值，觀察示波器輸出的結(jié)果，所得仿真結(jié)果如圖8所示。

圖8 系統(tǒng)（1）的Simulink仿真結(jié)果圖Fig.8 Simulation diagrams of system (1)

將圖8所示仿真結(jié)果與吸引子相軌跡圖進行對比，發(fā)現(xiàn)該超混沌電路的實驗結(jié)果與Matlab的數(shù)值仿真結(jié)果一致，這證明所設(shè)計的電路準確有效。

6 結(jié)語

相比于常見的三維混沌系統(tǒng)，超混沌系統(tǒng)有著更為復(fù)雜的混沌特征，更能滿足現(xiàn)代通信保密、圖像加密、人工神經(jīng)網(wǎng)絡(luò)、非線性電路等領(lǐng)域的需求。如何找到結(jié)構(gòu)簡單、動力學(xué)特性豐富的超混沌系統(tǒng)，并且將其應(yīng)用到實際生產(chǎn)與生活中，是混沌理論一直追求的目標。

本文采用狀態(tài)變量反饋法構(gòu)造了一個新的超混沌系統(tǒng)，通過根據(jù)平衡點對應(yīng)的雅可比矩陣計算出相應(yīng)的特征值判定平衡點的穩(wěn)定性等定性分析方法，以及計算Lyapunov指數(shù)和維數(shù)、繪制Poincaré截面圖等定量分析方法，對超混沌系統(tǒng)進行了基本的動力學(xué)分析。通過基本的運算放大電路，搭建了新超混沌系統(tǒng)的電路模型，并利用Matlab-Simulink仿真模塊庫驗證了電路實現(xiàn)的準確性，為超混沌系統(tǒng)的研究提供了新的思路。