曾尖尖 鄭遠(yuǎn)廣

(南昌航空大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院,南昌 330063)

引言

混沌是非線性系統(tǒng)中普遍存在的一種動(dòng)力學(xué)現(xiàn)象,其典型的特征是對(duì)擾動(dòng)的敏感性和不可預(yù)測(cè)性.在許多實(shí)際工程系統(tǒng)中,混沌常常是有害的.生態(tài)系統(tǒng)中的混沌會(huì)導(dǎo)致物種的突然消亡[1-2];被動(dòng)行走機(jī)器人中的混沌會(huì)降低行走效率,同時(shí)增加能量消耗[3-4];網(wǎng)絡(luò)信息傳輸中的混沌會(huì)導(dǎo)致網(wǎng)絡(luò)堵塞[5].因此,混沌控制是一項(xiàng)重要的研究課題[6-7].在早期的研究中,人們常常通過(guò)改變系統(tǒng)的分叉參數(shù)[8]來(lái)達(dá)到消除混沌的目的,或者用外部周期信號(hào)去同步混沌系統(tǒng),從而引導(dǎo)混沌系統(tǒng)到期望的周期運(yùn)動(dòng)上來(lái).采用這些控制方法控制混沌時(shí),受控系統(tǒng)的周期解軌線并不是原系統(tǒng)的混沌吸引子中包含的軌線,因此,通常需要比較大的控制力和比較大的能量消耗.

由于混沌吸引子中包含無(wú)數(shù)個(gè)不穩(wěn)定周期軌線[9],基于這一事實(shí),人們提出了OGY 控制方法[10].當(dāng)系統(tǒng)運(yùn)動(dòng)到預(yù)定的不穩(wěn)定周期軌道附近時(shí),對(duì)系統(tǒng)關(guān)鍵參數(shù)進(jìn)行微小的攝動(dòng)從而把系統(tǒng)穩(wěn)定到該周期軌道上.由于該方法是把混沌系統(tǒng)引導(dǎo)到混沌吸引子中原有的不穩(wěn)定周期軌道上,因此僅需要微小的控制力和少量的能量消耗.然而,OGY 控制方法對(duì)外界擾動(dòng)極度敏感,當(dāng)系統(tǒng)受到擾動(dòng)遠(yuǎn)離預(yù)定的周期軌道運(yùn)動(dòng)時(shí),該方法無(wú)法實(shí)施.為了克服這一缺點(diǎn),人們提出了時(shí)滯反饋控制(DFC)[11].DFC 用系統(tǒng)輸出狀態(tài)信號(hào)和輸出狀態(tài)信號(hào)的適當(dāng)滯后值之差生成反饋控制信號(hào),其中的時(shí)滯量設(shè)置成等于預(yù)定要控制的不穩(wěn)定周期軌道的周期.當(dāng)受控周期軌道穩(wěn)定后,DFC 的控制信號(hào)消失,因此DFC 不會(huì)改變?cè)到y(tǒng)的周期軌道,僅改變其穩(wěn)定性.在使用DFC 對(duì)混沌系統(tǒng)進(jìn)行控制時(shí),不需要知道系統(tǒng)精確的數(shù)學(xué)模型,只需要測(cè)得目標(biāo)周期軌道的周期.DFC 使用方便,且具有較好的控制效果,在實(shí)際工程中得到廣泛的應(yīng)用[12-14].

