概周期驅(qū)動二維正弦系統(tǒng)的奇異非混沌吸引子特性分析

摘要：研究了一類概周期驅(qū)動二維正弦系統(tǒng)奇異非混沌吸引子的存在性與度量特性。首先通過相圖描述了奇異非混沌吸引子的誕生過程，在給定參數(shù)下，奇異非混沌吸引子通過光滑環(huán)面倍化中斷而出現(xiàn)分形現(xiàn)象；其次，通過最大李雅普諾夫指數(shù)和相敏感指數(shù)證實了該奇異非混沌吸引子；最后利用功率譜、有限李雅普諾夫指數(shù)分布和遞歸圖表征了奇異非混沌吸引子特性。實驗結(jié)果表明，系統(tǒng)中存在奇異非混沌吸引子，并且具有較好的統(tǒng)計學特性。

關(guān)鍵詞：奇異非混沌吸引子；相敏感指數(shù)；有限李雅普諾夫指數(shù)分布

中圖分類號： O415.6;N941.7文獻標識碼： A

收稿日期：2023-05-11；修回日期：2023-09-19

基金項目：國家自然科學基金重點項目（11732014）；山東省自然科學基金面上項目（ZR2021MA095）

第一作者：張玉璽（2001-），男，山東菏澤人，碩士研究生，主要研究方向為非線性動力學。

通信作者：孫小淇（1986-），女，山東青島人，博士，副教授，主要研究方向為非線性微分方程及神經(jīng)網(wǎng)絡(luò)的動力學行為。

Analysis of Strange Nonchaotic Attractors of Periodic Two-dimensional Sinusoidal Systems

ZHANG Yuxi¹， SUN Xiaoqi²

（School of Mathematics and Statistics， Qingdao University， Qingdao 266071， China）

Abstract：The existence and metric characteristics of a class of quasiperiodic ally driven two-dimensional systems with sinusoidal functions are studied. Firstly， the birth process of strange nonchaotic attractors is described by phase diagram. Strange nonchaotic attractors are fractal due to the interruption of smooth torus doubling for the fixed parameters. Secondly， the strange nonchaotic attractors is confirmed by the largest Lyapunov exponents and phase sensitivity exponents. Finally， the power spectrum， finite Lyapunov exponents distribution and recurrence plots are used to characterize the strange nonchaotic attractors. According to the experimental results， there are strange nonchaotic attractors in the system and have good statistical characteristics.

Keywords： strange nonchaotic attractors; phase sensitivity exponent; the distribution of the finite-time Lyapunov exponents

0 引言

奇異非混沌吸引子幾何上具有分形特性（奇異性），但它的最大李雅普諾夫指數(shù)為負，體現(xiàn)了其非混沌特性。1984年，Grebogi在模擬非線性振蕩器的映射中首次發(fā)現(xiàn)并提出了奇異非混沌吸引子的概念^[1]。此后，研究奇異非混沌吸引子成為國內(nèi)外非線性動力學的重要課題之一，主要包括理論研究^[2-5]、數(shù)值研究^[6-7]和實驗研究^[8-11]，奇異非混沌吸引子也具有重要的應(yīng)用背景，如安全通信^[12-14]、電子電路^[15-18]、氣候動力學^[19]。

