王仲連,郭榮偉
齊魯工業(yè)大學(xué)(山東省科學(xué)院) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,山東 濟(jì)南 250353
混沌就是某些具有確定性的非線(xiàn)性系統(tǒng),在一定參數(shù)范圍內(nèi)能給出非周期的、表面看來(lái)很混亂的輸出,即來(lái)源于確定的體系的無(wú)規(guī)則運(yùn)動(dòng)?;煦绗F(xiàn)象揭示了在確定性和隨機(jī)性之間存在著由此及彼的橋梁,這在科學(xué)觀念上存在著深遠(yuǎn)的意義。越來(lái)越多的學(xué)者從不同的角度來(lái)研究混沌系統(tǒng)的同步問(wèn)題,從兩個(gè)相同系統(tǒng)的同步問(wèn)題推廣到兩個(gè)不同系統(tǒng)的同步問(wèn)題,甚至是兩個(gè)不同維數(shù)的混沌系統(tǒng)之間的同步問(wèn)題[1-4]。這里需要重點(diǎn)指出的是上述研究結(jié)果大部分沒(méi)有考慮系統(tǒng)所受到的外部擾動(dòng),或者考慮的只是在某種意義下的有界擾動(dòng)。事實(shí)上,對(duì)于一個(gè)系統(tǒng)來(lái)說(shuō),外部擾動(dòng)是不可能完全避免的,并且有時(shí)系統(tǒng)所受到的擾動(dòng)并不是有界的,比如在某一時(shí)刻系統(tǒng)受到了外部脈沖信號(hào)的干擾。因此,研究系統(tǒng)在外部干擾下的控制問(wèn)題不僅具有重要的理論意義,而且具有重大的應(yīng)用價(jià)值。UDE方法[5-7]是近年來(lái)處理系統(tǒng)中模型不確定性和外部擾動(dòng)的一個(gè)非常好的方法,對(duì)于一些有界的外部擾動(dòng),該方法能實(shí)現(xiàn)精確的估計(jì),從而實(shí)現(xiàn)精確跟蹤或者漸近穩(wěn)定。但是,對(duì)于某種類(lèi)型的無(wú)界擾動(dòng),該方法只能給出魯棒性的結(jié)果。受UDE思想的啟發(fā),通過(guò)尋找合適的函數(shù)進(jìn)行卷積的方法對(duì)兩種類(lèi)型的擾動(dòng)分別設(shè)計(jì)了相應(yīng)的濾波器,均實(shí)現(xiàn)了控制目標(biāo)。
四維超混沌金融系統(tǒng)[8]是一個(gè)的重要的經(jīng)濟(jì)學(xué)模型,其系統(tǒng)模型有四個(gè)變量,即貨幣、生產(chǎn)、勞動(dòng)力和股票。該模型自提出以來(lái),已經(jīng)取得了大量的研究成果[9-12]。然而,正如上面所介紹的,外部干擾的影響,特別是在無(wú)界的外部干擾的影響,并沒(méi)有被考慮。 因此,研究在無(wú)界的外部擾動(dòng)下的四維金融超混沌系統(tǒng)的完全同步問(wèn)題是非常重要和迫切的。
定義1[1]:考慮如下的混沌系統(tǒng)
其中,x∈Rn是系統(tǒng)的狀態(tài),g(x)是連續(xù)的向量函數(shù)。
設(shè)系統(tǒng)(1)是主系統(tǒng),則相應(yīng)的從系統(tǒng)為
其中,y∈Rn是系統(tǒng)的狀態(tài),g(y)是連續(xù)的向量函數(shù),B∈Rn×r是常矩陣,r≥1,u∈Rr是帶設(shè)計(jì)的控制器。
再令E=y-x,則誤差系統(tǒng)為
考慮在外部擾動(dòng)干擾下受控的四維超混沌金融系統(tǒng)[8]
其中,x∈R4是系統(tǒng)的狀態(tài),f(x)∈R4是連續(xù)的向量函數(shù),d(t)∈R4是外部擾動(dòng),即
外部擾動(dòng)d(t)有兩種情況:
情況1:
情況2:
其中A,C是未知非零常數(shù)。
令系統(tǒng)(4)為主系統(tǒng),則從系統(tǒng)為
其中,y∈R4是系統(tǒng)的狀態(tài),f(y)∈R4是連續(xù)的向量函數(shù)。
本文的主要目的是設(shè)計(jì)控制器u∈R使得
即主系統(tǒng)(4)和從系統(tǒng)(8)達(dá)到了完全同步。
為了實(shí)現(xiàn)系統(tǒng)對(duì)外部干擾的抑制,受UDE思想的啟發(fā),對(duì)常數(shù)類(lèi)型的有界擾動(dòng)和多項(xiàng)式乘以三角函數(shù)類(lèi)型的無(wú)界擾動(dòng)分別設(shè)計(jì)了相應(yīng)的濾波器,并實(shí)現(xiàn)了對(duì)上述干擾的精確抑制。
定理1 如果常數(shù)類(lèi)型的有界擾動(dòng)d1(t)=C,其中C是未知常數(shù),則設(shè)計(jì)的濾波器為
其中,k是一個(gè)正數(shù)。
證明:當(dāng)t→∞時(shí)
定理2 如果無(wú)界的外部擾動(dòng)為d1(t)=Atsin(t),其中A是未知常數(shù),則設(shè)計(jì)的濾波器為
證明:因?yàn)?/p>
而計(jì)算如下卷積得到
定理3 考慮主系統(tǒng)(4)和從系統(tǒng)(8)對(duì)于公式(6)所表示的有界擾動(dòng),設(shè)計(jì)的同步控制器如下:
u=us+ud, (11)
其中,
us=k(t)(1,0,0,0)E=k(t)E1, (12)
E=x-y,Gf(s)見(jiàn)公式(9),
證明:令E=x-y,則誤差系統(tǒng)為
其中,
注意到,若E1=0,則下面的系統(tǒng)
是漸近穩(wěn)定的。
再對(duì)下面的方程兩邊取Laplace變換
整理可得到公式(13)中的ud。
對(duì)于第二種情況的擾動(dòng),證明過(guò)程類(lèi)似,本文直接給出如下的結(jié)論。
定理4 考慮主系統(tǒng)(4)和從系統(tǒng)(8)對(duì)于公式(7)所表示的無(wú)界擾動(dòng),設(shè)計(jì)的同步控制器為
u=us+ud, (18)
其中us見(jiàn)公式(12),
E=x-y,Gf(s)見(jiàn)公式(10),k(t),F1(y,E)見(jiàn)公式(14)。
情況1:公式(6)所示的有界擾動(dòng)。
選擇主系統(tǒng)(4)的初值為x(0)=[1,2,-3,4],從系統(tǒng)(8)的初值為y(0)=[5,-2,3,2],反饋增益的初值k(0)=-1,d1(t)=100。MATLAB仿真結(jié)果見(jiàn)圖1、圖2、圖3。
圖1 誤差系統(tǒng)漸近穩(wěn)定
圖2 動(dòng)態(tài)反饋增益k(t)漸近收斂到一個(gè)合適的常數(shù)
圖漸近收斂到d1(t)
情況2:公式(7)所示的無(wú)界擾動(dòng)。
選擇主系統(tǒng)(4)的初值為x(0)=[3,2,5,2],從系統(tǒng)(8)的初值為y(0)=[3,2,5,2],反饋增益的初值k(0)=-1,擾動(dòng)d1(t)=10tsint。MATLAB仿真結(jié)果見(jiàn)圖4、圖5、圖6。
圖4 誤差系統(tǒng)漸近穩(wěn)定
圖5 動(dòng)態(tài)反饋增益k(t)漸近收斂到一個(gè)合適的常數(shù)
圖漸近收斂到d1(t)
本文用基于擾動(dòng)估計(jì)器的控制方法研究了四維混沌金融系統(tǒng)的完全同步問(wèn)題。受UDE思想的啟發(fā),對(duì)于常數(shù)類(lèi)型的有界外部擾動(dòng)和多項(xiàng)式乘以三角函數(shù)類(lèi)型的無(wú)界外部擾動(dòng),分別設(shè)計(jì)了相應(yīng)的濾波器,并設(shè)計(jì)了相應(yīng)的控制器,從而實(shí)現(xiàn)了該系統(tǒng)的完全同步問(wèn)題。首次解決了在某種無(wú)界的外部擾動(dòng)下四維超混沌金融系統(tǒng)的完全同步問(wèn)題。數(shù)值仿真結(jié)果驗(yàn)證了理論結(jié)果的正確性和有效性。
獲得的理論結(jié)果經(jīng)過(guò)適當(dāng)?shù)匦薷目赏茝V到研究一般的非線(xiàn)性系統(tǒng)的干擾抑制問(wèn)題,為研究非線(xiàn)性系統(tǒng)的各種控制問(wèn)題提供了理論支持。