一類分段非線性映射的混沌邊界分析

許宏飛，李群宏，寧 敏，商夢媛

（廣西大學 數(shù)學與信息科學學院，廣西 南寧 530004）

一類分段非線性映射的混沌邊界分析

許宏飛，李群宏*，寧 敏，商夢媛

（廣西大學 數(shù)學與信息科學學院，廣西 南寧 530004）

文章研究了一類具有非線性分支的分段映射的動力學行為.該模型可能應(yīng)用到物理科學、工程和醫(yī)學方面，也有助于一些經(jīng)濟模型的研究.以μ為分岔參數(shù)得到系統(tǒng)的分岔圖，發(fā)現(xiàn)在系統(tǒng)的不變吸引區(qū)間內(nèi)，周期軌道的每個周期點都有一定的存在范圍，這造成分岔結(jié)構(gòu)中出現(xiàn)迭代禁區(qū)現(xiàn)象.通過理論推導(dǎo)確定了周期軌道周期點的存在范圍和禁區(qū)邊界，進一步通過禁區(qū)邊界得到了混沌區(qū)域與周期n軌道區(qū)域的邊界的表達式，應(yīng)用Lyapunov指數(shù)對分析結(jié)果進行了驗證.

分岔；禁區(qū)邊界；混沌邊界；分段非線性不連續(xù)映射；Lyapunov指數(shù)

近些年，分段不連續(xù)映射動力系統(tǒng)在物理科學、工程、經(jīng)濟和醫(yī)學等方面已經(jīng)有了大量的研究[1-8].但以往的研究大多只考慮系統(tǒng)的邊界碰撞分岔，對系統(tǒng)的加周期現(xiàn)象和周期疊加現(xiàn)象進行理論分析和研究，對系統(tǒng)出現(xiàn)的混沌邊界的研究不多，尤其是對混沌區(qū)域與周期性區(qū)域邊界的研究較少.文獻[6]研究了一類具有可變禁區(qū)系統(tǒng)的不同類型的加周期分岔現(xiàn)象，發(fā)現(xiàn)系統(tǒng)的周期性區(qū)域與混沌區(qū)域的邊界線剛好處于左右邊界碰撞分岔同時發(fā)生的區(qū)域.文獻[9]建立了RC開關(guān)電路的一維分段映射模型，研究了系統(tǒng)參數(shù)變化誘導(dǎo)出的混沌吸引子的崩潰.現(xiàn)在已經(jīng)證明，魯棒混沌區(qū)接近由包含無數(shù)條分岔曲線的加周期結(jié)構(gòu)組成的周期性區(qū)域，這些分岔曲線形成加波段現(xiàn)象，這些現(xiàn)象標志著多波段混沌吸引子的發(fā)生[10].文獻[11]研究了TCP-RED離散反饋系統(tǒng)的周期穩(wěn)態(tài)經(jīng)過邊界碰撞分岔轉(zhuǎn)遷到混沌態(tài)的整個過程.文獻[12]研究了分段不連續(xù)映射模型，利用分岔解釋分段不連續(xù)映射誘導(dǎo)出的多波段混沌吸引子的結(jié)構(gòu).文獻[13]分析了分段不連續(xù)映射系統(tǒng)的整個混沌域的分岔結(jié)構(gòu)，發(fā)現(xiàn)系統(tǒng)的魯棒混沌吸引子嵌在兩個不同周期的周期性區(qū)域之間，給出了系統(tǒng)混沌吸引子的邊界表達形式.文獻[14]基于經(jīng)濟領(lǐng)域中的問題建立了一類分段不連續(xù)非線性映射模型，并分析了該系統(tǒng)的邊界碰撞分岔，給出了混沌區(qū)間的范圍.文獻[15]對如下不連續(xù)映射模型

進行了討論.研究結(jié)果表明，系統(tǒng)存在加周期現(xiàn)象和混沌現(xiàn)象.文獻[16]將模型（1）的線性分支變?yōu)閒（x）=ax+μ，主要討論了常數(shù)項對系統(tǒng)邊界碰撞分岔的影響，并對系統(tǒng)做了雙參數(shù)和三參數(shù)的分岔分析，但對混沌邊界沒有作相關(guān)研究.

本文主要研究文獻[16]中分段非線性映射系統(tǒng)分岔迭代結(jié)構(gòu)的禁區(qū)和系統(tǒng)出現(xiàn)的混沌的邊界.本文第1節(jié)通過理論分析系統(tǒng)的迭代表達式，得到系統(tǒng)的周期軌道中任意周期點的取值范圍，從而得到系統(tǒng)迭代點的禁區(qū)邊界曲線.第2節(jié)借助不連續(xù)映射的Lyapunov指數(shù)和分岔圖確定系統(tǒng)混沌區(qū)域與周期n軌道區(qū)域的邊界.在最后一節(jié)做出簡要的結(jié)論.

1 周期點的禁區(qū)

本文主要研究的模型為分段非線性映射，即

其中a＞1，0＜b＜1，0＜c＜1，μ＞-1.這里的x′代表映射的值，不代表x的導(dǎo)數(shù).當a+μ=c時，系統(tǒng)（2）為連續(xù)非線性映射，當a+μ≠c時，系統(tǒng)（2）為不連續(xù)非線性映射.記所研究的參數(shù)空間為