目前，研究奇異非混沌吸引子主要集中在概周期驅(qū)動的光滑動力學系統(tǒng)，且大量工作致力于研究其動力學轉(zhuǎn)遷和機制^[20-22]，如分形環(huán)面機制^[20]，間歇I型路徑^[21]、H-H路線^[22]等。近年來還有許多路線不斷被發(fā)現(xiàn)，如III型間歇路徑^[23]、鎖相路徑^[24]、擦邊分岔路徑^[25]。張在Wada域上描述了奇異非混沌吸引子的演化及分形性質(zhì)，通過不同的路徑描述了奇異非混沌吸引子的生成和演化，如環(huán)面分形路線、間歇-I路線、邊界路線和H-H路線^[26]、分岔誕生路徑^[27]。李^[28]通過理論和數(shù)值的方法證明了連續(xù)和離散系統(tǒng)中奇異非混沌吸引子的存在性，并進行了相關(guān)的數(shù)值分析和對多穩(wěn)態(tài)現(xiàn)象的控制。徐^[29]驗證了概周期驅(qū)動下二維分段線性范式系統(tǒng)中奇異非混沌吸引子的存在性并且用相軌跡圖、李雅普諾夫指數(shù)等方法進行了數(shù)值分析。在非光滑系統(tǒng)中奇異非混沌吸引子的產(chǎn)生機制和分形特性目前也引起廣泛關(guān)注^[30-32]。沈^[31]通過對概周期驅(qū)動的分段Logistic非光滑系統(tǒng)的研究，發(fā)現(xiàn)了存在著類型豐富的奇異非混沌吸引子，隨后，沈等^[32]研究了非光滑系統(tǒng)中環(huán)面邊界碰撞與奇異非混沌吸引子生成的機理，并通過不同路徑和機制識別了多種奇異非混沌吸引子。本文研究了一類概周期驅(qū)動二維正弦系統(tǒng)，探究了經(jīng)過環(huán)面倍化中斷后誕生奇異非混沌吸引子的過程，首先利用相圖來觀測奇異非混沌吸引子的產(chǎn)生過程，并展示如何用最大李雅普諾夫指數(shù)和相敏感指數(shù)來識別奇異非混沌吸引子，并通過功率譜、遞歸圖、譜分布函數(shù)、有限李雅普諾夫指數(shù)分布來進一步表征奇異非混沌吸引子。

1 系統(tǒng)模型

Lai等^[33]通過在Sprott B系統(tǒng)中引入正弦函數(shù)，得到了一個具有無限多個共存吸引子的混沌系統(tǒng)。根據(jù)Lai等的想法，Li等^[34]將正弦函數(shù)引入一類二維映射中，以生成無限多個共存吸引子，二維正弦系統(tǒng)描述如式（1）：

方程（1）中引入概周期驅(qū)動力后變?yōu)?/p>

其中，ɑ為控制參數(shù)，β為概周期驅(qū)動的振幅，無理數(shù)ω為系統(tǒng)頻率，ω=（ /5-1）/2。

為了研究奇異非混沌吸引子，首先利用相圖觀察奇異非混沌吸引子的產(chǎn)生過程，再利用最大李雅普諾夫指數(shù)和相敏感指數(shù)檢驗奇異非混沌吸引子。不妨取概周期驅(qū)動的振幅β=0.100，取參數(shù)b=1，c=0.1。

本節(jié)用一個典型的例子來證實奇異非混沌吸引子的產(chǎn)生，圖1是在取定振幅β的情況下，在α=0.456附近的吸引子相圖，表明從概周期吸引子到奇異非混沌吸引子的躍遷。能夠看出在概周期強迫下，不同周期為k的周期軌道變成了kT環(huán)面的準周期軌道，圖1給出了不同參數(shù)α下的相平面圖，當α=0.45時，系統(tǒng)的相圖表現(xiàn)為1T概周期環(huán)面，隨著α的增加到0.454時，系統(tǒng)演化成了2T概周期環(huán)面，當α增加到α=0.456前后時，吸引子表現(xiàn)出類似混沌吸引子的分形特性，但通過式（3）計算它的最大李雅普諾夫指數(shù)約為-0.008（小于零），因此，它是一個非混沌吸引子。當α增加到0.458時，吸引子的最大李雅普諾夫指數(shù)約為0.004（大于零），吸引子變?yōu)榛煦缥印?/p>

為了進一步證實并描述圖1c呈現(xiàn)的奇異非混沌吸引子的特性，本文通過相敏感函數(shù)與相敏感指數(shù)、功率譜與譜分布函數(shù)、遞歸圖、有限時間李雅普諾夫指數(shù)分布4種方法分析其特性。

2 奇異非混沌吸引子特性分析

2.1 相敏感函數(shù)與相敏感指數(shù)

相敏感函數(shù)是判斷概周期吸引子和奇異非混沌吸引子的重要方法之一^[17]，奇異非混沌吸引子常有如下性質(zhì)：Γ_N～N^μ，其中μ為相敏感指數(shù)，即相敏感函數(shù)Г_N隨著軌道Ν的長度而增長，并表現(xiàn)出冪律關(guān)系。相敏感函數(shù)計算方式為

