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

    一種改進(jìn)的移動(dòng)機(jī)器人自定位算法

    2014-10-22 02:44:08蒙禾青包琨超王宗玉
    關(guān)鍵詞:移動(dòng)機(jī)器人粒子角度

    張 琪,蒙禾青,周 嶸,包琨超,王宗玉

    (1.武漢理工大學(xué) 信息工程學(xué)院,湖北 武漢 430070;2.西安交通大學(xué)軟件學(xué)院,陜西 西安 710049)

    自定位是移動(dòng)機(jī)器人的核心問(wèn)題。通過(guò)自定位,移動(dòng)機(jī)器人可以在特定環(huán)境中完成人類(lèi)指定的各種任務(wù),例如傳遞指令、取送物品、危險(xiǎn)物排除、管理辦公室或家居和機(jī)器人足球等[1-2],因此這項(xiàng)研究具有重要價(jià)值。移動(dòng)機(jī)器人的定位,需要兩個(gè)條件:一是感知內(nèi)、外環(huán)境的傳感器,例如筆者采用經(jīng)典的測(cè)距傳感器來(lái)獲取環(huán)境信息;二是關(guān)于行走環(huán)境的地圖。地圖的獲取有兩種方式,一是機(jī)器人在行走過(guò)程中根據(jù)實(shí)際量測(cè)來(lái)自動(dòng)繪制環(huán)境地圖;二是把已經(jīng)標(biāo)定好的環(huán)境地圖存放在移動(dòng)機(jī)器人的芯片中。為了簡(jiǎn)化問(wèn)題,筆者采用第二種方式,即地圖在自定位之前已經(jīng)獲得。地圖采用離散柵格形式,如圖1所示。其中黑色柵格表示障礙物,用數(shù)字0表示,而白色部分表示機(jī)器人可以通過(guò)的柵格,用數(shù)字1表示。

    圖1 離散柵格地圖示例

    地圖建立以后,機(jī)器人通過(guò)移動(dòng),不斷獲取與墻壁或周邊物體的距離信息,經(jīng)過(guò)與地圖的匹配,逐步判斷機(jī)器人在當(dāng)前時(shí)刻的位置和方向。由于機(jī)器人定位是一個(gè)非線性、含有較多不確定性因素的問(wèn)題,因此當(dāng)前主要采用基于貝葉斯濾波的概率推導(dǎo)方法,如卡爾曼濾波器、多峰值估計(jì)法、馬爾科夫定位和蒙特卡洛法等[3-4]。在實(shí)例方面,卡內(nèi)基梅隆大學(xué)與德國(guó)伯恩大學(xué)基于蒙特卡洛方法,合作研發(fā)出博物館導(dǎo)游機(jī)器人Rhino及Minerva,可以從博物館中任意一個(gè)位置出發(fā),較快地實(shí)現(xiàn)自定位[5]。

    國(guó)內(nèi)對(duì)于自定位問(wèn)題的研究起步相對(duì)較晚,主要集中在近10年內(nèi)。如利用馬爾科夫算法對(duì)機(jī)器人進(jìn)行定位;利用改進(jìn)的自適應(yīng)粒子濾波算法,降低運(yùn)算復(fù)雜度,提升算法穩(wěn)健性;利用特征粒子提取思想,降低粒子數(shù)目;采用改進(jìn)的Hough密度譜和Fouder-Mellin變換方法來(lái)估計(jì)機(jī)器人的運(yùn)動(dòng)參數(shù)等[6-9]。

    筆者將基于粒子濾波器原理,對(duì)已知地圖的輪式移動(dòng)機(jī)器人的自定位展開(kāi)研究。在經(jīng)典粒子濾波算法的基礎(chǔ)上,提出了兩個(gè)改進(jìn)策略,可在保證跟蹤精度的同時(shí),顯著減少計(jì)算量:

    (1)粒子初始化時(shí),每個(gè)粒子的角度經(jīng)過(guò)一次匹配得到,可減少計(jì)算量。具體來(lái)說(shuō),首先估計(jì)每個(gè)粒子在所有可能方向上與周?chē)h(huán)境的距離,然后將這些距離估計(jì)值與移動(dòng)機(jī)器人的測(cè)距傳感器測(cè)量值進(jìn)行匹配,找出最接近該測(cè)量值的角度,作為該粒子的初始角度。這樣不僅能減少自定位所需要的粒子數(shù),而且能減少在計(jì)算移動(dòng)機(jī)器人角度方面的計(jì)算量。

    (2)采用查表式方法估計(jì)距離,可提高定位效率。首先對(duì)環(huán)境地圖進(jìn)行離散化,其后針對(duì)地圖中每個(gè)小格,進(jìn)行360度等間隔離散,并計(jì)算各角度上該小格與周?chē)h(huán)境的距離,將其存為表格形式。當(dāng)機(jī)器人在環(huán)境中行走,并獲得量測(cè)距離時(shí),可以直接調(diào)用上述表格文件,通過(guò)查表方法,判斷各個(gè)粒子在各方向上的概率,從而減少自定位的在線計(jì)算時(shí)間,增強(qiáng)機(jī)器人運(yùn)行及定位的實(shí)時(shí)性。

    1 機(jī)器人運(yùn)動(dòng)模型

    在對(duì)輪式機(jī)器人分析的過(guò)程中,常把機(jī)器人建模為輪子上的一個(gè)剛體,在水平面上運(yùn)動(dòng)。參見(jiàn)圖2,建立平面坐標(biāo)系XOY,則機(jī)器人的位置和姿態(tài)可以用如下?tīng)顟B(tài)向量表示:

    式中:x,y為機(jī)器人在二維直角坐標(biāo)系中的位置坐標(biāo);θ為機(jī)器人的朝向,即機(jī)器人運(yùn)動(dòng)正方向與X軸的夾角。

    圖2 機(jī)器人狀態(tài)示意圖

    假設(shè)在 t1時(shí)刻機(jī)器人處于狀態(tài) S1(x1,y1,θ1),若機(jī)器人左右輪都進(jìn)行勻速運(yùn)動(dòng),經(jīng)過(guò)短暫時(shí)間 Δt之后,機(jī)器人處于狀態(tài) S2(x2,y2,θ2),那么這段時(shí)間內(nèi)機(jī)器人的運(yùn)動(dòng)可近似看作圓周運(yùn)動(dòng),其運(yùn)動(dòng)模型如圖3所示。

    圖3 運(yùn)動(dòng)模型

    則當(dāng)前狀態(tài) S2(x2,y2,θ2)與上一狀態(tài) S1(x1,y1,θ1)之間的關(guān)系為:

    式中:Δsr和Δsl分別為右輪和左輪一次行走的距離;b為機(jī)器人同軸兩個(gè)輪子間的距離;Δs和Δθ分別為機(jī)器人一次運(yùn)動(dòng)行走的距離和方向的變化,這兩個(gè)量即為機(jī)器人在自身局部坐標(biāo)系中的移動(dòng)和轉(zhuǎn)動(dòng),式(2)即為局部坐標(biāo)系與全局坐標(biāo)系之間的轉(zhuǎn)換公式。

    2 移動(dòng)機(jī)器人自定位設(shè)計(jì)

    2.1 移動(dòng)機(jī)器人定位貝葉斯建模

    移動(dòng)機(jī)器人定位原理圖如圖4所示??梢圆捎酶怕蕡D來(lái)描述移動(dòng)機(jī)器人的運(yùn)動(dòng)過(guò)程,其中ut為作用于機(jī)器人的控制信息,pt為t時(shí)刻移動(dòng)機(jī)器人的狀態(tài),ct為機(jī)器人在t時(shí)刻得到的傳感器距離信息。

    圖4 移動(dòng)機(jī)器人定位原理圖

    可以將機(jī)器人定位理解為貝葉斯濾波問(wèn)題。機(jī)器人依據(jù)獲得的量測(cè)信息,以一定的置信度(Belief)確定它的位置:

    式(5)表示在已知0~t時(shí)刻的所有控制數(shù)據(jù)及觀測(cè)數(shù)據(jù) d0…t={ct,ut-1,ct-1,ut-2,…,c0}的情況下,機(jī)器人在t時(shí)刻位于狀態(tài)pt的概率。

    結(jié)合貝葉斯公式變換可得到:

    由于系統(tǒng)具有馬爾科夫性,且 P(ct|ut-1,ct-1,ut-2,…,c0)為一常數(shù),因此式(6)可進(jìn)一步簡(jiǎn)化為:

    通過(guò)對(duì)狀態(tài)空間離散化,則式(7)可以寫(xiě)為:

    式中,E為離散后的狀態(tài)空間。

    2.2 基于粒子濾波器的數(shù)值化實(shí)現(xiàn)

    式(8)給出了離散狀態(tài)空間中的定位模型,由于各種先驗(yàn)概率、似然度函數(shù)往往不具備解析形式,可能為非線性、非高斯的概率,直接由式(8)來(lái)解析求解狀態(tài)是不可能的。

    因此筆者采用經(jīng)典的粒子濾波器來(lái)數(shù)值化求解式(8)。粒子濾波器也稱(chēng)為序貫蒙特卡洛方法(sequential Monte Carlo,SMC),是一種在時(shí)間序列上的離散、非參數(shù)式的貝葉斯濾波方法,能夠處理各種非線性、非高斯、多模態(tài)問(wèn)題。

    2.2.1 狀態(tài)方程及量測(cè)方程

    機(jī)器人的定位問(wèn)題可以描述為:給定一組觀測(cè)向量Ct,估計(jì)機(jī)器人當(dāng)前所在狀態(tài)pt。其中:Ct={c0,c1,…,ct}為觀測(cè)序列;pt=[xt,yt,θt]T為機(jī)器人在t時(shí)刻的位置和朝向。

    設(shè)移動(dòng)機(jī)器人狀態(tài)方程為:

    式中:ut為機(jī)器人的運(yùn)動(dòng)控制輸入變量;wt為過(guò)程噪聲。函數(shù)f(x)由機(jī)器人的運(yùn)動(dòng)模型給定,即在式(2)的基礎(chǔ)上加上噪聲wt,則可得:

    式中,ws和wθ分別為機(jī)器人的運(yùn)動(dòng)距離和偏轉(zhuǎn)角度的噪聲。

    系統(tǒng)狀態(tài)量測(cè)方程為:

    式中:M為環(huán)境地圖;vt為量測(cè)噪聲;pt的估計(jì)由求解Bel(pt)得到。

    2.2.2 初始化粒子

    移動(dòng)機(jī)器人自定位程序首先是粒子的初始化。初始的粒子服從均勻分布,即每個(gè)粒子出現(xiàn)在圖像中任意一個(gè)點(diǎn)的概率都是1/N,其中N為移動(dòng)機(jī)器人在地圖上可能出現(xiàn)的位置總數(shù)。

    由式(1)可知,粒子的初始化除了要初始化粒子的位置之外,還要初始化粒子的朝向角度,即θ。因此需要大量的粒子來(lái)表示移動(dòng)機(jī)器人所有可能的位置和角度,計(jì)算量非常大。

    假定θ值按均勻分布生成粒子,即θ~U(0,2π),那么當(dāng)某個(gè)粒子的坐標(biāo)與機(jī)器人的真實(shí)2D位置吻合時(shí),該粒子隨機(jī)生成的θ與機(jī)器人真實(shí)的方向吻合的概率并不大,從而造成計(jì)算浪費(fèi)。

    實(shí)際上,在已獲得距離量測(cè)的情況下,狀態(tài)向量中角度分量θ是不獨(dú)立的,其與其他兩個(gè)分量,即機(jī)器人在二維空間的坐標(biāo)x、y是相關(guān)的。因此可以把狀態(tài)變量分為兩個(gè)部分,則有式(12):

    根據(jù)式(12),首先基于均勻分布p(xt,yt)~U(M)對(duì)二維坐標(biāo)x,y進(jìn)行采樣(M為環(huán)境地圖),然后對(duì)于每一個(gè)粒子,求出該粒子在自己位置上匹配權(quán)重最大的角度,將該值作為θ的初始值,這與文獻(xiàn)[10]提出的方法有些類(lèi)似,都是把狀態(tài)向量分為幾個(gè)部分,將聯(lián)合概率分解為條件概率的乘積形式。而筆者的創(chuàng)新之處在于,結(jié)合坐標(biāo)分量和量測(cè)來(lái)直接估計(jì)每個(gè)粒子的角度分量。其優(yōu)點(diǎn)在于既能有效提高粒子的采樣質(zhì)量,又能極大地減少用于確定機(jī)器人朝向的粒子數(shù)量,并且在以后的定位計(jì)算中也不需要再考慮移動(dòng)機(jī)器人的θ值問(wèn)題,只需根據(jù)式(10)將機(jī)器人運(yùn)動(dòng)產(chǎn)生的θ值的變化量Δθ加到每個(gè)粒子的θ值即可,不用再計(jì)算移動(dòng)機(jī)器人的θ值。

    2.2.3 粒子重采樣

    在完成粒子的初始化之后,每個(gè)粒子的權(quán)重也在初始化θ值時(shí)賦給了每個(gè)粒子。此時(shí)就可以找出權(quán)重最大的粒子作為本次機(jī)器人位置的最佳估計(jì)。下一步的工作就是對(duì)粒子進(jìn)行重采樣,為下一次估計(jì)機(jī)器人位置做準(zhǔn)備。

    重采樣的目的是去掉權(quán)重小的粒子,同時(shí)將權(quán)重大的粒子分解成若干個(gè)權(quán)重較小的新粒子。重采樣的思想可以用圖5來(lái)表示。

    圖5 重采樣示意圖

    圖5(a)為粒子權(quán)重的累計(jì)概率分布函數(shù),橫坐標(biāo)為粒子的序號(hào);圖5(b)為0~1間均勻分布的隨機(jī)變量,u為該隨機(jī)變量的一次采樣值重采樣過(guò)程中,首先計(jì)算出粒子概率的累積概率,然后在0~1間按均勻分布隨機(jī)采樣,投影到圖5(a)中,投影的位置向下指示出累積概率相應(yīng)定義域的位置,這就是新粒子的索引。這樣共產(chǎn)生L個(gè)粒子,其權(quán)重都重置為1/L。顯然,權(quán)重大的粒子就可能會(huì)被多次重復(fù)采樣。

    2.2.4 粒子狀態(tài)預(yù)測(cè)

    完成了粒子的重采樣之后,機(jī)器人運(yùn)動(dòng)一次。根據(jù)里程計(jì)的返回值由式(3)和式(4)計(jì)算出這一次運(yùn)動(dòng)后,機(jī)器人相對(duì)運(yùn)動(dòng)前的位置和角度的變化,并由該變化量按照式(10)系統(tǒng)狀態(tài)方程,更新每個(gè)粒子的位置和角度,使每個(gè)粒子都能相對(duì)于自己上一時(shí)刻的位置和角度,作與移動(dòng)機(jī)器人相同的運(yùn)動(dòng)。

    如機(jī)器人向前走了10 cm,每個(gè)粒子也按自己的方向向前走10 cm。實(shí)際環(huán)境中由于受車(chē)輪打滑、里程計(jì)測(cè)量有誤差等因素的影響,雖然里程計(jì)返回值是10 cm,但機(jī)器人很有可能不是走了剛好10 cm。因此式(10)引入了過(guò)程噪聲wt。

    2.2.5 粒子權(quán)值更新

    在所有粒子都隨移動(dòng)機(jī)器人運(yùn)動(dòng)后,下一步的工作就是為所有粒子重新計(jì)算權(quán)值。

    基于環(huán)境柵格地圖,對(duì)于每一個(gè)粒子,根據(jù)該粒子的狀態(tài)(包含位置及方向),計(jì)算出該粒子在該方向上與環(huán)境的距離d,并將距離值與移動(dòng)機(jī)器人測(cè)距傳感器的真實(shí)返回值d^相比較,從而求出式(8)中的P(ct|pt):

    式中:C為歸一化系數(shù);∑為由傳感器量測(cè)噪聲vt的方差。

    由于粒子的數(shù)目非常多,如果用實(shí)時(shí)計(jì)算的方式來(lái)求解每個(gè)粒子的測(cè)量值會(huì)需要很長(zhǎng)的時(shí)間。因此,設(shè)計(jì)了一種基于查表方式的權(quán)值更新方法。首先在離線形式下,把機(jī)器人在地圖中所有可能的位置和姿態(tài)列出來(lái),并計(jì)算出相應(yīng)各傳感器的距離信息,生成一個(gè)二進(jìn)制文件。

    通過(guò)查表方法,可以實(shí)時(shí)地將需要的數(shù)據(jù)d^從文件中找出來(lái),由式(13)計(jì)算量測(cè)概率,從而減少了機(jī)器人自定位的計(jì)算量,縮短了機(jī)器人兩次運(yùn)動(dòng)之間的時(shí)間,增加了自定位程序的實(shí)時(shí)性,但也一定程度地增加了機(jī)器人內(nèi)存的消耗。

    基于粒子濾波理論的移動(dòng)機(jī)器人自定位程序的流程圖如圖6所示。

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

    圖6 自定位流程圖

    筆者使用VC2005對(duì)算法進(jìn)行仿真。圖7中使用的粒子數(shù)為3500,灰色區(qū)域?yàn)檎系K物和墻壁,灰色的點(diǎn)為粒子。實(shí)線為仿真程序生成的機(jī)器人實(shí)際位置,機(jī)器人前方兩條線的長(zhǎng)度表示兩個(gè)測(cè)距傳感器的實(shí)際測(cè)量值;虛線為機(jī)器人最優(yōu)估計(jì)粒子的位置和方向,兩條射線的長(zhǎng)度表示在該位置和方向上,兩個(gè)測(cè)距傳感器對(duì)應(yīng)的測(cè)量估計(jì)值,可以通過(guò)查表獲得。

    仿真使用的地圖是由4個(gè)房間和1條走廊以及每個(gè)房間內(nèi)的障礙物組成簡(jiǎn)單結(jié)構(gòu)的環(huán)境。由于仿真環(huán)境地圖比較簡(jiǎn)單,環(huán)境中存在對(duì)稱(chēng)的點(diǎn),則機(jī)器人在自定位時(shí),粒子會(huì)集中到這些與機(jī)器人所在真實(shí)環(huán)境對(duì)稱(chēng)的點(diǎn)附近,只有當(dāng)機(jī)器人運(yùn)動(dòng)到獨(dú)特的位置時(shí),所有的粒子才會(huì)聚集到機(jī)器人的真實(shí)位置附近。

    圖7(b)為機(jī)器人仿真定位程序的第1次定位結(jié)果。與圖7(a)的初始狀態(tài)相比,粒子已經(jīng)開(kāi)始聚集,不過(guò)由于信息還比較少,因此最優(yōu)的粒子與實(shí)際的機(jī)器人位于不同房間。

    圖7(c)為仿真程序的第2次定位結(jié)果,與第1次定位的結(jié)果相比,粒子又進(jìn)一步集中,且最優(yōu)估計(jì)狀態(tài)也正確估計(jì)出機(jī)器人所在房間。

    隨著機(jī)器人的運(yùn)動(dòng),距離量測(cè)信息越來(lái)越多,第4次定位后,效果比較明顯,參見(jiàn)圖7(e),到第6次定位完成后,估計(jì)就非常準(zhǔn)確了,參見(jiàn)圖7(f)。

    將筆者提出的算法與經(jīng)典蒙特卡洛算法進(jìn)行比較可知,筆者采用了角度計(jì)算策略,而經(jīng)典算法直接對(duì)式(1)狀態(tài)進(jìn)行采樣;筆者采用了查表計(jì)算權(quán)值,而經(jīng)典算法是采用在線計(jì)算粒子權(quán)重的方法。

    采用的仿真實(shí)驗(yàn)粒子數(shù)為3500,運(yùn)行環(huán)境為Windows Xp下 VS2005,奔騰雙核 2.6 G,1 G 內(nèi)存。比較結(jié)果如表1所示,可以發(fā)現(xiàn)筆者算法在耗時(shí)及收斂速度上具有優(yōu)越性。

    圖7 筆者算法定位仿真

    表1 筆者算法與經(jīng)典蒙特卡洛算法比較

    4 結(jié)論

    筆者針對(duì)經(jīng)典蒙特卡洛定位算法,提出了兩種改進(jìn)策略。首先粒子初始化時(shí),每個(gè)粒子的角度經(jīng)過(guò)一次匹配得到,從而減小了計(jì)算量。其次,采用離線方式計(jì)算場(chǎng)景的權(quán)重表,在定位過(guò)程中,經(jīng)過(guò)直接查表,可極大地減少定位的時(shí)間消耗。

    通過(guò)與經(jīng)典蒙特卡洛定位算法的比較,筆者的算法具有計(jì)算量小,收斂速度快的優(yōu)點(diǎn)。

    [1]蔡自興.機(jī)器人學(xué)[M].北京:清華大學(xué)出版社,2000:16-25.

    [2]IEGWART R,NOURBAKHSH I R.自主移動(dòng)機(jī)器人導(dǎo)論[M].李仁厚,譯.西安:西安交通大學(xué)出版社,2006:21-120.

    [3]DURRANT-WHYTE H,BAILEY T.Simultaneous localization and mapping:part I[J].Robotics & Automation Magazine,IEEE,2006,13(2):99 -110.

    [4]BAILEY T,DURRANT -WHYTE H.Simultaneous localization and mapping(SLAM):part II[J].Robotics& Automation Magazine,IEEE,2006,13(3):108-117.

    [5]THURN S,BEETZ M,BENNEWITZ M,et al.Probabilistic algorithms and the interactive museum tourguide robot minerva[J].International Journal of Robotics Research,2000,19(11):972 -999.

    [6]方正,佟國(guó)峰,徐心和.基于貝葉斯濾波理論的自助機(jī)器人自定位方法研究[J].控制與決策,2006,21(8):841-862.

    [7]曾慧,吳福朝,胡占義.基于平面激光測(cè)量的移動(dòng)機(jī)器人自定位方法[J].自動(dòng)化學(xué)報(bào),2007,33(2):138 -144.

    [8]劉承俊,原魁,鄒偉,等.基于特征粒子的Monte Carlo自定位方法[J].機(jī)器人,2006,28(1):30 -35.

    [9]陳衛(wèi)東,張飛.移動(dòng)機(jī)器人的同步自定位與地圖創(chuàng)建研究進(jìn)展[J].控制理論與應(yīng)用,2005,22(3):455 -460.

    [10]OLIVIER C,SIMON J G,ERIC M.An overview of existing methods and recent advances in sequential Monte Carlo[J].IEEE,2007,95(5):654 -662.

    猜你喜歡
    移動(dòng)機(jī)器人粒子角度
    移動(dòng)機(jī)器人自主動(dòng)態(tài)避障方法
    神奇的角度
    一個(gè)涉及角度和的幾何不等式鏈的改進(jìn)
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    角度不同
    37°女人(2017年8期)2017-08-12 11:20:48
    人啊
    滇池(2017年7期)2017-07-18 19:32:42
    基于Twincat的移動(dòng)機(jī)器人制孔系統(tǒng)
    極坐標(biāo)系下移動(dòng)機(jī)器人的點(diǎn)鎮(zhèn)定
    基于引導(dǎo)角的非完整移動(dòng)機(jī)器人軌跡跟蹤控制
    午夜精品一区二区三区免费看| 亚洲内射少妇av| 成人三级黄色视频| 国产亚洲5aaaaa淫片| 亚洲中文字幕一区二区三区有码在线看| 老熟妇乱子伦视频在线观看| av在线亚洲专区| 精品国内亚洲2022精品成人| 啦啦啦韩国在线观看视频| 国产伦精品一区二区三区四那| 日本在线视频免费播放| 欧洲精品卡2卡3卡4卡5卡区| 一级av片app| 最近的中文字幕免费完整| 99国产精品一区二区蜜桃av| 深爱激情五月婷婷| 国产一区二区亚洲精品在线观看| 91麻豆精品激情在线观看国产| 精品久久久久久久久亚洲| 亚洲精品乱码久久久久久按摩| 国产一级毛片七仙女欲春2| 亚洲经典国产精华液单| 大香蕉久久网| 成人亚洲欧美一区二区av| 亚洲经典国产精华液单| 网址你懂的国产日韩在线| 欧美人与善性xxx| 99精品在免费线老司机午夜| 国产色婷婷99| 国产精品电影一区二区三区| 人人妻人人看人人澡| 99久久精品一区二区三区| 国产欧美日韩精品一区二区| 少妇猛男粗大的猛烈进出视频 | 亚洲成人精品中文字幕电影| 男人舔奶头视频| 亚洲精品亚洲一区二区| 久久人人精品亚洲av| 国产三级在线视频| 精品久久久久久久久亚洲| 2021天堂中文幕一二区在线观| 亚洲精品日韩在线中文字幕 | av在线播放精品| 欧美日韩乱码在线| 亚洲七黄色美女视频| 毛片一级片免费看久久久久| 91久久精品国产一区二区成人| 欧美一区二区亚洲| 综合色丁香网| 免费大片18禁| 日本av手机在线免费观看| 99九九线精品视频在线观看视频| 有码 亚洲区| 久久精品综合一区二区三区| 一级毛片久久久久久久久女| 日韩成人av中文字幕在线观看| 日本av手机在线免费观看| 亚洲av熟女| 一边摸一边抽搐一进一小说| 给我免费播放毛片高清在线观看| 久久精品影院6| 国产精品国产三级国产av玫瑰| 国产高清视频在线观看网站| 亚州av有码| 亚洲国产欧美在线一区| 男女边吃奶边做爰视频| a级毛片a级免费在线| 亚洲av一区综合| 夜夜爽天天搞| 精品一区二区三区视频在线| 久久99热6这里只有精品| 99久国产av精品| 亚洲国产精品久久男人天堂| 又爽又黄a免费视频| 亚洲av成人精品一区久久| 国产伦精品一区二区三区四那| 国产精品不卡视频一区二区| 日韩国内少妇激情av| 亚洲性久久影院| 日本在线视频免费播放| 中文在线观看免费www的网站| 午夜福利成人在线免费观看| 中文字幕制服av| 日韩欧美国产在线观看| 国内精品久久久久精免费| 精品久久久久久久人妻蜜臀av| 伦精品一区二区三区| 中文资源天堂在线| 国产单亲对白刺激| 久久韩国三级中文字幕| 亚洲av.av天堂| 国产爱豆传媒在线观看| 18+在线观看网站| 特大巨黑吊av在线直播| 亚洲精品久久国产高清桃花| 欧美精品国产亚洲| 蜜臀久久99精品久久宅男| 亚洲欧美成人精品一区二区| 免费看av在线观看网站| 一个人观看的视频www高清免费观看| 在线免费观看的www视频| 精品一区二区三区人妻视频| 直男gayav资源| 久久久久九九精品影院| 国产黄a三级三级三级人| 亚洲精品日韩在线中文字幕 | 午夜福利在线观看免费完整高清在 | 中文字幕精品亚洲无线码一区| 人妻少妇偷人精品九色| 国产麻豆成人av免费视频| 国产精品美女特级片免费视频播放器| 亚洲一区高清亚洲精品| 直男gayav资源| 久久久久久久亚洲中文字幕| 欧美高清成人免费视频www| 免费看光身美女| 亚洲电影在线观看av| 亚洲激情五月婷婷啪啪| 男女下面进入的视频免费午夜| 精品久久久久久久久av| 在线观看美女被高潮喷水网站| 少妇猛男粗大的猛烈进出视频 | 国产真实乱freesex| 婷婷色av中文字幕| 精品久久久久久久久久久久久| 99精品在免费线老司机午夜| kizo精华| 国产毛片a区久久久久| 在线观看免费视频日本深夜| 国产女主播在线喷水免费视频网站 | 丰满乱子伦码专区| 国产高清有码在线观看视频| 99久国产av精品国产电影| 成年av动漫网址| 国产高清激情床上av| 精华霜和精华液先用哪个| 91麻豆精品激情在线观看国产| a级毛片免费高清观看在线播放| 免费观看人在逋| 亚洲婷婷狠狠爱综合网| 一进一出抽搐动态| 又爽又黄无遮挡网站| 级片在线观看| 中文字幕av在线有码专区| 国产毛片a区久久久久| 国产精品久久久久久久电影| 99久久久亚洲精品蜜臀av| 精品人妻视频免费看| 日韩精品有码人妻一区| 国产高清视频在线观看网站| 深爱激情五月婷婷| 久久久久九九精品影院| 久久草成人影院| 成人一区二区视频在线观看| av天堂中文字幕网| 精品人妻偷拍中文字幕| 3wmmmm亚洲av在线观看| 一边亲一边摸免费视频| 欧美性感艳星| 一区二区三区免费毛片| 国产美女午夜福利| 亚洲人成网站在线播| 在线国产一区二区在线| h日本视频在线播放| 国产三级中文精品| 丰满的人妻完整版| 九九热线精品视视频播放| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美三级三区| 国产成人福利小说| 六月丁香七月| 美女 人体艺术 gogo| 国产精品永久免费网站| 亚洲性久久影院| 亚洲成人久久性| www日本黄色视频网| 国产一区二区亚洲精品在线观看| 亚洲精品久久国产高清桃花| 日本黄色片子视频| 亚洲在线自拍视频| 老熟妇乱子伦视频在线观看| 国产极品天堂在线| 精品一区二区免费观看| 99久国产av精品国产电影| 国产日本99.免费观看| 观看免费一级毛片| 欧美人与善性xxx| 小说图片视频综合网站| 亚洲精品乱码久久久久久按摩| 色哟哟哟哟哟哟| 国产精品国产高清国产av| 午夜福利高清视频| 丰满的人妻完整版| 69人妻影院| 久久精品夜色国产| 成人综合一区亚洲| 免费无遮挡裸体视频| 日韩人妻高清精品专区| 美女黄网站色视频| 日日干狠狠操夜夜爽| 中文字幕精品亚洲无线码一区| 老女人水多毛片| 少妇熟女aⅴ在线视频| 久久久午夜欧美精品| 深夜精品福利| 欧美激情在线99| 三级毛片av免费| 国产大屁股一区二区在线视频| 国产伦精品一区二区三区四那| 欧美日韩综合久久久久久| 国产三级在线视频| 亚洲av成人av| 亚洲五月天丁香| 日韩一区二区三区影片| 国产日韩欧美在线精品| 国产麻豆成人av免费视频| 美女黄网站色视频| 国产在线男女| 天美传媒精品一区二区| 一个人看的www免费观看视频| 舔av片在线| 丰满人妻一区二区三区视频av| 淫秽高清视频在线观看| 国产亚洲精品久久久久久毛片| 此物有八面人人有两片| 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 久久精品国产鲁丝片午夜精品| 免费看光身美女| 22中文网久久字幕| 久久人人爽人人爽人人片va| 国产成人影院久久av| 国产精品伦人一区二区| 在线国产一区二区在线| 国产一级毛片七仙女欲春2| 男人舔女人下体高潮全视频| 淫秽高清视频在线观看| 简卡轻食公司| 在线播放国产精品三级| 亚洲国产精品合色在线| av国产免费在线观看| a级毛片a级免费在线| 亚洲欧美中文字幕日韩二区| 男女那种视频在线观看| 不卡视频在线观看欧美| 久久欧美精品欧美久久欧美| 国产高清视频在线观看网站| 女人十人毛片免费观看3o分钟| 亚洲国产精品成人久久小说 | 日韩一本色道免费dvd| 国产蜜桃级精品一区二区三区| 国产亚洲91精品色在线| 亚洲高清免费不卡视频| 久久精品国产亚洲av天美| 久久久久久久久久成人| 国内少妇人妻偷人精品xxx网站| 久久精品国产清高在天天线| 天天躁日日操中文字幕| 国产一级毛片在线| 人人妻人人看人人澡| 欧美精品一区二区大全| 国产私拍福利视频在线观看| 一区二区三区四区激情视频 | 成人毛片a级毛片在线播放| 日本黄大片高清| 亚洲性久久影院| 中文字幕av成人在线电影| 长腿黑丝高跟| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 国产精品国产高清国产av| 自拍偷自拍亚洲精品老妇| 国产av不卡久久| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 中国美白少妇内射xxxbb| 2022亚洲国产成人精品| АⅤ资源中文在线天堂| 午夜精品在线福利| 国内精品久久久久精免费| 边亲边吃奶的免费视频| 国产成年人精品一区二区| 日本欧美国产在线视频| 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 国产私拍福利视频在线观看| 最后的刺客免费高清国语| 久久久a久久爽久久v久久| 精品不卡国产一区二区三区| 免费观看a级毛片全部| 99精品在免费线老司机午夜| 免费av观看视频| 亚洲乱码一区二区免费版| 国产精品爽爽va在线观看网站| 午夜激情福利司机影院| 亚洲人成网站在线观看播放| 午夜福利高清视频| 丰满乱子伦码专区| 免费在线观看成人毛片| 麻豆一二三区av精品| 男女啪啪激烈高潮av片| 丰满乱子伦码专区| 精品无人区乱码1区二区| 国产伦在线观看视频一区| 国产精品无大码| 国产高清有码在线观看视频| 蜜桃久久精品国产亚洲av| 成人特级黄色片久久久久久久| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 亚洲在线自拍视频| 欧美一区二区精品小视频在线| 淫秽高清视频在线观看| 久久久久性生活片| 日韩成人伦理影院| 日本色播在线视频| 亚洲国产欧美在线一区| 国产麻豆成人av免费视频| 国产探花极品一区二区| www.色视频.com| 亚洲成人久久爱视频| 国产精品乱码一区二三区的特点| 午夜精品一区二区三区免费看| 在线天堂最新版资源| 国产高清激情床上av| 国产精品美女特级片免费视频播放器| 免费av不卡在线播放| 日韩强制内射视频| 日本三级黄在线观看| 此物有八面人人有两片| 欧美性感艳星| 久久久久久久亚洲中文字幕| 国产午夜精品久久久久久一区二区三区| 人妻制服诱惑在线中文字幕| 中文在线观看免费www的网站| 99久久精品一区二区三区| 国产老妇女一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人妻久久中文字幕网| 亚洲av二区三区四区| 我的女老师完整版在线观看| 毛片女人毛片| 联通29元200g的流量卡| 国产精品综合久久久久久久免费| 中国国产av一级| 亚洲精品456在线播放app| av福利片在线观看| 国产男人的电影天堂91| 国产亚洲精品久久久com| 国产精品麻豆人妻色哟哟久久 | 一级毛片我不卡| 亚洲av成人精品一区久久| 22中文网久久字幕| 日本爱情动作片www.在线观看| 亚洲精品乱码久久久v下载方式| 不卡一级毛片| 久久久久久久午夜电影| 成人鲁丝片一二三区免费| 亚洲在久久综合| 国产精品免费一区二区三区在线| 春色校园在线视频观看| 日韩大尺度精品在线看网址| 成人毛片a级毛片在线播放| 欧美日韩乱码在线| 午夜福利成人在线免费观看| 国产老妇女一区| 国产一区二区在线av高清观看| 性欧美人与动物交配| 国产高清不卡午夜福利| 久久精品影院6| 久久精品久久久久久久性| 国产黄片美女视频| 亚洲国产欧洲综合997久久,| 蜜桃亚洲精品一区二区三区| 人人妻人人澡人人爽人人夜夜 | 久久99热6这里只有精品| 三级经典国产精品| 变态另类丝袜制服| 美女高潮的动态| 夫妻性生交免费视频一级片| 久久久久久久久久黄片| 热99re8久久精品国产| 18+在线观看网站| 亚洲精品久久久久久婷婷小说 | 久久综合国产亚洲精品| 又粗又硬又长又爽又黄的视频 | 欧美色视频一区免费| 大又大粗又爽又黄少妇毛片口| 亚洲成人久久性| 少妇的逼好多水| 能在线免费观看的黄片| 午夜视频国产福利| 精品人妻视频免费看| 综合色丁香网| 老司机福利观看| 欧美最新免费一区二区三区| 精品午夜福利在线看| 在现免费观看毛片| 欧美潮喷喷水| 此物有八面人人有两片| 免费大片18禁| 校园春色视频在线观看| 高清午夜精品一区二区三区 | 老师上课跳d突然被开到最大视频| 亚洲成人av在线免费| 国产精品一区二区三区四区久久| 成人性生交大片免费视频hd| 欧美日本视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美精品一区二区大全| 乱人视频在线观看| 99久国产av精品| 国产熟女欧美一区二区| 国产一区二区在线av高清观看| 色视频www国产| 国产精品一区二区在线观看99 | АⅤ资源中文在线天堂| 在线免费观看不下载黄p国产| 久久精品国产亚洲av涩爱 | 久久精品久久久久久久性| 国产精品.久久久| 亚洲自偷自拍三级| 国产精品一二三区在线看| 成人高潮视频无遮挡免费网站| 国产不卡一卡二| 亚洲欧美清纯卡通| 亚洲成av人片在线播放无| 日韩中字成人| 免费观看的影片在线观看| 国产激情偷乱视频一区二区| 久久人人精品亚洲av| 两个人视频免费观看高清| 久久人妻av系列| 亚洲激情五月婷婷啪啪| 欧美zozozo另类| 精品久久久久久久久久免费视频| 久久久久久九九精品二区国产| 超碰av人人做人人爽久久| 欧美精品国产亚洲| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 在线观看66精品国产| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 91久久精品电影网| 嫩草影院新地址| 久久人妻av系列| 欧美日本亚洲视频在线播放| 欧美成人一区二区免费高清观看| 国内精品久久久久精免费| 精品久久久久久久久亚洲| av女优亚洲男人天堂| 少妇高潮的动态图| 精品久久久久久久久av| 亚洲av电影不卡..在线观看| 午夜福利在线观看免费完整高清在 | 国产一区二区在线av高清观看| 国产精品久久电影中文字幕| 亚洲成人久久性| 搡老妇女老女人老熟妇| 亚洲在线自拍视频| 国产真实伦视频高清在线观看| 97热精品久久久久久| 国产精品野战在线观看| 日本色播在线视频| 久久久色成人| 欧美人与善性xxx| 国产精品免费一区二区三区在线| 国产美女午夜福利| 成年版毛片免费区| 一级毛片电影观看 | 黄片无遮挡物在线观看| 九九热线精品视视频播放| 婷婷色av中文字幕| 女同久久另类99精品国产91| 日韩强制内射视频| 尤物成人国产欧美一区二区三区| 亚洲一区二区三区色噜噜| 国产一区二区激情短视频| 欧美成人免费av一区二区三区| 国产成人精品久久久久久| 精品免费久久久久久久清纯| 亚洲国产精品国产精品| 男女啪啪激烈高潮av片| 69人妻影院| 国产精品电影一区二区三区| 看非洲黑人一级黄片| 最近的中文字幕免费完整| 人人妻人人澡人人爽人人夜夜 | 舔av片在线| 级片在线观看| 中文欧美无线码| 亚洲精品久久国产高清桃花| 久久久欧美国产精品| 99久久精品热视频| 麻豆精品久久久久久蜜桃| 亚洲av二区三区四区| 欧美性猛交╳xxx乱大交人| 久久久久久久亚洲中文字幕| 少妇人妻一区二区三区视频| 给我免费播放毛片高清在线观看| 久久久欧美国产精品| 免费av毛片视频| 晚上一个人看的免费电影| 大型黄色视频在线免费观看| 一个人观看的视频www高清免费观看| 免费一级毛片在线播放高清视频| 国产精品一及| 蜜桃久久精品国产亚洲av| 少妇人妻一区二区三区视频| 热99re8久久精品国产| 边亲边吃奶的免费视频| 亚洲av免费高清在线观看| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放| 日韩成人伦理影院| 亚洲国产欧美人成| 国产色婷婷99| 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 2021天堂中文幕一二区在线观| 日韩 亚洲 欧美在线| 成人毛片60女人毛片免费| 六月丁香七月| 白带黄色成豆腐渣| 一本一本综合久久| 91精品一卡2卡3卡4卡| 久久久久久国产a免费观看| 又黄又爽又刺激的免费视频.| 搡女人真爽免费视频火全软件| 欧美丝袜亚洲另类| 成人毛片60女人毛片免费| 日本一二三区视频观看| 国产黄片美女视频| 我要看日韩黄色一级片| 午夜视频国产福利| 亚洲av成人av| 男女那种视频在线观看| 久久久久久伊人网av| 天堂影院成人在线观看| 国产精品av视频在线免费观看| 国产精品蜜桃在线观看 | a级毛片a级免费在线| 成人一区二区视频在线观看| 黑人高潮一二区| 男女那种视频在线观看| 国产精品一及| 蜜桃久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 丝袜美腿在线中文| 日本熟妇午夜| 91在线精品国自产拍蜜月| 国产人妻一区二区三区在| 国内精品美女久久久久久| 午夜免费激情av| 日本黄色视频三级网站网址| 在线观看午夜福利视频| 欧美区成人在线视频| videossex国产| 网址你懂的国产日韩在线| 日韩av在线大香蕉| 国产黄片美女视频| 18禁在线无遮挡免费观看视频| 亚洲经典国产精华液单| 嫩草影院入口| 一个人看的www免费观看视频| 女人十人毛片免费观看3o分钟| 黄色一级大片看看| 久久6这里有精品| 看十八女毛片水多多多| 六月丁香七月| 99视频精品全部免费 在线| 91久久精品国产一区二区成人| 老熟妇乱子伦视频在线观看| av卡一久久| 国产老妇女一区| 99久久精品国产国产毛片| 精品不卡国产一区二区三区| 亚洲五月天丁香| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 中文亚洲av片在线观看爽| 亚洲四区av| 男的添女的下面高潮视频| 美女 人体艺术 gogo| 国产av在哪里看| 久久精品国产99精品国产亚洲性色| 亚洲激情五月婷婷啪啪| 午夜激情福利司机影院| 99国产极品粉嫩在线观看| 国产单亲对白刺激| 精华霜和精华液先用哪个| 日韩欧美三级三区| 老女人水多毛片| 久久人人爽人人爽人人片va| 欧美变态另类bdsm刘玥| 乱人视频在线观看| 日韩国内少妇激情av| av视频在线观看入口| av福利片在线观看| 最近手机中文字幕大全| av在线蜜桃| 国产中年淑女户外野战色| 国产精品无大码| 免费黄网站久久成人精品| 久久久久久伊人网av| 又黄又爽又刺激的免费视频.| 亚洲av成人av| 免费av观看视频| 一夜夜www| 国产精品.久久久| 禁无遮挡网站| 我的老师免费观看完整版|