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

    具有弱通訊的多智能體分布式自適應(yīng)協(xié)同跟蹤控制

    2021-01-26 10:48:50薛惠鋒
    關(guān)鍵詞:編隊障礙物通訊

    張 普, 薛惠鋒, 高 山, 左 軒

    (西北工業(yè)大學(xué)自動化學(xué)院, 陜西 西安 710129)

    0 引 言

    近年來,多智能體協(xié)同控制在現(xiàn)代工業(yè)、軍事以及民用領(lǐng)域的廣泛應(yīng)用,如智能機器人、無人機、電網(wǎng)系統(tǒng)、衛(wèi)星姿態(tài)等,吸引了大量學(xué)者的廣泛關(guān)注,如編隊控制[1-3]、協(xié)同追蹤控制[4-6]、協(xié)同容錯控制[7-9]、聚合一致性控制[10]、復(fù)雜動態(tài)網(wǎng)絡(luò)協(xié)同控制[11]、城市交通網(wǎng)絡(luò)擁堵協(xié)同控制[12]、傳感器網(wǎng)絡(luò)估計[13]、復(fù)雜系統(tǒng)協(xié)同控制等[14],其中多智能體協(xié)同跟蹤控制是該領(lǐng)域最典型的熱點研究問題之一。在已有多智能體跟蹤研究中,所采用的模型均具有零輸入狀態(tài)的同構(gòu)跟隨者,而且領(lǐng)航者的位置和速度信息是已知的;此外,所采用的參考軌跡均已規(guī)劃好,即給定已知參考軌跡。然而,在實際應(yīng)用中,通常參考軌跡是按照時間函數(shù)變化的,并通過分布式控制實現(xiàn)協(xié)同跟蹤。

    值得注意的是,與一般單系統(tǒng)的跟蹤控制相比,多智能體分布式協(xié)同跟蹤控制的局限性在于,多智能體編隊系統(tǒng)中只有個別智能體能夠獲得關(guān)于期望軌跡的全部信息,使得彼此之間的信息交互不全,導(dǎo)致跟蹤效果差。針對這一問題,已有文獻采用不同方法進行了相關(guān)研究。比如,基于行為法[15-16]、虛擬結(jié)構(gòu)法[17-18]、領(lǐng)導(dǎo)-跟蹤法[19-21]、基于圖論法[22]以及人工勢場法[23]等。然而,在實際應(yīng)用中,這些方法在解決高階非線性、系統(tǒng)參數(shù)未知或參數(shù)不確定、通信寬帶有限、有限時間收斂、輸入控制有界、網(wǎng)絡(luò)拓撲連通性差、智能體之間數(shù)據(jù)丟失以及自身通訊能力弱等問題時,仍具有一定的局限性。因此,研究多智能體編隊系統(tǒng)分布式自適應(yīng)控制策略,對于含有以上約束的該系統(tǒng)協(xié)同跟蹤研究具有一定的理論意義和工程應(yīng)用價值。

    針對線性多智能體編隊系統(tǒng)在協(xié)同跟蹤中的研究,文獻[24]采用一種分布式自適應(yīng)的方法,研究了帶有未知參數(shù)的高階線性多智能體系統(tǒng),并給出了該系統(tǒng)實現(xiàn)協(xié)同跟蹤的充要條件是網(wǎng)絡(luò)拓撲中大部分智能體之間的交互通訊信號可測量。Ren[25]等人基于低增益反饋理論和李雅普諾夫穩(wěn)定性理論,針對具有切換拓撲結(jié)構(gòu)的二階非線性多智能體系統(tǒng),設(shè)計了一種控制協(xié)同跟蹤控制律,能夠有效地解決該系統(tǒng)在特定約束下的一致性跟蹤問題。文獻[26]針對智能體之間缺乏交互信息而不能保持一致的問題,提出了一種分布式線性協(xié)議方法,設(shè)計了一種協(xié)同跟蹤控制律,使得智能體之間在缺乏信息的情況下,能夠?qū)崿F(xiàn)一致性跟蹤的目的。

    針對非線性多智能體系統(tǒng)的一致跟蹤研究,Yu[27]等人針對含有未知非線性部分的二階多智能體系統(tǒng),采用一種平均脈沖間隔和脈沖增益方法,有效地解決了二階非線性多智能體的耦合編隊跟蹤問題。Wang[28]等人針對追隨者信息未知的二階非線性多智能體系統(tǒng),采用分布式魯棒自適應(yīng)方法,有效地實現(xiàn)了該系統(tǒng)在有限時間的一致跟蹤。文獻[29]針對多智能體在拓撲信息交互動態(tài)時變下的一致性問題,基于圖論提出了多智能體系統(tǒng)實現(xiàn)一致性的條件是智能體之間的有向圖有足夠多的生成樹。文獻[30]針對通訊信號擁堵的二階非線性多智能體系統(tǒng),采用自適應(yīng)量化和控制通訊信號變化的方法,初步實現(xiàn)了對多智能體的輸入輸出量化反饋控制,有效地實現(xiàn)了該系統(tǒng)的一致性跟蹤。文獻[31]針對含有未知參數(shù)及干擾不確定的高階非線性多智能體系統(tǒng),提出一種自適應(yīng)模糊迭代學(xué)習(xí)算法,有效地實現(xiàn)了多智能體在某一區(qū)間內(nèi)能夠平穩(wěn)追蹤目標(biāo)的目的。文獻[32]針對高階非線性多智能體中存在集中信息干擾的問題,采用了一種分布式EOS拓展?fàn)顟B(tài)觀測器的新方法,有效地實現(xiàn)了該系統(tǒng)的一致性跟蹤。

    綜上所述,對于智能體之間信息可測且網(wǎng)絡(luò)通信狀況良好的多智能體系統(tǒng),采用分布式自適應(yīng)方法能夠有效地實現(xiàn)系統(tǒng)的協(xié)同跟蹤控制。然而,對于智能體之間存在信息失真、丟包、延遲、網(wǎng)絡(luò)通信較弱、間歇性通訊等問題的高階非線性多智能體系統(tǒng),目前仍然沒有成熟和統(tǒng)一的控制律設(shè)計規(guī)則,能夠保證智能體實現(xiàn)預(yù)設(shè)性能的協(xié)同跟蹤。

    因此,本文在文獻[27,32]研究的基礎(chǔ)上,針對多智能體編隊系統(tǒng)協(xié)同跟蹤的問題,研究其在弱通訊情況下的協(xié)同跟蹤控制,即跟蹤收斂性和跟蹤效果研究。通過對多智能體系統(tǒng)模型中的不同參數(shù)進行估計,使得智能體在平面內(nèi)整個追蹤過程中,對縱向、橫向以及任意兩個智能體之間的相對距離誤差進行補償控制,以保證其誤差快速收斂;同時引入時變參考軌跡,克服了常規(guī)的給定參考路徑的局限性,有利于工程實踐。相比文獻[27,32]的研究,本文的創(chuàng)新點為:① 基于反推技術(shù)和時變參考軌跡策略,能夠使多智能體在弱通訊情況下實現(xiàn)快速跟蹤,并保持一致性;② 能夠使在多智能體平面跟蹤過程中,其橫向跟蹤誤差、縱向跟蹤誤差、鎮(zhèn)定誤差以及任意兩個智能體之間的相對距離誤差均快速收斂;③ 針對系統(tǒng)中的不確定參數(shù)和非一致動態(tài)的非線性子系統(tǒng),所提方法能夠?qū)ζ溥M行準(zhǔn)確的估計,以便對該系統(tǒng)進行補償控制,實現(xiàn)有效的跟蹤。同時,該方法能夠應(yīng)用于多智能體編隊系統(tǒng)協(xié)同控制,實現(xiàn)以最優(yōu)軌跡避開障礙物,并保持穩(wěn)定跟蹤。

    1 建立數(shù)學(xué)模型

    1.1 圖論知識

    考慮由N個智能體組成的高階非線性系統(tǒng),這些智能體分別標(biāo)記為1,2,…,N。智能體之間的拓撲交互信息用圖論中的G=V×E表示,其中點集V={1,2,…,N}分別代表了N個智能體,邊集E=V×V。每個智能體#i可以與其鄰近的智能體進行信息交換。圖G的鄰接矩陣用A表示,其矩陣中元素aij=1表示智能體#i可以接受智能體#j的信息,反之a(chǎn)ij=0。在網(wǎng)絡(luò)拓撲結(jié)構(gòu)中,度矩陣采用D表示,其中D=diag{d1,d2,…,dN},di=card{j:(i,j)=E}。定義拉普拉斯矩陣為L=D-A。其中,D[dij]∈Rp×p稱為入度矩陣,當(dāng)i≠j時,dij=0;否則,

    1.2 弱通訊知識

    弱通訊是指智能體之間信息交互與共享過程中存在周期性間歇或者時延,會嚴重影響智能體系統(tǒng)的跟蹤或者協(xié)同控制。多智能體系統(tǒng)協(xié)同跟蹤過程中存在弱通訊,如圖1所示。

    圖1 智能體之間信息交互弱通訊示意圖Fig.1 Schematic diagram of weak communication for information interaction among agents

    在圖1中,弱通訊主要分為兩種類型:

    (1) 完全周期性間歇。在有向圖中,智能體與其鄰近的智能體之間的信息傳遞會出現(xiàn)周期性中斷,則稱之為完全周期性間歇。

    (2) 部分周期性間歇。在有向圖中相鄰僚機之間保持信息交互,但僚機與長機之間信息交互會出現(xiàn)周期性中斷,則稱之為部分周期性間歇。

    引理 1對于多智能體系統(tǒng),如果滿足

    (1)

    則其能夠?qū)崿F(xiàn)信息的一致性。其中,x和v分別表示智能體的位置和速度,(i,j)表示兩個不同的智能體?;谑?1),t同時滿足圖1弱通訊傳輸?shù)臈l件,這也是本文基于此設(shè)計控制律的依據(jù)。

    引理 2在圖1中,對于T>0,在間隔[mT,(m+1)T]內(nèi),定義信息的傳輸速率為ρ=δ/T。其中δ表示在有信號的情形下,傳輸過程所需要的時間,T表示在弱通訊情形下信號傳輸?shù)闹芷?包括有信號傳輸和無信號傳輸,其表達式為

    (2)

    式中,m∈Ν;k表示任意兩個智能體在信息傳遞過程的權(quán)重因子,其大小取決于其拓撲結(jié)構(gòu)中的拉普拉斯矩陣L;z表示跟蹤誤差;ui(t)表示第i個智能體的控制輸入。

    1.3 高階非線性模型

    針對含有弱通訊的多智能體系統(tǒng),建立其非線性模型:

    (3)

    Q(u)=p(u)u+Δ(u)

    (4)

    式中,Δ(u)表示任意兩個智能體在弱通訊情形下的控制信號誤差。在多智能體跟蹤過程中,視為一種干擾項。本文所設(shè)計的控制律是為了消除這種擾動對追蹤效果的影響。除此之外,p(u)的值取決于輸入信號的死區(qū)范圍:

    (5)

    式中,umin和umax分別表示控制信號u的死區(qū)上、下邊界。

    在任意兩個智能體之間通訊正常,但通訊信號強度較弱情況下,其信號誤差滿足

    Δ(u)≤λ|u|

    (6)

    式中,λ表示實常數(shù)。

    1.4 領(lǐng)航者動力學(xué)模型

    在慣性坐標(biāo)下,領(lǐng)航者的導(dǎo)航模型為

    (7)

    式(7)對時間求導(dǎo),可得

    (8)

    式中,(ax,ay,az)表示空間3個方向的加速度。

    根據(jù)式(8),可得領(lǐng)航者的動力學(xué)方程:

    (9)

    式中,(Fx,Fy,Fz)表示空間3個方向的力;m表示領(lǐng)航智能體的質(zhì)量。

    引理 3存在正實數(shù)a,b,p,q,滿足1/p+1/q=1,則ab≤1/pap+1/qbq, 等號成立當(dāng)且僅當(dāng)ap=bq,此為楊氏不等式理論。

    假設(shè) 1多智能體編隊系統(tǒng)的通信拓撲圖G是固定且有向的,確保多智能體編隊系統(tǒng)在追蹤過程中保持隊形恒定。

    在追蹤控制律設(shè)計過程中,為了消除非線部分Δ(u),根據(jù)引理3,有

    (10)

    于是,可得

    則式(10)成立。

    假設(shè) 2多智能體之間的通訊是正常且順暢,僅信號的強度不同。

    假設(shè) 3參考軌跡yr是光滑有界的,且對于整個編隊系統(tǒng)均已知,以此確保多智能體編隊系統(tǒng)均按照時變的參考軌跡跟蹤期望的位置信號,然后按照預(yù)設(shè)的隊形運動。

    1.5 避障模型

    為避免在此跟蹤過程中多智能體與障礙物之間發(fā)生碰撞,引入人工勢場法建立多智能體避障模型,即智能體和障礙物帶有同種電荷時相互排斥,帶異種電荷時相互吸引。多智能體中任意一個智能體所受的合力為

    fc i(t)=fi a(t)+fi r(t)+h(t)

    (11)

    式中,fc i表示第i個智能體所受的合力;fi a表示第i個智能體所受的引力;fi r(t)表示第i個智能體所受的斥力;h(t)表示調(diào)節(jié)函數(shù)。由于多智能體在追蹤過程中,易出現(xiàn)局部極小值點,為了打破這一平衡狀態(tài),設(shè)計函數(shù)

    (12)

    式中,sign(ε)表示自變量為ε的符號函數(shù),ε為很小的正常數(shù)。

    除此之外,每個智能體所受的合力為

    (13)

    式中,mi和mj分別表示第i和j個智能體的質(zhì)量;rij表示第i個智能體與第j個智能體之間的相對距離;bij表示拉普拉斯矩陣L中第ij項的數(shù)值。

    當(dāng)智能體靠近障礙物時,其相對距離變小,智能體所受的斥力大于引力,合力表現(xiàn)為斥力,智能體會遠離障礙物;當(dāng)智能體與障礙物遠離時,其相對距離減小,智能體所受的引力大于斥力,合力表現(xiàn)為引力,智能體會靠近障礙物;當(dāng)智能體與障礙物的相對距離為特定值時,斥力等于引力,此時會出現(xiàn)局部最小值點,通過函數(shù)h(t)微調(diào)其合力,表現(xiàn)為斥力或者引力,使智能體打破平衡狀態(tài)。同時,多智能體編隊系統(tǒng)以時變的參考估計跟蹤目標(biāo),因此多智能體會以最優(yōu)路徑避開障礙物,跟蹤運動目標(biāo),并保持隊形不變。

    2 分布式追蹤控制律設(shè)計

    在本文中,多智能體系統(tǒng)由N個智能體構(gòu)成“領(lǐng)航者-跟隨者”分布式模式的編隊結(jié)構(gòu),其中包含1個領(lǐng)航者,N-1個追隨者?;谙到y(tǒng)模型式(3)、領(lǐng)航者動力學(xué)模型式(9)和避障模型式(11),針對弱通訊中完全周期性間歇的多智能體系統(tǒng),采用自適應(yīng)反演控制方法進行協(xié)同跟蹤控制律設(shè)計,以避開障礙物實現(xiàn)協(xié)同跟蹤。

    定義智能體的追蹤誤差為

    ki(yi(t)-yr(t))

    (14)

    式中,i=1,2,…,n,aij表示智能體#i和智能體#j之間信息傳遞的權(quán)重系數(shù);ki表示正常系數(shù);yi-yj表示智能體#i和智能體#j之間的相對距離差;yr為參考軌跡。

    定義新誤差為

    zi,2(t)=xi,2(t)-αi,1(t)

    (15)

    式(14)對時間求導(dǎo),并將式(15)代入求導(dǎo)結(jié)果,整理化簡得

    (16)

    構(gòu)造李雅普諾夫函數(shù):

    (17)

    式(17)對時間求導(dǎo),并將式(14)代入求導(dǎo)結(jié)果,整理可得

    (18)

    由式(18)可得虛擬控制律和調(diào)節(jié)函數(shù),分別表示為

    (19)

    (20)

    將式(19)和式(20)代入式(18),可得

    (21)

    同時,將式(19)代入式(16),整理化簡可得

    (22)

    定義新誤差為

    zi,3(t)=xi,3(t)-αi,2(t)

    (23)

    式(15)對時間求導(dǎo),并結(jié)合式(23),可得

    (24)

    式(19)對時間求導(dǎo),可得

    (25)

    將式(25)代入式(24),可得

    (26)

    構(gòu)造李雅普諾夫函數(shù)

    (27)

    式(27)對時間求導(dǎo),可得

    (28)

    將式(21)和式(26)代入式(28),可得

    (29)

    同理,由式(29)可得調(diào)節(jié)函數(shù)和虛擬控制律,分別表示為

    (31)

    將式(30)和式(31)代入式(29),可得

    (32)

    將式(31)代入式(26),可得

    (33)

    定義新誤差為

    zi,q+1(t)=xi,q+1(t)-αi,q(t)

    (34)

    (35)

    構(gòu)造李雅普諾夫函數(shù):

    (36)

    式中,Vi,q -1是第q-1步迭代設(shè)計過程中的李雅普諾夫函數(shù)。

    (37)

    由式(37)可得調(diào)節(jié)函數(shù)和虛擬控制律,分別表示為

    (38)

    αi,q(t)=-zi,q -1(t)-cqzi,q(t)-

    (39)

    定義誤差為

    zi,n(t)=xi,n(t)-αi,n -1(t)

    (40)

    (41)

    構(gòu)造李雅普諾夫函數(shù):

    (42)

    式中,Vi,n -1(t)是第n-1步迭代設(shè)計過程中的李雅普諾夫函數(shù)。

    (43)

    根據(jù)式(3)和式(4),結(jié)合式(43),可得調(diào)節(jié)函數(shù)和最終控制律分別為

    (44)

    (45)

    因此,式(45)表示多智能體編隊系統(tǒng)處于弱通訊環(huán)境中的協(xié)同跟蹤控制律,確保系統(tǒng)穩(wěn)定。而多智能體編隊系統(tǒng)的模型為高階非線性系統(tǒng),于是所設(shè)計的控制律為n階非線性系統(tǒng)的追蹤控制律,包含給定的常系數(shù)、更新律以及跟蹤誤差。對于仿真驗證而言,給定二階系統(tǒng),即n=2,結(jié)合領(lǐng)航者動力學(xué)模型、初始設(shè)置的參數(shù),就可以得出具體的控制律:

    (46)

    將式(44)和式(45)代入式(43),可得

    (47)

    3 穩(wěn)定性分析

    針對系統(tǒng)模型式(3),所設(shè)計的分布式自適應(yīng)跟蹤控制律和參數(shù)自適應(yīng)更新律,能夠使多智能體編隊系統(tǒng)跟蹤期望的指定位置,同時保持一致。在滿足假設(shè)1和假設(shè)2的條件下,為了驗證所設(shè)計的控制律在多智能體系統(tǒng)跟蹤和保持協(xié)同過程中系統(tǒng)的穩(wěn)定性,構(gòu)造整個多智能體編隊系統(tǒng)的二次型函數(shù)為

    VN(t)=Vi,1(t)+Vi,2(t)+Vi,3(t)+…+Vi,n(t)=

    (48)

    式(48)對時間求導(dǎo),可得

    (49)

    (50)

    式中,ck>0(k=1,2,…,n)。

    圖2 障礙物、智能體以及目標(biāo)之間的受力關(guān)系圖Fig.2 Diagram of the force relationship among obstacle, agent and target

    (51)

    式中,多智能體編隊系統(tǒng)所受的合力為0,則多智能體之間的隊形保持不變,而且編隊系統(tǒng)避開障礙物,同時在完全周期性間歇情形下實現(xiàn)跟蹤的目的。在跟蹤和隊形保持過程中,智能體之間需要實時進行信息的傳遞,而在受限狀態(tài)(弱通訊)下,即信號時強時弱,導(dǎo)致通信受阻難以完成跟蹤與隊形保持任務(wù),而且極易發(fā)生碰撞。因此,本文設(shè)計了擾動補償控制律,能夠使多智能體編隊系統(tǒng)在受限狀態(tài)下依舊保持跟蹤效果和恒定隊形。

    4 仿真實驗

    在本節(jié)中,采用5個智能體作為被控對象構(gòu)成多智能體編隊系統(tǒng),其中包含1個領(lǐng)航者(編號為5)和4個跟隨者(編號為1, 2, 3, 4)。每個智能體結(jié)構(gòu)和尺寸都相等,即是同構(gòu)的。同時,每個智能體具有相同的運動特性,包含相同的偏航角和相對位置;這5個智能體在平面內(nèi)運動,且其通訊是順暢的,拓撲結(jié)構(gòu)示意圖如圖3所示。

    圖3 多智能體系統(tǒng)拓撲結(jié)構(gòu)示意圖Fig.3 Structure schematic diagram of the multi-agent system topology

    表1 5個智能體的初始位置,期望指定位置以及參數(shù)估計值

    由第1.1節(jié)圖論可知,鄰接矩陣A,度矩陣D以及拉普拉斯矩陣L,其關(guān)系為

    L=D-A

    (52)

    同時,結(jié)合第1.1節(jié),圖3以及式(52),可得

    4.1 收斂性分析

    為了驗證在弱通訊情況下所提方法在多智能體編隊追蹤過程中的收斂性,基于多智能體的系統(tǒng)模型和領(lǐng)航者的動力學(xué)模型,結(jié)合自適應(yīng)反推技術(shù),采用5個智能體進行仿真實驗,其結(jié)果如圖4~圖6所示。

    圖4 多智能體系統(tǒng)前行速度誤差曲線Fig.4 Forward speed error curve of the multi-agent system

    圖5 參數(shù)誤差曲線Fig.5 Curve of parameter error

    圖6 鎮(zhèn)定誤差曲線Fig.6 Curve of stabilization error

    圖4給出了多智能體系統(tǒng)前行速度誤差曲線。由圖4可知,隨著時間的變化,曲線大致呈現(xiàn)先緩慢減小后趨于穩(wěn)定的狀態(tài)。這是由于每個智能體從不同的位置出發(fā),各自以時變的參考軌跡運動,跟蹤各自期望的指定位置。在每個智能體跟蹤上各自期望的指定位置前,每個智能體與期望指定位置之間的相對速度逐漸縮小,相應(yīng)智能體前行速度誤差也相應(yīng)減小,當(dāng)智能體跟蹤上期望的指定位置時,此時指定期望位置也在運動,但是智能體的實際位置與指定的期望位置之間保持恒定的距離。因此,智能體的前行速度誤差基本接近0,并保持這種速度運動,協(xié)同達到一致運動的目的。在10 s之前,智能體#3和智能體#4曲線下降速度快,而智能體#5緩慢上升,這是由于其起點位置不同,而且初速度也不同,這使得每個智能體與期望的指定位置存在速度差。速度差最大為智能體#4,速度差最小的是智能體#5。在10 s之后,每個智能體前行速度誤差基本趨于0,也就是每個智能體已跟蹤上各自的期望指定位置,并保持一致運動。

    圖6給出了系統(tǒng)鎮(zhèn)定誤差曲線圖。隨著時間的變化,該曲線大致呈現(xiàn)先緩慢下降后趨于穩(wěn)定狀態(tài)的趨勢。在5個智能體跟蹤期望的指定位置過程中,所提方法能夠?qū)崿F(xiàn)快速有效的跟蹤,并使得5個智能體達到一致性。本文所采用的核心技術(shù)是反推技術(shù),而反推技術(shù)在求解控制律過程,每一步均會設(shè)定鎮(zhèn)定誤差,逐步遞推,然后得到最終的控制律。由圖6可知,智能體#1~智能體#3的鎮(zhèn)定誤差曲線逐漸緩慢下降,這說明其在跟蹤期望的指定位置的過程先于智能體#4和智能體#5。智能體#4以非常緩慢的趨勢逼近期望指定位置,鎮(zhèn)定誤差變化很小,說明該智能體在跟蹤期望的指定位置過程中很平穩(wěn)。然而,智能體#5曲線變化比較大,先是瞬間下降后緩慢上升至穩(wěn)定狀態(tài),因為智能體#5作為領(lǐng)航者,其速度比智能體#4稍微小點,但比其他智能體均大,在跟蹤過程中需要保持隊形的一致性,因此該曲線變化趨勢明顯。

    4.2 跟蹤效果分析

    為了進一步驗證所設(shè)計分布式自適應(yīng)一致跟蹤方法的有效性,采用同樣的仿真環(huán)境和相同的初始條件,對5個智能體編隊系統(tǒng)進行跟蹤效果分析,如圖7~圖10所示。

    圖7 多智能體系統(tǒng)在x方向的跟蹤誤差曲線Fig.7 Tracking error curve of multi-agent system in the x direction

    圖8 多智能體系統(tǒng)在y方向的追蹤誤差曲線Fig.8 Tracking error curve of the multi-agent system in the y direction

    圖9 多智能體系統(tǒng)在x方向的相對距離誤差曲線Fig.9 Relative distance exerror curve of the multi-agent system in the x direction

    圖10 多智能體系統(tǒng)在y方向的相對距離誤差曲線Fig.10 Relative distance error curve of the multi-agent system in the y direction

    圖7給出了多智能體系統(tǒng)在x方向的追蹤誤差曲線。隨著時間的變化,該曲線呈現(xiàn)正弦波的趨勢。5個智能體編隊系統(tǒng)從不同起點以時變的參考軌跡運動跟蹤期望的指定位置并保持一致運動,在此過程中分為跟蹤階段和保持階段。因此,在6 s之前,該曲線處于分離狀態(tài),這說明5個智能體以不同的初速度、偏航角以及控制參數(shù)進行運動,跟蹤指定的期望位置;通過調(diào)節(jié)系統(tǒng)的參數(shù)θ,進而實現(xiàn)不同的虛擬控制,使得每個智能體均能按照時變的參考軌跡跟蹤指定的期望位置。在6 s之后,該曲線呈現(xiàn)重合后正弦波的趨勢。這說明5個智能體在前進方向?qū)崿F(xiàn)了跟蹤后保持一致運動,采用的自適應(yīng)協(xié)同跟蹤算法可以實時控制每個智能體在x方向的速度、偏航角以及系統(tǒng)的參數(shù),實現(xiàn)對多智能體系統(tǒng)協(xié)同跟蹤目的。

    圖8給出了多智能體系統(tǒng)在y方向的跟蹤誤差曲線。隨著時間的變化,該曲線呈現(xiàn)從分離狀態(tài)到重合狀態(tài)的波動趨勢。同樣地,5個智能體編隊系統(tǒng)在平面中以不同的速度和偏航角開始運動,然后按照時變的參考軌跡運動,追蹤各自期望的指定位置,形成模型中給出的五邊形隊形,并保持該隊形繼續(xù)運動。在此過程中,同樣分為跟蹤和保持一致性兩個階段。在10 s之前,5個智能體從不同起點運動,由于參考的軌跡是正弦變化的路徑,在y方向有相應(yīng)的運動,5個智能體的追蹤誤差曲線有較大的變化,這是由于每個智能體在追蹤過程中出現(xiàn)通訊能力弱的情況,基于反推技術(shù)的自適應(yīng)控制算法能夠使得每個智能體較快調(diào)整各自在y方向的速度和偏航角,使其盡快按照時變參考軌跡追蹤期望的指定位置,實現(xiàn)跟蹤目的。在10 s之后,5個智能體已跟蹤到期望的指定位置并形成模型中的五邊形隊形,然后保持該隊形運動,其在y方向的路徑跟蹤誤差為0,實現(xiàn)了較好的隊形保持目的。

    圖9給出了多智能體系統(tǒng)在x方向的相對距離誤差曲線。隨著時間的變化,該曲線呈現(xiàn)先快速下降,然后緩慢下降至穩(wěn)定狀態(tài)的趨勢。5個智能體從不同起點出發(fā),以不同速度和偏航角跟蹤期望的指定位置,由于5個智能體的初速度均不相等,尤其領(lǐng)航者(智能體#5)速度較大,智能體#1~智能體#3的速度較小,智能體#4的速度最大,導(dǎo)致5個智能體跟蹤過程出現(xiàn)不等的相對距離差,智能體#2~智能體#4出現(xiàn)快速下降,而智能體#1出現(xiàn)較為緩慢的上升,之后5個智能體均跟蹤已確定的期望指定位置直至達到穩(wěn)定狀態(tài),這是由于基于反推技術(shù)的自適應(yīng)控制方法使得以參考軌跡快速跟蹤期望的指定位置,之后保持一致運動。因此,在10 s之前,5條曲線依次呈現(xiàn)不同的下降和上升的趨勢,直至35 s之后,曲線保持穩(wěn)定,實現(xiàn)5個智能體編隊系統(tǒng)跟蹤和保持一致性的目的。

    圖10給出了多智能體系統(tǒng)在y方向的相對距離誤差曲線。隨著時間的變化,該曲線呈現(xiàn)快速下降后快速達到穩(wěn)定的趨勢。5個智能體編隊系統(tǒng)開始從不同位置以不同速度按照時變的參考軌跡運動,5個智能體在完成編隊跟蹤期望的指定位置和保持五邊形的隊形之余,還需要使得5個智能體中任意兩個智能體在通訊能力差的情況下保持之間的相對距離大于安全距離,實現(xiàn)在跟蹤過程中不發(fā)生碰撞。因此,在y方向曲線會呈現(xiàn)不同速率的下降趨勢,而后在自適應(yīng)控制算法的控制下,實現(xiàn)穩(wěn)定狀態(tài),也就是曲線呈現(xiàn)水平趨勢。在10 s前,智能體#3的速度相對比較大(V3=1.0 m/s),快速跟蹤期望的指定位置,在控制參考的調(diào)整下,曲線變化趨勢大;智能體#5作為領(lǐng)航者,其速度較大(V5=1.2 m/s),在跟蹤過程保持較為平穩(wěn)的追蹤,曲線變化趨勢平緩。在10 s之后,5個智能體在y方向相對距離誤差均為0,實現(xiàn)一致跟蹤的目的。

    4.3 避障分析

    在本節(jié)中,基于避障模型式(11),采用相同的初始條件和相同的實驗環(huán)境,仿真結(jié)果如圖11和圖12所示。

    圖11 多智能體避障曲線Fig.11 Collision avoidance curve of the multi-agent

    圖12 多智能體避障過程相對距離誤差曲線圖Fig.12 Relative distance error curve of the multi-agent collision avoidance process

    圖11給出了多智能體在平面運動的避障圖。由圖11可知,隨著時間的變化,該曲線呈現(xiàn)以最優(yōu)軌跡避開障礙物波動的趨勢。在多智能體跟蹤過程中,每個智能體由不同的起點位置以時變的參考軌跡跟蹤期望的指定位置,然后保持既定的編隊隊形運動。在此過程中,會遇到障礙物,本文采用人工勢場方法進行避障,其核心為智能體與障礙物之間的距離決定彼此遠離或者靠近,從而達到避障目標(biāo)。在5 s之前,智能體逐漸靠近障礙物,其相對距離減小,產(chǎn)生的斥力大于引力,使得多智能體稍微遠離障礙物而進行運動。同理,當(dāng)智能體所示斥力小于引力,則相互之間靠近,呈現(xiàn)一種以最優(yōu)軌跡避開障礙物運動。

    圖12給出了多智能體避障過程相對距離誤差曲線。由圖12可知,隨著時間變化,該曲線首先快速上升,然后快速下降到很小,之后緩慢趨于穩(wěn)定狀態(tài)。這是由于多智能體編隊系統(tǒng)在剛開始逐漸靠近障礙物時,其合力為斥力,同時還繼續(xù)跟蹤期望的指定位置,從而使得其相對距離誤差變大;隨后,當(dāng)其距離逐漸增大,合力表現(xiàn)為引力;最后,智能體與障礙物之間的引力和斥力等大反向,達到動態(tài)平衡,其曲線處于水平狀態(tài)。

    4.4 討論

    為了驗證本文所提出算法的追蹤效果,將與滑模控制(sliding mode control, SMC)方法,自適應(yīng)控制方法(adaptive control method, ACM)進行仿真對比分析,采用相同的初始條件和假設(shè),其結(jié)果如圖13所示。

    圖13 多智能體跟蹤誤差曲線Fig.13 Tracking error curve of the multi-agent

    圖13給出了多智能體跟蹤誤差曲線。在所提方法、ACM以及SMC方法的作用下,隨著時間的變化,曲線先快速上升到最高峰值,然后迅速下降,最后緩慢趨于穩(wěn)定狀態(tài)。這是由于智能體在不同的初始位置以不同的角度開始運動,根據(jù)參考軌跡跟蹤期望的指定位置,在峰值之前,曲線迅速上升,這是由于每個智能體以時變的軌跡快速跟蹤期望的指定位置,智能體與期望的指定位置之間的距離相距較遠,這使得其跟蹤誤差較大。當(dāng)該曲線從峰值到水平位置時,3種方法均迅速下降后緩慢上升至水平位置,這是由于智能體即將接近各自的期望指定位置,速度減小,對應(yīng)的相對距離誤差也隨之減小至0,則實現(xiàn)一致跟蹤的目的。然后,所提方法使得多智能體系統(tǒng)的跟蹤誤差最小,達到穩(wěn)態(tài)所需要的時間最短,ACM對應(yīng)的跟蹤效果次之,SMC方法跟蹤效果最差。

    為了進一步直觀說明3種方法的跟蹤效果,對其追蹤過程中的穩(wěn)態(tài)誤差和穩(wěn)定時間進行對比分析,如表2所示。

    表2 3種方法跟蹤效果指標(biāo)對比

    由表2可知,本文所提方法能夠使多智能體在弱通訊情形下具有良好的追蹤效果。

    5 結(jié) 論

    本文基于反推技術(shù),提出分布式自適應(yīng)一致跟蹤控制方法,達到多智能體在弱通訊的情況能夠?qū)崿F(xiàn)有效的跟蹤和保持一致運動的目的。5個智能體作為被控對象構(gòu)成多智能體系統(tǒng),這5個智能體均攜帶有無線傳輸系統(tǒng),在復(fù)雜環(huán)境下,信號傳輸能力存在弱化情況。基于此,本文將5個智能體系統(tǒng)模型分解為5個子系統(tǒng)進行虛擬控制律設(shè)計,避免非線性系統(tǒng)的處理所帶來的困難,從而實現(xiàn)最終的控制律設(shè)計,實現(xiàn)對5個智能體系統(tǒng)進行有效控制。主要結(jié)論和存在問題如下:

    (1) 所提出的多智能體分布式自適應(yīng)方法能夠使5個智能體在平面x和y方向進行有效追蹤,同時保持鎮(zhèn)定誤差收斂于很小的領(lǐng)域內(nèi),使多智能體編隊系統(tǒng)在弱通訊情形下完成跟蹤任務(wù),并以最優(yōu)軌跡避開障礙物并保持隊形恒定;

    (2) 基于反推技術(shù),提出自適應(yīng)一致跟蹤控制律,能夠使多智能體在整個跟蹤過程具有較好的穩(wěn)定性,誤差收斂速度快;

    (3) 盡管所提出的方法能夠使多智能體編隊系統(tǒng)在弱通訊情形下具有較好的追蹤效果,但本文是將弱通訊情況作為系統(tǒng)的擾動進行處理的,沒有分析其內(nèi)部的通訊機理,也未能建立精確的多智能體模型。因此,下一步工作將對弱通訊情況下的傳輸機理進行深入分析,建立準(zhǔn)確的多智能體模型,進行跟蹤控制。

    猜你喜歡
    編隊障礙物通訊
    2023年1月25日,美軍一次演習(xí)期間,空軍正在進行編隊飛行
    軍事文摘(2023年5期)2023-03-27 08:56:26
    《茶葉通訊》簡介
    茶葉通訊(2022年2期)2022-11-15 08:53:56
    《茶葉通訊》簡介
    茶葉通訊(2022年3期)2022-11-11 08:43:50
    通訊報道
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計和處理
    基于事件驅(qū)動的多飛行器編隊協(xié)同控制
    通訊簡史
    基于預(yù)測控制的無人機編隊內(nèi)部避碰
    多彈編隊飛行控制技術(shù)研究
    国产精品99久久99久久久不卡| 成人特级av手机在线观看| 国产伦精品一区二区三区四那| 国产综合懂色| 十八禁网站免费在线| 国产精品一区二区三区四区免费观看 | 精品日产1卡2卡| 久久精品综合一区二区三区| h日本视频在线播放| 亚洲激情在线av| 淫妇啪啪啪对白视频| 最近最新免费中文字幕在线| 国产野战对白在线观看| 小说图片视频综合网站| 亚洲男人的天堂狠狠| 国产精品亚洲美女久久久| 久99久视频精品免费| 国产国拍精品亚洲av在线观看 | 男女下面进入的视频免费午夜| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站 | 日韩欧美 国产精品| 在线播放国产精品三级| 97人妻精品一区二区三区麻豆| 黄色日韩在线| 特大巨黑吊av在线直播| 高清毛片免费观看视频网站| 欧美一区二区精品小视频在线| 窝窝影院91人妻| 久久九九热精品免费| 有码 亚洲区| 美女高潮喷水抽搐中文字幕| 91久久精品电影网| 很黄的视频免费| 99精品欧美一区二区三区四区| 好看av亚洲va欧美ⅴa在| 亚洲成a人片在线一区二区| 欧美在线一区亚洲| 国产欧美日韩精品亚洲av| 久久久久亚洲av毛片大全| 99久久精品一区二区三区| 亚洲成人久久爱视频| 亚洲熟妇中文字幕五十中出| 悠悠久久av| 精品人妻一区二区三区麻豆 | 天堂网av新在线| 一边摸一边抽搐一进一小说| 99热只有精品国产| 亚洲专区国产一区二区| 亚洲久久久久久中文字幕| 成人午夜高清在线视频| 婷婷精品国产亚洲av在线| 精品久久久久久久人妻蜜臀av| 亚洲国产色片| 噜噜噜噜噜久久久久久91| 日本与韩国留学比较| 在线a可以看的网站| 亚洲av成人精品一区久久| 熟女人妻精品中文字幕| 男女之事视频高清在线观看| 日本 av在线| 国产精品亚洲一级av第二区| 国产高清videossex| 给我免费播放毛片高清在线观看| 97超视频在线观看视频| 神马国产精品三级电影在线观看| 亚洲国产精品合色在线| svipshipincom国产片| 一二三四社区在线视频社区8| 成人三级黄色视频| 9191精品国产免费久久| 国产中年淑女户外野战色| 舔av片在线| 久久精品夜夜夜夜夜久久蜜豆| 小蜜桃在线观看免费完整版高清| 97超级碰碰碰精品色视频在线观看| 欧美性猛交╳xxx乱大交人| avwww免费| 国内精品一区二区在线观看| 国产成人系列免费观看| 一个人看视频在线观看www免费 | bbb黄色大片| 亚洲性夜色夜夜综合| 午夜福利高清视频| 午夜福利成人在线免费观看| 小蜜桃在线观看免费完整版高清| 午夜精品在线福利| 欧美性猛交╳xxx乱大交人| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美免费精品| 午夜免费成人在线视频| 搞女人的毛片| 99久国产av精品| 哪里可以看免费的av片| 又黄又粗又硬又大视频| 51国产日韩欧美| 欧美成人性av电影在线观看| 国产真人三级小视频在线观看| 成人国产综合亚洲| 成人18禁在线播放| 99精品久久久久人妻精品| 19禁男女啪啪无遮挡网站| 嫩草影视91久久| 成人特级黄色片久久久久久久| 最新美女视频免费是黄的| 非洲黑人性xxxx精品又粗又长| 日本a在线网址| 淫秽高清视频在线观看| 成人特级av手机在线观看| 热99在线观看视频| 欧美一级a爱片免费观看看| 亚洲av熟女| 亚洲美女黄片视频| 一本综合久久免费| 欧美中文综合在线视频| 久久性视频一级片| 女生性感内裤真人,穿戴方法视频| 国产三级在线视频| 欧美日韩综合久久久久久 | 一a级毛片在线观看| 人人妻人人看人人澡| 中出人妻视频一区二区| 国内揄拍国产精品人妻在线| 精品国产超薄肉色丝袜足j| 波多野结衣巨乳人妻| 91麻豆av在线| 国模一区二区三区四区视频| 亚洲乱码一区二区免费版| 免费av不卡在线播放| 亚洲成av人片在线播放无| 日本三级黄在线观看| 午夜福利视频1000在线观看| 国产一区二区三区在线臀色熟女| 在线观看一区二区三区| 熟女电影av网| 高潮久久久久久久久久久不卡| 欧美日韩亚洲国产一区二区在线观看| 天天躁日日操中文字幕| 亚洲av电影在线进入| 亚洲精品成人久久久久久| 国产亚洲精品久久久久久毛片| 国产午夜精品论理片| 国产成人啪精品午夜网站| 十八禁人妻一区二区| 听说在线观看完整版免费高清| 俄罗斯特黄特色一大片| 日本黄大片高清| 色综合婷婷激情| 欧美日韩亚洲国产一区二区在线观看| 久久久国产成人免费| 色精品久久人妻99蜜桃| 一夜夜www| 波多野结衣高清作品| 日韩精品青青久久久久久| 免费av观看视频| 看免费av毛片| 黄色丝袜av网址大全| 亚洲成a人片在线一区二区| 欧美色视频一区免费| 狠狠狠狠99中文字幕| 九九热线精品视视频播放| 久久久久久久久久黄片| 久久亚洲精品不卡| 亚洲 国产 在线| 国产一级毛片七仙女欲春2| 亚洲午夜理论影院| 男人舔奶头视频| 一本综合久久免费| 小蜜桃在线观看免费完整版高清| 亚洲精品乱码久久久v下载方式 | 波多野结衣巨乳人妻| 午夜视频国产福利| 欧美最黄视频在线播放免费| 欧美日韩亚洲国产一区二区在线观看| 黄色片一级片一级黄色片| 在线观看免费午夜福利视频| 嫁个100分男人电影在线观看| 亚洲av免费在线观看| 男女床上黄色一级片免费看| 国产精品一及| 在线视频色国产色| 国产高清视频在线播放一区| 国产精华一区二区三区| 日本免费一区二区三区高清不卡| 搡老熟女国产l中国老女人| 日本三级黄在线观看| 久久精品国产99精品国产亚洲性色| 午夜福利免费观看在线| 一二三四社区在线视频社区8| 久久久久久大精品| 九色成人免费人妻av| 久久久久久久久大av| 又粗又爽又猛毛片免费看| 欧美日韩黄片免| 全区人妻精品视频| 精品电影一区二区在线| 搡老岳熟女国产| 天堂动漫精品| 啪啪无遮挡十八禁网站| 久久久国产成人免费| 成年女人毛片免费观看观看9| av黄色大香蕉| 亚洲国产日韩欧美精品在线观看 | 欧美bdsm另类| 高清日韩中文字幕在线| 五月伊人婷婷丁香| 国内揄拍国产精品人妻在线| 日韩欧美国产一区二区入口| 国产极品精品免费视频能看的| 亚洲中文字幕日韩| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 中文亚洲av片在线观看爽| 国产激情偷乱视频一区二区| 亚洲国产精品久久男人天堂| 久久草成人影院| 成年免费大片在线观看| 午夜福利在线观看吧| 制服丝袜大香蕉在线| 国产欧美日韩一区二区三| 精华霜和精华液先用哪个| 动漫黄色视频在线观看| 午夜日韩欧美国产| 国内精品久久久久久久电影| 男女视频在线观看网站免费| 18美女黄网站色大片免费观看| 听说在线观看完整版免费高清| 淫秽高清视频在线观看| 国产精品自产拍在线观看55亚洲| 香蕉av资源在线| 搡女人真爽免费视频火全软件 | 最近视频中文字幕2019在线8| 老司机在亚洲福利影院| 久久精品国产亚洲av涩爱 | 高潮久久久久久久久久久不卡| 在线播放国产精品三级| 国产aⅴ精品一区二区三区波| 国产又黄又爽又无遮挡在线| 欧美中文综合在线视频| 亚洲第一电影网av| 国产色婷婷99| 国产高清三级在线| 淫妇啪啪啪对白视频| www国产在线视频色| 日本撒尿小便嘘嘘汇集6| 久99久视频精品免费| 午夜精品在线福利| 久久精品综合一区二区三区| 欧美不卡视频在线免费观看| 无限看片的www在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 韩国av一区二区三区四区| 午夜免费男女啪啪视频观看 | 1000部很黄的大片| 国产高清视频在线播放一区| 国产三级中文精品| 中文字幕av在线有码专区| 欧美大码av| 最新美女视频免费是黄的| 成人18禁在线播放| 男女那种视频在线观看| 99精品在免费线老司机午夜| 欧美一级a爱片免费观看看| 怎么达到女性高潮| 人妻丰满熟妇av一区二区三区| 亚洲在线观看片| 成人欧美大片| 午夜日韩欧美国产| 中文资源天堂在线| 一区二区三区国产精品乱码| a级一级毛片免费在线观看| 老司机在亚洲福利影院| 国产精品亚洲美女久久久| 天天一区二区日本电影三级| 成人国产综合亚洲| 午夜免费男女啪啪视频观看 | 亚洲男人的天堂狠狠| 黄色丝袜av网址大全| 久久久久精品国产欧美久久久| 真实男女啪啪啪动态图| 最近最新中文字幕大全电影3| 日本三级黄在线观看| 级片在线观看| 国产在线精品亚洲第一网站| 国产免费男女视频| 激情在线观看视频在线高清| 看黄色毛片网站| 最近最新免费中文字幕在线| 国产aⅴ精品一区二区三区波| 脱女人内裤的视频| 日日夜夜操网爽| 日本一本二区三区精品| 最新美女视频免费是黄的| 九色国产91popny在线| 很黄的视频免费| 日韩欧美免费精品| 88av欧美| 97人妻精品一区二区三区麻豆| 亚洲精品成人久久久久久| 国产精品98久久久久久宅男小说| av在线蜜桃| 色综合亚洲欧美另类图片| 国内揄拍国产精品人妻在线| 啦啦啦免费观看视频1| 19禁男女啪啪无遮挡网站| 国产精品久久久久久人妻精品电影| 99精品久久久久人妻精品| 欧美日韩亚洲国产一区二区在线观看| 窝窝影院91人妻| 亚洲中文字幕一区二区三区有码在线看| 成人三级黄色视频| 亚洲熟妇中文字幕五十中出| 天美传媒精品一区二区| 亚洲国产色片| 中文字幕人妻熟人妻熟丝袜美 | 88av欧美| 亚洲av熟女| 欧美一级毛片孕妇| 国内精品久久久久精免费| 岛国在线观看网站| 日本与韩国留学比较| 久久久久久国产a免费观看| 男插女下体视频免费在线播放| 国产在线精品亚洲第一网站| 成人特级黄色片久久久久久久| 18禁黄网站禁片免费观看直播| 久久久久久国产a免费观看| 国产亚洲av嫩草精品影院| 亚洲久久久久久中文字幕| 亚洲成av人片在线播放无| 人妻夜夜爽99麻豆av| 国产精品一及| 俺也久久电影网| 日本一本二区三区精品| 三级国产精品欧美在线观看| 男女床上黄色一级片免费看| 一区二区三区免费毛片| 内地一区二区视频在线| 国产精品99久久久久久久久| 午夜影院日韩av| 18禁在线播放成人免费| 午夜激情欧美在线| 国产成+人综合+亚洲专区| 嫩草影视91久久| 舔av片在线| www国产在线视频色| 伊人久久精品亚洲午夜| 日本与韩国留学比较| 在线观看美女被高潮喷水网站 | 国产亚洲欧美在线一区二区| 午夜免费激情av| 精品久久久久久成人av| 内射极品少妇av片p| 久久国产乱子伦精品免费另类| 偷拍熟女少妇极品色| 亚洲精品456在线播放app | 美女大奶头视频| 丰满人妻一区二区三区视频av | 女警被强在线播放| 窝窝影院91人妻| 欧美一区二区精品小视频在线| 九色成人免费人妻av| 国产精品美女特级片免费视频播放器| 免费人成在线观看视频色| 午夜福利高清视频| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久末码| 久久久久性生活片| 啦啦啦观看免费观看视频高清| 国产视频一区二区在线看| 国产精品久久久久久精品电影| 美女cb高潮喷水在线观看| 日韩欧美国产在线观看| 一区二区三区高清视频在线| 少妇的逼好多水| 99精品久久久久人妻精品| 人人妻人人看人人澡| 一级a爱片免费观看的视频| 国产成人影院久久av| 国产综合懂色| 两个人的视频大全免费| 国产三级黄色录像| 欧美色欧美亚洲另类二区| 中文字幕av在线有码专区| 国产午夜精品久久久久久一区二区三区 | 成年人黄色毛片网站| 99热精品在线国产| 男女午夜视频在线观看| 久久6这里有精品| 亚洲男人的天堂狠狠| 网址你懂的国产日韩在线| 色噜噜av男人的天堂激情| 免费人成视频x8x8入口观看| 成熟少妇高潮喷水视频| 国产毛片a区久久久久| 3wmmmm亚洲av在线观看| 搡老妇女老女人老熟妇| 日本黄色片子视频| 好看av亚洲va欧美ⅴa在| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯| 成人亚洲精品av一区二区| 天堂av国产一区二区熟女人妻| 国产三级中文精品| 国产成人福利小说| 18禁黄网站禁片午夜丰满| 少妇的逼水好多| 久久精品国产自在天天线| 欧美一区二区亚洲| 俺也久久电影网| 天天一区二区日本电影三级| 波多野结衣高清无吗| 国产精品一及| 好男人在线观看高清免费视频| 最近最新中文字幕大全电影3| 动漫黄色视频在线观看| 成年女人永久免费观看视频| 久99久视频精品免费| 韩国av一区二区三区四区| 91av网一区二区| 嫩草影院入口| 国产激情欧美一区二区| 男女午夜视频在线观看| 国产精品久久久久久人妻精品电影| 午夜免费激情av| 99热精品在线国产| 久久精品人妻少妇| 男女视频在线观看网站免费| 黑人欧美特级aaaaaa片| 伊人久久大香线蕉亚洲五| 一夜夜www| 色尼玛亚洲综合影院| 大型黄色视频在线免费观看| 精品免费久久久久久久清纯| 久久伊人香网站| 国产精品久久久久久久电影 | 久久久久久久精品吃奶| 欧美中文综合在线视频| 中文亚洲av片在线观看爽| 在线观看免费视频日本深夜| 国内揄拍国产精品人妻在线| 日本一二三区视频观看| 在线国产一区二区在线| 国产精品免费一区二区三区在线| 成人鲁丝片一二三区免费| 亚洲中文字幕日韩| 免费看美女性在线毛片视频| 国产高清三级在线| 亚洲美女黄片视频| 黄片大片在线免费观看| 9191精品国产免费久久| 高清日韩中文字幕在线| 午夜激情欧美在线| 波野结衣二区三区在线 | 久9热在线精品视频| 欧美+日韩+精品| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 亚洲专区国产一区二区| 免费大片18禁| 最新在线观看一区二区三区| 国产一级毛片七仙女欲春2| 啪啪无遮挡十八禁网站| 国产欧美日韩精品一区二区| www日本黄色视频网| 国产欧美日韩一区二区精品| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 婷婷丁香在线五月| 中文字幕av成人在线电影| 给我免费播放毛片高清在线观看| 欧美日韩瑟瑟在线播放| 国产黄a三级三级三级人| e午夜精品久久久久久久| 国产一级毛片七仙女欲春2| 国内毛片毛片毛片毛片毛片| 国产又黄又爽又无遮挡在线| 久久久久亚洲av毛片大全| 久久久久国内视频| 伊人久久精品亚洲午夜| 色综合亚洲欧美另类图片| 99久久精品国产亚洲精品| 国产亚洲精品久久久久久毛片| 制服丝袜大香蕉在线| 亚洲在线观看片| 少妇丰满av| aaaaa片日本免费| 欧美最新免费一区二区三区 | 免费无遮挡裸体视频| 一个人免费在线观看的高清视频| 亚洲在线观看片| 国产欧美日韩一区二区三| 很黄的视频免费| 在线观看66精品国产| 噜噜噜噜噜久久久久久91| 欧美大码av| 国产一区在线观看成人免费| 成人av在线播放网站| 欧美区成人在线视频| 日韩欧美免费精品| 国产久久久一区二区三区| av福利片在线观看| 久久久久久久午夜电影| 国产高清videossex| 无遮挡黄片免费观看| 精品一区二区三区视频在线 | 俄罗斯特黄特色一大片| 最新美女视频免费是黄的| 欧美日本亚洲视频在线播放| 一级作爱视频免费观看| 手机成人av网站| 99国产综合亚洲精品| 久久久久久久亚洲中文字幕 | 最好的美女福利视频网| 性欧美人与动物交配| av欧美777| 国产精品永久免费网站| 午夜福利欧美成人| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 亚洲国产色片| 麻豆久久精品国产亚洲av| 在线观看日韩欧美| 1000部很黄的大片| 亚洲人与动物交配视频| 久久精品人妻少妇| 日本撒尿小便嘘嘘汇集6| 少妇熟女aⅴ在线视频| 叶爱在线成人免费视频播放| 观看免费一级毛片| 日韩欧美精品v在线| 热99re8久久精品国产| 午夜免费成人在线视频| 久久精品国产99精品国产亚洲性色| 男女下面进入的视频免费午夜| 熟女少妇亚洲综合色aaa.| 久久精品夜夜夜夜夜久久蜜豆| 一卡2卡三卡四卡精品乱码亚洲| 男女做爰动态图高潮gif福利片| 九色国产91popny在线| 在线免费观看的www视频| 国产成人a区在线观看| 欧美黑人欧美精品刺激| 可以在线观看毛片的网站| 九色成人免费人妻av| 国内久久婷婷六月综合欲色啪| 国产真人三级小视频在线观看| 99久久精品国产亚洲精品| 午夜日韩欧美国产| 国产一区二区在线av高清观看| 一卡2卡三卡四卡精品乱码亚洲| 在线视频色国产色| 成人鲁丝片一二三区免费| 色哟哟哟哟哟哟| 免费观看的影片在线观看| 国产视频一区二区在线看| 亚洲精华国产精华精| 网址你懂的国产日韩在线| 色播亚洲综合网| 一个人看的www免费观看视频| 69av精品久久久久久| 成人国产综合亚洲| 丝袜美腿在线中文| 少妇高潮的动态图| 女人十人毛片免费观看3o分钟| 欧美精品啪啪一区二区三区| а√天堂www在线а√下载| 国产午夜精品论理片| 听说在线观看完整版免费高清| 国产亚洲精品av在线| 国产精品自产拍在线观看55亚洲| 欧美高清成人免费视频www| 国产亚洲欧美在线一区二区| a级一级毛片免费在线观看| 19禁男女啪啪无遮挡网站| 可以在线观看毛片的网站| 少妇的逼好多水| 日日夜夜操网爽| 成年女人永久免费观看视频| 精品国内亚洲2022精品成人| 九色成人免费人妻av| 免费人成在线观看视频色| 床上黄色一级片| 99久久精品一区二区三区| 欧美一区二区国产精品久久精品| 人人妻人人看人人澡| 成年女人永久免费观看视频| 一级毛片高清免费大全| 午夜福利免费观看在线| 韩国av一区二区三区四区| 听说在线观看完整版免费高清| 亚洲无线观看免费| 一夜夜www| 手机成人av网站| 免费看十八禁软件| 亚洲专区国产一区二区| 男女视频在线观看网站免费| 婷婷精品国产亚洲av| 欧美成人一区二区免费高清观看| 九九在线视频观看精品| 俺也久久电影网| 欧美日韩黄片免| 18禁黄网站禁片午夜丰满| av中文乱码字幕在线| 蜜桃亚洲精品一区二区三区| 国产精品香港三级国产av潘金莲| 欧美中文综合在线视频| 他把我摸到了高潮在线观看| 制服人妻中文乱码| 国产精品,欧美在线| 一级毛片女人18水好多| 麻豆国产97在线/欧美| 18+在线观看网站| 欧美一区二区亚洲| 成人国产综合亚洲|