通過計算可得圖1c吸引子的相敏感函數(shù)，由圖2可以看出相敏感函數(shù)Г_N與Ν存在冪律關(guān)系，即Г_N～N^μ，通過計算得：相敏感指數(shù)α=0.456時，μ約為0.994，該吸引子明顯區(qū)別于概周期吸引子，又因其最大李雅普諾夫指數(shù)約為-0.008，因此，該吸引子為奇異非混沌吸引子。

2.2 功率譜與譜分布函數(shù)

功率譜也是驗證奇異非混沌吸引子的一個重要方法，功率譜特征通過用時間序列分析非線性動力系統(tǒng)的波動狀態(tài)來表示，對于周期吸引子，功率譜是由基頻諧波處σ-峰組成，混沌吸引子為連續(xù)譜，奇異非混沌吸引子的功率譜具有獨特的特征，即奇異連續(xù)性。

由傅里葉變換

得到功率譜

其中，X（ω，T）為傅里葉級數(shù)，ω∈[0，1]為功率譜頻率，n為傅里葉變換的迭代次數(shù)，x_n為系統(tǒng)第n次迭代時的輸入變量。圖3a顯示了在參數(shù)α=0.456時的功率譜，對于奇異非混沌吸引子，功率譜有分形特性并且具有自相似的峰。

譜分布函數(shù)也是驗證和表征奇異非混沌吸引子的重要手段，譜分布函數(shù)是通過計算功率譜中大于σ值的峰的個數(shù)而得出的，與周期吸引子和混沌吸引子不同，奇異非混沌吸引子滿足冪律關(guān)系G（σ）～σ^-β，1lt;βlt;2，β為擬合的斜率。得到了以lgσ為橫坐標，lgG（σ）為縱坐標的擬合圖，如圖3b所示，可以看到它們的關(guān)系可以擬合成一條直線，并且斜率約等于-1.92。

2.3 遞歸圖

通過引入遞歸圖，以可視化的時間，將二維動力系統(tǒng)的狀態(tài)顯現(xiàn)出來，通過遞歸矩陣上的不同結(jié)構(gòu)，能夠區(qū)分基礎(chǔ)軌跡所蘊含的各種運動狀態(tài)。圖4a顯示了參數(shù)α=0.456下的隊規(guī)矩陣的可視化，顯示了重復模式也沒有中斷的均勻紋理，從圖4b能夠看出，固定時間段的垂直和水平波段代表了兩組不同參數(shù)下不同時期狀態(tài)的遞歸，說明底層過程實際上還是復雜的、具有確定性的。但是最大李雅普諾夫指數(shù)小于0，這說明吸引子在圖1c參數(shù)下是奇異非混沌吸引子。

2.4 有限時間Lyapunov指數(shù)分布

由于不同種類的吸引子在有限時間李雅普諾夫指數(shù)分布上會表現(xiàn)出不同的特征，因此對于檢驗奇異非混沌吸引子是有幫助的。有限時間李雅普諾夫指數(shù)P（t，λ）的分布具有一個特殊性質(zhì)，奇異非混沌吸引子在有限時間間隔內(nèi)通常具有正的李雅普諾夫指數(shù)，但漸近指數(shù)為負。對于奇異非混沌吸引子，在有限時間李雅普諾夫指數(shù)大于零的區(qū)域會出現(xiàn)正值曲線。對于環(huán)面吸引子，沒有尾部延伸到局部李雅普諾夫指數(shù)大于零的區(qū)域。這里不妨取t=50，并計算有限時間李雅普諾夫的分布。圖5為奇異非混沌吸引子在ɑ=0.456，β=0.100時的有限時間李雅普諾夫指數(shù)分布圖，它具有一個顯著的特點是有限時間李雅普諾夫指數(shù)分布中的正尾呈指數(shù)（縱坐標為對數(shù)）衰減，在其負值中表現(xiàn)出細長的尾巴，這種特性不同于其它類型吸引子（周期吸引子、概周期吸引子和混沌吸引子）。

3 總結(jié)

