李文娟, 牛瀟萌, 李旭超, 俞元洪

(1.赤峰學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,內(nèi)蒙古 赤峰 024000;2.赤峰學(xué)院計(jì)算機(jī)與信息工程學(xué)院,內(nèi)蒙古 赤峰 024000);

時(shí)滯擾動類Lorenz系統(tǒng)的Hopf分岔

李文娟1, 牛瀟萌1, 李旭超2, 俞元洪2

(1.赤峰學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,內(nèi)蒙古 赤峰 024000;2.赤峰學(xué)院計(jì)算機(jī)與信息工程學(xué)院,內(nèi)蒙古 赤峰 024000);

通過非線性動力學(xué)理論,對時(shí)滯類 Lorenz系統(tǒng)在平衡點(diǎn)的穩(wěn)定性問題和發(fā)生Hopf分岔的條件進(jìn)行了研究.首先計(jì)算得到系統(tǒng)的平衡點(diǎn),然后通過分析系統(tǒng)在平衡點(diǎn)處的相應(yīng)特征方程根的分布,得到系統(tǒng)在平衡點(diǎn)局部漸近穩(wěn)定和產(chǎn)生Hopf分岔的時(shí)滯臨界點(diǎn).以時(shí)滯為分叉參數(shù),研究了時(shí)滯系統(tǒng)存在Hopf分岔的條件.最后,利用Matlab程序進(jìn)行仿真驗(yàn)證所得結(jié)論與理論分析一致.本文的結(jié)論是對一些已有文獻(xiàn)研究成果的推廣.

時(shí)滯;類Lorenz系統(tǒng);穩(wěn)定性;Hopf分岔

1 引言

自1963年,美國氣象學(xué)家Lorenz在研究區(qū)域小氣候時(shí)提出了第一個經(jīng)典的Lorenz系統(tǒng)[1]以來,混沌系統(tǒng)得到了更廣泛的研究和應(yīng)用.此系統(tǒng)在混沌學(xué)歷史上有著重要的地位,特別是對它的分析在了解非線性方程如何出現(xiàn)混沌解方面很有意義.隨著對Lorenz系統(tǒng)研究的不斷深入,各種新的類Lorenz混沌系統(tǒng)大量產(chǎn)生,這種新的類Lorenz混沌系統(tǒng)在保密通信、激光物理、控制工程、化學(xué)反應(yīng)和生物醫(yī)學(xué)等領(lǐng)域應(yīng)用廣泛[2-4].現(xiàn)有的各種新的類Lorenz混沌系統(tǒng)仍無法滿足各領(lǐng)域應(yīng)用的需要,于是人們開始尋找新的類Lorenz混沌系統(tǒng),在此期間,很多學(xué)者得到一些行為復(fù)雜的系統(tǒng),例如Chen系統(tǒng)、Lu系統(tǒng)、Liu系統(tǒng)、Qi系統(tǒng)、Chua系統(tǒng)等[5-9].但是上述研究主要是混沌的同步和控制,對分岔的研究相對較少.對于時(shí)滯類Lorenz系統(tǒng)的研究更是不多.

2015年,文獻(xiàn)[10]提出了一類新的類Lorenz系統(tǒng)

其中x,y,z是狀態(tài)變量，a,b,c,d為系統(tǒng)參數(shù)．2015年,文獻(xiàn)[11]提出了一類帶時(shí)滯的類Lorenz系統(tǒng)

其中x,y,z是狀態(tài)變量;a,b,c,d為系統(tǒng)參數(shù);τ＞0為系統(tǒng)常量時(shí)滯．

本文的主要結(jié)果是將文獻(xiàn)[11]的系統(tǒng)(2)中第一個非線性方程中的a(y?x)改為

第二個非線性方程中的bx(t?τ)改為bx,第三個非線性方程中的dx2改為

得到一類新的雙時(shí)滯類Lorenz系統(tǒng).通過對一類新的雙時(shí)滯類Lorenz系統(tǒng)在零平衡點(diǎn)的線性化系統(tǒng)的特征方程的分析給出該系統(tǒng)在零平衡點(diǎn)的穩(wěn)定性問題和發(fā)生Hopf分岔的條件.最后通過數(shù)值擬真驗(yàn)證所得結(jié)論的正確性.本文研究的是雙時(shí)滯類Lorenz系統(tǒng),文獻(xiàn)[10]研究的是沒有時(shí)滯的類Lorenz系統(tǒng),文獻(xiàn)[11]研究的是單時(shí)滯的類Lorenz系統(tǒng).當(dāng)本文中的時(shí)滯為零時(shí)即為文獻(xiàn)[10]研究的系統(tǒng).

