梁欣濤++高華強++康守強++朱建良++王玉靜++鄭勢

摘要：構(gòu)建了新的五維混沌系統(tǒng)，進行離散混沌模型的仿真，給出了系統(tǒng)的混沌吸引子相圖，對該系統(tǒng)的耗散性、吸引子的存在性、平衡點的穩(wěn)定性、Lyapunov指數(shù)及維數(shù)、功率譜、Poincare截面圖、Lyapunov指數(shù)譜、分岔圖特性進行分析，結(jié)果表明該系統(tǒng)具有混沌特性，有復雜的動力學行為，且該行為對系統(tǒng)參數(shù)具有敏感性.為了使混沌得到更廣泛應用，采用數(shù)字電路實現(xiàn)該系統(tǒng)，對離散化的五維混沌系統(tǒng)進行Modesim仿真，將VHDL程序配置到FPGA中，并利用數(shù)模轉(zhuǎn)換模塊在示波器上觀測到了該系統(tǒng)的混沌吸引子相圖.數(shù)字電路實驗結(jié)果與離散模型仿真分析是一致的，進一步從物理實現(xiàn)上說明了系統(tǒng)的混沌特性.

關(guān)鍵詞：混沌系統(tǒng)；分岔圖；Lyapunov指數(shù)；電路實現(xiàn)

DOI： 10.15938/j.jhust.2015.03.020

中圖分類號：TN911.73

文獻標志碼：A

文章編號：1007-2683（2015）03-0101-05

0 引 言

構(gòu)造全新的混沌系統(tǒng)或改進型混沌系統(tǒng)，在此基礎(chǔ)上，對其混沌特性及其應用進行分析研究，這是目前國內(nèi)外研究混沌的一個熱點課題.混沌系統(tǒng)主要有離散混沌系統(tǒng)和連續(xù)系統(tǒng)兩大類.離散混沌系統(tǒng)典型的有一維拋物映射（Logistic映射）和二維的Henon映射.連續(xù)混沌系統(tǒng)較多，典型的有廣義Lorenz系統(tǒng)族、Rossler系統(tǒng)、Chua系統(tǒng)、Sprott系統(tǒng)、Chen系統(tǒng)、Lu系統(tǒng)、Liu系統(tǒng)、Qi系統(tǒng)等.近年來，為了構(gòu)建新的復雜的混沌系統(tǒng)，學者們利用各種方式來構(gòu)建新的高維混沌系統(tǒng).對新混沌系統(tǒng)的構(gòu)建與分析進一步豐富和完善了混沌理論，為混沌應用提供了一些新的技術(shù)手段，從而促進了混沌在自然科學、電子、通信以及其他工程應用領(lǐng)域的發(fā)展.具體的應用比如：研制混沌信號發(fā)生器、高容量動態(tài)信息存儲器、信息加密、保密通信、信號檢測與處理等.在某些應用中，混沌系統(tǒng)如果采用模擬器件實現(xiàn)，由于元器件參數(shù)的離散性等因素，使應用系統(tǒng)的實現(xiàn)很困難，解決該問題的有效途徑是基于連續(xù)混沌系統(tǒng)離散化和數(shù)字化處理技術(shù)來實現(xiàn)混沌序列及算法，進而利用先進數(shù)字處理器件與技術(shù)來實現(xiàn).該方法為混沌的應用，尤其是在混沌保密通信領(lǐng)域中的應用提供了強大的技術(shù)支持.

本文構(gòu)造了一個新五維二次的混沌系統(tǒng).該系統(tǒng)每個方程中各含有一個二次的非線性交叉乘積項，所需乘法器數(shù)量少，實現(xiàn)簡單.對新五維混沌系統(tǒng)進行數(shù)值模擬，對系統(tǒng)的耗散性、吸引子的存在性、平衡點的穩(wěn)定性、Lyapunov指數(shù)及維數(shù)、功率譜、Poincare截面圖動力學特性進行研究，根據(jù)分岔圖和Lyapunov指數(shù)譜詳細分析了混沌行為的系統(tǒng)參數(shù)敏感性，其中部分參數(shù)在很大范圍內(nèi)呈現(xiàn)混沌.最后，利用FPGA實現(xiàn)了新五維混沌系統(tǒng)的硬件電路，在示波器上觀察到混沌吸引子相圖，證實了該系統(tǒng)的可實現(xiàn)性.

1 新的五維混沌系統(tǒng)

