含多項式取絕對值函數的混沌系統分析與應用

高正中　杜翔

摘要： 為獲得更復雜的動力學特性，設計了一種含有三次多項式取絕對值函數的混沌系統，該混沌系統的理論模型由含有3個狀態變量的非線性方程組描述。分析了該系統的基本屬性，以及相軌圖、時域波形、Lyapunov指數、Poincaré映射、分岔圖等動力學特性。系統在一定的參數條件下存在周期和混沌的特性，初值對稱時存在共存吸引子或聚合吸引子。此外，某些系統參數變化時，系統存在恒定的動力學特性，狀態變量初值變化時系統的動力學特性也是恒定的。電路仿真驗證了理論的正確性?；谛略O計的混沌系統設計了一種加密方案，對加密性能進行了分析，表明了加密方案的有效性。

關鍵詞： 混沌；動力學；共存吸引子；共存分岔；圖像加密

中圖分類號： TP271.62?????? 文獻標識碼： A

Analysis and Application of Chaotic System with Polynomial Absolute Valued Function

GAO Zhengzhong，DU Xiang

（College of Electrical Engineering and Automation， Shandong University of Science and Technology， Qingdao 266510， China）

Abstract：In order to get more complex dynamical properties， a chaotic system with cubic polynomial taking absolute value function is designed. The theoretical model of the chaotic system is described by a set of nonlinear equations with three state variables. The basic properties and the dynamic characteristics of the system like phase diagram， time domain waveform， Lyapunov exponent， Poincaré map and bifurcation diagram are analyzed. Under certain parameter conditions， the system has periodic and chaotic properties， and there are coexistence attractors or aggregation attractors when the initial values are symmetrical. In addition， when some system parameters change， the system has constant dynamic characteristics， when the initial value of the state variable changes， the dynamic characteristics of the system also remain unchanged. The correctness of the theory is verified by circuit simulation. Based on the newly designed chaotic system， an encryption scheme is designed， and the encryption performance is analyzed， which shows the effectiveness of the encryption scheme.

Keywords： chaos; dynamics; coexisting attractors; coexisting bifurcation; image encryption

0 引言

1963年，美國氣象學家Lorenz在研究天氣預報時，把大氣動力學方程組進行簡化并在計算機上進行模擬實驗，發現舍去小數點后幾位數字的模型與未舍去的模型迥然不同[1]，于是他得出了“短期的天氣預報可行，但長期不可行”的結論。Lorenz研究天氣預報時提出了“蝴蝶效應”與混沌理論。這一理論表明確定性系統可以產生復雜的、類似隨機的運動。自此，諸多學者紛紛開展了大量的研究?；煦缋碚撘褟V泛應用于圖像加密[24]、保密通信[57]和神經網絡[810]等領域。師東生等[11]構建了一個在兩個方向上具有四翼特征的五維超混沌系統。鮮永菊等[12]構造了一個四維超混沌系統，該系統有豐富的多穩態特性與多翼吸引子，最多有12種吸引子共存的情況。Yang等[13]從一個三維混沌系統出發，加入3個線性動態控制器得到了六維混沌系統，該系統有隱藏超混沌等復雜現象。Ma等[14]在三維自治混沌系統中加入平滑磁通控制的憶阻器和交叉項，設計了一種可以產生四翼混沌吸引子的新型四維超混沌系統。Qi等[15]從Lorenz系統出發，提出了一種每個狀態方程都含有一個二次交叉項的新型混沌系統，該混沌系統有不動點、周期環、混沌等多種吸引子。Atiyeh等[16]提出了一個具有平衡面的四維混沌系統。Bao等[17]在亞遺傳混沌Jerk電路里引入了一個憶阻器，得到了一個新型混沌電路，在系統狀態變量初值不同時，系統有完全不同的動力學行為。Chen等[18]利用切換函數生成了環形、球形混沌吸引子。臧鴻雁等[19]給出了一個分段線性混沌映射，并分析了其與一個三次多項式映射拓撲共軛的充分條件。高秉建[20]通過增加第二類鞍焦平衡點的數目提出了一種新的生成多翼混沌吸引子的方法。

