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

    基于自適應(yīng)聯(lián)邦濾波的AGV定位研究

    2022-10-23 14:16:38溫尊旺王堯堯陳柏姚佳烽
    機(jī)械制造與自動(dòng)化 2022年5期
    關(guān)鍵詞:單目里程計(jì)卡爾曼濾波

    溫尊旺,王堯堯,陳柏,姚佳烽

    (南京航空航天大學(xué) 機(jī)電學(xué)院,江蘇 南京 210016)

    0 引言

    隨著“工業(yè)4.0”、“中國制造2025”等戰(zhàn)略的提出,傳統(tǒng)工業(yè)逐漸向自動(dòng)化及智能化的方向發(fā)展。由于自動(dòng)引導(dǎo)運(yùn)輸車(AGV)智能化高,靈活性強(qiáng)的特點(diǎn),廣泛應(yīng)用于智能車間、倉儲(chǔ)運(yùn)輸?shù)阮I(lǐng)域[1]。常見的AGV室內(nèi)定位方式主要有WiFi定位[2]、激光雷達(dá)、UWB定位[3]、視覺[4]、動(dòng)作捕捉、電磁導(dǎo)航等。這些定位方式各有特點(diǎn),在不同場合均有使用。

    為了進(jìn)一步提高定位精度,便有了多傳感器的組合導(dǎo)航,即采用兩種或兩種以上的非相似性導(dǎo)航系統(tǒng)對(duì)同一導(dǎo)航信息做測量并解算以形成量測量[5]。COONEY J A等[6]介紹了麥克納姆輪(Mecanum wheel)中的打滑問題,使用視覺航位推算來處理和計(jì)算車輛的速度及方向,以此來完成防滑的功能。QIN H L等[7]提出了基于小波閾值去噪和函數(shù)系數(shù)回歸(FAR)建模的慣性組合數(shù)據(jù)處理方法,引入GPS姿態(tài)信息,提高了MEMS慣性傳感器誤差的估計(jì)精度。孫自飛等[8]提出一種RGB-D傳感器對(duì)激光傳感器進(jìn)行輔助的自定位方法,具有較好的準(zhǔn)確性和可靠性。李沛等[9]提出了一種基于視覺定位的AGV避障方案,將超聲波傳感器和紅外傳感器相結(jié)合,基于多傳感器信息融合技術(shù),實(shí)現(xiàn)避障算法。

    本文基于Mecanum輪全向機(jī)器人設(shè)計(jì)了一種多傳感器融合的低成本定位方案。結(jié)合各傳感器的特點(diǎn),采用自適應(yīng)聯(lián)邦卡爾曼濾波及多尺度融合的手段,提高了定位精度。

    1 Mecanum輪全向機(jī)器人運(yùn)動(dòng)學(xué)分析

    所設(shè)計(jì)的Mecanum輪全向機(jī)器人系統(tǒng)包含動(dòng)力系統(tǒng)、傳感器系統(tǒng)以及控制系統(tǒng)。動(dòng)力系統(tǒng)包括電機(jī)驅(qū)動(dòng)器、步進(jìn)電機(jī)以及麥克納姆輪;傳感器系統(tǒng)包括編碼器、加速度計(jì)、陀螺儀以及OpenMV,所用的傳感器成本都較低??刂葡到y(tǒng)由STM32 F103組成,主要用于傳感器數(shù)據(jù)讀取、解算。

    Mecanum輪全向機(jī)器人的安裝方式主要為X型和O型,本文采用X型方式,如圖1所示[10]。R為Mecanum輪半徑;a、b分別為平臺(tái)寬度、長度的一半;vx、vy、ωz分別為平臺(tái)的x方向速度、y方向的速度以及角速度;v1、v2、v3、v4分別為4個(gè)麥克納姆輪的分速度;α為輪轂軸和輥?zhàn)虞S的夾角。

    圖1 Mecanum輪全向機(jī)器人結(jié)構(gòu)圖

    Mecanum輪全向機(jī)器人的運(yùn)動(dòng)可以分解為3個(gè)獨(dú)立的變量,分別為:x軸的平移、y軸的平移以及繞z軸的旋轉(zhuǎn)。如圖1所示建立坐標(biāo)系,假設(shè)小車的運(yùn)動(dòng)速度為V=[vx,vy,ωz],對(duì)小車的運(yùn)動(dòng)速度進(jìn)行運(yùn)動(dòng)分解,根據(jù)疊加原理,各個(gè)輪子的角速度為

    (1)

    式中:kx在理想狀態(tài)=1/R;ky=1/(Rtanα);kz=(atanα+b)/(Rtanα)。

    2 傳感器定位原理

    為了獲取準(zhǔn)確的位姿,采用多種傳感器實(shí)時(shí)獲得與位姿相關(guān)的信息,包括:OpenMV單目相機(jī)通過AprilTag標(biāo)簽獲得x、y方向的位置、z軸的姿態(tài);通過4個(gè)電機(jī)的編碼器獲得當(dāng)前的速度和角速度;陀螺儀得到z軸角度及角速度;MPU6050獲得x、y方向的加速度。

    2.1 坐標(biāo)系規(guī)定及坐標(biāo)變換

    為了方便后續(xù)的計(jì)算,對(duì)各個(gè)傳感器的坐標(biāo)系進(jìn)行規(guī)定,均采用右手定則,如圖2所示。其中包括世界坐標(biāo)系、AprilTag坐標(biāo)系、OpenMv(單目相機(jī))坐標(biāo)系、AGV坐標(biāo)系。相對(duì)于世界坐標(biāo)系,AprilTag坐標(biāo)系是不變的,其余坐標(biāo)系都隨著AGV坐標(biāo)系而變化,假設(shè)初始位置AGV坐標(biāo)系和世界坐標(biāo)系重合。

    圖2 坐標(biāo)系約定

    (2)

    2.2 AprilTag定位原理

    AprilTag是一個(gè)二維平面目標(biāo)基準(zhǔn)系統(tǒng),被廣泛應(yīng)用于視覺定位,當(dāng)標(biāo)簽貼在目標(biāo)物上時(shí),便可識(shí)別出相機(jī)和標(biāo)簽之間的3D位姿關(guān)系[11]。為了保證在AGV運(yùn)動(dòng)過程中AprilTag信息不丟失,擴(kuò)大識(shí)別范圍,在天花板上貼4塊ID不一樣的標(biāo)簽。以天花板某一點(diǎn)O為原點(diǎn),創(chuàng)建AprilTag坐標(biāo)系,對(duì)稱分布4塊標(biāo)簽。其中AprilTag坐標(biāo)系可以由世界坐標(biāo)系繞x軸旋轉(zhuǎn)180°后再平移得到,可以計(jì)算出旋轉(zhuǎn)矩陣和平移向量為

    (3)

    式中:Cθ=cosθ;Sθ=sinθ;φ、θ、ψ分別為繞x軸、y軸、z軸的旋轉(zhuǎn)角度;t1、t2、t3為水平偏移。

    齊次變換矩陣為

    (4)

    在AGV運(yùn)動(dòng)過程中,若識(shí)別到多個(gè)標(biāo)簽時(shí),為了提高識(shí)別精度,取最近的標(biāo)簽識(shí)別值為測量值。

    2.3 里程計(jì)定位原理

    4個(gè)麥克納姆輪的電機(jī)上,均安裝編碼器,在AGV運(yùn)動(dòng)時(shí),讀取每個(gè)電機(jī)編碼器的數(shù)值,根據(jù)式(1)計(jì)算出AGV當(dāng)前的速度和角速度,通過航跡推算的方法來獲得當(dāng)前的位姿,若采樣間隔為Δt,計(jì)算公式如下[12]:

    (5)

    從公式(5)可以看出,航跡推算的方法是非線性的,將其轉(zhuǎn)化到世界坐標(biāo)系下進(jìn)行線性化。AGV在世界坐標(biāo)系下的速度為

    (6)

    其中,速度和角速度可以根據(jù)式(1)計(jì)算,而式(6)也可以寫成

    (7)

    2.4 加速度計(jì)定位原理

    加速度計(jì)可以獲得車體x、y、z方向的加速度,在靜止?fàn)顟B(tài)下,x、y方向的加速度為0,z方向的加速度值為-g。為了減小加速度的誤差,需要對(duì)加速度計(jì)進(jìn)行校準(zhǔn)和補(bǔ)償。當(dāng)加速度計(jì)靜止時(shí),誤差為當(dāng)前的測量值,表示為

    (8)

    校準(zhǔn)后的實(shí)際加速度值為

    (9)

    將所測量的加速度值轉(zhuǎn)換到世界坐標(biāo)系下,計(jì)算為

    (10)

    3 組合導(dǎo)航系統(tǒng)設(shè)計(jì)

    AGV的定位包括x、y方向的平移和z軸的旋轉(zhuǎn)。單一的傳感器難以獲得較準(zhǔn)確的位姿信息,單目相機(jī)獲得位姿存在較大誤差,且更新頻率較慢,而通過航跡推算計(jì)算的位姿存在明顯的累計(jì)誤差。因此,把第2節(jié)介紹的多種傳感系統(tǒng)進(jìn)行組合,充分利用各導(dǎo)航的信息進(jìn)行互補(bǔ),最終使整個(gè)系統(tǒng)的精度和可靠性得到提高。

    3.1 整體融合方案

    整個(gè)系統(tǒng)的狀態(tài)變量包括[xvxaxyvyayψωz],如果將所有狀態(tài)變量一起計(jì)算,每次進(jìn)行矩陣運(yùn)算次數(shù)是83,極大地增加了計(jì)算量。為了降低計(jì)算量,將狀態(tài)變量分為[xvxax]、[yvyay]、[ψωz]單獨(dú)估計(jì)。其中[ψωz]的估計(jì)采用陀螺儀和AprilTag標(biāo)簽偏航角的測量值進(jìn)行融合。[xvxax]、[yvyay]的估計(jì)采用單目相機(jī)、里程計(jì)、加速度計(jì),x、y方向的數(shù)據(jù)進(jìn)行卡爾曼濾波融合。

    根據(jù)傳感器的特點(diǎn),陀螺儀在一定時(shí)間內(nèi)可具有較高的精度,但長時(shí)間存在積分誤差,而單目相機(jī)在AGV速度較慢或者停止時(shí)具有較好的精度。因此,每次運(yùn)動(dòng)到任務(wù)點(diǎn)時(shí),AGV速度較慢,此時(shí)用單目相機(jī)去校準(zhǔn)陀螺儀。通過這樣的融合,[ψωz]具有較好的全局精度。

    3.2 經(jīng)典卡爾曼濾波

    以x方向?yàn)槔?,?duì)單目相機(jī)和加速度計(jì)的數(shù)據(jù)進(jìn)行卡爾曼濾波,估計(jì)最優(yōu)的系統(tǒng)狀態(tài),假設(shè)系統(tǒng)的當(dāng)前狀態(tài)為x1,t,包括了位移、速度、加速度3個(gè)狀態(tài)變量,x1,t-1為上一時(shí)刻狀態(tài),系統(tǒng)的狀態(tài)方程和量測方程為

    (11)

    式中:F為上一時(shí)刻到下一時(shí)刻的變換矩陣;B為外部輸入u1,t的輸入控制矩陣;w1,t為系統(tǒng)的過程噪聲;H1是系統(tǒng)的觀測矩陣;v1,t為觀測噪聲;zx1,t為單目相機(jī)獲得的x方向的位移量。其中:

    (12)

    卡爾曼濾波的步驟包括兩步:預(yù)測和更新。預(yù)測過程計(jì)算如下[13]:

    (13)

    更新過程計(jì)算如下:

    (14)

    式中:Kx1,t為最優(yōu)卡爾曼增益;Rx1為量測噪聲;Px1,t為當(dāng)前時(shí)刻協(xié)方差矩陣;I為單位陣。

    3.3 局部濾波器的自適應(yīng)卡爾曼濾波

    根據(jù)卡爾曼濾波計(jì)算過程可知,只有當(dāng)動(dòng)態(tài)系統(tǒng)的參數(shù)和噪聲特性參數(shù)都準(zhǔn)確已知的條件下,經(jīng)典卡爾曼濾波才能得到最優(yōu)估計(jì)。然而,在實(shí)際系統(tǒng)中,一些參數(shù)無法獲取或是實(shí)時(shí)變化的,可能導(dǎo)致濾波器的精度降低,甚至可能出現(xiàn)估計(jì)值和實(shí)際值的誤差越來越大,從而引起濾波發(fā)散。為了消除這種誤差,通過自適應(yīng)濾波的方法,采用噪聲估計(jì)器,實(shí)時(shí)修正和估計(jì)觀測噪聲的特性,從而達(dá)到提高濾波精度的目的。

    在卡爾曼濾波中,局部濾波器的量測預(yù)測誤差的公式計(jì)算為

    (15)

    量測噪聲方向的更新過程為

    (16)

    當(dāng)t趨于無窮大時(shí),1/t趨于0,此時(shí)濾波器的自適應(yīng)能力將逐漸減弱。為了始終保持自適應(yīng)能力,將1/t替換為

    (17)

    當(dāng)t趨于無窮大時(shí),γt趨于1-c,始終保持一定的自適應(yīng)能力,其中c取0.999。

    3.4 多尺度融合及容錯(cuò)設(shè)計(jì)

    在實(shí)際系統(tǒng)中,各個(gè)傳感器的更新頻率并不相同。如果卡爾曼濾波器的更新頻率較大,單目相機(jī)的測量值有著嚴(yán)重的滯后;如果更新頻率較小,則不能保證位置估計(jì)的實(shí)時(shí)性。因此,本文提出一種多尺度融合方法,兼有一定的容錯(cuò)能力。

    里程計(jì)更新頻率較快,以里程計(jì)的更新頻率作為基準(zhǔn),通過串口中斷檢測單目相機(jī)是否更新。若更新,則按照式(13)進(jìn)行計(jì)算;若不發(fā)生更新,則認(rèn)為測量值不準(zhǔn)確,完全信任預(yù)測值。整個(gè)過程可以表示為:

    (18)

    式中:當(dāng)單目相機(jī)更新時(shí)δ(t)=1;沒有更新時(shí)δ(t)=0。

    通過上述方法,不僅實(shí)現(xiàn)了非均勻采樣的多尺度系統(tǒng)的融合,也具有一定的容錯(cuò)能力,當(dāng)標(biāo)簽被遮擋或者AGV行駛范圍超出標(biāo)簽范圍而無法識(shí)別時(shí),卡爾曼濾波仍然能夠起作用。當(dāng)單目相機(jī)工作正常時(shí),狀態(tài)變量的精度會(huì)逐漸恢復(fù)到正常值。

    3.5 聯(lián)邦卡爾曼濾波

    多傳感器數(shù)據(jù)融合包括集中式和分散式濾波。集中式濾波即將所有的傳感器數(shù)據(jù)放到同一個(gè)濾波器中進(jìn)行濾波,分散濾波即利用子濾波器處理各系統(tǒng)的信息再經(jīng)過全局濾波實(shí)現(xiàn)融合。集中式濾波雖然精度較好,但由于狀態(tài)維數(shù)高,導(dǎo)致計(jì)算量較大,容錯(cuò)性低,且不便于故障排除。因此,本文采用聯(lián)邦濾波器,雖然輸出結(jié)果是次優(yōu)的,但容錯(cuò)性得到極大的提升,計(jì)算量顯著減少。

    所設(shè)計(jì)的整體濾波器結(jié)構(gòu)如圖3所示,加速度計(jì)數(shù)據(jù)為參考系統(tǒng),里程計(jì)和視覺定位數(shù)據(jù)為子系統(tǒng),經(jīng)自適應(yīng)局部濾波器1、濾波器2后,進(jìn)入主濾波器,最終得到系統(tǒng)的狀態(tài)變量。

    圖3 濾波器結(jié)構(gòu)圖

    將局部濾波器的信息融合為新的全局狀態(tài)估計(jì)的過程為

    (19)

    得到全局估計(jì)后,再經(jīng)反饋得到子濾波器的狀態(tài)、協(xié)方差矩陣:

    (20)

    式中β1、β2為各子系統(tǒng)分配系數(shù),這里均取為0.5。

    4 仿真和實(shí)驗(yàn)

    4.1 仿真與數(shù)據(jù)分析

    根據(jù)上一節(jié)設(shè)計(jì)的自適應(yīng)聯(lián)邦卡爾曼濾波器,進(jìn)行仿真驗(yàn)證。假設(shè)載體初始狀態(tài)為x=[0 m 1 m/s 0 m/s2],保持1 m/s的勻速運(yùn)動(dòng),0.1 s更新一次。假定Openmv的方均根誤差為1 m,里程計(jì)方均根誤差為0.5 m/s,陀螺儀方均根誤差為0.5 m/s2,經(jīng)濾波后偏航角方均根誤差為0.5°。Openmv未濾波前的誤差、聯(lián)邦卡爾曼濾波、自適應(yīng)聯(lián)邦卡爾曼濾波后誤差的仿真結(jié)果如圖4所示,并計(jì)算得到未濾波前的方均根誤差約為1 m,聯(lián)邦濾波后的方均根誤差約為0.49 m,自適應(yīng)聯(lián)邦濾波后的方均根誤差約為0.35 m。從仿真圖和統(tǒng)計(jì)數(shù)據(jù)可以明顯地看出,自適應(yīng)聯(lián)邦卡爾曼濾波器濾波后精度有明顯的提升。

    圖4 自適應(yīng)聯(lián)邦濾波器仿真圖

    為了驗(yàn)證多尺度融合效果,考慮傳感器的更新頻率并進(jìn)行仿真,仿真條件除更新頻率外與上面仿真一致。假定里程計(jì)0.1 s更新一次,單目相機(jī)1 s更新一次,均采用自適應(yīng)聯(lián)邦卡爾曼濾波,分別使用多尺度融合和不使用多尺度融合。濾波后誤差的仿真結(jié)果如圖5所示,可以發(fā)現(xiàn)沒有進(jìn)行多尺度融合的濾波效果明顯變差,這是因?yàn)楦螺^慢的單目相機(jī)存在滯后,通過多尺度融合減小了這種滯后。

    圖5 多尺度融合仿真圖

    4.2 實(shí)驗(yàn)與數(shù)據(jù)分析

    為了進(jìn)一步驗(yàn)證所設(shè)計(jì)的濾波算法,本節(jié)采用麥克納姆輪全向機(jī)器人平臺(tái)進(jìn)行實(shí)驗(yàn)驗(yàn)證,對(duì)2.6 m×2.35 m的矩形軌跡進(jìn)行跟蹤,小車由A點(diǎn)出發(fā),經(jīng)過B、C、D點(diǎn),最終回到A點(diǎn)。實(shí)驗(yàn)平臺(tái)如圖6所示,軌跡測量設(shè)備為激光跟蹤系統(tǒng),將靶球固定在AGV上,激光跟蹤儀實(shí)時(shí)追蹤靶球,獲取靶球的位置。

    圖6 實(shí)驗(yàn)及測量平臺(tái)

    一共進(jìn)行4組實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果如表1所示,由于里程計(jì)存在累計(jì)誤差,雖然一開始具有較好的精度,但隨著運(yùn)動(dòng)距離的增加,實(shí)際軌跡與期望軌跡越來越遠(yuǎn),無法形成完整的閉合軌跡。采用單目相機(jī)定位未濾波得到的軌跡如圖7所示,可以看出具有較大的跟蹤誤差且軌跡不平滑。采用聯(lián)邦卡爾曼濾波的軌跡如圖8所示,相比于沒有濾波之前的軌跡得到了一些改善。采用自適應(yīng)聯(lián)邦卡爾曼濾波得到的軌跡如圖9所示,通過對(duì)比軌跡圖,可以看出定位精度有明顯的改善。圖10為在擾動(dòng)下的軌跡,AGV初始位置并不在原點(diǎn),但隨著時(shí)間的推移,整個(gè)系統(tǒng)會(huì)趨于穩(wěn)定,并得到一個(gè)較為準(zhǔn)確的位姿估計(jì)。圖11是3種濾波方式的誤差比較,從圖中可以看出自適應(yīng)聯(lián)邦卡爾曼濾波具有較好的濾波效果。根據(jù)表1可知,經(jīng)自適應(yīng)聯(lián)邦卡爾曼濾波后最大誤差從50.162 9 mm減小為15.603 0 mm,方均根誤差從13.671 0 mm減小到5.922 4 mm,平均誤差從10.370 9 mm減小到4.698 7 mm,最大定位誤差減小了69%,方均根誤差減小了56.6%,可以看出所設(shè)計(jì)的濾波器可以提高定位精度。

    表1 幾種濾波方式的誤差統(tǒng)計(jì) 單位:mm

    圖7 未濾波軌跡

    圖8 聯(lián)邦濾波軌跡

    圖9 自適應(yīng)聯(lián)邦濾波軌跡

    圖10 初始位置干擾后軌跡

    圖11 幾種濾波方式的誤差圖

    根據(jù)上述的仿真和實(shí)驗(yàn)的結(jié)果,可以推斷出所設(shè)計(jì)的濾波器具有一定的容錯(cuò)能力。當(dāng)單目相機(jī)被遮擋或超出工作范圍時(shí),此時(shí)單目相機(jī)數(shù)據(jù)具有嚴(yán)重的滯后性,但是里程計(jì)仍然能夠進(jìn)行工作,濾波器此時(shí)只信任里程計(jì)。盡管由于里程計(jì)的特性會(huì)導(dǎo)致誤差逐漸增大,但是只要當(dāng)單目相機(jī)重新正常工作,AGV將逐漸減小誤差并回到正常的狀態(tài)。

    5 結(jié)語

    本文設(shè)計(jì)了一套Mecanum輪全向機(jī)器人系統(tǒng),推導(dǎo)了AGV的正逆運(yùn)動(dòng)學(xué)方程。建立了AGV及各傳感器之間的坐標(biāo)關(guān)系,基于單目相機(jī)的AprilTag定位、里程計(jì)定位、加速度計(jì)定位建立定位模型,并將傳感器的非線性問題轉(zhuǎn)換為線性問題。根據(jù)各傳感器的特點(diǎn),設(shè)計(jì)自適應(yīng)聯(lián)邦卡爾曼濾波器,并將不同更新頻率的傳感器進(jìn)行多尺度融合。通過仿真和實(shí)驗(yàn)表明:所設(shè)計(jì)濾波器能夠大幅提高AGV的精度,并具有良好的魯棒性。

    猜你喜歡
    單目里程計(jì)卡爾曼濾波
    室內(nèi)退化場景下UWB雙基站輔助LiDAR里程計(jì)的定位方法
    一種單目相機(jī)/三軸陀螺儀/里程計(jì)緊組合導(dǎo)航算法
    基于模板特征點(diǎn)提取的立體視覺里程計(jì)實(shí)現(xiàn)方法
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    單目SLAM直線匹配增強(qiáng)平面發(fā)現(xiàn)方法
    基于CAD模型的單目六自由度位姿測量
    大角度斜置激光慣組與里程計(jì)組合導(dǎo)航方法
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    基于擴(kuò)展卡爾曼濾波的PMSM無位置傳感器控制
    基于單目立體視覺的三坐標(biāo)在線識(shí)別技術(shù)研究
    大片电影免费在线观看免费| 嫩草影院入口| 五月天丁香电影| 乱人伦中国视频| 亚洲精华国产精华液的使用体验| 黄色视频不卡| 国产精品偷伦视频观看了| 亚洲精品国产一区二区精华液| 丁香六月欧美| 色婷婷久久久亚洲欧美| 亚洲国产精品一区三区| 国产探花极品一区二区| 悠悠久久av| 欧美 日韩 精品 国产| 日韩制服丝袜自拍偷拍| 男人舔女人的私密视频| 一个人免费看片子| 少妇被粗大的猛进出69影院| 欧美人与善性xxx| 日本黄色日本黄色录像| 国产男女超爽视频在线观看| 国产成人精品无人区| 免费人妻精品一区二区三区视频| 精品午夜福利在线看| 亚洲精品国产av成人精品| 伊人久久大香线蕉亚洲五| 捣出白浆h1v1| 制服诱惑二区| 黄色怎么调成土黄色| 亚洲熟女毛片儿| 午夜免费男女啪啪视频观看| 亚洲精品自拍成人| 夫妻午夜视频| 日韩 亚洲 欧美在线| 国产免费又黄又爽又色| 亚洲欧美一区二区三区黑人| 桃花免费在线播放| 国语对白做爰xxxⅹ性视频网站| 亚洲国产欧美日韩在线播放| 亚洲久久久国产精品| 91国产中文字幕| 搡老乐熟女国产| 久久久久国产精品人妻一区二区| 国产欧美日韩一区二区三区在线| 久久久久久人妻| 国产精品二区激情视频| 人人妻人人澡人人爽人人夜夜| 人成视频在线观看免费观看| 香蕉国产在线看| 日韩成人av中文字幕在线观看| 9热在线视频观看99| 国产色婷婷99| 纵有疾风起免费观看全集完整版| 成人国产av品久久久| 欧美在线黄色| 亚洲av日韩在线播放| 国产精品人妻久久久影院| 老司机深夜福利视频在线观看 | 丰满饥渴人妻一区二区三| 欧美变态另类bdsm刘玥| 色精品久久人妻99蜜桃| 成人影院久久| 啦啦啦啦在线视频资源| 久久精品国产亚洲av涩爱| 成人黄色视频免费在线看| 综合色丁香网| 女的被弄到高潮叫床怎么办| 伦理电影大哥的女人| 欧美日韩一级在线毛片| 91精品国产国语对白视频| 男人舔女人的私密视频| 色网站视频免费| 国产精品秋霞免费鲁丝片| netflix在线观看网站| 亚洲av在线观看美女高潮| 人人妻人人澡人人爽人人夜夜| 99香蕉大伊视频| 亚洲,一卡二卡三卡| 亚洲专区中文字幕在线 | 99久国产av精品国产电影| 丝袜喷水一区| 国产女主播在线喷水免费视频网站| 中国三级夫妇交换| 国产精品久久久久成人av| 国产又爽黄色视频| 午夜福利网站1000一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 日韩制服丝袜自拍偷拍| 色网站视频免费| 最新在线观看一区二区三区 | 亚洲综合色网址| 99久国产av精品国产电影| 日韩电影二区| 国产黄色免费在线视频| 色视频在线一区二区三区| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| 亚洲av成人精品一二三区| 大片免费播放器 马上看| 精品亚洲乱码少妇综合久久| 亚洲av日韩在线播放| 国产精品 欧美亚洲| 亚洲五月色婷婷综合| 亚洲男人天堂网一区| 日本午夜av视频| 中文字幕高清在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品久久二区二区91 | 女人高潮潮喷娇喘18禁视频| 制服丝袜香蕉在线| 黑人猛操日本美女一级片| 麻豆乱淫一区二区| 久久精品国产亚洲av涩爱| bbb黄色大片| 一区二区三区乱码不卡18| 在线观看一区二区三区激情| 麻豆av在线久日| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 免费观看人在逋| 国产精品免费大片| 久久精品久久久久久久性| 伦理电影大哥的女人| 亚洲一级一片aⅴ在线观看| 国产精品熟女久久久久浪| 午夜免费观看性视频| 十八禁人妻一区二区| 久久99热这里只频精品6学生| 国产日韩欧美视频二区| 久久久久久久国产电影| 国产深夜福利视频在线观看| 国产精品香港三级国产av潘金莲 | 国产又色又爽无遮挡免| 大陆偷拍与自拍| 黄片无遮挡物在线观看| 日韩中文字幕视频在线看片| 日日爽夜夜爽网站| 国产精品一二三区在线看| 欧美97在线视频| 午夜福利在线免费观看网站| 亚洲精品,欧美精品| 又粗又硬又长又爽又黄的视频| 在线观看免费高清a一片| 亚洲精品国产区一区二| 老司机影院毛片| 国产成人免费观看mmmm| 国产黄色视频一区二区在线观看| 免费在线观看完整版高清| 看非洲黑人一级黄片| 99九九在线精品视频| 欧美在线黄色| 18在线观看网站| 天美传媒精品一区二区| a 毛片基地| 欧美精品高潮呻吟av久久| 久久国产精品大桥未久av| 国产 一区精品| av国产久精品久网站免费入址| 91成人精品电影| 日韩av不卡免费在线播放| 亚洲精品久久午夜乱码| 午夜福利视频在线观看免费| 亚洲精品日韩在线中文字幕| 国产1区2区3区精品| 国产日韩欧美在线精品| 欧美精品一区二区免费开放| 婷婷色综合www| 蜜桃在线观看..| 精品国产乱码久久久久久小说| 好男人视频免费观看在线| 久久狼人影院| 亚洲色图 男人天堂 中文字幕| 在线观看国产h片| 免费黄网站久久成人精品| 在线看a的网站| 人人澡人人妻人| 国产日韩欧美亚洲二区| 亚洲精品视频女| 啦啦啦在线免费观看视频4| 国产野战对白在线观看| 丝袜人妻中文字幕| 亚洲国产欧美网| 日韩欧美一区视频在线观看| 亚洲天堂av无毛| 无限看片的www在线观看| 1024香蕉在线观看| 少妇人妻久久综合中文| 精品少妇一区二区三区视频日本电影 | 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 97在线人人人人妻| 精品国产露脸久久av麻豆| 国产福利在线免费观看视频| 国产亚洲av片在线观看秒播厂| 免费黄频网站在线观看国产| 国产成人精品在线电影| 一级,二级,三级黄色视频| 日韩大片免费观看网站| 亚洲精品第二区| 久久精品人人爽人人爽视色| www.熟女人妻精品国产| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 69精品国产乱码久久久| www日本在线高清视频| 精品少妇黑人巨大在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲熟女毛片儿| 国产精品av久久久久免费| av天堂久久9| av在线app专区| 汤姆久久久久久久影院中文字幕| 午夜福利视频精品| 少妇被粗大的猛进出69影院| 最近最新中文字幕免费大全7| 熟女av电影| 国产成人免费观看mmmm| 妹子高潮喷水视频| 亚洲av日韩在线播放| 如何舔出高潮| 制服丝袜香蕉在线| 国产免费视频播放在线视频| 国产成人午夜福利电影在线观看| 午夜av观看不卡| 咕卡用的链子| 中文字幕人妻丝袜一区二区 | 国产一级毛片在线| 欧美精品一区二区大全| netflix在线观看网站| 宅男免费午夜| 大片电影免费在线观看免费| 亚洲av在线观看美女高潮| 十八禁人妻一区二区| 亚洲欧美清纯卡通| 九草在线视频观看| 国产人伦9x9x在线观看| 国产xxxxx性猛交| 欧美日韩成人在线一区二区| 国产精品熟女久久久久浪| 岛国毛片在线播放| 久久久久人妻精品一区果冻| 美女午夜性视频免费| 丝袜人妻中文字幕| 咕卡用的链子| 大香蕉久久网| 久久影院123| 欧美精品高潮呻吟av久久| 免费观看a级毛片全部| 午夜激情av网站| 90打野战视频偷拍视频| 成人18禁高潮啪啪吃奶动态图| 亚洲av在线观看美女高潮| 久久久久久久久久久久大奶| 一区二区三区乱码不卡18| 亚洲,欧美,日韩| 午夜激情久久久久久久| 亚洲国产精品国产精品| 国产一区二区三区av在线| 下体分泌物呈黄色| 美女脱内裤让男人舔精品视频| 精品国产乱码久久久久久小说| 久久久久视频综合| 99精国产麻豆久久婷婷| 在线观看免费午夜福利视频| 欧美国产精品一级二级三级| 国产成人欧美| 这个男人来自地球电影免费观看 | 热99久久久久精品小说推荐| 成年女人毛片免费观看观看9 | 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 亚洲av综合色区一区| 男女免费视频国产| 一级毛片我不卡| 日本黄色日本黄色录像| 两个人看的免费小视频| 午夜福利一区二区在线看| 国产成人免费观看mmmm| 日韩制服丝袜自拍偷拍| 99香蕉大伊视频| 国产乱来视频区| 欧美精品一区二区大全| 国语对白做爰xxxⅹ性视频网站| 高清视频免费观看一区二区| 久久ye,这里只有精品| 一级爰片在线观看| av片东京热男人的天堂| 亚洲美女搞黄在线观看| 69精品国产乱码久久久| 超碰97精品在线观看| 亚洲一区二区三区欧美精品| 丝袜脚勾引网站| 美女脱内裤让男人舔精品视频| 精品一区二区三区av网在线观看 | 国产成人免费观看mmmm| 九九爱精品视频在线观看| 最新的欧美精品一区二区| 在现免费观看毛片| 国产av码专区亚洲av| 黄色毛片三级朝国网站| 久久鲁丝午夜福利片| 伦理电影免费视频| 久热这里只有精品99| 成年女人毛片免费观看观看9 | 秋霞伦理黄片| 精品久久蜜臀av无| 一个人免费看片子| 最近最新中文字幕免费大全7| 欧美 日韩 精品 国产| av在线观看视频网站免费| 九色亚洲精品在线播放| 亚洲av欧美aⅴ国产| 日本欧美国产在线视频| 肉色欧美久久久久久久蜜桃| 美国免费a级毛片| 亚洲av日韩精品久久久久久密 | 国语对白做爰xxxⅹ性视频网站| 国产探花极品一区二区| xxxhd国产人妻xxx| 深夜精品福利| 我的亚洲天堂| 亚洲男人天堂网一区| 欧美精品av麻豆av| 国产精品国产av在线观看| 亚洲国产中文字幕在线视频| 久久精品熟女亚洲av麻豆精品| 极品人妻少妇av视频| 久久久久久人妻| 中文字幕人妻熟女乱码| 纯流量卡能插随身wifi吗| 在线观看免费日韩欧美大片| 麻豆乱淫一区二区| 国产精品国产三级专区第一集| 一边摸一边抽搐一进一出视频| 亚洲成人av在线免费| 在线观看免费高清a一片| 日韩,欧美,国产一区二区三区| av一本久久久久| 精品一品国产午夜福利视频| 日韩av免费高清视频| 国产日韩欧美在线精品| 成人国语在线视频| 国产日韩欧美在线精品| av网站在线播放免费| 国产一区有黄有色的免费视频| 精品国产一区二区久久| 伊人久久国产一区二区| 亚洲国产av新网站| 丁香六月天网| 久久久久国产一级毛片高清牌| 国产深夜福利视频在线观看| 99久久精品国产亚洲精品| 中文字幕制服av| 欧美人与性动交α欧美精品济南到| 韩国av在线不卡| 午夜日本视频在线| 欧美成人午夜精品| 汤姆久久久久久久影院中文字幕| 在线免费观看不下载黄p国产| 亚洲av日韩在线播放| 亚洲四区av| 国产成人精品久久久久久| 亚洲综合精品二区| 亚洲欧美一区二区三区黑人| 下体分泌物呈黄色| 亚洲精品视频女| 韩国av在线不卡| 十八禁人妻一区二区| 亚洲欧美激情在线| 国产精品无大码| 国产精品久久久久久精品古装| 丝袜人妻中文字幕| 国产精品成人在线| 久热爱精品视频在线9| 午夜日本视频在线| 国产精品一二三区在线看| 国产精品偷伦视频观看了| 在线观看www视频免费| 制服诱惑二区| 免费黄网站久久成人精品| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 精品亚洲乱码少妇综合久久| av又黄又爽大尺度在线免费看| 精品亚洲成国产av| 一边亲一边摸免费视频| 国产av一区二区精品久久| 国产视频首页在线观看| 亚洲精品美女久久久久99蜜臀 | 91精品伊人久久大香线蕉| 人人妻人人爽人人添夜夜欢视频| 成人漫画全彩无遮挡| 国产日韩欧美亚洲二区| 国产极品粉嫩免费观看在线| 国产亚洲av片在线观看秒播厂| 99久久精品国产亚洲精品| 丝袜喷水一区| 午夜久久久在线观看| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 中文字幕人妻丝袜一区二区 | 国产有黄有色有爽视频| av女优亚洲男人天堂| 一级片免费观看大全| 侵犯人妻中文字幕一二三四区| 国产精品 国内视频| 亚洲欧美中文字幕日韩二区| 麻豆av在线久日| 国产亚洲一区二区精品| 欧美日韩一级在线毛片| 丁香六月天网| 男女边吃奶边做爰视频| 中文字幕人妻丝袜一区二区 | 免费女性裸体啪啪无遮挡网站| 99九九在线精品视频| 男女之事视频高清在线观看 | a 毛片基地| 十八禁人妻一区二区| 欧美亚洲日本最大视频资源| 久久久精品国产亚洲av高清涩受| 日日啪夜夜爽| 国产欧美日韩综合在线一区二区| 久久久久网色| 色94色欧美一区二区| 丰满乱子伦码专区| 国产在线一区二区三区精| 国产 一区精品| 狠狠精品人妻久久久久久综合| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 99久久精品国产亚洲精品| 免费少妇av软件| 欧美少妇被猛烈插入视频| 高清视频免费观看一区二区| 色播在线永久视频| 热99国产精品久久久久久7| 99国产综合亚洲精品| 国产老妇伦熟女老妇高清| 天天躁夜夜躁狠狠久久av| 性色av一级| 美女主播在线视频| 久久久国产欧美日韩av| 男女高潮啪啪啪动态图| 19禁男女啪啪无遮挡网站| 久久精品国产综合久久久| 久久人妻熟女aⅴ| 一二三四在线观看免费中文在| 亚洲精品视频女| 黄色 视频免费看| 99久久精品国产亚洲精品| 久久狼人影院| 午夜免费鲁丝| 欧美日韩综合久久久久久| 老司机亚洲免费影院| 亚洲av中文av极速乱| 超碰97精品在线观看| 精品一区二区三区四区五区乱码 | 久久久精品免费免费高清| 91老司机精品| 色网站视频免费| 亚洲一区中文字幕在线| 中文欧美无线码| 乱人伦中国视频| 2018国产大陆天天弄谢| 日韩av免费高清视频| 成年动漫av网址| 黄片播放在线免费| 天天添夜夜摸| 国产精品99久久99久久久不卡 | 日韩制服丝袜自拍偷拍| 高清av免费在线| 亚洲,欧美精品.| 一二三四中文在线观看免费高清| 丝袜美腿诱惑在线| 国产亚洲av高清不卡| 99热国产这里只有精品6| av在线老鸭窝| 精品第一国产精品| 欧美日韩综合久久久久久| a 毛片基地| 免费观看a级毛片全部| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 精品少妇黑人巨大在线播放| 高清欧美精品videossex| 日韩一本色道免费dvd| 国产精品久久久久久久久免| 中文字幕色久视频| 欧美黄色片欧美黄色片| 99九九在线精品视频| 天天躁夜夜躁狠狠久久av| 久久久久久久大尺度免费视频| 极品人妻少妇av视频| 久久精品aⅴ一区二区三区四区| 激情五月婷婷亚洲| 成人亚洲精品一区在线观看| 水蜜桃什么品种好| 操出白浆在线播放| avwww免费| 亚洲 欧美一区二区三区| 精品免费久久久久久久清纯 | 欧美黄色片欧美黄色片| 精品久久久精品久久久| 国产极品天堂在线| 亚洲国产精品国产精品| 欧美精品av麻豆av| 成人亚洲欧美一区二区av| 最近的中文字幕免费完整| 日韩精品有码人妻一区| 日韩精品免费视频一区二区三区| 黑人欧美特级aaaaaa片| 免费观看人在逋| 天堂俺去俺来也www色官网| 精品一区二区免费观看| 亚洲人成网站在线观看播放| av卡一久久| av网站免费在线观看视频| 国产欧美日韩一区二区三区在线| 亚洲精品乱久久久久久| 一本一本久久a久久精品综合妖精| 91精品伊人久久大香线蕉| 国产精品久久久久久人妻精品电影 | 国产男女超爽视频在线观看| 精品少妇一区二区三区视频日本电影 | 天堂8中文在线网| 国产有黄有色有爽视频| 精品国产超薄肉色丝袜足j| 日韩大片免费观看网站| 国产精品二区激情视频| 老司机亚洲免费影院| 日本午夜av视频| 美女午夜性视频免费| av国产精品久久久久影院| 亚洲一码二码三码区别大吗| 三上悠亚av全集在线观看| 日本猛色少妇xxxxx猛交久久| 精品少妇黑人巨大在线播放| 女性被躁到高潮视频| 99香蕉大伊视频| 久久精品熟女亚洲av麻豆精品| 欧美日韩视频高清一区二区三区二| 夫妻性生交免费视频一级片| 狂野欧美激情性bbbbbb| 亚洲国产av影院在线观看| 午夜福利一区二区在线看| 嫩草影院入口| 十八禁人妻一区二区| 91成人精品电影| 综合色丁香网| 国产成人精品无人区| 999久久久国产精品视频| 国产精品 欧美亚洲| 欧美最新免费一区二区三区| 亚洲熟女毛片儿| 一区二区日韩欧美中文字幕| 久久韩国三级中文字幕| 国产黄色免费在线视频| 午夜福利视频在线观看免费| 我的亚洲天堂| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区蜜桃| 成人国语在线视频| 老汉色av国产亚洲站长工具| 欧美日韩av久久| 精品一区二区三区四区五区乱码 | 久久久久精品人妻al黑| 亚洲国产精品成人久久小说| 国产精品成人在线| 男人添女人高潮全过程视频| 午夜福利视频精品| 国产精品国产三级专区第一集| 国产乱人偷精品视频| 久久热在线av| 天堂中文最新版在线下载| 女的被弄到高潮叫床怎么办| 免费av中文字幕在线| 美女脱内裤让男人舔精品视频| 曰老女人黄片| 视频区图区小说| 宅男免费午夜| 国产精品嫩草影院av在线观看| 国产精品无大码| 少妇精品久久久久久久| 啦啦啦在线免费观看视频4| 亚洲激情五月婷婷啪啪| 久久久国产欧美日韩av| 一区二区三区四区激情视频| 亚洲精品视频女| 少妇精品久久久久久久| 久久 成人 亚洲| 午夜日本视频在线| 亚洲欧美色中文字幕在线| 啦啦啦视频在线资源免费观看| av网站免费在线观看视频| 亚洲少妇的诱惑av| 黑人欧美特级aaaaaa片| 人成视频在线观看免费观看| 高清在线视频一区二区三区| 色视频在线一区二区三区| 亚洲国产欧美一区二区综合| 女人高潮潮喷娇喘18禁视频| xxxhd国产人妻xxx| 亚洲精品视频女| 午夜影院在线不卡| 亚洲成人手机| 国产一区二区三区av在线| 在线观看免费视频网站a站| 人妻人人澡人人爽人人| 99精品久久久久人妻精品| 国产一卡二卡三卡精品 | 美女高潮到喷水免费观看| 麻豆精品久久久久久蜜桃|