本文通過一類簡單概周期驅(qū)動二維正弦映射，通過多種方法驗證了系統(tǒng)中奇異非混沌吸引子的存在性，并描述了其分形特性和統(tǒng)計學特征。首先通過相圖，觀察到了奇異非混沌吸引子的形成過程，并通過最大李雅普諾夫指數(shù)和相敏感指數(shù)識別奇異非混沌吸引子，通過功率譜、遞歸圖、譜分布函數(shù)和有限時間李雅普諾夫指數(shù)對奇異非混沌吸引子進行進一步描述。研究結(jié)果表明，系統(tǒng)中存在奇異非混沌吸引子。通過大量實驗數(shù)據(jù)分析，此系統(tǒng)在其他參數(shù)區(qū)域也存在相似的奇異非混沌吸引子，關(guān)于奇異非混沌吸引子的實驗還有待進一步研究。

參考文獻：

[1]GREBOGI C， OTT E， PELIKAN S，et al. Strange attractors that are not chaotic[J]. Physica D， 1984， 13（1/2）：261-268.

[2]JAEGER T H.The creation of strange non-chaotic attractors in non-smooth saddle-node bifurcations[J]. Mem Amer Math Soc， 2009， 201（945）：1-106.

[3]THAMILMARAN K， SENTHILKUMAR D V， VENKATESAN A， et al. Experimental realization of strange nonchaotic attractors in a quasiperiodically-forced electronic circuit[J]. Phys Rev E， 2006， 74（3）：036205.

[4]ZHOU T， MOSS F， BULSARA A. Observation of a strange nonchaotic attractor in a multistable potential[J]. Phys Rev A， 1992， 45（8）：5394-5400.

[5]GRAGER M，JAEGER T H. Dimensions of attractors in pinched skew products[J]. Commun Math Phys， 2013， 320（1）：101-119.

[6]DING M， GREBOGI C， OTT E. Dimensions of strange nonchaotic attractors[J]. Phys Lett A， 1989， 137（4/5）：167-172.

[7]MITSUI T， UENOHARA S， MORIE T， et al. Torus-doubling process via strange nonchaotic attractors[J]. Physics Letters A， 2012， 376（24）：1907-1914.

[8]DITTO W L， SPANO M L， SAVAGE H T，et al. Experimental observation of a strange nonchaotic attractor[J]. Phys Rev Lett， 1990， 65（5）：533-536.

[9]ZHOU T， MOSS F， BULSARA A. Observation of a strange nonchaotic attractor in a multistable potential[J]. Phys Rev A， 1992， 45（8）：5394-5400.

[10] KETOJA JA， SATIJA II. Harper equation， the dissipative standard map and strange nonchaotic attractors： relationship between an eigenvalue problem and iterated maps[J]. Physica D， 1997， 109（1）：70-80.

[11] 劉劍波，葉春飛，張樹京.基于收縮映射的奇異非混沌系統(tǒng)同步[J].物理學報，2000，49（1）：20-23.

LIU J B， YE C F， ZHANG S J. Synchronization of strange nonchaotic systems basedon contraction mapping[J]. Acta Physica Sinica， 2000， 49（1）：20-23.

[12] ZHOU C S， CHEN T L. Robust communication via synchronization between nonchaotic strange attractors[J]. Europhys Lett， 1997， 38（4）：261-265.

[13] RAMASWAMY R. Synchronization of strange nonchaotic attractors[J]. Phys Rev E， 1997， 56：7294-7296.

[14] CHACON R， GRACIA-HOZ A M. Route to chaos via strange non-chaotic attractors by reshaping periodic excitations[J]. Europhys Lett， 2002， 57（1）：7-13.

[15] LIU Z， ZHU Z W. Strange nonchaotic attractors from quasiperiodically forced Ueda’s circuit[J]. Int J Bifurcation Chaos Appl Sci Eng， 1996， 6（7）：1383-1388.

[16] VENKATESAN A， MURALI K， LAKSHMANAN M. Birth of strange nonchaotic attractors through type III intermittency[J]. Phys Lett A， 1999， 259（3）：246-253.

[17] 張永祥，俞建寧，褚衍東，等.一類新電路系統(tǒng)的奇怪非混沌吸引子分析[J].河北師范大學學報（自然科學版），2008，32（6）：5.

ZHANG Y X， YU J N， CHU Y D， et al. Analysis of strange nonchaotic attractors for a class of new circuit systems[J]. Journal of Hebei Normal University （Natural Science）， 2008， 32（6）：5.

