王 琳，張慶杰，陳宏偉

(空軍航空大學(xué)， 吉林 長(zhǎng)春 130022)

群系統(tǒng)編隊(duì)控制問題[1-2]越來越受到學(xué)者的關(guān)注。傳統(tǒng)的編隊(duì)控制方法存在通信代價(jià)大，以及易單點(diǎn)失效導(dǎo)致整體癱瘓的缺點(diǎn)。伴隨一致性理論的出現(xiàn)，為編隊(duì)控制問題提供了新的解決思路。文獻(xiàn)[3]證明了一致性編隊(duì)控制方法可以解決編隊(duì)控制問題，且優(yōu)于傳統(tǒng)編隊(duì)控制方法。

文獻(xiàn)[4]解決了一階積分特性的群系統(tǒng)模型下的編隊(duì)跟蹤問題。文獻(xiàn)[5]解決了二階積分特性模型下的編隊(duì)跟蹤問題。文獻(xiàn)[6]在假設(shè)含有有向生成樹，且領(lǐng)導(dǎo)者須作為根節(jié)點(diǎn)的條件下，解決了二階積分特性模型下的編隊(duì)跟蹤問題。文獻(xiàn)[7]在通信拓?fù)鋱D為無向圖時(shí)，解決了編隊(duì)控制問題。文獻(xiàn)[8]在通信拓?fù)鋱D為有向圖時(shí)，解決了時(shí)變編隊(duì)問題。文獻(xiàn)[9]在通信拓?fù)鋱D為有向圖時(shí)，利用一致性方法，解決了無人機(jī)編隊(duì)問題。文獻(xiàn)[10]解決了高階積分特性的群系統(tǒng)模型下的編隊(duì)跟蹤問題。文獻(xiàn)[11]解決了最小化線性二次型調(diào)節(jié)器(linear quadratic regulator,LQR)性能指標(biāo)的一階積分特性模型下的編隊(duì)問題。文獻(xiàn)[12]研究了系統(tǒng)矩陣正定時(shí)領(lǐng)導(dǎo)者跟隨者系統(tǒng)跟蹤最優(yōu)控制。文獻(xiàn)[13]利用逆優(yōu)化方法設(shè)計(jì)時(shí)變編隊(duì)跟蹤協(xié)議，給出群系統(tǒng)實(shí)現(xiàn)時(shí)變編隊(duì)跟蹤的可行性條件的證明，并保證所提出的控制策略滿足LQR性能指標(biāo)。

文獻(xiàn)[4-6]的控制方法只適用于低階積分特性模型，模型條件過于苛刻。文獻(xiàn)[6]的控制方法針對(duì)的是領(lǐng)導(dǎo)者跟隨者系統(tǒng)，其拓?fù)錀l件要求含有有向生成樹且領(lǐng)導(dǎo)者須作為根節(jié)點(diǎn)，適用范圍較窄。文獻(xiàn)[7]要求通信拓?fù)鋱D切換時(shí)均為無向圖，針對(duì)有向拓?fù)錀l件，文獻(xiàn)[7]的結(jié)論并不適用。文獻(xiàn)[8-9]的通信拓?fù)鋱D雖為有向圖，但是解決的是編隊(duì)形成問題，其結(jié)論不適用于編隊(duì)跟蹤問題。文獻(xiàn)[10]考慮了高階積分特性的編隊(duì)跟蹤問題，但是沒有涉及編隊(duì)跟蹤性能方面。文獻(xiàn)[11]考慮了最小化LQR性能指標(biāo)的編隊(duì)控制問題，針對(duì)的是低階模型，且沒有涉及跟蹤軌跡方面。文獻(xiàn)[12]對(duì)領(lǐng)導(dǎo)者跟隨者系統(tǒng)達(dá)到跟蹤最優(yōu)控制的條件較為苛刻，要求系統(tǒng)矩陣正定，且不考慮時(shí)變編隊(duì)和軌跡跟蹤，應(yīng)用范圍較窄。文獻(xiàn)[13]雖然在通信拓?fù)鋱D為有向圖時(shí)解決了編隊(duì)跟蹤問題，但跟蹤軌跡的約束性強(qiáng)，其動(dòng)態(tài)特性表達(dá)式要滿足特定條件，不具有普遍性。在通信拓?fù)鋱D為一般有向圖且滿足LQR性能指標(biāo)的條件下，保證群系統(tǒng)編隊(duì)跟蹤的控制協(xié)議設(shè)計(jì)方法還比較少。