在經(jīng)典DFC 的使用中,人們發(fā)現(xiàn)該方法存在一定的局限性.比如,DFC 存在奇數(shù)限制(ONL)[15].ONL 表明DFC 不能穩(wěn)定具有奇數(shù)個(gè)實(shí)部大于零的條件Lyapunov 指數(shù)的周期軌線.此外,用DFC 穩(wěn)定某些不穩(wěn)定周期軌線時(shí),受控周期軌線常常具有比較小的穩(wěn)定區(qū)域[16].這里的穩(wěn)定區(qū)域是指參數(shù)空間的某個(gè)區(qū)域,當(dāng)控制參數(shù)在這個(gè)區(qū)域中取值時(shí),受控周期軌線是穩(wěn)定的.為了提高經(jīng)典DFC 的控制效果,人們對(duì)該方法進(jìn)行了多方面的改進(jìn)[16-18].比如,人們提出了半周期延遲反饋控制來(lái)避免DFC 的奇數(shù)限制[19];文獻(xiàn)[20]基于梯度下降法設(shè)計(jì)了一種自動(dòng)調(diào)整反饋增益的方法,該方法可以有效地確定穩(wěn)定目標(biāo)周期軌線的反饋增益系數(shù);為了穩(wěn)定高周期的不穩(wěn)定周期軌線,文獻(xiàn)[16]提出了廣義時(shí)滯反饋控制方法(EDFC).

切換控制策略在控制工程中得到廣泛應(yīng)用[21-23].Bang-bang 控制的控制信號(hào)在兩個(gè)常值之間切換[24],該方法能夠?qū)崿F(xiàn)最短時(shí)間控制.滑?？刂剖橇硪环N常用的切換控制技術(shù).滑?？刂频闹饕枷胧?通過(guò)不連續(xù)地切換控制信號(hào),使受控系統(tǒng)到達(dá)并保持在設(shè)定的滑模面上,從而把系統(tǒng)穩(wěn)定到期望的平衡狀態(tài)[25].在系統(tǒng)耦合同步中,切換策略也被用來(lái)改變系統(tǒng)的耦合方式以提高系統(tǒng)的同步能力[26].周期耦合是一種依賴于時(shí)間和切換頻率的耦合方式[27-28],只有在特征的時(shí)間段上,系統(tǒng)之間存在耦合,而在其他時(shí)間段上,系統(tǒng)之間不存在耦合.研究表明,具有適當(dāng)切換頻率的周期耦合系統(tǒng)同步性能要顯著優(yōu)于標(biāo)準(zhǔn)耦合系統(tǒng)的同步性能[26].為了理解周期耦合對(duì)同步流形穩(wěn)定性的影響,文獻(xiàn)[29]通過(guò)計(jì)算同步流形的局部條件Lyapunov 指數(shù),發(fā)現(xiàn)周期耦合以一種非線性的方式修正了同步流形的穩(wěn)定性.

在混沌控制的工程實(shí)踐中,受控系統(tǒng)常常存在某些不確定性因素,同時(shí)控制參數(shù)也會(huì)受到一定的限制.要實(shí)現(xiàn)對(duì)目標(biāo)周期軌線的穩(wěn)定控制,就需要受控周期軌線具有足夠大的穩(wěn)定區(qū)域.因此,采用適當(dāng)?shù)姆椒ㄟM(jìn)一步擴(kuò)大受控周期軌線的穩(wěn)定區(qū)域具有重要的工程實(shí)踐意義.本文基于時(shí)變切換策略[26,30]對(duì)經(jīng)典DFC 控制方法進(jìn)行改進(jìn),提出了時(shí)變切換時(shí)滯反饋控制方法(TSDFC).其目的是提高經(jīng)典DFC 的控制效果,擴(kuò)大受控周期軌線的穩(wěn)定區(qū)域.本文的結(jié)構(gòu)如下:在第二節(jié)中,提出時(shí)變切換時(shí)滯反饋控制方法.在第三節(jié)中,通過(guò)實(shí)例分析,研究TSDFC 的控制性能.最后在第四節(jié)中,對(duì)本文的研究工作進(jìn)行總結(jié)并對(duì)相關(guān)問(wèn)題進(jìn)行討論.

1 TSDFC

考慮受控的非線性系統(tǒng)

式中,G為反饋增益系數(shù)矩陣,τ為時(shí)滯.

選取系統(tǒng)中某個(gè)周期為T的不穩(wěn)定周期軌線X*(t) 作為待鎮(zhèn)定的目標(biāo)周期軌線.DFC 是將時(shí)滯τ設(shè)定等于目標(biāo)周期軌線X*(t) 的周期T,再通過(guò)調(diào)節(jié)反饋增益系數(shù)矩陣G,使受控不穩(wěn)定周期軌線X*(t)達(dá)到穩(wěn)定.當(dāng)受控周期軌線X*(t) 穩(wěn)定時(shí),反饋控制信號(hào)U(t) 趨于無(wú)窮小.

