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

    僅利用Crab脈沖星的徑向/橫向測(cè)距方法

    2023-08-22 05:05:16賈運(yùn)澤王奕迪
    宇航學(xué)報(bào) 2023年7期
    關(guān)鍵詞:脈沖星畸變測(cè)距

    劉 勁,賈運(yùn)澤,王奕迪,馬 辛

    (1. 武漢科技大學(xué)信息科學(xué)與工程學(xué)院,武漢 430081;2. 國(guó)防科技大學(xué)空天科學(xué)學(xué)院,長(zhǎng)沙 410003;3. 北京航空航天大學(xué)儀器科學(xué)與光電工程學(xué)院,北京 100191)

    0 引 言

    X射線脈沖星能為航天器提供徑向上的速度與位置等信息,對(duì)深空探測(cè)自主導(dǎo)航具有重要意義。因而,X射線脈沖星導(dǎo)航備受?chē)?guó)內(nèi)外學(xué)者關(guān)注[1-4]。

    自從1981年Chester等[5]提出X射線脈沖星導(dǎo)航的概念以來(lái),美國(guó)就率先開(kāi)展相關(guān)研究工作[6],并于2017年發(fā)射了中子星內(nèi)部組成探測(cè)器[7-8],完成了X射線脈沖星導(dǎo)航空間驗(yàn)證。

    從2006年開(kāi)始,中國(guó)緊隨其后,開(kāi)展了X射線脈沖星導(dǎo)航研究工作,至今取得了一系列研究成果[9-11]。中國(guó)先后發(fā)射了天宮二號(hào)[12]、脈沖星導(dǎo)航實(shí)驗(yàn)衛(wèi)星[13](XPAVN-1)、“慧眼”探測(cè)器[14],驗(yàn)證了脈沖星導(dǎo)航應(yīng)用的可行性。

    鑒于Crab脈沖星的流量遠(yuǎn)高于其他脈沖星,國(guó)內(nèi)外X射線脈沖星導(dǎo)航實(shí)驗(yàn)的觀測(cè)脈沖星都是Crab脈沖星(RSP B0531+21)。從Crab脈沖星累積輪廓中,不僅可以估計(jì)脈沖星到達(dá)時(shí)間,還可感知脈沖星累積輪廓的微弱畸變,以反演航天器在脈沖星徑向上的速度信息。鑒于其他脈沖星流量極低,即使航天器搭載大面積X射線敏感器,僅能從接收的X射線脈沖星信號(hào)中估計(jì)出低精度的脈沖星到達(dá)時(shí)間,并且無(wú)法感知低流量的脈沖星累積輪廓的極微弱畸變。

    1 脈沖星橫向測(cè)距的科學(xué)機(jī)理

    在非線性軌道中,由于位置誤差傳播,除徑向誤差外,橫向位置誤差也會(huì)引起徑向位置誤差,進(jìn)而引起X射線脈沖星累積輪廓的畸變。本節(jié)將對(duì)此進(jìn)行證明,并推導(dǎo)出位置誤差與脈沖星輪廓畸變之間的關(guān)系,以此證明基于脈沖星畸變輪廓的橫向測(cè)距方法的可行性。

    圖1給出了關(guān)于脈沖星徑向與橫面、軌道面、軌道投影、初始位置以及OTD的示意圖。其中,軌道投影是指軌道在橫面上的投影;位置角是以航天器飛行方向?yàn)檎?當(dāng)前位置相對(duì)于初始位置的轉(zhuǎn)動(dòng)角。如:初始位置角為0°;飛行1/2個(gè)軌道周期時(shí),位置角為180°;飛行一個(gè)軌道周期時(shí),位置角為360°。

    圖1 OTD示意圖

    1.1 非線性軌道下徑向上的位置誤差

    由于位置誤差的傳播性質(zhì),不同方向上的位置誤差會(huì)互相影響,即在脈沖星橫面上的位置誤差必然會(huì)引起脈沖星徑向上的位置誤差,從而引起脈沖星輪廓畸變。下面,將證明這一性質(zhì)。

    位置誤差傳播性質(zhì)1.在非線性軌道中,橫向位置誤差經(jīng)傳播導(dǎo)致徑向位置誤差。

    證.采用反證法。

    在地心坐標(biāo)系(J2000.0)下,軌道動(dòng)力學(xué)模型如式(1):

    (1)

    式中:X為狀態(tài)向量,w(t)為系統(tǒng)噪聲。

    將式(1)離散化可得到:

    Xi=Φi,i-1Xi-1+Wi-1

    (2)

    式中:Φi,i-1≈Fi-1·ΔT+I6×6;I為單位矩陣;ΔT為采樣周期;Fi-1為ti-1時(shí)刻的狀態(tài)轉(zhuǎn)移方程的雅可比矩陣,在忽略攝動(dòng)項(xiàng)后,其表達(dá)式為:

    (3)

    令X0為航天器的初始狀態(tài)向量,忽略噪聲項(xiàng)和非線性量化誤差,則有:

    Xi=Φi,i-1·Φi-1,i-2…Φ1,0·X0

    (4)

    (5)

    (6)

    其在t1時(shí)刻的具體表達(dá)式為:

    (7)

    式中:δr0為t0時(shí)刻的初始位置誤差。

    同理,可得t2,t3,t4時(shí)刻的具體表達(dá)式,如式(8)~(10)所示。

    (8)

    (9)

    S0S2ΔT4]·δr0

    (10)

    聯(lián)立式(7)~(10),可得方程組。

    若徑向位置誤差保持為零,即

    (11)

    脈沖星累積輪廓不會(huì)發(fā)生畸變。此時(shí),唯一解為δr0=03×1。這就說(shuō)明,若徑向位置誤差始終為0,初始位置誤差必為0。反之,初始位置誤差(包括初始橫向位置誤差)必會(huì)引起徑向位置誤差。

    位置誤差傳播性質(zhì)得證。

    證.采用反證法。

    式(8)減式(7),式(9)減式(8),式(10)減式(9),可得方程組(12):

    (12)

    若徑向位置誤差保持不變,即

    (13)

    推論得證。

    1.2 位置誤差引起脈沖星輪廓畸變的形成機(jī)理

    本節(jié)研究脈沖星輪廓畸變的形成機(jī)理。本質(zhì)上,脈沖星累積輪廓可視為多個(gè)脈沖星子輪廓的疊加,具體如下:

    為簡(jiǎn)化問(wèn)題,忽略噪聲干擾,僅考慮相移引起的脈沖星輪廓畸變。設(shè)歸一化的脈沖星標(biāo)準(zhǔn)輪廓為h(φ)。其中,φ為脈沖相位;h(φ)的周期為1,即h(φ)=h(φ+1)。在脈沖星信號(hào)累積過(guò)程中,利用預(yù)估的徑向位置對(duì)相移進(jìn)行補(bǔ)償。若存在徑向位置誤差,必然導(dǎo)致相位發(fā)生漂移。相移量Δφj為:

    (14)

    式中:T0為脈沖星固有周期;c為光速;Δφj為第j(j=0,1,2,…,J-1)個(gè)周期的徑向位置誤差引起的相位誤差;J為脈沖星累積時(shí)間內(nèi)的脈沖周期數(shù)。

    (15)

    根據(jù)1.1節(jié)中的位置誤差傳播性質(zhì)及其推論,由式(15)可知,脈沖星累積輪廓是由不同相位的脈沖星子輪廓疊加而成,必然發(fā)生畸變。

    若要使脈沖星累積輪廓不發(fā)生畸變,徑向位置誤差必須始終為零。根據(jù)1.1節(jié)的位置誤差傳播性質(zhì),徑向位置誤差若要始終為零,橫向位置誤差也必須為零;反之,橫向位置誤差必然引起徑向位置誤差。

    以上脈沖星輪廓畸變的形成機(jī)理是通過(guò)脈沖星輪廓畸變反演橫向位置誤差的理論基礎(chǔ)。這說(shuō)明脈沖星橫向測(cè)距是可行的。

    2 徑向/橫向測(cè)距方法

    本節(jié)研究徑向/橫向測(cè)距方法。相較于傳統(tǒng)的X射線脈沖星導(dǎo)航方法,徑向/橫向測(cè)距方法引入了可觀測(cè)橫向方向這一概念。徑向/橫向測(cè)距方法并非能夠估計(jì)脈沖星橫面上任意方向的位置,而是在可觀測(cè)橫向方向上提供較高精度的位置估計(jì)。因此,該方法也被稱(chēng)為徑向/OTD測(cè)距方法。

    2.1 迭代循環(huán)互相關(guān)法

    迭代循環(huán)互相關(guān)法是在循環(huán)互相關(guān)[22]的基礎(chǔ)上,利用位置估計(jì)值修正初始位置值。循環(huán)互相關(guān)法適合于勻速直線飛行條件下的脈沖星TOA估計(jì);而迭代循環(huán)互相關(guān)法適合于非線性軌道。這是因?yàn)樵诜蔷€性軌道下,初始徑向位置誤差會(huì)引起徑向位置誤差增量。通過(guò)多次迭代,消除徑向位置誤差自干擾。迭代循環(huán)互相關(guān)法的流程圖如圖2所示。具體的算法流程如下:

    圖2 迭代循環(huán)互相關(guān)法

    (16)

    (17)

    (18)

    2.2 面向OTD2法

    航天器的位置誤差可以分解為三個(gè)相互垂直的分量。位置誤差在脈沖星徑向上的投影是徑向位置誤差;位置誤差在脈沖星橫面上的投影也可再次分解為兩個(gè)垂直方向。根據(jù)第1節(jié)的位置誤差傳播性質(zhì)與脈沖星輪廓畸變?cè)?脈沖星橫面上的位置誤差引起徑向位置誤差,進(jìn)而導(dǎo)致脈沖星累積輪廓的畸變。經(jīng)研究發(fā)現(xiàn),橫面內(nèi)不同方向上誤差引起的脈沖星輪廓畸變程度不同。引起的脈沖星輪廓畸變?cè)酱?該橫向方向的可觀測(cè)度越好。主要目的是在脈沖星橫面內(nèi),確定引起脈沖星輪廓畸變最大的可觀測(cè)橫向方向。

    圖3 面向OTD2法

    (19)

    式中:δr表示采樣間隔;N表示單個(gè)方向上的采樣數(shù);k和l為變量,k=1,2,…,N,l=1,2,…,N。這樣可得N2個(gè)采樣值。

    (20)

    圖4 脈沖星輪廓畸變度模型

    徑向/OTD測(cè)距法僅在山脊垂直方向估計(jì),而非在山脊方向和山脊垂直方向同時(shí)開(kāi)展二維估計(jì)。究其原因,一是在山脊方向估計(jì),定位精度必然不高;二是脈沖星橫面上的二維搜索必然導(dǎo)致計(jì)算量增大,影響X射線脈沖星導(dǎo)航系統(tǒng)的實(shí)時(shí)性。綜上所述,本文僅在OTD上開(kāi)展測(cè)距。

    5) 可觀測(cè)橫向方向估計(jì)。構(gòu)建脈沖星輪廓畸變度模型后,根據(jù)脊線方向可計(jì)算出可觀測(cè)橫向方向。通過(guò)函數(shù)擬合該山脊模型,擬合函數(shù)如下:

    z=a0-b0|px-y+d|-e(py+x+c0)2

    (21)

    式中:a0為最大值;b0為當(dāng)前值與最大值之間的斜率;p為脊線斜率;d和e分別為脊線的截距和曲率;c0為常數(shù);x和y為n1和n2方向上的變量。通過(guò)擬合,即可估計(jì)出p。

    可觀測(cè)橫向方向nOTD可通過(guò)下式計(jì)算:

    nOTD=nr×(n1+p·n2)

    (22)

    式中:×表示矢量叉乘;n1+p·n2表示脊線方向矢量,其與脈沖星徑向矢量nr的叉乘可得到脊線垂直方向矢量。

    (23)

    式中:δr表示采樣間隔;N表示采樣數(shù);k為變量,k=1,2,…,N。

    (24)

    (25)

    (26)

    通過(guò)雙一維搜索,分別求出了徑向和OTD上的航天器位置矢量。與二維搜索相比,采樣數(shù)大幅下降,計(jì)算量也必將大幅下降。通過(guò)先估計(jì)脈沖星徑向位置,消除了徑向位置誤差的干擾,必然能有效提高OTD上的定位精度。

    此外,通過(guò)后續(xù)3.2節(jié)的仿真實(shí)驗(yàn)可以看出,OTD與初始位置矢量有關(guān)??筛鶕?jù)初始位置快速確定OTD,從而避免了本節(jié)關(guān)于OTD的復(fù)雜計(jì)算,進(jìn)一步減小了計(jì)算量。

    利用三顆脈沖星的TOA進(jìn)行解算,可得到航天器的三維位置。徑向/OTD測(cè)距法可以通過(guò)一顆脈沖星的輻射信號(hào)估計(jì)出OTD和徑向上的位置,再結(jié)合另一顆脈沖星的TOA實(shí)現(xiàn)三維定位。位置協(xié)方差矩陣C表示為

    [(HTH)-1HT]T

    (27)

    3 仿真實(shí)驗(yàn)與結(jié)果分析

    本節(jié)考察僅利用Crab脈沖星的徑向/OTD測(cè)距方法的準(zhǔn)確性和實(shí)時(shí)性。下面將研究初始位置、累積時(shí)間、軌道半長(zhǎng)軸、軌道傾角等因素對(duì)徑向/橫向測(cè)距方法的影響,并對(duì)比了雙一維搜索法與二維搜索法。為了體現(xiàn)徑向/OTD測(cè)距法的優(yōu)越性,將采用三顆脈沖星的傳統(tǒng)脈沖星測(cè)距法與采用兩顆脈沖星的徑向/OTD測(cè)距法作了對(duì)比。

    3.1 仿真條件

    導(dǎo)航脈沖星是Crab脈沖星PSR B0531+21,其基本參數(shù)如表1所示。研究發(fā)現(xiàn),脈沖星累積時(shí)間與軌道周期均對(duì)OTD位置誤差有影響。將脈沖星累積時(shí)間與軌道周期之比(RPATOP)作為徑向/OTD測(cè)距法的重要參數(shù)。地心坐標(biāo)系(J2000.0)下的地球衛(wèi)星軌道六要素如表2所示。仿真實(shí)驗(yàn)平臺(tái)是處理器為英特爾Core i5-8300H@2.30 GHz,內(nèi)存為8G的聯(lián)想筆記本。

    表1 Crab脈沖星及X射線敏感器相關(guān)參數(shù)

    表2 軌道六要素

    3.2 初始位置與累積時(shí)間

    為了探索脈沖星觀測(cè)初始位置與估計(jì)精度之間的關(guān)系,本文用角表示軌道面位置,脈沖星徑向矢量在軌道面上的投影為0°。圖5給出了初始位置角及RPATOP對(duì)位置估計(jì)誤差的影響。其中,圖5中的初始位置角是指初始位置矢量與脈沖星徑向矢量在軌道面上的投影之間的夾角。軌道周期約為5 215 s,RPATOP為1/2時(shí),脈沖星累積時(shí)間約為2 607 s。

    圖5 初始位置角與RPATOP對(duì)位置估計(jì)誤差的影響

    從圖5中可以看出,在不同RPATOP下,徑向和OTD上的位置誤差呈現(xiàn)兩個(gè)波峰。這表明位置誤差是呈周期變化的,且變化周期約為半個(gè)軌道周期。隨著脈沖星累積時(shí)間的延長(zhǎng),徑向和OTD上的位置誤差均減小。通過(guò)以上現(xiàn)象可以得出,當(dāng)RPATOP為1/2,且初始位置角為90°或270°時(shí),徑向/橫向測(cè)距方法達(dá)到最優(yōu)。原因有二:1) 精度高。此時(shí)的徑向和OTD上的位置誤差分別約為30 m和300 m,低于RPATOP為1/3時(shí)的位置誤差。2)穩(wěn)定。初始位置角為90°,累積完成后,下次的初始位置為270°,二者的精度相當(dāng)。而RPATOP為2/3時(shí),相鄰的兩次位置估計(jì)誤差相差較大,不利于導(dǎo)航Kalman濾波器的設(shè)計(jì)。后續(xù)實(shí)驗(yàn)中,本文設(shè)置RPATOP為1/2。90°或270°時(shí)的初始位置定義為垂直初始位置。3.3節(jié)、3.4節(jié)、3.5節(jié)的仿真實(shí)驗(yàn)均是在垂直初始位置上開(kāi)展的。

    圖6給出了RPATOP為1/2時(shí),初始位置角與OTD角之間的關(guān)系。其中,OTD角是指OTD與初始位置矢量投影之間的夾角。從圖6可以看出,OTD相對(duì)于初始位置矢量投影的位置角在小區(qū)間[-29°,-20°]內(nèi),且呈周期性變化。這一性質(zhì)為OTD的快速確定提供了依據(jù)。導(dǎo)航系統(tǒng)在軌工作時(shí),將初始位置矢量投影按照地面預(yù)先估算的結(jié)果,旋轉(zhuǎn)特定度數(shù)即可得到OTD。在RPATOP為1/3和2/3時(shí)也符合類(lèi)似規(guī)律,在此不再贅述。這樣避免了2.2節(jié)OTD確定過(guò)程的復(fù)雜計(jì)算。

    圖6 初始位置角與OTD角之間的關(guān)系

    此外,結(jié)合圖5可知,即使OTD方位存在誤差,其對(duì)應(yīng)的位置精度并不會(huì)大幅下降。即使快速確定OTD時(shí)存在一定誤差,也不會(huì)造成導(dǎo)航性能下降。因此,這是一種完全可行的方案。

    3.3 半長(zhǎng)軸

    圖7給出了軌道半長(zhǎng)軸對(duì)徑向和OTD位置估計(jì)誤差的影響。從圖7中可以看出,脈沖星累積時(shí)間隨軌道半長(zhǎng)軸延長(zhǎng),徑向和OTD位置估計(jì)誤差隨軌道半長(zhǎng)軸減小。

    圖7 半長(zhǎng)軸對(duì)位置誤差的影響

    理論上,脈沖星位置估計(jì)誤差應(yīng)與累積時(shí)間的開(kāi)方成反比。但是,OTD位置估計(jì)誤差僅略微下降。究其原因,軌道半長(zhǎng)軸越長(zhǎng),軌道非線性程度越弱。弱非線性軌道引起OTD上的位置估計(jì)誤差增大。長(zhǎng)累積時(shí)間和弱非線性軌道綜合作用,使得OTD位置估計(jì)誤差呈現(xiàn)圖7所示的結(jié)果。

    徑向位置估計(jì)誤差受非線性軌道的影響小,與累積時(shí)間的開(kāi)方近似成反比。

    3.4 X射線敏感器

    X射線敏感器的性能與徑向/OTD測(cè)距法息息相關(guān)。下面,考察X射線敏感器的時(shí)間分辨率和有效面積對(duì)徑向/OTD測(cè)距法性能的影響,如表3所示。從表3可以看出,隨著時(shí)間分辨率的下降和有效面積的減小,徑向/OTD測(cè)距法的精度均下降。此外,時(shí)間分辨率的下降有助于縮短計(jì)算時(shí)間。究其原因,時(shí)間分辨率的下降意味著脈沖星信號(hào)數(shù)據(jù)的減少。時(shí)間分辨率的提高意味著需提升傳感器工藝;而有效面積的增加則意味著載重的增加。在實(shí)際工程中,需綜合考慮以上因素,選擇合適的X射線敏感器指標(biāo)參數(shù)。

    表3 X射線敏感器性能

    3.5 雙一維搜索與二維搜索

    表4 雙一維搜索與二維搜索

    3.6 傳統(tǒng)脈沖星測(cè)距與徑向/OTD測(cè)距法

    若要實(shí)現(xiàn)三維定位,傳統(tǒng)脈沖星測(cè)距方法需要至少三顆脈沖星(RSP B0531+21、B1821-24和B1937+21)。而徑向/OTD測(cè)距法只需兩顆脈沖星(RSP B0531+21、B1937+21)。傳統(tǒng)脈沖星測(cè)距方法利用循環(huán)互相關(guān)法估計(jì)脈沖到達(dá)時(shí)間;徑向/OTD測(cè)距法對(duì)RSP B0531+21在兩個(gè)方向上估計(jì)位置信息,對(duì)B1937+21采用循環(huán)互相關(guān)法。估計(jì)完成后,直接利用式(27)求出協(xié)方差矩陣C。

    根據(jù)C可得脈沖星三維定位誤差,如表5所示。傳統(tǒng)脈沖星測(cè)距方法需三顆脈沖星,徑向/OTD測(cè)距法僅需兩顆脈沖星;在不同X射線敏感器面積和時(shí)間分辨率下,徑向/OTD測(cè)距法的定位精度均高于傳統(tǒng)脈沖星測(cè)距方法。究其原因,弱脈沖星流量很低,導(dǎo)致測(cè)距精度低。其測(cè)距精度甚至低于Crab脈沖星OTD上的定位精度。

    表5 兩種脈沖星測(cè)距方法的比較

    若采用傳統(tǒng)脈沖星測(cè)距法,航天器需安裝三個(gè)X射線敏感器分別對(duì)準(zhǔn)三顆脈沖星,或者用一個(gè)X射線敏感器輪流觀測(cè)三顆脈沖星。無(wú)論哪種方案都會(huì)增加導(dǎo)航系統(tǒng)的負(fù)擔(dān)。徑向/橫向測(cè)距方法只需安裝兩個(gè)X射線敏感器分別對(duì)準(zhǔn)兩顆脈沖星,或者用一個(gè)X射線敏感器輪流觀測(cè)兩顆脈沖星。這無(wú)疑大大減小了導(dǎo)航系統(tǒng)的負(fù)擔(dān)。綜上所述,徑向/OTD測(cè)距法優(yōu)于傳統(tǒng)脈沖星測(cè)距法。

    4 結(jié) 論

    傳統(tǒng)脈沖星測(cè)距方法僅能在脈沖星徑向上定位。針對(duì)這一問(wèn)題,本文利用脈沖星累積輪廓畸變,開(kāi)展OTD上的定位,并提出了僅使用Crab脈沖星的徑向/橫向測(cè)距方法。該方法分別使用循環(huán)互相關(guān)法和卡方法估計(jì)脈沖星徑向和OTD上的位置信息。實(shí)驗(yàn)結(jié)果表明,徑向/橫向測(cè)距法僅利用Crab脈沖星就能在兩個(gè)方向上實(shí)現(xiàn)定位,其徑向和OTD上的定位精度分別約為30 m和300 m,計(jì)算時(shí)間僅需1分鐘。與傳統(tǒng)X射線脈沖星導(dǎo)航相比,基于徑向/橫向測(cè)距的脈沖星導(dǎo)航所需脈沖星數(shù)量少,且定位精度高。綜上,該方法具有實(shí)時(shí)、高精度、強(qiáng)魯棒等優(yōu)點(diǎn)。若將徑向/橫向測(cè)距法應(yīng)用于X射線脈沖星導(dǎo)航領(lǐng)域,可在降低系統(tǒng)負(fù)擔(dān)的同時(shí)大幅提高定位精度,為未來(lái)X射線脈沖星導(dǎo)航系統(tǒng)設(shè)計(jì)提供新思路和理論參考。

    猜你喜歡
    脈沖星畸變測(cè)距
    “中國(guó)天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    發(fā)現(xiàn)脈沖星的女天文學(xué)家——貝爾
    科學(xué)(2022年4期)2022-10-25 02:43:42
    類(lèi)星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    在Lightroom中校正鏡頭與透視畸變
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    基于虛擬觀測(cè)值的X射線單脈沖星星光組合導(dǎo)航
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    相對(duì)差分單項(xiàng)測(cè)距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    輻射誘導(dǎo)染色體畸變的快速FISH方法的建立
    《癌變·畸變·突變》2014年第26卷索引
    91九色精品人成在线观看| 成人国语在线视频| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 久久久久久久大尺度免费视频| 人妻 亚洲 视频| 香蕉丝袜av| 亚洲人成电影免费在线| 可以免费在线观看a视频的电影网站| 国产激情久久老熟女| 精品第一国产精品| 好男人视频免费观看在线| 亚洲精品在线美女| 男女床上黄色一级片免费看| 一区在线观看完整版| 中国国产av一级| 亚洲精品一卡2卡三卡4卡5卡 | 欧美精品亚洲一区二区| 80岁老熟妇乱子伦牲交| 国产黄色免费在线视频| 熟女少妇亚洲综合色aaa.| av电影中文网址| av线在线观看网站| 久久女婷五月综合色啪小说| 日韩一卡2卡3卡4卡2021年| 啦啦啦视频在线资源免费观看| 久热爱精品视频在线9| 青春草视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 免费在线观看黄色视频的| 日本av手机在线免费观看| 国产精品一区二区精品视频观看| 激情视频va一区二区三区| 丰满少妇做爰视频| 国产一区二区三区av在线| 又粗又硬又长又爽又黄的视频| 国产成人精品久久二区二区免费| 亚洲精品第二区| 女人被躁到高潮嗷嗷叫费观| 男人添女人高潮全过程视频| 日韩熟女老妇一区二区性免费视频| 国产成人精品久久久久久| 免费在线观看黄色视频的| 久久ye,这里只有精品| 热re99久久国产66热| 一本久久精品| 女人被躁到高潮嗷嗷叫费观| 久久毛片免费看一区二区三区| 亚洲成色77777| 老鸭窝网址在线观看| 69精品国产乱码久久久| 日本vs欧美在线观看视频| xxx大片免费视频| 9色porny在线观看| 99九九在线精品视频| 男女免费视频国产| 亚洲成色77777| 十八禁人妻一区二区| 国产成人欧美在线观看 | 亚洲人成电影免费在线| 黄色片一级片一级黄色片| 国产成人免费观看mmmm| 色婷婷av一区二区三区视频| 首页视频小说图片口味搜索 | 啦啦啦啦在线视频资源| h视频一区二区三区| 高清视频免费观看一区二区| 99精国产麻豆久久婷婷| 国产一区亚洲一区在线观看| 亚洲欧美一区二区三区国产| 宅男免费午夜| 婷婷色综合大香蕉| 热re99久久国产66热| 亚洲av在线观看美女高潮| 波野结衣二区三区在线| 精品一区二区三区四区五区乱码 | 人人妻人人爽人人添夜夜欢视频| 久久九九热精品免费| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲高清精品| 叶爱在线成人免费视频播放| 午夜久久久在线观看| 日韩 亚洲 欧美在线| 亚洲人成网站在线观看播放| 国产黄色免费在线视频| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 亚洲av日韩精品久久久久久密 | 亚洲av综合色区一区| 精品少妇黑人巨大在线播放| 免费不卡黄色视频| 777久久人妻少妇嫩草av网站| 久久久精品免费免费高清| 1024视频免费在线观看| 久久99热这里只频精品6学生| 亚洲国产成人一精品久久久| 国产熟女欧美一区二区| 欧美精品av麻豆av| 午夜精品国产一区二区电影| 成年美女黄网站色视频大全免费| 国产片内射在线| 激情视频va一区二区三区| 捣出白浆h1v1| 成人国语在线视频| 秋霞在线观看毛片| av有码第一页| 精品欧美一区二区三区在线| 水蜜桃什么品种好| 黄色视频不卡| 婷婷丁香在线五月| 亚洲精品在线美女| 国产精品.久久久| 国产一区二区在线观看av| 午夜福利影视在线免费观看| 久久天堂一区二区三区四区| 国产成人av教育| 日日夜夜操网爽| 亚洲欧美日韩另类电影网站| 看十八女毛片水多多多| 赤兔流量卡办理| 黑人巨大精品欧美一区二区蜜桃| 男女午夜视频在线观看| av视频免费观看在线观看| 国产免费现黄频在线看| 婷婷色麻豆天堂久久| 成年美女黄网站色视频大全免费| 又大又爽又粗| 美女脱内裤让男人舔精品视频| 十八禁人妻一区二区| 免费不卡黄色视频| 啦啦啦中文免费视频观看日本| av不卡在线播放| 午夜免费成人在线视频| 欧美激情高清一区二区三区| √禁漫天堂资源中文www| 欧美成人午夜精品| 叶爱在线成人免费视频播放| 亚洲国产精品一区三区| 欧美日韩亚洲高清精品| 美女大奶头黄色视频| 久久午夜综合久久蜜桃| 精品人妻一区二区三区麻豆| 这个男人来自地球电影免费观看| 精品高清国产在线一区| tube8黄色片| 久久免费观看电影| 国产免费一区二区三区四区乱码| avwww免费| 中文欧美无线码| 亚洲中文日韩欧美视频| 久久毛片免费看一区二区三区| 麻豆国产av国片精品| 国产免费一区二区三区四区乱码| kizo精华| 成人国产av品久久久| 视频在线观看一区二区三区| 国精品久久久久久国模美| 亚洲av成人不卡在线观看播放网 | 99久久精品国产亚洲精品| 午夜福利影视在线免费观看| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 超碰成人久久| 欧美黑人欧美精品刺激| 午夜激情久久久久久久| 午夜精品国产一区二区电影| cao死你这个sao货| 99国产精品99久久久久| 亚洲人成网站在线观看播放| 国产成人一区二区三区免费视频网站 | 丁香六月欧美| 人妻 亚洲 视频| 观看av在线不卡| 女人爽到高潮嗷嗷叫在线视频| 久久精品aⅴ一区二区三区四区| 中文欧美无线码| 亚洲欧美一区二区三区久久| 人妻 亚洲 视频| 精品国产乱码久久久久久小说| 国产精品久久久久成人av| 国产不卡av网站在线观看| 亚洲人成网站在线观看播放| 成年人黄色毛片网站| 精品国产国语对白av| 欧美变态另类bdsm刘玥| 99久久精品国产亚洲精品| 久久国产精品影院| 成人免费观看视频高清| 国产精品免费大片| 狠狠精品人妻久久久久久综合| 久久精品国产亚洲av高清一级| 黄片播放在线免费| 久久精品aⅴ一区二区三区四区| 最黄视频免费看| 99国产精品一区二区三区| 18禁国产床啪视频网站| 久久久久视频综合| 18禁观看日本| 国产精品偷伦视频观看了| 国产免费一区二区三区四区乱码| 欧美黑人欧美精品刺激| 欧美国产精品一级二级三级| 亚洲国产成人一精品久久久| 亚洲黑人精品在线| 咕卡用的链子| 久久av网站| 国产成人欧美在线观看 | 亚洲人成电影免费在线| 欧美成人午夜精品| 永久免费av网站大全| 高潮久久久久久久久久久不卡| 在线观看一区二区三区激情| 国产国语露脸激情在线看| 亚洲午夜精品一区,二区,三区| 深夜精品福利| 久久久久视频综合| 99久久人妻综合| 又大又爽又粗| 亚洲午夜精品一区,二区,三区| 久久久久精品人妻al黑| 又粗又硬又长又爽又黄的视频| 丝袜在线中文字幕| 中国美女看黄片| 国产成人精品久久二区二区免费| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 黄色 视频免费看| 纵有疾风起免费观看全集完整版| 一区二区三区乱码不卡18| 亚洲第一青青草原| 最新在线观看一区二区三区 | 老鸭窝网址在线观看| 精品少妇内射三级| 亚洲精品美女久久久久99蜜臀 | 亚洲,欧美,日韩| 亚洲精品乱久久久久久| 亚洲国产精品国产精品| 婷婷色综合大香蕉| 亚洲精品在线美女| 国产黄色免费在线视频| 一本综合久久免费| videosex国产| 黄色a级毛片大全视频| 亚洲第一av免费看| 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 美女大奶头黄色视频| 99久久99久久久精品蜜桃| 国产不卡av网站在线观看| 十分钟在线观看高清视频www| 国产精品国产三级国产专区5o| 中文字幕人妻丝袜制服| 丝袜美足系列| 成年人免费黄色播放视频| 国产亚洲午夜精品一区二区久久| 在线观看人妻少妇| 久久久精品国产亚洲av高清涩受| 亚洲国产欧美日韩在线播放| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 一个人免费看片子| 国产欧美亚洲国产| 久久免费观看电影| 国产精品免费大片| 久久精品亚洲熟妇少妇任你| 日韩制服丝袜自拍偷拍| 国产亚洲欧美在线一区二区| 精品熟女少妇八av免费久了| 久久影院123| 丰满饥渴人妻一区二区三| 午夜福利一区二区在线看| 男女边摸边吃奶| 女警被强在线播放| 日日夜夜操网爽| 亚洲精品成人av观看孕妇| 国产亚洲一区二区精品| 亚洲精品美女久久久久99蜜臀 | 亚洲中文字幕日韩| 大香蕉久久网| 精品熟女少妇八av免费久了| 亚洲精品中文字幕在线视频| 一区二区三区精品91| 欧美精品啪啪一区二区三区 | 电影成人av| 无限看片的www在线观看| 老汉色av国产亚洲站长工具| 国产一区二区 视频在线| 在线观看免费视频网站a站| 老司机午夜十八禁免费视频| 自线自在国产av| 两性夫妻黄色片| 日韩 亚洲 欧美在线| 成年美女黄网站色视频大全免费| 夜夜骑夜夜射夜夜干| 黄色一级大片看看| 亚洲 欧美一区二区三区| 另类亚洲欧美激情| 多毛熟女@视频| 日本欧美视频一区| 国产成人欧美在线观看 | 国产成人啪精品午夜网站| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 好男人电影高清在线观看| 国产精品成人在线| 亚洲三区欧美一区| 嫁个100分男人电影在线观看 | 一区在线观看完整版| 一级黄片播放器| 高清黄色对白视频在线免费看| 女人爽到高潮嗷嗷叫在线视频| av有码第一页| 午夜福利视频精品| 国产在线免费精品| 国产免费视频播放在线视频| www.av在线官网国产| 欧美久久黑人一区二区| 黄片小视频在线播放| 男女高潮啪啪啪动态图| 赤兔流量卡办理| 高清av免费在线| av在线老鸭窝| 国产一区亚洲一区在线观看| 日本欧美视频一区| 高清av免费在线| 国产精品久久久av美女十八| 爱豆传媒免费全集在线观看| 在线观看一区二区三区激情| 777米奇影视久久| 女人高潮潮喷娇喘18禁视频| 亚洲免费av在线视频| 中文字幕制服av| 91国产中文字幕| 久热爱精品视频在线9| 在线看a的网站| 人妻人人澡人人爽人人| 中文字幕人妻丝袜制服| 热99国产精品久久久久久7| 精品国产一区二区三区久久久樱花| 日日爽夜夜爽网站| 久久ye,这里只有精品| 国产一区有黄有色的免费视频| 菩萨蛮人人尽说江南好唐韦庄| 18禁观看日本| 国产一区二区在线观看av| 成人国产一区最新在线观看 | 国产精品麻豆人妻色哟哟久久| 国产成人91sexporn| 好男人电影高清在线观看| 精品国产乱码久久久久久小说| 每晚都被弄得嗷嗷叫到高潮| 美女主播在线视频| www.自偷自拍.com| 亚洲av男天堂| 一个人免费看片子| 午夜激情久久久久久久| 久久久精品免费免费高清| 欧美中文综合在线视频| 免费看不卡的av| 激情五月婷婷亚洲| 亚洲,欧美精品.| 制服人妻中文乱码| 十八禁网站网址无遮挡| 国产1区2区3区精品| 国产欧美亚洲国产| 亚洲国产欧美网| 看十八女毛片水多多多| 国产成人a∨麻豆精品| videosex国产| 一本久久精品| 最近最新中文字幕大全免费视频 | 日本av免费视频播放| 亚洲精品在线美女| 在线观看免费午夜福利视频| 韩国精品一区二区三区| 成人国产av品久久久| 亚洲自偷自拍图片 自拍| 国产精品av久久久久免费| 一级毛片我不卡| 人体艺术视频欧美日本| 欧美人与性动交α欧美软件| 精品福利永久在线观看| 咕卡用的链子| 亚洲久久久国产精品| 悠悠久久av| 国产一区二区三区综合在线观看| 91麻豆av在线| 精品少妇一区二区三区视频日本电影| 天天添夜夜摸| 真人做人爱边吃奶动态| 欧美97在线视频| 最新的欧美精品一区二区| 亚洲情色 制服丝袜| 亚洲国产成人一精品久久久| 精品国产一区二区三区久久久樱花| av天堂久久9| 无遮挡黄片免费观看| 国产色视频综合| 男人爽女人下面视频在线观看| www.av在线官网国产| 欧美大码av| 国产真人三级小视频在线观看| 91国产中文字幕| 国产老妇伦熟女老妇高清| av有码第一页| 久久国产精品人妻蜜桃| 日韩视频在线欧美| 又大又黄又爽视频免费| 国产日韩欧美视频二区| 亚洲欧美精品自产自拍| 午夜福利视频精品| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 黑人巨大精品欧美一区二区蜜桃| 91国产中文字幕| 国产精品一区二区精品视频观看| 19禁男女啪啪无遮挡网站| 99国产精品99久久久久| 黄频高清免费视频| 黄色视频在线播放观看不卡| 欧美亚洲 丝袜 人妻 在线| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 精品亚洲乱码少妇综合久久| 中文精品一卡2卡3卡4更新| 亚洲精品久久成人aⅴ小说| 黑人猛操日本美女一级片| 熟女av电影| 午夜免费鲁丝| 国产免费现黄频在线看| 免费日韩欧美在线观看| 亚洲精品日韩在线中文字幕| 国产一级毛片在线| e午夜精品久久久久久久| 国产真人三级小视频在线观看| 欧美少妇被猛烈插入视频| 国产成人一区二区在线| 91精品伊人久久大香线蕉| 香蕉丝袜av| 精品福利观看| 如日韩欧美国产精品一区二区三区| 欧美成人精品欧美一级黄| 日韩一卡2卡3卡4卡2021年| 国产又爽黄色视频| 伦理电影免费视频| 九草在线视频观看| 夜夜骑夜夜射夜夜干| 电影成人av| 超碰97精品在线观看| 视频区图区小说| 91成人精品电影| 国产精品三级大全| 中文字幕人妻丝袜一区二区| 久热这里只有精品99| 国产男人的电影天堂91| 18禁国产床啪视频网站| 久久青草综合色| 最近中文字幕2019免费版| 手机成人av网站| 亚洲五月婷婷丁香| 人人妻,人人澡人人爽秒播 | a 毛片基地| 国产色视频综合| 亚洲图色成人| 大片电影免费在线观看免费| 国产伦理片在线播放av一区| 色播在线永久视频| 交换朋友夫妻互换小说| 蜜桃国产av成人99| 国产97色在线日韩免费| 热99久久久久精品小说推荐| 在线观看免费高清a一片| 国产精品.久久久| 不卡av一区二区三区| 日韩大码丰满熟妇| 亚洲av在线观看美女高潮| 一级片'在线观看视频| 午夜福利,免费看| 男男h啪啪无遮挡| 国产免费福利视频在线观看| 黄网站色视频无遮挡免费观看| 国产国语露脸激情在线看| 国产日韩欧美在线精品| 久久影院123| 午夜福利免费观看在线| 亚洲国产精品成人久久小说| 日本午夜av视频| 午夜福利视频在线观看免费| 在线观看国产h片| 看十八女毛片水多多多| 丝袜脚勾引网站| 日本色播在线视频| 搡老岳熟女国产| 成人影院久久| 性高湖久久久久久久久免费观看| 国产片特级美女逼逼视频| 欧美精品一区二区免费开放| 久久久国产欧美日韩av| 在线观看国产h片| 波野结衣二区三区在线| 黄色视频在线播放观看不卡| 日韩伦理黄色片| 国产精品麻豆人妻色哟哟久久| 制服人妻中文乱码| 最近中文字幕2019免费版| 久热爱精品视频在线9| 18禁黄网站禁片午夜丰满| 久久狼人影院| 亚洲色图综合在线观看| 亚洲av日韩在线播放| 亚洲中文字幕日韩| 十分钟在线观看高清视频www| 久久天堂一区二区三区四区| 精品国产国语对白av| 男女边摸边吃奶| 亚洲 国产 在线| 一边亲一边摸免费视频| 国产精品九九99| 中文乱码字字幕精品一区二区三区| 午夜av观看不卡| 久久久国产欧美日韩av| 黄片播放在线免费| 久久人妻熟女aⅴ| 妹子高潮喷水视频| 欧美黄色淫秽网站| 高清视频免费观看一区二区| 2021少妇久久久久久久久久久| 成人手机av| 午夜福利视频在线观看免费| 天堂8中文在线网| 免费观看人在逋| 午夜两性在线视频| 七月丁香在线播放| 亚洲av美国av| 一本大道久久a久久精品| 亚洲精品久久久久久婷婷小说| 久久久久久久精品精品| 51午夜福利影视在线观看| 亚洲精品国产区一区二| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品美女久久久久99蜜臀 | 青草久久国产| 精品亚洲成a人片在线观看| 制服人妻中文乱码| 一级毛片黄色毛片免费观看视频| 久久精品久久久久久久性| 亚洲精品乱久久久久久| 国产欧美日韩精品亚洲av| 中文乱码字字幕精品一区二区三区| 久9热在线精品视频| 男女无遮挡免费网站观看| 黄色毛片三级朝国网站| 一级,二级,三级黄色视频| 亚洲伊人色综图| 国产高清国产精品国产三级| 色网站视频免费| 国产亚洲一区二区精品| av又黄又爽大尺度在线免费看| 中文字幕人妻丝袜一区二区| 亚洲精品日本国产第一区| 欧美日韩一级在线毛片| 午夜福利乱码中文字幕| 男男h啪啪无遮挡| 亚洲av国产av综合av卡| 精品人妻一区二区三区麻豆| 久久综合国产亚洲精品| 久久精品国产亚洲av高清一级| 少妇粗大呻吟视频| 国产亚洲午夜精品一区二区久久| 夫妻性生交免费视频一级片| 丰满饥渴人妻一区二区三| 久久久久国产一级毛片高清牌| 亚洲精品美女久久久久99蜜臀 | 别揉我奶头~嗯~啊~动态视频 | 婷婷色麻豆天堂久久| 国产亚洲av片在线观看秒播厂| 我的亚洲天堂| 免费观看av网站的网址| 啦啦啦啦在线视频资源| 一级片'在线观看视频| videosex国产| 亚洲成人免费电影在线观看 | 国产成人一区二区在线| 男人爽女人下面视频在线观看| 少妇人妻 视频| 热99久久久久精品小说推荐| 高清欧美精品videossex| 后天国语完整版免费观看| 久久精品熟女亚洲av麻豆精品| 亚洲中文日韩欧美视频| 久久精品亚洲av国产电影网| 欧美精品高潮呻吟av久久| 亚洲精品av麻豆狂野| 色精品久久人妻99蜜桃| 一个人免费看片子| 黄色毛片三级朝国网站| 亚洲精品久久久久久婷婷小说| 午夜影院在线不卡| 男男h啪啪无遮挡| 亚洲五月婷婷丁香| 青春草视频在线免费观看| 99久久99久久久精品蜜桃| 日本av免费视频播放| 日韩欧美一区视频在线观看| 久久精品国产综合久久久| 男人操女人黄网站| 午夜福利一区二区在线看| 又黄又粗又硬又大视频| 免费在线观看黄色视频的| 宅男免费午夜| 美女国产高潮福利片在线看| 国产成人91sexporn| 日本黄色日本黄色录像| 亚洲五月色婷婷综合| 成年女人毛片免费观看观看9 |