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

    基于改進(jìn)對稱量測方程的多目標(biāo)跟蹤

    2016-01-21 08:00:13黃樹彩

    吳 瀟, 黃樹彩, 凌 強, 鐘 宇

    (空軍工程大學(xué)防空反導(dǎo)學(xué)院, 陜西 西安 710051)

    ?

    基于改進(jìn)對稱量測方程的多目標(biāo)跟蹤

    吳瀟, 黃樹彩, 凌強, 鐘宇

    (空軍工程大學(xué)防空反導(dǎo)學(xué)院, 陜西 西安 710051)

    摘要:采用傳統(tǒng)對稱量測方程對多維多目標(biāo)跟蹤會增加偽觀測量,在目標(biāo)航跡交叉點附近容易產(chǎn)生較大跟蹤誤差。針對這一問題,提出一種基于多項式因式分解思想的改進(jìn)對稱量測方程,通過建立同一目標(biāo)在不同坐標(biāo)系中狀態(tài)之間的聯(lián)系,將新量測集與原始量測集的對應(yīng)關(guān)系描述的更加準(zhǔn)確,利用泰勒展開公式推導(dǎo)出新量測的觀測誤差,求出相應(yīng)的觀測誤差協(xié)方差陣。與現(xiàn)有對稱量測方程方法和聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)算法仿真對比,驗證所提方法不僅延續(xù)了對稱量測方程時效性優(yōu)勢,而且提高了目標(biāo)在航跡交叉時的跟蹤精度。

    關(guān)鍵詞:對稱量測方程; 多目標(biāo)跟蹤; 多項式因式分解

    0引言

    與單目標(biāo)跟蹤相比,多目標(biāo)跟蹤不僅面臨更加復(fù)雜的噪聲背景,而且還要處理好量測-航跡關(guān)聯(lián)。傳統(tǒng)多目標(biāo)跟蹤算法[1]以數(shù)據(jù)關(guān)聯(lián)為核心,要先進(jìn)行關(guān)聯(lián)運算,然后再對多個目標(biāo)獨立跟蹤,如聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)算法(joint probability data association,JPDA)、多假設(shè)跟蹤算法(multiple hypothesis tracking,MHT)等,這些方法隨著目標(biāo)數(shù)目增多計算量會劇增,極大地制約了算法的時效性。

    20世紀(jì)90年代,Kamen[2]首創(chuàng)了一種基于對稱量測方程(symmetric measurement equation,SME)的多目標(biāo)跟蹤方法,通過構(gòu)造原始觀測值的對稱函數(shù)得到一組新量測,可以同時估計多個目標(biāo)的狀態(tài),而不需要數(shù)據(jù)關(guān)聯(lián)運算。起初,Kamen僅考慮了一維空間量測數(shù)目與目標(biāo)真實數(shù)目相等的情況,并且構(gòu)造的新量測也只有乘積和形式[2-3],隨著研究的深入,不僅將SME的適用范圍推廣到多維空間,以及存在虛警、漏警的情況[4],而且構(gòu)造了冪次和形式的SME[5]。

    在處理多維多目標(biāo)跟蹤問題時,Kamen做法[5]是各維獨立進(jìn)行SME跟蹤,Muder[6]指出該處理方法會增加偽觀測量,在航跡交叉時容易產(chǎn)生較大影響。Lee[4]在二維空間提出了一種基于復(fù)數(shù)運算的SME方法,將原始測量值的橫縱坐標(biāo)值分別表示為復(fù)數(shù)的實部和虛部,通過復(fù)數(shù)運算建立同一目標(biāo)在不同坐標(biāo)系中狀態(tài)之間的聯(lián)系,并且構(gòu)造出了基于復(fù)數(shù)運算的乘積和形式SME。Leven和Lanterman[7]在Lee研究的基礎(chǔ)上建立了基于復(fù)數(shù)運算的冪次和形式SME,并將擴展卡爾曼濾波(extended kalman filter,EKF)與無跡卡爾曼濾波(unscented kalman filter, UKF)對比,仿真發(fā)現(xiàn)UKF濾波可能會產(chǎn)生奇異矩陣,而選用EKF濾波效果更好。雖然基于復(fù)數(shù)運算的SME有效抑制了偽觀測量的產(chǎn)生,跟蹤效果與Kamen經(jīng)典的SME相比有了很大改善,但在航跡交叉點附近的跟蹤效果還有提升空間。近幾年,Baum等[8]針對SME誤差期望值和協(xié)方差難以有效計算和非線性濾波問題,提出了最優(yōu)高斯濾波器,由于考慮了新量測方程的高階項,因而能給出更優(yōu)估計效果,但該方法在目標(biāo)數(shù)目增多時也存在計算復(fù)雜度較高的缺點。

    本文利用多項式因式分解的思想構(gòu)造了一種改進(jìn)SME,在觀測誤差服從獨立同分布、零均值高斯白噪聲前提下,推導(dǎo)出相應(yīng)的觀測誤差協(xié)方差陣。與現(xiàn)有幾種典型SME和JPDA算法仿真對比,驗證新算法不僅延續(xù)了SME時效性優(yōu)勢,而且跟蹤精度基本達(dá)到與JPDA相近的水平,有效抑制了航跡交叉時跟蹤誤差較大現(xiàn)象。

    1對稱量測方程

    1.1基本概念

    (1)

    式中,T為采樣周期。

    (2)

    式中,IN為N×N的單位陣;0N為N×N的零矩陣;w(k)=[w1(k)w2(k)…wN(k)]T。

    (3)

    式中,xj(k)表示第j個目標(biāo)在x軸的坐標(biāo)值;uj(k)是第j個目標(biāo)的觀測誤差,為零均值高斯白噪聲;j=1,2,…,N。

    (4)

    式中,hi表示第i個對稱變換函數(shù);vi(k)為第i個新量測的觀測誤差。

    以一維空間N=2乘積和形式的SME為例,介紹具體構(gòu)造原理。假設(shè)獲得量測集為{z1(k),z2(k)},定義新量測

    (5)

    由式(4)可以得到

    (6)

    定義新量測矩陣Z(k)=[Z1(k),Z2(k)]T和觀測誤差矩陣V(k)=[v1(k),v2(k)]T,則觀測誤差協(xié)方差陣R(k)=E[V(k)V(k)T]。

    1.2基于多項式思想的SME

    以Kamen乘積和形式SME為例[3],在一維空間新量測可以表示為

    (7)

    式中,xi表示第i個目標(biāo)的在x軸上的坐標(biāo)值,i=1,…,N。

    這種構(gòu)造新量測的方法有以下特點:

    (1) 新量測通常具有對稱性,其對稱性表現(xiàn)為新量測與原始測量數(shù)據(jù)的輸入順序無關(guān),不會隨著原始測量數(shù)據(jù)順序的變化而發(fā)生改變。

    (2) 從原始量測變換到新量測過程中信息應(yīng)當(dāng)是完整而沒有丟失,即新量測集與原始量測集應(yīng)當(dāng)具有一一對應(yīng)關(guān)系。具體表現(xiàn)為:當(dāng)原始量測集已知時,可以確定與其對應(yīng)的新量測集;而當(dāng)新量測集作為已知條件時,原始量測集也應(yīng)當(dāng)可以求解得到。

    由式(7)可以看出,一維空間Kamen乘積和SME顯然具有對稱性,為驗證其原始量測集與新量測集具有一一對應(yīng)關(guān)系,構(gòu)造關(guān)于s的因式:pi=s+xi,則有一元多項式

    (8)

    Kamen乘積和SME即為多項式系數(shù)ai所組成的集合,由于一元多項式與其分解因式之間具有一一對應(yīng)關(guān)系[6],因此從原始測量集變換到新量測集的信息是完整而沒有丟失的。

    將SME方法推廣到二維空間時,較為簡單的處理思路是各維獨立作SME變換,與之對應(yīng)的Kamen乘積和SME可表示為

    (9)

    式中,xi、yi分別是第i個目標(biāo)的在x軸、y軸上的坐標(biāo)值,i=1,…,N。

    (10)

    當(dāng)新量測集已知時,可以得到x軸和y軸的觀測集分別為{x1,x2}和{y1,y2},但無法確定同一目標(biāo)不同坐標(biāo)之間的關(guān)系,即兩個目標(biāo)的組合形式不僅有{(x1,y1), (x2,y2)},而且還有可能存在{(x1,y2), (x2,y1)}的形式。在跟蹤過程中顯然不希望后一種組合方式出現(xiàn),而這種被稱之為“鬼點”的坐標(biāo)會對結(jié)果產(chǎn)生不良影響,尤其是在航跡交叉時表現(xiàn)更為明顯。因此,需要改進(jìn)SME方法,使這種“鬼點”坐標(biāo)越少越好。

    2改進(jìn)SME跟蹤算法

    2.1改進(jìn)SME

    在二維空間提出一種基于分解因式的改進(jìn)思路,以提高新量測集的信息完整性。構(gòu)造因式:qi=s+txi+yi,其中,xi、yi分別表示第i個目標(biāo)的在x軸、y軸上的坐標(biāo)值;i=1,…,N。于是有多項式

    (11)

    新量測集即為多項式系數(shù)bij所組成的集合,當(dāng)N=2時,新量測可以表示為

    (12)

    由于觀測量中包含Kamen乘積和SME,因此,當(dāng)新量測集已知時,容易得到x軸和y軸的觀測集分別為{x1,x2}和{y1,y2},并且新量測通過最后一項加強了不同坐標(biāo)之間的關(guān)系。假設(shè)不同坐標(biāo)之間的關(guān)系并不明確,{(x1,y2), (x2,y1)}的組合的形式仍然不能被消除,則有

    (x1y1+x2y2)-(x1y2+x2y1)=0

    (13)

    式(13)表明,當(dāng)x1=x2或y1=y2時,這種坐標(biāo)間關(guān)系不明確的情況才會發(fā)生。然而,一旦x1=x2或y1=y2,集合{(x1,y1), (x2,y2)}與{(x1,y2), (x2,y1)}并沒有區(qū)別。因此,改進(jìn)算法能夠有效抑制“鬼點”現(xiàn)象,提高新量測集的信息完整性。

    2.2觀測誤差協(xié)方差陣

    由于觀測方程發(fā)生改變,與之相應(yīng)的觀測誤差協(xié)方差陣也要做出調(diào)整。但是經(jīng)過對稱變換的觀測方程其觀測誤差往往不再服從高斯分布,在計算時卻通常視作高斯分布近似處理[9],難點是新量測誤差協(xié)方差陣的求解十分困難。

    假設(shè)在一維空間含誤差的原始量測值zi=xi+ui,其中,i=1,…,N;xi為真實值,ui為觀測誤差。則經(jīng)過hi對稱變換后的新量測誤差vi由文獻(xiàn)[3]式(28)泰勒公式展開

    vi=

    (14)

    設(shè)新量測的觀測誤差向量V=[v1,…,vn]T,其中,n為新量測方程個數(shù),通過R=E[VVT]可求得新觀測誤差協(xié)方差陣。根據(jù)式(14)可知,V為多項式向量,則VVT為多項式矩陣,文獻(xiàn)[10]定理1給出了當(dāng)變量服從高斯分布時,求解多項式期望的方法。

    定理 1假設(shè)u=[u1,…,um]T~N(μ,Σ),Σ為m×m的半正定矩陣,對于非負(fù)整數(shù)s1,…,sm,有

    (15)

    對于如式(12)所示的新觀測方程,設(shè)初始觀測誤差向量u=[ux1, ux2, uy1, uy2]T為零均值高斯白噪聲,且滿足獨立同分布,即

    (16)

    式中,σ為初始觀測噪聲標(biāo)準(zhǔn)差。

    則新的觀測誤差協(xié)方差陣有

    R=E[VVT]=

    (17)

    式中,H(x1,x2,y1,y2)是h(x1,x2,y1,y2)的Jacobian陣。

    3仿真分析

    考慮二維空間上兩個勻速直線運動的目標(biāo),設(shè)定初始位置和速度以保證兩目標(biāo)能夠發(fā)生航跡交叉。另外,在真實位置的基礎(chǔ)上附加獨立同分布、零均值高斯白噪聲作為原始測量數(shù)據(jù)。目標(biāo)1初始狀態(tài)為:X1=[500,20,0,30]T,目標(biāo)2初始狀態(tài)為:X2=[0,40,500,10]T,當(dāng)采樣周期T=1s,觀測噪聲標(biāo)準(zhǔn)差σ分別為50m、100m時,原始觀測值和真實值如圖1所示。

    圖1 目標(biāo)觀測值與真實值

    假設(shè)跟蹤時無差別地獲取兩個目標(biāo)的量測值,且量測與目標(biāo)的對應(yīng)關(guān)系未知,不考慮虛警和漏警的情況。經(jīng)過改進(jìn)SME算法濾波后的效果圖如圖2所示。

    圖2 跟蹤效果圖

    由圖2中可看出:改進(jìn)算法在不同觀測噪聲下都能基本實現(xiàn)對兩個目標(biāo)的跟蹤,當(dāng)觀測噪聲比較小時,跟蹤的比較穩(wěn)定,但隨著噪聲的增大,跟蹤的精度也逐漸降低。

    為了評估改進(jìn)算法的跟蹤效果,將改進(jìn)SME算法與Kamen乘積和SME[3]、復(fù)數(shù)乘積和SME[4]以及JPDA算法進(jìn)行比較。分別在觀測噪聲標(biāo)準(zhǔn)差σ=50m和σ=100m條件下,進(jìn)行100次蒙特卡羅仿真,以目標(biāo)1為例,各算法的位置均方根誤差如圖3所示。

    圖3 各算法對目標(biāo)1跟蹤誤差圖

    由圖3中可看出:當(dāng)觀測誤差較小時,幾種算法的跟蹤效果相差不大,僅在航跡交叉點附近有不同程度的差別;當(dāng)觀測誤差較大時,改進(jìn)SME算法跟蹤精度明顯比Kamen乘積和SME、復(fù)數(shù)乘積和SME算法的跟蹤精度要高,基本達(dá)到與JPDA相近的水平。當(dāng)然,評判一個算法的優(yōu)劣不能僅僅考慮跟蹤精度,其時效性也是極為重要的指標(biāo)。分別取不同次數(shù)的蒙特卡羅仿真,統(tǒng)計各復(fù)數(shù)運算SME和JPDA算法的總運行時間如表1所示。

    表1 各算法總運行時間比較 s

    由表1可以看出:隨著仿真次數(shù)的增加,各算法的運行時間逐漸增大,同一仿真次數(shù)下SME算法運行時間相差不多,改進(jìn)SME算法比基于復(fù)數(shù)運算SME算法的時間性能還要好一點,而JPDA算法的運行時間要比其他3種算法的運行時間大很多,這是由于JPDA涉及數(shù)據(jù)關(guān)聯(lián),而其他3種算法不需要數(shù)據(jù)關(guān)聯(lián)運算,因此基于SME算法在時效性方面比JPDA算法更有優(yōu)勢。

    4結(jié)論

    本文介紹了一種基于SME的多目標(biāo)跟蹤算法,該算法不需要復(fù)雜的數(shù)據(jù)關(guān)聯(lián)運算就能同時對多個目標(biāo)進(jìn)行跟蹤。針對現(xiàn)有SME的不足之處,構(gòu)造了一種基于多項式因式分解思想的改進(jìn)SME,在觀測誤差服從獨立同分布、零均值高斯白噪聲前提下,推導(dǎo)出相應(yīng)的觀測誤差協(xié)方差陣。

    在不同觀測噪聲下對存在航跡交叉的兩個目標(biāo)進(jìn)行跟蹤,取得了較好的跟蹤效果,驗證了新算法的有效性;與現(xiàn)有SME算法和JPDA算法對比,新算法不僅延續(xù)了SME時效性優(yōu)勢,而且跟蹤精度基本達(dá)到與JPDA相近的水平,有效抑制了航跡交叉時跟蹤誤差較大現(xiàn)象。然而新算法也存在一些不足,比如:構(gòu)造的新觀測量數(shù)目比現(xiàn)有SME要多一些,而且利用改進(jìn)算法進(jìn)行多目標(biāo)跟蹤時,目標(biāo)數(shù)目必須已知。

    目前SME算法發(fā)展并不成熟,許多問題還沒有很好的解決,比如:存在虛警漏警時計算較為復(fù)雜、目標(biāo)數(shù)目未知情況下應(yīng)用受限等。因此,需要投入更多的精力去了解、研究SME算法,使之成為一種成熟的理論,以便更好地解決多目標(biāo)跟蹤問題。

    參考文獻(xiàn):

    [1] Shalom Y B, Daum F, Huang J. The probabilistic data association filter[J].IEEETrans.onControlSystemsMagazine, 2009,29(6): 82-100.

    [2] Kamen E W. Multiple target tracking based on symmetric mea-surement equations[J].IEEETrans.onAutomaticControl, 1992, 37(3): 371-374.

    [3] Kamen E W, Sastry C R. Multiple target tracking using pro-ducts of position measurements[J].IEEETrans.onAerospaceandElectronicSystems, 1993, 29(2):476-493.

    [4] Lee Y. The SME filter approach to multiple target tracking with false and missing measurements[D]. Georgia: School of Electrical and Computer Engineering, 1994.

    [5] Kamen E W, Su J K.Introductiontooptimalestimation[M]. London: Springer-Verlag, 1999.

    [6] Muder D J, O’Neil S D. The multi-dimensional SME filter for multi-target tracking[C]∥Proc.oftheSPIESignalandDataProcessingofSmallTargets, 1993: 587-599.

    [7] Leven W F, Lanterman A D. Multiple target tracking with symmetric measurement equations using unscented Kalman and particle filters[C]∥Proc.ofthe36thSoutheasternSymposiumonSystemTheory, 2004:195-199.

    [8] Baum M, Noack B, Beutler F, et al. Optimal Gaussian filtering for polynomial systems applied to association-free multi-target tracking[C]∥Proc.ofthe14thInternationalConferenceonInformationFusion, 2011:1-8.

    [9] Leven W F. Approximate Cramer-Rao bounds for multiple target tracking[D]. Georgia: School of Electrical and Computer Engineering, 2006.

    [10] Raymond K. From moments of sum to moments of product[J].JournalofMultivariateAnalysis, 2008, 99(3): 542-554.

    吳瀟(1991-),男,碩士研究生,主要研究方向為多目標(biāo)跟蹤。

    E-mail:february27th@163.com

    黃樹彩(1967-),男,教授,博士,主要研究方向為系統(tǒng)辨識與模式識別、目標(biāo)檢測與跟蹤。

    E-mail:hsc67118@126.com

    凌強(1990-),男,碩士研究生,主要研究方向為紅外弱小目標(biāo)檢測與跟蹤。

    E-mail:lq910131@gmail.com

    鐘宇(1987-),男,博士研究生,主要研究方向為系統(tǒng)建模與仿真、天基紅外預(yù)警探測與跟蹤。

    E-mail:zhongyu257678@163.com

    網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20150818.1519.012.html

    Multiple targets tracking using modified symmetric measurement equations

    WU Xiao, HUANG Shu-cai, LING Qiang, ZHONG Yu

    (AirandMissileDefenseCollege,AirForceEngineeringUniversity,Xi’an710051,China)

    Abstract:The classical symmetric measurement equation (SME) approach will produce pseudo-measurements while tracking multiple targets in multi-dimensions, which induces large tracking errors around the coordinate switching areas. To solve this problem, a modified SME is proposed by polynomial factoring to connect states of one target with different coordinates. As a result, the mapping relationships of original observations to new mesurements are described more precisely. The new measurement noise covariance matrix is deduced corresponding to the new SME measure noises which are derived by Taylor series expansions. Finally, simulations are constructed to validate the new SME’s characteristics by comparing with the existing SME approaches and joint probability data association algorithm respectively. Results show that the modified SME approach not only maintains the advantage of real-time performances, but also does well in improving tracking accuracy during the coordinate switching areas.

    Keywords:symmetric measurement equation (SME); multiple targets tracking; polynomial factoring

    作者簡介:

    中圖分類號:V 249

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

    DOI:10.3969/j.issn.1001-506X.2016.01.04

    基金項目:航空科學(xué)基金(20130196004)資助課題

    收稿日期:2015-01-13;修回日期:2015-06-29;網(wǎng)絡(luò)優(yōu)先出版日期:2015-08-18。

    黑人高潮一二区| 成人美女网站在线观看视频| 97超级碰碰碰精品色视频在线观看| 简卡轻食公司| 能在线免费观看的黄片| 亚洲成a人片在线一区二区| 又黄又爽又刺激的免费视频.| 在线国产一区二区在线| 午夜福利在线观看免费完整高清在 | 国产亚洲91精品色在线| 国产真实伦视频高清在线观看| 久久婷婷人人爽人人干人人爱| 男女啪啪激烈高潮av片| av在线观看视频网站免费| 香蕉av资源在线| 卡戴珊不雅视频在线播放| 女人被狂操c到高潮| 在线播放无遮挡| 国产高潮美女av| 97超视频在线观看视频| 久久久国产成人免费| 真实男女啪啪啪动态图| 精品一区二区三区视频在线观看免费| 国产av不卡久久| 禁无遮挡网站| 久久久a久久爽久久v久久| 人人妻人人澡人人爽人人夜夜 | 精品人妻熟女av久视频| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 女生性感内裤真人,穿戴方法视频| 在现免费观看毛片| 欧美最黄视频在线播放免费| 日韩av不卡免费在线播放| 亚洲内射少妇av| 久久久色成人| 亚洲人成网站在线观看播放| 亚洲激情五月婷婷啪啪| 亚洲av美国av| 日韩欧美免费精品| 中出人妻视频一区二区| 国产精品国产高清国产av| 国产精品无大码| 听说在线观看完整版免费高清| 国产成人aa在线观看| 啦啦啦观看免费观看视频高清| 欧美最新免费一区二区三区| 一进一出好大好爽视频| 国产一区二区三区在线臀色熟女| 国产成人aa在线观看| a级一级毛片免费在线观看| 精品一区二区三区人妻视频| 精品久久国产蜜桃| 99热这里只有是精品在线观看| 赤兔流量卡办理| 又粗又爽又猛毛片免费看| 男女边吃奶边做爰视频| 久久精品夜色国产| 成人永久免费在线观看视频| 18禁黄网站禁片免费观看直播| 男女之事视频高清在线观看| 精品日产1卡2卡| 久久鲁丝午夜福利片| 97人妻精品一区二区三区麻豆| 午夜亚洲福利在线播放| 一卡2卡三卡四卡精品乱码亚洲| 搞女人的毛片| 一本一本综合久久| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 久久精品影院6| 麻豆成人午夜福利视频| 亚洲无线观看免费| 国产不卡一卡二| 久久韩国三级中文字幕| 久久久久久久亚洲中文字幕| 男人舔奶头视频| 亚洲最大成人av| 国产视频内射| 久久鲁丝午夜福利片| 久久精品久久久久久噜噜老黄 | 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器| 亚洲欧美成人精品一区二区| 一级a爱片免费观看的视频| 一区二区三区四区激情视频 | 欧美最新免费一区二区三区| 国产高清视频在线观看网站| 乱系列少妇在线播放| 赤兔流量卡办理| 国产av在哪里看| 综合色丁香网| 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| 91久久精品国产一区二区三区| 国产日本99.免费观看| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 99久久无色码亚洲精品果冻| 乱码一卡2卡4卡精品| 精品熟女少妇av免费看| or卡值多少钱| 国产91av在线免费观看| 精品熟女少妇av免费看| 成人漫画全彩无遮挡| 欧美成人精品欧美一级黄| 伦理电影大哥的女人| 99久久精品热视频| 国产人妻一区二区三区在| 真人做人爱边吃奶动态| 亚洲成人精品中文字幕电影| 欧美成人a在线观看| 亚洲天堂国产精品一区在线| 亚洲精品亚洲一区二区| а√天堂www在线а√下载| 精品日产1卡2卡| 欧美日韩在线观看h| 国产一区二区亚洲精品在线观看| 啦啦啦啦在线视频资源| 欧美最新免费一区二区三区| 亚洲18禁久久av| 亚洲欧美中文字幕日韩二区| 亚洲最大成人中文| 久久中文看片网| 国产人妻一区二区三区在| 亚洲成av人片在线播放无| 少妇人妻一区二区三区视频| 俺也久久电影网| 欧美最黄视频在线播放免费| 3wmmmm亚洲av在线观看| 亚洲精品久久国产高清桃花| 午夜日韩欧美国产| 久久这里只有精品中国| 日本一二三区视频观看| av在线亚洲专区| 最新中文字幕久久久久| 1024手机看黄色片| 不卡视频在线观看欧美| 亚洲18禁久久av| 成熟少妇高潮喷水视频| 欧美不卡视频在线免费观看| 亚洲熟妇熟女久久| 久久草成人影院| 网址你懂的国产日韩在线| 日韩三级伦理在线观看| 久久久国产成人精品二区| 国产精品一区二区三区四区免费观看 | 国产片特级美女逼逼视频| 久久久欧美国产精品| 97人妻精品一区二区三区麻豆| 亚洲精品影视一区二区三区av| 亚洲在线自拍视频| 五月伊人婷婷丁香| 欧美色视频一区免费| 中文字幕av在线有码专区| 国产午夜精品久久久久久一区二区三区 | 国产麻豆成人av免费视频| 91狼人影院| 婷婷六月久久综合丁香| 非洲黑人性xxxx精品又粗又长| 人妻制服诱惑在线中文字幕| 国产av不卡久久| 97超级碰碰碰精品色视频在线观看| 看十八女毛片水多多多| 一级毛片电影观看 | 最近中文字幕高清免费大全6| 精品日产1卡2卡| 日本-黄色视频高清免费观看| 淫秽高清视频在线观看| 亚洲国产精品久久男人天堂| 国产高潮美女av| 日韩成人av中文字幕在线观看 | 亚洲精品色激情综合| а√天堂www在线а√下载| 舔av片在线| 亚洲精品一区av在线观看| 国内揄拍国产精品人妻在线| 精品午夜福利在线看| 中国美女看黄片| 国产精品久久久久久久久免| 国产高潮美女av| 天美传媒精品一区二区| 亚洲丝袜综合中文字幕| eeuss影院久久| 日韩亚洲欧美综合| 亚洲av二区三区四区| 国产黄a三级三级三级人| 女人被狂操c到高潮| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 日韩欧美三级三区| 欧美另类亚洲清纯唯美| av在线观看视频网站免费| 在线观看午夜福利视频| 精品少妇黑人巨大在线播放 | 亚洲中文字幕一区二区三区有码在线看| 听说在线观看完整版免费高清| 亚洲第一电影网av| 国产男靠女视频免费网站| 精品久久久久久久末码| 欧美激情国产日韩精品一区| 国国产精品蜜臀av免费| 最近最新中文字幕大全电影3| 国产av在哪里看| 一级毛片我不卡| 国产精品人妻久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 我要搜黄色片| 啦啦啦韩国在线观看视频| 亚洲精品成人久久久久久| 日韩亚洲欧美综合| 黄色配什么色好看| 有码 亚洲区| 一卡2卡三卡四卡精品乱码亚洲| 国产高潮美女av| 欧美在线一区亚洲| 日韩在线高清观看一区二区三区| 精品午夜福利在线看| 免费看日本二区| 身体一侧抽搐| 岛国在线免费视频观看| 少妇的逼水好多| 日本 av在线| 亚洲人与动物交配视频| 日本黄色视频三级网站网址| 国产成人a∨麻豆精品| 中文字幕av在线有码专区| 国产成人91sexporn| 少妇熟女欧美另类| 色哟哟哟哟哟哟| 岛国在线免费视频观看| 亚洲av中文字字幕乱码综合| 美女黄网站色视频| 久久草成人影院| 国产高清视频在线播放一区| 国产成人a∨麻豆精品| 午夜日韩欧美国产| 在线免费观看的www视频| 欧美丝袜亚洲另类| 三级国产精品欧美在线观看| 性色avwww在线观看| 伦精品一区二区三区| 亚洲自偷自拍三级| 精品久久久久久成人av| 国产精品无大码| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 日本a在线网址| 亚洲丝袜综合中文字幕| 五月玫瑰六月丁香| 一本精品99久久精品77| 搡老岳熟女国产| 成人av在线播放网站| 黄色日韩在线| 欧美一区二区亚洲| 国产在线男女| 国产av不卡久久| 别揉我奶头 嗯啊视频| 亚洲精品在线观看二区| 如何舔出高潮| 亚洲一区高清亚洲精品| 欧美国产日韩亚洲一区| 欧美日韩综合久久久久久| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 一个人看的www免费观看视频| 91麻豆精品激情在线观看国产| 亚洲欧美日韩高清专用| 又黄又爽又刺激的免费视频.| 国产综合懂色| 22中文网久久字幕| 一个人观看的视频www高清免费观看| 成人综合一区亚洲| 人人妻,人人澡人人爽秒播| 日本免费a在线| 久久久久精品国产欧美久久久| 大型黄色视频在线免费观看| 91在线精品国自产拍蜜月| 国产成人a区在线观看| av在线观看视频网站免费| 嫩草影院新地址| 日本一本二区三区精品| 国产精品久久久久久久电影| 久久精品夜夜夜夜夜久久蜜豆| 激情 狠狠 欧美| 美女高潮的动态| 女同久久另类99精品国产91| 免费av观看视频| 97碰自拍视频| 一级毛片久久久久久久久女| 亚洲欧美精品综合久久99| 在线a可以看的网站| 免费高清视频大片| 中文字幕精品亚洲无线码一区| a级一级毛片免费在线观看| 日韩一本色道免费dvd| 亚洲国产欧美人成| 夜夜夜夜夜久久久久| 亚洲av五月六月丁香网| 波多野结衣高清作品| 日本撒尿小便嘘嘘汇集6| 大香蕉久久网| 亚洲av.av天堂| 性插视频无遮挡在线免费观看| 别揉我奶头 嗯啊视频| 国产免费一级a男人的天堂| 久久久精品94久久精品| 熟女人妻精品中文字幕| 免费一级毛片在线播放高清视频| 欧美一区二区精品小视频在线| 久久6这里有精品| 真实男女啪啪啪动态图| 精品人妻偷拍中文字幕| 男女之事视频高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 99久久中文字幕三级久久日本| 51国产日韩欧美| 精品少妇黑人巨大在线播放 | 男人狂女人下面高潮的视频| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 色视频www国产| 日韩精品中文字幕看吧| 精品午夜福利在线看| 韩国高清视频一区二区三区| 插逼视频在线观看| 亚洲情色 制服丝袜| videos熟女内射| 黄片无遮挡物在线观看| 久久午夜福利片| 亚洲人成网站在线观看播放| 少妇人妻 视频| 国国产精品蜜臀av免费| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 51国产日韩欧美| 免费看光身美女| 少妇被粗大猛烈的视频| 亚洲四区av| 亚洲国产精品999| 国产成人免费观看mmmm| 午夜激情福利司机影院| 国产成人freesex在线| 七月丁香在线播放| 51国产日韩欧美| 免费看光身美女| 精品一区二区免费观看| 一本大道久久a久久精品| 亚洲精华国产精华液的使用体验| 日本午夜av视频| 国模一区二区三区四区视频| 精品一品国产午夜福利视频| 亚洲欧美精品专区久久| 男男h啪啪无遮挡| 中文字幕久久专区| 精华霜和精华液先用哪个| 亚洲综合精品二区| 欧美一级a爱片免费观看看| 国产精品成人在线| 我要看日韩黄色一级片| 欧美 日韩 精品 国产| 久久久久久久久久人人人人人人| 亚洲精品乱久久久久久| 狂野欧美白嫩少妇大欣赏| 女性生殖器流出的白浆| 精品久久久久久久久av| 亚洲国产成人一精品久久久| 女性生殖器流出的白浆| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av成人精品| 国产伦精品一区二区三区四那| 国产精品成人在线| 国产免费视频播放在线视频| 又黄又爽又刺激的免费视频.| 中文欧美无线码| 久久 成人 亚洲| 一级毛片aaaaaa免费看小| 亚洲精品日韩在线中文字幕| 日韩成人av中文字幕在线观看| 久久精品国产亚洲av天美| 大话2 男鬼变身卡| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品久久久久久婷婷小说| 最近最新中文字幕免费大全7| av女优亚洲男人天堂| 亚洲无线观看免费| 亚洲一级一片aⅴ在线观看| 精品人妻熟女av久视频| 国产在线视频一区二区| 亚洲伊人久久精品综合| 在线 av 中文字幕| 亚洲欧美成人精品一区二区| 69精品国产乱码久久久| 免费观看无遮挡的男女| 蜜桃久久精品国产亚洲av| av.在线天堂| 寂寞人妻少妇视频99o| 黑人猛操日本美女一级片| 欧美xxxx性猛交bbbb| 欧美+日韩+精品| xxx大片免费视频| 亚洲av不卡在线观看| 在线观看av片永久免费下载| 欧美日韩国产mv在线观看视频| 欧美老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 欧美日韩av久久| 久久女婷五月综合色啪小说| √禁漫天堂资源中文www| 中文字幕制服av| 一级,二级,三级黄色视频| 性色avwww在线观看| 不卡视频在线观看欧美| 国产亚洲午夜精品一区二区久久| 免费黄网站久久成人精品| 一级,二级,三级黄色视频| 亚洲国产最新在线播放| 色5月婷婷丁香| 色吧在线观看| 热re99久久精品国产66热6| 亚洲国产精品一区二区三区在线| 日本av手机在线免费观看| 黄色欧美视频在线观看| 精品国产乱码久久久久久小说| 水蜜桃什么品种好| 国产深夜福利视频在线观看| 夫妻性生交免费视频一级片| a级毛片免费高清观看在线播放| 高清在线视频一区二区三区| 深夜a级毛片| 91精品伊人久久大香线蕉| 亚洲精品一二三| 久久久久久伊人网av| 国产美女午夜福利| 九九在线视频观看精品| 一级二级三级毛片免费看| 亚洲国产成人一精品久久久| 国产成人精品久久久久久| 免费av中文字幕在线| 人妻少妇偷人精品九色| 久久久国产精品麻豆| 国产一区亚洲一区在线观看| 麻豆乱淫一区二区| .国产精品久久| 各种免费的搞黄视频| 日韩一区二区三区影片| 99精国产麻豆久久婷婷| 成人国产麻豆网| 男女边吃奶边做爰视频| 亚洲美女搞黄在线观看| 欧美精品一区二区免费开放| a 毛片基地| 成年美女黄网站色视频大全免费 | 黄色毛片三级朝国网站 | 激情五月婷婷亚洲| 一个人免费看片子| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 男女边摸边吃奶| 国产亚洲精品久久久com| 人人妻人人澡人人爽人人夜夜| 欧美另类一区| 中文字幕制服av| 午夜福利在线观看免费完整高清在| 美女福利国产在线| 久久av网站| 啦啦啦在线观看免费高清www| 亚洲国产欧美日韩在线播放 | 老司机影院成人| av专区在线播放| 久久久久视频综合| 亚洲精品一区蜜桃| 日韩精品有码人妻一区| 美女视频免费永久观看网站| 一区二区三区精品91| 国产老妇伦熟女老妇高清| 熟女电影av网| 91在线精品国自产拍蜜月| 99国产精品免费福利视频| 18禁在线播放成人免费| 80岁老熟妇乱子伦牲交| 色哟哟·www| 成人综合一区亚洲| 国产一区二区在线观看av| 另类精品久久| 天天躁夜夜躁狠狠久久av| av不卡在线播放| 久久久久国产网址| 免费大片18禁| 欧美老熟妇乱子伦牲交| 亚洲内射少妇av| 成人综合一区亚洲| 七月丁香在线播放| 久久精品国产亚洲av天美| 91aial.com中文字幕在线观看| 一级黄片播放器| 日韩av不卡免费在线播放| 国产亚洲精品久久久com| 日本黄大片高清| 2021少妇久久久久久久久久久| 黄色日韩在线| 亚洲婷婷狠狠爱综合网| 国产精品麻豆人妻色哟哟久久| 丰满人妻一区二区三区视频av| 熟妇人妻不卡中文字幕| 久久精品国产自在天天线| 妹子高潮喷水视频| 91aial.com中文字幕在线观看| 久久99热6这里只有精品| 蜜臀久久99精品久久宅男| 国产中年淑女户外野战色| 久久免费观看电影| 国产欧美日韩一区二区三区在线 | 观看美女的网站| 一级毛片黄色毛片免费观看视频| 国产视频内射| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 国产精品久久久久久av不卡| a级毛片在线看网站| 99热这里只有精品一区| 一级a做视频免费观看| 看免费成人av毛片| 简卡轻食公司| 国内揄拍国产精品人妻在线| 精品视频人人做人人爽| 91久久精品国产一区二区成人| 成人毛片a级毛片在线播放| 欧美成人精品欧美一级黄| 自拍欧美九色日韩亚洲蝌蚪91 | 丰满少妇做爰视频| 男人舔奶头视频| 乱人伦中国视频| 丰满乱子伦码专区| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 一区二区av电影网| 另类亚洲欧美激情| 亚洲,一卡二卡三卡| 日韩欧美 国产精品| 中文字幕制服av| 一级a做视频免费观看| 亚洲中文av在线| 国产 精品1| 亚洲欧洲国产日韩| 欧美高清成人免费视频www| 亚洲高清免费不卡视频| 汤姆久久久久久久影院中文字幕| 中文在线观看免费www的网站| 国产高清国产精品国产三级| 国产av精品麻豆| 午夜激情福利司机影院| 久久久久久伊人网av| 久久av网站| 亚洲av免费高清在线观看| 哪个播放器可以免费观看大片| 青春草视频在线免费观看| 免费看日本二区| 亚洲精品乱久久久久久| 亚洲欧美日韩东京热| av女优亚洲男人天堂| 久久精品久久久久久噜噜老黄| 精品视频人人做人人爽| 精品少妇内射三级| 久久免费观看电影| 新久久久久国产一级毛片| 免费大片18禁| 在线看a的网站| 国产爽快片一区二区三区| 精品一区在线观看国产| 大陆偷拍与自拍| 啦啦啦啦在线视频资源| 啦啦啦在线观看免费高清www| www.av在线官网国产| 久久久久国产精品人妻一区二区| 亚洲综合色惰| 中文资源天堂在线| 在现免费观看毛片| 精品国产乱码久久久久久小说| 赤兔流量卡办理| 春色校园在线视频观看| 国产熟女午夜一区二区三区 | 男女边摸边吃奶| av在线app专区| 综合色丁香网| 久久久久久久久久久久大奶| 精品人妻一区二区三区麻豆| 综合色丁香网| 22中文网久久字幕| 91成人精品电影| 嫩草影院新地址| 日韩av不卡免费在线播放| 国产亚洲av片在线观看秒播厂| 黄色怎么调成土黄色| 在线观看av片永久免费下载| 涩涩av久久男人的天堂| 人人妻人人添人人爽欧美一区卜| 美女大奶头黄色视频| 中文字幕人妻熟人妻熟丝袜美| 午夜免费鲁丝| 国产永久视频网站| 又爽又黄a免费视频| 少妇裸体淫交视频免费看高清| 成人毛片a级毛片在线播放| 伊人久久国产一区二区| 视频中文字幕在线观看| 日韩人妻高清精品专区| 国产亚洲av片在线观看秒播厂| 国产av国产精品国产| 国产黄色免费在线视频| 在线观看一区二区三区激情| 亚洲成人手机| 少妇人妻一区二区三区视频|