朱澤飛,涂俐蘭,吳澤虎

(武漢科技大學(xué)理學(xué)院,湖北 武漢 430065)

異維混沌動力系統(tǒng)的有限時間廣義同步

朱澤飛,涂俐蘭,吳澤虎

(武漢科技大學(xué)理學(xué)院,湖北 武漢 430065)

本文研究了異維混沌動力系統(tǒng)的有限時間廣義同步的問題.利用有限時間 Lyapunov 穩(wěn)定性定理、Jensen 不等式等理論方法, 通過設(shè)置不同的控制器,從理論上提出了一般的異維驅(qū)動系統(tǒng)和響應(yīng)系統(tǒng)的有限時間廣義同步的兩種方案,并且對方案二中的影響同步時間因素做了理論分析和證明.最后,數(shù)值模擬驗證了提出理論的正確性和可行性.

異維混沌系統(tǒng);有限時間 Lyapunov 穩(wěn)定性定理;有限時間廣義同步;同步時間

1引言

混沌同步由于在人類大腦[1]、流體混合[2]、無線電通訊[3]等領(lǐng)域中的廣泛應(yīng)用, 而成為非線性科學(xué)研究領(lǐng)域中的熱點課題之一.早期人們在物理學(xué)、生物學(xué)、氣象學(xué)、工程學(xué)等眾多領(lǐng)域中對混沌同步進(jìn)行了深入地研究,得到了一系列實現(xiàn)混沌系統(tǒng)同步的方法,譬如,PC 同步法[4]、脈沖控制同步法[5]、廣義同步法[6]等. 這些方法又可以歸納為兩大類同步 – 完全同步和廣義同步.完全同步的最終目標(biāo)是兩個混沌系統(tǒng)的同步態(tài)完全相同,而廣義同步的同步態(tài)不同,往往呈現(xiàn)出某種函數(shù)關(guān)系.實際上,現(xiàn)實中的系統(tǒng)很難做到動力學(xué)完全相同.

近年來,混沌系統(tǒng)的同步主要集中在研究同維同結(jié)構(gòu)或同維異結(jié)構(gòu)混沌系統(tǒng)之間的漸近同步[4?9]. 然而到目前為止, 對異維異結(jié)構(gòu)混沌系統(tǒng)之間的理論研究結(jié)果還比較少[10?11]. 一方面,當(dāng)兩個混沌系統(tǒng)的維數(shù)不同時,則它們之間的結(jié)構(gòu)必然大相徑庭,而且在相空間中,它們的吸引域也有很大的差異性;另一方面,混沌系統(tǒng)對初值條件極端的敏感,初值條件的任何微小改變,最終必將導(dǎo)致系統(tǒng)之間動力學(xué)行為的巨大變化.所以,相比較同維同結(jié)構(gòu)或同維異結(jié)構(gòu)的混沌系統(tǒng)而言,實現(xiàn)異維異結(jié)構(gòu)混沌系統(tǒng)間的廣義同步就具有更大的挑戰(zhàn)性.

同時,現(xiàn)有的文獻(xiàn)研究的漸近同步是指同步時間趨于無窮大時,混沌系統(tǒng)能否趨于同步狀態(tài).這一要求在實踐中有時候并不現(xiàn)實,譬如,在保密通信中,如果混沌振子在有限時間里不能達(dá)到同步,加密信息不能在有限時間里被成功地恢復(fù)或發(fā)送,都將造成無法挽回的損失[12].