目前，含有多項式取絕對值的非線性函數并能產生共存吸引子且無憶阻器的混沌系統鮮有報道。不含有多穩態現象的混沌系統，動力學特性的復雜程度較低，實際應用時效果稍差，雖然憶阻系統可以產生共存吸引子，但憶阻器沒有完全商用化，制備成本高，因此這樣的混沌系統在物理實現時，成本自然很高。本文設計了一個表達式中含有三次多項式取絕對值的混沌系統，對系統進行了基本特性考察、動力學特性分析，發現系統存在共存吸引子或聚合吸引子，這降低了含共存吸引子的混沌系統的物理實現成本。并且進行了電路仿真，最后基于該系統設計了一種加密方案并驗證了其有效性。

1 混沌系統描述

本文設計的混沌系統表達式：

其中，狀態變量由x，y，z表示，系統參數由a，b，c，d，e，f，g表示。所有系統參數均大于0。

2 基本特性與動力學行為分析

2.1 系統相圖、時域波形、Poincaré映射圖

在系統狀態方程中，若取a=4，b=5，c=4，d=0.1，e=1，f=1，g=7，并且系統初始條件設定為（x0，y0，z0）=（-1，-1，-1）時，系統產生的吸引子如圖1所示。分析圖1可知：系統的運動軌跡是雜亂無章的，同時可以發現吸引子具有渦卷結構，其吸引子屬于混沌吸引子。

狀態變量x，y的時域波形如圖2 a、b所示，波形呈現類似隨機的信號；系統的Poincaré映射取截平面z=0時，如圖2 c所示，呈現連續、成片的密集點。計算系統的Lyapunov指數為LE1=0.242 7、LE2=-0.044 2、LE3=-1.848 0。Lyapunov之和為-1.649 5，該值小于0，保證了系統的穩定性；并且得出Lyapunov維數DL=2.107 4，大于2而小于3。再結合相軌跡圖等仿真圖可以判斷出系統處于混沌狀態。

2.2 對稱性、耗散性、平衡點與穩定性

對系統方程組做（x，y，z）→（-x，-y，-z）變換，在系統參數不變的前提下，方程組的解不變，于是可以得到結論：系統是關于原點對稱的。

顯然，式（4）總是有一個負實根λ=-e。其它兩個根及平衡點的性質取決于系統參數a，b，g的取值情況。在所有系統參數均大于0的前提下，若b2-4ag≥0，方程的其它兩個根均為正實數，平衡點為指標2的鞍點；若b2-4ag<0則為指標2的鞍焦點。系統參數按照前面取值時，平衡點處的特征值λ1=-1，λ2=2.5+4.663 7j，λ3=2.5-4.663 7j，于是平衡點為指標2的鞍焦點。

2.3 參數g對動力學行為的影響

本節令參數g在[4， 8]這一區間上變化，初始條件?。▁0，y0，z0）=（1，1，1）。圖3 a，b分別給出了z的分岔圖和Lyapunov指數譜。其中，第3個Lyapunov指數是一個絕對值相對較大的負數，為了更好地呈現其它兩個Lyapunov指數，本節與2.4節將其省略不在圖中呈現。在g∈[4，8]區間，周期與混沌態發生多次切換。具體為：g∈[4，4.78]∪[5.04，5.38]∪[6.32，6.88]∪[7.67，8]時，系統處于周期狀態；g∈[4.78，5.04]∪[5.38，6.32]∪[6.88，7.67]時，系統處于混沌狀態。

由圖3b分岔圖知，系統經歷了多次正向、反向的分岔。其中，當g∈[5.38，6.32]時，分岔圖顯示了一些狹窄的周期窗口。選取不同的參數g，系統在x-z平面上的相圖如圖4所示。

2.4 參數f對動力學特性的影響

與2.3節類似，令參數f在[2，14]區間上變化，進行分析。首先，從圖5a和b可以看出，首個Lyapunov指數一直大于0，第二根Lyapunov指數一直等于0，相圖中則總是有成片的密集點，表明系統一直處于混沌態。然后，僅改變一個系統參數，令g=5.2，圖5c和d重新給出了上述兩幅仿真圖，可以看出首個Lyapunov指數一直約等于0，第二根Lyapunov指數一直小于0，分岔圖中點分布得比較稀疏，意味著系統一直處于周期態，因此系統參數f變化時，根據其它系統參數取值的情況，系統具有持續混沌或持續周期的特性。

2.5 參數變化時持續動力學特性理論分析