本文提出的新五維混沌系統(tǒng)的數(shù)學模型為：其中， 為常數(shù)，當

時，系統(tǒng)存在典型的部分混沌吸引子，如圖l所示.

由圖1可以看出，提出的五維混沌系統(tǒng)所產(chǎn)生的混沌吸引子相圖清晰、飽滿.由于該混沌系統(tǒng)對初值極為敏感，它表現(xiàn)為局部不穩(wěn)定，而從相圖的形成過程可看出系統(tǒng)又是從暫態(tài)向漸進穩(wěn)態(tài)運動，尋求穩(wěn)態(tài)，系統(tǒng)的運動軌跡靠近又分開，分開又反折而靠近，來回折疊無數(shù)次，形成復雜吸引子結(jié)構(gòu).

2 基本動力學特性

2.1耗散性和吸引子的存在性

由于

即隨著時間的推移，包含系統(tǒng)軌跡的每個體積元以指數(shù)率 收縮到零.這種體積收縮作用將使相軌跡必須折回來，即產(chǎn)生折疊運動.拉伸運動和折疊運動兩者相互作用的結(jié)果，只能是形成具有分形和分維的混沌運動，因此，從該角度定性分析出了系統(tǒng)（1）可形成混沌吸引子.

2.2平衡點及穩(wěn)定性

系統(tǒng)（1）的平衡點可解下列代數(shù)方程組得到：

給定系統(tǒng)（1）中的參數(shù)值后，系統(tǒng)（1）的6個平衡點分別為：

考察系統(tǒng)的穩(wěn)定性，對系統(tǒng)在各平衡點處線性化，得到Jacobi矩陣并計算各平衡點對應的特征值，在平衡點So處的Jacobi矩陣為

根據(jù) ，得到其特征值為

.這里5個特征根的實部有正也有負，根據(jù)Routh-Hurwitz條件，可得平衡點SO是不穩(wěn)定的鞍點.同理，可得到其他平衡點對應的特征值，結(jié)果如表1所示.

從表1可以看出，每個平衡點對應的所有特征值中至少有一個實部為正，且至少有一個實部為負，因此系統(tǒng)（1）的所有平衡點均是不穩(wěn)定的鞍焦點.

2.3 Lyapunov指數(shù)與Lyapunov維數(shù)

Lyapunov指數(shù)（簡寫為LE）是混沌系統(tǒng)中定量描述狀態(tài)空間混沌吸引子軌線彼此排斥和吸引的且.本文利用LET程序包，計算得到系統(tǒng)（1）的所有Lyapunov指數(shù)分別為 .如圖2所示.可見，該系統(tǒng)具有正的Lyapunov指數(shù)，是混沌系統(tǒng).

新五維混沌系統(tǒng)Lyapunov指數(shù)的維數(shù)為：

這里， ，其中j是保證 的最大 值.因此可求得 的大小為：

即該系統(tǒng)LE的維數(shù)是分數(shù)維數(shù)，也就是所謂的分維，這點也證明混沌的存在.

2.4 時域波形、功率譜及Poincare截面圖

混沌系統(tǒng)的時域波形具有非周期性，以分量X3和 為例，從圖3的（a）和（b）可以看出系統(tǒng)（1）的時域波形具有這種特點.而從圖3的（c）和（d）可以看出，他們的頻譜存在連續(xù)寬頻帶特性，沒有明顯的波峰，并且峰值連續(xù)，說明系統(tǒng)（1）具有混沌特性，

利用Poincare截面圖進一步分析系統(tǒng)（l），給定系統(tǒng)（1）中的參數(shù)后，選擇既不包含系統(tǒng)的軌跡，也不與軌線相切的平面作為Poincare截面，通過觀察截面上截點的情況，判斷系統(tǒng)是否可產(chǎn)生混沌運動.如圖4所示，得到系統(tǒng)（1）在幾個截面上的Poin-care映像，可見，在Poincare截面上有無窮多個分形結(jié)構(gòu)的密集點，形成一段連續(xù)的曲線，進一步說明了此時系統(tǒng)的運動是混沌的.

2.5 Lyapunov指數(shù)譜、分岔圖

如系統(tǒng)（1）參數(shù)改變，系統(tǒng)平衡點的穩(wěn)定性將發(fā)生變化，其運行狀態(tài)也發(fā)生相應的改變.隨參數(shù)變化的Lyapunov指數(shù)譜和分岔圖可以直觀地分析出系統(tǒng)狀態(tài)變化情況.以系統(tǒng)（1）中的部分參數(shù)變化的情況為例進行討論.

