• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于傳感器攻擊的領(lǐng)導(dǎo)跟隨多旋翼蜂群驅(qū)離方法

    2024-02-22 00:00:00許兆勝王樂劉柟池趙良玉席建祥林浩申
    航空兵器 2024年6期

    摘 要:""""" 從攻擊者角度研究領(lǐng)導(dǎo)跟隨多旋翼蜂群的驅(qū)離反制問題, 結(jié)合領(lǐng)導(dǎo)者無(wú)人機(jī)根節(jié)點(diǎn)屬性, 提出一種針對(duì)領(lǐng)導(dǎo)者無(wú)人機(jī)的傳感器攻擊方法, 利用蜂群分布式網(wǎng)絡(luò)交互特點(diǎn), 有序穩(wěn)定地驅(qū)離蜂群至攻擊者期望的路徑。 首先結(jié)合無(wú)人機(jī)動(dòng)態(tài)特性構(gòu)建傳感器攻擊模型, 將攻擊信號(hào)引入蜂群編隊(duì)跟蹤控制過程; 其次考慮蜂群通信拓?fù)淝袚Q, 采用狀態(tài)空間分解方法將攻擊信號(hào)分解到編隊(duì)跟蹤參考子系統(tǒng)和編隊(duì)跟蹤誤差子系統(tǒng), 給出了傳感器攻擊下實(shí)現(xiàn)蜂群驅(qū)離反制的充分條件, 并通過兩步非奇異變換簡(jiǎn)化分析判據(jù); 最后通過數(shù)值仿真驗(yàn)證本文所提驅(qū)離條件的合理性。

    關(guān)鍵詞:"""" 傳感器攻擊; 領(lǐng)導(dǎo)跟隨; 蜂群反制; 多旋翼蜂群; 分布式網(wǎng)絡(luò)

    中圖分類號(hào):""""" V249.12

    文獻(xiàn)標(biāo)識(shí)碼:""" A

    文章編號(hào):"""" 1673-5048(2024)06-0104-08

    DOI: 10.12132/ISSN.1673-5048.2024.0117

    0 引" 言

    分布式多旋翼無(wú)人機(jī)蜂群(以下簡(jiǎn)稱無(wú)人機(jī)蜂群)因其優(yōu)秀的靈活性、 機(jī)動(dòng)性和可重構(gòu)性, 在災(zāi)害搜救、 低空運(yùn)輸、 遙感測(cè)繪、 戰(zhàn)場(chǎng)感知等經(jīng)濟(jì)民生領(lǐng)域和國(guó)防軍事領(lǐng)域得到廣泛應(yīng)用。 但同時(shí), 無(wú)人機(jī)蜂群技術(shù)的濫用帶來(lái)的安全威脅也日益突出[1, 如何反制驅(qū)離無(wú)人機(jī)蜂群成為當(dāng)下研究熱點(diǎn)。 為了對(duì)抗無(wú)人機(jī)蜂群帶來(lái)的威脅, 文獻(xiàn)[2-5]從策略手段和裝備功能等方面對(duì)無(wú)人機(jī)蜂群反制方法進(jìn)行了綜述, 鮮有詳細(xì)的反制技術(shù)分析。 文獻(xiàn)[6-8]通過欺騙無(wú)人機(jī)機(jī)載導(dǎo)航定位傳感器實(shí)現(xiàn)了對(duì)單個(gè)無(wú)人機(jī)的反制, 該技術(shù)理論上可以通過逐個(gè)擊破的方式實(shí)現(xiàn)對(duì)無(wú)人機(jī)蜂群整體的反制, 但是存在費(fèi)效比高以及效率低的缺陷。 文獻(xiàn)[9]通過重建目標(biāo)保護(hù)區(qū)內(nèi)導(dǎo)航定位信號(hào)環(huán)境, 直接實(shí)現(xiàn)了對(duì)整個(gè)無(wú)人機(jī)蜂群的空間欺騙, 但該方法同時(shí)也導(dǎo)致保護(hù)區(qū)內(nèi)其他依賴導(dǎo)航功能的設(shè)備失效。

    以上的反制技術(shù)在對(duì)無(wú)人機(jī)蜂群實(shí)施反制時(shí)均忽略了無(wú)人機(jī)蜂群中個(gè)體與個(gè)體之間的網(wǎng)絡(luò)交互特征。 無(wú)人機(jī)蜂群是采用分布式架構(gòu)的自主或半自主個(gè)體通過網(wǎng)絡(luò)交互構(gòu)成的群體智能系統(tǒng)。 文獻(xiàn)[10-12]利用分布式網(wǎng)絡(luò)交互的特點(diǎn)對(duì)群體智能系統(tǒng)實(shí)施攻擊, 分析了攻擊信號(hào)在網(wǎng)絡(luò)中的傳播影響。 其中, 文獻(xiàn)[10]針對(duì)分布式機(jī)器人群編隊(duì)控制問題, 提出了一種基于位置偏差的劫持攻擊策略, 驅(qū)使機(jī)器人編隊(duì)偏離原運(yùn)動(dòng)方向。 文獻(xiàn)[11]通過執(zhí)行器攻擊, 實(shí)現(xiàn)了無(wú)向拓?fù)錀l件下多智能體系統(tǒng)的編隊(duì)運(yùn)動(dòng)驅(qū)離。 文獻(xiàn)[12]分別研究了虛假偏差和給定狀態(tài)兩種欺騙信號(hào)對(duì)多領(lǐng)導(dǎo)者無(wú)人機(jī)蜂群編隊(duì)的攻擊問題, 實(shí)現(xiàn)了對(duì)無(wú)人機(jī)蜂群的軌跡欺騙。 可以發(fā)現(xiàn), 文獻(xiàn)[10-12]僅考慮了固定拓?fù)錀l件下的無(wú)人機(jī)蜂群驅(qū)離反制問題, 實(shí)際應(yīng)用中, 蜂群的通信網(wǎng)絡(luò)可能因干擾或者受通信距離影響出現(xiàn)切換的情況, 同時(shí)領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群執(zhí)行任務(wù)時(shí)往往需要保持特定的編隊(duì)隊(duì)形, 并且跟蹤參考航跡或任務(wù)目標(biāo)路徑, 目前鮮有文獻(xiàn)解決類似無(wú)人機(jī)蜂群的驅(qū)離反制問題。

    受文獻(xiàn)[10-12]等研究成果啟發(fā), 本文以領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群為反制對(duì)象, 意圖通過傳感器攻擊, 在不破壞無(wú)人機(jī)蜂群編隊(duì)隊(duì)形的情況下, 驅(qū)離蜂群偏離原有任務(wù)目標(biāo)路徑, 誘騙其按照攻擊者期望軌跡運(yùn)動(dòng), 避免次生災(zāi)害, 實(shí)現(xiàn)對(duì)無(wú)人機(jī)蜂群的安全迫降或者捕獲。

    本文在切換拓?fù)錀l件下, 從攻擊者角度研究針對(duì)領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群編隊(duì)的誘騙驅(qū)離反制方法。 與現(xiàn)有驅(qū)離反制方法相比, 本文主要貢獻(xiàn)點(diǎn)有三個(gè): 一是結(jié)合無(wú)人機(jī)動(dòng)力學(xué)特性, 構(gòu)建針對(duì)領(lǐng)導(dǎo)者無(wú)人機(jī)的傳感器攻擊模型, 分析傳感器攻擊對(duì)切換拓?fù)湎骂I(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群編隊(duì)的影響, 明確了切換拓?fù)涞挠绊憴C(jī)理; 二是給出了無(wú)人機(jī)蜂群在傳感器攻擊下能夠被誘騙驅(qū)離的分析判據(jù), 且被誘騙驅(qū)離的無(wú)人機(jī)蜂群編隊(duì)隊(duì)形沒有發(fā)生改變; 三是利用共同Lyapunov候選函數(shù)和線性矩陣不等式方法, 確保了攻擊后的無(wú)人機(jī)蜂群編隊(duì)收斂性, 通過將所有切換拓?fù)浼现蠰aplacian矩陣特征值壓縮到一個(gè)凸集, 確保所得線性矩陣不等式條件數(shù)量與無(wú)人機(jī)數(shù)量無(wú)關(guān), 降低了計(jì)算復(fù)雜度, 解決了文獻(xiàn)[12]判據(jù)中可能隨無(wú)人機(jī)數(shù)量增加帶來(lái)計(jì)算爆炸問題。

    在本文中, 代表實(shí)矩陣集合, 0代表維數(shù)適當(dāng)?shù)牧阆蛄浚?1N表示元素均為1的N維列向量, IN表示N維單位矩陣, 符號(hào)表示Kronecker積運(yùn)算。

    1 問題描述

    1.1 切換通信拓?fù)?/p>

    設(shè)G=(V, ε, W)表示具有N個(gè)節(jié)點(diǎn)的拓?fù)鋱D, 圖G由節(jié)點(diǎn)集合V={v1, v2, …vN}、 邊集合ε{(vi, vj), vi, vj∈V}和鄰接矩陣

    W=[wij]∈

    組成。 eij=(vi, vj)用于表示連接節(jié)點(diǎn)vi與節(jié)點(diǎn)vj的邊, wij表示邊eji的權(quán)值。 圖G的Laplacian矩陣表示為=[lij]∈瘙綆N×N, 其中l(wèi)ij=-wij, i≠j且lii=∑Nk=1, k≠iwik, i=1, 2, …, N。 如果圖G中有一個(gè)節(jié)點(diǎn)到其余所有節(jié)點(diǎn)間至少存在一條有向路徑, 稱該圖包含生成樹。 其他圖論知識(shí)參見文獻(xiàn)[13]。

    無(wú)人機(jī)蜂群個(gè)體之間的通信拓?fù)潢P(guān)系隨時(shí)間切換, 其所有可能的拓?fù)浼媳硎緸镚~={G1, G2, …, Gs}(s∈, 其中為自然數(shù)集合, 上標(biāo)sgt;1表示對(duì)應(yīng)的拓?fù)渚幪?hào)。 切換信號(hào)表示為σ(t):[0, ∞)→{1, 2, …, s}, 切換時(shí)刻為0=t0lt;t1lt;t2lt;…lt;tm=t, 假設(shè)信號(hào)σ(t)不能無(wú)限快速切換, 切換順序任意, 在時(shí)間區(qū)間[tk, tk+1)內(nèi)拓?fù)潢P(guān)系Gσ(t)是固定保持不變的, 其對(duì)應(yīng)的鄰接矩陣表示為Wσ(t), 其對(duì)應(yīng)的Laplacian矩陣表示為σ(t), t~k=tk+1-tk表示該時(shí)間區(qū)間內(nèi)的駐留時(shí)間。

    假設(shè)1: 無(wú)人機(jī)蜂群中領(lǐng)導(dǎo)者無(wú)人機(jī)不接收跟隨者無(wú)人機(jī)信息, 僅單向發(fā)送協(xié)同信息至跟隨者無(wú)人機(jī); 跟隨者無(wú)人機(jī)之間的信息交互是無(wú)向連通的。

    引理1[14: 如果圖Gσ(t)包含有向生成樹, 那么零是其對(duì)應(yīng)Laplacian矩陣σ(t)的單一特征值, 1N是其相應(yīng)的一個(gè)特征向量, 即σ(t)1N=0, 其余非零特征值均具有正實(shí)部。

    1.2 多旋翼無(wú)人機(jī)蜂群編隊(duì)跟蹤描述

    考慮由N個(gè)多旋翼無(wú)人機(jī)組成的領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群, 其中領(lǐng)導(dǎo)者無(wú)人機(jī)(下標(biāo)記為N)需要跟蹤期望的任務(wù)目標(biāo)路徑Td(t)∈瘙綆6×1, 其余跟隨者無(wú)人機(jī)(下標(biāo)集合記為F={1, 2, …, N-1})按照一定的時(shí)變編隊(duì)隊(duì)形跟蹤領(lǐng)導(dǎo)者無(wú)人機(jī)。 多旋翼無(wú)人機(jī)利用電機(jī)提供的升力組合實(shí)現(xiàn)空間中的外環(huán)位速控制和俯仰、 偏航、 橫滾的內(nèi)環(huán)姿態(tài)控制, 其中外環(huán)動(dòng)力學(xué)模型描述為

    p·i(t)=vi(t)

    v·i(t)=geZ-Fi(t)MiRi(t)eZ

    (1)

    式中: i=1, 2, …, N; pi(t)=[piX(t), piY(t), piZ(t)]T表示地理坐標(biāo)系下的無(wú)人機(jī)i的X, Y, Z三軸上的位置狀態(tài)向量; vi(t)=[viX(t), viY(t), viZ(t)]T表示地理坐標(biāo)系下無(wú)人機(jī)i的速度狀態(tài)向量; g是重力加速度常數(shù); eZ=[0, 0, 1]T表示Z軸正方向的單位向量; Mi表示第i個(gè)無(wú)人機(jī)的質(zhì)量; Fi(t)表示第i個(gè)無(wú)人機(jī)的總升力; Ri(t)表示機(jī)體坐標(biāo)系到地理坐標(biāo)系的旋轉(zhuǎn)矩陣:

    Ri(t)= cosθi(t)cosψi(t)sinθi(t)cosψi(t)sinφi(t)-sinψi(t)cosφi(t)sinθi(t)cosψi(t)cosφi(t)+sinψi(t)sinφi(t)

    cosθi(t)sinψi(t)sinθi(t)sinψi(t)sinφi(t)+cosψi(t)cosφi(t)sinθi(t)sinψi(t)cosφi(t)-cosψi(t)sinφi(t)

    -sinθi(t)cosθi(t)sinφi(t)cosθi(t)cosφi(t)

    (2)

    式中: θi(t)為地理坐標(biāo)系下的橫滾角; ψi(t)為偏航角; φi(t)為俯仰角。 定義-π/2lt;θi(t)lt;π/2, -π/2lt;φi(t)lt;π/2。 令ui(t)=[uiX(t), uiY(t), uiZ(t)]T=[v·iX(t), v·iY(t), v·iZ(t)]T表示外環(huán)位速控制輸入, 結(jié)合式(1)~(2)可得

    uiX(t)=-Fi(t)Mi(sinθi(t)cosψi(t)cosφi(t)+

    sinψi(t)sinφi(t))

    uiY(t)=-Fi(t)Mi(sinθi(t)sinψi(t)cosφi(t)-

    cosψi(t)sinφi(t))

    uiZ(t)=-Fi(t)Micosθi(t)cosφi(t)+g (3)

    式中: i=1, 2, …, N。 通過式(3)將動(dòng)力學(xué)模型(1)轉(zhuǎn)換為二階積分器形式:

    p·i(t)=vi(t)

    v·i(t)=ui(t)" (4)

    式中: i=1, 2, …, N。 除式(1)描述的外環(huán)位速動(dòng)力學(xué)模型外, 多旋翼無(wú)人機(jī)還具有表述空間轉(zhuǎn)動(dòng)性質(zhì)的內(nèi)環(huán)姿態(tài)動(dòng)力學(xué)方程, 在研究針對(duì)無(wú)人機(jī)蜂群編隊(duì)控制傳感器攻擊問題時(shí)不涉及無(wú)人機(jī)空間中轉(zhuǎn)動(dòng)情況的分析, 所以考慮將無(wú)人機(jī)視為質(zhì)點(diǎn)。 需要指出, 控制輸入ui(t)可以通過文獻(xiàn)[15]提出的位姿轉(zhuǎn)換方程得出期望的姿態(tài)角, 并利用姿態(tài)控制器轉(zhuǎn)化為無(wú)人機(jī)的控制力矩和升力。

    令xi=[pTi(t), vTi(t)]T∈6×1描述無(wú)人機(jī)的位置速度狀態(tài)信息, 式(4)可以表述為狀態(tài)空間形式:

    x·N(t)=AxN(t)+BuN(t)

    x·i(t)=Axi(t)+Bui(t), i∈F (5)

    式中: uN(t)∈×1為領(lǐng)導(dǎo)者無(wú)人機(jī)的控制輸入; ui(t)為跟隨者無(wú)人機(jī)控制輸入(i∈F); A和B分別表示蜂群系統(tǒng)的系統(tǒng)矩陣和輸入矩陣: A=0100I3,

    B=01I3。

    領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群系統(tǒng)期望的編隊(duì)隊(duì)形采用時(shí)變編隊(duì)向量f(t)=[fT1(t), fT2(t), …, fTN-1(t), 0]T表示, 其中fi(t)=[fTip(t), fTiv(t)]T∈6×1, i=1, 2, …, N為連續(xù)可微的時(shí)變向量, fi(t)(i∈F)表示跟隨者無(wú)人機(jī)i與編隊(duì)參考中心的相對(duì)偏移量, 記領(lǐng)導(dǎo)者無(wú)人機(jī)N編隊(duì)向量fN(t)=0。

    定義1: 對(duì)于任意有界初始狀態(tài), 如果領(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群系統(tǒng)滿足

    limt→∞(xi(t)-fi(t)-Td(t))=0 (i=1, 2, …, N)(6)

    則無(wú)人機(jī)蜂群系統(tǒng)實(shí)現(xiàn)領(lǐng)導(dǎo)跟隨結(jié)構(gòu)時(shí)變編隊(duì)跟蹤。

    本文以具有切換拓?fù)涞念I(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群為攻擊目標(biāo), 該無(wú)人機(jī)蜂群系統(tǒng)能夠?qū)崿F(xiàn)基于一致性的編隊(duì)跟蹤。 鑒于此, 可將領(lǐng)導(dǎo)者無(wú)人機(jī)控制輸入描述為

    uN(t)=K1(xN(t)-fN(t)-Td(t))(7)

    式中: K1∈3×6為增益矩陣。

    跟隨者無(wú)人機(jī)控制協(xié)議為[16

    ui(t)=K2∑N-1j=1wσ(t)ij(xi(t)-fi(t)-xj(t)+fj(t))+

    K2wσ(t)iN(xi(t)-fi(t)-xN(t))+δi(t)

    (8)

    式中: i∈F; K2∈ 6為增益矩陣; wσ(t)ij表示σ(t)對(duì)應(yīng)圖拓?fù)渲懈S者無(wú)人機(jī)j到i的通信權(quán)重; wσ(t)iN為領(lǐng)導(dǎo)者無(wú)人機(jī)N到跟隨者無(wú)人機(jī)i的通信權(quán)重。 控制協(xié)議中第一部分K2∑N-1j=1wσ(t)ij(xi(t)-fi(t)-xj(t)+fj(t))為跟隨者無(wú)人機(jī)之間編隊(duì)狀態(tài)差反饋?lái)?xiàng), 描述了相鄰無(wú)人機(jī)之間趨向編隊(duì)向量的吸引作用, 可使無(wú)人機(jī)狀態(tài)收斂至編隊(duì)向量。 第二部分K2wσ(t)iN(xi(t)-fi(t)-xN(t))為跟隨者無(wú)人機(jī)與領(lǐng)導(dǎo)者無(wú)人機(jī)之間編隊(duì)跟蹤狀態(tài)差的反饋?lái)?xiàng), 用于實(shí)現(xiàn)跟隨者無(wú)人機(jī)對(duì)領(lǐng)導(dǎo)者無(wú)人機(jī)的編隊(duì)跟蹤。 第三部分δi(t)∈3×1為編隊(duì)補(bǔ)償項(xiàng), 可以擴(kuò)大時(shí)變向量fi(t)的可行域。

    1.3 驅(qū)離反制方法描述

    為避免傳統(tǒng)無(wú)人機(jī)蜂群反制手段帶來(lái)的次生災(zāi)害, 本文通過設(shè)計(jì)傳感器攻擊策略, 將攻擊信號(hào)注入蜂群關(guān)鍵個(gè)體中, 利用無(wú)人機(jī)蜂群分布式網(wǎng)絡(luò)交互機(jī)制, 在不破壞無(wú)人機(jī)蜂群編隊(duì)隊(duì)形的前提下, 誘騙整個(gè)無(wú)人機(jī)蜂群脫離原有目標(biāo)路徑, 而沿攻擊者預(yù)設(shè)的軌跡運(yùn)動(dòng)。 攻擊者能夠通過調(diào)節(jié)攻擊信號(hào)對(duì)驅(qū)離后的無(wú)人機(jī)蜂群運(yùn)動(dòng)軌跡進(jìn)行設(shè)計(jì)。

    假設(shè)2: 攻擊者能夠獲取無(wú)人機(jī)蜂群的編隊(duì)隊(duì)形、 位置速度狀態(tài)、 拓?fù)浣换ズ涂刂茀?shù)信息K。

    注釋1: 通過光電探測(cè)設(shè)備獲取實(shí)時(shí)位置速度狀態(tài)等技術(shù)手段相對(duì)成熟。 文獻(xiàn)[17]給出一種基于紅外探測(cè)的無(wú)人機(jī)蜂群結(jié)構(gòu)特征感知方法來(lái)檢測(cè)識(shí)別無(wú)人機(jī)蜂群隊(duì)形。 對(duì)于異構(gòu)無(wú)人機(jī)蜂群, 可以通過外部特征識(shí)別無(wú)人機(jī)攜帶的載荷類型來(lái)區(qū)分領(lǐng)導(dǎo)者無(wú)人機(jī), 這是因?yàn)闊o(wú)人機(jī)任務(wù)分工不同, 領(lǐng)導(dǎo)者無(wú)人機(jī)一般攜帶光學(xué)鏡頭、 雷達(dá)探測(cè)等感知類載荷, 而跟隨者無(wú)人機(jī)一般攜帶執(zhí)行任務(wù)相關(guān)的載荷, 如投擲裝置等。 對(duì)于同構(gòu)無(wú)人機(jī)蜂群, 可以通過運(yùn)動(dòng)特征識(shí)別領(lǐng)導(dǎo)者無(wú)人機(jī), 這是因?yàn)轭I(lǐng)導(dǎo)者無(wú)人機(jī)與跟隨者無(wú)人機(jī)控制協(xié)議不同, 領(lǐng)導(dǎo)者無(wú)人機(jī)的控制協(xié)議僅受任務(wù)目標(biāo)路徑Td(t)影響, 其運(yùn)動(dòng)軌跡相對(duì)平滑, 而跟隨者無(wú)人機(jī)在保持編隊(duì)隊(duì)形的同時(shí)還需要跟蹤領(lǐng)導(dǎo)者無(wú)人機(jī), 故無(wú)人機(jī)蜂群整體在進(jìn)行轉(zhuǎn)彎或加速等運(yùn)動(dòng)時(shí), 跟隨者無(wú)人機(jī)會(huì)出現(xiàn)一定的跟蹤滯后現(xiàn)象, 攻擊者可以通過分析無(wú)人機(jī)蜂群外部觀測(cè)數(shù)據(jù)來(lái)識(shí)別領(lǐng)導(dǎo)者。 關(guān)于拓?fù)浣换ズ涂刂茀?shù)K等信息的獲取, 文獻(xiàn)[18]提出了利用辨識(shí)和估計(jì)的方法能夠獲取未知群系統(tǒng)的拓?fù)湫畔ⅲ?文獻(xiàn)[19-20]針對(duì)已知系統(tǒng)矩陣和輸入矩陣的情況下, K作為狀態(tài)反饋增益矩陣時(shí), 利用狀態(tài)空間模型參數(shù)相似變化關(guān)系, 構(gòu)造相似變換增益矩陣, 基于遞推子空間方法提出一種辨識(shí)和校驗(yàn)反饋增益矩陣的遞推算法, 利用最小二乘法計(jì)算得到反饋增益矩陣K。 因此假設(shè)2具有相對(duì)合理性, 具備一定的實(shí)際應(yīng)用意義。

    2 針對(duì)無(wú)人機(jī)蜂群編隊(duì)跟蹤的傳感器攻擊設(shè)計(jì)與分析

    2.1 針對(duì)領(lǐng)導(dǎo)者無(wú)人機(jī)的傳感器攻擊策略

    領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群在執(zhí)行任務(wù)中, 領(lǐng)導(dǎo)者無(wú)人機(jī)帶領(lǐng)跟隨者無(wú)人機(jī)形成編隊(duì)并沿任務(wù)目標(biāo)路徑運(yùn)動(dòng), 領(lǐng)導(dǎo)者無(wú)人機(jī)可作為關(guān)鍵攻擊目標(biāo)。 攻擊信號(hào)的注入方式可根據(jù)領(lǐng)導(dǎo)者無(wú)人機(jī)傳感器類型進(jìn)行調(diào)整。

    考慮向領(lǐng)導(dǎo)者無(wú)人機(jī)傳感器注入攻擊信號(hào), 其攻擊模型可以建立為

    x^N(t)=xN(t)+βxαN(t)(9)

    式中: x^N(t)∈表示領(lǐng)導(dǎo)者無(wú)人機(jī)被攻擊后的狀態(tài)信息; xαN(t)∈表示注入的攻擊信號(hào); β為攻擊信號(hào)注入標(biāo)識(shí), β=0表示未注入誘導(dǎo)攻擊信號(hào); β=1表示注入誘導(dǎo)攻擊信號(hào)。

    攻擊注入后, 式(7)轉(zhuǎn)換為

    u^N(t)=K1(xN(t)+βxαN(t)-Td(t))(10)

    式中: u^N(t)為攻擊后領(lǐng)導(dǎo)者無(wú)人機(jī)的實(shí)際控制輸入, 回顧上節(jié)可知fN(t)=0。

    2.2 針對(duì)無(wú)人機(jī)蜂群編隊(duì)跟蹤的傳感器攻擊分析

    如果假設(shè)1成立, 則領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群的通信拓?fù)銵aplacian矩陣可以表示為

    σ(t)=σ(t)ff+σ(t)lflσ(t)lf

    01×(N-1)0(11)

    式中: σ(t)ff∈-1)×(N-1)是由跟隨者無(wú)人機(jī)組成的交互拓?fù)鋵?duì)應(yīng)的Laplacian矩陣;" lσ(t)lf=-[w1N, w2N, …, w(N-1)N] T∈×1可以用來(lái)描述領(lǐng)導(dǎo)者無(wú)人機(jī)與跟隨者無(wú)人機(jī)之間的交互作用。 式(7)和式(8)表明領(lǐng)導(dǎo)者無(wú)人機(jī)不受任何跟隨者無(wú)人機(jī)的影響, 但領(lǐng)導(dǎo)者無(wú)人機(jī)的協(xié)調(diào)狀態(tài)信息可以通過lσ(t)lf傳遞給部分跟隨者無(wú)人機(jī)。 式(11)中, σ(t)lf=diag{w1N, w2N, …, w(N-1)N}表示對(duì)角線元素為wiN的對(duì)角矩陣, 其中i=1, 2, …, N-1。 需要注意, σ(t)ff+σ(t)lf是一個(gè)對(duì)稱正定矩陣21。

    因?yàn)檩斎刖仃嘊為列滿秩矩陣, 所以存在一個(gè)非奇異矩陣B=[B~T, B-T]T, 使得B~B=I3, B-B=0, 其中B~∈瘙綆3×6, B-∈令

    f(t)=[fTF(t), 0]T, fTF(t)∈

    x(t)=[xT1(t), xT2(t), …, xTN-1(t), xTN(t)]T=

    [xTF, xTN(t)]T

    δ(t)=[δT1(t), δT2(t), …, δTN-1(t), 0]T=

    [δTF(t), 0]T

    可得攻擊信號(hào)注入后無(wú)人機(jī)蜂群系統(tǒng)全局閉環(huán)形式:

    x·(t)=IN-1A00A+BK1+σ(t)BK2x(t)-

    σ(t)ff+σ(t)lf)BK2

    0fF(t)+

    0BK1(βxαN(t)-

    Td(t))+

    IN-1B0δF(t) (12)

    令θi(t)=xi(t)-fi(t), i=1, 2, …, N; θ(t)=[θTF(t), θTN(t)]T, θ(t)=[θT1(t), θT2(t), …, θTN-1(t), θTN(t)]T, 式(12)可以進(jìn)一步表示為

    θ·(t)=IN-1A0

    0A+BK1+σ(t)BK2θ(t)+

    IN-1A0fF(t)-I6(N-1)0f·F(t)+

    0BK1(βxαN(t)-

    Td(t))+IN-1B0δF(t)(13)

    攻擊作用下, 式(12)仍能保持時(shí)變編隊(duì)f(t)當(dāng)且僅當(dāng)式(13)實(shí)現(xiàn)狀態(tài)一致。

    令Η=IN-11N-1

    01×(N-1)1, 則

    Η-1=IN-1-1N-1

    01×(N-1)1。

    根據(jù)引理1可知σ(t)1N=0, 即

    σ(t)ff+σ(t)lf)1N-1+lσ(t)lf=0(N-1)×1

    可得

    Η-1σ(t)Η=σ(t)ff+σ(t)lf0

    00(14)

    式(14)表明領(lǐng)導(dǎo)者無(wú)人機(jī)的動(dòng)態(tài)可以通過非奇異變換矩陣Η與無(wú)人機(jī)蜂群系統(tǒng)的整個(gè)動(dòng)態(tài)解耦。

    令ζ(t)=[ξT1(t), ξT2(t), …, ξTN-1(t)]T, ξ(t)=

    [ζT(t), θTN(t)]T, 則跟隨者無(wú)人機(jī)與領(lǐng)導(dǎo)者無(wú)人機(jī)的跟蹤狀態(tài)差為

    ζ(t)=(IN-1I6)θF(t)-(1N-1I6)xN(t) (15)

    式(15)可由下式得到:

    ξ(t)=(Η-1I6)θ(t)

    將θ(t)=(ΗI6)ξ(t)代入式(13), 等式兩邊左乘Η-1I6并將其分解為編隊(duì)跟蹤參考子系統(tǒng):

    x·N(t)=(A+BK1)xN(t)+BK1(βxαN(t)-Td(t))(16)

    編隊(duì)跟蹤誤差子系統(tǒng):

    ζ·(t)=(IN-1A+(σ(t)ff+σ(t)lf)BK2)ζ(t)+

    (-1N-1BK1)(xN+βxαN(t)-Td(t))+

    (IN-1A)fF(t)-(IN-1I6)f·F(t)+

    (IN-1B)δF(t)(17)

    由假設(shè)1可知, 所有跟隨者無(wú)人機(jī)之間的通信交互是無(wú)向切換的, 所以σ(t)ff為對(duì)稱矩陣, σ(t)lf為非零矩陣, 故(σ(t)ff+σ(t)lf)T=σ(t)ff+σ(t)lfgt;0, 即其特征值均為正實(shí)數(shù)。 可知存在一個(gè)正交矩陣Uσ(t)1)×(N-1))使得

    (Uσ(t))T(σ(t)ff+σ(t)lf)Uσ(t)=diag{λσ(t)1, λσ(t)2, …, λσ(t)N-1},

    0lt;λσ(t)1≤λσ(t)2≤…≤λσ(t)N-1為σ(t)ff+σ(t)lf的特征值。

    令λ-1=min{λm1:k∈{1, 2, …, s}}表示所有拓?fù)浼现蟹橇闾卣髦档淖钚≈担?λ-N-1=max{λmN-1:k∈{1, 2, …, s}}表示所有拓?fù)浼现刑卣髦底畲笾怠?/p>

    令ζ(t)=(Uσ(t)I6)ζ(t), 進(jìn)一步地將式(17)兩邊左乘(Uσ(t))TI6轉(zhuǎn)化為

    ζ·i(t)=(A+λσ(t)iBK2)ζi(t)+(eTi(Uσ(t))TB)·

    δF(t)+(eTi(Uσ(t))T(-1N-1)BK1)·

    (xN+βxαN(t)-Td(t))+(eTi(Uσ(t))TA)·

    fF(t)-(eTi(Uσ(t))TI6)f·F(t)

    (18)

    式中: i=1, 2, …N-1; ei∈瘙綆(N-1)×1為第i個(gè)元素為1其余元素為0的N-1維列向量。

    注釋2: 上述分析利用狀態(tài)空間分解法將攻擊下的領(lǐng)導(dǎo)者無(wú)人機(jī)動(dòng)力學(xué)系統(tǒng)與跟隨者無(wú)人機(jī)動(dòng)力學(xué)系統(tǒng)解耦, 并采用兩步非奇異變換法將子系統(tǒng)變換為低維子系統(tǒng)。 具體的, 第一步利用非奇異變換將Laplacian矩陣轉(zhuǎn)換為對(duì)角塊矩陣, 第二步將矩陣σ(t)ff+σ(t)lf對(duì)角化, 使其可以被線性解耦。 需要指出, 子系統(tǒng)描述了跟隨者無(wú)人機(jī)之間的相對(duì)運(yùn)動(dòng)關(guān)系, 無(wú)人機(jī)蜂群若要實(shí)現(xiàn)定義1所描述的編隊(duì)跟蹤需要子系統(tǒng)是漸近穩(wěn)定的。

    下面給出領(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群編隊(duì)能夠被驅(qū)離反制的分析判據(jù)。

    定理1: 對(duì)于任意給定初始狀態(tài), 當(dāng)下列條件均成立時(shí), 具有切換拓?fù)涞臒o(wú)人機(jī)蜂群系統(tǒng)可以被驅(qū)離反制:

    ① 編隊(duì)可行性條件:

    limt→∞(B-AfF(t)-B-f·F(t))=0

    ② 存在矩陣P=PTgt;0使得

    PA+ATP+λ-k(PBK2+KT2BTP)lt;0, k∈{1, N-1}

    ③ 跟蹤可行性條件:

    limt→∞(ATd(t)-T·d(t))=0且A+BK1是Hurwitz的

    ④ 攻擊可行性條件:

    limt→∞β(AxαN(t)-x·αN(t))=0

    證明: 定義如下輔助函數(shù):

    C(t)=(HI6)[xTN(t), 0]T(19)

    C-(t)=(HI6)[0, ζT(t)]T(20)

    可知[xTN(t), 0]T=e1xN(t), 由式(19)可知

    C(t)=He1xN(t)=1NxN(t) (21)

    已知HI6是非奇異的, C(t)和C-(t)是線性無(wú)關(guān)的, 且θ(t)=C(t)+C-(t), 因此, 當(dāng)且僅當(dāng)limt→∞ζ(t)=0時(shí), 無(wú)人機(jī)蜂群中跟隨者能實(shí)現(xiàn)編隊(duì)fF(t)。 因?yàn)閁I6為非奇異矩陣, 上述條件等價(jià)于

    limt→∞ζi(t)=0, i=1, 2, …, N-1(22)

    當(dāng)由fF(t)確定的時(shí)變編隊(duì)可行, 即條件①被滿足, 結(jié)合B-B=0, 則有

    limt→∞(B-AfF(t)-B-f·F(t)+B-BδF(t))=0(23)

    同時(shí)可以找到δF(t)使得

    limt→∞(B~AfF(t)-B~f·F(t)+B~BδF(t))=0(24)

    結(jié)合式(23)~(24)可知

    limt→∞(BAfF(t)-Bf·F(t)+BBδF(t))=0 (25)

    將式(25)左乘B-1

    limt→∞(AfF(t)-f·F(t)+BδF(t))=0(26)

    接下來(lái)考慮子系統(tǒng)ζ-·i(t)=(A+λσ(t)iBK2)ζ-i(t), i∈F的穩(wěn)定性問題。

    選取如下共同Lyapunov候選函數(shù):

    Vi(t)=ζ-Ti(t)Pζ-i(t), i∈F(27)

    對(duì)式(27)兩邊求導(dǎo)得

    V·i(t)=ζ-Ti(t)(PA+ATP+λσ(t)i(PBK2+KT2BTP))ζ-i(t)(28)

    如果條件②被滿足, 即下式成立:

    PA+ATP+λ-k(PBK2+KT2BTP)lt;0, k∈{1, N-1}(29)

    根據(jù)線性矩陣不等式的凸屬性, 式(29)成立可以保證PA+ATP+λσ(t)i(PBK2+KT2BTP)lt;0, i∈F, 此時(shí)子系統(tǒng)ζ-·i(t)=(A+λσ(t)iBK2)ζ-i(t), i∈F是漸近穩(wěn)定的。

    考慮傳感器攻擊下領(lǐng)導(dǎo)者無(wú)人機(jī)跟蹤穩(wěn)定性問題。

    傳感器攻擊下領(lǐng)導(dǎo)者無(wú)人機(jī)相對(duì)于任務(wù)目標(biāo)路徑偏差記為

    μ(t)=xN(t)-βxαN(t)-Td(t)(30)

    將式(30)對(duì)時(shí)間求導(dǎo), 代入式(17)可得

    μ·(t)=(A+BK1)μ(t)+βx·αN(t)-AβxαN(t)+

    ATd(t)-T·d(t)(31)

    如果條件③和條件④被滿足, 則有

    limt→∞μ(t)=0(32)

    可知系統(tǒng)式(31)是漸近穩(wěn)定的。

    綜合式(26)、 (29)和式(32)可知,

    limt→∞ζi(t)=0, i=1, 2, …, N-1成立, 故子系統(tǒng)式(18)是漸近穩(wěn)定的。 定理1證畢。

    注釋3: 定理1給出了切換拓?fù)湎聼o(wú)人機(jī)蜂群能夠被驅(qū)離反制的充分條件, 條件①~③是傳感器攻擊的前提條件, 判斷了切換拓?fù)湎骂I(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群能夠形成編隊(duì)并跟蹤任務(wù)目標(biāo)路徑Td(t), 其中條件①編隊(duì)可行性條件約束了由fF(t)描述的編隊(duì)隊(duì)形可行域, 條件②利用線性矩陣不等式條件保證了切換系統(tǒng)的穩(wěn)定性, 條件③確保了誘騙領(lǐng)導(dǎo)者無(wú)人機(jī)跟蹤攻擊后的任務(wù)目標(biāo)路徑。 不受條件①~③的影響, 條件④給出了攻擊信號(hào)的可行性約束。 需要指出, 傳感器攻擊信號(hào)與無(wú)人機(jī)蜂群動(dòng)力學(xué)系統(tǒng)具有相同的模態(tài), 這表明攻擊信號(hào)注入后不會(huì)影響無(wú)人機(jī)蜂群系統(tǒng)編隊(duì)收斂性。

    注釋4:" 由式(7)和式(30)可知, 傳感器攻擊前后領(lǐng)導(dǎo)者無(wú)人機(jī)狀態(tài)軌跡由跟蹤任務(wù)目標(biāo)路徑Td(t)改變?yōu)楦櫐聏αN(t)+Td(t), 攻擊者可通過調(diào)節(jié)攻擊信號(hào)xαN(t)改變領(lǐng)導(dǎo)者無(wú)人機(jī)跟蹤軌跡, 來(lái)實(shí)現(xiàn)對(duì)無(wú)人機(jī)蜂群誘騙后狀態(tài)軌跡的設(shè)計(jì)。 需要說明, 領(lǐng)導(dǎo)跟隨結(jié)構(gòu)無(wú)人機(jī)蜂群的運(yùn)動(dòng)軌跡由領(lǐng)導(dǎo)者無(wú)人機(jī)決定, 當(dāng)領(lǐng)導(dǎo)者無(wú)人機(jī)被攻擊驅(qū)離時(shí), 整個(gè)無(wú)人機(jī)蜂群將沿著領(lǐng)導(dǎo)者的運(yùn)動(dòng)軌跡移動(dòng)。

    注釋5:" 以導(dǎo)航欺騙信號(hào)注入為例, 考慮領(lǐng)導(dǎo)者無(wú)人機(jī)帶領(lǐng)蜂群以“一”字形編隊(duì)襲擊機(jī)場(chǎng)、 車站等重要公共基礎(chǔ)設(shè)施。 蜂群飛行高度200 m, 間隔距離10 m, 蜂群中無(wú)人機(jī)均采用GPS衛(wèi)星導(dǎo)航定位系統(tǒng)。 在此場(chǎng)景下, 安保機(jī)構(gòu)可以通過分析無(wú)人機(jī)之間的間隔、 與發(fā)射天線的距離和發(fā)射波束寬度等幾何條件關(guān)系, 設(shè)計(jì)適當(dāng)驅(qū)離反制方案, 實(shí)現(xiàn)攻擊信號(hào)的單獨(dú)注入。 如圖1所示, 利用發(fā)射波束小于6°的定向天線, 能夠?qū)⒐粜盘?hào)單獨(dú)注入領(lǐng)導(dǎo)者無(wú)人機(jī)。 需要說明, 在實(shí)際應(yīng)用中攻擊信號(hào)的注入需要根據(jù)無(wú)人機(jī)傳感器類型進(jìn)行調(diào)整, 例如針對(duì)采用視覺傳感器的無(wú)人機(jī), 圖1的方案顯然不再適用。

    3 仿真分析

    3.1 7架無(wú)人機(jī)繞圓編隊(duì)仿真案例

    考慮由7架無(wú)人機(jī)組成的領(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群沿任務(wù)目標(biāo)路徑Td(t)運(yùn)動(dòng), 其中1架為領(lǐng)導(dǎo)者無(wú)人機(jī), 其余6架為跟隨者無(wú)人機(jī)。 無(wú)人機(jī)蜂群中領(lǐng)導(dǎo)者無(wú)人機(jī)單向傳輸協(xié)同信息至跟隨者無(wú)人機(jī), 且不接收任何跟隨者無(wú)人機(jī)信息, 跟隨者無(wú)人機(jī)之間采用無(wú)向交互方式, 無(wú)人機(jī)蜂群通信以隨機(jī)切換方式進(jìn)行, 其通信拓?fù)浼先鐖D2所示。

    無(wú)人機(jī)蜂群中領(lǐng)導(dǎo)者無(wú)人機(jī)跟蹤任務(wù)目標(biāo)路徑Td(t)滿足定理1中條件③, 取Td(0)=[0, 0, 8, 3, 0, 0]T。 領(lǐng)導(dǎo)者無(wú)人機(jī)狀態(tài)初始值為

    x1(0)=[0, 0, 0, 0, 0, 0]T

    其余跟隨者無(wú)人機(jī)狀態(tài)初始值為

    x2(0)=[2, 6, 0, 0, 0, 0]T

    x3(0)=[4, -6, 0, 0, 0, 0]T

    x4(0)=[5, -4, 0, 0, 0, 0]T

    x5(0)=[5, -3, 0, 0, 0, 0]T

    x6(0)=[3, 4, 0, 0, 0, 0]T

    x7(0)=[4, 7, 0, 0, 0, 0]T

    其中6架跟隨者無(wú)人機(jī)采用正六邊形隊(duì)形以領(lǐng)導(dǎo)者無(wú)人機(jī)為中心繞圓前進(jìn), 采用的時(shí)變編隊(duì)向量為

    fi(t)=[3cos((i-1)π/3+t)-3cos((i-1)π/3+t)·3sin((i-1)π/3+t)-3sin((i-1)π/3+t)·3sin((i-1)π/3+t)3cos((i-1)π/3+t)]T

    式中:i=2, 3, …, 7; tlt;50 s。 無(wú)人機(jī)蜂群在編隊(duì)跟蹤運(yùn)動(dòng)時(shí), 拓?fù)淝袚Q次序如圖3所示。

    當(dāng)注入的傳感器攻擊信號(hào)滿足定理1條件④, 初始值為xα1(0)=[0, 0, 0, 0, -0.2, -0.1]T。 仿真時(shí)長(zhǎng)為50 s, 攻擊信號(hào)注入時(shí)刻為第30 s, 攻擊效果如圖4所示。 圖4中, 紅色五角星表示領(lǐng)導(dǎo)者無(wú)人機(jī), 其運(yùn)動(dòng)路線由紅色實(shí)線表示, 跟隨者無(wú)人機(jī)成正六邊形圍繞領(lǐng)導(dǎo)者無(wú)人機(jī)繞圓運(yùn)動(dòng)。 未注入攻擊時(shí), 無(wú)人機(jī)蜂群整體跟蹤任務(wù)目標(biāo)路徑運(yùn)動(dòng), 在第30 s注入攻擊后, 領(lǐng)導(dǎo)者無(wú)人機(jī)運(yùn)動(dòng)偏離任務(wù)目標(biāo)路徑, 跟隨者無(wú)人機(jī)仍保持編隊(duì)隊(duì)形跟蹤領(lǐng)導(dǎo)者無(wú)人機(jī), 無(wú)人機(jī)蜂群整體運(yùn)動(dòng)軌跡被驅(qū)離。

    圖5~7為傳感器攻擊下無(wú)人機(jī)蜂群在三維空間X, Y, Z三軸的速度偏差變化情況。 可以看出, 攻擊信號(hào)初始值是分別在Y軸速度和Z軸速度方向施加-0.2 m/s和-0.1 m/s的速度攻擊, 圖6和圖7分別在攻擊注入時(shí)刻第30 s發(fā)生變化后又逐漸收斂, 表明無(wú)人機(jī)蜂群重新恢復(fù)原有編隊(duì)隊(duì)形, 驗(yàn)證了符合定理1中條件④約束的傳感器攻擊信號(hào)不會(huì)破壞無(wú)人機(jī)蜂群的編隊(duì)隊(duì)形。

    3.2 9架無(wú)人機(jī)收縮編隊(duì)仿真案例

    考慮由9架無(wú)人機(jī)組成的領(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群, 其中1架為領(lǐng)導(dǎo)者無(wú)人機(jī), 其余8架為跟隨者無(wú)人機(jī)。 無(wú)人機(jī)蜂群通信拓?fù)浼先鐖D8所示。

    本仿真案例中拓?fù)淝袚Q次序參照?qǐng)D3, 任務(wù)目標(biāo)路徑、 傳感器攻擊信號(hào)設(shè)計(jì)和攻擊注入時(shí)間等參數(shù)設(shè)置與3.1節(jié)仿真案例保持一致, 仿真時(shí)長(zhǎng)為50 s,

    增加的第8和第9架跟隨者無(wú)人機(jī)的狀態(tài)初始值設(shè)置為

    x8(0)=[2, 1, 0, 0, 0, 0]T

    x9(0)=[3, 3, 0, 0, 0, 0]T

    無(wú)人機(jī)蜂群中跟隨者無(wú)人機(jī)采用的時(shí)變編隊(duì)向量為

    fi(t)=[(5+3sin(t))×sin((i-1)π/4)(5+

    3sin(t))×sin((i-1)π/4)(5+3sin(t))×

    cos((i-1)π/4)3sin((i-1)π/4)×

    cos(t)3sin((i-1)π/4)×

    cos(t)3cos((i-1)π/4)×cos(t)]T

    式中: i=2, 3, …, 9; tlt;50 s。

    圖9展示了有無(wú)傳感器攻擊的無(wú)人機(jī)蜂群編隊(duì)運(yùn)動(dòng)對(duì)比圖。 從圖9(a)中可以看出, 在未向領(lǐng)導(dǎo)者無(wú)人機(jī)注入傳感器攻擊時(shí)無(wú)人機(jī)蜂群整體成正八邊形的編隊(duì)包圍領(lǐng)導(dǎo)者無(wú)人機(jī)運(yùn)動(dòng), 運(yùn)動(dòng)過程中編隊(duì)半徑不斷收縮。 在圖9(b)中領(lǐng)導(dǎo)者無(wú)人機(jī)被注入傳感器攻擊信號(hào)后, 領(lǐng)導(dǎo)者無(wú)人機(jī)帶領(lǐng)8架跟隨者無(wú)人機(jī)整體偏離了原有的任務(wù)目標(biāo)路徑, 同時(shí)可以發(fā)現(xiàn), 在無(wú)人機(jī)蜂群整體偏離運(yùn)動(dòng)過程中, 編隊(duì)中各無(wú)人機(jī)的相對(duì)方位, 以及由編隊(duì)函數(shù)確定的正八邊形半徑收縮速率都沒有發(fā)生變化, 說明注入的傳感器攻擊信號(hào)不會(huì)改變無(wú)人機(jī)蜂群的編隊(duì)隊(duì)形運(yùn)動(dòng), 即定理1中條件①與條件④是非耦合的獨(dú)立條件。

    4 結(jié)" 論

    本文研究了針對(duì)領(lǐng)導(dǎo)跟隨結(jié)構(gòu)分布式多旋翼無(wú)人機(jī)蜂群的驅(qū)離反制問題。 從攻擊者角度, 提出了一種基于傳感器攻擊的蜂群驅(qū)離反制方法, 給出了無(wú)人機(jī)蜂群編隊(duì)驅(qū)離的分析判據(jù)。 通過將攻擊信號(hào)注入領(lǐng)導(dǎo)者無(wú)人機(jī), 利用無(wú)人機(jī)蜂群分布式交互特點(diǎn), 在不改變無(wú)人機(jī)蜂群編隊(duì)隊(duì)形的情況下, 實(shí)現(xiàn)無(wú)人機(jī)蜂群的整體驅(qū)離。 仿真結(jié)果表明, 當(dāng)傳感器攻擊信號(hào)滿足攻擊可行條件時(shí), 攻擊者能夠?qū)崿F(xiàn)無(wú)人機(jī)蜂群整體的驅(qū)離反制。

    參考文獻(xiàn):

    [1] 陳麗蓉," 羅俊松," 寧進(jìn)," 等. 無(wú)人機(jī)集群智能反制技術(shù)的研究現(xiàn)狀及展望[J]. 無(wú)線電工程," 2024," 54(9): 2181-2192.

    Chen Lirong," Luo Junsong," Ning Jin," et al. Research Status and Prospects of Intelligent Countermeasure Technology for UAV Swarms[J]. Radio Engineering," 2024," 54(9): 2181-2192.(in Chinese)

    [2] 李威," 張艷紅," 李冰," 等. 分布式無(wú)人機(jī)蜂群導(dǎo)航欺騙反制方法研究[J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù)," 2023(3): 95-105.

    Li Wei," Zhang Yanhong," Li Bing," et al. Research on Navigation Deception Countermeasures for Distributed UAV Swarm[J]. Tactical Missile Technology," 2023(3): 95-105.(in Chinese)

    [3] 李姝," 裘昌利," 李晶. 攻防一體化無(wú)人機(jī)蜂群作戰(zhàn)體系研究[C]∥無(wú)人系統(tǒng)高峰論壇," 2022," 7-12.

    Li Shu," Qiu Changli," Li Jing. Research on Attack-Defense Integration of UAV Swarms[C]∥Unmanned Systems Summit," 2022," 7-12. (in Chinese)

    [4] Lyu C Y," Zhan R J. Global Analysis of Active Defense Technologies for Unmanned Aerial Vehicle[J]. IEEE Aerospace and Electronic Systems Magazine," 2022," 37(1): 6-31.

    [5] 張海龍," 楊榮強(qiáng)," 王玉茜. 智能集群對(duì)防空作戰(zhàn)影響及反制策略分析[J]. 空天防御," 2021," 4(1): 26-32.

    Zhang Hailong," Yang Rongqiang," Wang Yuqian. Analysis of the Influence of Intelligent Swarming on Air Defense Operations and Countermeasures[J]. Air amp; Space Defense," 2021," 4(1): 26-32.(in Chinese)

    [6] Kerns A J," Shepard D P," Bhatti J A," et al. Unmanned Aircraft Capture and Control via GPS Spoofing[J]. Journal of Field Robotics," 2014," 31(4): 617-636.

    [7] 郭妍," 唐康華," 張鷺. 面向無(wú)人機(jī)的逐點(diǎn)偏移式衛(wèi)星導(dǎo)航欺騙干擾方法[J]. 工程科學(xué)與技術(shù)," 2023," 55(2): 275-284.

    Guo Yan," Tang Kanghua," Zhang Lu. GNSS Navigation Spoofing Method of UAV Based on Point-by-Point Offset[J]. Advanced Engineering Sciences," 2023," 55(2): 275-284.(in Chinese)

    [8] Noh J," Kwon Y," Son Y," et al. Tractor Beam[J]. ACM Transactions on Privacy and Security," 2019," 22(2): 1-26.

    [9] Ceccato M," Formaggio F," Tomasin S. Spatial GNSS Spoofing Against Drone Swarms with Multiple Antennas and Wiener Filter[J]. IEEE Transactions on Signal Processing," 2020," 68: 5782-5794.

    [10] Yang Y," Xiao Y," Li T S. Attacks on Formation Control for Multiagent Systems[J]. IEEE Transactions on Cybernetics," 2022," 52(12): 12805-12817.

    [11] Li J L," Wang L," Xi J X," et al. Induced Attack on Formation Control of Multiagent Systems with Prescribed Reference Trajectories[J]. International Journal of Robust and Nonlinear Control," 2024," 34(12): 8374-8397.

    [12] 李威," 王樂," 高久安," 等. 面向領(lǐng)導(dǎo)跟隨無(wú)人機(jī)蜂群的定向驅(qū)離攻擊方法[J]. 航空學(xué)報(bào)," 2024," 45(12): 190-210.

    Li Wei," Wang Le," Gao Jiuan," et al. Directional Expelling Attack on Unmanned Aerial Vehicle Swarm with Leader-Follower Structure[J]. Acta Aeronautica et Astronautica Sinica," 2024," 45(12): 190-210.(in Chinese)

    [13] Godsil C D," Royle G. Algebraic Graph Theory[M]. New York: Springer Verlag," 2001.

    [14] Ren W," Beard R W. Consensus Seeking in Multiagent Systems under Dynamically Changing Interaction Topologies[J]. IEEE Transactions on Automatic Control," 2005," 50(5): 655-661.

    [15] Xi J X," Wang C," Wang L," et al. Optimal Formation Strategy for Interconnected Quadrotor Drones: Theory and Experiment[J]. Aerospace Science and Technology," 2024," 144: 108771.

    [16] Xi J X," Wang L," Zheng J F," et al. Energy-Constraint Formation for Multiagent Systems with Switching Interaction Topologies[J]. IEEE Transactions on Circuits and Systems I: Regular Papers," 2020," 67(7): 2442-2454.

    [17] 夏文新," 楊小岡," 席建祥," 等. 基于紅外探測(cè)的無(wú)人機(jī)群結(jié)構(gòu)特性感知方法[J]. 紅外與激光工程," 2024," 53(1): 257-268.

    Xia Wenxin," Yang Xiaogang," Xi Jianxiang," et al. Structure Chara-cteristics Sensing Method of Unmanned Aerial Vehicle Group Based on Infrared Detection[J]. Infrared and Laser Engineering," 2024," 53(1): 257-268.(in Chinese)

    [18] Zheng Y," Wu X Q," Fan Z Y," et al. Identifying Topology and System Parameters of Fractional-Order Complex Dynamical Networks[J]. Applied Mathematics and Computation," 2022," 414: 126666.

    [19] 倪智宇," 譚述君," 吳志剛. 線性時(shí)變反饋系統(tǒng)的狀態(tài)空間模型和模態(tài)參數(shù)遞推辨識(shí)[J]. 控制與決策," 2016," 31(2): 324-330.

    Ni Zhiyu," Tan Shujun," Wu Zhigang. Recursive Identification of State Space Model and Modal Parameter of Linear Time-Varying Feedback System[J]. Control and Decision," 2016," 31(2): 324-330.(in Chinese)

    [20] 倪智宇. 大型柔性航天器時(shí)變動(dòng)力學(xué)參數(shù)在軌辨識(shí)方法[D]. 大連: 大連理工大學(xué)," 2016.

    Ni Zhiyu. On-Orbit Identification Method of Time-Varying Dynamic Parameters of Large Flexible Spacecraft[D]. Dalian: Dalian University of Technology," 2016. (in Chinese)

    [21] Li Z K, Duan Z S. Cooperative Control of Multi-Agent Systems: A Consensus Region Approach[M]. Boca Raton: CRC Press," 2015.

    Leader-Following Multirotor UAV Swarm Expulsion

    Method Based on Sensor Attacks

    Xu Zhaosheng1, 2," Wang Le1, 3*," Liu Nanchi1," Zhao Liangyu3," Xi Jianxiang1," Lin Haoshen1

    (1. Rocket Force University of Engineering," Xi’an 710025," China; 2. Unit 63768 of PLA," Xi’an 710025," China;

    3. Beijing Institute of Technology," Beijing 100081," China)

    Abstract: Studying countermeasures for leader-following multirotor UAV swarms from an attacker’s perspective. Proposing a sensor attack method against the leader UAV based on the root node attributes," and using swarm distributed network interaction characteristics of the swarm to expulse it to the path desired by the attacker. Firstly," the sensor attack model is constructed by using the dynamic characteristics of the UAV," and the attack signal is induced into the swarm formation tracking control process. Secondly," considering the swarm communication topology switching," a state space decomposition method is used to decompose the attack signal into a formation tracking reference subsystem and a formation tracking error subsystem. Sufficient conditions for realizing swarm repellent countermeasures under sensor attacks are given," and the analytical criterion is simplified by a two-step non-singular transformation. Finally," the effectiveness of the proposed countermeasures is verified by numerical simulation.

    Key words: sensor attacks; leader-following; UAVs countermeasures; multirotor swarm; distributed network

    成年人午夜在线观看视频| 亚洲色图综合在线观看| 日韩制服丝袜自拍偷拍| 啦啦啦视频在线资源免费观看| 久久青草综合色| 看免费av毛片| 精品第一国产精品| 国产1区2区3区精品| 亚洲欧美一区二区三区久久| 午夜免费鲁丝| 久久青草综合色| 脱女人内裤的视频| 日本欧美国产在线视频| 一级毛片电影观看| 午夜激情久久久久久久| 美女国产高潮福利片在线看| 久久亚洲精品不卡| 国产精品 国内视频| 日韩中文字幕视频在线看片| 免费日韩欧美在线观看| 日韩 欧美 亚洲 中文字幕| 成人亚洲欧美一区二区av| 18禁黄网站禁片午夜丰满| 久久99热这里只频精品6学生| 99热全是精品| 深夜精品福利| 亚洲美女黄色视频免费看| 丰满人妻熟妇乱又伦精品不卡| 女性被躁到高潮视频| 91成人精品电影| 午夜福利免费观看在线| 色婷婷av一区二区三区视频| 日韩免费高清中文字幕av| 极品少妇高潮喷水抽搐| 欧美xxⅹ黑人| 国产成人啪精品午夜网站| 一本大道久久a久久精品| 精品亚洲乱码少妇综合久久| 1024香蕉在线观看| 久久亚洲精品不卡| 久久久久久久久免费视频了| 各种免费的搞黄视频| 精品一区二区三区四区五区乱码 | 一本色道久久久久久精品综合| 国产三级黄色录像| 人体艺术视频欧美日本| 日本午夜av视频| 欧美日韩视频高清一区二区三区二| 侵犯人妻中文字幕一二三四区| h视频一区二区三区| 纵有疾风起免费观看全集完整版| 丰满人妻熟妇乱又伦精品不卡| 九色亚洲精品在线播放| 免费在线观看完整版高清| 日韩大码丰满熟妇| 9色porny在线观看| 成年av动漫网址| 五月开心婷婷网| 操出白浆在线播放| 少妇猛男粗大的猛烈进出视频| 少妇粗大呻吟视频| 人体艺术视频欧美日本| 国产精品免费大片| 在线观看www视频免费| 两个人免费观看高清视频| 国产精品久久久人人做人人爽| 成人三级做爰电影| 在线观看免费午夜福利视频| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美在线一区| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精| 久久精品成人免费网站| 国精品久久久久久国模美| 在线亚洲精品国产二区图片欧美| 久久久国产精品麻豆| 国产国语露脸激情在线看| 中文字幕高清在线视频| 十八禁高潮呻吟视频| 99国产综合亚洲精品| 国产成人91sexporn| 国语对白做爰xxxⅹ性视频网站| 涩涩av久久男人的天堂| 亚洲国产精品一区二区三区在线| 成人影院久久| 男人舔女人的私密视频| 亚洲国产精品一区二区三区在线| 男女边吃奶边做爰视频| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 亚洲伊人久久精品综合| 97人妻天天添夜夜摸| av视频免费观看在线观看| 亚洲精品日韩在线中文字幕| 久久久久久久久免费视频了| 日本五十路高清| 国产有黄有色有爽视频| av一本久久久久| 精品国产超薄肉色丝袜足j| 高清视频免费观看一区二区| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 午夜免费男女啪啪视频观看| 99热国产这里只有精品6| 精品国产国语对白av| 国产一区有黄有色的免费视频| 亚洲一码二码三码区别大吗| 我的亚洲天堂| av在线老鸭窝| 国产精品国产三级专区第一集| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三区在线| 欧美xxⅹ黑人| 校园人妻丝袜中文字幕| 亚洲图色成人| 精品久久久精品久久久| 日韩电影二区| 成人黄色视频免费在线看| 91精品国产国语对白视频| 亚洲精品一区蜜桃| 久久免费观看电影| 国产成人av激情在线播放| 久久久久久久大尺度免费视频| 欧美日本中文国产一区发布| 色综合欧美亚洲国产小说| 久久久久网色| 老汉色av国产亚洲站长工具| 亚洲欧美成人综合另类久久久| 丰满少妇做爰视频| 蜜桃在线观看..| 欧美成狂野欧美在线观看| 亚洲成人手机| 两个人免费观看高清视频| 国产精品一国产av| 不卡av一区二区三区| 少妇被粗大的猛进出69影院| 69精品国产乱码久久久| 国产精品久久久久久精品电影小说| 亚洲第一青青草原| 最黄视频免费看| 91九色精品人成在线观看| 免费少妇av软件| 这个男人来自地球电影免费观看| 天天躁日日躁夜夜躁夜夜| 亚洲第一av免费看| 成人黄色视频免费在线看| 午夜91福利影院| avwww免费| 只有这里有精品99| 2021少妇久久久久久久久久久| 少妇被粗大的猛进出69影院| 人人妻人人澡人人看| 老司机亚洲免费影院| kizo精华| 桃花免费在线播放| 九色亚洲精品在线播放| 国产高清不卡午夜福利| 欧美大码av| 午夜精品国产一区二区电影| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 激情视频va一区二区三区| 亚洲中文av在线| 精品欧美一区二区三区在线| 亚洲欧美激情在线| 狠狠婷婷综合久久久久久88av| 多毛熟女@视频| 纯流量卡能插随身wifi吗| 国产有黄有色有爽视频| 脱女人内裤的视频| 两个人免费观看高清视频| 亚洲一卡2卡3卡4卡5卡精品中文| 丝袜脚勾引网站| 亚洲精品美女久久久久99蜜臀 | 黄色 视频免费看| 99国产综合亚洲精品| 国产精品欧美亚洲77777| 欧美 日韩 精品 国产| 国产精品av久久久久免费| 午夜激情久久久久久久| 免费看av在线观看网站| 永久免费av网站大全| 黄片播放在线免费| 飞空精品影院首页| 国产免费又黄又爽又色| 国产成人免费无遮挡视频| 曰老女人黄片| 国产淫语在线视频| 欧美精品高潮呻吟av久久| 国产在线观看jvid| 青春草视频在线免费观看| 19禁男女啪啪无遮挡网站| 电影成人av| 少妇的丰满在线观看| 啦啦啦视频在线资源免费观看| 秋霞在线观看毛片| 亚洲精品久久成人aⅴ小说| 大陆偷拍与自拍| 亚洲一卡2卡3卡4卡5卡精品中文| 男女高潮啪啪啪动态图| 亚洲午夜精品一区,二区,三区| 女人高潮潮喷娇喘18禁视频| 在线天堂中文资源库| 国产一区二区三区av在线| 亚洲综合色网址| 亚洲午夜精品一区,二区,三区| 免费少妇av软件| 在线看a的网站| 亚洲av电影在线观看一区二区三区| 自线自在国产av| 久久久亚洲精品成人影院| 少妇猛男粗大的猛烈进出视频| 午夜免费观看性视频| 中文字幕人妻熟女乱码| 一边摸一边做爽爽视频免费| 国产女主播在线喷水免费视频网站| 亚洲,一卡二卡三卡| 亚洲国产看品久久| bbb黄色大片| 成年人午夜在线观看视频| 看免费成人av毛片| 国产精品一区二区在线不卡| av在线app专区| 久久99精品国语久久久| 曰老女人黄片| 国产伦人伦偷精品视频| 精品少妇内射三级| 亚洲精品美女久久av网站| 色综合欧美亚洲国产小说| 中文字幕制服av| 欧美在线一区亚洲| 亚洲国产最新在线播放| www.av在线官网国产| 熟女少妇亚洲综合色aaa.| 日韩制服丝袜自拍偷拍| 国产精品成人在线| 侵犯人妻中文字幕一二三四区| 成在线人永久免费视频| 国精品久久久久久国模美| 久热这里只有精品99| 国产成人91sexporn| 在线观看www视频免费| 大话2 男鬼变身卡| 一本久久精品| 亚洲中文av在线| 国产视频首页在线观看| 巨乳人妻的诱惑在线观看| 99精品久久久久人妻精品| 水蜜桃什么品种好| 少妇裸体淫交视频免费看高清 | 亚洲国产精品一区三区| 久热这里只有精品99| 色婷婷av一区二区三区视频| 男人舔女人的私密视频| 超色免费av| 国产精品久久久久成人av| 日韩,欧美,国产一区二区三区| 亚洲国产中文字幕在线视频| 一个人免费看片子| 亚洲免费av在线视频| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人| 大片免费播放器 马上看| 亚洲男人天堂网一区| 成人三级做爰电影| 国产在线观看jvid| 亚洲伊人久久精品综合| 晚上一个人看的免费电影| 亚洲一区二区三区欧美精品| 手机成人av网站| 老司机影院毛片| 人人妻人人添人人爽欧美一区卜| 91麻豆精品激情在线观看国产 | 亚洲欧美清纯卡通| 午夜福利在线免费观看网站| 久久 成人 亚洲| 国产淫语在线视频| 亚洲精品乱久久久久久| 91成人精品电影| 精品熟女少妇八av免费久了| 亚洲精品国产区一区二| 亚洲人成电影观看| 精品少妇一区二区三区视频日本电影| 欧美精品一区二区大全| 超碰成人久久| 久久人妻熟女aⅴ| 99国产精品99久久久久| 日本黄色日本黄色录像| 高潮久久久久久久久久久不卡| 成人午夜精彩视频在线观看| av网站在线播放免费| 新久久久久国产一级毛片| 婷婷色综合大香蕉| 老司机靠b影院| 1024视频免费在线观看| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 91精品伊人久久大香线蕉| 黄片小视频在线播放| 免费久久久久久久精品成人欧美视频| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 新久久久久国产一级毛片| 欧美日韩福利视频一区二区| 婷婷成人精品国产| 中文欧美无线码| 老司机在亚洲福利影院| 国产成人精品久久二区二区91| 天天躁日日躁夜夜躁夜夜| 亚洲男人天堂网一区| 男女高潮啪啪啪动态图| 国产精品久久久人人做人人爽| 色综合欧美亚洲国产小说| 91成人精品电影| 亚洲av成人精品一二三区| 免费看av在线观看网站| 制服诱惑二区| 黄色 视频免费看| 久久99精品国语久久久| 建设人人有责人人尽责人人享有的| 永久免费av网站大全| 婷婷色综合www| 波野结衣二区三区在线| 亚洲欧美一区二区三区久久| 99国产综合亚洲精品| 热re99久久国产66热| 国产欧美日韩精品亚洲av| 亚洲欧洲精品一区二区精品久久久| 国产精品香港三级国产av潘金莲 | 1024香蕉在线观看| 亚洲国产精品一区三区| 深夜精品福利| 一级毛片女人18水好多 | 咕卡用的链子| 叶爱在线成人免费视频播放| av电影中文网址| 波野结衣二区三区在线| 丰满少妇做爰视频| 肉色欧美久久久久久久蜜桃| 一本综合久久免费| 欧美xxⅹ黑人| 一边摸一边抽搐一进一出视频| 日本a在线网址| 一级毛片我不卡| 视频区欧美日本亚洲| 国精品久久久久久国模美| 国产欧美日韩精品亚洲av| 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 精品一区二区三区av网在线观看 | 日韩一区二区三区影片| 久久精品国产亚洲av高清一级| 亚洲 欧美一区二区三区| av欧美777| 美女主播在线视频| 久久 成人 亚洲| 久久国产精品男人的天堂亚洲| 国产亚洲精品久久久久5区| 99久久精品国产亚洲精品| 麻豆乱淫一区二区| 亚洲国产欧美一区二区综合| 亚洲精品第二区| 视频在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 色精品久久人妻99蜜桃| 亚洲精品成人av观看孕妇| a级毛片在线看网站| 国产精品久久久久久精品电影小说| 大型av网站在线播放| 老司机影院成人| 国产成人精品久久二区二区91| 欧美精品啪啪一区二区三区 | 一区二区三区精品91| 亚洲熟女毛片儿| 十八禁人妻一区二区| 黄网站色视频无遮挡免费观看| 婷婷色麻豆天堂久久| 97精品久久久久久久久久精品| 中国美女看黄片| 亚洲精品第二区| 汤姆久久久久久久影院中文字幕| 建设人人有责人人尽责人人享有的| 只有这里有精品99| 每晚都被弄得嗷嗷叫到高潮| 国产成人一区二区三区免费视频网站 | 精品人妻在线不人妻| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 黑人欧美特级aaaaaa片| 看免费成人av毛片| 亚洲欧美一区二区三区久久| 香蕉丝袜av| 免费日韩欧美在线观看| 日韩av在线免费看完整版不卡| 成年人午夜在线观看视频| 伊人亚洲综合成人网| 久热爱精品视频在线9| 1024视频免费在线观看| 亚洲av综合色区一区| 日本av免费视频播放| 晚上一个人看的免费电影| 国产欧美日韩综合在线一区二区| 精品亚洲乱码少妇综合久久| 少妇的丰满在线观看| 青春草视频在线免费观看| 欧美黄色淫秽网站| 国产欧美日韩一区二区三区在线| 中文字幕人妻熟女乱码| 极品少妇高潮喷水抽搐| 国产成人av教育| 亚洲精品成人av观看孕妇| 亚洲精品自拍成人| 飞空精品影院首页| 国产成人免费观看mmmm| 日韩精品免费视频一区二区三区| 亚洲国产精品一区三区| 亚洲中文av在线| 欧美日韩国产mv在线观看视频| 狠狠婷婷综合久久久久久88av| 中文乱码字字幕精品一区二区三区| 国产成人精品无人区| 97精品久久久久久久久久精品| 国产黄频视频在线观看| 91字幕亚洲| 免费不卡黄色视频| 男女无遮挡免费网站观看| 国产视频首页在线观看| 在线天堂中文资源库| 99久久99久久久精品蜜桃| 久久人妻福利社区极品人妻图片 | 成人黄色视频免费在线看| 亚洲专区中文字幕在线| 狂野欧美激情性xxxx| 日韩欧美一区视频在线观看| 欧美久久黑人一区二区| 成人亚洲精品一区在线观看| 欧美精品啪啪一区二区三区 | 午夜老司机福利片| tube8黄色片| 国产男女内射视频| 亚洲,欧美,日韩| 国产免费视频播放在线视频| 丰满人妻熟妇乱又伦精品不卡| 男女边摸边吃奶| 国产精品一区二区在线不卡| 天天影视国产精品| 日韩制服骚丝袜av| 色婷婷久久久亚洲欧美| 国产精品.久久久| 啦啦啦在线免费观看视频4| 90打野战视频偷拍视频| 欧美国产精品va在线观看不卡| 久久久亚洲精品成人影院| 亚洲av综合色区一区| 久久久精品国产亚洲av高清涩受| 男人添女人高潮全过程视频| 黄片播放在线免费| 成在线人永久免费视频| 日韩 亚洲 欧美在线| 91国产中文字幕| 纯流量卡能插随身wifi吗| 成人国产av品久久久| 久久人人爽人人片av| 麻豆国产av国片精品| 女人高潮潮喷娇喘18禁视频| svipshipincom国产片| 精品人妻熟女毛片av久久网站| 99久久99久久久精品蜜桃| 国产成人精品无人区| 丰满饥渴人妻一区二区三| 亚洲欧洲国产日韩| 亚洲国产看品久久| 亚洲情色 制服丝袜| 久久精品人人爽人人爽视色| 曰老女人黄片| 麻豆国产av国片精品| 国产亚洲av片在线观看秒播厂| 国产成人啪精品午夜网站| 亚洲精品美女久久久久99蜜臀 | 亚洲一区二区三区欧美精品| 18禁观看日本| 熟女少妇亚洲综合色aaa.| 高清视频免费观看一区二区| 女人爽到高潮嗷嗷叫在线视频| 一级片免费观看大全| 国产一级毛片在线| 可以免费在线观看a视频的电影网站| 久热这里只有精品99| 啦啦啦中文免费视频观看日本| 一级毛片黄色毛片免费观看视频| 一边摸一边做爽爽视频免费| 亚洲国产av新网站| av视频免费观看在线观看| 亚洲 欧美一区二区三区| 男女之事视频高清在线观看 | 亚洲av欧美aⅴ国产| 啦啦啦啦在线视频资源| 首页视频小说图片口味搜索 | 亚洲av片天天在线观看| 老鸭窝网址在线观看| 久久久久久久大尺度免费视频| 久久久久视频综合| 国产伦理片在线播放av一区| 精品卡一卡二卡四卡免费| 男女下面插进去视频免费观看| 亚洲av综合色区一区| 久久久久久亚洲精品国产蜜桃av| 亚洲精品一二三| 亚洲成人手机| 麻豆av在线久日| 欧美+亚洲+日韩+国产| 国产男人的电影天堂91| 一二三四社区在线视频社区8| 夜夜骑夜夜射夜夜干| 亚洲精品久久久久久婷婷小说| 国产爽快片一区二区三区| 高清av免费在线| 看免费成人av毛片| 曰老女人黄片| bbb黄色大片| 久久精品aⅴ一区二区三区四区| 一区二区三区四区激情视频| 精品国产乱码久久久久久男人| 国产亚洲av高清不卡| 菩萨蛮人人尽说江南好唐韦庄| 赤兔流量卡办理| 99久久综合免费| 日韩人妻精品一区2区三区| 中文字幕制服av| 美女扒开内裤让男人捅视频| 亚洲中文av在线| 水蜜桃什么品种好| 精品国产一区二区久久| 另类亚洲欧美激情| 9191精品国产免费久久| 日韩大码丰满熟妇| 日本91视频免费播放| 黑人欧美特级aaaaaa片| 肉色欧美久久久久久久蜜桃| 妹子高潮喷水视频| 一本色道久久久久久精品综合| 人妻一区二区av| 亚洲自偷自拍图片 自拍| 欧美精品一区二区免费开放| 天天躁狠狠躁夜夜躁狠狠躁| 日本av手机在线免费观看| 亚洲一码二码三码区别大吗| 亚洲精品国产av蜜桃| 夫妻午夜视频| 色网站视频免费| 亚洲欧美激情在线| 各种免费的搞黄视频| 精品人妻熟女毛片av久久网站| 亚洲精品日本国产第一区| 精品熟女少妇八av免费久了| 午夜免费观看性视频| 丝袜美足系列| 亚洲人成电影观看| 国产精品麻豆人妻色哟哟久久| 久久久久久久精品精品| 亚洲av欧美aⅴ国产| 黑人巨大精品欧美一区二区蜜桃| 国产人伦9x9x在线观看| 狠狠婷婷综合久久久久久88av| 999精品在线视频| 国产97色在线日韩免费| 男人添女人高潮全过程视频| www.精华液| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 精品人妻熟女毛片av久久网站| 99国产精品99久久久久| 天堂中文最新版在线下载| 亚洲五月色婷婷综合| 美女主播在线视频| 丝袜在线中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 日韩中文字幕欧美一区二区 | 国产黄色视频一区二区在线观看| 少妇人妻 视频| 欧美激情 高清一区二区三区| av不卡在线播放| 18禁黄网站禁片午夜丰满| 国产不卡av网站在线观看| 一本大道久久a久久精品| 99精品久久久久人妻精品| 少妇被粗大的猛进出69影院| 亚洲欧美清纯卡通| 欧美日韩福利视频一区二区| 精品人妻在线不人妻| 美国免费a级毛片| 国产一区二区三区av在线| 一区福利在线观看| 欧美乱码精品一区二区三区| 两性夫妻黄色片| 欧美黑人精品巨大| 人妻人人澡人人爽人人| 欧美人与性动交α欧美精品济南到| 国产淫语在线视频| 亚洲国产最新在线播放| 日本色播在线视频| 午夜福利一区二区在线看| 中文字幕精品免费在线观看视频| 99精品久久久久人妻精品| 国产不卡av网站在线观看| 热99国产精品久久久久久7| 亚洲av成人精品一二三区| 人妻一区二区av| 午夜激情av网站| 一本色道久久久久久精品综合| 亚洲中文av在线| 性色av一级| 国产精品 欧美亚洲|