2 主要結(jié)果及證明

本文考慮時(shí)滯類Lorenz系統(tǒng):

其中x,y,z是狀態(tài)變量;a,b,c,d,h為系統(tǒng)參數(shù);τ＞0為系統(tǒng)常量時(shí)滯,可理解為信號傳輸?shù)淖璧K時(shí)間等.設(shè)系統(tǒng)參數(shù)

系統(tǒng)(3)的平衡點(diǎn)滿足下列條件:

由(4)式知系統(tǒng)(3)的平衡點(diǎn)有三個,它們分別為

2.1 平衡點(diǎn)穩(wěn)定性分析

定理 2.1假設(shè)τ=0,則系統(tǒng)(3)在平衡點(diǎn)

處是不穩(wěn)定的.

證明在平衡點(diǎn)

處,系統(tǒng)(3)的Jacobian矩陣為:

矩陣(5)對應(yīng)的特征方程為:

即

其中

當(dāng)τ=0時(shí),(6)式轉(zhuǎn)化為

因?yàn)橄到y(tǒng)參數(shù) a＞0,b＜0,c＞0,d＞0,h＞0,所以 2abch＜0根據(jù)羅斯 -霍維茲(Routh-Hurwitz)判據(jù)可知,(7)式的特征根具有正實(shí)部.所以當(dāng)τ=0時(shí),系統(tǒng)(3)在平衡點(diǎn)

處是不穩(wěn)定的.

注 2.1因?yàn)樵谧儞QS:(x,y,z)?→(?x,?y,z)的作用下具有不變性,即關(guān)于z具有對稱性.所以只討論系統(tǒng)(3)在平衡點(diǎn)

的穩(wěn)定性.

下面討論在平衡點(diǎn)O(0,0,0)處的穩(wěn)定性.在平衡點(diǎn)O(0,0,0)處易求得線性化系統(tǒng)

線性化系統(tǒng)(8)對應(yīng)的特征方程為

特征方程(9)可化為

其中

引理 2.1假設(shè)τ=0,則系統(tǒng)(3)在平衡點(diǎn)O(0,0,0)處是漸近穩(wěn)定的.

證明當(dāng)τ=0時(shí),(10)式轉(zhuǎn)化為

因?yàn)橄到y(tǒng)參數(shù)a＞0,b＜0,c＞0,d＞0,h＞0,所以(a+c)＞0,?abch＞0且有

根據(jù)羅斯 -霍維茲 (Routh-Hurwitz)判據(jù)可知,(11)式的所有特征根都具有負(fù)實(shí)部.所以當(dāng)τ=0時(shí),系統(tǒng)(3)在平衡點(diǎn)O(0,0,0)處是局部漸近穩(wěn)定的.

2.2 Hopf分岔問題

接下來討論在平衡點(diǎn)O(0,0,0)處的Hopf分岔存在的條件.

當(dāng)τ＞0時(shí),由于只考慮虛根,方程(10)等價(jià)于下面方程

設(shè)λ=iω(ω是大于零的常數(shù)）是(10)式的一個純虛根,則虛部ω滿足

根據(jù)復(fù)數(shù)相等可得

由三角函數(shù)的等式推得(13)式可化為

對(14)式有以下結(jié)論.

引理 2.2(14)式至少有一個正實(shí)根.

證明令u=ω2,則式(14)可化為

設(shè)

由(16)式得

根據(jù)函數(shù)零點(diǎn)存在定理,至少存在一個實(shí)數(shù)u0∈(0,+∞),使得f(uo)=0.所以(15)式至少有一個正實(shí)根．因?yàn)閡=ω2,從而(14)式至少有一個正實(shí)根.

設(shè)ω0為(14)式的正實(shí)根,則(10)式有一純虛根iω0.又由(13)式得

將ω=ω0代入方程(17),則時(shí)滯τ的值為

因此 (iω0,τk)是 (10)式的解,即當(dāng)時(shí)滯 τ=τk時(shí),λ=±iω0是 (10)式的一對共軛的純虛根.

設(shè)τ0=min{τk},則時(shí)滯τ=τ0是(10)式出現(xiàn)純虛根λ=±iω0時(shí)τ的最小值.故可得如下的引理3.

