一個超混沌類Lorenz系統的非線性動力學行為及計算機仿真

徐鴻鵬，尹社會，張勇

（河南工業職業技術學院河南南陽473000）

一個超混沌類Lorenz系統的非線性動力學行為及計算機仿真

徐鴻鵬，尹社會，張勇

（河南工業職業技術學院河南南陽473000）

利用Mat1ab軟件和數學微分方程理論分析給出了一個新五維超混沌類Lorenz系統的非線性動力學特性。通過定性分析和定量分析相結合的手段探討了主要包括對稱性、耗散性、平衡點的穩定性、空間相圖、時序波形圖等方面的非線性動力學行為，并運用Wo1f方法計算出了系統的Lyapunov指數和Lyapunov維數，結合系統的不同運動狀態的分岔圖、Poincaré映射圖和功率譜圖等手段進一步表明該系統具有復雜的動力學特性，為進一步的混沌控制、同步與加密通信等工程應用提供了理論依據。

類Lorenz超混沌系統；分岔；耗散性；功率譜

著名的Lorenz系統是研究混沌現象的典型范例［1_2］。自從Lorenz在一個簡單的三維自治系統中首先發現了蝴蝶混沌吸引子之后，又有新的混沌吸引子不斷被發現，1999年，Chen系統和Lü系統相繼被提出［3_4］，隨后一個統一的混沌系統也被陳關榮和呂金虎提出并研究［5］。此后，國內外不少學者相繼提出新的混沌或超混沌系統［6_7］，新系統的提出促進了人們對混沌現象的深入研究和認識，提高了混沌理論在工程上的應用能力。文獻［8］提出了構造新超混沌系統的必要條件并構造了一類新的五維超混沌類Lorenz系統［8］，但并未進行動力學方面的研究，文中通過數值仿真給出了系統的相圖、時序波形圖、分岔圖、Poincaré映射、功率譜圖等，結果驗證了該系統屬于超混沌系統及其豐富的動力學行為。

1　數學模型及其吸引子的存在性

章秀君等構造了一個五維超混沌類Lorenz系統的方程為［：

其中：（x，y，z，w，u）∈R5為狀態變量，a，b，c，d為系統實參數。當a=10，b=8/3，c=28，d=2時，系統相空間的時間平均散度為

可知系統為一個耗散非線性動力系統，其軌線隨著時間不斷演化到一個不變的吸引子集合中，其初值?。?.1，0.4，0.1，0.1，0.1），系統（1）的吸引子軌線的相圖如圖1所示。

2　系統的LyaPunoV指數和LyaPunoV維數

由數值計算可得系統的5個Lyapunov指數分別為：λ1=1.863 9，λ2=0.043 589，λ3=_0.733 93，λ4=_4.612 6，λ5=_12.117 6，其中有兩個Lyapunov指數大于0，說明系統處于超混沌狀態。進一步通過Kap1an_Yorke猜想可以得到其分形維數為：

可見系統具有分數維數，進一步驗證了系統處于混沌狀態。

圖1　吸引子相圖

3　系統在超混沌態的時序波形、功率譜和Polncaré映射

對混沌系統的吸引子分析，其空間結構十分復雜，軌線無窮延伸、壓縮和折疊，其軌線在特定的吸引域內具有遍歷性。時域分析如圖2所示也表明，序列具有典型的非周期性，且具有對初值敏感依賴性。各個變量隨時間演化的時序波形圖如圖2所示。如果初值發生細微的變化，系統的行為會發生明顯的變化，如圖3所示為初值?。?.1，0.4，0.1，0.1，0.1）和（0.1，0.4+0.00001，0.1，0.1，0.1）的時序圖形。吸引子的非周期性也可以通過系統的連續功率譜表現出來，采用快速Fourier變換法，對系統的第一個狀態變量x的時間序列進行功率譜分析如圖4所示。由于是非周期運動，所以表現出連續功率譜，又由于存在分叉現象，所以有峰值出現。

圖2　時序波形圖

混沌系統的Poincaré截面上是沿幾條曲線弧分布著一些具有分形結構的密集點，如圖5所示，圖中給出了x_y平面上的Poincaré映射，進一步說明了系統處于混沌態。

4　系統的平衡點分析

令系統（1）右邊等于0，當a=10，b=8/3，c=28，d=2時，經計算系統有無窮多個平衡點：S0（0，0，0，0，0），S1（0，c，0，0，_10c），S2（4.32，0，14.0，0，43.2），S3（_4.32，0，14.0，0，_43.2），S4（_ 13.3，25.3，1.33，16.9，_284.0）。

圖3　初值敏感性時序波形圖

圖4　功率譜圖

圖5　Poincaré映射

在平衡點S0（0，0，0，0，0）處線性化系統（1）可得其Jacobi矩陣，

由|λI_J|=0可得5個特征值為λ1=_1，λ2=_10，λ3=_2.6667，λ4=_2，λ5=0。其中4個負的實根表示在這4個方向收縮，因此平衡點S0（0，0，0，0，0）為穩定的結點。

在系列平衡點S1（0，c，0，0，_10c）處線性化系統（1）可得其Jacobi矩陣（為分析方便取c=1），

由|λI_J|=0可得5個特征值為λ1=_1，λ2=_10.099，λ3= 0.099，λ4=_2，λ5=_2.6667。其中4個負的實根表示在這2個方向收縮，1個正的實根表示在這1個方向擴張，因此平衡點S1（0，c，0，0，_10c）為不穩定的指標1的鞍點。

