【Paper】2007_Consensus control for a class of networks of dynamic agents 二阶静态一致性
文章目錄
- 3. Control Protocol and Network Dynamics
- 4. Network with Fixed Topology
- 定理 1
- Remark 1
- Lemma 2
- 證明 Lemma 2
- 證明 定理 1
3. Control Protocol and Network Dynamics
在本節中,我們將介紹解決上述平均一致性問題的控制協議。我們將使用固定拓撲和無通信時延的線性協議:
ui=ui1+ui2(5)u_i = u_{i1} + u_{i2} \tag{5}ui?=ui1?+ui2?(5)
其中
ui1=kvi(6)u_{i1} = k v_i \tag{6}ui1?=kvi?(6)
是有著反饋增益 kkk 的本地速度反饋,反饋增益之后會來設計,并且
ui2=∑j∈Niaij(xj?xi)(7)u_{i2} = \sum_{j\in N_i} a_{ij} (x_j - x_i) \tag{7}ui2?=j∈Ni?∑?aij?(xj??xi?)(7)
是節點 viv_ivi? 的鄰居對應的部分,在拓撲固定的網絡中是恒定的,在拓撲切換的網絡中是可變的。
通過使用協議(5),智能體動力學可寫成如下形式:
x˙i=vimv˙i=kvi+∑j∈Niaij(xj?xi)(8)\begin{aligned} \dot{x}_i & = v_i \\ m \dot{v}_i & = k v_i + \sum_{j\in N_i} a_{ij} (x_j - x_i) \end{aligned}\tag{8}x˙i?mv˙i??=vi?=kvi?+j∈Ni?∑?aij?(xj??xi?)?(8)
定義
ξi=[xi,vi]T\xi_i = [x_i, v_i]^\text{T}ξi?=[xi?,vi?]T
我們有
ξ˙i=Aξi+BKξi+BF∑j∈Niaij(ξj?ξi)(9)\dot{\xi}_i = A \xi_i + B K \xi_i + B F \sum_{j\in N_i} a_{ij} (\xi_j - \xi_i) \tag{9}ξ˙?i?=Aξi?+BKξi?+BFj∈Ni?∑?aij?(ξj??ξi?)(9)
其中 A=[0100]A = \left[\begin{matrix} 0 & 1 \\ 0 & 0 \\ \end{matrix}\right]A=[00?10?],B=[01]B = \left[\begin{matrix} 0 \\ 1 \\ \end{matrix}\right]B=[01?],K=[0k]K = \left[\begin{matrix} 0 & k \end{matrix}\right]K=[0?k?],F=[10]F = \left[\begin{matrix} 1 & 0 \end{matrix}\right]F=[1?0?]。
進一步,定義
ξ=[ξ1ξ2?ξM]\xi = \left[\begin{matrix} \xi_1 \\ \xi_2 \\ \vdots \\ \xi_M \\ \end{matrix}\right]ξ=??????ξ1?ξ2??ξM????????
動力學模型描述為
ΦG(t)=IM?(A+BK)?LG(t)?BF(10)\Phi_{G(t)} = I_M \otimes (A+BK) - L_{G(t)} \otimes B F \tag{10}ΦG(t)?=IM??(A+BK)?LG(t)??BF(10)
其中 LG(t)L_{G(t)}LG(t)? 是之前提到的關于 ttt 時刻圖 G(t)G(t)G(t) 的拉普拉斯矩陣。
如果針對任意時間 ttt 都有圖 G(t)≡GG(t) \equiv GG(t)≡G,網絡的動力學是一個線性時不變系統。否則的話,網絡的動力學通常混合了連續狀態 ξ∈R2M\xi \in \R^{2M}ξ∈R2M 和離散狀態 GGG。結構就是,網絡成為一個典型的線性時間切換系統。
4. Network with Fixed Topology
本節,我們提供了具有固定拓撲網絡的平均一致性的分析。
定理 1
考慮一個固定拓撲網絡 GGG 是連通圖,假設所有初始速度都是 000,即 vi(0)=0v_i(0) = 0vi?(0)=0,那么針對任意的負增益 k<0k<0k<0,協議(5)全局漸進解決了平均一致性問題。
Remark 1
一個關于定理 1 更加通用的形式已經在文獻 Tanner HG, Jadbabaie A, Pappas GJ. Stable ?ocking of mobile agents, part I: ?xed topology. Proceedings of the IEEE Conference on Decision and Control, Hyatt Regency Maui, HI, U.S.A., vol. 2. 2003; 2010–2015. 中用一個不同的方式進行了證明。
在給出定理 1 的證明前,我們首先考慮 lim?t→∞exp?(Φt)\lim_{t\rightarrow \infty} \exp(\Phi t)limt→∞?exp(Φt) 的解。
Lemma 2
假設圖 GGG 是連通的,那么對于任意 k<0k<0k<0,
lim?t→∞exp?(Φt)=wrwlT(11)\lim_{t\rightarrow \infty} \exp(\Phi t) = w_r w_l^\text{T} \tag{11}t→∞lim?exp(Φt)=wr?wlT?(11)
其中 wl,wrw_l, w_rwl?,wr? 分別是 Φ\PhiΦ 關于 000 特征值的左特征向量和右特征向量。進一步地,
wr=1M1M?[10]T(12)w_r = \frac{1}{\sqrt{M}} \textbf{1}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T} \tag{12}wr?=M?1?1M??[1?0?]T(12)
wl=1?kM1M?[?k1]T(13)w_l = \frac{1}{-k \sqrt{M}} \textbf{1}_M \otimes [\begin{matrix} -k & 1 \end{matrix}]^\text{T} \tag{13}wl?=?kM?1?1M??[?k?1?]T(13)
并且 wrTwl=1w_r^\text{T} w_l = 1wrT?wl?=1。
證明 Lemma 2
首先,我們有
Φ1M1M?[10]T=1M(IM?(A+BK)?L?BF)1M?[10]T=1M(IM1M)?((A+BK)[10]T)?1M(L1M)?(BF[10]T)=1M1M?02?1M0M?[10]T=02M\begin{aligned} \Phi \frac{1}{\sqrt{M}} \textbf{1}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T} &= \frac{1}{\sqrt{M}} (I_M \otimes (A + B K) - L \otimes B F) \textbf{1}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T} \\ &= \frac{1}{\sqrt{M}} (I_M \textbf{1}_M) \otimes ((A+BK) [\begin{matrix} 1 & 0 \end{matrix}]^\text{T}) -\frac{1}{\sqrt{M}} (L \textbf{1}_M) \otimes (BF [\begin{matrix} 1 & 0 \end{matrix}]^\text{T}) \\ &= \frac{1}{\sqrt{M}} \textbf{1}_M \otimes \textbf{0}_2 - \frac{1}{\sqrt{M}} \textbf{0}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T} \\ &= \textbf{0}_{2M} \end{aligned}ΦM?1?1M??[1?0?]T?=M?1?(IM??(A+BK)?L?BF)1M??[1?0?]T=M?1?(IM?1M?)?((A+BK)[1?0?]T)?M?1?(L1M?)?(BF[1?0?]T)=M?1?1M??02??M?1?0M??[1?0?]T=02M??
這意味著 wr=1/M1M?[10]Tw_r = 1 / \sqrt{M} \textbf{1}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T}wr?=1/M?1M??[1?0?]T。
接下來
現在我們給出定理 1 的證明。
證明 定理 1
根據引理 2,我們有
ξ(t)=exp?(Φt)ξ(0)\xi(t) = \exp (\Phi t) \xi(0)ξ(t)=exp(Φt)ξ(0)
由此可見
lim?t→∞ξ(t)=lim?t→∞exp?(Φt)ξ(0)=wrwlTξ(0)=(wlTξ(0))wr=1M(wlTξ(0))1M?[10]T\begin{aligned} \lim_{t\rightarrow \infty} \xi(t) &= \lim_{t\rightarrow \infty} \exp(\Phi t) \xi(0) \\ &= w_r w_l^\text{T} \xi(0) \\ &= (w_l^\text{T} \xi(0)) w_r \\ &= \frac{1}{\sqrt{M}} (w_l^\text{T} \xi(0)) \textbf{1}_M \otimes [\begin{matrix} 1 & 0 \end{matrix}]^\text{T} \end{aligned}t→∞lim?ξ(t)?=t→∞lim?exp(Φt)ξ(0)=wr?wlT?ξ(0)=(wlT?ξ(0))wr?=M?1?(wlT?ξ(0))1M??[1?0?]T?
由于
1M(wlTξ(0))=1M(1MT?[1?1k])[x1(0)v1(0)?xM(0)vM(0)]=1M∑i=1M(xi(0)?vi(0)k)=Ave(x0)\begin{aligned} \frac{1}{\sqrt{M}} (w_l^\text{T} \xi(0)) &= \frac{1}{M} (\textbf{1}_M^\text{T} \otimes [\begin{matrix} 1 & -\frac{1}{k} \end{matrix}]) [\begin{matrix} x_1(0) & v_1(0) & \cdots & x_M(0) & v_M(0) \end{matrix}] \\ &= \frac{1}{M} \sum_{i=1}^M (x_i(0) - \frac{v_i(0)}{k}) \\ &= \text{Ave}(x_0) \end{aligned}M?1?(wlT?ξ(0))?=M1?(1MT??[1??k1??])[x1?(0)?v1?(0)???xM?(0)?vM?(0)?]=M1?i=1∑M?(xi?(0)?kvi?(0)?)=Ave(x0?)?
并且能夠知道的是
lim?t→∞vi(t)=0,i=1,2,?,M\lim_{t\rightarrow \infty} v_i(t) = 0, \quad i=1,2,\cdots,Mt→∞lim?vi?(t)=0,i=1,2,?,M
這表明協議(5)全局漸進解決了平均一致性問題。完成了證明。
總結
以上是生活随笔為你收集整理的【Paper】2007_Consensus control for a class of networks of dynamic agents 二阶静态一致性的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 【Paper】2013_Cooperat
- 下一篇: 【Paper】2021_Analysis