雖然DFC 使用方便,但受控周期軌線X*(t) 常常僅有較小的穩(wěn)定區(qū)域.為了提高DFC 的控制效果,擴(kuò)大受控周期軌線X*(t) 的穩(wěn)定區(qū)域,本文基于時(shí)變切換策略[26,30]對(duì)DFC 進(jìn)行改進(jìn),提出了TSDFC.TSDFC 的反饋控制信號(hào)為

式中,反饋增益系數(shù)矩陣為

式中,s gn(x)=G為定常反饋增益系數(shù)矩陣,ω為切換頻率,其決定了切換速度的快慢.TSDFC的控制信號(hào)僅在特征的時(shí)段上存在,而在其他時(shí)段控制信號(hào)不存在,即控制信號(hào)按照給定的頻率 ω 進(jìn)行周期性切換.

為確定受控周期軌線X*(t) 的穩(wěn)定區(qū)域,求得主穩(wěn)定方程為[28]

式中,η(t)=X(t)-X*(t)為偏差變量,DF(t,X*(t))為F(t,X(t))在X*(t) 上的雅可比矩陣.X*(t) 的穩(wěn)定性由主穩(wěn)定方程(5)的條件Lyapunov 指數(shù)確定.當(dāng)主穩(wěn)定方程(5)的所有條件Lyapunov 指數(shù)均具有負(fù)實(shí)部時(shí),受控周期軌線X*(t) 穩(wěn)定.

主穩(wěn)定方程(5)是時(shí)滯系統(tǒng),其初值狀態(tài)空間是無(wú)窮維Banach 空間C=C([-T,0],Rm).在Banach 空間C中的初始函數(shù)定義為ηt0(θ)=η(t0+θ),其中,t0為初始時(shí)間,θ ∈[-T,0].系統(tǒng)(5)的最大條件Lyapunov指數(shù)定義為其中

由于系統(tǒng)(5)的初值狀態(tài)空間是無(wú)窮維的,從而無(wú)法通過(guò)對(duì)系統(tǒng)(5) 進(jìn)行積分來(lái)計(jì)算最大條件Lyapunov 指數(shù) λm[31].為克服這一困難,本文采用Runge-Kutta 方法,將系統(tǒng)(5)轉(zhuǎn)化成具有有限維狀態(tài)空間的離散系統(tǒng),進(jìn)而求得最大條件Lyapunov 指數(shù) λm[32],詳細(xì)的計(jì)算過(guò)程請(qǐng)參考附錄A.

將TSDFC 應(yīng)用于系統(tǒng)(1)得到

當(dāng)ω=0時(shí),TSDFC 退化為DFC.

2 TSDFC 鎮(zhèn)定不穩(wěn)定周期軌線

本節(jié)用TSDFC 控制兩個(gè)典型的非線性混沌系統(tǒng),考察該方法鎮(zhèn)定不穩(wěn)定周期軌線的具體效果.首先,通過(guò)求得受控周期軌線X*(t) 的穩(wěn)定區(qū)域與切換頻率 ω 的關(guān)系圖,找到最優(yōu)切換頻率,使受控周期軌線X*(t) 的穩(wěn)定區(qū)域達(dá)到最大.其次,對(duì)比TSDFC和經(jīng)典DFC 的控制效果,闡述TSDFC 在鎮(zhèn)定混沌系統(tǒng)的不穩(wěn)定周期軌線的有效性.

2.1 TSDFC 控制 R?ssler 系統(tǒng)

受TSDFC 控制的 R?ssler 系統(tǒng)方程為

式中,x,y,z為系統(tǒng)狀態(tài)變量,a,b,c為系統(tǒng)參數(shù),為反饋增益系數(shù),其在定常值g和0之間周期切換.不失一般性,這里假設(shè)狀態(tài)變量x為唯一可獲得的系統(tǒng)狀態(tài)變量,則反饋控制信號(hào)表示為Ut(t)=gt(x(t-T)-x(t)).