本文主要研究了一種針對(duì)群系統(tǒng)滿足LQR指標(biāo)的編隊(duì)跟蹤優(yōu)化方法。

符號(hào)RN×N表示N×N維矩陣。IN是N×N維單位矩陣。A>B意味著A-B是正定的，A≥B意味著A-B是半正定的。A?B表示矩陣A和矩陣B的克羅內(nèi)克積。

1 圖論知識(shí)和相關(guān)引理

1.1 圖論知識(shí)

圖G=(V,ε,W)中V是節(jié)點(diǎn)集合，ε是邊集合，W是鄰接矩陣。其中，V={ξ1,ξ2,…,ξN}，ε={(ξi,ξj) ∶ξi,ξj∈V}，W=[wij]∈RN×N。用eij=(ξi,ξj)表示節(jié)點(diǎn)ξi延伸到節(jié)點(diǎn)ξj的邊。如果對(duì)任意的eij都存在eji∈ε，則圖G是無向圖。其他情況下，稱圖G是有向圖。用wij表示邊eji的連接權(quán)重。當(dāng)wij>0，節(jié)點(diǎn)ξi可以接收來自節(jié)點(diǎn)ξj的信息，而對(duì)i={1,2,…,N}，有wij=0。如果存在一個(gè)節(jié)點(diǎn)ξi，可以傳遞信息到其他任意節(jié)點(diǎn)，則稱圖G含有一個(gè)有向生成樹。定義節(jié)點(diǎn)ξi的鄰居節(jié)點(diǎn)集合為Ni={ξj∈V∶(ξj,ξi)∈ε}。定義節(jié)點(diǎn)ξi的入度為degin(ξi)，則圖的入度矩陣為D=diag{degin(ξ1),degin(ξ2),…,degin(ξN)}。通常用L=D-W表示圖G的拉普拉斯矩陣。

1.2 相關(guān)引理

引理1[14]有向圖含有向生成樹，則0是拉普拉斯矩陣L的單個(gè)特征值且L1=0，其他非零特征值均具有正實(shí)部。如果無向圖含有向生成樹，其他非零特征值均為正數(shù)。

引理2[15]如果矩陣W∈RN×N每一特征根均有正的實(shí)部，存在正定矩陣Q>0使得WTQ+QW>0。

引理3[16]對(duì)于矩陣IN+L，L是圖G的拉普拉斯矩陣，存在正定矩陣Q，使得

(IN+L)TQ+Q(IN+L)>2αQ

其中，0<α<1。

引理4[17]對(duì)于系統(tǒng)

如果：

1)存在正定函數(shù)V(x,t)；

則平衡狀態(tài)xe=0是漸近穩(wěn)定的。

引理5[18]考慮連續(xù)時(shí)間線性定常系統(tǒng)

選取二次型性能指標(biāo)函數(shù)

其中，A和B為輸入矩陣，(A,B)可控且B是列滿秩矩陣。x∈Rn是第i個(gè)主體的狀態(tài)，u∈Rm是控制輸入。

選取加權(quán)陣Q>0和R>0，控制變量u=Kx，其中，K=-R-1BTP且P>0滿足里卡蒂方程ATP+PA+Q-PBR-1BTP=0。則反饋增益矩陣K使二次型性能指標(biāo)函數(shù)J最小化且矩陣-KB是可對(duì)角化且正定的。

引理6[19]對(duì)于逆最優(yōu)問題，反饋增益矩陣K=-R-1BTP對(duì)于二次型性能指標(biāo)函數(shù)J是最優(yōu)的且P>0滿足里卡蒂方程

ATP+PA+Q-PBR-1BTP=0

(Q>0和R>0)的條件是：

1)-KB是可對(duì)角化且正定的，而且滿足rank(KB)=rank(K)。即存在對(duì)稱矩陣P>0和R>0使得K=-R-1BTP。

2)K穩(wěn)定到矩陣P的核空間。

2 問題描述

2.1 編隊(duì)跟蹤問題

考慮群系統(tǒng)滿足如式(1)所示動(dòng)態(tài)特性。

(1)

式中，ξi(t)∈Rl是第i個(gè)主體的狀態(tài)，ui(t)∈Rm是第i個(gè)主體的控制輸入。

假設(shè)1B是列滿秩矩陣，通信拓?fù)鋱DG包含有向生成樹。

定義1在控制輸入ui(t)下，群系統(tǒng)式(1)的主體狀態(tài)能夠滿足