由2.3和2.4的分析知道，系統參數g變化能影響系統的動力學行為，而系統參數f變化不影響。將系統的平衡點代入特征方程，參數c，d，f會完全消去。具體來說，當x=y=z=0時，特征方程的系數a2=e-b，a1=ag-be，a0=aeg，即特征方程的系數與參數c，d，f無關。從理論上來講，當這3個參數改變時，系統平衡點類型相同，擁有相同的動力學特性。

另外，由圖5b，d知，當系統參數f改變時，分岔一直在下降，且分岔圖或其上方輪廓類似反比例函數曲線，因此該系統參數改變時系統具有調幅特性。從理論上來說，系統參數c，d改變時，系統同樣具有調幅特性[2122]。

2.6 狀態變量初值對動力學行為的影響

本節中，令系統參數按照3.1取值，每次固定兩個狀態變量的初值等于1，第3個狀態變量初值在[-30，30]之間變化，從圖6所示的Lyapunov指數譜可以發現，3個狀態變量初值在[-30，30]這一區間上變化時，Lyapunov指數譜符號總為（+，0，-），因此，在一定的參數條件下，系統初值變化時存在持續混沌的特性。

2.7 系統的共存吸引子與共存分岔

共存吸引子、共存分岔也是混沌系統的重要特征。共存吸引子是指系統參數取相同數值而初值不相等的時候，所產生的不一樣的吸引子。其它系統參數按照3.1取值，改變g的數值，共存吸引子在x-y平面上的投影如圖7a～e所示。其中，紅色的軌跡為系統初值?。?1，-1，-2）時，系統的吸引子投影，藍色的軌線則對應系統初值?。?，1，2）。x的共存分岔圖如圖7 f所示，紅色的部分對應初值?。?1，-1，-2），藍色對應系統初值?。?，1，2）。由圖7各圖可見，系統初值取對稱值時，g∈[5.38，6.32]系統產生聚合混沌吸引子，g∈[4，5.38]∪[6.32，8]時產生對稱共存吸引子。

3 系統的模擬電路仿真

模擬電路仿真通常在Multisim軟件平臺進行。首先，運算放大器的線性范圍為±13.5V，模擬乘法器線性范圍為±10V，因此需要對超過上述范圍的狀態變量z進行壓縮變換。設z→2z，變換后，式（1）轉化為

在電路的狀態方程中，首先取電容C1=C2=C3= 10nF，R4=R5=R8=R10=R11=R12=R13=R15=R16=R17=R19 =10kΩ，對比系統方程和電路狀態方程，得

1R1C1=1 000，1R2C1=2 000，1R3C2=8 000，1R6C2=5 000，1R18C2=4 000，1100R14C2=100，1R20C3=3 500

進而可求得R1=100kΩ，R2=50kΩ，R3=12.5kΩ，R6=20kΩ，R14=10kΩ，R18=25kΩ，R20=28.571kΩ。取R20=28.7kΩ，其它電路元件按照計算結果取值，令3個電容的初始電壓均取-1V，即設定系統初值為（x0，y0，z0）=（-1，-1，-1），y-x平面、y-z平面的仿真實驗結果如圖9所示，證明了理論的正確性。

4 系統在圖像加密中的應用

本加密解密算法基于新設計的混沌系統，結合了像素置亂和擴散的算法。像素置亂可以使得像素點偏離原始位置，使相鄰像素相關程度下降；而擴散可以改變像素的灰度值，影響圖像的統計特性。

本加密解密算法分為4部分——一個混沌密碼發生器、兩個明文無關的擴散算法和一個明文關聯的置亂算法。加密解密過程如圖10所示。

加密詳細過程：

1）混沌密碼發生器。借助新混沌系統生成6個M×N維的偽隨機矩陣。記密鑰為K={x0，y0，z0，r1，r2}，具體步驟：

第1步：由輸入密鑰K的前3個值作為新混沌系統的初始值，共迭代r1+r2+MN次，其中前r1+r2次是為了跳過系統的暫態。得到的3個偽隨機序列記作{xi}、{yi}、{zi}，i=1，2，…，MN。

2）明文敏感性分析。該性能測試方法為：借助密鑰K先后加密兩幅圖像P1和P2，其中P2從P1中隨機選取一個像素點（i，j），將該像素點的值改變1得到。相應的密文圖像分別記為C1和C2。計算C1和C2之間的NPCR（像素改變率）、UACI（歸一化平均變化強度）值。上述方法共重復100次。最后計算上述2個指標的平均值。