將系統(tǒng)參數(shù)設(shè)定為a=0.2,b=0.2,c=5.7,此時(shí)R?ssler系統(tǒng)呈現(xiàn)混沌動(dòng)力學(xué)行為,其混沌吸引子如圖1 所示.

圖1 R?ssler 系統(tǒng)的混沌吸引子Fig.1 The chaotic attractor of R?ssler system

選擇 R?ssler 混沌吸引子中兩個(gè)典型的不穩(wěn)定周期軌線Xk(t)=(xk(t),yk(t),zk(t))T(k=1,2) 作為待鎮(zhèn)定的目標(biāo)周期軌線.圖2 展示了目標(biāo)周期軌線Xk(t)在xoy平面上的投影.

圖2 R?ssler 系統(tǒng)的不穩(wěn)定周期軌線在 x oy 平面上的投影Fig.2 The projection of the unstable periodic trajectory ofR?ssler system on the x oy plane

為確定受控周期軌線Xk(t) 的穩(wěn)定區(qū)域,求得主穩(wěn)定方程為

式中,時(shí)滯 τk等于目標(biāo)周期軌線Xk(t) 的周期(T1=5.880,T2=11.760).當(dāng)受控周期軌線Xk(t) 的最大條件Lyapunov 指數(shù) λm的實(shí)部小于零,受控周期軌線Xk(t) 穩(wěn)定.受控周期軌線Xk(t) 的穩(wěn)定區(qū)域是反饋增益系數(shù)g的某個(gè)區(qū)間,在這個(gè)區(qū)間上受控周期軌線穩(wěn)定,該區(qū)間的長(zhǎng)度記為穩(wěn)定區(qū)域?qū)挾萕s.受控周期軌線Xk(t) 的穩(wěn)定區(qū)域越大表明控制效果越好.

當(dāng)ω=0時(shí),TSDFC 退化成經(jīng)典DFC.通過(guò)計(jì)算受控周期軌線Xk(t) 的最大條件Lyapunov 指數(shù),得到受控周期軌線Xk(t) 的穩(wěn)定區(qū)域.圖3 中紅色虛線是應(yīng)用DFC 鎮(zhèn)定不穩(wěn)定周期軌線Xk(t) 時(shí)的最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g變化的曲線圖.同時(shí)在表1 中給出了受控周期軌線Xk(t) 的穩(wěn)定區(qū)域和穩(wěn)定區(qū)域?qū)挾萕s.結(jié)果表明受DFC 控制的周期軌線Xk(t) 的穩(wěn)定區(qū)域較小,這表明DFC 的控制效果比較差.

為提 高經(jīng)典DFC 的控制性能,令 ω >0,這時(shí)TSDFC 的控制信號(hào)將依據(jù)時(shí)段進(jìn)行切換,其切換過(guò)程如圖4 所示,在灰色時(shí)段上,控制信號(hào)存在,而在灰色時(shí)段之外,控制信號(hào)消失.TSDFC 的控制效果與切換頻率 ω 密切相關(guān).為確定最佳切換頻率 ω,對(duì)應(yīng)于周期軌線X1(t)和X2(t),分別求得各自的穩(wěn)定區(qū)域?qū)挾萕s隨切換頻率ω 的變化曲線(見圖5).從圖5 中可見,切換過(guò)程以非線性方式修改了受控不穩(wěn)定周期軌線的穩(wěn)定性和穩(wěn)定區(qū)域[33],穩(wěn)定區(qū)域?qū)挾萕s隨切換頻率 ω 的增大而非平滑地變化.對(duì)應(yīng)于周期軌線X1(t)和X2(t) 的最佳切換頻率分別為ω1=3.9和ω2=2.2.采用最佳切換頻率,求得對(duì)應(yīng)于周期軌線X1(t)和X2(t) 的最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g的變化曲線如圖3 所示(藍(lán)色實(shí)線),同時(shí)在表1 中給出了相應(yīng)的穩(wěn)定區(qū)域和穩(wěn)定區(qū)域?qū)挾萕s.