基于以上所述,本文將探討一般的異維混沌系統(tǒng)的有限時間廣義同步.文獻(xiàn) [10]研究了異維混沌系統(tǒng)的廣義同步問題, 但是同步時間仍然是趨于無窮大時的情形; 文獻(xiàn) [13] 研究了隨機(jī)擾動下統(tǒng)一混沌系統(tǒng)的有限時間同步,但是討論的混沌系統(tǒng)仍然是同維混沌系統(tǒng);文獻(xiàn)[11] 研究了異維混沌系統(tǒng)的有限時間廣義同步, 提出了一種控制器的設(shè)計方案, 但是設(shè)計方案還有待提高,而且文中并沒有給出同步時間影響因素的結(jié)論證明. 本文在文獻(xiàn) [11]的基礎(chǔ)之上給出了控制器的兩種設(shè)計方案,其中方案一提出了與文獻(xiàn) [11]的控制器完全不同的設(shè)計思路,而方案二把文獻(xiàn) [11]中定理 1 中的控制器推廣到更一般的情形. 另外,本文把文獻(xiàn) [11]中作者給出的同步時間影響因素的結(jié)論修改得更加合理,而且給出了理論證明.

2預(yù)備定義和引理

考察一般混沌系統(tǒng)的驅(qū)動系統(tǒng)為

響應(yīng)系統(tǒng)為

其中 x=(x1,x2,···,xn)T∈ Rn是驅(qū)動系統(tǒng) (2.1) 的狀態(tài)變量,aij∈ R 是不全為零的常數(shù), fi(x):Rn?→ R(i=1,2,···,n) 是 非 線 性連續(xù)函數(shù);y=(y1,y2,···,ym)T∈ Rm是響應(yīng)系統(tǒng)(2.2) 的 狀態(tài)變量,bij∈ R 是不全為 零 的常數(shù),gj(y):Rm?→ R(j=1,2,···,m) 是非線性連續(xù)函數(shù).

對于響應(yīng)系統(tǒng) (2.2) 式,若施加控制器, 則有

其中 u=(u1,u2,···,um)T∈ Rm是待設(shè)計的控制器.

下面給出有限時間廣義同步的相關(guān)定義以及文中需要用到的引理.

定義1對于驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 中的狀態(tài)變量, 設(shè) E(t)=y(t) ? φ(x(t)),稱 E(t) 為 驅(qū) 動 系 統(tǒng) (2.1) 和 響 應(yīng) 系 統(tǒng) (2.3) 的 廣 義 同 步 誤 差, 其 中 φ(x):Rn?→ Rm, 即φ(x)=(φ1(x),φ2(x),···,φm(x))T∈ Rm是任意一個給定的連續(xù)可微的向量函數(shù).

定義2對于定義 1 中定義的廣義同步誤差 E(t), 如果存在一個時間常量 T ＞ 0, 使得t?→limT?||E(t)||=0, 且當(dāng) t ≥ T 時有 ||E(t)||≡ 0, 則稱驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 在控制器 u 的作用下關(guān)于向量函數(shù) φ(x) 在 T 時刻達(dá)到有限時間廣義同步,其中 T 稱為同步時間, ||.|| 是 2- 范數(shù).

定義3[14]考察下面的非線性動力系統(tǒng)

其中 x ∈ Rn是 系統(tǒng)的狀 態(tài)變量,f(x):Rn?→ Rn是光滑 的 非線性函 數(shù). 如果存在 一 個常數(shù)T ＞ 0(T 依賴于初始條件 x(0) 的值) 使得

且當(dāng) t ≥ T 時有 ||x(t)||≡ 0, 則稱系統(tǒng) (2.4) 是有限時間穩(wěn)定的.

引理1[14](有限時間 Lyapunov 穩(wěn)定性定理) 假設(shè)存在連續(xù)、正定的函數(shù) V(t) 滿足如下微分不等式

其中 c ＞ 0,0 ＜ η＜ 1 為正常數(shù). 那么對于任意給定的 t0,都有

其中

引理2[15](Jensen 不等式) 對于任意的實數(shù) ai(i=1,2,···,n) 以及 0 ＜ p ＜ 1, 有

注1本文的主要目標(biāo)是設(shè)計合適的控制器 u, 在定義 2 的意義下使得驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 達(dá)到有限時間廣義同步.