(2)

式中，r(t)∈Rl為給定的軌跡，則稱群系統(tǒng)式(1)能夠形成時(shí)變編隊(duì)h(t)，同時(shí)可以跟蹤軌跡r(t)。

2.2 控制協(xié)議框架

基于一致性理論，考慮如下編隊(duì)控制協(xié)議：

ui(t)=ui1(t)+ui2(t)+ui3(t),i=1,2,…N

(3)

式中，ui1(t)為自身反饋控制輸入，ui2(t)為輔助函數(shù)輸入，ui3(t)為鄰居反饋控制輸入，其具體表達(dá)式為

(4)

其中，Ni表示拓?fù)鋱DG第i個(gè)節(jié)點(diǎn)的鄰居集合，K1和K2是待設(shè)計(jì)的增益矩陣，vi(t)∈Rm是輔助函數(shù)。

3 問題分析和協(xié)議設(shè)計(jì)

3.1 充分必要條件

由定義1可知，群系統(tǒng)式(1)的主體狀態(tài)與相應(yīng)的時(shí)變編隊(duì)和軌跡的差值需要收斂到零，才能實(shí)現(xiàn)編隊(duì)跟蹤，這樣對(duì)編隊(duì)跟蹤控制協(xié)議的設(shè)計(jì)會(huì)有一定的限制條件，定理1給出了控制協(xié)議框架下群系統(tǒng)實(shí)現(xiàn)編隊(duì)跟蹤的充要條件。

定理1有界初始條件下，群系統(tǒng)式(1)通過ui(t)的控制可以形成時(shí)變編隊(duì)h(t)并跟蹤軌跡r(t)的充分必要條件是：

1)對(duì)于控制函數(shù)vi(t)∈Rm有

(5)

2)如下系統(tǒng)是漸近穩(wěn)定的

(6)

證明：將式(3)代入式(1)，并令

得到群系統(tǒng)的閉環(huán)方程為

(IN?BK1)h(t)+(L?BK2)h(t)

(7)

式中，L為通信拓?fù)鋱DG的拉普拉斯矩陣。令

ζi(t)=ξi(t)-hi(t)-r(t)，i=1,2,…,N

則群系統(tǒng)式(7)可以轉(zhuǎn)換為

(8)

(9)

且式(10)所示閉環(huán)系統(tǒng)是漸近穩(wěn)定的，則群系統(tǒng)式(1)能夠形成時(shí)變編隊(duì)h(t)，同時(shí)可以跟蹤軌跡r(t)。由式(9)可得到定理1的條件1，由式(10)可得到定理1的條件2。根據(jù)上述線性變換的證明過程，條件1和條件2是充分必要條件。

(10)

□

注1：定理1的適用范圍較廣，不同于文獻(xiàn)[20]，除了考慮群系統(tǒng)實(shí)現(xiàn)時(shí)變編隊(duì)，還加以考慮跟蹤軌跡。其控制函數(shù)v(t)的引入是為了補(bǔ)償時(shí)變編隊(duì)h(t)和軌跡r(t)帶來的多余項(xiàng)，將群系統(tǒng)式(1)轉(zhuǎn)化為閉環(huán)自治系統(tǒng)，便于討論系統(tǒng)的穩(wěn)定性。其控制函數(shù)v(t)可以通過式(5)求解得到，條件1較易滿足。而條件2可通過李雅普諾夫穩(wěn)定性理論予以證明，在理論證明中，難點(diǎn)在于閉環(huán)系統(tǒng)式(6)中增益矩陣K1和K2的設(shè)計(jì)方法，下面定理2給出了增益矩陣的設(shè)計(jì)方法。

3.2 穩(wěn)定性分析

定理2有界初始條件下，如果K=-R-1BTP，若K1=cK，K2=-cK，c是增益常數(shù)，群系統(tǒng)式(1)在編隊(duì)跟蹤控制協(xié)議式(4)下能夠?qū)崿F(xiàn)時(shí)變編隊(duì)并跟蹤軌跡。

證明：考慮分段連續(xù)的李雅普諾夫函數(shù)

V=ζT(t)(Qσ?P)ζ(t)

(11)

P是里卡蒂方程ATP+PA+Q-PBR-1BTP=0的正定解。V是連續(xù)的，對(duì)其求導(dǎo)并將式(6)、K1=cK、K2=-cK、K=-R-1BTP代入式(11)可得