圖3 受控周期軌線 Xk(t) 最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g的變化曲線對(duì)比圖,紅色虛線和藍(lán)色實(shí)線分別對(duì)應(yīng)于DFC和TSDFCFig.3 The comparison of the maximum Lyapunov exponent of periodic trajectory Xk(t) as a function of the feedback strength,the red dotted line and the blue solid line correspond to DFC and TSDFC,respectively

圖4 系統(tǒng)(9)的時(shí)間序列和TSDFC 控制信號(hào)在不同時(shí)段上切換示意圖,在灰色時(shí)段中控制信號(hào)存在,而在灰色時(shí)段外控制信號(hào)消失Fig.4 The time series of system(9) and the schematic diagram of the switching of TSDFC control signals in different time intervals,where the control signal exists in the gray time intervals and disappears outside the gray time intervals

圖5 受控周期軌線 Xk(t) 的穩(wěn)定區(qū)域?qū)挾?Ws隨切換頻率 ω 的變化曲線Fig.5 The stable region width Wsof controlled periodic trajectory Xk(t) as a function of the switching frequencyω

表1 受控周期軌線 Xk(t) 的穩(wěn)定區(qū)域Table 1 The stability regions of the controlled periodic trajectory Xk(t)

通過(guò)圖3和表1 中對(duì)比DFC和TSDFC 的控制效果,可見應(yīng)用TSDFC 時(shí)受控周期軌線Xk(t) 的穩(wěn)定區(qū)域明顯大于應(yīng)用DFC 時(shí)受控周期軌線Xk(t) 的穩(wěn)定區(qū)域.因此,采用恰當(dāng)?shù)那袚Q頻率,TSDFC 的控制性能明顯優(yōu)于DFC 的控制性能.

為進(jìn)一步檢驗(yàn)上述結(jié)果的正確性,考察誤差指數(shù)

式中,X(tj)(j=1,2,···) 是系統(tǒng)狀態(tài)變量的時(shí)間序列,〈·〉表示足夠長(zhǎng)的時(shí)間內(nèi)平均值.當(dāng)受控周期軌線X*(t) 穩(wěn)定時(shí),誤差指數(shù)Sdex趨近于零.誤差指數(shù)Sdex隨反饋增益系數(shù)的變化曲線如圖6 所示.從圖6 中可見,圖中Sdex趨近于零的區(qū)域與圖3和表1 中得到的穩(wěn)定區(qū)域是一致的,從而驗(yàn)證了圖3和表1 中結(jié)果的正確性.

圖6 受控周期軌線 Xk(t) 的誤差指數(shù)隨反饋增益系數(shù) g 的變化曲線對(duì)比圖,紅色虛線和藍(lán)色實(shí)線分別對(duì)應(yīng)于DFC和TSDFCFig.6 The comparison of the error index of the controlled periodic trajectory Xk(t) as a function of the feedback strength g,the red dotted line and the blue solid line correspond to DFC and TSDFC,respectively

2.2 TSDFC 控制Van der Pol 系統(tǒng)

受TSDFC 控制的Van der Pol 系統(tǒng)方程為

式中,x和y為系統(tǒng)狀態(tài)變量,μ為系統(tǒng)參數(shù),a和分別表示外激勵(lì)幅值和頻率.反饋增益系數(shù)gt={sgn[cos(ωt)]+1}·g在定常值g和0之間周期切換.

將系統(tǒng)參數(shù)設(shè)定為 μ=1,a=0.998,=0.45.此時(shí)Van der Pol 系統(tǒng)呈現(xiàn)混沌動(dòng)力學(xué)行為,其混沌吸引子如圖7(a)所示.

選擇Van der Pol 混沌吸引子中一個(gè)典型的不穩(wěn)定周期軌線X3(t)=(x3(t),y3(t))T,作為待鎮(zhèn)定的目標(biāo)周期軌線.圖7(b)展示了目標(biāo)周期軌線X3(t)在xoy平面上的投影.