3主要結(jié)果

3.1控制器的設(shè)計方案

由定義 1,可以給出驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 的廣義同步誤差為 E=y ? φ(x) 且

其中 Dφ(x) 是 φ(x) 的 Jacobi矩陣,即

把 (2.1),(2.3),(3.2) 式代入 (3.1) 式中, 得到驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 的誤差系統(tǒng)為

其中 E=(e1,e2,···,em)T∈ Rm,

根據(jù)定義 3, 對于驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 的有限時間廣義同步問題的研究, 可以等價地轉(zhuǎn)換為研究誤差系統(tǒng) (3.3) 在零點的有限時間穩(wěn)定性問題. 接下來的目標(biāo)是設(shè)計合適的控制器在定義 3 的意義下,誤差系統(tǒng) (3.3) 能夠達(dá)到有限時間穩(wěn)定. 本文設(shè)計了兩種方案來施加控制器u.

注2以下設(shè)文中的 α=qp是合適的有理數(shù),p,q 為正奇數(shù)且 p ＞ q.

方案一分 m 個步驟依次地設(shè)置控制器 u=(u1,u2,···,um)T∈ Rm中的 m 個分量.

第1步對 (3.3) 式的第一個方程設(shè)置控制器

其中 c1,d1為任意給定的正常數(shù).

將 (3.5) 式代入誤差系統(tǒng) (3.3) 中的第一個方程中得

構(gòu)造 Lyapunov 函數(shù)為

對 (3.7) 式關(guān)于時間 t 求導(dǎo)數(shù)得

把 (3.6) 式代入 (3.8) 式得

由引理 1 可知在某個 T1時刻,誤差 e1趨于零,且當(dāng) t≥ T1時,有 e1≡ 0.

第2步類似地, 當(dāng) t ≥ T1時, 由第一步可知 e1≡ 0,于是誤差系統(tǒng) (3.3) 的子系統(tǒng)變?yōu)?/p>

對 (3.10) 式的第一個方程設(shè)置控制器

其中 c2,d2為任意給定的正常數(shù). 將 (3.11) 式代入誤差子系統(tǒng) (3.10) 式中的第一個方程得

類似地,構(gòu)造一個 Lyapunov 函數(shù)為

仍然可以得到

由引理 1 可知在某個 T2＞ T1時刻, 誤差 e2趨于零, 且當(dāng) t ≥ T2時, 有 e2≡ 0.

如此類似地,在上一步的基礎(chǔ)之上,一步一步地設(shè)置控制器 u=(u1,u2,···,um)T∈ Rm中剩下的 m ? 2 個分量 u3,u4,···,um?1,um, 其中

其中 cj,dj(j=3,4,···,m) 為任意給定的正常數(shù). 仍然由引理 1 可知分別在某個 T3＜ T4＜···＜ Tm?1＜ Tm時刻,誤差 ej趨于零,且當(dāng) t≥ Tj時,有 ej≡ 0(j=3,4,···,m). 從而在最后的 Tm時刻有 ‖E(t) ‖ 趨于零, 且當(dāng) t ≥ Tm時有 ‖ E(t) ‖≡ 0. 即誤差系統(tǒng) (3.3) 是有限時間穩(wěn)定的. 也就是說在控制器 u=(u1,u2,···,um)T∈ Rm的作用下,其中驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 關(guān)于向量函數(shù) φ(x) 在 Tm時刻達(dá)到有限時間廣義同步.

方案二一次性地把控制器 u=(u1,u2,···,um)T∈ Rm中的 m 個分量 u1,u2,···,um同時施加上去,只求最后的同步時間.

考察誤差系統(tǒng) (3.3), 若設(shè)施加的控制器 u=(u1,u2,···,um)T∈ Rm中的分量滿足如下條件

其中 pi(i=1,2,···,m) 為任意給定的正常數(shù),那么誤差系統(tǒng) (3.3) 是有限時間穩(wěn)定的.