仿真時，明文和密文圖像分別選取2 000對相鄰的像素點。由表1中的相關系數計算結果與圖13的相鄰像素灰度值示意圖可得出結論：明文圖像在各個方向都具有很強的相關性，而密文圖像則無相關性。

3）密文敏感性分析。該性能測試方法與明文敏感性類似，不同在于，改變的不是明文圖像的像素點，而是密文圖像的像素點。并且需要用同一密鑰對兩個不同的密文圖像解密分別得到P1和P2，計算P1和P2之間的NPCR、UACI值。對“Cameraman”加密后的圖像，計算得NPCR=99.603 1%，UACI=31.173 9%，結果接近理論值，因此該加密方案具有較好的密文敏感性。

4）密鑰靈敏度。正確的解密結果應如圖14 a所示。但進行解密操作時若取密鑰K={x0，y0，z0，r1，r2}={-1，-1，-1-10-15，35，201}，解密結果將如圖14 b所示。雖然密鑰僅有一個數值有10-15的偏差，但直觀上可以發現解密圖像仍然呈現噪聲形式，無法傳達有效信息。這表明該加密方案的密鑰靈敏度極優。

5 結論

本文設計了一種含有三次多項式取絕對值函數的混沌系統。通過數值分析和MATLAB仿真分析了基本屬性和動力學特性。結果表明，雖然這一系統在一定的參數條件下僅有一個指標2的鞍焦平衡點，但存在混沌現象，且狀態變量初值對稱時存在共存吸引子或聚合吸引子。系統的某幾個參數或系統初值變化時，具有持續周期或混沌的動力學特性，也有調幅特性。電路仿真驗證了理論的正確性。最后基于新設計的混沌系統設計了加密方案，仿真與分析表明了加密方案的有效性。

參考文獻：

[1]LORENZ E N. Deterministic non-periodic flows [J]. Journal of the Atmospheric Sciences， 1963， 20（3）： 130－141.

[2]田興華， 張紀會， 李陽. 基于混沌映射的自適應退火型粒子群算法 [J]. 復雜系統與復雜性科學， 2020， 17（1）： 45－54.

TIAN X H， ZHANG J H， LI Y. An adaptive annealing particle swarm optimization based on chaotic mapping [J]. ComplexSystems and Complexity Science， 2020， 17（1）： 45－54.

[3]LUO Y L， ZHOU R L， LIU J X， et al. A novel image encryption scheme based on Keplers third law and random Hadamard transform [J]. Chinese Physics B， 2017， 26（12）： 120504.

[4]SUN J W， YANG Q F， WANG Y F. Dynamical analysis of novel memristor chaotic system and DNA encryption application [J]. Iranian Journal of Science and Technology， Transactions of Electrical Engineering， 2020， 44： 449－460.

[5]劉林芳， 芮國勝， 張洋，等. 基于相空間對稱Lorenz陣子群的混沌通信保密研究 [J]. 通信學報， 2019， 40（5）： 32－38.

LIU L F， RUI G S， ZHANG Y， et al. Research on the chaotic secure communication of the phase-space symmetric Lorenz oscillator group [J]. Journal on Communications， 2019， 40（5）： 32－38.

[6]SAHIN M， TASJIRAN Z， HASAN G， et al. Application and modeling of a novel 4D memristive chaotic system for communication systems[J]. Circuits， Systems， and Signal Processing， 2020， 39： 3320－3349.

[7]YU F， LIU L， HE B Y， et al. Analysis and FPGA Realization of a novel 5D hyperchaotic four-wing memristive system， active control synchronization， and secure communication application [DB/OL].[2022-02-02].http：//doi.org/10.1155/2019/4047957.

[8]盧輝斌， 王麗佳. 基于Hopfield網絡的彩色圖像混沌加密算法 [J]. 吉林大學學報（信息科學版）， 2014， 32（2）： 131－137.

LU H B， WANG L J. Color image encryption algorithm of chaotic based on the hopfield network [J]. Journal of Jilin University （Information Science Edition）， 2014， 32（2）： 131－137.

