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

    基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法

    2015-10-21 01:50:28杜開(kāi)君茅云生向祖權(quán)周永清宋利飛
    船海工程 2015年3期
    關(guān)鍵詞:航速航向海事

    杜開(kāi)君,茅云生,向祖權(quán),周永清,宋利飛,劉 彬

    (武漢理工大學(xué)a.高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室;b.交通學(xué)院,武漢430063)

    基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法

    杜開(kāi)君a,b,茅云生a,b,向祖權(quán)a,b,周永清a,b,宋利飛b,劉 彬b

    (武漢理工大學(xué)a.高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室;b.交通學(xué)院,武漢430063)

    針對(duì)水面無(wú)人艇對(duì)探測(cè)到的電子海圖上沒(méi)有標(biāo)示的動(dòng)態(tài)障礙物的規(guī)避問(wèn)題,提出符合國(guó)際海上避碰規(guī)則公約的規(guī)避方法,將動(dòng)態(tài)障礙物某一運(yùn)動(dòng)時(shí)刻轉(zhuǎn)換為相對(duì)無(wú)人艇瞬時(shí)靜止的狀態(tài),對(duì)避障模型進(jìn)行實(shí)時(shí)計(jì)算;并通過(guò)粒子群優(yōu)化算法解算出無(wú)人艇進(jìn)行避障所需的最優(yōu)航速和航向改變量。仿真實(shí)驗(yàn)表明,該方法可得到合理的規(guī)避策略,有效地完成動(dòng)態(tài)障礙物的規(guī)避任務(wù)。

    水面無(wú)人艇;動(dòng)態(tài)障礙物;粒子群優(yōu)化算法;路徑規(guī)劃;海事規(guī)則

    水面無(wú)人艇(unmanned surface vehicle,USV)作為一種無(wú)人海上航行器,是惟一一種可以在水上、水面、水下同時(shí)發(fā)揮作用的無(wú)人運(yùn)載平臺(tái)[1]。在航行過(guò)程中,水面無(wú)人艇處在不可預(yù)見(jiàn)的動(dòng)態(tài)的海洋環(huán)境中,障礙物很可能出現(xiàn)在已經(jīng)規(guī)劃好的全局路徑上,而AIS或航海雷達(dá)并沒(méi)有探測(cè)到,電子海圖上也沒(méi)有標(biāo)示障礙物及其運(yùn)動(dòng)信息。為了保證無(wú)人艇能夠快速、安全地到達(dá)目標(biāo),需要對(duì)使其處在危險(xiǎn)狀況下的動(dòng)態(tài)障礙物進(jìn)行規(guī)避,這正是決定無(wú)人艇智能化程度高低的關(guān)鍵技術(shù)之一[1-5]。在分析各種無(wú)人艇障礙規(guī)避方法的基礎(chǔ)上,提出基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法。該方法根據(jù)國(guó)際海上航行規(guī)則公約將碰撞情形劃分為追越、正面相遇、交叉相遇3種局面;引入極坐標(biāo)系下的粒子群優(yōu)化算法,通過(guò)改變無(wú)人艇航速和航向,實(shí)現(xiàn)規(guī)避動(dòng)態(tài)障礙物的目的。

    1 極坐標(biāo)系下的建模

    常用的兩種平面坐標(biāo)系為笛卡爾坐標(biāo)系和極坐標(biāo)系。笛卡爾坐標(biāo)系在建模時(shí)具有局限性,如只能向前走,容易掉入陷阱;而極坐標(biāo)空間則對(duì)長(zhǎng)度和角度敏感,可以將各維分量與路徑處的期望航速和航向關(guān)聯(lián),并且便于描述障礙物信息。所以,采取極坐標(biāo)系進(jìn)行建模。

    極坐標(biāo)系下的起點(diǎn)和終點(diǎn)分別為在全局規(guī)劃路徑中選取的子目標(biāo)段的起點(diǎn)和終點(diǎn)。這兩點(diǎn)的連線為極坐標(biāo)系的極軸,起始點(diǎn)為極點(diǎn),逆時(shí)針?lè)较驗(yàn)檎?/p>

    真實(shí)環(huán)境下,障礙物的形狀各種各樣,通過(guò)雷達(dá)掃描出來(lái)的雷達(dá)圖像或通過(guò)聲吶、傳感器得到的障礙物形狀會(huì)很復(fù)雜,直接對(duì)這些未經(jīng)處理的障礙物形狀進(jìn)行避障,其計(jì)算過(guò)程往往會(huì)代價(jià)高昂,導(dǎo)致計(jì)算時(shí)間增加。為了縮短計(jì)算時(shí)間,在進(jìn)行路徑規(guī)劃之前,需要對(duì)障礙物進(jìn)行包圍體處理,使障礙物形狀得到合理的簡(jiǎn)化。

    國(guó)內(nèi)大部分學(xué)者對(duì)障礙物包圍體的處理僅采取圓形包圍處理,未充分考慮到海面上運(yùn)行的船舶均為長(zhǎng)寬比較大的障礙,僅采取圓形包圍會(huì)損失大量可行區(qū)間,從而造成避障策略不合理。這里對(duì)障礙物的包圍體處理包括兩種方式。

    1)對(duì)長(zhǎng)寬比較小(取小于2∶1)的障礙物來(lái)說(shuō),采用圓形進(jìn)行障礙包圍,不會(huì)損失太多的可行路徑;且極大地簡(jiǎn)少了計(jì)算量,可有效縮短計(jì)算時(shí)長(zhǎng)。外包圓具體取值方法為:其直徑取障礙物所有頂點(diǎn)連線的最大值,即

    式中:P——障礙頂點(diǎn)且(i,j∈N,i≠j),圓心即為最大連線的中點(diǎn)。

    2)對(duì)長(zhǎng)寬比較大(大于2∶1)的障礙物來(lái)說(shuō),將障礙物進(jìn)行更貼合實(shí)際形狀的模擬(即有向包圍盒,OBB[6])。雖然增大了計(jì)算量,但采用此方法,可得到更多的有效路徑區(qū)間,產(chǎn)生更優(yōu)的避障路徑,使得避障能耗更小、續(xù)航力更久。

    無(wú)人艇在海面航行過(guò)程中,會(huì)存在許多不確定性因素,如障礙物運(yùn)動(dòng)的不確定性、無(wú)人艇本身運(yùn)動(dòng)的不確定性、執(zhí)行機(jī)構(gòu)運(yùn)行時(shí)的不確定性、無(wú)人艇自身定位和對(duì)障礙物進(jìn)行探測(cè)的系統(tǒng)也不可避免地存在誤差等。為了簡(jiǎn)化計(jì)算,將水面無(wú)人艇作為一個(gè)質(zhì)心處理,所以還需要對(duì)障礙物進(jìn)行膨脹建模。

    2 動(dòng)態(tài)障礙物避障模型

    無(wú)人艇和圓形動(dòng)態(tài)障礙物當(dāng)前時(shí)刻的幾何模型見(jiàn)圖1。圖中:VUSV為無(wú)人艇的運(yùn)動(dòng)速度;VObs為障礙物的運(yùn)動(dòng)速度;α為極軸到VUSV的角度,可表示為α=∠(VUSV,ex),其中:ex為極軸;同理β=∠(VObs,ex)為極軸到VObs的角度;θ=∠(LRO,ex)為極軸到無(wú)人艇與障礙物圓心連線(LRO)的角度;φ=∠(VUSV,ΔV)為ΔV到VUSV的角度;γ=∠(ΔV,LRO),μ=∠(LRO,tan L),tan L為障礙圓切線。

    圖1 避障模型示意

    為使無(wú)人艇能在下一時(shí)刻避開(kāi)障礙區(qū)間,γ應(yīng)該取在(LRO-μ,LRO+μ)范圍之外的角度。下面將通過(guò)求解γ,找到無(wú)人艇航速和轉(zhuǎn)向需要滿足的條件,以實(shí)現(xiàn)無(wú)人艇避障。

    如圖2所示,將ΔV分解為指向障礙物圓心的速度分量ΔV0,和其垂直速度分量ΔVγ。其中ΔV0促使無(wú)人艇朝向障礙物運(yùn)動(dòng),ΔVγ促使無(wú)人艇逃離障礙物。

    圖2 對(duì)速度ΔV分解示意圖

    利用VUSV、VObs、ΔV構(gòu)成的速度三角形求解ΔVr和γ。

    式中:

    式(5)中,假設(shè)圓形動(dòng)態(tài)障礙物運(yùn)動(dòng)不發(fā)生突變,其運(yùn)動(dòng)速度VObs和運(yùn)動(dòng)方向β在很短的時(shí)間變化很小,可以忽略不計(jì),d VObs=0,dβ=0,則有

    所以式(3)可以轉(zhuǎn)化為用差分形式表示為

    在VUSV、VObs、ΔV構(gòu)成的速度三角形(如圖2)中,存在式(11)關(guān)系,其中φ為VUSV、VObs夾角,將其帶入式(10),可得到式(12)。

    要使無(wú)人艇逃離危險(xiǎn)區(qū)域(LRO-μ,LRO+μ)范圍,Δγ需要滿足下面的不等式。

    式(12)說(shuō)明Δγ同無(wú)人艇的航速改變量VUSV和航向改變量Δα存在關(guān)聯(lián),于是避障任務(wù)就轉(zhuǎn)化成調(diào)整無(wú)人艇的航速VUSV和航向Δα。

    3 國(guó)際海上避碰規(guī)則公約的約束

    水面無(wú)人艇在海面航行過(guò)程中,為了避免對(duì)海上其他船舶及無(wú)人艇自身的安全造成威脅,在進(jìn)行避障的過(guò)程中必須滿足相關(guān)規(guī)則。國(guó)際海事組織制定了國(guó)際海上避碰規(guī)則公約(International Regulations for Preventing Collisions At Sea,1972/ COLREGS),但現(xiàn)在還沒(méi)有專(zhuān)門(mén)的法律法規(guī)對(duì)無(wú)人艇的航行行為做出明確的規(guī)定,較合理的解決方法就是使無(wú)人艇遵守國(guó)際海上避碰規(guī)則公約[7]。

    國(guó)際海上避碰規(guī)則公約的13、14、15條分別對(duì)無(wú)人艇在航行過(guò)程中可能出現(xiàn)的碰撞局面:追越、正面相遇和交叉相遇,做出了規(guī)定。

    圖3給出了4種情況的避障模型,障礙物A狀態(tài)為無(wú)人艇在起始點(diǎn)“Start”處探測(cè)到該障礙物的位置,箭頭為其運(yùn)動(dòng)方向;障礙物B狀態(tài)則為無(wú)人艇對(duì)該障礙物避障結(jié)束時(shí),障礙物所處的位置。

    以下無(wú)人艇和運(yùn)動(dòng)障礙物之間的航向角度差,均以無(wú)人艇航向方向?yàn)榛鶞?zhǔn)(即航向角為0°)進(jìn)行描述。海事規(guī)則沖突情景定義見(jiàn)圖4。

    1)追越。若無(wú)人艇高速航行與障礙物航向角度差在[0°,45°]與 [315°,360°)之間時(shí),接近危險(xiǎn)距離,則無(wú)人艇左側(cè)航行。

    圖3 避障模型示意

    圖4 海事規(guī)則沖突情景定義

    2)正面相遇。若夾角在[165°,195°]之間,接近危險(xiǎn)距離,則無(wú)人艇應(yīng)右側(cè)航行,見(jiàn)圖3b)。

    3)交叉相遇。若夾角在(45°,165°)之間,接近危險(xiǎn)距離,則障礙物在無(wú)人艇右側(cè)交叉相遇,無(wú)人艇應(yīng)右側(cè)航行,見(jiàn)圖3c);若夾角在(195°, 315°)之間,接近危險(xiǎn)距離,則障礙物在無(wú)人艇左側(cè)交叉相遇,此時(shí)無(wú)人艇應(yīng)左側(cè)航行,見(jiàn)圖3d)。

    4 動(dòng)態(tài)障礙規(guī)避

    4.1 粒子群優(yōu)化算法及其適應(yīng)度函數(shù)的確定

    粒子群優(yōu)化算法是一種集群優(yōu)化算法。假設(shè)有一群隨機(jī)粒子以一定的速度飛行,每個(gè)粒子當(dāng)前自己找到的最優(yōu)解為Pi,d,全局最優(yōu)解為Pg,d。其中:d=1,2,…,N。N為粒子群迭代次數(shù)[8]。粒子在找到這兩個(gè)最優(yōu)解后,按式(14)、(15)來(lái)分別更新速度和位置。

    式中:r1和r2為0-1之間均勻分布的隨機(jī)數(shù)。w是慣性權(quán)值,文中w的取值選取廣泛采用的線性遞減慣性權(quán)重[9];c1和c2為學(xué)習(xí)因子的取值,經(jīng)多次試驗(yàn),均取值為2.0。

    無(wú)人艇在航行過(guò)程中,較多使用的是不改變航速,僅改變航向;但在某些情況下還需要同時(shí)改變航速和航向的策略,進(jìn)行避障。航速和航向改變量即為粒子群優(yōu)化算法的兩個(gè)維度值,無(wú)人艇路徑規(guī)劃問(wèn)題就可以表達(dá)成一個(gè)多條件下的目標(biāo)優(yōu)化問(wèn)題。

    式中:φ、μ、γ、ΔV已知,f(ΔVUSV,Δα)為要實(shí)現(xiàn)最優(yōu)的目標(biāo),即為粒子群優(yōu)化算法中的適應(yīng)度函數(shù), m1、m2分別為航速改變量和航向改變量的權(quán)值。計(jì)算得到最優(yōu)解后還需判斷是否在可行域內(nèi),即需要滿足國(guó)際海上避碰規(guī)則公約的約束,無(wú)人艇最小回轉(zhuǎn)直徑的約束,以及不與全局障礙物碰撞的約束。若該段路徑中存在多個(gè)動(dòng)態(tài)障礙物,則只需對(duì)每個(gè)障礙物加入式(16)中第2、3個(gè)約束即可。

    4.2 避障結(jié)束點(diǎn)的確定

    在對(duì)某一子目標(biāo)段內(nèi)的動(dòng)態(tài)障礙物實(shí)現(xiàn)成功規(guī)避后,無(wú)人艇需要沿著全局規(guī)劃路徑點(diǎn)繼續(xù)前行;采用上述的算法得到實(shí)現(xiàn)避障后的最優(yōu)航速和航向,所以還需確定無(wú)人艇成功避障時(shí)刻的避障結(jié)束點(diǎn)。

    若障礙物是朝向無(wú)人艇運(yùn)行,在未到達(dá)避障結(jié)束點(diǎn)之前,障礙物在全局中的極徑在不斷縮小,見(jiàn)圖5a)。圖中ˉVusv為經(jīng)過(guò)優(yōu)化后的無(wú)人艇速度、ΔˉV為優(yōu)化后無(wú)人艇和障礙物的合速度。此時(shí)避障結(jié)束點(diǎn)的選取,可根據(jù)在避障模型中,將障礙物視為靜止障礙物,以合速度ΔˉV到達(dá)障礙物最遠(yuǎn)端(圖中1/4圓弧段處)所耗費(fèi)的時(shí)間t,作為避障結(jié)束所耗費(fèi)的時(shí)間;進(jìn)而轉(zhuǎn)換到全局坐標(biāo)系下,以此時(shí)間段t計(jì)算避障結(jié)束點(diǎn)的極徑值和極角。

    若障礙物正遠(yuǎn)離起始點(diǎn)運(yùn)行,在未到達(dá)避障結(jié)束點(diǎn)之前,障礙物在全局中的極徑在不斷增大,見(jiàn)圖5b)。此情況下若還按上述方法進(jìn)行處理,無(wú)疑會(huì)造成避障結(jié)束點(diǎn)遠(yuǎn)離本段子目標(biāo)點(diǎn)。所以在此情況下,設(shè)置避障結(jié)束點(diǎn)即為本段起始點(diǎn)時(shí)障礙物所處位置的極徑處為避障結(jié)束點(diǎn)處的極徑。

    綜合以上兩種求避障結(jié)束點(diǎn)的情形,在計(jì)算時(shí)選取兩者極徑較小的那一點(diǎn)作為避障結(jié)束點(diǎn)。

    以上為針對(duì)動(dòng)態(tài)圓形障礙物的避障結(jié)束點(diǎn)選取方法,動(dòng)態(tài)矩形障礙物的避障結(jié)束點(diǎn)選取方法相同。

    圖5 避障結(jié)束點(diǎn)示意

    4.3 算法描述

    結(jié)合文中對(duì)基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法,給出該算法的流程,見(jiàn)圖6。

    圖6 基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法流程

    在粒子群優(yōu)化算法中,可行域判斷是指要同時(shí)滿足海事規(guī)則、無(wú)人艇最小回轉(zhuǎn)直徑及不與全局障礙物碰撞的約束。為了避免算法陷入局部最優(yōu),降低解算時(shí)間,在可行域的判斷中設(shè)置海事規(guī)則判斷最大次數(shù);若判斷次數(shù)大于設(shè)置的最大次數(shù),則不對(duì)海事規(guī)則進(jìn)行考慮。

    5 仿真實(shí)驗(yàn)結(jié)果及分析

    為了驗(yàn)證基于海事規(guī)則的水面無(wú)人艇動(dòng)態(tài)障礙規(guī)避方法,采用在Visual C++平臺(tái)上開(kāi)發(fā)的水面無(wú)人艇危險(xiǎn)避障算法仿真實(shí)驗(yàn)平臺(tái)進(jìn)行仿真實(shí)驗(yàn)分析。設(shè)立了如圖7所示的避障環(huán)境,該避障環(huán)境包括了海事規(guī)則中正面相遇、交叉相遇、追越的3種情形,以及多障礙路徑段的避障。

    圖7 仿真分析示意

    在該仿真實(shí)驗(yàn)中,假定水面無(wú)人艇已經(jīng)得到了全局最優(yōu)路徑節(jié)點(diǎn)Path1、Path2、Path3,“Start”為起始點(diǎn),“Goal”為終止點(diǎn),相距3 000 m;無(wú)人艇運(yùn)行速度為40 kn,最大航速為50 kn。粒子群優(yōu)化算法中慣性權(quán)重w采用線型遞減策略,wmax= 0.95,wmin=0.2;c1=c2=2;m1=1,m2=70;設(shè)置粒子群迭代次數(shù)50,粒子數(shù)50。

    該環(huán)境下包括一個(gè)全局靜態(tài)圓形障礙物SOC,和一個(gè)全局靜態(tài)矩形障礙物SOR。經(jīng)過(guò)全局路徑規(guī)劃方法解算,得到全局規(guī)劃路徑節(jié)點(diǎn)見(jiàn)表1。包括1個(gè)動(dòng)態(tài)圓形障礙物DOC1(見(jiàn)表2)和3個(gè)動(dòng)態(tài)矩形障礙物DOR1、DOR2、DOR3(見(jiàn)表3)(分別選取上海外高橋制造的85 000 m3液化氣船,18 000 TEU集裝箱船和106 kDWT阿芙拉型原油船)。為了方便驗(yàn)證,將動(dòng)態(tài)障礙物的起始位置均定為無(wú)人艇起始點(diǎn)的位置。動(dòng)態(tài)障礙物起始時(shí)刻的位置用Ai表示,避障結(jié)束點(diǎn)時(shí)刻障礙物的位置用Bi表示,箭頭方向?yàn)檎系K物運(yùn)動(dòng)方向。

    表1 全局規(guī)劃路徑節(jié)點(diǎn)

    表2 圓形動(dòng)態(tài)障礙物DOC參數(shù)

    表3 矩形動(dòng)態(tài)障礙物DOR參數(shù)

    表4 動(dòng)態(tài)障礙規(guī)避路徑節(jié)點(diǎn)

    無(wú)人艇在起始時(shí)刻時(shí)沿子目標(biāo)點(diǎn)Path1航行時(shí),探測(cè)到已知?jiǎng)討B(tài)矩形障礙物DOR1,并與無(wú)人艇呈185.17°相向運(yùn)行,即夾角在165°和195°之間,為海事規(guī)則中的正面相遇情形。若無(wú)人艇繼續(xù)向前行進(jìn),會(huì)同障礙物相遇碰撞,所以需要無(wú)人艇做出相應(yīng)的調(diào)整,從障礙物左舷即無(wú)人艇右轉(zhuǎn)向通行。通過(guò)該算法解算,耗時(shí)0.081 s;得到的避障策略為向右舷轉(zhuǎn)向11.34°,行進(jìn)航速不變。避障結(jié)束點(diǎn)為P-AD0,其極坐標(biāo)值見(jiàn)表4。圖7中較粗的實(shí)線即為無(wú)人艇改變航速、航向進(jìn)行避障的航線,成功避障后逐漸恢復(fù)航速至執(zhí)行任務(wù)航速,朝向本段的全局子目標(biāo)終點(diǎn)航行。

    無(wú)人艇繼續(xù)航行到Path1位置時(shí),探測(cè)到動(dòng)態(tài)圓形障礙物DOC1,和動(dòng)態(tài)矩形障礙物DOR2,并分別從無(wú)人艇兩側(cè)駛來(lái),角度分別為75.17°和254.17°;為海事規(guī)則中的交叉相遇情形,需要分別從兩個(gè)障礙物的后方通過(guò),即對(duì)DOC1而言為右側(cè)通行,對(duì)DOR2而言為左側(cè)通行。無(wú)人艇在Path1時(shí)刻進(jìn)行避障,算法解算耗時(shí)0.119 s;得到的避障策略為航向不改變,航速?gòu)?0 kn減少到27.1 kn行進(jìn)。無(wú)人艇運(yùn)行到避障結(jié)束點(diǎn)P-AD1后,就已成功完成了對(duì)障礙物DOC1和DOR2的規(guī)避,此時(shí)便逐漸加速到40 kn繼續(xù)朝向本段的子目標(biāo)點(diǎn)Path2航行。

    無(wú)人艇運(yùn)行到Path3位置時(shí),探測(cè)到已知?jiǎng)討B(tài)矩形障礙物DOR3,并與無(wú)人艇呈32.06°運(yùn)行,即角度差在(0°,45°)之間,且無(wú)人艇航速(40 kn)大于障礙物速度(15 kn),為海事規(guī)則中的追越情形。此時(shí)需要無(wú)人艇做出相應(yīng)的調(diào)整,從障礙物左舷即無(wú)人艇左轉(zhuǎn)向通行。所以無(wú)人艇在Path3時(shí)刻通過(guò)該算法解算,耗時(shí)0.078 s;得到的避障策略為向左舷轉(zhuǎn)向7.62°,航速增加到42.0 kn行進(jìn)。無(wú)人艇運(yùn)行到避障結(jié)束點(diǎn)P-AD3后,就已成功完成了對(duì)障礙物DOR3的規(guī)避。此時(shí)便逐漸恢復(fù)至執(zhí)行任務(wù)航速40 kn,轉(zhuǎn)向本段的子目標(biāo)點(diǎn)即終點(diǎn)航行。

    6 結(jié)束語(yǔ)

    對(duì)障礙物進(jìn)行包圍體處理,使障礙物形狀得到合理簡(jiǎn)化,減少了計(jì)算消耗,縮短了計(jì)算時(shí)間。且根據(jù)障礙物長(zhǎng)寬比,采用圓形包圍體或有向包圍盒對(duì)障礙物進(jìn)行簡(jiǎn)化;較好地解決了國(guó)內(nèi)學(xué)者對(duì)無(wú)人艇會(huì)遇障礙物簡(jiǎn)化時(shí),僅采用單一的圓形包圍體處理而導(dǎo)致的損失大量可行區(qū)間、產(chǎn)生不合理避障策略的問(wèn)題。

    將國(guó)際海上避碰行為規(guī)則作為避障約束條件。通過(guò)仿真實(shí)驗(yàn)可得出,該算法不僅對(duì)單一障礙物避障適用,同時(shí)針對(duì)多障礙物規(guī)避問(wèn)題,也可得到合理的避障策略。

    將粒子群優(yōu)化算法引入無(wú)人艇動(dòng)態(tài)障礙規(guī)避中,很好地發(fā)揮了該算法在解決多目標(biāo)優(yōu)化問(wèn)題中操作流程簡(jiǎn)單、參數(shù)簡(jiǎn)潔、收斂效果好的特點(diǎn)。通過(guò)仿真實(shí)驗(yàn)得出,避障解算時(shí)間在0.15 s之內(nèi);且根據(jù)避障情景可解算出僅改變航速、僅改變航向、同時(shí)改變航速和航向的3種可行策略。

    仿真實(shí)驗(yàn)結(jié)果表明,該算法能夠很好地完成對(duì)動(dòng)態(tài)障礙物的規(guī)避,具有較強(qiáng)的環(huán)境適應(yīng)能力、較快的解算速度、較好的避障效果,滿足無(wú)人艇對(duì)障礙物規(guī)避的要求。然而由于部分障礙物存在運(yùn)動(dòng)不確定性,預(yù)測(cè)到的障礙物運(yùn)動(dòng)信息難免出現(xiàn)不準(zhǔn)確的情況。針對(duì)這些障礙物若僅采用文中提出的算法,不能保障無(wú)人艇成功進(jìn)行避障。所以下一步需要解決的問(wèn)題是:無(wú)人艇根據(jù)視覺(jué)傳感器或避障聲吶等障礙識(shí)別系統(tǒng)實(shí)時(shí)探測(cè)并反饋障礙物當(dāng)前位置信息,在文中提出算法的基礎(chǔ)上,對(duì)阻礙無(wú)人艇航行的障礙物做出敏捷的應(yīng)激反應(yīng)式避障。

    [1]YAN R,PANG S,S-UN H,et al.Development and missions of unmanned surface vehicle[J].Journal ofMarine Science and Application,2010,9(4):451-457.

    [2]BENJAMIN M R,CURCIO J A,LEONARD J J,et al. Navigation of unmanned marine vehicles in accordance with the rules of the road[C]∥Proceedings of the 2006 IEEE International Conference on Robotics and Automation.Orlando,Florida,USA,2006:3581-3587.

    [3]莊佳園,蘇玉民,廖煜雷,等.基于航海雷達(dá)的水面無(wú)人艇局部路徑規(guī)劃[J].上海交通大學(xué)學(xué)報(bào),2012,46 (9):1371-1375.

    [4]盧艷爽.水面無(wú)人艇路徑規(guī)劃算法研究[D].哈爾濱:哈爾濱工程大學(xué),2010.

    [5]馬 闖,殷 波,馬文帥.水上機(jī)器人三維實(shí)時(shí)避障算法研究[J].微計(jì)算機(jī)信息,2009,25(3-2):235-237.

    [6]CASALINO G,TURETTA A,SIMETTI E.A three-layered architecture for real time path planning and obstacle avoidance for surveillance USVs operating in harbour fields[C]∥OCEANS 2009-EUROPE.Bremen,German,2009:1-8.

    [7]LEE S,KWON K,JOH J.A fuzzy logic for autonomous navigation ofmarine vehicles satisfying COLREG guidelines[J].International Journal of Control Automation and Systems,2004,2(2):171-181.

    [8]蘇紹娟,王麗錚,王呈方.不確定性航線配船數(shù)學(xué)模型建模方法[J].船海工程,2007,36(4):100-103.

    [9]張麗平.粒子群優(yōu)化算法的理論及實(shí)踐[D].杭州:浙江大學(xué),2005.

    A Dynamic Obstacle Avoidance Method of USV Based on COLREGS

    DU Kai-juna,b,MAO Yun-shenga,b,XIANG Zu-quana,b,ZHOU Yong-qinga,b,SONG Li-feib,LIU Binb
    (a.Key Laboratory of High Performance Ship Technology; b.School of Transportation,Wuhan University of Technology,Wuhan 430063,China)

    To solve the problem for the unmanned surface vehicle(USV)avoiding the dynamic obstaclewhich is notmarked on an electronic chart,a dynamic obstacle avoidancemethod complies the international regulations for preventing collisions atsea (COLREGS)are presented.The dynamic obstacles are converted into an instantaneous static state related to USV in a motion moment,so as to compute the optimal speed and heading in real-time for the USV to avoid the obstacles by particle swarm optimization algorithm.Simulation experiments demonstrate that the dynamic obstacle avoidancemethod can calculate a reasonable avoidance strategy and accomplish the dynamic obstacle avoidance task effectively.

    USV;dynamic obstacles;particle swarm optimization;path planning;maritime rule

    10.3963/j.issn.1671-7953.2015.03.028

    U675.85

    A

    1671-7953(2015)03-0119-06

    2014-12-16

    修回日期:2015-01-19

    杜開(kāi)君(1987-),男,碩士生

    研究方向:船舶先進(jìn)制造技術(shù)

    E-m ail:dkjkls@126.com

    猜你喜歡
    航速航向海事
    VLCC在波浪中的航速優(yōu)化與能效優(yōu)化分析
    提升全回轉(zhuǎn)港作拖輪航速的有效途徑
    水上消防(2022年1期)2022-06-16 08:06:56
    信息精要與海事掃描
    水上消防(2021年5期)2022-01-18 05:33:42
    信息精要與海事掃描
    水上消防(2021年3期)2021-08-21 03:12:22
    知坐標(biāo),明航向
    信息精要與海事掃描
    水上消防(2020年5期)2020-12-14 07:16:30
    信息精要與海事掃描
    水上消防(2020年3期)2020-07-25 02:36:32
    考慮幾何限制的航向道模式設(shè)計(jì)
    低速水面目標(biāo)航速精度分析及精確解算
    基于干擾觀測(cè)器的船舶系統(tǒng)航向Backstepping 控制
    電子制作(2017年24期)2017-02-02 07:14:16
    两个人看的免费小视频| 色视频在线一区二区三区| 大码成人一级视频| 亚洲一区二区三区欧美精品| 视频在线观看一区二区三区| av在线观看视频网站免费| 老鸭窝网址在线观看| 叶爱在线成人免费视频播放| 国产精品久久久久成人av| 大话2 男鬼变身卡| 久久精品国产鲁丝片午夜精品| 久久97久久精品| 国产av精品麻豆| 日韩视频在线欧美| 男男h啪啪无遮挡| 亚洲综合色网址| 亚洲美女视频黄频| 欧美亚洲 丝袜 人妻 在线| 亚洲第一青青草原| 亚洲少妇的诱惑av| 国产 精品1| 又粗又硬又长又爽又黄的视频| 成人毛片60女人毛片免费| 国产精品国产av在线观看| 国产精品国产三级国产专区5o| 啦啦啦中文免费视频观看日本| 大片电影免费在线观看免费| 只有这里有精品99| 国产免费现黄频在线看| 免费少妇av软件| 激情视频va一区二区三区| 亚洲,一卡二卡三卡| 亚洲欧洲国产日韩| 亚洲国产精品一区二区三区在线| 九色亚洲精品在线播放| 国产 一区精品| 国产精品二区激情视频| 精品国产超薄肉色丝袜足j| 精品卡一卡二卡四卡免费| 亚洲av.av天堂| 免费观看a级毛片全部| 国产精品一国产av| 夫妻性生交免费视频一级片| 日韩一区二区视频免费看| 丝袜脚勾引网站| 两个人看的免费小视频| 亚洲欧洲国产日韩| 制服诱惑二区| 两性夫妻黄色片| 亚洲激情五月婷婷啪啪| 亚洲av成人精品一二三区| 一级毛片电影观看| 菩萨蛮人人尽说江南好唐韦庄| 91久久精品国产一区二区三区| 久久97久久精品| 亚洲欧美清纯卡通| 一级a爱视频在线免费观看| 韩国av在线不卡| 久热久热在线精品观看| 午夜福利乱码中文字幕| 久久久久久久久久久久大奶| 日本-黄色视频高清免费观看| 亚洲欧美一区二区三区黑人 | 日韩伦理黄色片| 视频在线观看一区二区三区| 黄片无遮挡物在线观看| 日韩一卡2卡3卡4卡2021年| 日本爱情动作片www.在线观看| 97人妻天天添夜夜摸| 精品久久久精品久久久| 国产精品熟女久久久久浪| 在现免费观看毛片| 国产日韩欧美亚洲二区| 国产精品国产av在线观看| 九九爱精品视频在线观看| 一级毛片电影观看| 久久久精品94久久精品| 久久午夜福利片| 亚洲国产精品成人久久小说| 老司机亚洲免费影院| 黄片播放在线免费| 侵犯人妻中文字幕一二三四区| 黑人巨大精品欧美一区二区蜜桃| 捣出白浆h1v1| 久久久国产精品麻豆| 精品少妇黑人巨大在线播放| 1024视频免费在线观看| 亚洲欧洲国产日韩| www.熟女人妻精品国产| 啦啦啦在线观看免费高清www| 毛片一级片免费看久久久久| 香蕉精品网在线| 欧美精品高潮呻吟av久久| 午夜福利视频在线观看免费| 国产精品一区二区在线观看99| 国产成人精品一,二区| 岛国毛片在线播放| av免费观看日本| 激情五月婷婷亚洲| 天美传媒精品一区二区| 色视频在线一区二区三区| 成年人午夜在线观看视频| 国产欧美日韩一区二区三区在线| 国产精品国产av在线观看| 成人免费观看视频高清| 国产在线免费精品| 欧美精品高潮呻吟av久久| 性色av一级| 考比视频在线观看| 在线观看免费日韩欧美大片| freevideosex欧美| 精品久久久久久电影网| 日本黄色日本黄色录像| 91精品国产国语对白视频| 又大又黄又爽视频免费| av网站免费在线观看视频| 精品国产露脸久久av麻豆| 黑丝袜美女国产一区| 亚洲欧美一区二区三区黑人 | 成人黄色视频免费在线看| 哪个播放器可以免费观看大片| 午夜福利网站1000一区二区三区| 大陆偷拍与自拍| 久久精品国产亚洲av高清一级| 亚洲av国产av综合av卡| 日本免费在线观看一区| 青草久久国产| 久久这里只有精品19| 日本av手机在线免费观看| 中文字幕精品免费在线观看视频| 热re99久久精品国产66热6| 黄片小视频在线播放| 天天躁夜夜躁狠狠久久av| 天美传媒精品一区二区| 大香蕉久久网| 亚洲欧美色中文字幕在线| 免费黄色在线免费观看| 一级毛片我不卡| 久久国产亚洲av麻豆专区| av国产久精品久网站免费入址| 黄色毛片三级朝国网站| av卡一久久| 老司机亚洲免费影院| 免费看av在线观看网站| 涩涩av久久男人的天堂| 午夜福利,免费看| 国产精品av久久久久免费| 国产在线视频一区二区| 少妇的逼水好多| 18在线观看网站| 国产视频首页在线观看| 亚洲精品成人av观看孕妇| 成年动漫av网址| 亚洲国产毛片av蜜桃av| 色网站视频免费| 人人妻人人澡人人看| 777米奇影视久久| 久久97久久精品| 最近的中文字幕免费完整| 亚洲成人一二三区av| 一区二区三区激情视频| 精品国产一区二区三区四区第35| 亚洲av国产av综合av卡| 黑人欧美特级aaaaaa片| 不卡av一区二区三区| 热99国产精品久久久久久7| 亚洲图色成人| 成人亚洲精品一区在线观看| 一本—道久久a久久精品蜜桃钙片| 国产一区二区在线观看av| 午夜福利视频精品| 老女人水多毛片| 人人妻人人澡人人爽人人夜夜| 精品福利永久在线观看| 高清av免费在线| 春色校园在线视频观看| 国产精品香港三级国产av潘金莲 | 一级毛片电影观看| 丝袜脚勾引网站| 午夜精品国产一区二区电影| 丝袜美腿诱惑在线| 久久99精品国语久久久| av卡一久久| 极品少妇高潮喷水抽搐| 寂寞人妻少妇视频99o| 亚洲一级一片aⅴ在线观看| 春色校园在线视频观看| h视频一区二区三区| 国产av码专区亚洲av| 精品一区二区免费观看| 在现免费观看毛片| 中文字幕人妻丝袜一区二区 | 成人毛片60女人毛片免费| 18禁裸乳无遮挡动漫免费视频| 国产一区二区 视频在线| 不卡av一区二区三区| 亚洲精品国产av蜜桃| 国产日韩欧美亚洲二区| 在线天堂最新版资源| 亚洲成人av在线免费| 91aial.com中文字幕在线观看| 亚洲av电影在线观看一区二区三区| 水蜜桃什么品种好| 丝袜喷水一区| 边亲边吃奶的免费视频| 欧美国产精品一级二级三级| 婷婷色综合www| 久久精品久久久久久久性| 久久毛片免费看一区二区三区| 青草久久国产| 亚洲国产最新在线播放| 少妇人妻精品综合一区二区| 国产探花极品一区二区| 日韩中文字幕视频在线看片| 最近手机中文字幕大全| 日本91视频免费播放| 在线看a的网站| 久久精品久久久久久噜噜老黄| 人人澡人人妻人| 久久久久国产精品人妻一区二区| 在线观看一区二区三区激情| 人人妻人人澡人人爽人人夜夜| 久久久欧美国产精品| 久久综合国产亚洲精品| 啦啦啦在线观看免费高清www| 有码 亚洲区| 毛片一级片免费看久久久久| 9191精品国产免费久久| 大片免费播放器 马上看| 卡戴珊不雅视频在线播放| 不卡av一区二区三区| 99久国产av精品国产电影| 欧美精品亚洲一区二区| 国产精品久久久久久精品古装| videos熟女内射| 亚洲,欧美,日韩| 在线看a的网站| av片东京热男人的天堂| 午夜免费鲁丝| 天天躁夜夜躁狠狠躁躁| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 午夜福利视频精品| 97在线人人人人妻| 国产精品成人在线| 大片免费播放器 马上看| 熟妇人妻不卡中文字幕| 乱人伦中国视频| 亚洲国产精品一区二区三区在线| 欧美国产精品一级二级三级| 国产 精品1| 一区福利在线观看| 国产精品二区激情视频| 中文字幕色久视频| 婷婷成人精品国产| 久久久久久久久久久免费av| 男人操女人黄网站| 亚洲四区av| 久久毛片免费看一区二区三区| 亚洲国产精品一区二区三区在线| 极品人妻少妇av视频| 成人黄色视频免费在线看| 26uuu在线亚洲综合色| www.精华液| 一二三四在线观看免费中文在| 久热这里只有精品99| 永久网站在线| 色吧在线观看| 久久人人爽av亚洲精品天堂| 啦啦啦啦在线视频资源| 久久女婷五月综合色啪小说| 亚洲欧洲日产国产| av电影中文网址| 久久这里有精品视频免费| 日韩中字成人| 精品少妇久久久久久888优播| 国精品久久久久久国模美| 国产成人aa在线观看| 欧美日韩国产mv在线观看视频| 91精品伊人久久大香线蕉| 久久久久国产精品人妻一区二区| 成人影院久久| 国语对白做爰xxxⅹ性视频网站| 男的添女的下面高潮视频| 亚洲av.av天堂| 在线观看免费日韩欧美大片| 人人澡人人妻人| 亚洲精品成人av观看孕妇| 久久久久久久久久人人人人人人| 亚洲伊人久久精品综合| 精品卡一卡二卡四卡免费| av国产久精品久网站免费入址| 一个人免费看片子| 欧美xxⅹ黑人| 国产一区亚洲一区在线观看| 中国国产av一级| 欧美成人午夜精品| 久久国产精品大桥未久av| 伦理电影免费视频| 狂野欧美激情性bbbbbb| 激情视频va一区二区三区| av.在线天堂| 国产精品女同一区二区软件| 9色porny在线观看| 国产精品久久久久久久久免| 欧美黄色片欧美黄色片| 国产精品偷伦视频观看了| 69精品国产乱码久久久| 咕卡用的链子| 亚洲伊人色综图| 日韩电影二区| 久久99蜜桃精品久久| 黑人猛操日本美女一级片| 97精品久久久久久久久久精品| 欧美亚洲日本最大视频资源| 亚洲第一青青草原| 90打野战视频偷拍视频| 亚洲婷婷狠狠爱综合网| 男女国产视频网站| 久久久久久久久免费视频了| 欧美日韩精品网址| 久久婷婷青草| 五月开心婷婷网| 精品亚洲成国产av| 精品国产一区二区久久| 日产精品乱码卡一卡2卡三| 亚洲av男天堂| 有码 亚洲区| 春色校园在线视频观看| 在线免费观看不下载黄p国产| 韩国精品一区二区三区| 欧美av亚洲av综合av国产av | 一级毛片 在线播放| 亚洲av电影在线观看一区二区三区| 精品国产国语对白av| 婷婷色av中文字幕| 一本大道久久a久久精品| 超碰97精品在线观看| av在线老鸭窝| 又大又黄又爽视频免费| 男女高潮啪啪啪动态图| 免费观看无遮挡的男女| 久久久久网色| 久久国产亚洲av麻豆专区| 亚洲欧洲精品一区二区精品久久久 | 一级爰片在线观看| 日本免费在线观看一区| 国产精品一国产av| 国产精品三级大全| av国产久精品久网站免费入址| av电影中文网址| 麻豆乱淫一区二区| 亚洲三区欧美一区| 97精品久久久久久久久久精品| 99热网站在线观看| 两性夫妻黄色片| 久久久久久久亚洲中文字幕| 久久久久精品人妻al黑| 免费观看a级毛片全部| 久久久久久久久免费视频了| 在线免费观看不下载黄p国产| 欧美黄色片欧美黄色片| 天天躁夜夜躁狠狠久久av| 久久久久久伊人网av| 国产亚洲午夜精品一区二区久久| 日韩在线高清观看一区二区三区| 久久亚洲国产成人精品v| 色94色欧美一区二区| 99热网站在线观看| 一本大道久久a久久精品| 美女中出高潮动态图| 街头女战士在线观看网站| 亚洲欧洲国产日韩| 欧美日韩精品网址| 国产成人精品在线电影| 9色porny在线观看| 亚洲精品国产av蜜桃| 90打野战视频偷拍视频| 爱豆传媒免费全集在线观看| 丝瓜视频免费看黄片| 这个男人来自地球电影免费观看 | 26uuu在线亚洲综合色| 在线免费观看不下载黄p国产| 国产免费福利视频在线观看| 欧美人与性动交α欧美软件| 国产精品免费视频内射| 在线看a的网站| 一本久久精品| 欧美成人午夜免费资源| 国产成人午夜福利电影在线观看| 2018国产大陆天天弄谢| av在线播放精品| 午夜91福利影院| 久久久久久免费高清国产稀缺| 亚洲欧美成人精品一区二区| 观看av在线不卡| 高清不卡的av网站| 亚洲男人天堂网一区| 国产精品人妻久久久影院| 少妇猛男粗大的猛烈进出视频| 欧美黄色片欧美黄色片| 国产精品亚洲av一区麻豆 | 国产激情久久老熟女| 亚洲三区欧美一区| 交换朋友夫妻互换小说| 久久久精品94久久精品| 久久精品夜色国产| 国产精品av久久久久免费| 久久综合国产亚洲精品| 三上悠亚av全集在线观看| 人妻一区二区av| 国产男女内射视频| 少妇人妻精品综合一区二区| 国产成人精品福利久久| 国产精品免费大片| 久久久久久人妻| 不卡视频在线观看欧美| 精品一区二区免费观看| 天天躁夜夜躁狠狠躁躁| 国产熟女午夜一区二区三区| 亚洲四区av| 久久久国产精品麻豆| 国产精品一区二区在线不卡| 久久精品国产亚洲av天美| 黑丝袜美女国产一区| 精品少妇内射三级| 久久久久国产网址| 日本午夜av视频| 日韩中文字幕欧美一区二区 | 日韩,欧美,国产一区二区三区| 在线观看国产h片| 国产成人午夜福利电影在线观看| 伦理电影免费视频| 国产午夜精品一二区理论片| 激情五月婷婷亚洲| 宅男免费午夜| 久久人人爽av亚洲精品天堂| 深夜精品福利| 日本爱情动作片www.在线观看| 99国产精品免费福利视频| 青春草亚洲视频在线观看| 国产又色又爽无遮挡免| 国产亚洲午夜精品一区二区久久| 侵犯人妻中文字幕一二三四区| 亚洲精品aⅴ在线观看| 午夜日韩欧美国产| 亚洲精品久久成人aⅴ小说| 一级片'在线观看视频| 晚上一个人看的免费电影| 国产成人精品久久久久久| www.熟女人妻精品国产| 伊人久久大香线蕉亚洲五| 熟女电影av网| 丝袜美足系列| 午夜福利在线免费观看网站| 美女国产视频在线观看| 精品亚洲乱码少妇综合久久| 亚洲av中文av极速乱| 99九九在线精品视频| 久久国产精品男人的天堂亚洲| 久久久国产精品麻豆| 青青草视频在线视频观看| 久久婷婷青草| 国产成人91sexporn| 久久热在线av| 日韩精品有码人妻一区| av国产精品久久久久影院| 丝袜喷水一区| 哪个播放器可以免费观看大片| 亚洲,欧美,日韩| 欧美精品一区二区免费开放| 亚洲欧美一区二区三区久久| 欧美精品高潮呻吟av久久| 亚洲男人天堂网一区| 久久久久久久久久久久大奶| 一级片'在线观看视频| 亚洲av综合色区一区| 久久99一区二区三区| 9191精品国产免费久久| 国产av码专区亚洲av| www.精华液| av国产精品久久久久影院| 飞空精品影院首页| 国产一区二区激情短视频 | 久久精品久久久久久噜噜老黄| 国产精品一二三区在线看| 亚洲精品自拍成人| 日韩免费高清中文字幕av| 免费少妇av软件| 亚洲人成77777在线视频| 精品少妇内射三级| 国产一区二区在线观看av| 另类精品久久| 高清黄色对白视频在线免费看| 亚洲av电影在线观看一区二区三区| 天天影视国产精品| 国产精品成人在线| 久久久久久人人人人人| 国产精品秋霞免费鲁丝片| 一级a爱视频在线免费观看| 高清在线视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 午夜久久久在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲熟女精品中文字幕| 97精品久久久久久久久久精品| 亚洲国产精品999| 9色porny在线观看| 午夜免费观看性视频| av免费在线看不卡| 激情视频va一区二区三区| 成人黄色视频免费在线看| 日韩 亚洲 欧美在线| 最近中文字幕2019免费版| 最近中文字幕高清免费大全6| 国产又爽黄色视频| 男女边吃奶边做爰视频| 欧美在线黄色| 国产成人91sexporn| 我要看黄色一级片免费的| 国产av一区二区精品久久| 黄片播放在线免费| 一个人免费看片子| 777米奇影视久久| 欧美在线黄色| 久久久久视频综合| 晚上一个人看的免费电影| 日本wwww免费看| 国产亚洲一区二区精品| 亚洲av欧美aⅴ国产| 最近手机中文字幕大全| 国产毛片在线视频| 日韩成人av中文字幕在线观看| 巨乳人妻的诱惑在线观看| 欧美日韩视频精品一区| 国产视频首页在线观看| 少妇猛男粗大的猛烈进出视频| 老女人水多毛片| 国产成人精品无人区| 国产成人a∨麻豆精品| 亚洲少妇的诱惑av| 亚洲,欧美精品.| 亚洲国产精品一区二区三区在线| 最新中文字幕久久久久| 丰满少妇做爰视频| 1024香蕉在线观看| 中文字幕人妻丝袜制服| 日韩不卡一区二区三区视频在线| 性色av一级| 国产毛片在线视频| 久久久久久伊人网av| 一边摸一边做爽爽视频免费| 久久精品亚洲av国产电影网| 久久国内精品自在自线图片| 丝瓜视频免费看黄片| 青草久久国产| 久久精品人人爽人人爽视色| 欧美国产精品va在线观看不卡| 免费大片黄手机在线观看| 大片免费播放器 马上看| 国产欧美亚洲国产| 亚洲精品国产色婷婷电影| 欧美人与善性xxx| 丰满乱子伦码专区| 秋霞伦理黄片| 成人国产麻豆网| 少妇 在线观看| 亚洲av国产av综合av卡| 2021少妇久久久久久久久久久| 亚洲激情五月婷婷啪啪| 日韩 亚洲 欧美在线| 制服丝袜香蕉在线| 美女福利国产在线| 国产精品久久久久久久久免| 国产麻豆69| 午夜免费观看性视频| 国产一区二区三区av在线| 自线自在国产av| 91国产中文字幕| 丰满饥渴人妻一区二区三| 成人亚洲精品一区在线观看| 多毛熟女@视频| 亚洲欧美成人精品一区二区| 欧美另类一区| 欧美日韩精品网址| 香蕉国产在线看| 日韩制服丝袜自拍偷拍| 欧美日韩一区二区视频在线观看视频在线| 黄色毛片三级朝国网站| 日本-黄色视频高清免费观看| 国产综合精华液| 亚洲少妇的诱惑av| 欧美日本中文国产一区发布| 精品第一国产精品| 精品人妻在线不人妻| 夫妻午夜视频| 一区福利在线观看| 亚洲一级一片aⅴ在线观看| 亚洲av在线观看美女高潮| 国产男女内射视频| 99re6热这里在线精品视频| 欧美精品国产亚洲| 久久av网站| www.精华液| 免费少妇av软件| 黄色怎么调成土黄色| 午夜福利一区二区在线看| xxxhd国产人妻xxx| 久久久久久久久久久免费av| 91精品伊人久久大香线蕉| 午夜福利,免费看| 少妇人妻 视频| a级毛片黄视频| 久久久久久免费高清国产稀缺| 飞空精品影院首页|