構(gòu)造 Lyapunov 函數(shù)為

對 (3.18) 式沿著軌線 (3.3) 對時間 t 求導(dǎo)數(shù)得

把 (3.17) 式代入 (3.19) 式中得

由引理2可知

從而

結(jié)合 (3.19)–(3.22) 式可知

由引理 1 可知存在 T ＞ 0 使得在 T 時刻有 ‖ E(t)‖ 趨于零,且當(dāng) t≥ T 時有 ‖E(t) ‖≡ 0.即誤差系統(tǒng) (3.3) 是有限時間穩(wěn)定的. 所以驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 關(guān)于向量函數(shù)φ(x) 在 T 時刻達(dá)到有限時間廣義同步.

注3顯然, 利用方案一在進(jìn)行控制器的設(shè)計時, 誤差系統(tǒng) (3.3) 的子系統(tǒng)變得更加簡單,而且其復(fù)雜程度是低于方案二的.

注4特別地,在方案二的控制器的分量 (3.17) 式中,當(dāng) p1=p2= ···=pm=1 時,此時的控制器就是文獻(xiàn) [11]中定理 1 的結(jié)果. 也就是說, 文獻(xiàn) [11] 中的定理 1 只是本文方案二中(3.17) 式的特殊情形.

3.2同步時間的影響因素分析

在上一節(jié)中用兩種方案設(shè)計出了達(dá)到有限時間廣義同步的控制器,在這一節(jié)中本文將分析方案二中達(dá)到廣義同步所需時間的影響因素.

對于 (3.23) 式,結(jié)合引理 1 可知

把 (3.24) 式代入引理 1 中的 (2.7) 式,整理得

其中 E(t0)=y(t0) ? φ(x(t0)),t0∈ [0,+∞),α ∈ (0,1).

由 (3.25) 式可知驅(qū)動系統(tǒng) (2.1) 和響應(yīng)系統(tǒng) (2.3) 關(guān)于向量函數(shù) φ(x) 的同步時間 T, 可以通過控制初始誤差狀態(tài) E(t0) 和參數(shù) α ∈ (0,1) 來實現(xiàn). 下面討論 E(t0) 與參數(shù) α 分別變化時,對同步時間T 的影響結(jié)果.

注5下文中的e是自然對數(shù)的底數(shù).

定理1假設(shè)給定初始狀態(tài) E(t0),

(1) 若 0 ＜ ||E(t0)||≤ e, 則 T 為 α ∈ (0,1) 上的嚴(yán)格遞增函數(shù);

(2) 若 ||E(t0)||＞ e,

證(1) 由題設(shè)可知 t0和 ||E(t0)|| 都是已知的, 則 T 是關(guān)于 α ∈ (0,1) 上的一元函數(shù),對T關(guān)于α求導(dǎo)數(shù)得

(ii) 若 ?∞ ＜ ln ‖ E(t0) ‖＜ 0, 則從而于是所以 T 為 α ∈ (0,1) 上的嚴(yán)格遞增函數(shù);

綜上 (i),(ii),(iii) 所述,T 為 α ∈ (0,1) 上的嚴(yán)格遞增函數(shù).

上的嚴(yán)格遞增函數(shù);

定理2假設(shè)給定參數(shù) α,則 T 為關(guān)于初始狀態(tài)上的嚴(yán)格遞增函數(shù).

該定理顯然成立,證明從略.

注6定理1從理論上嚴(yán)格地證明了同步時間T和參數(shù) α 之間的關(guān)系.

4數(shù)值模擬

為了證明本文提出方案的有效性和正確性,本節(jié)將對上節(jié)方案二的結(jié)果做數(shù)值模擬,并且只討論驅(qū)動系統(tǒng) (2.1) 的維數(shù) n 大于響應(yīng)系統(tǒng) (2.3) 的維數(shù) m 時的情形.