cζT(t){[(IN+L)TQσ]?PBR-1BTP}ζ(t)-

cζT(t){[Qσ(IN+L)]?PBR-1BTP}ζ(t)

(12)

由引理3得

(IN+L)TQσ+Qσ(IN+Lσ(t))>2αQσ

(13)

則式(12)轉(zhuǎn)換為

(14)

由里卡蒂方程可得

ATP+PA+Q-PBR-1BTP=0

(15)

則式(14)轉(zhuǎn)換為

(16)

(17)

□

注2：由于跟蹤控制協(xié)議式(4)包含了自身反饋控制輸入，從式(12)～(16)判斷，假設(shè)1并不是群系統(tǒng)實(shí)現(xiàn)時(shí)變編隊(duì)和跟蹤軌跡的必要條件，即使各主體之間不連通，時(shí)變編隊(duì)和軌跡跟蹤仍然可以實(shí)現(xiàn)。根據(jù)文獻(xiàn)[16]可知，通信拓?fù)鋱D中含有有向生成樹，能夠增強(qiáng)群系統(tǒng)的魯棒性。

在實(shí)際應(yīng)用中，往往還需要考慮系統(tǒng)能否滿足某種性能指標(biāo)，下面定理3給出了編隊(duì)跟蹤控制協(xié)議式(4)可以最小化LQR性能指標(biāo)的充分條件。

3.3 最優(yōu)性分析

考慮如下LQR性能指標(biāo)：

(18)

定理3有界初始條件下，選取編隊(duì)跟蹤控制協(xié)議式(4)可以最小化LQR性能指標(biāo)的充分條件是：矩陣IN+L是可對(duì)角化且正定的。

證明：令K1=cK，K2=-cK，K=-R-1BTP，代入式(10)，得到全局誤差動(dòng)態(tài)特性

=IN?Aζ(t)+c(IN+L)?BKζ(t)

=IN?Aζ(t)+(IN?B)c(IN+L)?Kζ(t)

=IN?Aζ(t)+(IN?B)v

(19)

由引理5可知，矩陣-KB是可對(duì)角化且正定的，而且滿足rank(KB)=rank(K)，則有

rank{[c(IN+L)?K](IN?B)}

=rank[c(IN+L)]rank(KB)

=rank[c(IN+L)]rank(K)

=rank[c(IN+L)?K]

(20)

□

=-(IN?R)-1(IN?B)T(Lσ?P)

(21)

(22)

(23)

注4：在滿足條件式(5)和閉環(huán)系統(tǒng)式(6)漸近穩(wěn)定的情況下，編隊(duì)跟蹤控制協(xié)議式(4)可以使群系統(tǒng)形成指定時(shí)變編隊(duì)并跟蹤軌跡，但不一定可以最小化二次型性能指標(biāo)，通信拓?fù)溥€需滿足特定條件，即IN+L是可對(duì)角化且正定的，才可以實(shí)現(xiàn)最小化二次型性能指標(biāo)。根據(jù)引理1，當(dāng)拓?fù)鋱D為無向圖時(shí)，矩陣IN+L是實(shí)對(duì)稱矩陣，必定可以對(duì)角化，且特征值都為正的實(shí)數(shù)，滿足正定性。

注5：傳統(tǒng)最小化性能指標(biāo)的方法，是先選定加權(quán)陣，確定性能指標(biāo)之后，得出控制輸入的最優(yōu)形式，進(jìn)而決定控制協(xié)議的設(shè)計(jì)結(jié)構(gòu)。而傳統(tǒng)方法大都是針對(duì)低階積分特性模型，如文獻(xiàn)[11]，或者不考慮時(shí)變編隊(duì)和跟蹤軌跡，如文獻(xiàn)[12]。對(duì)于高階群系統(tǒng)模型，考慮時(shí)變編隊(duì)和跟蹤軌跡的影響，選定性能指標(biāo)之后，設(shè)計(jì)控制協(xié)議是十分困難的，基于此，本文設(shè)計(jì)的控制協(xié)議使得群系統(tǒng)轉(zhuǎn)化為自治系統(tǒng)，通過變量代換，得到性能指標(biāo)的加權(quán)陣，能夠最小化性能指標(biāo)。

4 數(shù)值仿真與分析

4.1 數(shù)值仿真條件

圖1給出了各主體之間的拓?fù)鋱D。可以看出，圖1中的拓?fù)鋱DG包含有向生成樹，G相應(yīng)的拉普拉斯矩陣為