不連續(xù)點x=1將狀態(tài)空間分割為兩個區(qū)域

1.1 不動點分析

性質(zhì) 當(a,b,μ,c)∈P時，系統(tǒng)（2）的迭代點在(-∞,]內(nèi)隨著迭代次數(shù)單調(diào)遞減，在(,1]內(nèi)隨著迭代次數(shù)單調(diào)遞增.

證明 設(shè)x0為迭代的初始點，在S1區(qū)域內(nèi)迭代k次，則xk=f(xk-1),(k≥1)，經(jīng)遞推可得

當(a,b,μ,c)∈P時，有a＞1，由于為S1區(qū)域內(nèi)容許不動點，所以當x0≤時，xk≤xk-1；當1≥x0>時，xk>xk-1；即系統(tǒng)（2）的迭代點在(-∞,]內(nèi)單調(diào)遞減，在(,1]內(nèi)單調(diào)遞增.因此，S1區(qū)域內(nèi)的點經(jīng)過有限次迭代一定會到達S2區(qū)域內(nèi).

1.2 禁區(qū)分析

周期點的禁區(qū)表示周期軌道的周期點不存在的區(qū)間.若用Rk表示非線性函數(shù)g（x）迭代k次，用L表示線性函數(shù)（fx）迭代一次，形成的周期k+1解的符號序列記為RkL.

選取參數(shù)a=1.6,b=0.1,c=0.8,-0.4≤μ≤1得到全局分岔圖.圖1為添加了迭代禁區(qū)邊界線的全局分岔圖.觀察圖1容易得到，在μ>-0.2區(qū)域存在周期3軌道區(qū)域，周期4軌道區(qū)域，且在周期3軌道區(qū)域與周期4軌道區(qū)域之間夾有周期7軌道區(qū)域，也存在更復(fù)雜的周期n軌道區(qū)域.還觀察到在直線x=f(g(1))和直線x=g(f(1))之間沒有迭代點，直線x=g(f(g(1)))和直線x=g2(f(1))之間也為空白區(qū)域；同時還發(fā)現(xiàn)這些禁區(qū)邊界線將分岔結(jié)構(gòu)分割為若干“長帶”.若每一組禁區(qū)邊界線都相交，即分岔結(jié)構(gòu)的“長帶”之間沒有空白區(qū)域，迭代點鋪滿整個不變吸引區(qū)間，也就是說每一個μ所對應(yīng)的這組參數(shù)，系統(tǒng)（2）有無數(shù)個迭代點，這是不可能形成周期軌道的，此時系統(tǒng)（2）應(yīng)該出現(xiàn)混沌現(xiàn)象.由以上分析可得，系統(tǒng)（2）禁區(qū)邊界線的交點應(yīng)該落在周期n軌道區(qū)域與混沌區(qū)域的邊界線上.下一節(jié)，通過Lyapunov指數(shù)譜來驗證這個結(jié)論.

圖1 全局分岔圖Fig.1 Global bifurcation diagram

2 混沌邊界

Lyapunov指數(shù)是判別一個動力系統(tǒng)做周期運動還是混沌運動的一個重要物理量.根據(jù)系統(tǒng)（2），其Ly?apunov指數(shù)λ為

由（9）式可以得到系統(tǒng)（2）的Lyapunov指數(shù)圖.為了便于比較，我們將兩幅圖繪制在同一坐標下，見圖2（a）.

圖2（a）是禁區(qū)邊界曲線與Lyapunov指數(shù)圖合并圖，圖中最下方為Lyapunov指數(shù)曲線，圖2（b）為圖2（a）Lyapunov指數(shù)圖的放大圖.由圖2（a）知，在直線μ=μ0的左側(cè)，Lyapunov指數(shù)大于零，在右側(cè)，Lyapunov指數(shù)小于零，因此豎直線μ=μ0即為混沌區(qū)域與周期n軌道區(qū)域的分界線.圖2（a）的虛線由上到下分別是：

實線由上到下分別是：

邊界線x=f(g(1))與x=g(f(1))相交，x=g(f(g(1)))與x=g2(f(1))相交，等等.經(jīng)過觀察發(fā)現(xiàn)每一組禁區(qū)邊界線gi(f(1))和gi-1(f(g(1))),i=1,2,…,k都相交且交點都落在直線μ=μ0上；因此，禁區(qū)邊界線的交點可以確定直線μ=μ0.由式（8）可得：

圖2 (a)禁區(qū)邊界曲線與Lyapunov指數(shù)圖合并圖,(b)Lyapunov指數(shù)圖Fig.2(a)Combination of forbidden boundary curve and Lyapunov exponent diagram;(b)Lyapunov exponent diagram

（10）式即為混沌區(qū)域與周期n軌道區(qū)域的分界線解析式，注意這里計算得到的是兩個μ值，具體的選擇需要參照做分岔圖時μ所取的范圍.因為選取的參數(shù)為a=1.6,b=0.1,c=0.8,-0.4≤μ≤1，所以圖1中混沌區(qū)域與周期n軌道區(qū)域的邊界線為μ=-0.2092.