本節(jié)以三維 Lorenz 系統(tǒng)和二維 Duffi ng 系統(tǒng)分別作為驅(qū)動系統(tǒng)和響應(yīng)系統(tǒng)進(jìn)行驗證. 驅(qū)動系統(tǒng)為

響應(yīng)系統(tǒng)為

由于 φ(x):R3→ R2是任意給定的連續(xù)可微的向量函數(shù),所以為了簡便起見不妨設(shè)

那么

且廣義同步誤差為

由公式 (3.3) 可得誤差系統(tǒng)為

根據(jù) (3.17) 式可得所求的控制器為

因為 p1,p2是任意給定的正常數(shù),所以為了簡便起見,不妨取 p1=p2=1, 那么 (4.7) 式可化為

取驅(qū)動系統(tǒng) (4.1) 和響應(yīng)系統(tǒng) (4.2) 的初始狀態(tài)分別為

下面對參數(shù) α 分別取57和35時做數(shù)值模擬.

利用 MATLAB,可獲得驅(qū)動系統(tǒng)、響應(yīng)系統(tǒng)以及誤差系統(tǒng)的軌跡圖. 圖 1 是驅(qū)動系統(tǒng)(4.1) 的軌跡 圖, 圖 2 是 響 應(yīng)系統(tǒng) (4.2) 的軌 跡 圖, 圖 3 是 當(dāng)參數(shù) α =57時 在控制器 (4.8) 作用下的誤差系統(tǒng)軌跡圖,圖4是當(dāng)參數(shù) α=35時在控制器 (4.8) 作用下的誤差系統(tǒng)軌跡圖.

從模擬圖中可以發(fā)現(xiàn)圖 3 與圖 4 顯示了誤差系統(tǒng)在有限時間收斂到零,也就是說,驅(qū)動系統(tǒng) (4.1) 和響應(yīng)系統(tǒng) (4.2) 在控制器 (4.8) 的作用下, 在有限時間里達(dá)到了廣義同步, 而且圖 3的同步時間比圖 4的同步時間要長一些,這與定理 1的理論結(jié)果是一致的.事實上, 經(jīng) 過計算可知,而又 因 為所以, 于是由定理 1 中的結(jié)論①可知同步時間 T 是關(guān)于參數(shù) α 的嚴(yán)格遞增函數(shù),所以圖3的同步時間比圖4的同步時間要長一些.

5結(jié)論

本文通過設(shè)置不同的控制器,從理論上提出了一般的異維驅(qū)動系統(tǒng)與響應(yīng)系統(tǒng)的有限時間廣義同步的兩種方案,其中方案一給出了與文獻(xiàn) [11]完全不同的設(shè)計思路,方案二把文獻(xiàn)[11]中的結(jié)論推廣到一般情形. 進(jìn)一步地,本文給出了方案二中的參數(shù)和誤差系統(tǒng)的初始狀態(tài)對同步時間影響的理論分析和證明,從而可以通過適當(dāng)?shù)馗淖儏?shù)和系統(tǒng)的初始狀態(tài)來控制同 步速度. 最 后, 利用 三維 的 Lorenz 系 統(tǒng) 和二 維的 Duffi ng 系 統(tǒng) 進(jìn)行 了數(shù)值 仿真 實 驗, 驗證了文中理論的有效性和正確性.

[1]Schiff S J,Jerger K,Duong D H,Chang T,Spano ML,Ditto W L.Controlling chaos in the brain[J]. Nature,1994,370(6491):615–620.

[2]Ottino J M,Muzzio F J,Tjahjadi M,Franjione J G,Jana S C,Kusch H A.Chaos,symmetry,and self-similarity:exploiting order and disorder in mixing processes[J].Sci.,1992,257(5071):754–760.

