徐昌彪, 鐘 德, 夏 誠, 黎 周
(1.重慶郵電大學(xué) 光電工程學(xué)院,重慶 400065; 2.重慶郵電大學(xué) 通信與信息工程學(xué)院,重慶 400065)
混沌是確定性系統(tǒng)的內(nèi)在隨機(jī)性,是非線性動(dòng)力學(xué)系統(tǒng)所特有的一種形式,并且廣泛地存在于自然界諸如物理[1]、化學(xué)[2]、地質(zhì)學(xué)[3]以及生物學(xué)[4]等領(lǐng)域。自從1963年Lorenz[5]在三維自治系統(tǒng)中發(fā)現(xiàn)了第一個(gè)混沌吸引子以來,人們不斷地發(fā)現(xiàn)新的混沌系統(tǒng),如Chen等[6]系統(tǒng)、Lü等[7]系統(tǒng)、Chua等[8]系統(tǒng)、Liu等[9]系統(tǒng)等。但大多數(shù)混沌系統(tǒng)的參數(shù)變化范圍較小[10-14],它們的混沌特性受到系統(tǒng)任意一個(gè)參數(shù)很大影響,參數(shù)的變更與誤差會(huì)改變系統(tǒng)的動(dòng)力學(xué)行為,使得系統(tǒng)相軌收斂于不動(dòng)點(diǎn),或是處于周期、擬周期甚至混沌振蕩狀態(tài)等。而具有超大范圍混沌系統(tǒng)的混沌特性在某個(gè)參數(shù)的變動(dòng)下始終處于混沌狀態(tài),其最大Lyapunov指數(shù)始終大于零,這使得該系統(tǒng)具有更豐富的動(dòng)力學(xué)特性,在混沌雷達(dá)[15],激光振蕩[16]以及信息處理[17]等領(lǐng)域中具有重大的應(yīng)用價(jià)值。
但總的來說,關(guān)于超大范圍混沌系統(tǒng)的研究還處于起步階段,Liu等[18]提出了一個(gè)具有較大參數(shù)范圍的新型三維混沌系統(tǒng),參數(shù)范圍為[-100,100]。賈紅艷等[19]在已有的超混沌系統(tǒng)研究的基礎(chǔ)上提出了一個(gè)大范圍超混沌系統(tǒng),參數(shù)范圍為[10,1 000],梅蓉等[20]構(gòu)建了一類新的超大范圍超混沌系統(tǒng),參數(shù)范圍為[0,105],但仿真只在[0,2 000]進(jìn)行。
基于以上研究,本文提出了一個(gè)超大范圍的混沌系統(tǒng),其中參數(shù)b的取值為[0,107],比本文中提到的其他系統(tǒng)的參數(shù)范圍都要大。理論分析了系統(tǒng)的動(dòng)力學(xué)特性,進(jìn)行了數(shù)值仿真和電路仿真。設(shè)計(jì)了一個(gè)自適應(yīng)控制器和一個(gè)自適應(yīng)滑??刂破?,分別用于具有未知參數(shù)混沌系統(tǒng)的全局穩(wěn)定和給定信號(hào)的追蹤與未知參數(shù)的辨識(shí)。仿真結(jié)果表明了所設(shè)計(jì)控制器的可行性和有效性。
新混沌系統(tǒng)的數(shù)學(xué)模型為
(1)
式中:x,y,z為系統(tǒng)變量。取a=3,b=10,初值為[1,1,1]時(shí),系統(tǒng)存在一個(gè)典型的吸引子,如圖1所示。此時(shí)系統(tǒng)的三個(gè)Lyapunov指數(shù)為0.122 7,-0.001 7,-3.121 0。
(a) x-y相圖(b) x-z相圖(c) y-z相圖(d) x-y-z相圖
圖1 系統(tǒng)相圖
Fig.1 Phase diagram of the system
系統(tǒng)在變換(x,y,z)→(-x,-y,z)下具有不變性,即系統(tǒng)關(guān)于z軸對(duì)稱,且這種對(duì)稱性對(duì)所有系統(tǒng)參數(shù)均成立。
根據(jù)式(1),有
(2)
可見系統(tǒng)是耗散的,則系統(tǒng)的軌線最終會(huì)被限制在一個(gè)體積為0的極限點(diǎn)集上,并且它的漸近動(dòng)力學(xué)行為會(huì)被固定在一個(gè)吸引子上,這說明了吸引子的存在性。
令式(1)的左邊等于0,即
(3)
其中a=3,b=10時(shí),系統(tǒng)有3個(gè)平衡點(diǎn),分別為:
S1=(0,0,0),
S2=(-0.223 6,0.223 6,1),
S3=(0.223 6,-0.223 6,1)。
在平衡點(diǎn)S1=(0,0,0)處線性化系統(tǒng),得其Jacobi矩陣為
(4)
特征值方程為
|λI-J|=0
(5)
計(jì)算出特征值為λ1=-3,λ2=1,λ3=-1。由于特征值λ1和λ3是負(fù)實(shí)數(shù),λ2是正實(shí)數(shù),故平衡點(diǎn)S1是不穩(wěn)定的鞍點(diǎn)。
同理,在平衡點(diǎn)S2和S3處分別線性化系統(tǒng),得到它們的特征值均為λ1=-3.434 8,λ2=0.217 4+1.856 4i,λ3=0.217 4-1.856 4i。由于特征值λ1是負(fù)實(shí)數(shù),λ2和λ3均有正實(shí)部,故平衡點(diǎn)S2和S3都是不穩(wěn)定的鞍點(diǎn)。
取a=3,b=10,初值為[1,1,1], 采用四階龍格-庫塔(ODE45)算法對(duì)式(1)進(jìn)行求解,可得系統(tǒng)狀態(tài)變量x的時(shí)域波形圖,如圖2所示,可以看出系統(tǒng)為非周期系統(tǒng)。如圖3所示,系統(tǒng)的功率譜是連續(xù)譜,沒有明顯的波峰,并且峰值連成一片,說明了系統(tǒng)是混沌系統(tǒng)。選取的Poincare截面如圖4所示,可以看出Poincare截面上有一些成片的具有分形結(jié)構(gòu)的密集點(diǎn),吸引子的葉片清晰可見,也表明了系統(tǒng)具有混沌特性。
圖2 狀態(tài)變量的時(shí)域波形
圖3 系統(tǒng)的功率譜
圖4 系統(tǒng)的Poincare截面:y=0
采用Jacobi法計(jì)算出系統(tǒng)的Lyapunov指數(shù)譜如圖5所示,在參數(shù)b從0~107變化時(shí),系統(tǒng)的三個(gè)Lyapunov指數(shù)都趨于固定常數(shù),為恒Lyapunov指數(shù)譜。圖6為系統(tǒng)變量y隨參數(shù)b從0~107變化時(shí)的分岔圖,可以看出系統(tǒng)都處于混沌狀態(tài)。
系統(tǒng)的Lyapunov維數(shù)為:
(6)
當(dāng)a=3,b=10時(shí),系統(tǒng)的三個(gè)Lyapunov指數(shù)分別為L1=0.122 7,L2=-0.001 7,L3=-3.121 0,故可計(jì)算出DL=2.038 8。由此可見,系統(tǒng)的Lyapunov維數(shù)為分?jǐn)?shù)維數(shù),驗(yàn)證了系統(tǒng)為混沌系統(tǒng)。
采用線性電阻、電容、LM2924N運(yùn)算放大器、MULTIPLIER模擬乘法器(乘法器的輸出增益為0.1),設(shè)計(jì)出系統(tǒng)的模擬電子電路,如圖7所示。
圖5 系統(tǒng)的Lyapunov指數(shù)譜Fig.5 Lyapunov exponent spectrum of the system圖6 系統(tǒng)的分岔圖Fig.6 Bifurcation diagram of the system
圖7 電路原理圖
根據(jù)系統(tǒng)參數(shù)以及電路理論,由乘法器的輸出增益為0.1,得電路方程為:
(7)
當(dāng)取R1=R6=R9=10 kΩ,R3=R8=R11=100 kΩ,R2=3 kΩ,R4=R5=R7=R10=R12=R13=R14=R15=1 kΩ,R16=R17=0.1 kΩ,C1=C2=C3=1 μF時(shí),采用Multisim軟件對(duì)電路進(jìn)行了仿真實(shí)驗(yàn),仿真結(jié)果如圖8所示,可以看出仿真實(shí)驗(yàn)結(jié)果與數(shù)值仿真結(jié)果是一致的。
(a) x-y
(b) x-z
(c) y-z
受控系統(tǒng)為
(8)
式中:x,y,z為系統(tǒng)變量,u1,u2,u3為自適應(yīng)控制器。
設(shè)計(jì)自適應(yīng)控制器為
(9)
參數(shù)估計(jì)誤差定義為
(10)
對(duì)式(8)求導(dǎo),得
(11)
將式(9)和式(10)代入式(8),得閉環(huán)系統(tǒng)
(12)
構(gòu)造以下Lyapunov函數(shù)
(13)
其中λ1,λ2均大于零,顯然V是一個(gè)正定函數(shù)。對(duì)V求導(dǎo)可得
(14)
根據(jù)式(14),設(shè)定參數(shù)更新定律如下
(15)
定理1具有未知參數(shù)的受控系統(tǒng)(8)通過自適應(yīng)控制器(9)和參數(shù)更新定律(15)在所有初始條件下全局穩(wěn)定,其中k1,k2,k3是正增益常數(shù)。
(16)
顯然,式(16)是半負(fù)定的,狀態(tài)變量x,y,z和參數(shù)估計(jì)誤差ea(t),eb(t)是全局有界的。
根據(jù)式(16)可得
(17)
因此,可得:
(18)
將不等式(18)兩端從0到t積分,得到
(19)
根據(jù)式(19)可得x,y,z,ea,eb∈L2。
此外,由圖9,10可知,在自適應(yīng)控制器(9)和參數(shù)更新定律(15)的作用下受控系統(tǒng)(8)很快穩(wěn)定到平衡點(diǎn)原點(diǎn),并且實(shí)現(xiàn)了對(duì)系統(tǒng)未知參數(shù)的識(shí)別。
圖9 受控狀態(tài)x,y,z
考慮如下受控系統(tǒng)
圖10 參數(shù)和的收斂曲線
(20)
定義滑模函數(shù)
(21)
式中:r為理想位置信號(hào),e=y-r為跟蹤誤差,k>0。
設(shè)計(jì)自適應(yīng)控制器為
(22)
(23)
對(duì)參數(shù)估計(jì)誤差(23)求導(dǎo)得到
(24)
定義Lyapunov函數(shù)為
(25)
對(duì)V求導(dǎo),將式(22)和式(23)代入,得
(26)
取參數(shù)更新定律為
(27)
則
(28)
根據(jù)Lyapunov穩(wěn)定性定理,可知受控系統(tǒng)式(20)的平衡狀態(tài)是全局漸進(jìn)穩(wěn)定的。
圖11 跟蹤sin(t)的曲線
圖12 參數(shù)和的收斂曲線
本文提出了一個(gè)超大范圍的混沌系統(tǒng),分析了系統(tǒng)的動(dòng)力學(xué)特性,設(shè)計(jì)了系統(tǒng)的模擬電路。最后設(shè)計(jì)了一個(gè)自適應(yīng)控制器和一個(gè)自適應(yīng)滑模控制器,分別用于具有未知參數(shù)的混沌系統(tǒng)的全局穩(wěn)定和給定信號(hào)的追蹤與未知參數(shù)的辨識(shí),數(shù)值仿真結(jié)果表明所設(shè)計(jì)控制器是有效的。
圖13 跟蹤cos(t)的曲線
圖14 參數(shù)和的收斂曲線