[18] 謝帆，楊汝，張波.電流反饋型Buck變換器二維分段光滑系統(tǒng)邊界碰撞和分岔研究[J].物理學報，2010，59（12）：8393-8406.

XIE F， YANG R， ZHANG B. Research on boundary collisions and bifurcations in a two dimensional segmental smooth system of current feedback Buck converter[J]. Acta Physica Sinica， 2010， 59（12）：8393-8406.

[19] MITSUI T， AIHARA K. Dynamics between order and chaos in conceptual models of glacial cycles[J]. Clim Dynam， 2014， 42（11）：3087-3099.

[20] NISHIKAWA T， KANEKO K. Fractalization of a torus as a strange nonchaotic attractor[J]. Physical Review E， 1996， 54（6）：6114-6124.

[21] PRASAD A， MEHRA V， RAMASWAMY R. Intermittency route to strange nonchaotic attractors[J]. Physical Review Letters， 1997， 79（21）：4127-4130.

[22] HEAGY J F， HAMMEL S M. The birth of strange nonchaotic attractors[J]. Physica D： Nonlinear Phenomena， 1994， 70（1/2）：140-153.

[23] VENKATESAN A， MURALI K， LAKSJMANAN M. Birth of strange nonchaotic attractors through type III intermittency[J]. Physics Letters A， 1999， 259（3）：246-253.

[24] ALSTR?M P， CHRISTIANSEN B， LEVINSEN M T. Nonchaotic transition from quasiperiodicity to complete phase locking[J]. Physical Review Letters， 1988， 61（15）：1679-1682.

[25] WEN G L， YIN S， XU H D， et al. Analysis of grazing bifurcation from periodic motion to quasi-periodic motion in impact-damper systems[J]. Chaos， Solitons amp; Fractals， 2016， 83：112-118.

[26] ZHANG Y. Strange nonchaotic attractors with Wada basins[J]. Physica D： Nonlinear Phenomena， 2013， 259：26-36.

[27] PIKOVSKYA， FEUDEL U. Characterizing strange nonchaotic attractors[J]. Chaos， 1995， 5（1）：253-260.

[28] 李高磊.動力系統(tǒng)中的奇異非混沌吸引子及多穩(wěn)態(tài)動力學研究[D].成都：西南交通大學， 2022.

LI G L. Study on Strange NonChaotic Attractors and Multistable Dynamics in Dynamic Systems[D]. Chengdu： Southwest Jiaotong University， 2022.

[29] 徐震，沈云柱.概周期驅(qū)動的二維分段線性范式系統(tǒng)的奇異非混沌吸引子[J].濟南大學學報（自然科學版），2022，36（2）：231-236.

XU Z， SHEN Y Z. Strange nonchaotic attractors for almost periodic driven two dimensional piecewise linear normal form systems[J]. Journal of University of Jinan（Science and Technology）， 2022， 36（2）：231-236.

[30] LI G L， YUE Y， XIE J H， et al. Strange nonchaotic attractors in a nonsmooth dynamical system[J]. Communications in Nonlinear Science and Numerical Simulation， 2019. 78：104858.

[31] 沈云柱，張凡輝，東廣霞.概周期驅(qū)動分段Logistic系統(tǒng)的奇異非混沌吸引子[J].河北師范大學學報：自然科學版，2019，43（3）：6.

SHEN Y Z， ZHANG F H， DONG G X. Strange nonchaotic attractors for almost periodic driven piecewise logistic systems[J]. Journal of Hebei Normal University （Natural Science）， 2019， 43（3）：6.

[32] SHEN Y， ZHANG Y. Mechanisms of strange nonchaotic attractors in a nonsmooth system with border-collision bifurcations[J]. Nonlinear Dynamics， 2019， 96（2）：1405-1428.

[33] LAI Q， KUATE P D K， LIU F， et al. An extremely simple chaotic system with infinitely many coexisting attractors[J]. IEEE Trans Circuits Syst II， 2019， 67（6）：1129-1133.

[34] ZHANG L，LIU Y， WEI Z C， et al. A novel class of two-dimensional chaotic maps with infinitely many coexisting attractors[J]. Chinese Physics B， 2020， 29（6）：060501.

（責任編輯 耿金花）