圖1 群系統(tǒng)各主體之間的通信拓?fù)鋱DFig.1 Communication topologies between agents of swarm systems

基于多無人機(jī)偵察任務(wù)場(chǎng)景，無人機(jī)須從不同角度偵察目標(biāo)，以及被偵察目標(biāo)軌跡會(huì)發(fā)生變化，定義時(shí)變編隊(duì)h(t)和軌跡r(t)分別為

其中，ω=0.25 rad/s。由式(5)可以得到控制函數(shù)v(t)。

當(dāng)0≤t<30時(shí)，

當(dāng)30≤t≤70時(shí)，

取α=0.5，得c=1。各主體的初始狀態(tài)分別為：xi1(0)=4(δ-0.5)，xi2(0)=3(δ-0.5)，xi3(0)=2(δ-0.5)(i=1,2,…,8)。其中，δ為(0,1)之間的隨機(jī)數(shù)，設(shè)置仿真時(shí)間70 s。

4.2 結(jié)果討論與對(duì)比分析

4.2.1 結(jié)果討論

圖2給出了8個(gè)主體的狀態(tài)在0.5 s、10 s、30 s、45 s、60 s和70 s時(shí)的狀態(tài)演化過程和編隊(duì)構(gòu)型。在初始階段，8個(gè)主體的構(gòu)型為不規(guī)則圖形，隨著時(shí)間的推移，群系統(tǒng)的8個(gè)主體形成了指定的時(shí)變編隊(duì)并可以保持穩(wěn)定。

(a) t=0.5 s

(b) t=10 s

(c) t=30 s

(d) t=45 s

(e) t=60 s

(f) t=70 s圖2 不同時(shí)刻群系統(tǒng)狀態(tài)演化過程Fig.2 State evolution process of the swarm systems at different times

圖3為20～30 s和50～60 s之間各主體的軌跡圖，其中方框表示主體的起點(diǎn)，五角星表示主體的終點(diǎn)。在20～30 s和50～60 s的時(shí)間段里，8個(gè)主體已經(jīng)形成了規(guī)則的構(gòu)型，并可以跟蹤預(yù)先設(shè)定的軌跡。當(dāng)軌跡發(fā)生變化時(shí)，群系統(tǒng)可以保持對(duì)軌跡的跟蹤且編隊(duì)的構(gòu)型并沒有受到影響。

(a) 20～30 s

(b) 50～60 s圖3 不同時(shí)間段的軌跡截圖Fig.3 Snapshot of motion trajectories in different time

圖4(a)～(c)分別給出了編隊(duì)跟蹤誤差在3個(gè)不同方向上的差值曲線。從圖中不難看出，不同方向上編隊(duì)跟蹤誤差都可以趨于零。說明各主體的三個(gè)狀態(tài)與編隊(duì)和軌跡相應(yīng)狀態(tài)的差值趨于零，這也說明群系統(tǒng)形成了指定的時(shí)變編隊(duì)，并可以保持穩(wěn)定。

(a) ξi1(t)-hi1(t)-r(t)

(b) ξi2(t)-hi2(t)-r(t)

(c) ξi3(t)-hi3(t)-r(t)圖4 3個(gè)不同方向上的編隊(duì)跟蹤誤差Fig.4 Formation tracking errors in three different directions

從增益常數(shù)c分析，隨著增益常數(shù)c取值的不同，編隊(duì)跟蹤誤差趨于零的時(shí)間長(zhǎng)短不同，即群系統(tǒng)形成時(shí)變編隊(duì)并跟蹤軌跡的時(shí)間長(zhǎng)短不同。以x方向?yàn)槔?，圖5(a)～(c)分別給出了增益常數(shù)c在取值為0.6、1和2時(shí)的編隊(duì)跟蹤誤差曲線。這里認(rèn)為差值小于10-3時(shí)，編隊(duì)跟蹤誤差趨于零。圖中豎直虛線表示差值趨于零的時(shí)間?？梢钥闯鲭S著增益常數(shù)c的增大，編隊(duì)跟蹤誤差趨于零的時(shí)間越短。圖6(a)～(b)分別給出了增益常數(shù)c在取值為0.6和2時(shí)，8個(gè)主體的狀態(tài)在0.5 s時(shí)的狀態(tài)截圖，與圖2(a)中c取值為1對(duì)比，可以看出隨著增益常數(shù)c的增大，8個(gè)主體的狀態(tài)在0.5 s時(shí)的編隊(duì)構(gòu)型越規(guī)則，也說明8個(gè)主體的編隊(duì)形成速度越快。