[3]Argyris A,Syvridis D,Larger L,Annovazzi L V,Colet P,Fischer I,Garcia O J,Mirasso C R, Pesquera L,Shore K A.Chaos-based communications at high bit rates using commercial fi bre-optic links[J].Nature,2005,438(7066):343–346.

[4]Pecora L M,Carroll T L.Synchronization in chaotic systems[J].Phys.Rev.Lett.,1990,64(8): 821–824.

[5]Liu J G.A novel study on the impulsive synchronization of fractional-order chaotic systems[J]. Chinese Phys.B,2013,22(6):285–288.

[6]Rulkov N F,Sushchik M M,Tsimring L S.Generalized synchronization of chaos in directionally coupled chaotic systems[J].Phys.Rev.E,1995,51(2):980–994.

[7] 于娜, 丁群, 陳紅. 異結(jié)構(gòu)系統(tǒng)混沌同步及其在保密通信中的應(yīng)用 [J]. 通信學(xué)報,2007,28(10):73–78.

[8] 李爽, 徐偉, 李瑞紅, 李玉鵬. 異結(jié)構(gòu)系統(tǒng)混沌同步的新方法 [J]. 物理學(xué)報,2006,55(11):5681–5687.

[9]Jia Z,Lu J A,Deng G M.Hybrid projective synchronization of uncertain chaotic systems[J].J. Math.,2011,31(2):275–283.

[10]Ge Z M,Yang C H.The generalized synchronization of a Quantum-CNN chaotic oscillator with diff erent order systems[J].Chaos Sol.Frac.,2008,35(5):980–990.

[11]Cai N,Li W Q,Jing Y W.Finite-time generalized synchronization of chaotic systems with diff erent order[J].Nonl.Dyn.,2011,64(4):385–393.

[12]Bowong S,Kakmeni M,Koina R.Chaos synchronization and duration time of a class of uncertain chaotic systems[J].Math.Comp.Simul.,2006,71(3):212–228.

[13] 王嬌, 涂俐蘭, 朱澤飛. 隨機(jī)擾動下統(tǒng)一混沌系統(tǒng)的有限時間同步 [J]. 數(shù)學(xué)雜志,2014,34(1):1–8.

[14]Wang H,Han Z Z,Xie Q Y,Zhang W.Finite-time chaos control of unifi ed chaotic systems with uncertain parameters[J].Nonl.Dyn.,2009,55(4):323–328.

[15]Gao F,Yuan F.Adaptive fi nite-time stabilization for a class of uncertain high order nonholonomic systems[J].Isa Trans.,2014,54:75–82.

FINITE-TIME GENERALIZED SYNCHRONIZATION OF CHAOTIC DYNAMICAL SYSTEMS WITH DIFFERENT DIMENSIONS

ZHU Ze-fei,TU Li-lan,WU Ze-hu
(School of Science,Wuhan University of Science and Technology,Wuhan 430065,China)

In this paper,fi nite-time generalized synchronization of chaotic dynamical systems with diff erent dimensions is investigated.Based on the theoretical approaches of fi nite-time Lyapunov stability theorem and Jensen inequality etc.,and by means of setting up diff erent controllers,two kinds of schemes are proposed in theory so as to achieve fi nite-time generalized synchronization between general drive system and response system with diff erent dimensions. Furthermore,the factors in the second scheme have been theoretically analyzed and proved which have an impact on the synchronization time.Finally,some numerical simulations are presented to verify that the proposed theories are correct and feasible.

chaotic system with diff erent dimensions;fi nite-time Lyapunov stability theorem;finite-time generalized synchronization;synchronization time

tion:93C15;93D99

3C15;93D99

O231.2

A

0255-7797(2017)02-0365-11

2015-11-16 接收日期:2016-02-26

國家自然科學(xué)基金資助 (61473338);國家自然科學(xué)基金資助 (61304164).

朱澤飛 (1989–), 男, 湖北十堰,碩士, 主要研究方向: 非線性分析與控制