圖7 Van der Pol 系統(tǒng)的混沌吸引子和不穩(wěn)定周期軌線在 x oy 平面上的投影Fig.7 The projection of the chaotic attractors and unstable periodic trajectory of Van der Pol systems on the x oy plane

為確定受控周期軌線X3(t) 的穩(wěn)定性區(qū)域,求得主穩(wěn)定方程為

式中,時(shí)滯 τ3等于目標(biāo)周期軌線X3(t)的周期T3=2π/.當(dāng)受控周期軌線X3(t) 的最大條件Lyapunov 指數(shù)λm的實(shí)部小于零,受控周期軌線X3(t) 穩(wěn)定.

當(dāng)ω=0時(shí),TSDFC 退化成經(jīng)典DFC.圖8 中紅色虛線是應(yīng)用DFC 鎮(zhèn)定不穩(wěn)定周期軌線X3(t) 時(shí)的最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g變化的曲線圖,同時(shí)求得受控周期軌線X3(t) 的穩(wěn)定區(qū)域?qū)挾萕s為0.88.結(jié)果表明受DFC 控制時(shí)周期軌線X3(t)的穩(wěn)定區(qū)域較小,這表明DFC 的控制效果比較差.

圖8 受控周期軌線 X3(t) 最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g的變化曲線對(duì)比圖,紅色虛線和藍(lán)色實(shí)線分別對(duì)應(yīng)于DFC和TSDFCFig.8 The comparison of the maximum conditional Lyapunov exponent of the controlled periodic trajectory X3(t) as a function of the feedback strength g,the red dotted line and the blue solid line correspond to DFC and TSDFC,respectively

為了擴(kuò)大受控周期軌線X3(t) 穩(wěn)定區(qū)域,令 ω >0,這時(shí)TSDFC 的控制信號(hào)將依據(jù)時(shí)段進(jìn)行切換.為確定最佳切換頻率 ω,對(duì)于周期軌線X3(t),求得穩(wěn)定區(qū)域?qū)挾萕s隨切換頻率 ω 的變化曲線(見圖9).從圖9中可見,切換過(guò)程以非線性方式修改了受控不穩(wěn)定周期軌線的穩(wěn)定性和穩(wěn)定區(qū)域[33],穩(wěn)定區(qū)域?qū)挾萕s隨切換頻率 ω 的增大而非平滑地變化,對(duì)于周期軌線X3(t) 的最佳切換頻率為 ω3=5.4.采用最佳切換頻率 ω3,求得周期軌線X3(t) 的最大條件Lyapunov 指數(shù)隨反饋增益系數(shù)g的變化曲線如圖8 所示(藍(lán)色實(shí)線),同時(shí)也求得受控周期軌線X3(t) 的穩(wěn)定區(qū)域?qū)挾萕s為4.95.

圖9 受控周期軌線 X3(t) 的穩(wěn)定區(qū)域?qū)挾?Ws隨切換頻率 ω 的變化曲線Fig.9 The width of the stable region of the controlled periodic trajectory X3(t) as a function of the switching frequency ω

通過(guò)圖8 中對(duì)比DFC和TSDFC 的控制效果,可見應(yīng)用TSDFC 時(shí)受控周期軌線X3(t) 的穩(wěn)定區(qū)域明顯大于應(yīng)用DFC 時(shí)受控周期軌線X3(t) 的穩(wěn)定區(qū)域.因此采用恰當(dāng)?shù)那袚Q頻率,TSDFC 的控制性能明顯優(yōu)于DFC 的控制性能.

為進(jìn)一步檢驗(yàn)上述結(jié)果的正確性,考察誤差指數(shù)

誤差指數(shù)Sdex隨反饋增益系數(shù)g的變化曲線如圖10所示.從圖10中可見,圖中Sdex趨近于零的區(qū)域與圖8 中得到的穩(wěn)定區(qū)域是一致的,從而驗(yàn)證了圖8 中結(jié)果的正確性.