[9]LIN H R， WANG C H， TAN Y M， et al. Hidden extreme multistability with hyperchaos and transient chaos in a Hopfield neural network affected by electromagnetic radiation [J]. Nonlinear Dynamics， 2020， 99： 2369－2386.

[10] 孫亮， 羅佳， 喬印虎.局部有源憶阻器電路及其在HR耦合神經元網絡中的應用[J].電子與信息學報，2021，43（11）： 3374－3383.

SUN L， LUO J， QIAO Y H. A Locally active memristor circuit and its application to a coupled hindmarsh-rose neuron network [J]. Journal of Electronics & Information Technology， 2021， 43（11）： 3374－3383.

[11] 師東生， 石煒. 一個新四翼高維超混沌系統的復雜動力學行為研究與仿真 [J]. 現代電子技術， 2020， 43（19）： 10－13.

SHI D S， SHI W. Research and simulation of complex dynamic behavior of new four-wing high-dimensional hyperchaotic system [J]. Modern Electronics Technique， 2020， 43（19）： 10－13.

[12] 鮮永菊， 扶坤榮， 徐昌彪. 一個具有多翼吸引子的四維多穩態超混沌系統 [J]. 振動與沖擊， 2020， 40（1）： 15－22， 38.

XIAN Y J， FU K R， XU C B. A 4-D multi-stable hyper-chaotic system with multi-wing attractors [J]. Journal of Vibration and Shock， 2020， 40 （1）： 15－22， 38.

[13] YANG L B， YANG Q G， CHEN G R. Hidden attractors， singularly degenerate heteroclinic orbits， multistability and physical realization of a new 6D hyperchaotic system [J]. Communications in Nonlinear Science and Numerical Simulation， 2020， 90： 105362.

[14] MA J， CHEN Z Q， WANG Z L， et al. A four-wing hyper-chaotic attractor generated from a 4-D memristive system with a line equilibrium [J]. Nonlinear Dynamics， 2015， 81： 1275－1288.

[15] QI G Y， CHEN G R， DU S Z， et al. Analysis of a new chaotic system [J]. Physica A， 2005， 352： 295－308.

[16] ATIYEH B， KARTHIKEYAN R， ABDUL J， et al. Dynamical analysis of a new multistable chaotic system with hidden attractor： antimonotonicity， coexisting multiple attractors， and offset boosting [J]. Physics Letters A， 2019， 383： 1450－1456.

[17] BAO H， WANG N， BAO B C， et al. Initial condition-dependent dynamics and transient period in memristor-based hypogenetic jerk system with four line equilibria [J]. Communications in Nonlinear Science and Numerical Simulation， 2018， 57： 264－275.

[18] CHEN Q F， HONG Y G， CHEN G R. Chaotic behaviors and toroidal/spherical attractors generated by discontinuous dynamics[J]. Physica A， 2006， 371： 293－302.

[19] 臧鴻雁， 韋心元， 袁悅. 一類三次多項式混沌映射的判定及性能分析 [J]. 電子與信息學報， 2021， 43（2）： 454－460.

ZANG H Y， WEI X Y， YUAN Y. Determination and properties analysis of a cubic polynomial chaotic map [J]. Journal of Electronics & Information Technology， 2021， 43（2）： 454－460.

[20] 高秉建. 基于Liu混沌系統生成的多翅膀蝴蝶吸引子 [J]. 復雜系統與復雜性科學， 2016， 13（1）： 91－94.

GAO B J. Multi-wing butterfly attractor from a modified chaotic system [J].Complex Systems and Complexity Science， 2016， 13（1）： 91－94.

[21] WANG Z， SUN C， JIN F， et al. A widely amplitude-adjustable chaotic oscillator based on a physical model of HP memristor [J]. IEICE Electronics Express， 2018， 15（8）： 1－6.

[22] WANG R， LI C B， CICEK S， et al. Amemristive hyperjerk chaotic system： amplitude control， FPGA design， and prediction with artificial neural network [DB/OL].[2022-02-01].http：//doi.org/10.1155/2021/6636813.

（責任編輯 李 進）

收稿日期：2022－05－18；修回日期：2022－07－08

基金項目： 中國博士后科學基金（2015T80279）；山東省自然科學基金（ZR2020MF071）

第一作者： 高正中（1971－），男，山東濟寧人，博士，教授，主要研究方向為計算機檢測與控制、機器人技術等。