3 結(jié)論

本文針對具有非線性右側(cè)分支的分段一維映射模型（2）展開研究.在系統(tǒng)（2）的不變吸引區(qū)間內(nèi)，討論了周期軌道各周期點的存在范圍.得到了分岔圖中分岔結(jié)構(gòu)的迭代禁區(qū)的邊界線.對禁區(qū)邊界線進行理論分析發(fā)現(xiàn)，分段映射系統(tǒng)的禁區(qū)邊界線的交點落在混沌區(qū)域與周期n軌道區(qū)域的分界線上.進而給出了混沌區(qū)域邊界線的解析表達式.采用Lyapunov指數(shù)譜驗證了該結(jié)論的可靠性.

[1]沙金，包伯成，許建平，等.脈沖序列控制電流斷續(xù)模式Buck變換器的動力學建模與邊界碰撞分岔[J].物理學報，2012，61（12）：120501-1-120501-13.

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

[3]戴棟，馬西奎，李小峰.一類具有兩個邊界的分段光滑系統(tǒng)中邊界碰撞分岔現(xiàn)象及混沌[J].物理學報，2003，52（11）：2729-2736.

[4]Qu Z，Shiferaw Y，Weiss J N.Nonlinear dynamics of cardiac excitation-contraction coupling:an iterated map study[J].Physical Review E，2007，75（1）：112-118.

[5]Sushko I，Gardini L.Degenerate bifurcations and border collisions in piecewise smooth 1D and 2D maps[J].International Journal of Bifurcation&Chaos，2010，20（7）：2045-2070.

[6]楊科利.一類可變禁區(qū)的不連續(xù)系統(tǒng)的加周期分岔[J].物理學報，2015，64（12）：91-97.

[7]Gardini L，Tramontana F，Avrutin V.Border-collision bifurcations in 1D piecewise-linear maps and LEONOV's approach[J]. International Journal of Bifurcation&Chaos，2011，20（10）：3085-3104.

[8]秦志英，李群宏.一類非光滑映射的邊界碰撞分岔[J].力學學報，2013，45（1）：25-29.

[9]Avrutin V，Clüver M，Mahout V，et al.Bandcount adding structure and collapse of chaotic attractors in a piecewise linear bimodal map[J].Physica D:Nonlinear Phenomena，2015，309：37-56.

[10]Avrutin V，Schanz M，Gardini L.Self-similarity of the bandcount adding structures:Calculation by map replacement[J].Regular &Chaotic Dynamics，2010，15（15）：685-703.

[11]楊談，金躍輝，程時端.TCP-RED離散反饋系統(tǒng)中的邊界碰撞分岔及混沌控制[J].Acta Physica Sinica，2009，58（8）：5224-5237.

[12]Avrutin V，Schanz M.On the fully developed bandcount adding scenario[J].Nonlinearity，2008，21（5）：1077-1103.

[13]Avrutin V，Schanz M，Schenke B.Coexistence of the Bandcount-Adding and Bandcount-Increment Scenarios[J].Discrete Dy?namics in Nature&Society，2011，2011：701-716.

[14]Gardini L，Sushko I，Naimzada A K.Growing through chaotic intervals[J].Journal of Economic Theory，2008，143（1）：541-557.

[15]Jain P，Banerjee S.Border-collision bifurcations in one-dimensional discontinuous maps[J].International Journal of Bifurcation &Chaos,2011,13(11):3341-3351.

[16]韋麗梅.分段映射的邊界碰撞分岔及混沌控制[D].南寧：廣西大學，2014.

[17]Bernardo M di，Budd C J，Champneys A R，et al.Piecewise-smooth Dynamical Systems:Theory and Applications[M].London: Springer-Verlag，2008.

責任編輯：吳興華

Analysis of the Chaotic Boundary in a Class of Piecewise Nonlinear Mapping

XU Hongfei，LI Qunhong*，NING Min，SHANG Mengyuan
（School of Mathematics and Information Sciences，Guangxi University，Nanning530004，China）

In this paper the dynamical behaviors of a class of discontinuous one-dimensional mappings with a nonlinear branch are studied.This kind of models can be used in physical science,engineering,and medical science,and is also help?ful to the study of economics models.Takingμas a bifurcation parameter to draw the bifurcation diagram of the system,we find that in the invariant attracting region of the system there is an existence range for each point of the periodic orbit,and it leads to iteration forbidden region appearing in the bifurcation structure.By theoretical derivation,the article determines the existence ranges of the periodic orbits and the boundary of the forbidden region,obtains the boundary expression of the chaot?ic region and the period-norbits region by the boundary of the forbidden region,and finally verifies the analytic results by the Lyapunov exponents.

bifurcation；forbidden boundary；chaotic boundary；piecewise nonlinear mapping；Lyapunov exponent

O 317；O 193

：A

：1674-4942（2016）04-0363-06

10.12051/j.issn.1674-4942.2016.04.002

2016-09-17

廣西自然科學基金（2013GXNSFAA019017，2014GXNSFBA118024）

*通訊作者