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

    彈道目標(biāo)多傳感器低重合率航跡關(guān)聯(lián)技術(shù)研究

    2014-09-08 11:53:18朱冠宇李昌璽
    雷達(dá)與對(duì)抗 2014年3期
    關(guān)鍵詞:重合航跡彈道

    朱冠宇,周 焰,郭 戈,徐 穎,李昌璽

    (空軍預(yù)警學(xué)院,武漢 430019)

    彈道目標(biāo)多傳感器低重合率航跡關(guān)聯(lián)技術(shù)研究

    朱冠宇,周 焰,郭 戈,徐 穎,李昌璽

    (空軍預(yù)警學(xué)院,武漢 430019)

    彈道目標(biāo)伴隨著運(yùn)動(dòng)速度快、區(qū)間跨度廣、測(cè)量誤差大等因素,多傳感器在協(xié)同跟蹤時(shí)相互間因重疊率低產(chǎn)生航跡中斷是一個(gè)比較普遍的現(xiàn)象。這種現(xiàn)象因?qū)桔E關(guān)聯(lián)得到連貫完整的航跡產(chǎn)生影響而成為一個(gè)亟待解決的問題。針對(duì)該問題,通過對(duì)現(xiàn)有航跡的延續(xù)和反推,本文提出一種利用目標(biāo)狀態(tài)進(jìn)行航跡關(guān)聯(lián)的方法。實(shí)驗(yàn)證明,該方法可以有效解決彈道目標(biāo)在低重合率條件下的航跡關(guān)聯(lián)問題。

    彈道目標(biāo);航跡預(yù)測(cè);航跡關(guān)聯(lián);多傳感器

    0 引 言

    傳統(tǒng)的防空系統(tǒng)中,同一地區(qū)有多部雷達(dá)進(jìn)行監(jiān)視,傳感器之間相互重疊率高,對(duì)于同一目標(biāo)有多部雷達(dá)進(jìn)行跟蹤監(jiān)視,可以獲得不同傳感器對(duì)于同一目標(biāo) 足夠長(zhǎng)度的同一時(shí)間段內(nèi)的航跡。在對(duì)來(lái)襲目標(biāo)進(jìn)行關(guān)聯(lián)判斷時(shí),主要的根據(jù)是這些同一時(shí)間段內(nèi)的航跡之間的距離、相似程度。但是,在彈道導(dǎo)彈預(yù)警體系中,對(duì)彈道目標(biāo)的跟蹤采用分段協(xié)同探測(cè)的模式,各傳感器關(guān)于同一目標(biāo)在同一時(shí)段里航跡重疊的部分較少,甚至有時(shí)會(huì)因波束無(wú)法交接而產(chǎn)生沒有重疊的情況。這時(shí),傳統(tǒng)的航跡關(guān)聯(lián)方法就失效了。因此,需要專門尋找一種對(duì)于彈道目標(biāo)的航跡關(guān)聯(lián)的方法。

    文獻(xiàn)[4]對(duì)彈道導(dǎo)彈運(yùn)行過程中的狀態(tài)估計(jì)算法和定軌預(yù)報(bào)方法進(jìn)行了討論分析,并指出了其在遠(yuǎn)程相控陣?yán)走_(dá)與多功能地基相控陣?yán)走_(dá)之間進(jìn)行交接預(yù)報(bào)的作用。文獻(xiàn)[5]討論了在多彈道目標(biāo)發(fā)生航跡中斷之后對(duì)中斷前后的航跡片段進(jìn)行優(yōu)化的關(guān)聯(lián)方法,但基于同一傳感器的數(shù)據(jù)。本文所提供方法與前述相比是一種基于航跡預(yù)測(cè)的航跡關(guān)聯(lián)方法。將目標(biāo)航跡在兩傳感器之間的部分通過遞推來(lái)延長(zhǎng)公共部分,并在延長(zhǎng)之后利用航跡關(guān)聯(lián)方法進(jìn)行關(guān)聯(lián)判斷。通過仿真實(shí)驗(yàn)表明,該方法可以有效解決彈道目標(biāo)運(yùn)動(dòng)過程中的航跡關(guān)聯(lián)問題。

    1 彈道目標(biāo)航跡關(guān)聯(lián)基本思想

    彈道導(dǎo)彈飛行過程如圖1所示。目標(biāo)經(jīng)過不同傳感器的掃描范圍時(shí)會(huì)分別產(chǎn)生一段航跡。由于彈道導(dǎo)彈目標(biāo)飛行跨度大等特殊因素,不同傳感器產(chǎn)生的航跡之間可能會(huì)有空白或交接部分較少。大致可以分為圖1中的兩種情況。

    圖1 彈道目標(biāo)軌跡情況

    圖1中,情況1表示目標(biāo)航跡通過傳感器之間的重合區(qū)域,但處在公共部分的航跡較少;情況2表示目標(biāo)航跡在兩部傳感器間通過時(shí)有部分處在兩部傳感器的公共監(jiān)視范圍之外。

    當(dāng)情況1中傳感器之間的公共區(qū)域足夠大,則可以利用傳統(tǒng)的航跡關(guān)聯(lián)方法進(jìn)行關(guān)聯(lián);若情況2中兩部傳感器之間的空隙過大,則會(huì)影響關(guān)聯(lián)的精度。因此,需要對(duì)具體使用情況進(jìn)行判斷。

    傳統(tǒng)的航跡關(guān)聯(lián)方法只需要對(duì)有限個(gè)點(diǎn)進(jìn)行關(guān)聯(lián)處理(例如采用序貫法最少只需要6個(gè)點(diǎn))就可以完成關(guān)聯(lián)判斷。所以,對(duì)于情況1和情況2,采取不同的判斷標(biāo)準(zhǔn)。因?yàn)榍闆r1中已經(jīng)存在一部分公共航跡點(diǎn),這部分公共點(diǎn)可以利用來(lái)直接進(jìn)行關(guān)聯(lián)判斷;而情況2中,目標(biāo)在兩部傳感器之間缺失的點(diǎn)需要通過航跡遞推來(lái)得到。這樣,針對(duì)情況1,當(dāng)目標(biāo)通過的公共區(qū)域的點(diǎn)數(shù)小于或等于6個(gè)點(diǎn)時(shí)(因?yàn)樵诤桔E關(guān)聯(lián)算法中最少只需要6個(gè)點(diǎn)就可以判定兩條航跡完全關(guān)聯(lián));針對(duì)情況2,則需要缺失的點(diǎn)數(shù)小于20(缺失點(diǎn)過多將導(dǎo)致遞推誤差加大,發(fā)生漏關(guān)聯(lián)的概率增大)。

    對(duì)情況1,設(shè)n1為公共區(qū)域內(nèi)點(diǎn)的個(gè)數(shù),則判斷標(biāo)準(zhǔn)為:n1≤6。

    對(duì)情況2,設(shè)n2為兩傳感器之間缺失的點(diǎn)的個(gè)數(shù),則判斷標(biāo)準(zhǔn)為:n2<20。

    n2的確定方法如下:設(shè)t1為前航跡最后一個(gè)點(diǎn)出現(xiàn)的時(shí)刻,t2為后航跡第1個(gè)點(diǎn)出現(xiàn)的時(shí)刻,T為傳感器掃描周期,則

    (1)

    2 航跡遞推長(zhǎng)度選取

    航跡遞推的目的是為了延長(zhǎng)航跡中時(shí)間重合部分的長(zhǎng)度(情況1)或產(chǎn)生時(shí)間重合部分(情況2),并利用時(shí)間重合部分通過航跡關(guān)聯(lián)方法進(jìn)行關(guān)聯(lián)判斷。為了保證關(guān)聯(lián)正確率和運(yùn)算速度,需要確定航跡遞推的長(zhǎng)度。為方便敘述,可定義:(1)前航跡:即前一傳感器濾波得到的航跡片段;(2)后航跡:后一傳感器濾波得到的航跡片段。情況1與情況2的基本思路都是使得時(shí)間重合部分的航跡點(diǎn)數(shù)達(dá)到20個(gè),這樣可以有充分的時(shí)間重合部分的航跡來(lái)保證關(guān)聯(lián)的質(zhì)量,又可以使得航跡關(guān)聯(lián)部分不至于過于冗長(zhǎng)而影響運(yùn)行速度。

    T1、T2分別為前、后航跡,為了減小系統(tǒng)誤差,其向前、后遞推的點(diǎn)數(shù)應(yīng)相同。T1、T2遞推增加的部分航跡分別記為T10、T20,遞推的點(diǎn)數(shù)分別記為n10、n20。n1、n2定義如前所述保持不變。前、后航跡每增加一個(gè)遞推點(diǎn)即進(jìn)行一次判定, 當(dāng)判定結(jié)果第一次達(dá)到判斷標(biāo)準(zhǔn)時(shí)即停止遞推。針對(duì)情況1,判斷標(biāo)準(zhǔn)為

    n10+n20+n1≥20

    (2)

    遞推后情況1航跡如圖2所示,虛線部分為遞推航跡。圖示中,n1=4,n10=n20=8。

    圖2 情況1遞推關(guān)系示意圖

    針對(duì)情況2,判斷標(biāo)準(zhǔn)為

    n10+n20-n1≥20

    (3)

    遞推后情況2航跡如圖3所示,虛線部分為遞推航跡。圖示中,n2=4,n10=n20=12。

    圖3 情況2遞推關(guān)系示意圖

    3 航跡遞推

    進(jìn)行航跡遞推時(shí),情況1與情況2、前航跡與后航跡所使用的方法和基本原理是一致的。這里以前航跡遞推為例來(lái)進(jìn)行說明。

    由于考慮到彈道導(dǎo)彈在中段所處時(shí)間較長(zhǎng),傳感器間的跟蹤接力主要是在中段進(jìn)行,所以采用中段時(shí)的狀態(tài)進(jìn)行推測(cè)。

    設(shè)彈道導(dǎo)彈目標(biāo)在中段運(yùn)行時(shí)的方程為

    X(k+1)=FX(k)+Gu(k)+V(k)

    (4)

    其中,u(k)為目標(biāo)運(yùn)動(dòng)過程中所受到的地球加速度,V(k)為受到的其他擾動(dòng)。因?yàn)樵谥卸芜\(yùn)行時(shí)重力加速度是目標(biāo)所受到的主要影響,故可以將目標(biāo)視為僅受重力作用[1],可將此時(shí)運(yùn)動(dòng)方程簡(jiǎn)化為

    X(k+1)=FX(k)+Gu(k)

    (5)

    利用式(5)可以對(duì)目標(biāo)狀態(tài)進(jìn)行推測(cè)。對(duì)于前航跡向后遞推,可以直接使用式(5);對(duì)于后航跡的前向遞推則需要將方向進(jìn)行反向后進(jìn)行,即

    X(k)=F-1[X(k+1)-Gu(k)]

    (6)

    目標(biāo)所處位置不同,受到的引力不同。根據(jù)地球平均半徑R和導(dǎo)彈距地心位置r可得受到加速度的大小為

    (7)

    4 航跡關(guān)聯(lián)

    在這里,采用序貫法對(duì)時(shí)間重合條件下的目標(biāo)航跡(包括遞推得到的航跡)進(jìn)行航跡關(guān)聯(lián)處理。

    航跡關(guān)聯(lián)算法主要分為兩大類,即傳統(tǒng)的基于統(tǒng)計(jì)的算法和基于模糊數(shù)學(xué)的算法。其中序貫法是基于統(tǒng)計(jì)的航跡關(guān)聯(lián)算法的一種。序貫法是把航跡當(dāng)前時(shí)刻的關(guān)聯(lián)與其歷史聯(lián)系起來(lái),并賦予良好的航跡關(guān)聯(lián)質(zhì)量管理和多義性處理技術(shù),能夠應(yīng)對(duì)一定的密集目標(biāo)環(huán)境、交叉、分岔航跡較多的場(chǎng)合。同時(shí),相比于模糊數(shù)學(xué)的方法,其原理簡(jiǎn)單,運(yùn)算速度快,且由于傳感器精度的提高其關(guān)聯(lián)精度已經(jīng)完全能夠滿足需要[2]。

    將時(shí)間重合部分航跡重新編號(hào)。假設(shè)該部分共有n個(gè)點(diǎn),則前后航跡該部分的狀態(tài)估計(jì)之差記為

    (8)

    兩航跡關(guān)聯(lián)實(shí)際上是一種假設(shè)檢驗(yàn)過程。設(shè)H0和H1是下列事件:

    H0:前、后航跡關(guān)聯(lián);

    H1:前、后航跡不關(guān)聯(lián)。

    序貫法采取的假設(shè)檢驗(yàn)量為

    (9)

    其中

    C12(k|k)=P1(k|k)+P2(k|k)

    (10)

    其中P1(k|k)與P2(k|k)分別為前、后航跡的狀態(tài)協(xié)方差。由于在航跡外推過程中兩部雷達(dá)沒有獲得測(cè)量值,這里在應(yīng)用到外推航跡時(shí)對(duì)P1(k|k)與P2(k|k)計(jì)算利用線性卡爾曼濾波方法進(jìn)行近似:協(xié)方差一步預(yù)測(cè)為

    P(k+1|k)=F(k)P(k|k)F′(k)+Q(k)

    (11)

    信息協(xié)方差為

    (12)

    這里

    R(k+1)由X(k+1)轉(zhuǎn)化出斜距、俯仰、方位信息后近似求得,繼而

    K(k+1)=P(k+1|k)H′(k+1)S-1(k+1)

    (13)

    (14)

    λ12(k)是具有knx自由度的χ2分布隨即變量,nx為狀態(tài)向量的維數(shù)。這樣即可對(duì)H0或H1進(jìn)行假設(shè)檢驗(yàn),即如果

    λ12(k)≤δ(k)

    (15)

    則接受H0,否則接受H1。其中閾值滿足

    Pr{λ12(k)>δ(k)|H0}=α

    (16)

    式中,α為檢驗(yàn)的顯著水平,這里一般取0.05。于是,對(duì)于給定的α,δ(k)隨著k的增加而變大,且需實(shí)時(shí)修正。對(duì)于某一nx值,當(dāng)k大于某一值時(shí),δ(k)近似為一遞推式。例如,nx=4,α=0.05,當(dāng)k≥8時(shí):

    δ(k)=δ(k-1)+5

    (17)

    于是,便可以利用其判斷出前、后航跡是否關(guān)聯(lián)。再對(duì)時(shí)間重合部分的航跡進(jìn)行航跡融合,并與前、后航跡一起組成了一條完整、連貫的航跡。

    5 仿真實(shí)驗(yàn)

    實(shí)驗(yàn)?zāi)M某型射程為5000 km的彈道導(dǎo)彈發(fā)射,導(dǎo)彈發(fā)射點(diǎn)位于147°E, 18°N。

    濾波器采用的是無(wú)味卡爾曼濾波器(Unscented Kalman Filter, UKF),主要原因,一是不需要對(duì)非線性狀態(tài)方程和觀測(cè)方程進(jìn)行線性化處理,直接利用非線性狀態(tài)方程來(lái)估算目標(biāo)狀態(tài)向量的概率密度函數(shù),適用于任意非線性模型,狀態(tài)估計(jì)精度較高;二是UKF在使用時(shí)不需要計(jì)算雅克比矩陣,實(shí)現(xiàn)簡(jiǎn)單,時(shí)效性高[3]。

    坐標(biāo)系采用地球固連坐標(biāo)系,目標(biāo)的三維坐標(biāo)在該坐標(biāo)系下產(chǎn)生。

    5.1 單目標(biāo)跟蹤仿真

    兩部雷達(dá)均采用遠(yuǎn)程相控陣?yán)走_(dá),測(cè)距精度為50 m,方位角和俯仰角測(cè)量噪聲標(biāo)準(zhǔn)差為0.2°[4],均服從均值為0的正態(tài)分布。其中,雷達(dá)還存在1×10-3m與1×10-3rad的系統(tǒng)誤差。采樣時(shí)間T=1 s。

    導(dǎo)彈在發(fā)射之后共歷時(shí)1538 s后命中目標(biāo)。兩部雷達(dá)在其中的探測(cè)區(qū)間分別是[1100 s,1200 s],[1190 s,1290 s]。

    兩部雷達(dá)對(duì)目標(biāo)跟蹤濾波后的圖像見圖4。經(jīng)過航跡遞推后,得到新航跡與之前航跡結(jié)合,如圖5所示。

    圖4 遞推前航跡

    圖5 遞推后航跡

    圖6中“*”表示為前航跡向后遞推得到的部分,“+”表示后航跡前向遞推得到的部分,其中遞推部分放大之后的圖像可見圖7。再利用航跡關(guān)聯(lián)的方法,即可判斷出兩部分航跡屬于同一目標(biāo)。

    圖6 遞推部分詳細(xì)表示

    5.2 群目標(biāo)仿真

    在5.1的基礎(chǔ)上,假設(shè)同時(shí)有10批次彈道導(dǎo)彈目標(biāo)進(jìn)入兩部雷達(dá)威力區(qū)域。進(jìn)行多次關(guān)聯(lián)判斷仿真,其中一次的目標(biāo)在雷達(dá)威力區(qū)域軌跡圖如圖7所示。

    圖7 群目標(biāo)運(yùn)行軌跡

    此時(shí),兩部雷達(dá)對(duì)群目標(biāo)探測(cè)的航跡圖像及航跡遞推后關(guān)聯(lián)圖像如圖8、圖9所示。

    圖8 兩部雷達(dá)分別測(cè)得目標(biāo)軌跡

    圖9 遞推后航跡

    從上述圖中可以看出目標(biāo)在雷達(dá)威力區(qū)域明顯的航跡中斷現(xiàn)象,再通過遞推關(guān)聯(lián)判斷等方法處理后可以得出多目標(biāo)航跡的關(guān)聯(lián)情況。

    針對(duì)上述多目標(biāo)情況進(jìn)行30次仿真,結(jié)果如表1所示。

    表1 多目標(biāo)關(guān)聯(lián)效果

    6 結(jié)束語(yǔ)

    本文提出了一種利用目標(biāo)狀態(tài)進(jìn)行航跡關(guān)聯(lián)的方法。仿真結(jié)果證實(shí),該方法關(guān)聯(lián)正確率較高,能夠較好地解決在多傳感器低重合率情況下的彈道目標(biāo)航跡關(guān)聯(lián)問題。這對(duì)后續(xù)的航跡識(shí)別、平滑處理以及航跡管理具有較大的參考價(jià)值。

    [1] 唐宏斌,鄭健,李駿,姜文利,周宇.基于主動(dòng)段狀態(tài)估計(jì)的彈道估計(jì)方法研究[J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù),2007(3):27-31.

    [2] 何友,修建娟,張晶煒,等.雷達(dá)數(shù)據(jù)處理及應(yīng)用[M].北京:電子工業(yè)出版社,2005:26-28,139-144.

    [3] 吳紅,等.彈道導(dǎo)彈自由段跟蹤初始狀態(tài)估計(jì)方法研究[J]. 計(jì)算機(jī)仿真,2006(6).

    [4] 張榮濤.多雷達(dá)跟蹤彈道導(dǎo)彈交接預(yù)報(bào)技術(shù)研究[J]. 現(xiàn)代雷達(dá),2010(8).

    [5] 俞建國(guó),劉梅,陳錦海, 等.彈道目標(biāo)航跡片段關(guān)聯(lián)及優(yōu)化[J]. 航空學(xué)報(bào),2011,32(10).

    [6] William J F.Interacting multiple model filter for tactical ballistic missile tracking[J].IEEE Trans. on Aerospace and Electronic Systems,2008,44(2):418-426.

    Research on track correlation technology of ballistic targets for multiple sensors with low overlapping rate

    ZHU Guan-yu, ZHOU Yan, GUO Ge, XU Ying, LI Chang-xi

    (Air Force Early Warning Academy, Wuhan 430019)

    On account of such factors as high moving speed, wide span space and big measurement error for ballistic targets, the track interrupt caused by low overlapping rate is a relatively common phenomenon when multiple sensors are tracking the targets cooperatively and it is a serious problem that needs to be solved. Aimed at this problem, a method of using the target state to perform track correlation is proposed through extending the existing tracks. The test results show that the method can effectively solve the problem of track correlation of ballistic targets with low overlapping rate.

    ballistic target; track prediction; track correlation; multiple sensors

    2014-05-03;

    2014-05-16

    朱冠宇(1990-),男,碩士研究生,研究方向:信息融合;周焰(1966-),男,教授,研究方向:信息融合與圖像處理;郭戈(1985-),男,博士研究生,研究方向:信息融合;徐穎(1985-),女,博士研究生,研究方向:信息融合與圖像處理;李昌璽(1986-),男,博士研究生,研究方向:信息融合。

    TN957.522

    A

    1009-0401(2014)03-0056-05

    猜你喜歡
    重合航跡彈道
    彈道——打勝仗的奧秘
    夢(mèng)的航跡
    青年歌聲(2019年12期)2019-12-17 06:32:32
    一維彈道修正彈無(wú)線通信系統(tǒng)研制
    電子制作(2019年7期)2019-04-25 13:17:48
    自適應(yīng)引導(dǎo)長(zhǎng)度的無(wú)人機(jī)航跡跟蹤方法
    電力系統(tǒng)單回線自適應(yīng)重合閘的研究
    電子制作(2017年10期)2017-04-18 07:23:07
    視覺導(dǎo)航下基于H2/H∞的航跡跟蹤
    基于PID控制的二維彈道修正彈仿真
    考慮暫態(tài)穩(wěn)定優(yōu)化的自適應(yīng)重合閘方法
    消除彈道跟蹤數(shù)據(jù)中伺服系統(tǒng)的振顫干擾
    基于航跡差和航向差的航跡自動(dòng)控制算法
    久久久久久人人人人人| 变态另类成人亚洲欧美熟女| 久久国产精品影院| 日本 av在线| 中文字幕熟女人妻在线| 色视频www国产| 精品国内亚洲2022精品成人| 身体一侧抽搐| 哪里可以看免费的av片| 18禁裸乳无遮挡免费网站照片| 久99久视频精品免费| 亚洲成人免费电影在线观看| 午夜精品一区二区三区免费看| 国产精品一区二区免费欧美| 丁香六月欧美| 老汉色∧v一级毛片| 日本免费一区二区三区高清不卡| 一二三四在线观看免费中文在| 天天一区二区日本电影三级| 亚洲午夜理论影院| 国产激情久久老熟女| 成人三级做爰电影| 国产精品美女特级片免费视频播放器 | 日韩欧美国产在线观看| 国内精品久久久久久久电影| 99riav亚洲国产免费| 中文字幕精品亚洲无线码一区| 最近最新中文字幕大全电影3| 欧美精品啪啪一区二区三区| 日本五十路高清| 天堂动漫精品| 18禁裸乳无遮挡免费网站照片| 女人被狂操c到高潮| 丁香六月欧美| 亚洲熟女毛片儿| 日韩 欧美 亚洲 中文字幕| 男女做爰动态图高潮gif福利片| 亚洲中文av在线| 亚洲精品美女久久久久99蜜臀| 麻豆国产97在线/欧美| 欧美日韩乱码在线| 国产淫片久久久久久久久 | 久久久精品大字幕| 淫妇啪啪啪对白视频| 久久午夜综合久久蜜桃| 很黄的视频免费| 免费看美女性在线毛片视频| 成人性生交大片免费视频hd| 舔av片在线| 国产亚洲欧美在线一区二区| 黄色女人牲交| 久久久久久久久免费视频了| 欧美成人免费av一区二区三区| 久久久久性生活片| 国产高清三级在线| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 亚洲精品在线观看二区| 热99在线观看视频| 综合色av麻豆| 香蕉国产在线看| 久久天躁狠狠躁夜夜2o2o| 免费观看的影片在线观看| 18美女黄网站色大片免费观看| 中亚洲国语对白在线视频| 午夜福利在线观看吧| 欧美乱妇无乱码| 久久性视频一级片| 国产成人欧美在线观看| 精品一区二区三区av网在线观看| 亚洲成av人片在线播放无| 又黄又爽又免费观看的视频| 99久久久亚洲精品蜜臀av| 少妇丰满av| 国产成年人精品一区二区| 观看美女的网站| 国产高清有码在线观看视频| 99精品久久久久人妻精品| 久久国产精品人妻蜜桃| 天天一区二区日本电影三级| 久久久国产欧美日韩av| 久久久国产精品麻豆| 日本 欧美在线| 欧美3d第一页| 日本一本二区三区精品| 亚洲无线观看免费| 看免费av毛片| 搡老熟女国产l中国老女人| 在线播放国产精品三级| 一个人观看的视频www高清免费观看 | 免费大片18禁| 91在线精品国自产拍蜜月 | 99久久国产精品久久久| 国产人伦9x9x在线观看| 一本综合久久免费| 亚洲专区中文字幕在线| 欧美黄色片欧美黄色片| 两个人视频免费观看高清| 亚洲欧美日韩高清专用| 中国美女看黄片| 国产精品亚洲一级av第二区| 色哟哟哟哟哟哟| 村上凉子中文字幕在线| 在线永久观看黄色视频| 精品国产乱码久久久久久男人| 男人和女人高潮做爰伦理| 成年免费大片在线观看| 中文字幕最新亚洲高清| 特级一级黄色大片| 人妻夜夜爽99麻豆av| 久久婷婷人人爽人人干人人爱| 99国产综合亚洲精品| 日韩欧美国产在线观看| 一本精品99久久精品77| 亚洲av片天天在线观看| 小说图片视频综合网站| 精华霜和精华液先用哪个| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 人妻丰满熟妇av一区二区三区| 亚洲人成网站高清观看| 老熟妇乱子伦视频在线观看| 一级a爱片免费观看的视频| 成年女人永久免费观看视频| 99精品久久久久人妻精品| 午夜精品一区二区三区免费看| 亚洲成人免费电影在线观看| 国产精品av视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲七黄色美女视频| 国产亚洲av嫩草精品影院| 亚洲欧美激情综合另类| 在线播放国产精品三级| 中文字幕高清在线视频| 国内精品美女久久久久久| 美女午夜性视频免费| 国产蜜桃级精品一区二区三区| 成人一区二区视频在线观看| 亚洲天堂国产精品一区在线| 超碰成人久久| 黄色 视频免费看| 人人妻,人人澡人人爽秒播| 色吧在线观看| 香蕉丝袜av| aaaaa片日本免费| 国产欧美日韩精品一区二区| 在线a可以看的网站| 人人妻人人澡欧美一区二区| 国产不卡一卡二| 国产成年人精品一区二区| 欧美日韩中文字幕国产精品一区二区三区| 日本精品一区二区三区蜜桃| 999精品在线视频| 亚洲成人精品中文字幕电影| 亚洲欧美日韩无卡精品| 欧美日韩中文字幕国产精品一区二区三区| 成人午夜高清在线视频| 成年人黄色毛片网站| 97超级碰碰碰精品色视频在线观看| 麻豆国产av国片精品| 午夜激情福利司机影院| 好男人在线观看高清免费视频| 最新在线观看一区二区三区| 免费一级毛片在线播放高清视频| 一区二区三区高清视频在线| 亚洲精品美女久久久久99蜜臀| 男女视频在线观看网站免费| 成熟少妇高潮喷水视频| 欧美日本视频| 亚洲自拍偷在线| 国产精品九九99| 国产亚洲精品久久久com| 法律面前人人平等表现在哪些方面| 亚洲最大成人中文| 久久久久亚洲av毛片大全| 色吧在线观看| 淫秽高清视频在线观看| 国产精品久久久久久人妻精品电影| 亚洲人与动物交配视频| 狂野欧美白嫩少妇大欣赏| 91老司机精品| 欧美乱码精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | ponron亚洲| 国产精品久久久人人做人人爽| 观看美女的网站| 在线观看免费午夜福利视频| 日韩欧美在线二视频| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 88av欧美| 国产成人影院久久av| 亚洲avbb在线观看| 欧美乱色亚洲激情| 身体一侧抽搐| 99久久精品热视频| 69av精品久久久久久| 欧美色欧美亚洲另类二区| 国产aⅴ精品一区二区三区波| 白带黄色成豆腐渣| 亚洲精品在线观看二区| 国产精品九九99| 青草久久国产| 日韩欧美国产在线观看| 日韩国内少妇激情av| 亚洲自拍偷在线| 欧美日韩乱码在线| av女优亚洲男人天堂 | 亚洲真实伦在线观看| 十八禁网站免费在线| 欧美激情在线99| www国产在线视频色| 国产欧美日韩精品一区二区| 欧美三级亚洲精品| 国产欧美日韩一区二区三| 欧美在线一区亚洲| 国产日本99.免费观看| 人妻丰满熟妇av一区二区三区| 亚洲人成网站在线播放欧美日韩| 日本与韩国留学比较| 琪琪午夜伦伦电影理论片6080| 亚洲精品久久国产高清桃花| 精品无人区乱码1区二区| 国产激情欧美一区二区| 亚洲电影在线观看av| 色噜噜av男人的天堂激情| 日本免费a在线| 两人在一起打扑克的视频| 亚洲国产中文字幕在线视频| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 欧美中文日本在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 久久草成人影院| 精品无人区乱码1区二区| 男人和女人高潮做爰伦理| 性欧美人与动物交配| 淫妇啪啪啪对白视频| 啦啦啦观看免费观看视频高清| 激情在线观看视频在线高清| 国产精品1区2区在线观看.| 日韩免费av在线播放| 久久人人精品亚洲av| 久久久国产欧美日韩av| 美女免费视频网站| 岛国视频午夜一区免费看| 性色av乱码一区二区三区2| 99久久综合精品五月天人人| 日日摸夜夜添夜夜添小说| 不卡av一区二区三区| 动漫黄色视频在线观看| 国产乱人视频| 巨乳人妻的诱惑在线观看| 99国产综合亚洲精品| 高清毛片免费观看视频网站| 婷婷丁香在线五月| www.熟女人妻精品国产| 天堂av国产一区二区熟女人妻| 看免费av毛片| 国产精品永久免费网站| 久久草成人影院| 成年免费大片在线观看| 国内精品久久久久精免费| 日韩国内少妇激情av| 成人永久免费在线观看视频| 亚洲自拍偷在线| 国产男靠女视频免费网站| 国产精品98久久久久久宅男小说| 成人高潮视频无遮挡免费网站| 亚洲va日本ⅴa欧美va伊人久久| 黑人欧美特级aaaaaa片| 亚洲av五月六月丁香网| 亚洲欧美日韩东京热| 少妇的逼水好多| 精品福利观看| 中文亚洲av片在线观看爽| 老汉色av国产亚洲站长工具| bbb黄色大片| 一级毛片精品| 欧美又色又爽又黄视频| 亚洲精品中文字幕一二三四区| 午夜福利成人在线免费观看| 手机成人av网站| 九九久久精品国产亚洲av麻豆 | 真实男女啪啪啪动态图| 一夜夜www| 国产激情久久老熟女| 欧美激情久久久久久爽电影| 在线永久观看黄色视频| 丝袜人妻中文字幕| 久久久精品大字幕| 亚洲一区二区三区色噜噜| 亚洲激情在线av| 91在线精品国自产拍蜜月 | 男插女下体视频免费在线播放| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影| 国产av在哪里看| 国产精品影院久久| 91av网一区二区| 久久久久久久久久黄片| 久久香蕉国产精品| 亚洲国产欧美人成| 51午夜福利影视在线观看| 一边摸一边抽搐一进一小说| 免费一级毛片在线播放高清视频| 国产乱人伦免费视频| 亚洲熟女毛片儿| 曰老女人黄片| 在线永久观看黄色视频| 久久久久国内视频| 黄色丝袜av网址大全| 亚洲av片天天在线观看| 欧美日韩综合久久久久久 | 国产精品乱码一区二三区的特点| 久久久久国产精品人妻aⅴ院| 97超视频在线观看视频| 日韩欧美免费精品| 色av中文字幕| 全区人妻精品视频| 亚洲片人在线观看| 毛片女人毛片| 免费一级毛片在线播放高清视频| 大型黄色视频在线免费观看| 久久久久久久午夜电影| 精品日产1卡2卡| 国内精品一区二区在线观看| 国产精品久久久久久精品电影| 久99久视频精品免费| 激情在线观看视频在线高清| 波多野结衣高清作品| 成年女人毛片免费观看观看9| 熟女电影av网| 中文字幕熟女人妻在线| 蜜桃久久精品国产亚洲av| 成人av一区二区三区在线看| 亚洲精品中文字幕一二三四区| 两个人看的免费小视频| 亚洲美女黄片视频| 十八禁人妻一区二区| 搡老熟女国产l中国老女人| 久久天堂一区二区三区四区| 精品一区二区三区四区五区乱码| 白带黄色成豆腐渣| 国产精品九九99| 亚洲国产精品久久男人天堂| 国产一区二区三区视频了| 日韩高清综合在线| 熟女电影av网| 成年人黄色毛片网站| 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 麻豆成人av在线观看| 中亚洲国语对白在线视频| 欧美日韩精品网址| 老鸭窝网址在线观看| 老熟妇乱子伦视频在线观看| 久久天躁狠狠躁夜夜2o2o| h日本视频在线播放| 欧美zozozo另类| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人 | 麻豆成人av在线观看| 免费看十八禁软件| 欧美三级亚洲精品| 麻豆成人午夜福利视频| ponron亚洲| 日韩有码中文字幕| 99国产极品粉嫩在线观看| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 女人被狂操c到高潮| 久久国产乱子伦精品免费另类| 美女高潮的动态| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 巨乳人妻的诱惑在线观看| 婷婷丁香在线五月| 久久精品国产亚洲av香蕉五月| 欧美国产日韩亚洲一区| 97超视频在线观看视频| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 欧美乱码精品一区二区三区| 亚洲精品中文字幕一二三四区| 757午夜福利合集在线观看| 国产成人精品久久二区二区91| 日韩精品中文字幕看吧| a级毛片在线看网站| 黄色日韩在线| 一本久久中文字幕| 国产高清视频在线播放一区| 亚洲欧美激情综合另类| 精华霜和精华液先用哪个| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久久久精品电影| 久久午夜综合久久蜜桃| 九色成人免费人妻av| 亚洲av中文字字幕乱码综合| 国产精品自产拍在线观看55亚洲| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人 | 成人av一区二区三区在线看| x7x7x7水蜜桃| 巨乳人妻的诱惑在线观看| 亚洲男人的天堂狠狠| 久久久成人免费电影| 观看免费一级毛片| 女生性感内裤真人,穿戴方法视频| 中文字幕高清在线视频| 一夜夜www| 成年版毛片免费区| 久久性视频一级片| 久久精品人妻少妇| 桃色一区二区三区在线观看| avwww免费| 伦理电影免费视频| 亚洲熟妇熟女久久| 亚洲无线观看免费| 母亲3免费完整高清在线观看| 小说图片视频综合网站| 久久久国产成人免费| 又粗又爽又猛毛片免费看| 亚洲人成伊人成综合网2020| 五月伊人婷婷丁香| 欧美大码av| 亚洲熟妇中文字幕五十中出| 亚洲av片天天在线观看| 亚洲av熟女| 一二三四在线观看免费中文在| 国产精华一区二区三区| 精品电影一区二区在线| 久久精品国产99精品国产亚洲性色| 午夜福利18| 非洲黑人性xxxx精品又粗又长| 激情在线观看视频在线高清| 国产精品亚洲美女久久久| 九九久久精品国产亚洲av麻豆 | 亚洲avbb在线观看| 黄色 视频免费看| 一边摸一边抽搐一进一小说| 看片在线看免费视频| 国产私拍福利视频在线观看| 国产v大片淫在线免费观看| 国产乱人视频| 香蕉国产在线看| 欧美日韩综合久久久久久 | 91av网站免费观看| 丁香六月欧美| 国产高清视频在线观看网站| 高清在线国产一区| 中文字幕人妻丝袜一区二区| 精品国产乱子伦一区二区三区| 欧美一区二区国产精品久久精品| 精品国产乱子伦一区二区三区| 国产私拍福利视频在线观看| 黄色丝袜av网址大全| 国产伦精品一区二区三区视频9 | 久久精品亚洲精品国产色婷小说| 一个人观看的视频www高清免费观看 | 久久婷婷人人爽人人干人人爱| 国产成人av激情在线播放| 国产一区二区在线av高清观看| 9191精品国产免费久久| av视频在线观看入口| 国产亚洲精品久久久久久毛片| 亚洲人成电影免费在线| 亚洲成人中文字幕在线播放| 久久热在线av| 露出奶头的视频| 欧美激情久久久久久爽电影| 一本一本综合久久| 村上凉子中文字幕在线| 国产人伦9x9x在线观看| 久久久久久大精品| 老熟妇乱子伦视频在线观看| 99在线视频只有这里精品首页| 2021天堂中文幕一二区在线观| 午夜成年电影在线免费观看| 一二三四社区在线视频社区8| 搡老熟女国产l中国老女人| 夜夜爽天天搞| ponron亚洲| 精品久久蜜臀av无| 精品乱码久久久久久99久播| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 亚洲成人久久性| 欧美丝袜亚洲另类 | 精品国产乱子伦一区二区三区| 天天添夜夜摸| 国产精品99久久99久久久不卡| 久9热在线精品视频| 亚洲成人精品中文字幕电影| 最近最新中文字幕大全电影3| 久久久久久大精品| 国语自产精品视频在线第100页| 后天国语完整版免费观看| 日本与韩国留学比较| 国内久久婷婷六月综合欲色啪| 日韩三级视频一区二区三区| 欧美日韩瑟瑟在线播放| 两性夫妻黄色片| 亚洲乱码一区二区免费版| 国产高清videossex| 久久精品综合一区二区三区| 亚洲色图av天堂| 后天国语完整版免费观看| 色在线成人网| 中文在线观看免费www的网站| 综合色av麻豆| 免费av不卡在线播放| 精品熟女少妇八av免费久了| 国产不卡一卡二| 国内精品久久久久精免费| 免费观看人在逋| 夜夜看夜夜爽夜夜摸| 亚洲人成伊人成综合网2020| 在线永久观看黄色视频| 亚洲国产欧美网| 夜夜爽天天搞| 又粗又爽又猛毛片免费看| 国产精品久久电影中文字幕| 国产三级黄色录像| 激情在线观看视频在线高清| 熟女人妻精品中文字幕| 日本 av在线| 看片在线看免费视频| 夜夜爽天天搞| 嫩草影视91久久| 精品99又大又爽又粗少妇毛片 | 亚洲最大成人中文| 欧美色视频一区免费| 巨乳人妻的诱惑在线观看| 91久久精品国产一区二区成人 | 欧美av亚洲av综合av国产av| 天堂网av新在线| 亚洲七黄色美女视频| 黄色 视频免费看| 人妻久久中文字幕网| 九色成人免费人妻av| 日韩精品中文字幕看吧| 国产成人aa在线观看| 久久亚洲精品不卡| 在线国产一区二区在线| 午夜福利在线观看吧| 国产精品影院久久| 最新美女视频免费是黄的| 免费大片18禁| 身体一侧抽搐| 亚洲avbb在线观看| 午夜福利免费观看在线| 小说图片视频综合网站| 国产69精品久久久久777片 | 国产精品99久久久久久久久| 久久久色成人| av欧美777| 亚洲中文日韩欧美视频| 十八禁人妻一区二区| 久久国产乱子伦精品免费另类| 夜夜躁狠狠躁天天躁| av视频在线观看入口| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 久久久精品欧美日韩精品| 国产精品久久久久久久电影 | 国产精品99久久久久久久久| 欧美色视频一区免费| 最新中文字幕久久久久 | 亚洲国产欧洲综合997久久,| 国产毛片a区久久久久| 黄色片一级片一级黄色片| 色综合亚洲欧美另类图片| 好看av亚洲va欧美ⅴa在| 亚洲成av人片免费观看| 色噜噜av男人的天堂激情| 中文字幕高清在线视频| 高清毛片免费观看视频网站| 91av网一区二区| 日日摸夜夜添夜夜添小说| 欧美xxxx黑人xx丫x性爽| 久久人人精品亚洲av| 人妻丰满熟妇av一区二区三区| 亚洲国产精品成人综合色| 久久久久亚洲av毛片大全| 日本一本二区三区精品| 99国产精品一区二区蜜桃av| 神马国产精品三级电影在线观看| 日本一本二区三区精品| 69av精品久久久久久| 此物有八面人人有两片| 黑人欧美特级aaaaaa片| 69av精品久久久久久| 亚洲精品在线观看二区| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 2021天堂中文幕一二区在线观| 人妻丰满熟妇av一区二区三区| 一个人观看的视频www高清免费观看 | 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩综合久久久久久 | 午夜福利在线在线| 九色国产91popny在线| 成人午夜高清在线视频| 欧美又色又爽又黄视频| 久久久精品大字幕| a级毛片在线看网站| 亚洲自拍偷在线| 国产av麻豆久久久久久久| 一本久久中文字幕| 国产不卡一卡二| 国产熟女xx| 日韩有码中文字幕| 精品乱码久久久久久99久播| 曰老女人黄片|