引理2.3如果

那么(10)式有一對純虛根

設(shè)(10)式的特征根

滿足

引理 2.4對于任意的τ=τk/=0,則

證明對(10)式兩邊關(guān)于τ求導(dǎo)，可得

由(10)式可得

易得

由(21)式和

可得

由于

同號.證畢.

根據(jù)引理2.1-引理2.4和Hopf分岔理論可得下面結(jié)論．

定理2.2如果

那么

(1)當(dāng)τ∈[0,τ0)時(shí),系統(tǒng)(3)在平衡點(diǎn)O(0,0,0)是漸近穩(wěn)定的;

(2)當(dāng)τ＞τ0時(shí),系統(tǒng)(3)在平衡點(diǎn)O(0,0,0)是不穩(wěn)定的;

(3)當(dāng)

時(shí),系統(tǒng)(3)在平衡點(diǎn)O(0,0,0)處發(fā)生Hopf分岔,產(chǎn)生極限環(huán).

注 2.2定理的結(jié)論可以由下部分的數(shù)值仿真來驗(yàn)證.文獻(xiàn)[1]的系統(tǒng)是系統(tǒng)(3)中

時(shí)的特例.此外,定理也推廣了文獻(xiàn)[2]的結(jié)果.

3 數(shù)值仿真

時(shí)滯類Lorenz系統(tǒng)(3)的參數(shù)a＞0,b＜0,c＞0,d＞0,令

這時(shí)系統(tǒng)(3)可化為

利用Matlab軟件計(jì)算得(11)式的正實(shí)根

和(15)式中

故,由定理可得下面的推論.

推論3.1如果

則

(1)當(dāng)τ∈[0,0.7111)時(shí),系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)是漸近穩(wěn)定的;

(2)當(dāng)τ＞0.7111時(shí),系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)是不穩(wěn)定的;

(3)當(dāng)

時(shí),系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)處發(fā)生Hopf分岔,產(chǎn)生極限環(huán).

下面用Matlab軟件繪出時(shí)滯τ取不同值時(shí),系統(tǒng)(23)的狀態(tài)變量隨時(shí)間t的軌線圖和相圖,驗(yàn)證所得結(jié)論的正確性.

圖1-1 時(shí)滯τ=0.7時(shí),系統(tǒng)(23)的狀態(tài) 變量x隨時(shí)間t的軌線圖

圖1-2 時(shí)滯τ=0.7時(shí),系統(tǒng)(23)的狀態(tài) 變量y隨時(shí)間t的軌線圖

圖1-3 時(shí)滯τ=0.7時(shí),系統(tǒng)(23)的狀態(tài) 變量z隨時(shí)間t的軌線圖

圖1-4 滯τ=0.7時(shí),系統(tǒng)(23)的狀態(tài) 在xoy平面內(nèi)的相圖

從圖 1-1、圖1-2、圖1-3可以看出，當(dāng)時(shí)滯τ=0.7時(shí),系統(tǒng)(23)的狀態(tài)變量x,y,z的值隨時(shí)間t的增大而趨于平衡點(diǎn)O(0,0,0),所以系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)是漸近穩(wěn)定的．

圖2-1 時(shí)滯τ=0.7111時(shí),系統(tǒng)(23)的狀態(tài)變量x隨時(shí)間t的軌線圖

圖2-2 時(shí)滯τ=0.7111時(shí),系統(tǒng)(23)的狀態(tài) 變量y隨時(shí)間t的軌線圖

圖2-3 時(shí)滯τ=0.7111時(shí),系統(tǒng)(23)的狀態(tài)}變量z隨時(shí)間t的軌線圖

圖2-4 時(shí)滯τ=0.7111時(shí),系統(tǒng)(23)的狀態(tài) 在xoy平面內(nèi)的相圖

從圖2-1、圖 2-2、圖 2-3可以看出,當(dāng)時(shí)滯 τ=0.7111時(shí),系統(tǒng)(23)的狀態(tài)變量 x,y,z的值隨時(shí)間t的增大保持等周期震蕩,所以系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)處發(fā)生Hopf分岔,產(chǎn)生極限環(huán).

圖3-1 時(shí)滯 τ=0.72時(shí),系統(tǒng)(23)的狀態(tài)變量x隨時(shí)間t的軌線圖