圖10 受控周期軌線 X3(t) 的誤差指數(shù)隨反饋增益系數(shù) g 的變化曲線對(duì)比圖,紅色虛線和藍(lán)色實(shí)線分別對(duì)應(yīng)于DFC和TSDFCFig.10 The comparison of the error index of the controlled periodic trajectory X3(t) as a function of the feedback strength g,the red dotted line and the blue solid line correspond to DFC and TSDFC,respectively

總結(jié)前面的兩個(gè)算例分析,我們得到以下結(jié)論.TSDFC 的控制性能與切換頻率密切相關(guān).受控周期軌線的穩(wěn)定區(qū)域隨切換頻率的變化呈現(xiàn)強(qiáng)非線性變化.當(dāng)采用恰當(dāng)?shù)那袚Q頻率時(shí),受TSDFC控制時(shí)周期軌線的穩(wěn)定區(qū)域明顯大于受經(jīng)典DFC控制時(shí)周期軌線的穩(wěn)定區(qū)域.

備注:TSDFC 可看成是DFC 的拓展.當(dāng)ω=0時(shí),TSDFC 退化成經(jīng)典DFC.雖然在本文中TSDFC僅用于穩(wěn)定兩個(gè)混沌系統(tǒng)的不穩(wěn)定周期軌線,但與經(jīng)典DFC 一樣,TSDFC 可以應(yīng)用于廣泛的混沌系統(tǒng)控制中.

3 結(jié)論

本文首次在經(jīng)典DFC 的基礎(chǔ)上,加入時(shí)變切換策略來(lái)提高其控制性能,建立了TSDFC 方法.實(shí)例研究表明,當(dāng)選擇恰當(dāng)?shù)那袚Q頻率時(shí),TSDFC 可以穩(wěn)定不穩(wěn)定周期軌線,且受控周期軌線的穩(wěn)定區(qū)域明顯大于應(yīng)用經(jīng)典DFC 時(shí)受控周期軌線的穩(wěn)定區(qū)域.從實(shí)際應(yīng)用的角度來(lái)看,相比經(jīng)典DFC,TSDFC具有顯著的優(yōu)點(diǎn).例如,在一些實(shí)際受控系統(tǒng)中,反饋增益系數(shù)被限制在某些有限的區(qū)域內(nèi).利用TSDFC方法時(shí),可以選擇恰當(dāng)?shù)那袚Q頻率,改變和擴(kuò)大受控周期軌線的穩(wěn)定區(qū)域,使反饋增益系數(shù)落入穩(wěn)定區(qū)域,從而實(shí)現(xiàn)對(duì)周期軌線的鎮(zhèn)定.此外,提出的TSDFC方法,開辟了另一種利用不連續(xù)控制信號(hào)來(lái)控制混沌系統(tǒng)的可能性.TSDFC 的控制性能取決于控制時(shí)段的選擇.本文僅研究了一種比較簡(jiǎn)單的時(shí)變切換策略.對(duì)于其他更復(fù)雜的時(shí)變切換策略,TSDFC 控制效果如何？特別是哪種時(shí)變切換策略對(duì)提高經(jīng)典DFC 的控制性能是最好的,目前還不清楚.這些問(wèn)題值得進(jìn)一步研究.

附錄A 計(jì)算主穩(wěn)定方程的最大條件Lyapunov 指數(shù)

對(duì)于主穩(wěn)定方程

式中,η(t)∈Rm為偏差變量,T為時(shí)滯,X*(t+T)=X*(t)為待鎮(zhèn)定的目標(biāo)周期軌線,DF(t+T,X*(t+T))=DF(t,X*(t))∈Rm×m為系統(tǒng)矩陣.設(shè)時(shí)區(qū)[t0,t0+T] 上的時(shí)間序列為tk=t0+kh,其中k=0,1,···,kT,kT=T/h,h為步長(zhǎng).為減少計(jì)算最大條件Lyapunov指數(shù)的計(jì)算量,用二階Runge-Kutta 方法將主穩(wěn)定方程(A1)離散化得到

在式(A2)中