• <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ù)公開更直接
    男女国产视频网站| 国产精品一区二区精品视频观看| 久久精品国产亚洲av香蕉五月 | 国产一区二区三区综合在线观看| 亚洲av男天堂| 婷婷成人精品国产| 亚洲激情五月婷婷啪啪| 欧美国产精品va在线观看不卡| 91av网站免费观看| 一区二区三区四区激情视频| 国产伦理片在线播放av一区| 午夜影院在线不卡| 国产主播在线观看一区二区| 日本五十路高清| 色婷婷久久久亚洲欧美| 免费久久久久久久精品成人欧美视频| 电影成人av| 亚洲欧美精品自产自拍| 九色亚洲精品在线播放| 如日韩欧美国产精品一区二区三区| 免费在线观看影片大全网站| a在线观看视频网站| 一边摸一边做爽爽视频免费| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 日韩有码中文字幕| 99久久人妻综合| 久久久久久人人人人人| 在线十欧美十亚洲十日本专区| 美女扒开内裤让男人捅视频| 汤姆久久久久久久影院中文字幕| 在线天堂中文资源库| 日韩精品免费视频一区二区三区| 国产免费av片在线观看野外av| 正在播放国产对白刺激| 亚洲精品美女久久av网站| 久久精品国产综合久久久| 桃红色精品国产亚洲av| 亚洲av日韩精品久久久久久密| 丝瓜视频免费看黄片| 一本色道久久久久久精品综合| 日韩 欧美 亚洲 中文字幕| 男女国产视频网站| 国产一区二区三区综合在线观看| 9191精品国产免费久久| 国产伦人伦偷精品视频| 老司机影院成人| 国产色视频综合| 日韩中文字幕欧美一区二区| 欧美日韩视频精品一区| 女人高潮潮喷娇喘18禁视频| 亚洲成国产人片在线观看| 国产精品香港三级国产av潘金莲| 伦理电影免费视频| 岛国毛片在线播放| 亚洲欧美一区二区三区久久| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲熟妇少妇任你| 大型av网站在线播放| 一区二区三区激情视频| 青草久久国产| 久久久久久人人人人人| 最新在线观看一区二区三区| 亚洲国产看品久久| 熟女少妇亚洲综合色aaa.| 巨乳人妻的诱惑在线观看| 中文精品一卡2卡3卡4更新| 一个人免费看片子| 国产高清国产精品国产三级| 亚洲中文av在线| av不卡在线播放| 伊人亚洲综合成人网| 97精品久久久久久久久久精品| 国产色视频综合| av在线老鸭窝| 一级片'在线观看视频| 国产一卡二卡三卡精品| 99久久综合免费| 狠狠狠狠99中文字幕| 激情视频va一区二区三区| 亚洲国产av新网站| 巨乳人妻的诱惑在线观看| 精品国内亚洲2022精品成人 | 成年人免费黄色播放视频| 在线十欧美十亚洲十日本专区| 亚洲熟女精品中文字幕| 免费少妇av软件| 50天的宝宝边吃奶边哭怎么回事| 国产伦人伦偷精品视频| 婷婷色av中文字幕| 各种免费的搞黄视频| 午夜福利,免费看| 色婷婷久久久亚洲欧美| av福利片在线| 人妻 亚洲 视频| 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 爱豆传媒免费全集在线观看| 一本一本久久a久久精品综合妖精| 我要看黄色一级片免费的| 香蕉国产在线看| 男女无遮挡免费网站观看| 亚洲av片天天在线观看| 最新的欧美精品一区二区| 精品亚洲成国产av| av不卡在线播放| 青草久久国产| 久热这里只有精品99| 国产成+人综合+亚洲专区| 香蕉国产在线看| 国产精品秋霞免费鲁丝片| 亚洲精品国产色婷婷电影| 18禁裸乳无遮挡动漫免费视频| 国产成人av教育| 1024视频免费在线观看| 精品国产乱码久久久久久男人| 伊人久久大香线蕉亚洲五| 久久国产精品人妻蜜桃| 精品福利永久在线观看| 国产亚洲av片在线观看秒播厂| 在线观看免费高清a一片| 黄片大片在线免费观看| 日韩视频在线欧美| 色播在线永久视频| 中文欧美无线码| 欧美在线黄色| 国产亚洲av片在线观看秒播厂| 91麻豆精品激情在线观看国产 | 一区二区三区精品91| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av片天天在线观看| 美女视频免费永久观看网站| 国产成人系列免费观看| 99精国产麻豆久久婷婷| 精品国产乱码久久久久久小说| 国产高清videossex| 国产成+人综合+亚洲专区| 亚洲成人免费av在线播放| 99国产精品免费福利视频| 亚洲av欧美aⅴ国产| 亚洲欧美成人综合另类久久久| 丰满人妻熟妇乱又伦精品不卡| 97人妻天天添夜夜摸| 精品人妻1区二区| 人人妻人人爽人人添夜夜欢视频| 国产精品久久久久久人妻精品电影 | 曰老女人黄片| 一区福利在线观看| 国产福利在线免费观看视频| 精品久久久久久电影网| 日韩中文字幕视频在线看片| 亚洲专区字幕在线| 每晚都被弄得嗷嗷叫到高潮| 欧美激情 高清一区二区三区| 亚洲国产欧美一区二区综合| 狂野欧美激情性xxxx| 国产一区二区 视频在线| 亚洲国产精品成人久久小说| 男女午夜视频在线观看| 亚洲av电影在线进入| 国产国语露脸激情在线看| 亚洲avbb在线观看| 少妇猛男粗大的猛烈进出视频| 婷婷色av中文字幕| 国产精品.久久久| 热99re8久久精品国产| 亚洲 欧美一区二区三区| 成年人午夜在线观看视频| 黑人巨大精品欧美一区二区mp4| 欧美精品高潮呻吟av久久| 久久人妻福利社区极品人妻图片| av片东京热男人的天堂| netflix在线观看网站| 乱人伦中国视频| 天天添夜夜摸| 精品久久久久久电影网| 精品国产国语对白av| 99国产极品粉嫩在线观看| 欧美黑人欧美精品刺激| 亚洲精品久久久久久婷婷小说| 欧美另类一区| 亚洲专区国产一区二区| 久久毛片免费看一区二区三区| 母亲3免费完整高清在线观看| 91麻豆av在线| 日日夜夜操网爽| 免费在线观看影片大全网站| 亚洲成人国产一区在线观看| 一区二区三区精品91| 欧美中文综合在线视频| 少妇 在线观看| 热99久久久久精品小说推荐| 亚洲国产欧美一区二区综合| 欧美日韩亚洲国产一区二区在线观看 | 亚洲avbb在线观看| 天堂8中文在线网| 国产精品一区二区精品视频观看| 日韩欧美国产一区二区入口| 高清视频免费观看一区二区| 久久免费观看电影| 在线av久久热| 国产成人a∨麻豆精品| 97人妻天天添夜夜摸| 美女午夜性视频免费| 精品久久久久久久毛片微露脸 | 精品第一国产精品| 建设人人有责人人尽责人人享有的| 欧美变态另类bdsm刘玥| 久久99一区二区三区| 欧美中文综合在线视频| 极品少妇高潮喷水抽搐| 中文欧美无线码| 一级毛片女人18水好多| 在线观看人妻少妇| 亚洲天堂av无毛| 又黄又粗又硬又大视频| 国产精品九九99| 老司机午夜福利在线观看视频 | 午夜精品国产一区二区电影| 色精品久久人妻99蜜桃| 国产亚洲精品一区二区www | 丝袜喷水一区| 欧美在线一区亚洲| 国产成人av激情在线播放| 国产日韩欧美亚洲二区| 日韩有码中文字幕| 一区二区三区激情视频| 国产精品久久久久久人妻精品电影 | 午夜福利视频精品| 蜜桃在线观看..| 国产精品秋霞免费鲁丝片| av免费在线观看网站| 国产免费现黄频在线看| 午夜免费鲁丝| 80岁老熟妇乱子伦牲交| 老司机午夜福利在线观看视频 | 精品人妻1区二区| 18在线观看网站| 美女福利国产在线| 久久亚洲国产成人精品v| 人人妻人人澡人人看| 777久久人妻少妇嫩草av网站| 99国产精品免费福利视频| 婷婷丁香在线五月| 精品福利观看| 久久中文字幕一级| 亚洲精品av麻豆狂野| 亚洲国产看品久久| 成人手机av| av不卡在线播放| 欧美老熟妇乱子伦牲交| 久久久久久久久久久久大奶| 王馨瑶露胸无遮挡在线观看| 国产欧美日韩一区二区三 | 亚洲精品在线美女| 久久天躁狠狠躁夜夜2o2o| 午夜日韩欧美国产| 18禁国产床啪视频网站| 一级片'在线观看视频| 久久久久久久国产电影| 麻豆av在线久日| 王馨瑶露胸无遮挡在线观看| a级片在线免费高清观看视频| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 国产色视频综合| 男人爽女人下面视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品粉嫩美女一区| 久久ye,这里只有精品| 亚洲少妇的诱惑av| 99精品久久久久人妻精品| 国产精品免费视频内射| 中文字幕另类日韩欧美亚洲嫩草| netflix在线观看网站| 亚洲男人天堂网一区| 亚洲成国产人片在线观看| 飞空精品影院首页| 亚洲精品久久久久久婷婷小说| 欧美 亚洲 国产 日韩一| 99精品欧美一区二区三区四区| 日韩中文字幕视频在线看片| 人人妻,人人澡人人爽秒播| 人成视频在线观看免费观看| 欧美 亚洲 国产 日韩一| 国产亚洲欧美在线一区二区| 正在播放国产对白刺激| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 老熟妇乱子伦视频在线观看 | www.精华液| 国产野战对白在线观看| 免费不卡黄色视频| 老汉色∧v一级毛片| 亚洲av欧美aⅴ国产| 一本综合久久免费| 国产免费现黄频在线看| 中文字幕av电影在线播放| 亚洲精品一二三| 精品久久蜜臀av无| 精品人妻一区二区三区麻豆| 国产麻豆69| 婷婷成人精品国产| 亚洲美女黄色视频免费看| 男人爽女人下面视频在线观看| 黄色视频,在线免费观看| 精品国内亚洲2022精品成人 | 亚洲午夜精品一区,二区,三区| 人人妻人人澡人人爽人人夜夜| 黑人欧美特级aaaaaa片| 搡老岳熟女国产| 精品国产一区二区久久| 黑人巨大精品欧美一区二区蜜桃| 国产av又大| 亚洲人成电影观看| 精品一区二区三区av网在线观看 | 久久久久网色| 久久亚洲国产成人精品v| 国产国语露脸激情在线看| 99国产精品一区二区三区| 亚洲第一青青草原| 中文字幕人妻熟女乱码| 老司机深夜福利视频在线观看 | 精品一区二区三区四区五区乱码| 搡老熟女国产l中国老女人| 一个人免费看片子| 国产亚洲精品久久久久5区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美激情久久久久久爽电影 | av天堂久久9| 久久久精品区二区三区| 亚洲视频免费观看视频| 波多野结衣一区麻豆| 日本av手机在线免费观看| 亚洲全国av大片| 免费观看av网站的网址| 欧美另类一区| 激情视频va一区二区三区| 美女脱内裤让男人舔精品视频| 狂野欧美激情性xxxx| 国产高清视频在线播放一区 | 国产成人av激情在线播放| 国产欧美日韩一区二区三 | av在线app专区| 热99re8久久精品国产| 成人av一区二区三区在线看 | 新久久久久国产一级毛片| 热99re8久久精品国产| 国产成人欧美| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 波多野结衣一区麻豆| 一区二区三区激情视频| 两性夫妻黄色片| 天堂俺去俺来也www色官网| 老鸭窝网址在线观看| 热99国产精品久久久久久7| 日韩视频一区二区在线观看| 亚洲av片天天在线观看| 99国产综合亚洲精品| 99久久99久久久精品蜜桃| 日韩三级视频一区二区三区| 啦啦啦视频在线资源免费观看| 日本精品一区二区三区蜜桃| 久久影院123| 蜜桃在线观看..| 国产片内射在线| 国产淫语在线视频| 欧美另类亚洲清纯唯美| 青春草视频在线免费观看| 老司机福利观看| 极品少妇高潮喷水抽搐| 一本色道久久久久久精品综合| videos熟女内射| 黄色片一级片一级黄色片| 成人av一区二区三区在线看 | 天天躁狠狠躁夜夜躁狠狠躁| 首页视频小说图片口味搜索| 天天躁狠狠躁夜夜躁狠狠躁| av网站在线播放免费| 国产一区有黄有色的免费视频| av片东京热男人的天堂| 国产精品久久久av美女十八| 午夜久久久在线观看| 老汉色∧v一级毛片| 国产一区二区三区av在线| 不卡一级毛片| 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美| 国产麻豆69| 亚洲少妇的诱惑av| 婷婷丁香在线五月| 精品人妻在线不人妻| 欧美另类一区| 欧美午夜高清在线| av天堂久久9| 亚洲,欧美精品.| 中文字幕最新亚洲高清| 亚洲久久久国产精品| 黄片播放在线免费| 国产在线视频一区二区| 国产成人精品久久二区二区免费| 日韩欧美免费精品| 一本综合久久免费| 国产成人精品无人区| 欧美日韩av久久| 国产精品久久久久久人妻精品电影 | 97精品久久久久久久久久精品| 亚洲精品自拍成人| 狂野欧美激情性bbbbbb| 中国美女看黄片| 少妇 在线观看| 亚洲人成电影观看| 亚洲国产成人一精品久久久| 少妇人妻久久综合中文| 亚洲精品一区蜜桃| 日韩大片免费观看网站| av欧美777| 午夜激情久久久久久久| 久久精品aⅴ一区二区三区四区| 精品一区在线观看国产| 18禁国产床啪视频网站| 精品国内亚洲2022精品成人 | 国产高清视频在线播放一区 | 国产精品香港三级国产av潘金莲| 亚洲欧洲精品一区二区精品久久久| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆| 男人舔女人的私密视频| 曰老女人黄片| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 亚洲国产欧美网| 国产亚洲精品一区二区www | 又大又爽又粗| 国产精品久久久久久精品电影小说| 国产精品一区二区在线观看99| netflix在线观看网站| 亚洲va日本ⅴa欧美va伊人久久 | 久久性视频一级片| 国产免费一区二区三区四区乱码| 80岁老熟妇乱子伦牲交| 亚洲av男天堂| 黄网站色视频无遮挡免费观看| 成在线人永久免费视频| 久久ye,这里只有精品| 精品人妻在线不人妻| 亚洲欧洲日产国产| 亚洲国产看品久久| 精品国产乱码久久久久久小说| 好男人电影高清在线观看| e午夜精品久久久久久久| 欧美在线黄色| 黑人欧美特级aaaaaa片| 91成年电影在线观看| 亚洲自偷自拍图片 自拍| 电影成人av| 精品国产超薄肉色丝袜足j| 成人黄色视频免费在线看| 国产高清视频在线播放一区 | av视频免费观看在线观看| 国产成人av激情在线播放| 成人免费观看视频高清| 色婷婷久久久亚洲欧美| 久久精品亚洲av国产电影网| 久久久久久久大尺度免费视频| 久久久精品国产亚洲av高清涩受| 国产真人三级小视频在线观看| 日本猛色少妇xxxxx猛交久久| 超色免费av| 亚洲精品第二区| 午夜精品久久久久久毛片777| 日韩制服丝袜自拍偷拍| 2018国产大陆天天弄谢| 精品乱码久久久久久99久播| 亚洲免费av在线视频| 自线自在国产av| 少妇粗大呻吟视频| www.av在线官网国产| 伦理电影免费视频| 久久狼人影院| 99久久精品国产亚洲精品| 欧美成人午夜精品| 视频在线观看一区二区三区| 在线亚洲精品国产二区图片欧美| 久久久久国产一级毛片高清牌| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 欧美激情久久久久久爽电影 | 色婷婷av一区二区三区视频| 国产真人三级小视频在线观看| 青春草亚洲视频在线观看| 性高湖久久久久久久久免费观看| 男男h啪啪无遮挡| 免费观看人在逋| 欧美精品亚洲一区二区| 两性夫妻黄色片| 国产在线一区二区三区精| 国产成人欧美在线观看 | 99久久人妻综合| 亚洲综合色网址| 蜜桃在线观看..| 亚洲中文日韩欧美视频| 大型av网站在线播放| 不卡一级毛片| 一本—道久久a久久精品蜜桃钙片| avwww免费| 久久综合国产亚洲精品| 90打野战视频偷拍视频| 国产精品99久久99久久久不卡| 香蕉国产在线看| 国产日韩欧美亚洲二区| 男女国产视频网站| 久久久国产精品麻豆| 夜夜骑夜夜射夜夜干| 亚洲一区二区三区欧美精品| 老汉色∧v一级毛片| 亚洲伊人久久精品综合| 日本欧美视频一区| 一级毛片电影观看| 亚洲成人国产一区在线观看| 欧美国产精品一级二级三级| 一个人免费看片子| 男人操女人黄网站| 曰老女人黄片| 天天操日日干夜夜撸| a 毛片基地| 69av精品久久久久久 | 黑人操中国人逼视频| 美女脱内裤让男人舔精品视频| 女人爽到高潮嗷嗷叫在线视频| 99国产综合亚洲精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲av美国av| 多毛熟女@视频| 久久狼人影院| 久久久久视频综合| 三上悠亚av全集在线观看| 亚洲人成电影观看| 老司机在亚洲福利影院| 精品免费久久久久久久清纯 | 国产有黄有色有爽视频| 人妻人人澡人人爽人人| 大陆偷拍与自拍| 午夜日韩欧美国产| 国产免费现黄频在线看| 91老司机精品| 在线观看一区二区三区激情| 十分钟在线观看高清视频www| 免费黄频网站在线观看国产| 亚洲精品国产一区二区精华液| 亚洲视频免费观看视频| av天堂久久9| 欧美激情极品国产一区二区三区| 曰老女人黄片| 黄色片一级片一级黄色片| 两个人免费观看高清视频| 精品高清国产在线一区| 免费看十八禁软件| 在线观看人妻少妇| 黄色a级毛片大全视频| 免费高清在线观看日韩| 亚洲精品第二区| 免费不卡黄色视频| 黄片小视频在线播放| 搡老熟女国产l中国老女人| 桃红色精品国产亚洲av| 国产在线观看jvid| 日本av手机在线免费观看| 国产免费福利视频在线观看| 少妇粗大呻吟视频| 老司机午夜福利在线观看视频 | 国产精品久久久av美女十八| 亚洲欧洲日产国产| 国产成人精品在线电影| 国产精品国产三级国产专区5o| 国产极品粉嫩免费观看在线| 国产亚洲精品第一综合不卡| 美女扒开内裤让男人捅视频| 久久久国产精品麻豆| 国产视频一区二区在线看| 久久久久久亚洲精品国产蜜桃av| 亚洲,欧美精品.| 热99国产精品久久久久久7| 99香蕉大伊视频| 国产免费av片在线观看野外av| 咕卡用的链子| 久久久久国内视频| 国产精品久久久av美女十八| 中文字幕av电影在线播放| 亚洲精品自拍成人| 午夜福利视频在线观看免费| 久久青草综合色| av网站免费在线观看视频| 国产深夜福利视频在线观看| 久久久久久亚洲精品国产蜜桃av| 中文字幕人妻丝袜一区二区| 男人爽女人下面视频在线观看| 热99re8久久精品国产| 91麻豆精品激情在线观看国产 | 天天添夜夜摸| 另类精品久久| 色综合欧美亚洲国产小说| 久久久国产精品麻豆| 午夜两性在线视频| e午夜精品久久久久久久| 捣出白浆h1v1| www日本在线高清视频| 久久久久久久久久久久大奶| 黄网站色视频无遮挡免费观看| 国产一区二区三区av在线| 国产精品免费大片|