• <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。

    免费av毛片视频| 午夜激情av网站| 黄色a级毛片大全视频| 无人区码免费观看不卡| а√天堂www在线а√下载| 免费高清视频大片| 国产精品久久久av美女十八| 欧美成人午夜精品| 又黄又粗又硬又大视频| 啦啦啦韩国在线观看视频| 亚洲欧美精品综合久久99| 啪啪无遮挡十八禁网站| 天堂√8在线中文| av在线播放免费不卡| 纯流量卡能插随身wifi吗| 九色亚洲精品在线播放| 女人爽到高潮嗷嗷叫在线视频| 日韩大尺度精品在线看网址 | 成人三级做爰电影| 中文字幕人妻熟女乱码| 午夜福利视频1000在线观看 | 亚洲在线自拍视频| 久久久久国产一级毛片高清牌| 国产精品久久久久久人妻精品电影| 久久久水蜜桃国产精品网| 宅男免费午夜| 男女午夜视频在线观看| 757午夜福利合集在线观看| 精品国内亚洲2022精品成人| 国产亚洲欧美精品永久| 可以在线观看的亚洲视频| 69av精品久久久久久| 在线永久观看黄色视频| 久久人人精品亚洲av| 国产亚洲精品久久久久久毛片| 美国免费a级毛片| 一级毛片精品| 九色国产91popny在线| 成人亚洲精品av一区二区| 国产国语露脸激情在线看| 中亚洲国语对白在线视频| 成人亚洲精品一区在线观看| 久久久久久人人人人人| 少妇熟女aⅴ在线视频| 看免费av毛片| 在线播放国产精品三级| 99久久精品国产亚洲精品| 久久天躁狠狠躁夜夜2o2o| 免费高清在线观看日韩| 最好的美女福利视频网| 久久国产精品影院| 国产成+人综合+亚洲专区| 精品国内亚洲2022精品成人| 黄色丝袜av网址大全| 欧美黄色淫秽网站| 欧美人与性动交α欧美精品济南到| 亚洲黑人精品在线| 两性夫妻黄色片| 黄色视频不卡| 国产精品,欧美在线| 国产不卡一卡二| 国产野战对白在线观看| 欧美+亚洲+日韩+国产| 欧美一区二区精品小视频在线| 午夜视频精品福利| 欧美乱色亚洲激情| 人人澡人人妻人| 欧美大码av| 久久人人爽av亚洲精品天堂| 老司机福利观看| 婷婷丁香在线五月| 黄色片一级片一级黄色片| 国产亚洲av嫩草精品影院| 日韩大码丰满熟妇| 国产精品久久视频播放| 99国产精品一区二区蜜桃av| 久久精品国产亚洲av高清一级| 午夜亚洲福利在线播放| 精品久久久久久久毛片微露脸| 国产精品日韩av在线免费观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲五月婷婷丁香| 久久午夜亚洲精品久久| 后天国语完整版免费观看| 人成视频在线观看免费观看| 国内毛片毛片毛片毛片毛片| 极品教师在线免费播放| 99国产综合亚洲精品| 国产精品av久久久久免费| 一个人免费在线观看的高清视频| 99在线人妻在线中文字幕| 日本精品一区二区三区蜜桃| 亚洲视频免费观看视频| 国产精品久久久久久精品电影 | 色在线成人网| 亚洲精品中文字幕在线视频| 久久久久国产精品人妻aⅴ院| 精品国产国语对白av| 亚洲中文日韩欧美视频| 精品一区二区三区av网在线观看| 黑人欧美特级aaaaaa片| 精品一区二区三区av网在线观看| www.精华液| 熟女少妇亚洲综合色aaa.| 欧美色欧美亚洲另类二区 | 精品福利观看| 国产亚洲精品综合一区在线观看 | 国产av一区二区精品久久| 美女大奶头视频| 久久久久久人人人人人| 亚洲精品国产区一区二| 久久精品91蜜桃| 亚洲一区二区三区色噜噜| 亚洲欧美日韩另类电影网站| 免费无遮挡裸体视频| 女性被躁到高潮视频| 中文字幕精品免费在线观看视频| 色播在线永久视频| 波多野结衣巨乳人妻| 久久 成人 亚洲| 亚洲精品国产一区二区精华液| 十分钟在线观看高清视频www| 精品无人区乱码1区二区| x7x7x7水蜜桃| 黄色片一级片一级黄色片| 亚洲人成伊人成综合网2020| 久久久精品国产亚洲av高清涩受| 18禁黄网站禁片午夜丰满| 欧美乱妇无乱码| 制服诱惑二区| 麻豆国产av国片精品| 午夜福利视频1000在线观看 | 真人做人爱边吃奶动态| 亚洲国产日韩欧美精品在线观看 | 久久久久久久午夜电影| 美女高潮到喷水免费观看| 欧美激情极品国产一区二区三区| 久久香蕉国产精品| 久久久久久久午夜电影| 国产精品日韩av在线免费观看 | 色尼玛亚洲综合影院| 日韩欧美在线二视频| 欧美另类亚洲清纯唯美| 99久久精品国产亚洲精品| 热99re8久久精品国产| 老司机午夜十八禁免费视频| 婷婷丁香在线五月| 免费在线观看完整版高清| 亚洲av美国av| 视频区欧美日本亚洲| 久久久久久大精品| 亚洲人成77777在线视频| 亚洲性夜色夜夜综合| 狠狠狠狠99中文字幕| 在线观看午夜福利视频| 18美女黄网站色大片免费观看| 午夜精品在线福利| 亚洲人成电影免费在线| 很黄的视频免费| www国产在线视频色| 国产精品久久电影中文字幕| 丰满人妻熟妇乱又伦精品不卡| 欧美久久黑人一区二区| 精品国产乱子伦一区二区三区| 国内精品久久久久精免费| 91精品国产国语对白视频| 久久 成人 亚洲| 国产亚洲精品av在线| av电影中文网址| 超碰成人久久| 757午夜福利合集在线观看| 日韩国内少妇激情av| 黄色成人免费大全| 一本综合久久免费| 亚洲一区二区三区不卡视频| netflix在线观看网站| 日韩国内少妇激情av| 又紧又爽又黄一区二区| 国产精品亚洲一级av第二区| 妹子高潮喷水视频| 久久婷婷成人综合色麻豆| 大型黄色视频在线免费观看| 日韩精品青青久久久久久| 国产男靠女视频免费网站| 美女午夜性视频免费| avwww免费| 韩国精品一区二区三区| 成人特级黄色片久久久久久久| 免费一级毛片在线播放高清视频 | ponron亚洲| 午夜影院日韩av| 一a级毛片在线观看| 欧美午夜高清在线| 欧美黑人欧美精品刺激| 一个人免费在线观看的高清视频| 久久人妻熟女aⅴ| 久久久久久久久久久久大奶| 久久九九热精品免费| 中文字幕人妻熟女乱码| 国产熟女午夜一区二区三区| 久久九九热精品免费| 亚洲欧美激情综合另类| 亚洲av电影不卡..在线观看| 国产aⅴ精品一区二区三区波| 国产99久久九九免费精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲一区中文字幕在线| 动漫黄色视频在线观看| 伊人久久大香线蕉亚洲五| 日韩成人在线观看一区二区三区| 国产主播在线观看一区二区| 免费在线观看完整版高清| 90打野战视频偷拍视频| 久久人妻熟女aⅴ| 18禁裸乳无遮挡免费网站照片 | 国产精品香港三级国产av潘金莲| 国产成年人精品一区二区| 国产欧美日韩精品亚洲av| 视频区欧美日本亚洲| 亚洲欧美一区二区三区黑人| 视频在线观看一区二区三区| 99香蕉大伊视频| 国内精品久久久久久久电影| 成人亚洲精品av一区二区| av天堂久久9| 欧美中文日本在线观看视频| 成人亚洲精品av一区二区| svipshipincom国产片| 淫妇啪啪啪对白视频| 成人精品一区二区免费| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区mp4| 国产高清视频在线播放一区| 成人三级黄色视频| 国产欧美日韩一区二区三| 久久天堂一区二区三区四区| 久久精品人人爽人人爽视色| 久久亚洲真实| 女警被强在线播放| 成人三级做爰电影| 级片在线观看| 中文字幕久久专区| 国产熟女xx| 国产亚洲精品久久久久久毛片| 伊人久久大香线蕉亚洲五| 婷婷精品国产亚洲av在线| 国产不卡一卡二| 久久香蕉国产精品| 母亲3免费完整高清在线观看| 如日韩欧美国产精品一区二区三区| 深夜精品福利| 99精品欧美一区二区三区四区| 日本a在线网址| 99riav亚洲国产免费| 男人操女人黄网站| 国产精品久久电影中文字幕| 一级毛片高清免费大全| 老司机在亚洲福利影院| 此物有八面人人有两片| 女生性感内裤真人,穿戴方法视频| 精品国产乱码久久久久久男人| 69av精品久久久久久| 波多野结衣一区麻豆| 亚洲欧美精品综合久久99| 人人澡人人妻人| 超碰成人久久| 免费在线观看影片大全网站| 国产精品久久久av美女十八| 亚洲精品粉嫩美女一区| 十八禁人妻一区二区| 97超级碰碰碰精品色视频在线观看| 最近最新中文字幕大全电影3 | 1024视频免费在线观看| 欧美日韩一级在线毛片| www.精华液| 亚洲国产欧美网| 精品国产超薄肉色丝袜足j| 亚洲成人精品中文字幕电影| 变态另类丝袜制服| 亚洲黑人精品在线| 日本撒尿小便嘘嘘汇集6| 最新美女视频免费是黄的| 欧美性长视频在线观看| 午夜福利高清视频| 久久婷婷人人爽人人干人人爱 | 亚洲第一电影网av| 国产激情欧美一区二区| 国产精品国产高清国产av| 久久久久久人人人人人| 亚洲精品久久成人aⅴ小说| 亚洲精品粉嫩美女一区| 国产精品久久视频播放| 狂野欧美激情性xxxx| 国产1区2区3区精品| 国产视频一区二区在线看| 亚洲美女黄片视频| 亚洲中文av在线| 丝袜人妻中文字幕| 大码成人一级视频| 一区二区三区国产精品乱码| 中文亚洲av片在线观看爽| 97人妻精品一区二区三区麻豆 | 国产精品1区2区在线观看.| 亚洲精品在线观看二区| 啦啦啦 在线观看视频| 丝袜人妻中文字幕| 熟女少妇亚洲综合色aaa.| 色综合婷婷激情| 精品一区二区三区视频在线观看免费| 一a级毛片在线观看| 亚洲av日韩精品久久久久久密| av超薄肉色丝袜交足视频| 成人国语在线视频| 亚洲 欧美一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久,| 成人亚洲精品一区在线观看| 两性夫妻黄色片| 午夜久久久在线观看| 欧美绝顶高潮抽搐喷水| 一进一出抽搐gif免费好疼| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 午夜福利成人在线免费观看| 一边摸一边抽搐一进一出视频| 亚洲国产毛片av蜜桃av| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲色图av天堂| 精品久久久久久,| or卡值多少钱| 如日韩欧美国产精品一区二区三区| 天堂影院成人在线观看| 9热在线视频观看99| 欧美日本亚洲视频在线播放| 午夜福利免费观看在线| 亚洲人成电影观看| 色尼玛亚洲综合影院| 岛国在线观看网站| 97人妻精品一区二区三区麻豆 | 少妇熟女aⅴ在线视频| 日韩成人在线观看一区二区三区| 成在线人永久免费视频| 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 国产精品免费视频内射| 亚洲伊人色综图| 亚洲精品美女久久久久99蜜臀| www.www免费av| 久久久久精品国产欧美久久久| 久久久久国产精品人妻aⅴ院| 美女免费视频网站| 国产高清有码在线观看视频 | 级片在线观看| 一边摸一边抽搐一进一出视频| 国产野战对白在线观看| 香蕉丝袜av| 久久草成人影院| 亚洲午夜精品一区,二区,三区| 精品国产一区二区三区四区第35| 最近最新中文字幕大全免费视频| 久久久国产成人精品二区| 亚洲成a人片在线一区二区| 国产精品电影一区二区三区| 老熟妇仑乱视频hdxx| 国产精品一区二区三区四区久久 | 中文字幕人成人乱码亚洲影| 久久人妻熟女aⅴ| 亚洲欧洲精品一区二区精品久久久| 精品不卡国产一区二区三区| 欧美老熟妇乱子伦牲交| 69精品国产乱码久久久| 欧美最黄视频在线播放免费| 国产真人三级小视频在线观看| 亚洲精品中文字幕一二三四区| 久久草成人影院| 曰老女人黄片| av视频免费观看在线观看| 精品欧美国产一区二区三| 国产在线精品亚洲第一网站| 一区二区三区高清视频在线| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 国产1区2区3区精品| 天天躁夜夜躁狠狠躁躁| 色综合站精品国产| 在线观看日韩欧美| 看片在线看免费视频| 亚洲av成人不卡在线观看播放网| 国产成人精品在线电影| av片东京热男人的天堂| 国产免费男女视频| 男女下面进入的视频免费午夜 | 18禁观看日本| 黄色视频不卡| 国产国语露脸激情在线看| netflix在线观看网站| 人人妻,人人澡人人爽秒播| 中文字幕另类日韩欧美亚洲嫩草| 91麻豆av在线| 一级黄色大片毛片| 免费av毛片视频| 成年人黄色毛片网站| 中文字幕av电影在线播放| 神马国产精品三级电影在线观看 | 日韩欧美在线二视频| 久热爱精品视频在线9| 男女床上黄色一级片免费看| 久久性视频一级片| 久久国产亚洲av麻豆专区| 好男人电影高清在线观看| 香蕉丝袜av| 多毛熟女@视频| 在线观看www视频免费| 99精品欧美一区二区三区四区| 国产欧美日韩综合在线一区二区| 91九色精品人成在线观看| 日韩 欧美 亚洲 中文字幕| 99riav亚洲国产免费| 中文字幕av电影在线播放| 人人妻人人澡欧美一区二区 | 成人精品一区二区免费| 久久国产乱子伦精品免费另类| 免费搜索国产男女视频| 黄片播放在线免费| 午夜免费成人在线视频| 男女下面插进去视频免费观看| 国产精品一区二区在线不卡| 纯流量卡能插随身wifi吗| a在线观看视频网站| 欧美日韩乱码在线| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 日韩三级视频一区二区三区| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲| 日本 av在线| 桃红色精品国产亚洲av| 一区二区三区国产精品乱码| 国产在线观看jvid| 桃红色精品国产亚洲av| 精品高清国产在线一区| 午夜福利欧美成人| 久久午夜亚洲精品久久| 9191精品国产免费久久| 国产又色又爽无遮挡免费看| 久久久久九九精品影院| 免费一级毛片在线播放高清视频 | 又大又爽又粗| 高潮久久久久久久久久久不卡| 亚洲av成人不卡在线观看播放网| cao死你这个sao货| xxx96com| 精品久久久久久,| a在线观看视频网站| 韩国av一区二区三区四区| 久久亚洲精品不卡| 中文字幕人妻熟女乱码| 女性被躁到高潮视频| 91精品三级在线观看| 在线av久久热| 亚洲电影在线观看av| 亚洲av电影在线进入| 校园春色视频在线观看| 午夜视频精品福利| 一区二区三区精品91| 黑人操中国人逼视频| avwww免费| 女警被强在线播放| 国产精品99久久99久久久不卡| 午夜福利,免费看| 亚洲免费av在线视频| 国产精品亚洲一级av第二区| 91成年电影在线观看| www.熟女人妻精品国产| 国产精品影院久久| 国产精品久久视频播放| 热re99久久国产66热| 久热这里只有精品99| 波多野结衣一区麻豆| 欧美性长视频在线观看| 日本精品一区二区三区蜜桃| 成年女人毛片免费观看观看9| 久久性视频一级片| 国产精品电影一区二区三区| avwww免费| 午夜老司机福利片| 精品少妇一区二区三区视频日本电影| videosex国产| 国产精品永久免费网站| 亚洲精品在线观看二区| 桃红色精品国产亚洲av| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 欧美中文日本在线观看视频| 国产一级毛片七仙女欲春2 | 淫秽高清视频在线观看| 淫妇啪啪啪对白视频| 村上凉子中文字幕在线| 禁无遮挡网站| 不卡av一区二区三区| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 一进一出抽搐gif免费好疼| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 一进一出好大好爽视频| 日韩欧美国产在线观看| 欧美日本亚洲视频在线播放| 禁无遮挡网站| 亚洲国产精品合色在线| 亚洲中文av在线| 欧美大码av| 一a级毛片在线观看| 精品国产一区二区久久| 成人特级黄色片久久久久久久| 午夜福利视频1000在线观看 | 日韩欧美国产在线观看| 欧美性长视频在线观看| 国产欧美日韩精品亚洲av| 欧美一级毛片孕妇| 国产精品98久久久久久宅男小说| 免费高清视频大片| 自拍欧美九色日韩亚洲蝌蚪91| 午夜视频精品福利| 亚洲人成电影观看| 精品国产国语对白av| 国产蜜桃级精品一区二区三区| 操出白浆在线播放| 精品福利观看| 淫秽高清视频在线观看| 午夜成年电影在线免费观看| 一进一出抽搐gif免费好疼| 一级片免费观看大全| 麻豆一二三区av精品| 美女午夜性视频免费| 久99久视频精品免费| 黄色成人免费大全| 成在线人永久免费视频| 亚洲情色 制服丝袜| 日本欧美视频一区| 欧美人与性动交α欧美精品济南到| www日本在线高清视频| 黄色视频不卡| 免费观看精品视频网站| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到| 欧美不卡视频在线免费观看 | 母亲3免费完整高清在线观看| 亚洲色图综合在线观看| 99re在线观看精品视频| 国产精品美女特级片免费视频播放器 | 午夜精品久久久久久毛片777| 国产精品免费视频内射| 久久久久久久精品吃奶| 搞女人的毛片| 精品国产乱码久久久久久男人| 咕卡用的链子| 露出奶头的视频| 中文字幕高清在线视频| 久久精品影院6| 午夜免费鲁丝| 美女高潮喷水抽搐中文字幕| 久久午夜综合久久蜜桃| 午夜福利一区二区在线看| 亚洲国产精品久久男人天堂| 男女下面进入的视频免费午夜 | 欧美激情极品国产一区二区三区| 欧美一级a爱片免费观看看 | 欧美中文综合在线视频| 在线视频色国产色| 老司机午夜福利在线观看视频| 99久久久亚洲精品蜜臀av| 成年人黄色毛片网站| 午夜两性在线视频| 丝袜人妻中文字幕| 变态另类成人亚洲欧美熟女 | 99久久99久久久精品蜜桃| 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 国产精品二区激情视频| 黄片播放在线免费| 国产精品二区激情视频| 亚洲情色 制服丝袜| 国产欧美日韩精品亚洲av| 男男h啪啪无遮挡| 国产av精品麻豆| 咕卡用的链子| а√天堂www在线а√下载| 99久久99久久久精品蜜桃| 国产精品日韩av在线免费观看 | av欧美777| 在线播放国产精品三级| 在线永久观看黄色视频| 国产精品一区二区在线不卡| 亚洲视频免费观看视频| 91大片在线观看| 黄片小视频在线播放| 亚洲情色 制服丝袜| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 人人妻人人爽人人添夜夜欢视频| 丰满人妻熟妇乱又伦精品不卡| 两个人免费观看高清视频| 日本三级黄在线观看| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 亚洲第一青青草原| 两个人看的免费小视频| 88av欧美| 丁香六月欧美| 一个人免费在线观看的高清视频| 午夜福利影视在线免费观看| 亚洲一码二码三码区别大吗| 亚洲av五月六月丁香网| 免费看十八禁软件|