在平衡點S2（4.32，0，14.0，0，43.2）處線性化系統（1）可得其Jacobi矩陣，

由|λI_J|=0可得5個特征值為λ1=_19.4534，λ2=5.5193+ 3.1942i，λ3=5.5193_3.1942i，λ4=_6.1224，λ5=_1.1295。其中3個負的實根表示在這3個方向收縮，一對共軛復根的實部為正說明平衡點S2（4.32，0，14.0，0，43.2）為指標2的不穩定鞍點。

在平衡點S3（_4.32，0，14.0，0，_43.2）處線性化系統（1）可得其Jacobi矩陣，

由|λI_J|=0可得5個特征值為λ1=_6.9549+12.1504i，λ2=_ 6.9549_12.1504i，λ3=_0.4746+4.2467i，λ4=_0.4746_4.2467i，λ5=_0.8076。其中1個負的實根表示在這1個方向收縮，2對共軛復根的實部為負說明平衡點S3（_4.32，0，14.0，0，_43.2）為穩定焦點。

在平衡點S4（_13.3，25.3，1.33，16.9，_284.0）處線性化系統（1）可得其Jacobi矩陣，

由|λI_J|=0可得5個特征值為λ1=_9.1149，λ2=_2.6059+ 4.5315i，λ3=_2.6059_4.5315i，λ4=_0.67+0.3003i，λ5=_0.67_ 0.3003i。其中1個負的實根表示在這1個方向收縮，2對共軛復根的實部為負說明平衡點S4（_13.3，25.3，1.33，16.9，_284.0）為穩定焦點。

通過以上分析，在參數a=10，b=8/3，c=28，d=2下，系統產生混沌現象。

5　系統參數的影響

非線性動力系統的動力學行為主要由系統參數決定，隨著系統參數的變化，系統表現出極限環（周期軌或擬周期軌）和奇怪吸引子等不同的非線性行為，即出現Hopf分叉和混沌現象。分岔圖顯示，當固定其余參數只有一個參數發生變化時，系統表現出不同的動力學行為，為敘述方便，這里只列出改變參數的區間范圍。在a∈［10，16］，b∈［2.4，4.5］，c∈［8，28.5］，d∈［0，20］，區間內變化，系統表現出混沌或超混沌行為，如圖6所示。下面通過吸引子相圖可以進一步驗證分岔圖中的結論，如圖7所示。

圖6　參數對變量x的分岔圖

6　結論

文中研究了一種五維超混沌類Lorenz系統，對系統的基本動力學行為進行了深入的研究，包括功率譜、Poincare截面、分岔圖等。分析表明該系統具有豐富的動力學行為，系統在參數變化時的動力學行為的演變呈現出周期、復雜周期（擬周期）、混沌以及超混沌運動，這些結論為系統的電子振蕩電路的實現和通信工程設計等應用提供了理論依據。

圖7　周期軌或擬周期軌相圖

［1］Lorenz E.N.Deterministic non_periods f1ows［J］.Journa1 of Atmosphere Science，1963，20（2）：130_141.

［2］尹社會，張勇，張付臣，等.基于Lorenz系統的強迫Lorenz混沌系統的動力學研究［J］.東北師大學報：自然科學版，2014，46（1）：42_47.

［3］Guan_rong CHEN，Tetsushi UETA.Yet another chaotic attractor［J］.Internationa1 Journa1 of Bifurcation and Chaos，1999，9（7）：1465_1466.

［4］Jin_hu Lü，Guan_rong Chen.A new chaotic attractor coined［J］.Internationa1 Journa1 of Bifurcation and Chaos，2002，12 （3）：659_661.

［5］Jin_hu Lü，Guan_rong Chen，Dai_zhan CHENG，etc.Bridge the gap between the Lorenz system and the Chen system［J］. Internationa1 Journa1 of Bifurcation and Chaos，2002，12 （12）：2917_2926.

［6］尹社會，張勇，皮小力.自治混沌系統的動力學行為及計算機仿真［J］.廣西物理，2015，36（1）：32_37.

［7］高智中.一個新超混沌系統及其線性反饋同步［J］.中山大學學報：自然科學版，2012，51（6）：30_34.

［8］章秀君，吳志強，方正.超混沌系統的構造方法研究［J］.計算機工程與應用，2014，50（2）：92_98，151.

Nonllnear dynamlc behaVlor of a lorenz-llke hyPerchaotlc system and lts comPuter slmulatlon

XU Hong-peng，YIN She-hui，ZHANG Yong
（Henan Polytechnic Institute，Nanyang 473000，China）

The non1inear characteristic properties of a nove1 five-dimensiona1 Lorenz-1ike hyper-chaotic system is further investigated by differentia1 equation theoretica1 ana1ysis and simu1ative ana1ysis based on Mat1ab software.The symmetry of system，dissipation，the stabi1ity of the equi1ibrium points，the phase diagram，time domain waveform，Lyapunov exponents and Lyapunov dimension，bifurcation diagram，Poincaré mapping diagram and the power spectrum are given by Mat1ab software.The resu1ts showed the nove1 chaotic system has rich dynamic behavior.The theoretica1 basis is provided for chaotic contro1，sychronization and encryption communications engineering and so on.

Lorenz-1ike hyper-chaotic system；bifurcation；dissipation；power spectrum

TN911.6；TN918.1

A

1674_6236（2016）10_0038_04

2015_07_13稿件編號：201507087

河南省科學技術發展計劃項目（142300410416）

徐鴻鵬（1980—），男，河南南陽人，講師。研究方向：計算物理和大學物理教學。