分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源反步控制器設(shè)計(jì)

趙飛，鄭永愛(ài)

（揚(yáng)州大學(xué)信息工程學(xué)院，江蘇 揚(yáng)州 225127）

最近，學(xué)者們提出了多種混沌系統(tǒng)的控制和同步方法[1-10]。其中，反步控制是一種將李雅普諾夫函數(shù)的選擇與反饋控制器的設(shè)計(jì)相結(jié)合的系統(tǒng)控制方法。文獻(xiàn)[11-13]利用反步控制法實(shí)現(xiàn)了整數(shù)階和分?jǐn)?shù)階混沌系統(tǒng)的控制和同步。無(wú)源控制是從系統(tǒng)的能量角度入手，使閉環(huán)系統(tǒng)滿足無(wú)源性的一種控制方法，從而實(shí)現(xiàn)系統(tǒng)的穩(wěn)定性。利用無(wú)源控制，文獻(xiàn)[14]和文獻(xiàn)[15-16]分別實(shí)現(xiàn)了超混沌系統(tǒng)和分?jǐn)?shù)階系統(tǒng)的控制和同步。結(jié)合反步控制和無(wú)源控制技巧，該文提出了實(shí)現(xiàn)分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源反步控制法。

1 預(yù)備知識(shí)

Caputo 微分定義為[17]：

其中，n-1 ＜α＜n,n∈N，Γ(·)為伽馬函數(shù)。非自治分?jǐn)?shù)階系統(tǒng)如式（2）：

其中，α∈(0,1)是分?jǐn)?shù)階，x=(x1,x2,…,xn)T是系統(tǒng)的狀態(tài)變量，f=(f1,f2,…,fn)T是可微的非線性向量函數(shù)。

引理1[18]設(shè)x=0 是式（2）系統(tǒng)的平衡點(diǎn)，若存在李雅普諾夫函數(shù)V(x(t),t) 和K 類(lèi)函數(shù)γi(i=1,2,3)滿足式（3）：

那么式（2）系統(tǒng)是漸近穩(wěn)定的。

引理2[19]若x(t)為連續(xù)可微函數(shù)向量，則有：

分?jǐn)?shù)階系統(tǒng)如式（5）：

其中，x∈Rn是狀態(tài)變量，u∈Rm是輸入變量，y∈Rm是輸出變量，f(x)和g(x)是光滑向量場(chǎng)，h(x)是光滑映射，且f(0)=h(0)=0。

假設(shè)式（5）系統(tǒng)可以用式（6）表示：

其中，(zT,yT)T∈Rn是系統(tǒng)的新坐標(biāo)，z∈Rn-m是狀態(tài)變量。

定義1如果存在一個(gè)正定李雅普諾夫函數(shù)V(z,y)(稱為存儲(chǔ)函數(shù)或能量函數(shù))，滿足對(duì)任意的t≥0有≤uT y，則式（6）系統(tǒng)是無(wú)源的。如果存在一個(gè)正定李雅普諾夫函數(shù)V(z,y)和一個(gè)K 類(lèi)函數(shù)γ4，滿足對(duì)任意的t≥0 有：

則式（6）系統(tǒng)是嚴(yán)格無(wú)源的。

由引理1 易得下面的引理3：

引理3假設(shè)式（6）系統(tǒng)是嚴(yán)格無(wú)源的，則式（6）系統(tǒng)的平衡點(diǎn)x=0 在輸入u=0 的情況下是漸近穩(wěn)定的。

2 分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源反步控制設(shè)計(jì)

將利用反步控制研究分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源控制問(wèn)題，系統(tǒng)描述如下：

其中，α是分?jǐn)?shù)階階次，z∈Rn-r，(x1,x2,…,xr)T∈Rr是狀態(tài)，u∈R是輸入，函數(shù)f(·)、g(·)和φi(·)(i=1,2…,r)都是光滑函數(shù)，且f(0,0)=0,g(0,0)=0,φi(0,0)=0(i=1,2,…,r)。

定理1若函數(shù)f(z,0)滿足zT f(z,0)≤-|z|2l(z)，其中l(wèi)(z)是光滑函數(shù)，且對(duì)于任意的z，都有l(wèi)(z)＞0，則存在如下反饋控制器（式（9）），使得式（8）系統(tǒng)是嚴(yán)格無(wú)源的，且在外部輸入信號(hào)v=0 的情況下是漸近穩(wěn)定的。

其中，α1=-e1-zT g(z,x1)-φ1(z,x1)與αi=-ei-ei-1-φi(z,x1,…,xm)+(i=2,…,r-1) 為設(shè)計(jì)的虛擬控制器，且轉(zhuǎn)換變量e1=x1、ei=xi-αi-1(i=2,…,r)。

證明設(shè)計(jì)步驟共有r步，第1 步：設(shè)x2為控制器，α1為虛擬控制器。由于e1=x1、e2=x2-α1，則有子系統(tǒng)：

第i(2 ≤i＜r)步：設(shè)xi+1為實(shí)際控制器，αi是虛擬控制器。由于ei+1=xi+1-αi，則有子系統(tǒng)：

從而由定義1 知，式（12）子系統(tǒng)是嚴(yán)格無(wú)源的。

第r步：式（8）系統(tǒng)變?yōu)椋?/p>

所以由定義1 知，式（14）系統(tǒng)是嚴(yán)格無(wú)源的。即式（8）系統(tǒng)是嚴(yán)格無(wú)源的。進(jìn)一步由引理3 知，在輸入v=0 的情況下，式（8）系統(tǒng)也是漸近穩(wěn)定的。

3 數(shù)值仿真

