蔡 萍, 唐駕時(shí), 李震波
(1. 閩南師范大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 福建 漳州 363000; 2.湖南大學(xué) 機(jī)械與運(yùn)載工程學(xué)院, 湖南 長(zhǎng)沙 410082)
分岔和混沌是非線性科學(xué)的一個(gè)重要研究領(lǐng)域,也是非線性微分方程研究的重要組成部分。雖然, 關(guān)于分岔理論的研究近年來(lái)十分活躍, 取得了眾多成果[1~3],但在控制方面,分岔遠(yuǎn)不及混沌成熟。分岔控制是指設(shè)計(jì)一個(gè)控制器,它能夠修改一個(gè)非線性系統(tǒng)的分岔特性,從而實(shí)現(xiàn)人們所要的理想動(dòng)力學(xué)行為。Hopf分岔是一種動(dòng)態(tài)分岔,具有重要的理論意義和實(shí)際價(jià)值。關(guān)于Hopf分岔的控制已引起了許多科學(xué)工作者的極大關(guān)注。劉爽等[4]設(shè)計(jì)非線性反饋控制器對(duì)一類(lèi)耦合非線性相對(duì)轉(zhuǎn)動(dòng)系統(tǒng)的Hopf分岔進(jìn)行了控制。郁培等[5]利用非線性狀態(tài)反饋,給出了推遲、消除Hopf 分岔以及改變分岔解穩(wěn)定性的控制方法。劉曾榮等[6]結(jié)合反饋與非反饋控制,對(duì)一類(lèi)三維系統(tǒng)的Hopf分岔進(jìn)行了控制。近年來(lái),Hopf分岔控制理論也成功地應(yīng)用于神經(jīng)網(wǎng)絡(luò)模型中[7,8],成為神經(jīng)網(wǎng)絡(luò)非線性動(dòng)力學(xué)的一個(gè)重要組成部分。然而,在沒(méi)有辦法或沒(méi)有必要消除Hopf分岔的情形下,理想的目標(biāo)是對(duì)極限環(huán)的振幅與頻率進(jìn)行控制,這方面也有很多的工作[9~11]。
van der Pol 振子是一個(gè)著名的自激振動(dòng)系統(tǒng),最早出現(xiàn)在非線性電路中,后來(lái)在機(jī)械、生物、化學(xué)等很多工程和科學(xué)領(lǐng)域中發(fā)現(xiàn)了它的存在。以往對(duì)van der Pol 振子的分岔控制研究,主要集中在van der Pol振子、廣義的van der Pol振子和耦合的van der Pol 振子這三類(lèi)系統(tǒng)。對(duì)參數(shù)激勵(lì)作用下的廣義van der Pol 振子的研究并不多見(jiàn)。本文考慮如下含有參數(shù)激勵(lì)的廣義van der Pol方程
(1)
含參數(shù)激勵(lì)的非線性系統(tǒng)有著豐富的動(dòng)力學(xué)性態(tài),1/2亞諧共振是其中之一。唐駕時(shí)等[12]做了該Mathieu-van der Pol系統(tǒng)在強(qiáng)非線性下的分岔分析,但相關(guān)的分岔控制并未有文獻(xiàn)給出。本文研究的目的是設(shè)計(jì)線性位移和速度時(shí)滯反饋控制器,對(duì)系統(tǒng)(1)的1/2亞諧共振下的分岔響應(yīng)進(jìn)行控制。受控的系統(tǒng)可以描述為
(2)
其中,τ1和τ2為時(shí)滯量,g1和g2為增益系數(shù)。首先,通過(guò)多尺度法得到受控系統(tǒng)的分岔響應(yīng)方程,給出了極限環(huán)幅值與控制參數(shù)之間的關(guān)系;其次,通過(guò)零解穩(wěn)定性分析,得到Hopf分岔產(chǎn)生的條件以及控制參數(shù)的取值范圍;最后,通過(guò)數(shù)值模擬說(shuō)明理論結(jié)果的正確性和控制的有效性。
采用多尺度法[13],設(shè)方程(2)的攝動(dòng)解形式為
x=x0(T0,T1)+εx1(T0,T1)+…
(3)
其中Ti=εit(i=0,1,…) .討論1/2 亞諧共振情況,即Ω=2ω0+εσ,則
(4)
其中,σ為調(diào)諧參數(shù)。將(3),(4)式代入(2)式,比較方程兩邊ε的同次冪系數(shù),得到下列攝動(dòng)方程組
(5)
cosΩT0x0+g1x0(t-τ1)+g2D0x0(t-τ2)
(6)
(7)
(8)
令方程(8)中的久期項(xiàng)為零得
(9)
令A(yù)=aeiφ,其中a和φ是關(guān)于T1的實(shí)函數(shù),代入(9)式,得系統(tǒng)極坐標(biāo)形式的平均方程
(10)
(11)
這里
(12)
在方程(9)中令A(yù)=u+iv,u和v為關(guān)于T1的實(shí)函數(shù),則可得直角坐標(biāo)形式的平均方程
(13)
系統(tǒng)(13)的平衡點(diǎn)在原點(diǎn)(u,v)=(0,0) ,平均方程(13)的Jacobian矩陣為
(14)
這里w=(u,v)T.對(duì)應(yīng)零解的特征方程為
(15)
由Routh-Huriwigz (R-H)判據(jù)可知,當(dāng)滿(mǎn)足條件(16)時(shí)
(16)
原點(diǎn)是漸近穩(wěn)定的平衡點(diǎn)。并且當(dāng)M=0時(shí),方程(15)存在一對(duì)純虛根,根據(jù)Hopf 分岔定理,該情況系統(tǒng)發(fā)生Hopf分岔。由于M 是時(shí)滯項(xiàng)系數(shù)g1,g2和時(shí)滯參數(shù)τ1,τ2的函數(shù),因此,調(diào)節(jié)g1,g2和τ1,τ2可使系統(tǒng)的平衡點(diǎn)由穩(wěn)定變?yōu)椴环€(wěn)定并從中生長(zhǎng)出極限環(huán),即可以控制系統(tǒng)Hopf 分岔的發(fā)生。由式(11)可知,參數(shù)g1,g2和τ1,τ2同樣影響系統(tǒng)極限環(huán)幅值大小,即參數(shù)可以控制極限環(huán)的大小。由上述分析可知, 系統(tǒng)的線性部分能夠影響Hopf分岔的產(chǎn)生, 因此非線性控制器僅僅能夠改變極限環(huán)的幅值, 而無(wú)法控制Hopf分岔. 線性控制器雖然既能改變極限環(huán)幅值也能控制Hopf分岔的發(fā)生, 但顯然時(shí)滯反饋控制器的調(diào)節(jié)幅度和空間更大,在實(shí)際工程中的應(yīng)用也會(huì)更加廣泛。這是時(shí)滯反饋控制器的優(yōu)越性所在,也是其它控制器無(wú)法做到的。
(a)x-t 圖圖
(a)x-t 圖圖
(a)x-t 圖圖
建立了含參數(shù)激勵(lì)的廣義Van der Pol方程的時(shí)滯動(dòng)力系統(tǒng),主要討論了在1/2亞諧共振下控制參數(shù)對(duì)系統(tǒng)分岔的控制和極限環(huán)幅值的影響。通過(guò)攝動(dòng)分析,得到了時(shí)滯動(dòng)力系統(tǒng)的分岔響應(yīng)方程,給出極限環(huán)幅值和控制參數(shù)的關(guān)系。借助R-H判據(jù)對(duì)系統(tǒng)平衡點(diǎn)的穩(wěn)定性做了分析,從而獲得系統(tǒng)產(chǎn)生Hopf分岔的條件。數(shù)值模擬結(jié)果表明,時(shí)滯控制器不僅能控制極限環(huán)的幅值,也能控制系統(tǒng)Hopf分岔的產(chǎn)生,與理論分析一致,得到滿(mǎn)意的結(jié)果。因此,時(shí)滯系統(tǒng)在分岔控制領(lǐng)域有著廣泛的應(yīng)用前景。
參考文獻(xiàn):
[1]陳予恕, 唐 云. 非線性動(dòng)力學(xué)的現(xiàn)代分析方法[M]. 北京: 科學(xué)出版社, 1992.
[2]劉延柱, 陳立群. 非線性振動(dòng)[M]. 北京: 高等教育出版社, 2001.
[3]聞邦椿, 李以農(nóng), 韓清凱. 非線性振動(dòng)理論中的解析方法及工程應(yīng)用[M]. 沈陽(yáng): 東北大學(xué)出版社, 2001.
[4]劉 爽, 劉浩然. 一類(lèi)耦合非線性相對(duì)轉(zhuǎn)動(dòng)系統(tǒng)的Hopf分岔控制[J]. 物理學(xué)報(bào), 2010, 8(59): 5223~5228.
[5]Yu P, Chen G R. Hopf bifurcation control using nonlinear feedback with polynomial functions [J]. International Journal of Bifurcation and Chaos, 2004, 14 (5): 1683~1704.
[6]Liu Z R, Chung K W. Hybrid control of bifurcation in continuous nonlinear dynamical systems coupled chaotic systems [J]. International Journal of Bifurcation and Chaos, 2005, 15 (12): 3895~3903.
[7]張 洪, 陳天麒. 時(shí)延神經(jīng)元的Hopf分岔控制[J]. 控制理論與應(yīng)用, 2006, 4(23): 181~186 .
[8]趙洪涌, 陳 凌. 一類(lèi)慣性神經(jīng)網(wǎng)絡(luò)的分岔與控制[J]. 物理學(xué)報(bào), 2011, 60(7): 070202.
[9]Tang J S, Chen Z L. Amplitude control of limit cycle in van der Pol system [J]. International Journal of Bifurcation and Chaos, 2006, 16 (2): 487~495.
[10]Lee Y S.Suppression of limit cycle oscillations in the van der Pol oscillator by means of passive nonlinear energy sinks [J]. Structural of Control and Health Monitoring, 2006, 13: 41~75.
[11]沈建和, 陳樹(shù)輝. 非線性振動(dòng)系統(tǒng)極限環(huán)振幅與頻率的控制 [J]. 振動(dòng)與沖擊, 2009, 28 (6): 90~92 .
[12]Tang J S, Qin J Q. Bifurcations of a generalized van der Pol oscillator with strong nonlinearity [J]. Journal of Sound and Vibration, 2007, 306: 890~896.
[13]Nayfeh A H, Mook D T.Nonlinear Oscillations[M]. New York: Wiley, 1978.
湖北師范大學(xué)學(xué)報(bào)(自然科學(xué)版)2013年4期