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

    多目標(biāo)無序量測(cè)OOSM-GMPHD濾波算法

    2017-11-23 02:10:26胡建旺
    傳感器與微系統(tǒng) 2017年11期
    關(guān)鍵詞:無序高斯時(shí)刻

    趙 凱, 胡建旺, 吉 兵, 劉 鋼

    (1.軍械工程學(xué)院 信息工程系,河北 石家莊 050003; 2.71697部隊(duì),河南 新鄉(xiāng)453600)

    多目標(biāo)無序量測(cè)OOSM-GMPHD濾波算法

    趙 凱1, 胡建旺1, 吉 兵1, 劉 鋼2

    (1.軍械工程學(xué)院信息工程系,河北石家莊050003;2.71697部隊(duì),河南新鄉(xiāng)453600)

    針對(duì)無序量測(cè)條件下多目標(biāo)跟蹤問題,提出了一種適用于線性系統(tǒng)的單步滯后無序量測(cè)濾波算法(OOSM-GMPHD)。在前向預(yù)測(cè)框架內(nèi),以高斯混合概率假設(shè)密度(GMPHD)濾波器為基礎(chǔ)濾波算法,對(duì)每一高斯分量分別用延遲到達(dá)的量測(cè)與等價(jià)量測(cè)進(jìn)行預(yù)測(cè)、更新,經(jīng)剪枝與合并等步驟獲得最終的目標(biāo)數(shù)量與狀態(tài)估計(jì)。仿真結(jié)果表明:算法可有效消除無序量測(cè)的影響,準(zhǔn)確估計(jì)多目標(biāo)數(shù)目和狀態(tài)。

    多目標(biāo)跟蹤; 無序量測(cè); 單步滯后; 高斯混合; 前向預(yù)測(cè)

    0 引 言

    多目標(biāo)跟蹤技術(shù)一直是跟蹤領(lǐng)域研究的熱點(diǎn),被廣泛應(yīng)用于軍事、民用領(lǐng)域[1~4]。然而在多目標(biāo)跟蹤中,存在著目標(biāo)數(shù)目未知且隨時(shí)間變化、探測(cè)概率未知、虛警概率大、目標(biāo)雜波多等不利因素,嚴(yán)重影響到多目標(biāo)的有效跟蹤。傳統(tǒng)的多目標(biāo)跟蹤技術(shù)主要采用數(shù)據(jù)關(guān)聯(lián)技術(shù),需要建立“量測(cè)—航跡”對(duì)應(yīng)關(guān)系,從而使多目標(biāo)跟蹤問題轉(zhuǎn)化為多個(gè)簡單的單目標(biāo)跟蹤問題。但當(dāng)目標(biāo)數(shù)量過大、雜波較多時(shí),存在NP難問題。

    為避免復(fù)雜的數(shù)據(jù)關(guān)聯(lián)難題,文獻(xiàn)[5]提出了基于隨機(jī)有限集理論的概率假設(shè)密度(probability hypothesis density,PHD)濾波,但PHD求解過程中存在著復(fù)雜的多重積分運(yùn)算,難以計(jì)算。因此,文獻(xiàn)[6]給出了PHD濾波器的序貫蒙特卡羅(sequential Monte Carlo,SMC)實(shí)現(xiàn);文獻(xiàn)[7]給出了PHD濾波器的高斯混合(Gaussian mixture,GM)實(shí)現(xiàn)方式。由于GM實(shí)現(xiàn)方式的計(jì)算量相對(duì)較小且精度較高,更易于實(shí)現(xiàn)狀態(tài)與誤差的估計(jì),因此,其被廣泛用于近似化求解PHD濾波。

    在多目標(biāo)跟蹤系統(tǒng)中,融合中心須通過網(wǎng)絡(luò)傳輸接收各傳感器量測(cè)進(jìn)行處理[8]。由于各傳感器采樣頻率、預(yù)處理時(shí)間、通信延遲不同,各量測(cè)傳輸?shù)饺诤现行臅r(shí),常會(huì)產(chǎn)生無序量測(cè)(out-of-sequence measurement,OOSM)問題[9~11],即同一目標(biāo)的較早時(shí)刻的量測(cè)在較晚時(shí)刻的量測(cè)之后到達(dá),然而現(xiàn)有的PHD濾波無法直接處理此類負(fù)時(shí)間更新問題。文獻(xiàn)[12]提出了處理OOSM的前向預(yù)測(cè)方法,并指出該方法與過程噪聲模型無關(guān),精度較高,且計(jì)算量與儲(chǔ)存量較小,利用滯后的量測(cè)數(shù)據(jù)對(duì)最新時(shí)刻的目標(biāo)狀態(tài)進(jìn)行再更新,以獲得更精確的狀態(tài)估計(jì)與協(xié)方差矩陣。這一方法為研究OOSM問題提供了框架。

    對(duì)此,在前向預(yù)測(cè)框架內(nèi)提出了OOSM-GMPHD算法,可處理單個(gè)量測(cè)單步滯后情形下的多目標(biāo)跟蹤問題。

    1 無線射頻系統(tǒng)與PHD濾波原理

    1)k時(shí)刻多目標(biāo)狀態(tài)模型用RFS描述為

    Xk=Sk|k-1(Xk-1)∪Bk|k-1(Xk-1)∪Γk

    (1)

    式中Xk-1為k-1時(shí)刻的目標(biāo)狀態(tài)RFS;Sk|k-1為從k-1~k時(shí)刻存活目標(biāo)狀態(tài)RFS;Bk|k-1為k-1~k時(shí)刻由Xk-1衍生(Spawn Birth)的目標(biāo)狀態(tài)RFS;Γk為k時(shí)刻新生(New Birth)目標(biāo)狀態(tài)RFS。

    2)k時(shí)刻目標(biāo)量測(cè)模型用RFS描述為

    Zk=Ck(Xk-1)∪Θk(Xk-1)

    (2)

    式中Ck為雜波與虛警量測(cè)RFS;Θk為Xk的量測(cè)RFS。

    根據(jù)經(jīng)典Bayes估計(jì)理論,基于RFS的多目標(biāo)Bayes估計(jì)表示為

    pk|k-1(Xk|Z1:k-1)=

    (3)

    pk(Xk|Z1:k)=

    (4)

    式中pk|k-1與pk為多目標(biāo)先驗(yàn)、后驗(yàn)概率;fk|k-1為狀態(tài)轉(zhuǎn)移概率;gk為量測(cè)似然函數(shù);μ為某一空間近似Lebesgue測(cè)度。

    由于Bayes估計(jì)中積分運(yùn)算的存在,且量測(cè)似然函數(shù)gk的計(jì)算量隨著目標(biāo)數(shù)目增加而呈指數(shù)級(jí)增長,因此直接用其對(duì)目標(biāo)狀態(tài)與數(shù)目進(jìn)行估計(jì)十分困難。Ronald Mahler利用多目標(biāo)后驗(yàn)分布的概率假設(shè)密度(PHD)來近似后驗(yàn)概率密度,極大減小了計(jì)算量[13]。

    2 多目標(biāo)OOSM問題描述

    假定目標(biāo)i的狀態(tài)向量為xi,量測(cè)向量為Zi

    (5)

    (6)

    假設(shè)tk時(shí)刻,已有目標(biāo)i的狀態(tài)估計(jì)

    (7)

    來自t=td時(shí)刻的量測(cè)

    (8)

    圖1 單步滯后OOSM情形

    3 基于前向預(yù)測(cè)的OOSM-GMPHD濾波算法

    3.1 GMPHD多目標(biāo)濾波算法

    在應(yīng)用GMPHD濾波算法時(shí),除PHD濾波器常用假設(shè)之外,還需要滿足[14]:

    1)單目標(biāo)馬爾科夫轉(zhuǎn)移密度與似然函數(shù)均為線性高斯;

    2)目標(biāo)存活概率PS和為常數(shù);檢測(cè)概率PD為常數(shù);

    3)目標(biāo)新生和衍生過程的PHD均為高斯混合形式。

    GMPHD濾波器大致可分為預(yù)測(cè)與更新2步進(jìn)行。假設(shè)已知tk-1時(shí)刻的目標(biāo)后驗(yàn)強(qiáng)度Dk-1(x),且其形式滿足

    (9)

    式中Jk-1為高斯混合分量個(gè)數(shù)。

    a.預(yù)測(cè)步

    tk時(shí)刻預(yù)測(cè)強(qiáng)度為

    Dk|k-1(x)=γk(x)+Ds,k|k-1(x)+Db,k|k-1(x)

    (10)

    式中γk(x),Ds,k|k-1(x)與Db,k|k-1(x)分別為新生目標(biāo)密度,存活目標(biāo)和衍生目標(biāo)的PHD,即

    (11)

    (12)

    (13)

    各分量具體形式可表示為

    更進(jìn)一步地,預(yù)測(cè)步PHD可表示為

    (14)

    式中

    Jk|k-1=Jb,k|k-1+Jγ,k+Jk-1

    (15)

    b.更新步

    經(jīng)過更新,tk時(shí)刻后驗(yàn)強(qiáng)度為

    (16)

    (17)

    由此可知,經(jīng)過PHD更新后,k時(shí)刻的高斯分量個(gè)數(shù)為Jk=(1+|Zk|)Jk|k-1(|Zk|為集合Zk的勢(shì),即元素的個(gè)數(shù))。這樣,后驗(yàn)概率密度的高斯項(xiàng)隨時(shí)間變化將會(huì)無限制地增加,為了解決這一問題,需要采用剪枝、合并等措施將其控制在一定的數(shù)目范圍內(nèi)。

    3.2 等價(jià)量測(cè)

    等價(jià)量測(cè)方法在基于卡爾曼濾波與粒子濾波的OOSM更新算法中應(yīng)用較多,其原理是將該OOSM發(fā)生時(shí)刻之后的所有量測(cè)用一個(gè)等價(jià)量測(cè)代替,用來更新到最新時(shí)刻的狀態(tài)估計(jì),解決多步滯后OOSM問題,而不儲(chǔ)存該段時(shí)間內(nèi)的量測(cè)數(shù)據(jù)、狀態(tài)估計(jì)與協(xié)方差矩陣。

    (18)

    式中

    (19)

    (20)

    (21)

    (22)

    (23)

    3.3 前向預(yù)測(cè)框架內(nèi)的OOSM-GMPHD濾波算法

    基于前向預(yù)測(cè)方法的基本思想為:當(dāng)系統(tǒng)接收到OOSM時(shí),將濾波過程返回其產(chǎn)生時(shí)刻的前一時(shí)刻,順序地用該OOSM與晚于其產(chǎn)生時(shí)刻的等價(jià)量測(cè)分別進(jìn)行濾波更新,得到最新時(shí)刻的狀態(tài)估計(jì)。具體過程如圖2所示。在接收到OOSM后,根據(jù)時(shí)戳,得到該OOSM在量測(cè)序列中的位置,按照?qǐng)D2中①,②,③所示流程進(jìn)行運(yùn)算。

    圖2 前向預(yù)測(cè)方法流程

    從實(shí)質(zhì)上看,前向預(yù)測(cè)法是一種不完全的順序更新,在濾波精度與計(jì)算量之間取得了平衡。這一過程還可表述為

    (24)

    為解決線性系統(tǒng)無序量測(cè)條件下的GMPHD濾波問題,提出基于前向預(yù)測(cè)框架的OOSM-GMPHD濾波算法,將最新時(shí)刻每一個(gè)高斯分量納入前向預(yù)測(cè)框架內(nèi)。融合中心須將tk-1時(shí)刻的PHDDk-1(x)進(jìn)行儲(chǔ)存,當(dāng)tk時(shí)刻的系統(tǒng)得到后驗(yàn)PHDDk(x)后,收到滯后到達(dá)的OOSMZd,則須將濾波過程返回至tk-1時(shí)刻,具體處理過程如下:

    1)預(yù)測(cè)并更新tk-1~td時(shí)刻的分量狀態(tài)。利用OOSMZd對(duì)tk-1時(shí)刻的每一個(gè)高斯分量進(jìn)行濾波,得到td時(shí)刻的PHDDd(x)。具體過程如式(13)~式(20)所示。

    2)對(duì)td時(shí)刻的PHDDd(x)進(jìn)行剪枝與合并處理。設(shè)定一個(gè)權(quán)值門限Tprun,權(quán)重小于合并距離的高斯項(xiàng)將會(huì)被剔除,該步驟稱為剪枝。Clark,Panta和Vo提出采用式(25)定義的合并距離[14]

    (25)

    另外設(shè)定門限Tmerg,當(dāng)兩個(gè)分量合并距離小于此門限時(shí),將會(huì)合并成單個(gè)分量。經(jīng)過此步,高斯分量數(shù)量會(huì)有所減少。

    3)提取td時(shí)刻目標(biāo)數(shù)量與狀態(tài)。將Jk個(gè)分量權(quán)值求和,即得到目標(biāo)數(shù)量;權(quán)值大于0.5的高斯分量的均值對(duì)應(yīng)為目標(biāo)的狀態(tài)。其中,目標(biāo)數(shù)量可表示為

    (26)

    5)對(duì)Dk(x)進(jìn)行剪枝與合并、狀態(tài)提取,得到最新時(shí)刻的目標(biāo)數(shù)量與狀態(tài)。

    4 仿真驗(yàn)證

    以二維平面內(nèi)多目標(biāo)運(yùn)動(dòng)情形為例,通過仿真實(shí)驗(yàn)綜合驗(yàn)證OOSM-GMPHD算法在無序條件下的跟蹤性能,并忽略O(shè)OSM的處理方式,即丟棄滯后量測(cè)法與真實(shí)值做對(duì)比,以驗(yàn)證算法的有效性。

    4.1 參數(shù)設(shè)置

    假設(shè)目標(biāo)的監(jiān)測(cè)區(qū)域?yàn)閇-100 m,100 m]×[-100 m,100 m]的二維平面,對(duì)應(yīng)的目標(biāo)狀態(tài)方程為

    (27)

    空間中有4個(gè)目標(biāo)在運(yùn)動(dòng),相繼出現(xiàn)或消失。傳感器對(duì)4個(gè)目標(biāo)進(jìn)行觀測(cè),得到的量測(cè)為其位置的二維坐標(biāo)。假設(shè)每次掃描的雜波數(shù)服從泊松分布,且該過程中忽略衍生目標(biāo)出現(xiàn),則傳感器量測(cè)方程可表示為

    (28)

    4.2 仿真一

    在跟蹤過程的前7 ,傳感器獲得7個(gè)量測(cè)。融合中心最終收到的量測(cè)序列為{z1,z2,z3,z5,z4,z6,z7}。比較為t=6時(shí)刻的位置分量均方根誤差(RMSE)及目標(biāo)數(shù)量。

    表1 仿真一結(jié)果

    分別用OOSM-GMPHD算法、丟棄滯后量測(cè)法進(jìn)行處理,將得到的結(jié)果與真實(shí)值做比較,結(jié)果如表1所示。可以看出,OOSM-GMPHD的RMSE遠(yuǎn)低于丟棄滯后量測(cè)法,即跟蹤精度較高;且丟棄滯后量測(cè)法在目標(biāo)數(shù)量估計(jì)上出現(xiàn)了偏差,此項(xiàng)的RMSE數(shù)值忽略了誤差較大的目標(biāo)。

    4.3 仿真二

    現(xiàn)將跟蹤時(shí)間延長至35,即融合中心收到35個(gè)量測(cè)。同樣假設(shè)每7個(gè)量測(cè)中包含一個(gè)單步滯后的OOSM。分別使用OOSM-GMPHD算法與丟棄滯后量測(cè)法處理上述量測(cè),將濾波結(jié)果與目標(biāo)真實(shí)運(yùn)動(dòng)軌跡、數(shù)量做比較,結(jié)果如圖3、圖4所示。

    圖3 多目標(biāo)真實(shí)航跡與不同算法狀態(tài)估計(jì)

    圖4 不同算法多目標(biāo)數(shù)目

    圖3為多目標(biāo)真實(shí)航跡與不同算法狀態(tài)估計(jì)的對(duì)比,可以看出:盡管多目標(biāo)的量測(cè)數(shù)據(jù)中含有多個(gè)OOSM,但在使用OOSM-GMPHD算法濾波后,這些OOSM被應(yīng)用到了狀態(tài)更新中,得到的多目標(biāo)的狀態(tài)估計(jì)量可以很好地跟蹤目標(biāo)的真實(shí)運(yùn)動(dòng)狀態(tài),表明了該算法的有效性 。

    圖4給出了在不同時(shí)刻量測(cè)數(shù)據(jù)經(jīng)過OOSM-GMPHD算法處理得到的目標(biāo)數(shù)量與丟棄滯后量測(cè)法及真實(shí)值的對(duì)比,算法對(duì)目標(biāo)數(shù)目的估計(jì)基本接近真實(shí)值,優(yōu)于丟棄滯后量測(cè)法,目標(biāo)的數(shù)目估計(jì)更加準(zhǔn)確。

    5 結(jié) 論

    對(duì)于無序量測(cè)條件下的多目標(biāo)跟蹤問題,給出的OOSM-GMPHD濾波算法能夠有效地跟蹤多個(gè)目標(biāo),可對(duì)目標(biāo)的數(shù)目及狀態(tài)做出有效估計(jì)。算法在前向預(yù)測(cè)框架內(nèi),對(duì)每一個(gè)高斯分量分別使用OOSM與等價(jià)量測(cè)進(jìn)行更新,而后利用剪枝與合并算法,提取出目標(biāo)數(shù)量與狀態(tài)估計(jì)。仿真結(jié)果表明:該算法能準(zhǔn)確跟蹤無序量測(cè)條件下的多個(gè)目標(biāo)。

    [1] Bar-Shalom Y,Fortmann T E.Tracking and data association[M].San Diego:Academic,1988.

    [2] 郝燕玲,孟凡彬,張崇猛,等.多傳感器多目標(biāo)跟蹤的粒子PHD濾波算法[J].傳感器與微系統(tǒng),2010,29(4):18-21.

    [3] 趙 斌,胡建旺,吉 兵.改進(jìn)匹配方法的BFG-GMPHD濾波算法[J].傳感器與微系統(tǒng),2016,35(7):136-139.

    [4] 劉麗娟,劉國棟.基于ET-GM-PHD的多傳感器多目標(biāo)跟蹤算法[J].傳感器與微系統(tǒng),2013,32(10):126-128.

    [5] Mahler R P S.Multitarget Bayes filtering via first-order multitarget moments[J].IEEE Transaction on Aerospace and Electronic Systems,2003,39(4):1152-1178.

    [6] Vo B N,Singh S,Doucet A.Sequential Monte Carlo methods for multi-target filtering with random finite sets[J].IEEE Transactions on Aerospace and Electronic Systems,2005,41(4):1224-1245.

    [7] Vo B N,Ma W K.The Gaussian mixture probability hypothesis density filter[J].IEEE Transactions on Signal Processing,2006,54(11):4091-4104.

    [8] 袁 丁,胡建旺,吉 兵,等.快速邊緣粒子濾波在無序量測(cè)問題中的應(yīng)用[J].傳感器與微系統(tǒng),2014,33(6):157-160.

    [9] 韓崇昭,朱紅艷,段戰(zhàn)勝,等.多源信息融合 [M].2版.北京:清華大學(xué)出版社,2010:368-379.

    [10] 王洪峰,周 磊,單甘霖.國外軍事信息融合理論與應(yīng)用的研究進(jìn)展[J].電光與控制,2007,14(4):13-17.

    [11] 吳衛(wèi)華,江 晶.無序量測(cè)問題研究綜述[J].傳感器與微系統(tǒng),2014,33(12):5-9.

    [12] 周文輝,李 琳,陳國海,等.單步延遲無序量測(cè)濾波算法的最優(yōu)性分析[J].中國科學(xué),2007,37(4):564-580.

    [13] Mahler R.A theory of PHD filters of higher order in target numb-r[C]∥Proceedings of SPIE—The International Society for Optical Engineering,2006:62350K—1—62350K-12.

    [14] Daniel Clark,Ba-Ngu Vo.Convergence analysis of the Gaussian mixture PHD filter[J].IEEE Transactions on Singnal Processing,2007,55(4):1204-1212

    Multi-targetOOSM-GMPHDfilteringalgorithmwithout-of-sequencemeasurement

    ZHAO Kai1, HU Jian-wang1, JI Bing, LIU Gang2

    (1.DepartmentofInformationEngineering,OrdnanceEngineeringCollege,Shijiazhuang050003,China;2.Unit71697ofPLA,Xinxiang453600,China)

    For addressing multi-target tracking problem with out-of-sequence measurement(OOSM),a filtering algorithm is proposed with multiple single-step-lag OOSM for linear system,which denotes as OOSM-Gaussian mixture probability hypothesis density(GMPHD).Within the forward prediction framework,taking GMPHD as basis filtering algorithm,predicting and updating each Gaussian component,and then can obtained the target number and state estimation after pruning and merging.etc.Simulation results show that the algorithm can effectively filter out influence of OOSM and can accurately estimate the multi-target number and state.

    multi-target tracking; out-of-sequence measurement(OOSM); one-step-lag; Gaussian mixture(GM); forward prediction

    10.13873/J.1000—9787(2017)11—0154—04

    TP 391

    A

    1000—9787(2017)11—0154—04

    2016—10—27

    趙 凱(1991-),男,碩士研究生,主要研究方向?yàn)槟繕?biāo)跟蹤與多源信息融合,E—mail:609587194@qq.com。

    猜你喜歡
    無序高斯時(shí)刻
    小高斯的大發(fā)現(xiàn)
    車身無序堆疊零件自動(dòng)抓取系統(tǒng)
    冬“傲”時(shí)刻
    捕獵時(shí)刻
    天才數(shù)學(xué)家——高斯
    張博庭:煤電不能再這么無序發(fā)展下去了
    能源(2017年11期)2017-12-13 08:12:30
    高速路上右行規(guī)則與無序行駛規(guī)則的比較研究
    無序體系中的國際秩序
    街拍的歡樂時(shí)刻到來了
    有限域上高斯正規(guī)基的一個(gè)注記
    亚洲五月天丁香| 9191精品国产免费久久| 欧美日韩黄片免| 最新中文字幕久久久久 | 国产成人啪精品午夜网站| 香蕉久久夜色| 大型黄色视频在线免费观看| 淫秽高清视频在线观看| 精品一区二区三区四区五区乱码| 国产精品电影一区二区三区| 国内少妇人妻偷人精品xxx网站 | 国产成人精品久久二区二区免费| 亚洲av成人不卡在线观看播放网| 久久久国产精品麻豆| 国产欧美日韩一区二区精品| 亚洲在线观看片| 国产精品1区2区在线观看.| 97碰自拍视频| 国产一区二区在线av高清观看| 禁无遮挡网站| 制服人妻中文乱码| 国语自产精品视频在线第100页| 一级a爱片免费观看的视频| 国产精品av久久久久免费| www.www免费av| 亚洲人成伊人成综合网2020| 亚洲国产色片| a在线观看视频网站| x7x7x7水蜜桃| 欧美黄色淫秽网站| 国产精品 欧美亚洲| svipshipincom国产片| 看片在线看免费视频| 看片在线看免费视频| 在线观看舔阴道视频| 超碰成人久久| av天堂中文字幕网| 国产野战对白在线观看| 九九久久精品国产亚洲av麻豆 | 51午夜福利影视在线观看| 免费看十八禁软件| av女优亚洲男人天堂 | 麻豆成人av在线观看| 一进一出抽搐gif免费好疼| 99热只有精品国产| 一级作爱视频免费观看| 亚洲午夜理论影院| 日韩欧美精品v在线| 亚洲九九香蕉| 人妻久久中文字幕网| www日本在线高清视频| 天天躁狠狠躁夜夜躁狠狠躁| cao死你这个sao货| 亚洲欧洲精品一区二区精品久久久| 婷婷精品国产亚洲av| 亚洲片人在线观看| 这个男人来自地球电影免费观看| 高潮久久久久久久久久久不卡| 国产极品精品免费视频能看的| 美女黄网站色视频| 国产成人aa在线观看| 美女被艹到高潮喷水动态| 欧美一区二区精品小视频在线| 亚洲熟妇熟女久久| 国产日本99.免费观看| 中文资源天堂在线| 国产不卡一卡二| 欧美zozozo另类| 精品99又大又爽又粗少妇毛片 | 麻豆国产av国片精品| 淫妇啪啪啪对白视频| 亚洲精品在线美女| 午夜福利免费观看在线| 中出人妻视频一区二区| 国产亚洲av高清不卡| 色吧在线观看| 成人av在线播放网站| 国内精品久久久久精免费| 91在线观看av| 国产极品精品免费视频能看的| 国产三级中文精品| 国产精品久久久av美女十八| 日韩欧美三级三区| 极品教师在线免费播放| 久久天躁狠狠躁夜夜2o2o| bbb黄色大片| 法律面前人人平等表现在哪些方面| 欧美最黄视频在线播放免费| 白带黄色成豆腐渣| 亚洲专区中文字幕在线| 国产精品久久久久久精品电影| 中文字幕熟女人妻在线| 真人一进一出gif抽搐免费| 午夜成年电影在线免费观看| 1000部很黄的大片| 极品教师在线免费播放| 亚洲第一电影网av| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 亚洲精品456在线播放app | 一本一本综合久久| 99久久99久久久精品蜜桃| 欧美日本亚洲视频在线播放| 久久亚洲真实| 怎么达到女性高潮| 天天躁狠狠躁夜夜躁狠狠躁| 中亚洲国语对白在线视频| 中文在线观看免费www的网站| 久久精品91蜜桃| 国产av不卡久久| 99久久国产精品久久久| 又黄又爽又免费观看的视频| 久久久成人免费电影| 久久久久九九精品影院| 中文字幕av在线有码专区| 女人高潮潮喷娇喘18禁视频| 757午夜福利合集在线观看| 婷婷精品国产亚洲av在线| av在线天堂中文字幕| 精品久久蜜臀av无| 手机成人av网站| 久久天堂一区二区三区四区| 99久久国产精品久久久| 成人国产综合亚洲| 国产精品久久久久久精品电影| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 女人高潮潮喷娇喘18禁视频| 国产成人精品无人区| 欧美中文综合在线视频| 国产一级毛片七仙女欲春2| 日韩欧美一区二区三区在线观看| 亚洲欧美一区二区三区黑人| www.精华液| 久久久久久久久中文| 波多野结衣高清无吗| 91在线精品国自产拍蜜月 | 国产视频一区二区在线看| 亚洲国产欧洲综合997久久,| 欧美中文日本在线观看视频| 成人特级黄色片久久久久久久| 亚洲精品国产精品久久久不卡| 97超级碰碰碰精品色视频在线观看| 在线十欧美十亚洲十日本专区| 九九久久精品国产亚洲av麻豆 | 久久国产精品影院| 亚洲精品美女久久av网站| 99在线视频只有这里精品首页| 国产亚洲精品av在线| 中文字幕人成人乱码亚洲影| 亚洲国产色片| 国产精品亚洲美女久久久| 国产精品香港三级国产av潘金莲| 在线观看日韩欧美| 免费看光身美女| 亚洲片人在线观看| 国产精品久久久久久人妻精品电影| 国产69精品久久久久777片 | 国产成年人精品一区二区| 岛国在线免费视频观看| 亚洲欧美激情综合另类| 亚洲 欧美一区二区三区| 国产亚洲精品久久久久久毛片| 91av网一区二区| 亚洲av日韩精品久久久久久密| 禁无遮挡网站| 99久久成人亚洲精品观看| 国产精品亚洲美女久久久| 日韩人妻高清精品专区| 亚洲av美国av| 一本久久中文字幕| 欧美另类亚洲清纯唯美| 少妇熟女aⅴ在线视频| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 精品久久久久久成人av| 亚洲欧美精品综合一区二区三区| 久久这里只有精品中国| 国产亚洲欧美在线一区二区| 91麻豆精品激情在线观看国产| 九九在线视频观看精品| 欧美黄色片欧美黄色片| 99久久国产精品久久久| 在线十欧美十亚洲十日本专区| 欧美午夜高清在线| 1024手机看黄色片| 午夜福利成人在线免费观看| 岛国在线观看网站| 久久精品aⅴ一区二区三区四区| 男女视频在线观看网站免费| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费电影在线观看| 亚洲18禁久久av| 一本一本综合久久| 亚洲欧美日韩高清在线视频| 亚洲 国产 在线| 国产精品av久久久久免费| 熟女少妇亚洲综合色aaa.| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧洲综合997久久,| 午夜激情欧美在线| 少妇熟女aⅴ在线视频| 91麻豆精品激情在线观看国产| 国产美女午夜福利| 色老头精品视频在线观看| 18禁国产床啪视频网站| 好男人在线观看高清免费视频| 99国产精品99久久久久| av中文乱码字幕在线| 69av精品久久久久久| 久久午夜综合久久蜜桃| 国产精品综合久久久久久久免费| 亚洲av日韩精品久久久久久密| 夜夜夜夜夜久久久久| 好看av亚洲va欧美ⅴa在| 亚洲精品一区av在线观看| av在线天堂中文字幕| 99久国产av精品| av福利片在线观看| 亚洲精品色激情综合| 日本 欧美在线| 国产aⅴ精品一区二区三区波| 99久久久亚洲精品蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 午夜免费成人在线视频| 中文字幕人妻丝袜一区二区| 好男人在线观看高清免费视频| 欧美中文综合在线视频| 精品久久久久久久人妻蜜臀av| 国产黄片美女视频| 神马国产精品三级电影在线观看| 天堂动漫精品| 亚洲国产精品成人综合色| 国产美女午夜福利| 日日摸夜夜添夜夜添小说| av天堂中文字幕网| 日韩成人在线观看一区二区三区| 床上黄色一级片| 99riav亚洲国产免费| 又爽又黄无遮挡网站| 国产在线精品亚洲第一网站| 搞女人的毛片| 国产高潮美女av| 亚洲国产精品sss在线观看| 国产午夜精品久久久久久| 亚洲黑人精品在线| 国产精品一区二区免费欧美| 精品久久久久久久久久免费视频| 国产亚洲精品av在线| 日韩中文字幕欧美一区二区| 精品人妻1区二区| 久久久久久久久久黄片| 性色av乱码一区二区三区2| 国产亚洲精品久久久com| 给我免费播放毛片高清在线观看| 久久久久国产一级毛片高清牌| 日本一本二区三区精品| 免费看光身美女| 午夜免费成人在线视频| 国产成年人精品一区二区| 亚洲av熟女| 日本黄色片子视频| 精品一区二区三区av网在线观看| 亚洲国产精品合色在线| av中文乱码字幕在线| 欧美黄色片欧美黄色片| av视频在线观看入口| 国产成人一区二区三区免费视频网站| 午夜久久久久精精品| 亚洲va日本ⅴa欧美va伊人久久| 日韩免费av在线播放| 久久天堂一区二区三区四区| 午夜亚洲福利在线播放| 老司机午夜福利在线观看视频| h日本视频在线播放| 亚洲专区国产一区二区| 亚洲精品在线美女| 制服丝袜大香蕉在线| 两个人看的免费小视频| 亚洲一区二区三区色噜噜| 19禁男女啪啪无遮挡网站| 麻豆成人午夜福利视频| 真人做人爱边吃奶动态| 男人舔奶头视频| 国产黄色小视频在线观看| 亚洲国产精品久久男人天堂| 夜夜看夜夜爽夜夜摸| 色综合站精品国产| 人妻久久中文字幕网| 亚洲 欧美 日韩 在线 免费| 亚洲精品国产精品久久久不卡| 欧美黄色淫秽网站| 亚洲欧美日韩高清专用| 黄色丝袜av网址大全| 99久久久亚洲精品蜜臀av| 精品国产乱码久久久久久男人| 亚洲国产精品久久男人天堂| 国产一区二区三区视频了| 听说在线观看完整版免费高清| 午夜成年电影在线免费观看| 日本黄大片高清| 俄罗斯特黄特色一大片| 国产 一区 欧美 日韩| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av在线| 91久久精品国产一区二区成人 | 桃红色精品国产亚洲av| 国产视频一区二区在线看| 欧美激情在线99| 久久精品亚洲精品国产色婷小说| 99国产综合亚洲精品| 日日摸夜夜添夜夜添小说| 中文字幕人妻丝袜一区二区| 亚洲在线观看片| 欧美午夜高清在线| 欧美大码av| 黄色女人牲交| 国产人伦9x9x在线观看| 色吧在线观看| 国产精品一区二区免费欧美| 国产精品一及| 99久久成人亚洲精品观看| 在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产伦精品一区二区三区视频9 | 免费一级毛片在线播放高清视频| 国产三级黄色录像| 少妇的逼水好多| 亚洲av成人av| 一级毛片精品| 亚洲无线在线观看| 在线观看舔阴道视频| 淫秽高清视频在线观看| 久久精品国产综合久久久| 视频区欧美日本亚洲| 中文字幕久久专区| 床上黄色一级片| 中文字幕久久专区| 国产一区二区激情短视频| 怎么达到女性高潮| 麻豆国产97在线/欧美| 久久久精品大字幕| 亚洲精品美女久久久久99蜜臀| xxx96com| 欧美黄色片欧美黄色片| 天堂网av新在线| 久久久久久久久久黄片| 久久久精品大字幕| 麻豆国产97在线/欧美| 99热精品在线国产| 欧美性猛交黑人性爽| 欧美最黄视频在线播放免费| 哪里可以看免费的av片| 可以在线观看毛片的网站| 一二三四社区在线视频社区8| 国产在线精品亚洲第一网站| АⅤ资源中文在线天堂| 亚洲av成人精品一区久久| 97超视频在线观看视频| 国产综合懂色| 亚洲激情在线av| 国产av一区在线观看免费| 国产蜜桃级精品一区二区三区| 午夜成年电影在线免费观看| 1024香蕉在线观看| 欧美一级a爱片免费观看看| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添小说| 国产三级黄色录像| 99国产精品99久久久久| 俄罗斯特黄特色一大片| 免费在线观看视频国产中文字幕亚洲| 黄色视频,在线免费观看| 9191精品国产免费久久| 性色av乱码一区二区三区2| 日韩国内少妇激情av| 国产高清有码在线观看视频| 欧美日韩精品网址| 午夜免费激情av| 日韩中文字幕欧美一区二区| 99在线视频只有这里精品首页| 人妻夜夜爽99麻豆av| 国产真实乱freesex| 一进一出抽搐gif免费好疼| 欧美日韩一级在线毛片| 久久久久亚洲av毛片大全| 欧美色欧美亚洲另类二区| 国产野战对白在线观看| 午夜福利成人在线免费观看| 他把我摸到了高潮在线观看| 夜夜躁狠狠躁天天躁| 长腿黑丝高跟| 亚洲成av人片在线播放无| 偷拍熟女少妇极品色| 成人永久免费在线观看视频| 国产精品久久久久久人妻精品电影| 免费无遮挡裸体视频| 18美女黄网站色大片免费观看| 小说图片视频综合网站| 日本一二三区视频观看| 一区二区三区高清视频在线| 亚洲欧美激情综合另类| 99re在线观看精品视频| 午夜久久久久精精品| 99精品欧美一区二区三区四区| 免费在线观看日本一区| 九九热线精品视视频播放| 好看av亚洲va欧美ⅴa在| 一夜夜www| 亚洲色图av天堂| www.www免费av| 女同久久另类99精品国产91| 给我免费播放毛片高清在线观看| 大型黄色视频在线免费观看| 又爽又黄无遮挡网站| 一a级毛片在线观看| 精品久久久久久成人av| 欧美黄色片欧美黄色片| 成人精品一区二区免费| 国产精品爽爽va在线观看网站| 中文资源天堂在线| 啪啪无遮挡十八禁网站| 成年女人看的毛片在线观看| 热99re8久久精品国产| 2021天堂中文幕一二区在线观| 久久久久国产一级毛片高清牌| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久久久久久久| 国产成年人精品一区二区| 好男人电影高清在线观看| 免费av毛片视频| 日本黄色片子视频| 亚洲精品在线美女| а√天堂www在线а√下载| 日韩欧美国产在线观看| 国产又色又爽无遮挡免费看| 色视频www国产| 国内精品美女久久久久久| 天天添夜夜摸| 日本一本二区三区精品| 色在线成人网| 国产美女午夜福利| 中文资源天堂在线| 亚洲国产欧美人成| 99热这里只有精品一区 | 五月玫瑰六月丁香| 伦理电影免费视频| 国产乱人视频| 搡老妇女老女人老熟妇| 国产av麻豆久久久久久久| 岛国视频午夜一区免费看| 国产精品永久免费网站| 日本五十路高清| 亚洲国产色片| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 丁香六月欧美| 国产午夜精品久久久久久| 两个人看的免费小视频| 天天躁狠狠躁夜夜躁狠狠躁| 日本一二三区视频观看| 久久久久国产精品人妻aⅴ院| 亚洲熟妇中文字幕五十中出| 国产真实乱freesex| 曰老女人黄片| 日韩欧美免费精品| 国产熟女xx| 国产私拍福利视频在线观看| 国产日本99.免费观看| 国产午夜福利久久久久久| 国内毛片毛片毛片毛片毛片| 美女被艹到高潮喷水动态| 国产欧美日韩精品亚洲av| bbb黄色大片| 久久久久久人人人人人| 亚洲av第一区精品v没综合| 国产在线精品亚洲第一网站| 久久久成人免费电影| 精品久久久久久成人av| 久久人人精品亚洲av| 国产成人影院久久av| 欧美精品啪啪一区二区三区| 亚洲18禁久久av| 特大巨黑吊av在线直播| 亚洲熟女毛片儿| 岛国在线免费视频观看| www日本在线高清视频| 久久中文字幕人妻熟女| 成人一区二区视频在线观看| 午夜免费成人在线视频| 欧美色欧美亚洲另类二区| 亚洲熟女毛片儿| 男女视频在线观看网站免费| 脱女人内裤的视频| 午夜免费成人在线视频| 国产麻豆成人av免费视频| 欧美一区二区精品小视频在线| 久久九九热精品免费| 99视频精品全部免费 在线 | 成人欧美大片| 麻豆国产97在线/欧美| 国产高清三级在线| 757午夜福利合集在线观看| 国产精品自产拍在线观看55亚洲| 国产精品野战在线观看| 久久人人精品亚洲av| 日本撒尿小便嘘嘘汇集6| 男人舔奶头视频| 午夜激情欧美在线| 国产精品影院久久| 国产97色在线日韩免费| www日本黄色视频网| 欧美3d第一页| 国产真人三级小视频在线观看| 国产精品久久电影中文字幕| 久久这里只有精品中国| 99久久成人亚洲精品观看| 国产伦精品一区二区三区视频9 | 精品国产乱码久久久久久男人| 亚洲av片天天在线观看| 波多野结衣高清作品| 99精品在免费线老司机午夜| 一a级毛片在线观看| 狂野欧美白嫩少妇大欣赏| 18禁观看日本| 俺也久久电影网| 亚洲成av人片在线播放无| 舔av片在线| 国内精品久久久久精免费| 亚洲18禁久久av| 看免费av毛片| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久久久久久久| 九九久久精品国产亚洲av麻豆 | 国内久久婷婷六月综合欲色啪| 91麻豆精品激情在线观看国产| 天堂网av新在线| 老汉色av国产亚洲站长工具| 久久久精品大字幕| 久久久久久久久久黄片| 国产亚洲精品久久久com| 夜夜看夜夜爽夜夜摸| 日本 欧美在线| 久久精品91蜜桃| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 亚洲激情在线av| 99热这里只有是精品50| 99精品欧美一区二区三区四区| 久久久精品大字幕| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| av片东京热男人的天堂| 欧美又色又爽又黄视频| 麻豆国产av国片精品| 欧美日韩综合久久久久久 | 天堂√8在线中文| 99热这里只有是精品50| 变态另类丝袜制服| 91九色精品人成在线观看| 狠狠狠狠99中文字幕| 窝窝影院91人妻| 久久中文字幕人妻熟女| 色老头精品视频在线观看| 久久久国产精品麻豆| 国产精品电影一区二区三区| 国产高潮美女av| 欧美xxxx黑人xx丫x性爽| 国产乱人伦免费视频| 狠狠狠狠99中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 亚洲七黄色美女视频| 亚洲专区字幕在线| 亚洲国产精品成人综合色| 久久精品亚洲精品国产色婷小说| 亚洲七黄色美女视频| 欧美不卡视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| av视频在线观看入口| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区久久| 久久午夜综合久久蜜桃| 久久精品国产亚洲av香蕉五月| 亚洲国产精品久久男人天堂| 色噜噜av男人的天堂激情| 亚洲国产中文字幕在线视频| 午夜a级毛片| av片东京热男人的天堂| 国产爱豆传媒在线观看| 欧美日韩福利视频一区二区| 国产精品综合久久久久久久免费| 给我免费播放毛片高清在线观看| 亚洲电影在线观看av| 丝袜人妻中文字幕| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久末码| 免费av不卡在线播放| 日韩av在线大香蕉| 色综合婷婷激情| 国产精品爽爽va在线观看网站| 色综合婷婷激情| 欧美性猛交黑人性爽| 99久久99久久久精品蜜桃| 国产三级黄色录像| 小蜜桃在线观看免费完整版高清| 99在线人妻在线中文字幕| 国产精品电影一区二区三区|