例1 考慮如下分?jǐn)?shù)階Shimizu-Morioka 混沌系統(tǒng)[20]：

其中，x、y、z是狀態(tài)變量，u∈R是輸入，0＜α≤1是分?jǐn)?shù)階階次，a∈R、b∈R+是參數(shù)。當(dāng)α=0.98、a=0.75、b=0.45，初始條件為x=0、y=0.25、z=1 時(shí)，系統(tǒng)呈現(xiàn)混沌現(xiàn)象。分?jǐn)?shù)階Shimizu-Morioka系統(tǒng)的混沌吸引子如圖1 所示。

圖1 分?jǐn)?shù)階Shimizu-Morioka系統(tǒng)的混沌吸引子

設(shè)z=z、x1=x、x2=y，式（16）系統(tǒng)可表示為如式（8）的形式，如式（17）：

對(duì)于式（17）系統(tǒng)，n=3,r=2,f(z,0)=-bz,zf(z,0)=-bz2,l(z)=b＞0 滿足定理1 的條件。根據(jù)定理1，在虛擬控制律α1=-e1-ze1與控制律u=-2e1+ae2+aα1+e1z+CDα t α1-e2+v作用下，式（17）系統(tǒng)是嚴(yán)格無(wú)源的，且在輸入v=0 的情況下是漸近穩(wěn)定的。

對(duì)于式（17）系統(tǒng)，當(dāng)α=0.98、a=0.75、b=0.45，初始條件為z=-2、x1=-0.2、x2=-6 時(shí)，受控分?jǐn)?shù)階Shimizu-Moriok 系統(tǒng)的狀態(tài)響應(yīng)如圖2 所示，結(jié)果表明，所設(shè)計(jì)的控制器能將分?jǐn)?shù)階Shimizu-Moriok 系統(tǒng)穩(wěn)定到原點(diǎn)。

圖2 受控Shimizu-Moriok系統(tǒng)的狀態(tài)響應(yīng)

例2分?jǐn)?shù)階Lü混沌系統(tǒng)[21]如式（18）：

其中，x、y、z是狀態(tài)變量，u∈R是輸入，0 ＜α≤1是分?jǐn)?shù)階階次，a、b、c是正常參數(shù)。當(dāng)α=0.96、a=36、b=3、c=20，初始條件為x=1、y=2、z=2 時(shí)，系統(tǒng)呈現(xiàn)混沌現(xiàn)象。分?jǐn)?shù)階Lü系統(tǒng)的混沌吸引子如圖3 所示。

圖3 分?jǐn)?shù)階Lü系統(tǒng)的混沌吸引子

設(shè)z1=x,z2=z,x1=y，式（18）系統(tǒng)可表示為如式（8）的形式，如式（19）：

對(duì)于式（19）系統(tǒng)，n=3,r=1,f(z,0)=(-az1,-bz2)T,(z1,z2)f(z,0)=≤-max{a,b}||z||2,l(z)=max{a,b}＞0 滿足定理1 的條件。根據(jù)定理1，在控制律u=-ce1-az1-e1+v作用下，式（19）系統(tǒng)是嚴(yán)格無(wú)源的。且在輸入v=0 的情況下是漸近穩(wěn)定的。

對(duì)于式（19）系統(tǒng)，當(dāng)α=0.96、a=36、b=3、c=20，初始條件為z1=-22、z2=1、x1=16 時(shí)，受控分?jǐn)?shù)階Lü系統(tǒng)的狀態(tài)響應(yīng)如圖4 所示，結(jié)果表明，所設(shè)計(jì)的控制器能將分?jǐn)?shù)階Lü系統(tǒng)穩(wěn)定到原點(diǎn)。

圖4 受控分?jǐn)?shù)階Lü系統(tǒng)的狀態(tài)響應(yīng)

針對(duì)下面的整數(shù)階Lü混沌系統(tǒng)（式（20）），通過(guò)將該文所提方法與文獻(xiàn)[11]的反步控制方法做對(duì)比，進(jìn)一步來(lái)說(shuō)明該文所提方法的優(yōu)越性和研究動(dòng)機(jī)。

該文和文獻(xiàn)[11]設(shè)計(jì)的控制器分別為u=-cx1-az1-x1和u=-(1+c)x1-az1+z1z2，很明顯該文設(shè)計(jì)的控制器是線性控制器，但也具有和文獻(xiàn)[11]一樣的控制效果，所以該文方法更簡(jiǎn)單有效。當(dāng)a=36、b=3、c=20，初始條件為x1=10、x2=10、z=-10 時(shí)，受控整數(shù)階Lü系統(tǒng)分別在該文和文獻(xiàn)[11]所設(shè)計(jì)的控制器下的狀態(tài)響應(yīng)如圖5 和圖6 所示。

圖5 在該文控制器下的狀態(tài)響應(yīng)

圖6 在文獻(xiàn)[11]控制器下的狀態(tài)響應(yīng)

4 結(jié)論

該文研究了分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源反步控制問(wèn)題，利用反步控制的遞推過(guò)程，獨(dú)立設(shè)計(jì)每個(gè)子系統(tǒng)的能量函數(shù)和虛擬控制律，最終設(shè)計(jì)出整個(gè)系統(tǒng)的實(shí)際控制律，實(shí)現(xiàn)對(duì)分?jǐn)?shù)階混沌系統(tǒng)的無(wú)源反步控制。以分?jǐn)?shù)階Shimizu-Morioka 系統(tǒng)和分?jǐn)?shù)階Lü系統(tǒng)為例進(jìn)行數(shù)值仿真，理論分析和數(shù)值仿真結(jié)果驗(yàn)證了該方法的有效性。