(a) c=0.6

(b) c=1

(c) c=2圖5 ξi1(t)-hi1(t)-r(t)在不同取值c的誤差曲線Fig.5 The error curves of ξi1(t)-hi1(t)-r(t) in different values of c

(a) t=0.5 s，c=0.6

(b) t=0.5 s，c=2圖6 群系統(tǒng)狀態(tài)在c不同取值時(shí)的截圖Fig.6 Snapshot of states of the swarm systems in different values of c

注6：基于多無人機(jī)偵察任務(wù)場(chǎng)景，可以發(fā)現(xiàn)，當(dāng)其中一個(gè)無人機(jī)可以感知到危險(xiǎn)源時(shí)，需要聯(lián)動(dòng)整個(gè)無人機(jī)群系統(tǒng)突然改變方向。很明顯這種情況下的無人機(jī)群系統(tǒng)不僅需要位置和速度的誤差收斂于零，還要實(shí)現(xiàn)加速度的誤差收斂于零，才可以實(shí)現(xiàn)更加精準(zhǔn)的飛行。因此，針對(duì)高階系統(tǒng)模型，設(shè)此仿真案例，以驗(yàn)證編隊(duì)控制算法。

4.2.2 對(duì)比分析

為便于比較和考慮到控制協(xié)議最小化性能指標(biāo)需要滿足特定的拓?fù)錀l件，選取拓?fù)鋱D為無向圖，如圖7所示。

圖7 無向圖Fig.7 Undirected graph

其他仿真條件采用本文的仿真條件。

從不同控制方法分析，編隊(duì)跟蹤誤差趨于零的時(shí)間長(zhǎng)短不同。以x方向?yàn)槔瑘D8給出了文獻(xiàn)[16]方法和本文方法下x方向的收斂時(shí)間。圖中豎直虛線是增益常數(shù)c選取的邊界，即參數(shù)c的選取要大于0.5。利用文獻(xiàn)[16]方法得到的增益矩陣，使得編隊(duì)跟蹤誤差趨于零的時(shí)間為6.3 s，而利用本文方法得到的增益矩陣，可以通過調(diào)節(jié)增益常數(shù)c，改變編隊(duì)跟蹤誤差趨于零的快慢。當(dāng)c取0.5時(shí)，收斂時(shí)間為6.3 s；當(dāng)c取5時(shí)，收斂時(shí)間為2.3 s。從圖中可以看出，隨著增益常數(shù)c的增大，收斂時(shí)間從6.3 s減小至2.3 s，而相應(yīng)的最小性能指標(biāo)J*越來越大，在實(shí)際應(yīng)用中，可以綜合考慮收斂時(shí)間與性能指標(biāo)對(duì)系統(tǒng)的影響，適當(dāng)選取增益常數(shù)c。在實(shí)際應(yīng)用中，基于多無人機(jī)偵察任務(wù)場(chǎng)景，當(dāng)一個(gè)無人機(jī)收到被偵察目標(biāo)的軌跡信息，需要無人機(jī)群系統(tǒng)迅速做出反應(yīng)，為達(dá)到良好的偵察效果，無人機(jī)群系統(tǒng)需要迅速形成編隊(duì)隊(duì)形實(shí)施偵察任務(wù)，映射到仿真案例中，調(diào)節(jié)增益常數(shù)c的大小以改變編隊(duì)的形成速度，從而適應(yīng)偵察任務(wù)的需要。

圖8 收斂時(shí)間和性能指標(biāo)與參數(shù)c的關(guān)系Fig.8 Relationship between convergence time performance index and the parameter c

5 結(jié)論

本文提出了一種滿足LQR性能指標(biāo)的群系統(tǒng)編隊(duì)跟蹤優(yōu)化方法，結(jié)論為：

1)建立了編隊(duì)跟蹤問題數(shù)學(xué)描述，設(shè)計(jì)分布式控制協(xié)議框架，給出了編隊(duì)跟蹤的充要條件。

2)通過變量代換，將群系統(tǒng)編隊(duì)跟蹤控制問題轉(zhuǎn)化為閉環(huán)系統(tǒng)的穩(wěn)定性問題，并利用李雅普諾夫第二方法，分析閉環(huán)系統(tǒng)的穩(wěn)定性。

3)給出了滿足LQR性能指標(biāo)下通信拓?fù)錀l件。