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

    X射線脈沖星TOA數(shù)量對(duì)計(jì)時(shí)精度和導(dǎo)航的影響分析

    2021-04-26 12:35:00韓孟納童明雷朱鴻旭趙成仕朱幸芝
    時(shí)間頻率學(xué)報(bào) 2021年1期
    關(guān)鍵詞:脈沖星計(jì)時(shí)航天器

    韓孟納,童明雷,朱鴻旭,趙成仕,朱幸芝

    X射線脈沖星TOA數(shù)量對(duì)計(jì)時(shí)精度和導(dǎo)航的影響分析

    韓孟納1,2,3,童明雷1,2,朱鴻旭4,趙成仕1,2,朱幸芝1,2

    (1. 中國科學(xué)院 國家授時(shí)中心,西安 710600;2. 中國科學(xué)院 時(shí)間頻率基準(zhǔn)重點(diǎn)實(shí)驗(yàn)室,西安 710600;3. 中國科學(xué)院大學(xué),北京 100049;4. 山東航天電子技術(shù)研究所,煙臺(tái) 264000)

    針對(duì)XPNAV-1衛(wèi)星發(fā)布的Crab脈沖星首批數(shù)據(jù),采用Taylor頻域相關(guān)法得到了脈沖到達(dá)時(shí)間(TOA)及其對(duì)應(yīng)的測量精度。在此基礎(chǔ)上,去掉光子數(shù)量過少的第5組數(shù)據(jù),得到34組數(shù)據(jù)的擬合前計(jì)時(shí)殘差的均方根(RMS)約27.1 μs,擬合后計(jì)時(shí)殘差RMS約22.5 μs。將光子數(shù)量最多的第24組觀測數(shù)據(jù)按照不同的光子數(shù)量分組,得到了脈沖輪廓信噪比、TOA測量精度及計(jì)時(shí)殘差與光子數(shù)量之間的關(guān)系??紤]空間飛行器實(shí)時(shí)定位需求,本文分析了在保障一定信噪比的前提下增加TOA數(shù)量對(duì)定位精度的影響。將觀測數(shù)據(jù)按20000個(gè)光子重新分組,得到了214個(gè)TOA,其擬合前計(jì)時(shí)殘差RMS約為61.8 μs,擬合后計(jì)時(shí)殘差RMS約61.3 μs。這不但大大提高了導(dǎo)航的實(shí)時(shí)性,而且將定位精度提高了約10%。

    X射線脈沖星;脈沖星計(jì)時(shí);光子個(gè)數(shù);脈沖輪廓

    0 引言

    現(xiàn)代衛(wèi)星導(dǎo)航系統(tǒng)能夠?yàn)榈厍虮砻婧徒乜臻g的各類用戶提供全天候、高精度的自主導(dǎo)航信息服務(wù),但是地球衛(wèi)星導(dǎo)航卻無法對(duì)執(zhí)行深空探測和星際飛行任務(wù)的航天器導(dǎo)航。在深空探測與星際飛行任務(wù)中,航天器距離地球很遠(yuǎn),相對(duì)于地球衛(wèi)星來說,對(duì)深空航天器的位置確定以及姿態(tài)控制難度較大。考慮到在深空探測器跟蹤、導(dǎo)航中,陸基系統(tǒng)由于頻繁的地面操作以及導(dǎo)航誤差隨測控距離的增加而線性增大等問題,利用地面測控系統(tǒng)對(duì)在軌航天器進(jìn)行軌道監(jiān)測與位置修正是不方便的,航天器最好能夠自主導(dǎo)航,而其中一種具有很大發(fā)展前景的航天器自主導(dǎo)航方法就是基于X射線脈沖星計(jì)時(shí)的導(dǎo)航方法[1]。X射線脈沖星導(dǎo)航的自主性強(qiáng),信號(hào)源為宇宙中的自然天體—脈沖星,可以減少航天器對(duì)地面測控系統(tǒng)的依賴,將廣泛應(yīng)用于航天器深空探測和星際飛行任務(wù)[2-4]。

    脈沖星為高速自轉(zhuǎn)的致密天體,其磁極產(chǎn)生強(qiáng)磁場,磁場隨星體的旋轉(zhuǎn)產(chǎn)生感應(yīng)電場,電子被加速離開星體表面,產(chǎn)生輻射光束,如同燈塔的光束[5]。當(dāng)輻射束掃過地球,可被地面射電望遠(yuǎn)鏡接收或空間X射線探測器探測。脈沖星的自轉(zhuǎn)具有高度穩(wěn)定性,尤其是毫秒脈沖星。脈沖星的多方面應(yīng)用,如建立脈沖星時(shí)間尺度[6-8]、探測低頻引力波[9]、計(jì)算太陽系行星質(zhì)量[10]以及脈沖星導(dǎo)航[11-12]等,都是基于它自轉(zhuǎn)的高度穩(wěn)定性的[13]。目前發(fā)現(xiàn)的眾多脈沖星中,有不少具有從射電到X射線的全波段輻射,這將使我們更清楚地認(rèn)識(shí)脈沖星的輻射機(jī)制及其內(nèi)部結(jié)構(gòu)。脈沖星輻射的X射線光子易于被小型化探測器接收,這有利于減少航天器有效載荷的尺寸。脈沖星按照旋轉(zhuǎn)能量來源可分為三類,分別為旋轉(zhuǎn)供能脈沖星、吸積供能脈沖星和磁星。吸積供能脈沖星多位于雙星系統(tǒng),主要通過吸取伴星的物質(zhì)來提供輻射所需的能量,自轉(zhuǎn)極不穩(wěn)定[14]。X射線脈沖星導(dǎo)航所選用的脈沖星須為旋轉(zhuǎn)供能脈沖星,且同時(shí)具有X射線波段與射電波段的輻射。之所以要具有射電波段的輻射,是為了得到脈沖星星歷并建立脈沖星計(jì)時(shí)模型。蟹狀星云(Crab)脈沖星(B0531+21)正符合脈沖星篩選條件,可應(yīng)用于X射線脈沖星導(dǎo)航。Crab脈沖星位于超新星1054 AD中心處,是恒星超新星爆發(fā)后在蟹狀星云中的殘骸,其自轉(zhuǎn)周期約為33 ms,在X射線波段具有很強(qiáng)的流量。但由于Crab脈沖星比較年輕,為得到高精度的計(jì)時(shí)模型,其歷表需要頻繁更新。

    脈沖星導(dǎo)航試驗(yàn)衛(wèi)星(XPNAV-1)是我國首顆X射線脈沖星導(dǎo)航試驗(yàn)衛(wèi)星,該衛(wèi)星由中國航天科技集團(tuán)第五研究院研制,并于2016年11月10日在酒泉衛(wèi)星發(fā)射中心發(fā)射升空[15]。2017年5月9日,北斗衛(wèi)星導(dǎo)航系統(tǒng)網(wǎng)站發(fā)布了35組首批由星載Wolter-I聚焦型探測器探測到的Crab脈沖星數(shù)據(jù),探測器的各項(xiàng)性能參數(shù)如表1所示[16]。XPNAV-1衛(wèi)星主要用于驗(yàn)證星載X射線探測器的性能和其對(duì)空間環(huán)境的適應(yīng)力,并通過積累在軌觀測數(shù)據(jù)完善脈沖星導(dǎo)航算法,為未來航天器深空導(dǎo)航奠定基礎(chǔ)。本文主要采用XPNAV-1衛(wèi)星的發(fā)布數(shù)據(jù)進(jìn)行相應(yīng)的處理與分析。

    表1 Wolter-I聚焦型X射線探測器性能參數(shù)

    1 觀測數(shù)據(jù)處理流程與方法

    1.1 TOA的獲取

    XPNAV-1衛(wèi)星的發(fā)布數(shù)據(jù)給出了衛(wèi)星的軌道與光子事件兩個(gè)文件。其中,軌道文件包括軌道遙測歷元與遙測時(shí)刻衛(wèi)星的狀態(tài)向量兩部分,光子事件文件包括光子TOA(光子到達(dá)探測器的時(shí)刻)以及每個(gè)光子的能量。X射線探測器探測到Crab脈沖星輻射的X射線光子,并由星載原子時(shí)鐘記錄光子到達(dá)探測器的時(shí)刻。光子TOA采用UTC秒的記錄形式,即光子被捕獲時(shí)刻相對(duì)UTC 2008年1月1日0時(shí)0分0秒的秒數(shù)[16]。X射線光子到達(dá)衛(wèi)星的時(shí)刻與衛(wèi)星狀態(tài)向量的遙測時(shí)刻并不一致,因此,需要在光子到達(dá)探測器的時(shí)刻內(nèi)插出航天器的位置與速度,內(nèi)插的方法采用一維線性插值。

    衛(wèi)星在觀測過程中不斷繞地運(yùn)動(dòng),其所處的引力場在不斷變化,星載時(shí)鐘記錄的光子到達(dá)時(shí)刻需轉(zhuǎn)化到慣性系中的某一參考點(diǎn)。太陽系質(zhì)心天球參考系(barycentric celestial reference system,BCRS)是一個(gè)準(zhǔn)慣性系,可用于描述XPNAV-1衛(wèi)星的軌道運(yùn)動(dòng)和脈沖星的觀測信息(如角位置)。我們將Crab脈沖星光子到達(dá)探測器的時(shí)刻轉(zhuǎn)換到BCRS的坐標(biāo)原點(diǎn)——太陽系質(zhì)心(solar system barycentre,SSB)。此過程主要涉及相應(yīng)參考系與時(shí)間系統(tǒng)的轉(zhuǎn)換[16]。X射線波段觀測相對(duì)于射電觀測的一個(gè)優(yōu)點(diǎn)在于:X射線波段觀測無需考慮信號(hào)在星際介質(zhì)中的色散效應(yīng)以及地球大氣傳播時(shí)延。由于發(fā)布的Crab脈沖星星歷未考慮視差項(xiàng),為自洽本文也不予考慮,只進(jìn)行Roemer時(shí)延、Shapiro時(shí)延以及Einstein時(shí)延的改正。

    Roemer時(shí)延項(xiàng)為航天器相對(duì)于SSB的距離在Crab脈沖星視線方向上的投影所產(chǎn)生的時(shí)間延遲量,也將其稱為真空傳播時(shí)延。Shapiro時(shí)延項(xiàng)是由于太陽系中的大質(zhì)量天體引起時(shí)空彎曲,與平直時(shí)空相比這將引起額外的時(shí)間延遲[17-18]。由于地球公轉(zhuǎn)軌道并非嚴(yán)格的圓軌道,且太陽的引力場效應(yīng)不可忽略,在地球上實(shí)現(xiàn)的原子時(shí)從太陽系全局來看并不是均勻的,必須扣除太陽系引力場與地球相對(duì)太陽系質(zhì)心運(yùn)動(dòng)速度的影響,將時(shí)間基準(zhǔn)定義到太陽系質(zhì)心參考系。脈沖星計(jì)時(shí)觀測確定的脈沖星星歷參數(shù)必須明確說明所參考的時(shí)間尺度。BCRS對(duì)應(yīng)的時(shí)間尺度為太陽系質(zhì)心坐標(biāo)時(shí)(TCB),而由于歷史原因,美國NASA噴氣推進(jìn)實(shí)驗(yàn)室(JPL)給出的、用于描述太陽系天體的位置與速度的行星歷表采用的時(shí)間尺度是太陽系質(zhì)心力學(xué)時(shí)(TDB),TDB與TCB之間不存在周期性差異,僅存在線性速率差。由于Crab脈沖星的星歷歸算采用的是TDB,因此須將光子到達(dá)時(shí)刻最終轉(zhuǎn)換為TDB。Einstein時(shí)延項(xiàng)涉及協(xié)調(diào)世界時(shí)(coordinated universal time,UTC)到TDB的時(shí)間系統(tǒng)轉(zhuǎn)換[19]。轉(zhuǎn)換過程需考慮兩個(gè)效應(yīng),即狹義相對(duì)論中的“時(shí)鐘變慢”效應(yīng)(運(yùn)動(dòng)的時(shí)鐘會(huì)變慢)與廣義相對(duì)論中的“引力紅移”效應(yīng)(所處的引力場越強(qiáng),原子鐘的走時(shí)速率越慢)。各種時(shí)間尺度轉(zhuǎn)換公式如下[20-22]:

    在計(jì)算光子到達(dá)SSB處的時(shí)刻時(shí),需要使用太陽系歷表給出太陽與大行星的質(zhì)量以及觀測時(shí)刻各天體到SSB的距離等,一般使用JPL發(fā)布的DE系列歷表。由于Crab脈沖星星歷[23]是使用DE200歷表[24]得到的,為保持一致,本文做光子到達(dá)SSB的時(shí)刻轉(zhuǎn)換時(shí),太陽系歷表依然選用DE200。至此,經(jīng)過一系列的計(jì)算后,得到了光子到達(dá)SSB的時(shí)刻。利用光子到達(dá)SSB的時(shí)刻,在SSB處進(jìn)行脈沖輪廓折疊,分別得到積分脈沖輪廓與標(biāo)準(zhǔn)脈沖輪廓,然后將二者做互相關(guān)處理得到TOA。積分脈沖輪廓是以起始觀測時(shí)刻為起點(diǎn)折疊得到,至于利用多少光子數(shù)據(jù)能得到相對(duì)穩(wěn)定的積分脈沖輪廓,將在下文討論。標(biāo)準(zhǔn)脈沖輪廓需要更長積分時(shí)間的數(shù)據(jù),為提高信噪比,這里采用全部的35組數(shù)據(jù)折疊形成。本文Crab脈沖星的標(biāo)準(zhǔn)脈沖輪廓是以英國Jodrell Bank天文臺(tái)發(fā)布的Crab脈沖星星歷定義的參考?xì)v元為起點(diǎn)折疊得到,星歷定義的參考?xì)v元是射電觀測得到的TOA。由于衛(wèi)星觀測數(shù)據(jù)為X射線波段觀測得到,與脈沖星星歷所采用的射電數(shù)據(jù)屬于不同的波段,因此折疊出來的標(biāo)準(zhǔn)脈沖輪廓存在0.198 ms的零點(diǎn)相位差。標(biāo)準(zhǔn)脈沖輪廓如圖3所示。因脈沖輪廓經(jīng)離散傅里葉變換(DFT)之后在頻域互相關(guān)[25-26],與簡單的進(jìn)行時(shí)域互相關(guān)相比,其得到的相位偏移量更加精確,TOA測量精度更高,所以本文的TOA是采用積分脈沖輪廓與標(biāo)準(zhǔn)脈沖輪廓經(jīng)DFT之后在頻域互相關(guān)的方法得到的。最終,TOA通過起始觀測時(shí)刻加上由相位偏移量轉(zhuǎn)化的時(shí)間偏移量來確定。

    計(jì)時(shí)殘差為測量的TOA與脈沖星計(jì)時(shí)模型在SSB處預(yù)報(bào)的TOA之差,即[5]:

    式(15)中,A為探測器有效面積,T為總的觀測時(shí)間,P為脈沖周期,s為脈沖星輻射流量,b為背景X射線光子到達(dá)速率,參數(shù)與分別表示X射線光子信號(hào)本身的固有噪聲和背景噪聲影響下脈沖形狀對(duì)TOA測量精度的影響。由公式(16)可以看出,若探測器飛行環(huán)境近似一致,用相同的探測器觀測同一顆脈沖星時(shí),TOA不確定度與觀測時(shí)間之間存在關(guān)系。觀測時(shí)間越長,探測器積累的光子數(shù)越多,觀測時(shí)間與光子數(shù)近似為線性關(guān)系。因此,TOA測量的不確定度與光子數(shù)之間近似滿足平方反比關(guān)系。圖2中曲線是TOA不確定度擬合曲線,擬合曲線方程為,其中=5 492.3。擬合方程中的、分別為光子數(shù)與。由圖2可知,光子數(shù)越多,測量的TOA不確定度越小,平方反比關(guān)系成立。

    圖2 光子個(gè)數(shù)與TOA不確定度的關(guān)系

    1.2 觀測數(shù)據(jù)分組方案

    航天器的初始軌道位置由軌道力學(xué)模型給出,航天器在軌運(yùn)動(dòng)時(shí)并不會(huì)嚴(yán)格按照預(yù)定軌道運(yùn)行,真實(shí)位置與估計(jì)位置會(huì)有偏差,需要對(duì)偏差值進(jìn)行實(shí)時(shí)修正。而確定航天器的位置估計(jì)誤差主要是通過將測量的脈沖到達(dá)SSB時(shí)間與脈沖星計(jì)時(shí)模型在SSB處預(yù)報(bào)的TOA的差值來加以改正[12]。由于測量的TOA誤差與積分脈沖輪廓的信噪比成反比,因此使用多長的觀測數(shù)據(jù)可以得到相對(duì)穩(wěn)定的脈沖輪廓是一個(gè)值得研究的問題。雖然較長的積分時(shí)間會(huì)形成高信噪比的脈沖輪廓以及較小的統(tǒng)計(jì)測量誤差,但是如果積分時(shí)間過長,就會(huì)得到較少的積分脈沖輪廓,從而得到較少的TOA,降低了導(dǎo)航效率,無法對(duì)航天器的位置誤差做出實(shí)時(shí)修正。

    下面針對(duì)至少需要多少個(gè)光子可以形成相對(duì)穩(wěn)定的積分脈沖輪廓這一問題來開展研究。由于第22組與24組觀測數(shù)據(jù)的光子數(shù)最多,數(shù)據(jù)量分別為287 446與300 627,先使用這兩組數(shù)據(jù)試驗(yàn)。將兩組數(shù)據(jù)依次按照每組150 000、100 000、75 000、60 000、50 000、37 500、30 000、25 000、20 000、18 750、15 000、10 000個(gè)光子數(shù)分組。假設(shè)觀測過程中探測器所處的空間環(huán)境近似相同,我們定義信噪比(SNR):歸一化的平滑標(biāo)準(zhǔn)脈沖輪廓與噪聲的比值。這里的噪聲是積分脈沖輪廓與平滑標(biāo)準(zhǔn)脈沖輪廓光子數(shù)經(jīng)歸一化后的對(duì)應(yīng)值之差。將標(biāo)準(zhǔn)脈沖輪廓平滑處理是為了更加真實(shí)地還原脈沖信號(hào)。我們采用Savitzky-Golay平滑方法,這種方法通過平滑局部數(shù)據(jù)多項(xiàng)式回歸來實(shí)現(xiàn)曲線的平滑,能夠有效地保留數(shù)據(jù)的原始特征。在平滑過程中,平滑值不能設(shè)置的過大,否則會(huì)導(dǎo)致脈沖輪廓波形失真。平滑之后的標(biāo)準(zhǔn)脈沖輪廓如圖3所示。脈沖輪廓?dú)w一化方法為每個(gè)bin對(duì)應(yīng)的光子數(shù)減去其最小值(具有最少光子bin內(nèi)的光子數(shù))之后除以極差(最大值-最小值),這樣得到的脈沖輪廓位于[0,1]之間。將第22組數(shù)據(jù)與第24組數(shù)據(jù)依次按不同光子數(shù)進(jìn)行分組之后,得到每組光子個(gè)數(shù)與對(duì)應(yīng)積分脈沖輪廓SNR均值的關(guān)系如圖4所示。隨著光子個(gè)數(shù)的增多,積分脈沖輪廓SNR均值有一定的起伏,但整體呈上升趨勢。說明折疊所用的光子數(shù)越多,積分脈沖輪廓的信噪比越高,脈沖輪廓越穩(wěn)定。

    圖3 未平滑與經(jīng)過平滑之后的標(biāo)準(zhǔn)脈沖輪廓

    圖4 光子個(gè)數(shù)與SNR均值的關(guān)系

    表2給出第24組數(shù)據(jù)分組之后每組光子個(gè)數(shù)與其對(duì)應(yīng)的平均TOA不確定度的關(guān)系,而圖5給出了更直觀的展示??梢钥闯觯跀?shù)據(jù)總量相等的情況下,隨著每組光子個(gè)數(shù)的增加,TOA不確定度逐漸下降。由此可知光子個(gè)數(shù)越多,折疊出的脈沖輪廓越穩(wěn)定,與圖2得到的結(jié)論一致。此外,從圖5可以看出,計(jì)時(shí)殘差的RMS值隨著光子個(gè)數(shù)的增加也有減小的趨勢。

    表2 第24組數(shù)據(jù)再次分組后的平均TOA不確定度

    圖5 不同光子數(shù)量對(duì)應(yīng)的TOA不確定度及計(jì)時(shí)殘差特征

    圖6展示了TOA不確定度的變化率隨光子個(gè)數(shù)的增加不斷變小,當(dāng)光子數(shù)足夠多時(shí),TOA不確定度的變化率趨于零。說明隨著積分時(shí)間的不斷增加,脈沖輪廓穩(wěn)定度最終將達(dá)到飽和。在分組折疊過程中,每組光子個(gè)數(shù)達(dá)到20 000個(gè)以上時(shí),脈沖輪廓的波峰才能體現(xiàn)出來。因此,為形成高質(zhì)量積分脈沖輪廓,光子個(gè)數(shù)最好能超過20 000個(gè)。

    圖6 光子個(gè)數(shù)與TOA不確定度變化率的關(guān)系

    2 分組前后計(jì)時(shí)殘差及測量精度比較與分析

    在基于脈沖星計(jì)時(shí)的航天器自主導(dǎo)航中,利用星載X射線探測器接收脈沖星輻射的X射線光子,并將光子數(shù)據(jù)按前文介紹進(jìn)行相應(yīng)處理之后得到脈沖到達(dá)SSB的時(shí)間。如果航天器的位置不準(zhǔn),則將產(chǎn)生額外的計(jì)時(shí)殘差。此時(shí)可以通過測量的脈沖到達(dá)時(shí)間(或相位)信息來給航天器定位,或者說對(duì)航天器做軌道修正以防偏離預(yù)設(shè)軌道太大。XPNAV-1衛(wèi)星發(fā)布的數(shù)據(jù)每分鐘可以獲得800多個(gè)Crab脈沖星X射線光子數(shù)據(jù)。在獲得積分脈沖輪廓時(shí),若增加所需的光子數(shù)據(jù)量,則積分時(shí)間就要增加,同時(shí)數(shù)據(jù)處理過程耗時(shí)也會(huì)增加。若航天器處于磁場與引力場較強(qiáng)的極端飛行環(huán)境,在觀測時(shí)間以及數(shù)據(jù)處理時(shí)間內(nèi),航天器的實(shí)際飛行軌道會(huì)嚴(yán)重偏離預(yù)定軌道。由于航天器的位置估計(jì)誤差會(huì)隨時(shí)間累積,若采用較長積分時(shí)間的觀測數(shù)據(jù)作為一次有效觀測,反而增加位置修正的時(shí)間與難度,航天器的軌道無法實(shí)時(shí)修正。因此,在積分脈沖輪廓折疊時(shí),光子數(shù)據(jù)量的選取需要權(quán)衡觀測積分時(shí)間與脈沖輪廓信噪比,在保證脈沖輪廓信噪比的同時(shí),減少觀測積分時(shí)間。在分析過程中發(fā)現(xiàn),利用20 000個(gè)光子數(shù)據(jù)作為一次有效觀測是合適的。

    表3 分組前后計(jì)時(shí)殘差及不確定度的比較

    圖7 XPNAV-1衛(wèi)星214組觀測數(shù)據(jù)的擬合前、后計(jì)時(shí)殘差

    3 結(jié)語

    本文利用XPNAV-1衛(wèi)星發(fā)布的Crab脈沖星X射線觀測數(shù)據(jù),以及英國Jodrell Bank天文臺(tái)發(fā)布的Crab脈沖星的射電星歷,分析了計(jì)時(shí)結(jié)果。首先得到了35組數(shù)據(jù)對(duì)應(yīng)的TOA和擬合前計(jì)時(shí)殘差。將光子數(shù)過少的第5組數(shù)據(jù)剔除后,得到34組數(shù)據(jù)擬合前計(jì)時(shí)殘差的RMS值為27.117 6ms,擬合后計(jì)時(shí)殘差的RMS值為22.525 6ms。將數(shù)據(jù)量最多的第22、24組數(shù)據(jù)按光子個(gè)數(shù)分組,并對(duì)分組數(shù)據(jù)處理,得到了光子個(gè)數(shù)與積分脈沖輪廓SNR及TOA不確定度的關(guān)系,驗(yàn)證了TOA測量精度隨著光子個(gè)數(shù)的增加不斷提高的結(jié)論。通過分析脈沖輪廓波峰的有無與波形的失真情況,發(fā)現(xiàn)20 000個(gè)光子可以形成相對(duì)穩(wěn)定的積分脈沖輪廓。將去掉第5組后的34組數(shù)據(jù)按每組20 000個(gè)光子重新分組,得到了214個(gè)TOA,求得擬合前計(jì)時(shí)殘差的RMS值為61.836 7μs,擬合后計(jì)時(shí)殘差的RMS值為61.314 1 μs。與分組之前的結(jié)果相比,雖然擬合后的計(jì)時(shí)殘差變大,但在損失一定脈沖輪廓信噪比的情況下獲得了更多的TOA數(shù)據(jù)。若將此結(jié)果應(yīng)用于脈沖星導(dǎo)航,則定位精度會(huì)提升約10%。

    在應(yīng)用X射線脈沖星的深空探測飛行器導(dǎo)航中,采用較長積分時(shí)間,可以提高TOA測量精度。但適當(dāng)減少TOA觀測積分時(shí)間,可以增加TOA觀測數(shù)量,以便實(shí)時(shí)性地修正飛行器軌道誤差。合理斟酌TOA的測量精度與觀測數(shù)量是非常必要的。本文脈沖到達(dá)時(shí)間是通過將脈沖輪廓折疊后形成的積分脈沖輪廓與標(biāo)準(zhǔn)脈沖輪廓進(jìn)行互相關(guān)得到的,在脈沖輪廓折疊過程中會(huì)損失單個(gè)光子相位信息。而利用最大似然估計(jì)法可以在航天器處直接獲得TOA且無需進(jìn)行脈沖輪廓折疊[29-30]。后續(xù)的工作中,可以考慮用這種方法來確定TOA,提高導(dǎo)航的實(shí)時(shí)性。

    [1] 楊廷高. X射線脈沖星脈沖到達(dá)航天器時(shí)間測量[J]. 空間科學(xué)學(xué)報(bào), 2008, 4(4): 60-64.

    [2] 鄭偉, 王奕迪, 湯國建, 等. X射線脈沖星導(dǎo)航理論與應(yīng)用[M]. 北京: 科學(xué)出版社, 2015.

    [3] SHEMAR S, FRASER G, HEIL L, et al. Towards practical autonomous deep-space navigation using X-Ray pulsar timing[J]. Experimental Astronomy, 2016, 42(2): 101-138.

    [4] RAD A, AZARI L. Determining attitude and position in deep space missions using X ray pulsars[J]. International Journal of Astronomy and Astrophysics, 2014, 4(4): 628-648.

    [5] LORIMER D R, KRAMER M. Handbook of Pulsar Astronomy[M]. Cambridge: Cambridge University Press, 2005.

    [6] HOBBS G, COLES W, MANCHESTER R N, et al. Development of a pulsar-based time-scale[J]. Monthly Notices of the Royal Astronomical Society, 2012(427): 2780-2787.

    [7] PETIT G, TAVELLA P. Pulsars and time scales[J]. Astronomy and Astrophysics, 1996, 308(1): 290-298.

    [8] 童明雷, 楊廷高, 趙成仕, 等. 脈沖星計(jì)時(shí)模型參數(shù)的測量精度分析與估計(jì)[J]. 中國科學(xué)(物理學(xué) 力學(xué) 天文學(xué)), 2017, 47(9): 9503.

    [9] HOBBS G, ARCHIBALD A, ARZOUMANIAN Z, et al. The international pulsar timing array project: using pulsars as a gravitational wave detector[J]. Classical and Quantum Gravity, 2010, 27(8): 4013.

    [10] CHAMPION D J, HOBBS G B, MANCHESTER R N, et al. Measuring the mass of solar system planets using pulsar timing[J]. The Astrophysical Journal Letters, 2010, 720(2): 201-205.

    [11] SHEIKH S I. The use of variable celestial X-ray sources for spacecraft navigation[D]. Maryland: University of Maryland, 2005.

    [12] BECKER W, BERNHARDT M G, JESSNER A. Autonomous spacecraft navigation with pulsars[J]. Acta Futura, 2013(7): 11-28.

    [13] TAYLOR J H. Millisecond pulsars: nature's most stable clocks[J]. Proceedings of the IEEE, 2002, 79(7): 1054-1062.

    [14] GALLOWAY D, BASSA C, WANG Z, et al. Accreting neutron star spins and the equation of state[J]. AIP ConferenceProceedings, 2008(983): 510-518.

    [15] GUO Qian. LM-11 successfully Launched XPNAV-1[J]. Aerospace China, 2016, 17(4): 61-61.

    [16] 北斗網(wǎng). 脈沖星試驗(yàn)01星在軌試驗(yàn)數(shù)據(jù)[EB/OL]. (2017-05-09) [2020-06-01]. http://www.beidou.gov.cn/yw/xwzx/201710/t20171010_824.html.

    [17] EDWARDS R T, HOBBS G B, MANCHESTER R N. TEMPO2, a new pulsar timing package-II. The timing model and precision estimates[J]. Monthly Notices of the Royal Astronomical Society, 2006, 372(4): 1549-1574.

    [18] SHAPIRO I I. Fourth test of general relativity[J]. Physical Review Letters, 1964, 13(26): 789-791.

    [19] PETIT G. From atomic clocks to coordinate times[J]. Proceedings of the International Astronomical Union, 2007, 2(14):478-478.

    [20] International Astronomical Union. IAU 2006 Resolution B3[EB/OL]. (2020-03-12) [2020-06-01]. http://www.iau.org/static/resolutions/IAU2006 Resol3.pdf.

    [21] International Astronomical Union. IAU 2000 Resolution B1[EB/OL]. (2020-03-12) [2020-06-01]. https://www.iau.org/static/resolutions/IAU2000_French.pdf.

    [22] FUKUSHIMA T. Time ephemeris and general relativistic scale factor[J]. Proceedings of the International Astronomical Union, 2009, 5(261): 89-94.

    [23] LYNE A G, PRITCHARD R S, GRAHAM-SMITH F. Jodrell Bank Crab Pulsar Monthly Ephemeris[EB/OL]. (2020-03-21) [2020-06-01]. http://www.jb.man.ac.uk/pulsar/crab.html.

    [24] Jet Propulsion Laboratory. Planetary Ephemeris DE200[EB/OL]. (2020-02-15) [2020-06-01]. https://naif.jpl.nasa.gov/pub/naif/generic_kernels/.

    [25] 楊廷高, 童明雷, 趙成仕, 等. Crab脈沖星X射線計(jì)時(shí)觀測數(shù)據(jù)處理與分析[J]. 天文學(xué)報(bào), 2018, 59(2): 10-16.

    [26] TAYLOR J H. Pulsar timing and relativistic gravity[J]. Philosophical Transactions: Physical Sciences and Engineering, 1992, 341(1660): 117-134.

    [27] 朱鴻旭, 童明雷, 楊廷高, 等. XPNAV-1衛(wèi)星先期發(fā)布數(shù)據(jù)的計(jì)時(shí)分析[J]. 宇航學(xué)報(bào), 2019, 40(12): 1492-1500.

    [28] 帥平, 李明, 陳紹龍, 等. X射線脈沖星導(dǎo)航系統(tǒng)原理與方法[M]. 北京: 中國宇航出版社, 2009: 431-434.

    [29] RAY P S, KERR M, PARENT D, et al. Precise gamma-ray timing and radio observations of 17 fermi gamma-ray pulsars[J]. Astrophysical Journal Supplement, 2011, 194(2): 17.

    [30] LIVINGSTONE M A, SCOTT S M, CAMILO F, et al. X-ray and radio timing of the pulsar in 3C 58[J]. Astrophysical Journal, 2009, 706(2): 1163-1173.

    Analysis of the influences of TOA number on timing precision and navigation

    HAN Meng-na1,2,3, TONG Ming-lei1,2, ZHU Hong-xu4,ZHAO Cheng-shi1,2, ZHU Xing-zhi1,2

    (1. National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China;2. Key Laboratory of Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600, China;3. University of Chinese Academy of Sciences, Beijing 100049, China;4. Shandong Institute of Space Electronic Technology, Yantai 264000, China)

    Using the first public data set of Crab pulsar observed by XPNAV-1 satellite, we obtained the resulting pulse times of arrival (TOA) and the uncertainty using the Taylor frequency domain correlation method. After removing the fifth group of data for its too few photons, we obtained the root-mean-square (RMS) of the pre-fit and post-fit timing residuals of the 34 groups left to be about 27.1 μs and 22.5 μs, respectively. The data in the 24th group, which is with the most photons, were divided according to various photon quantities. The relations between the signal to noise ratio (SNR) of the pulse profiles, and the measurement uncertainties of TOA and the timing residuals, and the number of photons were obtained. On the real-time positioning requirements from spacecrafts, the effect of increasing the number of TOAs while maintaining a given SNR on the positioning precision is analyzed. By regrouping the whole data set by 20 000 photons, 214 TOAs were obtained, and the RMS of the pre-fit and post-fit residuals are about 61.8 μs and 61.3 μs. It not only largely enhances the real-time of the navigation, but also improves the positioning precision by about 10 percent.

    X-ray pulsar; pulsar timing; number of photons; pulse profile

    10.13875/j.issn.1674-0637.2021-01-0066-11

    韓孟納, 童明雷, 朱鴻旭, 等. X射線脈沖星TOA數(shù)量對(duì)計(jì)時(shí)精度和導(dǎo)航的影響分析[J]. 時(shí)間頻率學(xué)報(bào), 2021, 44(1): 66-76.

    2020-06-12;

    2020-07-20

    國家自然科學(xué)基金資助項(xiàng)目(U1831130;U1531112;11873050);中國科學(xué)院青年創(chuàng)新促進(jìn)會(huì)資助項(xiàng)目(2017450)

    猜你喜歡
    脈沖星計(jì)時(shí)航天器
    “中國天眼”已發(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
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    國際太空(2022年7期)2022-08-16 09:52:50
    暢游計(jì)時(shí)天地
    車迷(2022年1期)2022-03-29 00:50:24
    腕表計(jì)時(shí)2.0
    中國化妝品(2020年9期)2020-10-09 08:56:56
    12時(shí)計(jì)時(shí)法與24時(shí)計(jì)時(shí)法的互化
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    國際太空(2019年9期)2019-10-23 01:55:34
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    國際太空(2018年12期)2019-01-28 12:53:20
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    國際太空(2018年9期)2018-10-18 08:51:32
    24時(shí)計(jì)時(shí)法
    av视频免费观看在线观看| www.精华液| 王馨瑶露胸无遮挡在线观看| 女人高潮潮喷娇喘18禁视频| 黄色一级大片看看| 中文字幕人妻熟女乱码| 久久青草综合色| 秋霞在线观看毛片| 午夜福利,免费看| 精品福利永久在线观看| 美女扒开内裤让男人捅视频| 久久精品亚洲熟妇少妇任你| 无遮挡黄片免费观看| 国产精品一区二区在线观看99| 欧美精品一区二区免费开放| 午夜免费鲁丝| 激情视频va一区二区三区| 考比视频在线观看| av在线观看视频网站免费| 成人三级做爰电影| 777久久人妻少妇嫩草av网站| 国产日韩欧美在线精品| 中文字幕人妻丝袜制服| 啦啦啦在线观看免费高清www| 一二三四中文在线观看免费高清| 久久精品久久久久久久性| 国产又色又爽无遮挡免| 91aial.com中文字幕在线观看| 国产高清国产精品国产三级| 国产精品一二三区在线看| 亚洲视频免费观看视频| 亚洲美女黄色视频免费看| 国产成人啪精品午夜网站| 精品久久久精品久久久| 久久精品亚洲熟妇少妇任你| 激情视频va一区二区三区| 女的被弄到高潮叫床怎么办| videos熟女内射| 日本欧美视频一区| 宅男免费午夜| 国产精品.久久久| 欧美日韩一区二区视频在线观看视频在线| 亚洲av国产av综合av卡| 亚洲欧洲国产日韩| 午夜久久久在线观看| 亚洲国产精品一区三区| 亚洲国产欧美网| 亚洲一区中文字幕在线| 国产精品久久久久久精品古装| 女人高潮潮喷娇喘18禁视频| 人妻 亚洲 视频| 大香蕉久久网| 精品免费久久久久久久清纯 | 国产精品国产av在线观看| 欧美黄色片欧美黄色片| 亚洲第一区二区三区不卡| 国产精品成人在线| 桃花免费在线播放| 国产片特级美女逼逼视频| 欧美97在线视频| 制服丝袜香蕉在线| av福利片在线| 搡老乐熟女国产| 国产视频首页在线观看| 蜜桃在线观看..| 精品国产露脸久久av麻豆| 欧美激情高清一区二区三区 | 久久精品国产综合久久久| 少妇人妻精品综合一区二区| 中文字幕亚洲精品专区| 国产日韩欧美在线精品| 久久女婷五月综合色啪小说| 久久人妻熟女aⅴ| 十八禁高潮呻吟视频| 99re6热这里在线精品视频| 最近最新中文字幕大全免费视频 | 一级毛片电影观看| 久久精品亚洲熟妇少妇任你| 欧美日韩亚洲高清精品| 免费看不卡的av| 18禁国产床啪视频网站| 日韩视频在线欧美| 日韩熟女老妇一区二区性免费视频| a级毛片在线看网站| 亚洲av成人精品一二三区| 制服人妻中文乱码| 国产日韩一区二区三区精品不卡| 亚洲国产最新在线播放| 久久久精品区二区三区| av线在线观看网站| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 国产av码专区亚洲av| 午夜日本视频在线| 男女边吃奶边做爰视频| 日韩免费高清中文字幕av| 黄色 视频免费看| 午夜福利乱码中文字幕| 久久国产精品大桥未久av| 51午夜福利影视在线观看| 亚洲欧美成人综合另类久久久| 操出白浆在线播放| 久久久久国产一级毛片高清牌| 精品一区二区三区av网在线观看 | 一边摸一边抽搐一进一出视频| 飞空精品影院首页| 女人高潮潮喷娇喘18禁视频| 十八禁网站网址无遮挡| 大陆偷拍与自拍| 久久久久精品久久久久真实原创| 热re99久久国产66热| 成年av动漫网址| 亚洲精品美女久久久久99蜜臀 | 一区二区av电影网| 天堂8中文在线网| 在现免费观看毛片| 男人舔女人的私密视频| 亚洲少妇的诱惑av| 国产精品香港三级国产av潘金莲 | 亚洲成av片中文字幕在线观看| 久久久久精品人妻al黑| www.av在线官网国产| 亚洲国产欧美一区二区综合| 我的亚洲天堂| 嫩草影院入口| 悠悠久久av| 国产精品麻豆人妻色哟哟久久| 尾随美女入室| 男女免费视频国产| 国产精品久久久久久精品古装| 国产在线视频一区二区| 19禁男女啪啪无遮挡网站| 中文欧美无线码| 可以免费在线观看a视频的电影网站 | 一区二区av电影网| √禁漫天堂资源中文www| 久久久久精品国产欧美久久久 | 黄频高清免费视频| 99热国产这里只有精品6| 国产一区二区在线观看av| 国产精品一区二区在线不卡| av在线播放精品| 亚洲成人一二三区av| 亚洲视频免费观看视频| av有码第一页| 免费高清在线观看视频在线观看| 制服人妻中文乱码| 老司机亚洲免费影院| 久久精品国产a三级三级三级| 人妻一区二区av| 男女边摸边吃奶| 久久这里只有精品19| 久久青草综合色| 9热在线视频观看99| 亚洲国产成人一精品久久久| 欧美97在线视频| 久久国产精品大桥未久av| 黄网站色视频无遮挡免费观看| 亚洲国产欧美日韩在线播放| 一级a爱视频在线免费观看| 国产日韩欧美在线精品| 一级片免费观看大全| 国产成人免费观看mmmm| netflix在线观看网站| 中文字幕制服av| 欧美变态另类bdsm刘玥| 精品国产露脸久久av麻豆| 91aial.com中文字幕在线观看| 欧美日韩亚洲国产一区二区在线观看 | 男女床上黄色一级片免费看| 熟女av电影| 1024视频免费在线观看| 中文字幕最新亚洲高清| 9191精品国产免费久久| 国产精品久久久久久精品古装| 国产毛片在线视频| 国产黄色免费在线视频| 搡老乐熟女国产| 久久人人爽人人片av| avwww免费| 美女大奶头黄色视频| 日韩不卡一区二区三区视频在线| 9色porny在线观看| 亚洲精品国产区一区二| 99re6热这里在线精品视频| 久久久亚洲精品成人影院| 搡老乐熟女国产| 曰老女人黄片| 在线观看一区二区三区激情| 亚洲欧美精品自产自拍| 日本猛色少妇xxxxx猛交久久| 黄色视频在线播放观看不卡| tube8黄色片| 中文字幕人妻丝袜制服| 一本色道久久久久久精品综合| 老司机影院成人| 无遮挡黄片免费观看| 欧美激情极品国产一区二区三区| 18在线观看网站| www.精华液| 国产精品嫩草影院av在线观看| 久久久久久人妻| 久久久久久久久久久免费av| 欧美精品一区二区免费开放| 久久99热这里只频精品6学生| 亚洲精品第二区| 久久婷婷青草| 人成视频在线观看免费观看| 一区在线观看完整版| 欧美国产精品一级二级三级| 精品卡一卡二卡四卡免费| 黄网站色视频无遮挡免费观看| avwww免费| 99精国产麻豆久久婷婷| 国产极品粉嫩免费观看在线| 啦啦啦在线免费观看视频4| 999久久久国产精品视频| 人妻一区二区av| 亚洲美女黄色视频免费看| 午夜免费鲁丝| 婷婷色av中文字幕| 精品酒店卫生间| 99热网站在线观看| 午夜日韩欧美国产| 久久久国产欧美日韩av| 久久女婷五月综合色啪小说| av片东京热男人的天堂| 高清欧美精品videossex| 青春草视频在线免费观看| 日韩一本色道免费dvd| 免费av中文字幕在线| 一区二区三区激情视频| 日本黄色日本黄色录像| 七月丁香在线播放| 精品卡一卡二卡四卡免费| av线在线观看网站| 最近中文字幕2019免费版| 七月丁香在线播放| 高清欧美精品videossex| 亚洲欧美色中文字幕在线| 久久精品国产亚洲av高清一级| av网站在线播放免费| 老司机在亚洲福利影院| 久久av网站| 一级片免费观看大全| 黄色视频在线播放观看不卡| 欧美97在线视频| 男女之事视频高清在线观看 | 欧美日韩一区二区视频在线观看视频在线| 极品人妻少妇av视频| 亚洲精品av麻豆狂野| 丝瓜视频免费看黄片| 大陆偷拍与自拍| 色吧在线观看| 国产精品熟女久久久久浪| 国产xxxxx性猛交| 亚洲欧美一区二区三区黑人| www.精华液| 在线观看www视频免费| 久久av网站| av福利片在线| 王馨瑶露胸无遮挡在线观看| 亚洲精品一区蜜桃| 成年美女黄网站色视频大全免费| 日本猛色少妇xxxxx猛交久久| 国产 精品1| 黄片小视频在线播放| 91aial.com中文字幕在线观看| 一级片'在线观看视频| 久久久久精品性色| 成年人免费黄色播放视频| 日本av免费视频播放| 亚洲欧美清纯卡通| 国产一卡二卡三卡精品 | 男男h啪啪无遮挡| 在线免费观看不下载黄p国产| 午夜福利乱码中文字幕| 这个男人来自地球电影免费观看 | 欧美日韩综合久久久久久| 男女高潮啪啪啪动态图| 亚洲美女搞黄在线观看| 老鸭窝网址在线观看| 男女边吃奶边做爰视频| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 国产免费视频播放在线视频| 巨乳人妻的诱惑在线观看| 成人手机av| 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 欧美在线一区亚洲| 久久热在线av| 精品亚洲成a人片在线观看| 国产成人精品无人区| 黄色怎么调成土黄色| 日日啪夜夜爽| 最近最新中文字幕大全免费视频 | 大片免费播放器 马上看| 秋霞伦理黄片| bbb黄色大片| 国产在视频线精品| 日韩 欧美 亚洲 中文字幕| 狠狠婷婷综合久久久久久88av| 国产一区二区 视频在线| 亚洲av在线观看美女高潮| 婷婷色麻豆天堂久久| 校园人妻丝袜中文字幕| 欧美精品高潮呻吟av久久| 欧美精品一区二区免费开放| 亚洲欧美一区二区三区久久| 一区二区日韩欧美中文字幕| 精品一区二区免费观看| 建设人人有责人人尽责人人享有的| 99国产精品免费福利视频| 日日撸夜夜添| 亚洲国产欧美日韩在线播放| 亚洲欧美精品自产自拍| 啦啦啦视频在线资源免费观看| av在线老鸭窝| 久久精品亚洲av国产电影网| 亚洲婷婷狠狠爱综合网| 深夜精品福利| 亚洲国产中文字幕在线视频| 人人妻,人人澡人人爽秒播 | 日本黄色日本黄色录像| av又黄又爽大尺度在线免费看| 丰满饥渴人妻一区二区三| 黄色视频不卡| 2018国产大陆天天弄谢| 超碰97精品在线观看| 国产伦人伦偷精品视频| 最新在线观看一区二区三区 | 日韩欧美精品免费久久| 两个人免费观看高清视频| 别揉我奶头~嗯~啊~动态视频 | 在线观看免费午夜福利视频| 中文字幕精品免费在线观看视频| 两个人免费观看高清视频| 操出白浆在线播放| 国语对白做爰xxxⅹ性视频网站| 考比视频在线观看| 国产97色在线日韩免费| 丝袜脚勾引网站| 亚洲成色77777| 国产亚洲欧美精品永久| 久久人人97超碰香蕉20202| 亚洲国产欧美一区二区综合| 少妇猛男粗大的猛烈进出视频| 日韩电影二区| 日本一区二区免费在线视频| 日韩 亚洲 欧美在线| 香蕉丝袜av| 综合色丁香网| 国产又爽黄色视频| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 成年人免费黄色播放视频| 激情五月婷婷亚洲| 欧美日韩亚洲高清精品| 亚洲一区二区三区欧美精品| av不卡在线播放| 亚洲精品美女久久av网站| 国产成人欧美在线观看 | 国产男人的电影天堂91| 久久精品久久久久久久性| 亚洲欧美色中文字幕在线| videosex国产| 一本—道久久a久久精品蜜桃钙片| 欧美国产精品va在线观看不卡| 美女中出高潮动态图| 蜜桃在线观看..| 可以免费在线观看a视频的电影网站 | 国产免费又黄又爽又色| 69精品国产乱码久久久| 秋霞在线观看毛片| 久久人妻熟女aⅴ| 中文字幕av电影在线播放| 在线观看免费视频网站a站| 少妇人妻 视频| 黄网站色视频无遮挡免费观看| 大话2 男鬼变身卡| 久久国产亚洲av麻豆专区| 尾随美女入室| 久久久久精品久久久久真实原创| 最近的中文字幕免费完整| 丝袜美腿诱惑在线| 最近中文字幕高清免费大全6| 久久精品aⅴ一区二区三区四区| 中文字幕人妻熟女乱码| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 久久久久久久久免费视频了| 欧美黑人精品巨大| 十分钟在线观看高清视频www| 久久久久久免费高清国产稀缺| 中文天堂在线官网| 日韩欧美精品免费久久| 19禁男女啪啪无遮挡网站| 午夜激情久久久久久久| 亚洲国产精品国产精品| 欧美人与性动交α欧美精品济南到| 国产精品一区二区精品视频观看| 成年美女黄网站色视频大全免费| 日日爽夜夜爽网站| 久久久久精品人妻al黑| 欧美另类一区| 天天躁狠狠躁夜夜躁狠狠躁| av免费观看日本| 色婷婷久久久亚洲欧美| 亚洲成人av在线免费| 搡老乐熟女国产| 在线观看人妻少妇| 婷婷色综合大香蕉| 国产精品久久久久成人av| 久久国产精品大桥未久av| 午夜免费鲁丝| 国产精品香港三级国产av潘金莲 | e午夜精品久久久久久久| 一区二区三区激情视频| 国产成人精品无人区| 咕卡用的链子| 黄色一级大片看看| 久久综合国产亚洲精品| 丝瓜视频免费看黄片| 婷婷色麻豆天堂久久| 亚洲国产最新在线播放| 蜜桃在线观看..| 啦啦啦在线观看免费高清www| 99久久人妻综合| 亚洲色图 男人天堂 中文字幕| 日韩 亚洲 欧美在线| 男女之事视频高清在线观看 | 国产激情久久老熟女| videos熟女内射| bbb黄色大片| 亚洲成人一二三区av| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 国产高清不卡午夜福利| 亚洲,一卡二卡三卡| 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区久久| 国产欧美日韩一区二区三区在线| 街头女战士在线观看网站| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人| 黄色视频在线播放观看不卡| 色精品久久人妻99蜜桃| 欧美日韩亚洲高清精品| 亚洲精品乱久久久久久| 热99久久久久精品小说推荐| av在线播放精品| 丝袜美腿诱惑在线| 亚洲五月色婷婷综合| 婷婷色综合大香蕉| 少妇精品久久久久久久| 制服人妻中文乱码| 久久女婷五月综合色啪小说| av女优亚洲男人天堂| 只有这里有精品99| 国产在线一区二区三区精| 一边亲一边摸免费视频| 国产av一区二区精品久久| 两个人看的免费小视频| 亚洲av日韩精品久久久久久密 | 欧美国产精品va在线观看不卡| 亚洲精品自拍成人| 天天躁夜夜躁狠狠躁躁| 日韩中文字幕视频在线看片| 人妻 亚洲 视频| 综合色丁香网| 国产av码专区亚洲av| 亚洲欧美一区二区三区国产| 欧美日韩视频精品一区| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 国产免费现黄频在线看| a级毛片黄视频| 91aial.com中文字幕在线观看| 国产在线一区二区三区精| 2018国产大陆天天弄谢| 久久久久精品人妻al黑| av电影中文网址| 亚洲精品视频女| 18禁观看日本| 在线观看三级黄色| 深夜精品福利| 一本大道久久a久久精品| 婷婷色麻豆天堂久久| 欧美日韩亚洲高清精品| 自拍欧美九色日韩亚洲蝌蚪91| 免费在线观看完整版高清| 大片免费播放器 马上看| 男女免费视频国产| 婷婷色av中文字幕| 中文字幕色久视频| 亚洲成av片中文字幕在线观看| 精品国产国语对白av| 欧美精品一区二区大全| 国产成人欧美| 成人国语在线视频| 赤兔流量卡办理| 亚洲,一卡二卡三卡| 精品国产一区二区三区久久久樱花| 国产精品嫩草影院av在线观看| 大香蕉久久成人网| 精品亚洲成国产av| 久久热在线av| 国产av码专区亚洲av| 一二三四在线观看免费中文在| av免费观看日本| 国产亚洲av片在线观看秒播厂| 国产激情久久老熟女| 亚洲在久久综合| 丰满少妇做爰视频| 亚洲精品久久久久久婷婷小说| 日韩制服丝袜自拍偷拍| 啦啦啦在线观看免费高清www| 国产不卡av网站在线观看| 国产免费视频播放在线视频| 国产成人精品无人区| 一区二区三区激情视频| 国产麻豆69| 亚洲精品久久成人aⅴ小说| 久久热在线av| 国产av精品麻豆| 亚洲欧美精品综合一区二区三区| 91aial.com中文字幕在线观看| 九九爱精品视频在线观看| 精品午夜福利在线看| 99国产精品免费福利视频| 亚洲国产欧美网| 麻豆av在线久日| 黑人巨大精品欧美一区二区蜜桃| 一边亲一边摸免费视频| 大香蕉久久网| 亚洲欧美中文字幕日韩二区| 观看av在线不卡| 亚洲第一av免费看| av国产精品久久久久影院| 成年美女黄网站色视频大全免费| 免费黄色在线免费观看| 日本vs欧美在线观看视频| tube8黄色片| 两个人看的免费小视频| 最近中文字幕2019免费版| 美女国产高潮福利片在线看| 亚洲精品自拍成人| 亚洲av在线观看美女高潮| 午夜久久久在线观看| 久久精品国产亚洲av涩爱| 日本欧美国产在线视频| 国产精品成人在线| 丝袜喷水一区| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 久久久精品区二区三区| 国产男人的电影天堂91| 黄色视频不卡| 国产 精品1| 人人妻人人澡人人看| 国产又色又爽无遮挡免| 亚洲精品国产av成人精品| 91国产中文字幕| 丝袜喷水一区| av国产精品久久久久影院| 亚洲国产欧美一区二区综合| 亚洲av成人不卡在线观看播放网 | 又粗又硬又长又爽又黄的视频| 日韩,欧美,国产一区二区三区| 夜夜骑夜夜射夜夜干| 最新的欧美精品一区二区| 精品亚洲成国产av| 国产精品一区二区在线不卡| 卡戴珊不雅视频在线播放| 国产精品 国内视频| 中文字幕亚洲精品专区| 黑人猛操日本美女一级片| 国产高清不卡午夜福利| 精品少妇一区二区三区视频日本电影 | 国产免费视频播放在线视频| 1024视频免费在线观看| 极品少妇高潮喷水抽搐| 日本av手机在线免费观看| 大话2 男鬼变身卡| 嫩草影视91久久| 午夜久久久在线观看| 国产精品一区二区在线观看99| 一级黄片播放器| 精品午夜福利在线看| videosex国产| 性色av一级| 国产激情久久老熟女| 国产av精品麻豆| 国产精品久久久久久久久免| 一级毛片电影观看| 国产伦人伦偷精品视频| 亚洲专区中文字幕在线 | 欧美最新免费一区二区三区| 亚洲精品在线美女| 日本猛色少妇xxxxx猛交久久| 黑丝袜美女国产一区| 婷婷色av中文字幕| 又粗又硬又长又爽又黄的视频| 久久女婷五月综合色啪小说| 妹子高潮喷水视频| 国产精品免费大片| 丝袜美腿诱惑在线| 中文字幕制服av| 成人18禁高潮啪啪吃奶动态图| 国产成人精品无人区| 高清在线视频一区二区三区| 成人漫画全彩无遮挡|