圖3-2 時(shí)滯τ=0.72時(shí),系統(tǒng)(23)的狀態(tài) 變量y隨時(shí)間t的軌線圖

圖3-3 時(shí)滯 τ=0.72時(shí),系統(tǒng)(23)的狀態(tài) 變量z隨時(shí)間t的軌線圖

圖3-4 時(shí)滯τ=0.72時(shí),系統(tǒng)(23)在xoy平面里的相圖

從圖3-1、圖3-2、圖3-3可以看出,當(dāng)時(shí)滯τ=0.72時(shí),系統(tǒng)(23)的狀態(tài)變量x,y,z的值隨時(shí)間t的增大而逐漸遠(yuǎn)離平衡點(diǎn)O(0,0,0),說明系統(tǒng)(23)在平衡點(diǎn)O(0,0,0)是不穩(wěn)定的.

圖 1-4、圖 2-4、圖 3-4可以看出,當(dāng)時(shí)滯 τ＜τ0,τ=τ0,τ＞τ0時(shí),系統(tǒng) (23)xoy平面內(nèi)的相圖.

[1]Lorenz E N.Deterministic non-periodic fl ows[J].J.Atmos.Sci.,1963,20:130-141.

[2]劉式適,劉式達(dá).物理學(xué)中的非線性方程[M].北京:北京大學(xué)出版社,2012:341-345.

[3]王興元.混沌系統(tǒng)的同步及在保密通訊中的應(yīng)用[M].北京:科學(xué)出版社,2012:356-385.

[4]劉揚(yáng)正.超混沌Lü系統(tǒng)的電路實(shí)現(xiàn)[J].物理學(xué)報(bào),2008,57(3):1439-1443.

[5]Wenxin Q,Gangrong C.On the boundedness of solutions of the Chen system[J].J.Math.Anal.Appl.,2007,329(1):445-451.

[6]Mkaouar H,Boubaker O.Chaos synchronization for master slave piecewise linear systems application to Chua’s circuit[J].Commun Nonlinear Sci.Numer Simulat.,2012,17(3):1292-1302.

[7]Chongxin L,Tao L,Ling L,et al.A new chaotic attractor[J].Chaos,Solitons and Fractals,2004,22:1031-1038.

[8]Yongguang Y,Suochun Z.Hopf bifurcation in the Lü system[J].Chaos,Solitons and Fractals,2003,17(5):901-906.

[9]Guoyuan Q,Gangrong C,Shengzhi D,et al.Analysis of a new chaotic system[J].Physica A,2005,352:295-308.

[10]官國榮,吳成茂,賈倩.一種改進(jìn)的高性能Lorenz系統(tǒng)構(gòu)造及其應(yīng)用[J].物理學(xué)報(bào),2015,64(2):1-14.

[11]李德奎,連玉平.單時(shí)滯類Lorenz系統(tǒng)的Hopf分岔分析[J].數(shù)學(xué)雜志,2015,35(3):633-642.

Hopf bifurcation analysis of the disturbed Lorenz-like System with the delayed

Li Wenjuan1,Niu Xiaomeng1,Li Xvchao2,Yu Yuanhong3
(1.School of Mathematics and Statistics,Chifeng University,Chifeng 024000,China;2.Computer and Information Technology Institute,Chifeng University,Chifeng 024000,China;3.Academy of Mathematics System Sciences,Chinese Academy of Sciences,Beijing 100190,China)

In this paper,we analyze the stability and Hopf bifurcation condition of the system at its balance by using the nonlinear dynamics theory.The equilibrium point is obtained and then critical point of the delay for local stability of the equilibrium and existence of local Hopf bifurcation is also obtained by analyzing distribution of roots of the corresponding characteristic equation.the condition for the existence of Hopf bifurcation of this delay system is studied by taking delay as bifurcation parameter.Furthermore,by Matlab program some numerical simulations were given to show the correctness of the obtained conclusion.This paper extends the research achievements in the cited literature.

time-delay,Lorenz-like system,stability,Hopf bifurcation

O175

A

1008-5513(2017)05-0475-11

10.3969/j.issn.1008-5513.2017.05.005

2017-03-27.

國家自然科學(xué)基金(11561001);內(nèi)蒙古自然科學(xué)基金(2014MS0101);內(nèi)蒙古高等學(xué)?？蒲谢穑∟JZY17301).

李文娟(1981-),碩士,講師,研究方向：穩(wěn)定性理論及應(yīng)用.

2010 MSC:34C10,34C15