1）參數(shù)k變化情況：其他參數(shù)不變，改變k，

令參數(shù)k在[O，5.9]范圍內(nèi)變化，圖5的（a）和（b）給出了隨k變化時的Lyapunov指數(shù)譜和分岔圖.可以看出二者具有很好的一致性，當k在[0，1.1）時，所有LE均小于零，系統(tǒng)（1）處于穩(wěn)定狀態(tài).當k在[1.1，1.95），且除去1.2附近的值時，LE1大于零，其他LE2等于零，系統(tǒng)處于混沌狀態(tài).當k取1.2附近的值時，系統(tǒng)處于周期狀態(tài)，當k在[1. 95，2.9）時，系統(tǒng)又處于周期分岔狀態(tài)，當k在[2.9，5.9）時，系統(tǒng)又處于混沌狀態(tài)，

2）參數(shù)d變化情況：其他參數(shù)不變，改變d.

令參數(shù)d在[0，1000]范圍內(nèi)變化，圖5的（c）和（d）給出了隨著d變化時的Lyapunov指數(shù)譜和分岔圖，當d在[0，8]時，系統(tǒng)（1）處于穩(wěn)定狀態(tài)；當d在（8，20]時，系統(tǒng)出現(xiàn)倍周期分岔；當d在（20，50]時，系統(tǒng)處于混沌狀態(tài)；當d在（50，91]時，系統(tǒng)出現(xiàn)倍周期分岔；當d在（91，1000]時，系統(tǒng)又處于混沌狀態(tài)，同時在整個混沌帶內(nèi)存在著數(shù)個周期窗口.因此，系統(tǒng)（1）當參數(shù)d在[0，1000]內(nèi)變化時，LE可得到較大值，最大LE可達到5，而且相比其他系統(tǒng)呈現(xiàn)混沌的參數(shù)范圍較大.

3 系統(tǒng)的離散化仿真及FPGA實現(xiàn)

3.1 混沌系統(tǒng)的Modelsim仿真

為了采用數(shù)字電路實現(xiàn)新五維混沌系統(tǒng)，對該系統(tǒng)模型進行離散化，得到VHDL語言程序文件.利用Test Bench生成.tcl文件用于Modelsim進行RTL門級仿真，系統(tǒng)的xl、x2、x3、x4和x5變量的Modelsim仿真波形如圖6所示.可見，該離散化仿真結(jié)果與圖3中時域波形的Matlab仿真結(jié)果完全一致，說明新五維混沌系統(tǒng)離散模型正確，并可以在FPGA中實現(xiàn).

3.2 系統(tǒng)的FPGA實現(xiàn)

用Modelsim進行功能仿真后，將VHDL語言程序配置到FPGA中，本文選用型號為EP3C25E144C8的Cyclone系列FPCA構(gòu)建系統(tǒng)，以驗證混沌吸引子的存在，通過高速數(shù)模轉(zhuǎn)換芯片DAC904E，利用示波器觀察到模擬混沌吸引子相圖.為了和數(shù)值仿真結(jié)果做比較，本文在圖7中給出了五維混沌系統(tǒng)的部分吸引子相圖，這些相圖分別對應于圖1中給出的數(shù)值仿真相圖，可見，通過示波器觀測到的相軌跡圖同數(shù)值仿真分析是一致的，從物理意義上進一步驗證了新五維混沌系統(tǒng)的混沌特性.

4 結(jié) 論

本文提出一個新的五維混沌系統(tǒng)，對其進行了數(shù)值仿真，分析了基本的混沌動力學特性，定性地確定了系統(tǒng)混沌的存在性.并且該系統(tǒng)呈現(xiàn)混沌的d參數(shù)范圍很大，當該參數(shù)變化時，系統(tǒng)的最大Lya-punov指數(shù)可達到5，存在復雜的混沌動力學行為.最后，設(shè)計實現(xiàn)了該系統(tǒng)的硬件數(shù)字電路，在示波器上觀測到的吸引子相圖同數(shù)值仿真分析一致，從物理意義上進一步驗證了五維系統(tǒng)的混沌特性.提出的五維混沌系統(tǒng)及實現(xiàn)的電路可為隨參數(shù)變化的混沌信號發(fā)生器提供依據(jù)，也可為提高混沌保密通信安全性、混沌數(shù)值通信的可靠性及其他應用提供新的信號源.