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

    融入二維碼信息的自適應(yīng)蒙特卡洛定位算法

    2019-08-01 01:54:12胡章芳曾林全羅元羅鑫趙立明
    計算機(jī)應(yīng)用 2019年4期
    關(guān)鍵詞:二維碼

    胡章芳 曾林全 羅元 羅鑫 趙立明

    摘 要:蒙特卡洛定位(MCL)算法存在計算量大、定位精度差的問題,由于二維碼具有攜帶信息的多樣性、二維碼識別的方便性與易用性的特點(diǎn),提出一種融入二維碼信息的自適應(yīng)蒙特卡洛定位算法。首先,利用二維碼提供的絕對位置信息修正里程計模型的累計誤差后進(jìn)行采樣;然后,采用激光傳感器提供的觀測模型確定粒子的重要性權(quán)重;最后,因?yàn)橹夭蓸硬糠植捎霉潭颖炯瘯?dǎo)致大計算量,所以利用Kullback-Leibler距離(KLD)進(jìn)行重采樣,根據(jù)粒子在狀態(tài)空間的分布情況自適應(yīng)調(diào)整下一次迭代所需粒子數(shù),從而減小計算量?;谝苿訖C(jī)器人進(jìn)行的實(shí)驗(yàn)結(jié)果表明,改進(jìn)算法與傳統(tǒng)蒙特卡洛算法相比定位精度提高了15.09%,時間縮短了15.28%。

    關(guān)鍵詞:蒙特卡洛定位;里程計運(yùn)動模型;觀測模型;二維碼;Kullback-Leibler距離采樣

    中圖分類號:TP242.6

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1001-9081(2019)04-0989-05

    Abstract: Monte Carlo Localization (MCL) algorithm has many problems such as large computation and poor positioning accuracy. Because of the diversity of information carried by two-dimensional code and usability and convenience of two-dimensional code recognition, an adaptive MCL algorithm integrated with two-dimensional code information was proposed. Firstly, the cumulative error of odometer model was corrected by absolute position information provided by two-dimensional code and then sampling was performed. Sencondly, the measurement model provided by laser sensor was used to determine the importance weights of the particles. Finally, as fixed sample set used in the resampling part caused large computation, Kullback-Leibler Distance (KLD) was utilized in resampling to reduce the computation by adaptively adjusting the number of particles required for the next iteration according to the distribution of particles in state space. Experimental result on the mobile robot show that the proposed algorithm improves the localization accuracy by 15.09% and reduces the localization time by 15.28% compared to traditional Monte-Carlo algorithm.

    Key words: Monte-Carlo Localization (MCL); odometer motion model; measurement model; two-dimensional coding; Kullback-Leibler Distance (KLD) sampling

    0?引言

    在環(huán)境地圖已知的前提下,移動機(jī)器人根據(jù)環(huán)境感知和自身運(yùn)動確定其在環(huán)境中的位姿問題稱為定位問題。蒙特卡洛定位(Monte Carlo Localization, MCL)算法[1]以運(yùn)動模型采樣,并結(jié)合觀測模型評估每個粒子的重要性權(quán)重,得到系統(tǒng)狀態(tài)的后驗(yàn)信度分布。

    該算法成功應(yīng)用于移動機(jī)器人領(lǐng)域,適用于局部定位和全局定位兩類問題。里程計運(yùn)動模型通過整合車輪上的光電編碼器信息,進(jìn)而獲得機(jī)器人相對于上一采樣時刻位姿的相對差,在固定的時間間隔可進(jìn)行位姿估計。但是由于漂移或打滑等因素的影響,導(dǎo)致運(yùn)動模型的精度隨時間的增大而下降,MCL算法的先驗(yàn)分布依賴于里程計運(yùn)動模型的精度,若預(yù)測粒子分布與實(shí)際位姿相差太大易導(dǎo)致定位定位誤差增大;另外,重采樣后粒子會產(chǎn)生退化效應(yīng),粒子多樣性降低,固定大樣本粒子會導(dǎo)致計算資源浪費(fèi),因此學(xué)者們一直在研究如何解決這兩類問題。Borenstein等[2]在20世紀(jì)90年代就將里程計誤差分為系統(tǒng)誤差和非系統(tǒng)誤差兩個部分,并提出了一種里程計標(biāo)定方法UMBmark消除系統(tǒng)誤差對里程計精度影響,機(jī)器人按預(yù)定軌跡運(yùn)動幾次便可標(biāo)定出差動輪模型參數(shù)。Yap等[3]利用最大期望(Expectation Maximization, EM)算法并結(jié)合環(huán)境建圖來同時計算里程計運(yùn)動模型和激光器觀測模型的參數(shù),最后實(shí)現(xiàn)了在線自適應(yīng)標(biāo)定。Alhashimi等[4]改進(jìn)蒙特卡洛算法的觀測模型,通過設(shè)置的閾值來確定粒子樣本集合的大小,有效減小了計算量。黃露等[5]設(shè)計人工路標(biāo)并建立了路標(biāo)庫來修正里程計的累計誤差;

    但是由于人工路標(biāo)設(shè)計復(fù)雜且路標(biāo)庫建立需要大量數(shù)據(jù)導(dǎo)致定位精度低、計算量大,而且準(zhǔn)確性也無法滿足要求。

    因此,本文提出一種融入二維碼信息的自適應(yīng)蒙特卡洛定位算法(Adaptive MCL algorithm Integrated with two-dimensional code information, I-AMCL),在采樣過程中融入二維碼攜帶的絕對位置信息,激光傳感器信息用于更新粒子權(quán)重,并采用Kullback-Leibler距離(Kullback-Leibler Distance, KLD)[6]重采樣,根據(jù)采樣在狀態(tài)空間的近似分布的統(tǒng)計界限來確定粒子數(shù)。最后,通過實(shí)驗(yàn)驗(yàn)證在保證最優(yōu)位姿估計和穩(wěn)定性不受影響的情況下,本文算法減小計算量,并提高定位精度。

    假設(shè)當(dāng)前的粒子群St只依賴于上一時刻的粒子群St-1,初始置信度bel(x0)由先驗(yàn)分布 p(x0)隨機(jī)產(chǎn)生的m個這樣的粒子得到,并為每個粒子分配相同的重要性因子m-1。該算法第3)行利用運(yùn)動模型進(jìn)行采樣,以當(dāng)前置信度為起點(diǎn)分配粒子;第4)行觀測模型對每個粒子的權(quán)重進(jìn)行更新;第5)~11)行根據(jù)粒子權(quán)重重采樣得到新的粒子群,并不斷迭代。

    2?本文算法

    將視覺傳感器提供的二維碼的絕對位置信息融入到里程計運(yùn)動模型中,可修正里程計的累計誤差;而每次迭代過程中的重采樣會帶來繁重的計算量,

    因此本文提出一種融入二維碼信息的自適應(yīng)蒙特卡洛定位算法(I-AMCL),以提高定位精度,減少計算量。

    2.1?運(yùn)動模型采樣

    2.1.1?里程計運(yùn)動模型

    里程計使用的機(jī)器人內(nèi)部坐標(biāo)和全局世界坐標(biāo)之間不存在固定的坐標(biāo)變換,因此里程計模型的建立需要機(jī)器人內(nèi)部編碼器觀測的相對運(yùn)動信息。在時間間隔(t-1,t]內(nèi),給定運(yùn)動信息ut為:

    這就是里程計誤差,在機(jī)器人定位過程中容易導(dǎo)致定位誤差。

    2.1.2?二維碼信息

    由于直接將視覺傳感器識別的離散坐標(biāo)融入里程計模型誤差較大,所以應(yīng)提供機(jī)器人位置的直接觀測信息,因此需要計算出二維碼中心點(diǎn)在車體坐標(biāo)系中的精確位置[7-9]。

    為了計算三維空間中物體的幾何信息,需要對攝像機(jī)進(jìn)行標(biāo)定,由此得出攝像機(jī)的內(nèi)參數(shù)矩陣A和外參數(shù)矩陣R,確定攝像機(jī)相對于世界坐標(biāo)系的方位。

    由于攝像頭焦距不為0,有矩陣A和R可逆??赏ㄟ^下式計算拍攝照片中的像素點(diǎn)在車體坐標(biāo)系中的實(shí)際坐標(biāo)。

    當(dāng)知道二維碼中心在圖像中的坐標(biāo)(uc,vc),就可以計算出其在車體坐標(biāo)系中的坐標(biāo)(Xc,Yc,Zc),因此本文使用一個簡單易用的二維碼識別庫(Zbar庫),獲取DM碼三個位置探測圖形在圖像中的坐標(biāo)(u1,v1)、(u2,v2)、(u3,v3)。

    二維碼中心位于對角探測圖形連線中心上。由于已知探測圖形的位置,可計算二維碼的中心(uc,vc)在圖形中的坐標(biāo)為:(uc,vc)=(u1+u32,v1+v32)(14)

    將(uc,vc)代入就可以獲取二維碼在車體坐標(biāo)系中的坐標(biāo)xc=(Xc,Yc,Zc)。

    2.1.3?融入二維碼信息的采樣

    運(yùn)動模型采樣將初始姿態(tài)xt-1、里程計讀數(shù)ut和相機(jī)讀數(shù)xc作為輸入,輸出一個隨機(jī)xt服從p(xtut,xt-1)的分布。機(jī)器人在運(yùn)動過程中沒有掃描到二維碼時,首先由式(2)~(8)推測一個位姿xt,當(dāng)相機(jī)獲取二維碼信息,與此時里程計采樣輸出的位姿xt進(jìn)行比較,兩者誤差大于臨界值τ,則采樣算法輸出的值為二維碼坐標(biāo)的絕對值xc,并令當(dāng)前時刻的位姿xt為二維碼信息提供的位姿,繼續(xù)進(jìn)行采樣直到掃描下一個二維碼。

    2.2?激光傳感器的似然域

    其中:xt是機(jī)器人的位姿;zt是t時刻的觀測;M是環(huán)境地圖。假設(shè)每個觀測束噪聲之間獨(dú)立。

    機(jī)器人在定位時需將傳感器掃描的終點(diǎn)zt映射到地圖的全局坐標(biāo)空間。令xt=(x, y,θ)T表示機(jī)器人在t時刻的位姿。用(xk,sens, yk,sens)T表示安裝在機(jī)器人上的激光傳感器局部坐標(biāo)系統(tǒng)的位置,用θk,sens表示激光光束相對于機(jī)器人航向的角度。觀測zkt的終點(diǎn)經(jīng)過三角變換映射到全局坐標(biāo)系:

    傳感器檢測到障礙物時這些坐標(biāo)才具有意義,因此最大距離讀數(shù)將被丟棄。

    三個參數(shù)zhit、zrand和zmax進(jìn)行加權(quán)平均混合,并且zhit+zrand+zmax=1。

    2.3?基于KLD的重采樣

    用于表示置信度的采樣集的大小對粒子濾波器的效率影響程度很大。為了使移動機(jī)器人能夠處理全局定位和位置跟蹤問題,必須選擇大樣本集合來避免采樣消耗引起的發(fā)散,若采用固定大小的樣本集合,會導(dǎo)致非常大的計算量。KLD采樣能隨時間自適應(yīng)改變粒子數(shù)。

    KLD采樣對每次粒子濾波迭代都以概率1-δ確定樣本數(shù),使得真實(shí)的后驗(yàn)分布與基于采樣的近似分布之間的誤差小于ε,由此確定重采樣樣本集合的大小[10-11]。當(dāng)粒子數(shù)m滿足一定值時,可以保證概率的真實(shí)值與估計值之間的K-L距離小于閾值ε,此時m的值為:

    2.4?融入二維碼信息的自適應(yīng)MCL算法

    與蒙特卡洛算法相比,本文算法I-AMCL除了將t-1時刻的采樣集合St-1、地圖M和最新的控制ut及觀測zkt作為輸入,KLD采樣還把一個加權(quán)后采樣集合和二維碼提供的全局坐標(biāo)信息xc與臨界值τ作為輸入,且要求輸入統(tǒng)計誤差界限ε和δ。

    算法I-AMCL在滿足第27)行統(tǒng)計界限之前,KLD采樣將一直產(chǎn)生粒子。這個界限基于粒子所覆蓋的狀態(tài)空間的體積,粒子覆蓋的體積由覆蓋在三維狀態(tài)空間的位來衡量,直方圖H的位有兩種狀態(tài):空或被一個粒子占據(jù)(第3)~5)行),這個粒子是從St-1樣本集合中提取出來的(第7)行)?;谶@個粒子,一個新粒子被預(yù)測、加權(quán),并插入到新采樣集合,這里預(yù)測采用的是二維碼修正過的里程計模型(第8)~14)行),加權(quán)采用的激光傳感器的似然域(第15)~21)行)。h是至少填充了一個粒子的直方圖的位數(shù),如果新產(chǎn)生的粒子落入直方圖的一個空位中,那么非空位數(shù)目h是增加的,并且該位被標(biāo)記為非空(第22~30)行)。本文算法I-AMCL如算法2所示。

    3?實(shí)驗(yàn)與結(jié)果分析

    3.1?仿真研究

    首先通過仿真實(shí)驗(yàn)進(jìn)行研究。使用Matlab 2012a在計算機(jī)上進(jìn)行仿真,I-AMCL參數(shù)設(shè)置如下:c=0.7,ε=0.1,τ=0.01,δ=0.01,實(shí)驗(yàn)所用初始粒子數(shù)為1000。從圖1可看出,I-AMCL與真實(shí)狀態(tài)值更接近,有相近的收斂速度,由此說明本文算法估計精度更高。

    3.2?系統(tǒng)實(shí)現(xiàn)

    3.2.1?實(shí)驗(yàn)平臺

    為驗(yàn)證本文算法I-AMCL的有效性,采用如圖2所示的室內(nèi)移動機(jī)器人來進(jìn)行。

    在倉庫物流場景中,機(jī)器人根據(jù)訂單任務(wù)在粘貼有二維碼的位置進(jìn)行裝貨卸貨,該機(jī)器人采用光電編碼器進(jìn)行相對位置觀測,并配有CPU為1.86GHz運(yùn)行內(nèi)存為2GB的筆記本電腦,機(jī)器人的運(yùn)動速度為V=0.2m/s,采用灰度攝像頭,在1280×1024的分辨率下圖像傳輸速度達(dá)100fps。為保證旋轉(zhuǎn)過程中二維碼始終位于攝像頭視野之內(nèi),并給攝像頭提供穩(wěn)定的照明,本文將攝像頭安裝在驅(qū)動輪之間,使攝像頭的光心和機(jī)器人的旋轉(zhuǎn)中心重合。運(yùn)行時,安裝LED陣列在車輛底部,可以提供穩(wěn)定照明。

    高速二維碼的識別對攝像頭的拍攝幀率、曝光時間等都有要求。由于攝像頭安裝在車底,距離地面的高度有限制,實(shí)測其視野的邊長為0.15m。本文使用的二維碼為DM碼,用于存儲坐標(biāo)信息,坐標(biāo)信息計算方法為,利用精度為0.1cm的卷尺測量二維碼之間的距離,并通過幾何關(guān)系求出坐標(biāo)(X,Y),航向角θ值由機(jī)器人正方向和讀取二維碼的方向之間的夾角確定。然后利用二維碼的位置探測圖形來確認(rèn)機(jī)器人和二維碼間的相對坐標(biāo)。因此,可利用網(wǎng)頁工具或二維碼批量生成工具等,將實(shí)際的坐標(biāo)信息/坐標(biāo)編號存入到DM碼中,由此,當(dāng)機(jī)器人行駛到二維碼上方時,可以通過掃描二維碼來確認(rèn)自身的當(dāng)前坐標(biāo)。二維碼粘貼示例如圖3所示。

    3.2.2?實(shí)驗(yàn)結(jié)果分析

    利用機(jī)器人操作系統(tǒng)(Robot Operating System, ROS)中的gmapping節(jié)點(diǎn)構(gòu)建室內(nèi)環(huán)境地圖[12],初始位置已知,MCL和I-AMCL的初始粒子數(shù)都設(shè)置為1000,其相關(guān)參數(shù)設(shè)置同步仿真實(shí)驗(yàn)。完成機(jī)器人在先驗(yàn)地圖中的精確定位。

    機(jī)器人從A點(diǎn)到B點(diǎn)的規(guī)劃路徑如圖4(a)所示(其中①~④表示二維碼粘貼位置)。圖4(b)~(e)表示機(jī)器人在運(yùn)動過程中通過二維碼粘貼位置時的粒子變化情況,中心大圓點(diǎn)表示機(jī)器人,周圍的小點(diǎn)表示粒子,箭頭表示機(jī)器人移動方向。從圖4可看出移動機(jī)器人在利用I-AMCL進(jìn)行定位時,在初始位置,粒子數(shù)較多并分布在機(jī)器人周圍,隨著機(jī)器人運(yùn)動,觀測值和控制值在不斷矯正機(jī)器人的位姿,其定位所需粒子數(shù)在不斷減少,且機(jī)器人在二維碼粘貼位置的定位精度較高。

    圖5為兩種算法的粒子數(shù)隨機(jī)器人運(yùn)動距離的變化情況,進(jìn)一步說明MCL采用固定粒子數(shù),而I-AMCL在機(jī)器人運(yùn)動過程中可根據(jù)自由空間狀態(tài)來實(shí)時調(diào)整定位所需粒子數(shù),結(jié)合二維碼位姿信息和激光傳感器的觀測信息,移動機(jī)器人定位所需粒子數(shù)為110左右。

    如圖6所示的室內(nèi)環(huán)境為重慶郵電大學(xué)自動化學(xué)院一樓。為模擬倉庫物流環(huán)境,將二維碼粘貼在①~⑥位置,如圖6(a)所示。實(shí)驗(yàn)過程中6個二維碼均能被快速識別,分別記錄機(jī)器人采用兩種算法的運(yùn)動軌跡的實(shí)測值。定位數(shù)據(jù)分析對比如圖6所示:圖6(b)為MCL算法定位路徑,與規(guī)劃路徑相比誤差明顯;圖6(c)為I-AMCL定位路徑,可以明顯看出定位精度得到了明顯的提高,使得定位路徑明顯趨近于規(guī)劃路徑,整個導(dǎo)航系統(tǒng)的誤差大幅度降低。

    表1對比了兩種定位算法的性能,可以說明本文算法I-AMCL有更精確的定位,可根據(jù)機(jī)器人的狀態(tài)實(shí)時調(diào)整采樣粒子數(shù)達(dá)到提高計算效率的目的,實(shí)現(xiàn)機(jī)器人室內(nèi)導(dǎo)航。

    4?結(jié)語

    本文提出了一種融入二維碼信息的自適應(yīng)蒙特卡洛定位(I-AMCL)算法,利用攝像機(jī)提供的二維碼信息修正里程計模型后再采樣,激光傳感器提供的觀測模型用以確定粒子的重要性權(quán)重;并采用KLD重采樣調(diào)整所需粒子數(shù),減小了計算量同時提高了定位精度;最后通過仿真和實(shí)驗(yàn)驗(yàn)證了本文算法的可行性。但是,本文人為地改變了部分環(huán)境地圖信息,導(dǎo)致機(jī)器人的起始點(diǎn)必須在二維碼正上方,在一定程度上缺乏靈活性。

    參考文獻(xiàn)(References)

    [1] THRUN S, FOX D, BURGARD W, et al. Robust Monte Carlo localization for mobile robots[J]. Artificial Intelligence, 2001, 128(1/2): 99-141.

    [2] BORENSTEIN J, FENG L. Measurement and correction of systematic odometry errors in mobile robots[J]. IEEE Transactions on Robotics and Automation, 1997, 12(6): 869-880.

    [3] YAP T N, SHELTON C R. Simultaneous learning of motion and sensor model parameters for mobile robots[C]// Proceedings of the 2008 IEEE International Conference on Robotics and Automation. Piscataway, NJ: IEEE, 2008: 2091-2097.

    [4] ALHASHIMIl A W, HOSTETTLER R, GUSTAFSSON T. An improvement in the observation model for Monte Carlo localization[C]// Proceedings of the 2014 11th International Conference on Informatics in Control, Automation and Robotics. Piscataway, NJ: IEEE, 2014: 498-505.

    [5] 黃露, 朱明. 單目視覺人工路標(biāo)輔助的移動機(jī)器人導(dǎo)航方法[J]. 計算機(jī)系統(tǒng)應(yīng)用, 2018, 27(1): 106-112. (HUANG L, ZHU M. Mobile robot navigation algorithm assisted by monocular vision artificial road sign[J]. Computer Systems & Applications, 2018, 27(1): 106-112.)

    [6] ZHOU T, PENG D, XU C, et al. An adaptive particle filter based on Kullback-Leibler distance for underwater terrain aided navigation with multi-beam sonar[J]. IET Radar, Sonar & Navigation, 2018, 12(4): 433-441.

    [7] 尚明超, 楊斌, 張翠芳. 二維碼修正EKF-SLAM定位的室內(nèi)無人駕駛小車[J]. 單片機(jī)與嵌入式系統(tǒng)應(yīng)用, 2017, 17(7): 63-66. (SHANG M C, YANG B, ZHANG C F. Indoor driverless car with EKF-SLAM localization modified by QR code[J]. Microcontroller and Embedded System Applications, 2017, 17(7): 63-66.)

    [8] LI Z Q, HUANG J D. Study on the use of QR codes as landmarks for indoor positioning: preliminary results[C]// Proceedings of the 2018 IEEE/ION Position, Location and Navigation Symposium. Piscataway, NJ: IEEE, 2018: 1270-1276.

    [9] TRUC N T, KIM Y T. Navigation method of the transportation robot using fuzzy line tracking and QR code recognition[J]. International Journal of humanoid Robotics, 2016, 14(2): 456-460.

    [10] WANG F, LIN Y. Improving particle filter with a new sampling strategy[C]// Proceedings of the 2009 4th International Conference on Computer Science and Convergence Information Technology. Piscataway, NJ: IEEE, 2009: 408-412.

    [11] WANG F, ZHAO Q, ZHANG Y, et al. Particle filtering with multi proposal distributions[J]. International Journal of Communications Network & System Sciences, 2008, 1(1): 22-28.

    [12] 張建偉, 張立新, 胡穎, 等. 開源機(jī)器人操作系統(tǒng)——ROS[M]. 北京: 科學(xué)出版社, 2012: 1-6. (ZHANG J W, ZHANG L X, HU Y, et al. The Open Source Robot Operating System — ROS[M]. Beijing: Science Press, 2012: 1-6.)

    猜你喜歡
    二維碼
    可以吃的二維碼
    二維碼通行券自助發(fā)行系統(tǒng)在新興市場的應(yīng)用
    二維碼迷宮
    二維碼
    童話世界(2020年32期)2020-12-18 18:15:41
    小康二維碼
    心聲歌刊(2020年5期)2020-11-27 01:56:06
    從條形碼到二維碼
    從條形碼到二維碼
    會賺錢的收款二維碼
    讓嚴(yán)肅的二維碼呆萌起來
    二維碼讓政務(wù)公開更直接
    欧美精品一区二区免费开放| 久久精品国产综合久久久| 黄色视频,在线免费观看| 国产精品免费一区二区三区在线 | 日韩欧美国产一区二区入口| 亚洲欧美激情在线| 黄色成人免费大全| 搡老熟女国产l中国老女人| 香蕉久久夜色| 超碰成人久久| 九色亚洲精品在线播放| 侵犯人妻中文字幕一二三四区| 丝袜在线中文字幕| 天堂俺去俺来也www色官网| 亚洲专区国产一区二区| 国产精品久久久久成人av| 99国产精品一区二区三区| 亚洲专区国产一区二区| 亚洲精品成人av观看孕妇| 激情视频va一区二区三区| 国产在线精品亚洲第一网站| 69精品国产乱码久久久| 黄色女人牲交| 成年版毛片免费区| www.熟女人妻精品国产| 麻豆国产av国片精品| 国产欧美日韩一区二区三| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品久久久久久毛片 | 国产亚洲精品久久久久5区| 亚洲九九香蕉| 亚洲av第一区精品v没综合| 亚洲一码二码三码区别大吗| 亚洲黑人精品在线| 1024视频免费在线观看| 女人被狂操c到高潮| 中文字幕人妻熟女乱码| 99国产精品一区二区蜜桃av | 国产熟女午夜一区二区三区| 一a级毛片在线观看| 国产高清视频在线播放一区| 亚洲中文日韩欧美视频| www.精华液| 操出白浆在线播放| 日韩人妻精品一区2区三区| 亚洲男人天堂网一区| 啦啦啦免费观看视频1| 少妇 在线观看| 王馨瑶露胸无遮挡在线观看| 国产精品电影一区二区三区 | a级毛片黄视频| 啦啦啦视频在线资源免费观看| 女人精品久久久久毛片| 国产麻豆69| 国产麻豆69| av网站在线播放免费| 精品国内亚洲2022精品成人 | 女人被躁到高潮嗷嗷叫费观| 999久久久国产精品视频| 久久久久久久国产电影| 久久久国产欧美日韩av| 乱人伦中国视频| tube8黄色片| avwww免费| av视频免费观看在线观看| 亚洲美女黄片视频| 在线观看www视频免费| 黄色成人免费大全| 丝袜在线中文字幕| tube8黄色片| 深夜精品福利| 久久中文看片网| 久久久精品免费免费高清| 欧美日本中文国产一区发布| 黄色视频,在线免费观看| 亚洲精品自拍成人| 国产精品秋霞免费鲁丝片| 人妻久久中文字幕网| 夫妻午夜视频| 一区二区三区国产精品乱码| 亚洲va日本ⅴa欧美va伊人久久| av免费在线观看网站| 免费在线观看影片大全网站| 男女午夜视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品免费一区二区三区在线 | 一边摸一边抽搐一进一小说 | 岛国在线观看网站| 日日摸夜夜添夜夜添小说| av电影中文网址| 在线观看午夜福利视频| 久久久久国内视频| 99香蕉大伊视频| 国产成人av教育| 久久国产精品大桥未久av| 精品一区二区三卡| 村上凉子中文字幕在线| 在线十欧美十亚洲十日本专区| 午夜久久久在线观看| 国产伦人伦偷精品视频| 国产精品二区激情视频| 香蕉丝袜av| 51午夜福利影视在线观看| 亚洲精华国产精华精| 最近最新免费中文字幕在线| 国产熟女午夜一区二区三区| 亚洲一区高清亚洲精品| 中文字幕最新亚洲高清| 韩国av一区二区三区四区| 激情视频va一区二区三区| 亚洲伊人色综图| 亚洲av美国av| 久久人妻福利社区极品人妻图片| 国产精品成人在线| 操美女的视频在线观看| 久久国产精品人妻蜜桃| 精品少妇久久久久久888优播| 欧美精品av麻豆av| 精品欧美一区二区三区在线| 无人区码免费观看不卡| 一区在线观看完整版| 最近最新免费中文字幕在线| 亚洲人成77777在线视频| 日本黄色日本黄色录像| 精品乱码久久久久久99久播| 国产欧美日韩一区二区三| 欧美午夜高清在线| 欧美老熟妇乱子伦牲交| 他把我摸到了高潮在线观看| 91精品三级在线观看| 久久久久精品人妻al黑| 色精品久久人妻99蜜桃| 狠狠狠狠99中文字幕| 最近最新中文字幕大全免费视频| 亚洲欧美一区二区三区黑人| 午夜福利在线观看吧| а√天堂www在线а√下载 | 色婷婷av一区二区三区视频| 黄色片一级片一级黄色片| 久9热在线精品视频| 精品熟女少妇八av免费久了| 国产欧美日韩综合在线一区二区| 亚洲熟女精品中文字幕| 制服诱惑二区| 久99久视频精品免费| 99热国产这里只有精品6| 国产精品美女特级片免费视频播放器 | 超碰97精品在线观看| 18禁观看日本| 热99久久久久精品小说推荐| 香蕉丝袜av| 女警被强在线播放| 男女高潮啪啪啪动态图| 午夜福利影视在线免费观看| 美女视频免费永久观看网站| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品一区在线观看| 国产成人精品在线电影| 国产黄色免费在线视频| 国产成人一区二区三区免费视频网站| 18禁裸乳无遮挡动漫免费视频| 色尼玛亚洲综合影院| av一本久久久久| 国产人伦9x9x在线观看| 成人精品一区二区免费| 一边摸一边做爽爽视频免费| 亚洲欧美色中文字幕在线| 久久香蕉国产精品| 1024视频免费在线观看| 国产区一区二久久| 中文字幕色久视频| 国产aⅴ精品一区二区三区波| 国产亚洲精品一区二区www | 亚洲av日韩在线播放| 国产熟女午夜一区二区三区| 亚洲avbb在线观看| 亚洲自偷自拍图片 自拍| 亚洲精品在线美女| 欧美黑人精品巨大| 免费高清在线观看日韩| 亚洲熟女精品中文字幕| 久久久久视频综合| 看黄色毛片网站| 变态另类成人亚洲欧美熟女 | 99在线人妻在线中文字幕 | 精品熟女少妇八av免费久了| 人人妻人人澡人人看| 一级片免费观看大全| 国产精品自产拍在线观看55亚洲 | 国产精品久久电影中文字幕 | av不卡在线播放| 丰满饥渴人妻一区二区三| 日韩成人在线观看一区二区三区| 黄色视频不卡| 男人操女人黄网站| 午夜成年电影在线免费观看| 热re99久久国产66热| 999久久久精品免费观看国产| 色在线成人网| 亚洲精品美女久久av网站| 亚洲av熟女| 欧美黄色淫秽网站| 精品国内亚洲2022精品成人 | 亚洲三区欧美一区| 中文亚洲av片在线观看爽 | 亚洲精品美女久久av网站| av在线播放免费不卡| 亚洲国产精品sss在线观看 | 亚洲午夜理论影院| 欧美精品亚洲一区二区| 亚洲熟女毛片儿| 欧美亚洲日本最大视频资源| 国产成人系列免费观看| 亚洲视频免费观看视频| 捣出白浆h1v1| 久久精品国产亚洲av高清一级| 亚洲人成电影观看| 嫁个100分男人电影在线观看| 欧美成狂野欧美在线观看| 午夜老司机福利片| 搡老乐熟女国产| 韩国av一区二区三区四区| 91字幕亚洲| 亚洲精品美女久久久久99蜜臀| 亚洲av日韩在线播放| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 午夜视频精品福利| 午夜福利,免费看| 久久精品国产99精品国产亚洲性色 | 99re6热这里在线精品视频| 成在线人永久免费视频| 久久香蕉激情| 国产精品电影一区二区三区 | 母亲3免费完整高清在线观看| 国产成人av教育| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久成人av| 久久精品国产a三级三级三级| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区免费开放| 男女高潮啪啪啪动态图| 亚洲午夜精品一区,二区,三区| 午夜免费成人在线视频| 欧美激情极品国产一区二区三区| 在线观看免费日韩欧美大片| 夫妻午夜视频| 国产精品免费视频内射| 亚洲欧美色中文字幕在线| 69av精品久久久久久| 久久久国产一区二区| 老司机在亚洲福利影院| av天堂在线播放| 亚洲三区欧美一区| 国产不卡一卡二| 搡老岳熟女国产| 亚洲国产精品sss在线观看 | 精品国产乱子伦一区二区三区| 国产精品1区2区在线观看. | 天堂动漫精品| 亚洲第一青青草原| 男女之事视频高清在线观看| 欧美日韩成人在线一区二区| 国产区一区二久久| 午夜福利影视在线免费观看| 亚洲精品在线观看二区| 一个人免费在线观看的高清视频| 国产熟女午夜一区二区三区| 亚洲在线自拍视频| 欧美成人午夜精品| 欧美黑人精品巨大| 欧美乱妇无乱码| 精品一区二区三区av网在线观看| 黄网站色视频无遮挡免费观看| 国产淫语在线视频| 国产精品久久视频播放| 18禁裸乳无遮挡动漫免费视频| 黄色怎么调成土黄色| 动漫黄色视频在线观看| 男人舔女人的私密视频| 搡老岳熟女国产| 欧美激情 高清一区二区三区| 欧美激情极品国产一区二区三区| 久久精品人人爽人人爽视色| 欧美精品人与动牲交sv欧美| 99香蕉大伊视频| 黑人猛操日本美女一级片| 香蕉久久夜色| 丝袜在线中文字幕| 国产蜜桃级精品一区二区三区 | 9色porny在线观看| 91在线观看av| 水蜜桃什么品种好| 久热爱精品视频在线9| 涩涩av久久男人的天堂| 在线天堂中文资源库| 最新的欧美精品一区二区| 在线av久久热| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区精品视频观看| 女人被躁到高潮嗷嗷叫费观| 日韩欧美一区二区三区在线观看 | 午夜精品在线福利| 黑人猛操日本美女一级片| 好看av亚洲va欧美ⅴa在| 丝袜在线中文字幕| 国产蜜桃级精品一区二区三区 | 亚洲国产欧美网| 老司机靠b影院| 国产免费av片在线观看野外av| 大型av网站在线播放| 电影成人av| 国产区一区二久久| 亚洲一区二区三区不卡视频| 18禁观看日本| 91精品三级在线观看| 国产成人一区二区三区免费视频网站| 咕卡用的链子| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美免费精品| 婷婷丁香在线五月| 老司机在亚洲福利影院| 国产精品二区激情视频| 免费在线观看黄色视频的| av电影中文网址| 午夜福利欧美成人| 在线看a的网站| 女人精品久久久久毛片| 免费观看精品视频网站| 久热爱精品视频在线9| 亚洲精品一二三| 青草久久国产| 久久午夜综合久久蜜桃| 欧美日韩亚洲高清精品| 人人妻人人澡人人爽人人夜夜| 18禁黄网站禁片午夜丰满| 精品第一国产精品| 天堂俺去俺来也www色官网| 国产一卡二卡三卡精品| 一级,二级,三级黄色视频| 国产高清国产精品国产三级| 亚洲成人免费电影在线观看| 91精品三级在线观看| av福利片在线| 韩国av一区二区三区四区| 精品高清国产在线一区| 好看av亚洲va欧美ⅴa在| 午夜免费成人在线视频| 狠狠婷婷综合久久久久久88av| 成年女人毛片免费观看观看9 | 涩涩av久久男人的天堂| 精品免费久久久久久久清纯 | 欧美精品亚洲一区二区| 午夜免费成人在线视频| 夜夜爽天天搞| 亚洲免费av在线视频| 一级毛片女人18水好多| a级毛片黄视频| 天天躁夜夜躁狠狠躁躁| 80岁老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人| 69精品国产乱码久久久| 日日爽夜夜爽网站| 日韩免费av在线播放| 五月开心婷婷网| 日韩有码中文字幕| 咕卡用的链子| 欧美午夜高清在线| 18禁黄网站禁片午夜丰满| 高清欧美精品videossex| 国产精品香港三级国产av潘金莲| 国产色视频综合| av不卡在线播放| 国产精品二区激情视频| 精品一区二区三区四区五区乱码| 不卡av一区二区三区| a级片在线免费高清观看视频| 久久九九热精品免费| 看黄色毛片网站| 久久狼人影院| 午夜精品久久久久久毛片777| 国产精品乱码一区二三区的特点 | 国产精品成人在线| 免费久久久久久久精品成人欧美视频| 久久国产精品影院| 丝瓜视频免费看黄片| 在线永久观看黄色视频| 人成视频在线观看免费观看| 亚洲成a人片在线一区二区| 国产精品电影一区二区三区 | 国产精品一区二区在线观看99| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 国产成人欧美| 男女床上黄色一级片免费看| 国产精华一区二区三区| 久久国产精品人妻蜜桃| 精品少妇一区二区三区视频日本电影| 黄色女人牲交| 国产精品 国内视频| 高清av免费在线| 黑人操中国人逼视频| 纯流量卡能插随身wifi吗| avwww免费| 免费高清在线观看日韩| 久久久久久亚洲精品国产蜜桃av| 亚洲成人手机| 亚洲精品在线观看二区| 国产男女超爽视频在线观看| 久久这里只有精品19| 在线国产一区二区在线| 亚洲伊人色综图| 亚洲av电影在线进入| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 99国产综合亚洲精品| aaaaa片日本免费| 男男h啪啪无遮挡| 国产深夜福利视频在线观看| 日韩欧美免费精品| 国产成人啪精品午夜网站| 激情在线观看视频在线高清 | 很黄的视频免费| av国产精品久久久久影院| 天堂俺去俺来也www色官网| 国产免费av片在线观看野外av| 国产黄色免费在线视频| 窝窝影院91人妻| 涩涩av久久男人的天堂| 啪啪无遮挡十八禁网站| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 51午夜福利影视在线观看| svipshipincom国产片| 亚洲精品国产区一区二| 久久中文字幕人妻熟女| 99久久综合精品五月天人人| 亚洲色图av天堂| 在线永久观看黄色视频| 亚洲色图 男人天堂 中文字幕| 国产成人免费无遮挡视频| videos熟女内射| 成人黄色视频免费在线看| av不卡在线播放| 99re6热这里在线精品视频| 久久久久视频综合| 国产亚洲精品久久久久久毛片 | 美女扒开内裤让男人捅视频| 12—13女人毛片做爰片一| 黄色视频,在线免费观看| 高清在线国产一区| 欧美日韩黄片免| 欧美日韩亚洲综合一区二区三区_| 黑人欧美特级aaaaaa片| av视频免费观看在线观看| 久热这里只有精品99| 亚洲黑人精品在线| 啦啦啦在线免费观看视频4| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 欧美日韩成人在线一区二区| 波多野结衣av一区二区av| 99国产精品一区二区蜜桃av | 精品人妻1区二区| 日韩大码丰满熟妇| av电影中文网址| 在线免费观看的www视频| 亚洲欧美激情综合另类| 女同久久另类99精品国产91| 亚洲九九香蕉| 国产高清国产精品国产三级| av中文乱码字幕在线| 在线观看午夜福利视频| 久久精品国产综合久久久| 十八禁网站免费在线| 精品久久蜜臀av无| 亚洲av成人不卡在线观看播放网| 亚洲av日韩在线播放| av一本久久久久| 男人舔女人的私密视频| 久久久久久久久久久久大奶| 午夜成年电影在线免费观看| 天天躁夜夜躁狠狠躁躁| 狠狠婷婷综合久久久久久88av| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 黑丝袜美女国产一区| 人妻久久中文字幕网| 免费在线观看影片大全网站| 亚洲人成77777在线视频| 女人精品久久久久毛片| 国产精品免费视频内射| 啦啦啦 在线观看视频| 777米奇影视久久| 国产在视频线精品| 亚洲第一青青草原| 日本vs欧美在线观看视频| 在线观看66精品国产| 大陆偷拍与自拍| 伦理电影免费视频| 午夜影院日韩av| 午夜精品在线福利| 欧美乱码精品一区二区三区| 国产91精品成人一区二区三区| 黄色片一级片一级黄色片| 九色亚洲精品在线播放| 亚洲中文av在线| 午夜久久久在线观看| 大型黄色视频在线免费观看| 交换朋友夫妻互换小说| 黄色毛片三级朝国网站| 色在线成人网| 久久香蕉激情| 变态另类成人亚洲欧美熟女 | 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 久久久久精品人妻al黑| 成人三级做爰电影| 999精品在线视频| 丝袜美腿诱惑在线| 伦理电影免费视频| 精品亚洲成国产av| 欧美乱妇无乱码| 热99国产精品久久久久久7| 成年人免费黄色播放视频| 美女扒开内裤让男人捅视频| 极品少妇高潮喷水抽搐| av福利片在线| 亚洲欧美激情在线| 女人久久www免费人成看片| 黑人猛操日本美女一级片| 99久久人妻综合| 女性生殖器流出的白浆| 成年女人毛片免费观看观看9 | 妹子高潮喷水视频| 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 丝袜美腿诱惑在线| 真人做人爱边吃奶动态| 人妻一区二区av| 欧美精品av麻豆av| 国产不卡一卡二| 久久草成人影院| 纯流量卡能插随身wifi吗| 精品国产乱码久久久久久男人| 欧美国产精品一级二级三级| 国产精品综合久久久久久久免费 | 少妇的丰满在线观看| www.自偷自拍.com| 如日韩欧美国产精品一区二区三区| 黄色成人免费大全| 欧美激情高清一区二区三区| 欧美最黄视频在线播放免费 | 国产精品久久久久久精品古装| 黄片大片在线免费观看| 精品少妇一区二区三区视频日本电影| 99国产精品一区二区蜜桃av | 中文字幕人妻丝袜制服| 亚洲专区字幕在线| videosex国产| 免费一级毛片在线播放高清视频 | 日日爽夜夜爽网站| 99国产精品免费福利视频| 精品第一国产精品| 亚洲av美国av| 9色porny在线观看| 国产一卡二卡三卡精品| 黄色丝袜av网址大全| 国产又色又爽无遮挡免费看| 国产欧美日韩综合在线一区二区| 久热爱精品视频在线9| 久久精品国产亚洲av高清一级| 91精品国产国语对白视频| 91麻豆精品激情在线观看国产 | 国产极品粉嫩免费观看在线| 亚洲精品国产区一区二| 国产精品久久视频播放| 国产精品乱码一区二三区的特点 | 我的亚洲天堂| 欧美黄色片欧美黄色片| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三| 黄片播放在线免费| 久久久水蜜桃国产精品网| 丁香六月欧美| 一个人免费在线观看的高清视频| 欧美乱码精品一区二区三区| 丝袜美腿诱惑在线| 午夜福利影视在线免费观看| 日韩欧美免费精品| 亚洲自偷自拍图片 自拍| 久久久久久久精品吃奶| 欧美不卡视频在线免费观看 | 国产成人系列免费观看| 成人国语在线视频| 日韩成人在线观看一区二区三区| 成人国语在线视频| 大型黄色视频在线免费观看| 中文欧美无线码| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 中文亚洲av片在线观看爽 | 老熟妇乱子伦视频在线观看| 69av精品久久久久久| av电影中文网址| 五月开心婷婷网| 9191精品国产免费久久| 桃红色精品国产亚洲av| 国产精品国产高清国产av | 国产精品一区二区免费欧美| a级片在线免费高清观看视频| 伦理电影免费视频| 久热这里只有精品99|