• <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í)法
    成年女人在线观看亚洲视频 | 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| av网站免费在线观看视频 | 最后的刺客免费高清国语| freevideosex欧美| 欧美 日韩 精品 国产| 国产成人福利小说| 黄色配什么色好看| 欧美丝袜亚洲另类| 日韩制服骚丝袜av| 日日摸夜夜添夜夜添av毛片| 亚洲国产精品成人久久小说| 少妇人妻一区二区三区视频| 日本三级黄在线观看| 久久精品综合一区二区三区| 亚洲电影在线观看av| 久久人人爽人人片av| 97在线视频观看| 精品国产一区二区三区久久久樱花 | 亚洲国产成人一精品久久久| 91久久精品国产一区二区三区| 免费观看av网站的网址| 九九爱精品视频在线观看| 校园人妻丝袜中文字幕| 亚洲av国产av综合av卡| 国产日韩欧美在线精品| 在线观看免费高清a一片| 全区人妻精品视频| 婷婷色av中文字幕| 建设人人有责人人尽责人人享有的 | 国产午夜精品论理片| 99久久精品一区二区三区| 中文在线观看免费www的网站| 久久精品国产亚洲av天美| 亚洲精品乱码久久久v下载方式| 嘟嘟电影网在线观看| 亚洲成人一二三区av| 国产精品熟女久久久久浪| 蜜臀久久99精品久久宅男| 小蜜桃在线观看免费完整版高清| 久久韩国三级中文字幕| 91av网一区二区| 国产欧美日韩精品一区二区| 一区二区三区免费毛片| 亚洲美女视频黄频| 伊人久久精品亚洲午夜| 免费看av在线观看网站| 国产精品日韩av在线免费观看| av网站免费在线观看视频 | 国产人妻一区二区三区在| 国产精品av视频在线免费观看| 精品少妇黑人巨大在线播放| 免费看日本二区| 噜噜噜噜噜久久久久久91| 久久久久久九九精品二区国产| 精品人妻一区二区三区麻豆| 久久久久九九精品影院| 九九久久精品国产亚洲av麻豆| 看十八女毛片水多多多| 欧美日韩亚洲高清精品| 国产高清国产精品国产三级 | 少妇人妻精品综合一区二区| 欧美精品国产亚洲| 国产熟女欧美一区二区| 极品少妇高潮喷水抽搐| 99久国产av精品国产电影| 日韩制服骚丝袜av| 美女脱内裤让男人舔精品视频| 亚洲成人一二三区av| 色播亚洲综合网| 久久精品国产自在天天线| 国产成人一区二区在线| 在现免费观看毛片| 男女边吃奶边做爰视频| 日本黄大片高清| 免费黄频网站在线观看国产| 国产高清有码在线观看视频| 最近最新中文字幕免费大全7| 熟女电影av网| 老女人水多毛片| 男插女下体视频免费在线播放| 亚洲熟妇中文字幕五十中出| 91午夜精品亚洲一区二区三区| 中文字幕久久专区| 精品人妻熟女av久视频| 久久久精品欧美日韩精品| 国产一区二区三区av在线| av在线亚洲专区| 亚洲精品中文字幕在线视频 | 成年av动漫网址| 国产在视频线精品| 亚洲国产精品专区欧美| 人体艺术视频欧美日本| 又黄又爽又刺激的免费视频.| 麻豆成人午夜福利视频| 国产亚洲午夜精品一区二区久久 | 校园人妻丝袜中文字幕| 亚洲伊人久久精品综合| 精品不卡国产一区二区三区| 亚洲丝袜综合中文字幕| 国产成人aa在线观看| 欧美3d第一页| 高清av免费在线| 少妇丰满av| 亚洲成色77777| 秋霞伦理黄片| 建设人人有责人人尽责人人享有的 | 一级av片app| 水蜜桃什么品种好| 久久99热这里只频精品6学生| 日本午夜av视频| 国内精品一区二区在线观看| 高清在线视频一区二区三区| 欧美xxxx黑人xx丫x性爽| 精品午夜福利在线看| 黄片wwwwww| 国产三级在线视频| 国产永久视频网站| 色网站视频免费| 人人妻人人澡人人爽人人夜夜 | 欧美最新免费一区二区三区| 天天一区二区日本电影三级| 人人妻人人看人人澡| 国产视频内射| 69人妻影院| 亚洲精品乱久久久久久| 日韩人妻高清精品专区| 亚洲精华国产精华液的使用体验| 欧美日韩亚洲高清精品| 一区二区三区乱码不卡18| 亚洲最大成人av| 国产三级在线视频| 国产高清不卡午夜福利| 精品久久久久久久人妻蜜臀av| 91精品国产九色| 天天躁日日操中文字幕| 国产成年人精品一区二区| 汤姆久久久久久久影院中文字幕 | 性插视频无遮挡在线免费观看| 777米奇影视久久| 少妇人妻一区二区三区视频| 精品国产露脸久久av麻豆 | 国产高清有码在线观看视频| 免费电影在线观看免费观看| 日韩av在线免费看完整版不卡| 日韩成人伦理影院| 国产久久久一区二区三区| 99久国产av精品国产电影| or卡值多少钱| 亚洲精品色激情综合| 91午夜精品亚洲一区二区三区| 成人亚洲精品一区在线观看 | 老司机影院毛片| 久久久精品免费免费高清| 日韩欧美国产在线观看| 美女cb高潮喷水在线观看| 美女内射精品一级片tv| 国精品久久久久久国模美| av又黄又爽大尺度在线免费看| 亚洲自拍偷在线| 国产精品1区2区在线观看.| 久久这里有精品视频免费| 亚洲欧美中文字幕日韩二区| 观看美女的网站| 欧美精品国产亚洲| 2021天堂中文幕一二区在线观| 亚洲国产高清在线一区二区三| 亚洲av成人精品一区久久| 特大巨黑吊av在线直播| 午夜爱爱视频在线播放| 日本av手机在线免费观看| 久久久久久久久久黄片| 国产午夜精品久久久久久一区二区三区| 免费电影在线观看免费观看| 亚洲av国产av综合av卡| 日本熟妇午夜| 免费看不卡的av| 内地一区二区视频在线| 91久久精品国产一区二区成人| 99热全是精品| 国产一区二区亚洲精品在线观看| 亚洲av.av天堂| 亚洲激情五月婷婷啪啪| 人妻夜夜爽99麻豆av| 国产探花极品一区二区| 又黄又爽又刺激的免费视频.| 久久久久久久久久久丰满| 国产精品美女特级片免费视频播放器| 一区二区三区乱码不卡18| 18禁在线无遮挡免费观看视频| 欧美性感艳星| 日韩强制内射视频| 99热网站在线观看| 欧美日韩精品成人综合77777| 国产午夜精品久久久久久一区二区三区| 国产v大片淫在线免费观看| 国产人妻一区二区三区在| 青春草国产在线视频| 大又大粗又爽又黄少妇毛片口| 亚洲第一区二区三区不卡| 国产乱人视频| 久久国产乱子免费精品| 亚洲av免费在线观看| 美女xxoo啪啪120秒动态图| 国精品久久久久久国模美| 色播亚洲综合网| 午夜福利视频1000在线观看| 久久久久精品性色| 国产亚洲一区二区精品| 丝瓜视频免费看黄片| 成人午夜精彩视频在线观看| 99re6热这里在线精品视频| 免费av不卡在线播放| 一级毛片黄色毛片免费观看视频| 午夜福利在线在线| 日日撸夜夜添| 我要看日韩黄色一级片| 最后的刺客免费高清国语| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久久免| 亚洲av在线观看美女高潮| 啦啦啦韩国在线观看视频| 久久久久久久久久成人| 欧美三级亚洲精品| 久久99蜜桃精品久久| 高清日韩中文字幕在线| 久久精品综合一区二区三区| 日日干狠狠操夜夜爽| 国产淫片久久久久久久久| 99热这里只有是精品在线观看| 久久精品国产亚洲av天美| 午夜免费观看性视频| 欧美区成人在线视频| 亚洲乱码一区二区免费版| 菩萨蛮人人尽说江南好唐韦庄| 日本猛色少妇xxxxx猛交久久| 搡老乐熟女国产| av卡一久久| 水蜜桃什么品种好| 亚洲国产精品国产精品| 国产伦在线观看视频一区| 三级国产精品欧美在线观看| 精华霜和精华液先用哪个| 五月玫瑰六月丁香| 一二三四中文在线观看免费高清| 精品一区二区三区视频在线| 国产单亲对白刺激| 高清在线视频一区二区三区| 丝袜喷水一区| 国产伦理片在线播放av一区| 亚洲成色77777| 91午夜精品亚洲一区二区三区| 国产精品久久久久久精品电影| 国产精品精品国产色婷婷| videossex国产| 麻豆精品久久久久久蜜桃| 精品久久久久久久久久久久久| 久久久久免费精品人妻一区二区| 日本-黄色视频高清免费观看| 日本与韩国留学比较| 国产伦精品一区二区三区视频9| 久久久久久国产a免费观看| 欧美日韩国产mv在线观看视频 | 内射极品少妇av片p| 国产免费福利视频在线观看| 日本wwww免费看| 搡老妇女老女人老熟妇| 午夜爱爱视频在线播放| 直男gayav资源| 黄色一级大片看看| 亚洲精品自拍成人| 日韩精品有码人妻一区| 中国国产av一级| 一区二区三区免费毛片| 亚洲国产精品成人久久小说| 又爽又黄a免费视频| 特级一级黄色大片| 男女边摸边吃奶| 亚洲国产欧美在线一区| 亚洲欧美日韩卡通动漫| av女优亚洲男人天堂| 国产乱人偷精品视频| 久久鲁丝午夜福利片| 国产永久视频网站| 我的女老师完整版在线观看| 日本av手机在线免费观看| 青春草国产在线视频| 丝袜喷水一区| 九九久久精品国产亚洲av麻豆| a级毛片免费高清观看在线播放| 亚洲四区av| 亚洲精品乱码久久久v下载方式| 伊人久久精品亚洲午夜| 国产黄频视频在线观看| 亚洲欧洲日产国产| 欧美成人一区二区免费高清观看| 中文字幕免费在线视频6| 日韩电影二区| 在线a可以看的网站| 婷婷色综合大香蕉| 九色成人免费人妻av| 国产在视频线精品| 亚洲不卡免费看| 黄色一级大片看看| 亚洲国产精品sss在线观看| 一级a做视频免费观看| 国产成人aa在线观看| 成人一区二区视频在线观看| 亚洲av二区三区四区| 国产精品国产三级专区第一集| 亚洲熟妇中文字幕五十中出| 日本欧美国产在线视频| 一级av片app| 少妇熟女aⅴ在线视频| 国产 一区精品| 亚洲最大成人av| 久99久视频精品免费| 亚洲不卡免费看| 内射极品少妇av片p| 神马国产精品三级电影在线观看| 一级片'在线观看视频| 男插女下体视频免费在线播放| 亚洲不卡免费看| 国产探花极品一区二区| 国产av在哪里看| 听说在线观看完整版免费高清| 免费黄色在线免费观看| a级一级毛片免费在线观看| 亚洲欧美一区二区三区黑人 | 午夜精品国产一区二区电影 | 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 日韩av不卡免费在线播放| 欧美另类一区| 亚洲欧美中文字幕日韩二区| av在线老鸭窝| 免费观看av网站的网址| 亚洲精品色激情综合| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 国产午夜福利久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久av不卡| 尤物成人国产欧美一区二区三区| 亚洲成人av在线免费| 国产伦精品一区二区三区视频9| 国产高清有码在线观看视频| 成人二区视频| 久久久久久久久久成人| 精品国产三级普通话版| 国产亚洲精品久久久com| 国产伦理片在线播放av一区| 婷婷色麻豆天堂久久| 亚洲伊人久久精品综合| 亚洲av电影在线观看一区二区三区 | 亚洲国产精品sss在线观看| 成人特级av手机在线观看| 国产亚洲精品久久久com| 久久精品国产鲁丝片午夜精品| 高清视频免费观看一区二区 | 五月天丁香电影| 国产亚洲av片在线观看秒播厂 | 国产三级在线视频| 亚洲欧洲日产国产| 美女cb高潮喷水在线观看| 亚洲av电影不卡..在线观看| 高清视频免费观看一区二区 | 欧美成人a在线观看| 晚上一个人看的免费电影| 赤兔流量卡办理| 午夜福利高清视频| 在线免费观看的www视频| 国产精品久久久久久av不卡| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 亚洲精品日韩在线中文字幕| 亚洲av成人av| 午夜免费激情av| 国产91av在线免费观看| 亚洲欧美日韩东京热| 亚洲精品,欧美精品| 国产精品爽爽va在线观看网站| 熟女人妻精品中文字幕| 国产淫语在线视频| 国产成人精品久久久久久| 国产精品久久久久久精品电影| 国产免费视频播放在线视频 | ponron亚洲| 日韩av免费高清视频| 国内精品宾馆在线| 精品久久久久久久久久久久久| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频 | 观看免费一级毛片| 国产精品一区二区性色av| 偷拍熟女少妇极品色| 国产日韩欧美在线精品| 男女啪啪激烈高潮av片| 一个人观看的视频www高清免费观看| 国产 一区 欧美 日韩| 日韩av在线免费看完整版不卡| 精品久久久久久久久av| 成人鲁丝片一二三区免费| 国产欧美日韩精品一区二区| 伊人久久精品亚洲午夜| 亚洲av中文av极速乱| 一区二区三区免费毛片| 男女那种视频在线观看| 最近中文字幕2019免费版| 床上黄色一级片| 男人和女人高潮做爰伦理| 欧美潮喷喷水| av播播在线观看一区| 免费大片18禁| 日日啪夜夜爽| 在线a可以看的网站| 亚洲精品色激情综合| 简卡轻食公司| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 99九九线精品视频在线观看视频| 久热久热在线精品观看| 色综合色国产| 五月伊人婷婷丁香| 亚洲综合精品二区| 真实男女啪啪啪动态图| 在线天堂最新版资源| 97在线视频观看| 亚洲人与动物交配视频| www.av在线官网国产| 丰满人妻一区二区三区视频av| 国产人妻一区二区三区在| 国产午夜精品一二区理论片| 久久热精品热| 91在线精品国自产拍蜜月| 久久草成人影院| 国产精品av视频在线免费观看| 午夜福利在线观看吧| 日本三级黄在线观看| 三级国产精品欧美在线观看| 亚州av有码| 日韩亚洲欧美综合| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 91精品伊人久久大香线蕉| 日本一本二区三区精品| 白带黄色成豆腐渣| 天堂√8在线中文| 精品欧美国产一区二区三| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 国产精品女同一区二区软件| 如何舔出高潮| 高清午夜精品一区二区三区| 精品午夜福利在线看| 国产单亲对白刺激| 亚洲欧美日韩东京热| av又黄又爽大尺度在线免费看| 成人av在线播放网站| 国产成人freesex在线| 久久精品国产自在天天线| 美女内射精品一级片tv| 国产色爽女视频免费观看| 国产精品国产三级专区第一集| 欧美一区二区亚洲| 丝瓜视频免费看黄片| 国产永久视频网站| 久久久久久久久久黄片| a级毛色黄片| 午夜免费观看性视频| 国产不卡一卡二| 国内精品宾馆在线| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 国模一区二区三区四区视频| 青春草视频在线免费观看| 日韩av不卡免费在线播放| 91在线精品国自产拍蜜月| 啦啦啦啦在线视频资源| 九草在线视频观看| 黑人高潮一二区| 亚洲精品第二区| 国产一区二区三区综合在线观看 | 午夜福利高清视频| 一级a做视频免费观看| 中文字幕亚洲精品专区| 天堂中文最新版在线下载 | 免费观看在线日韩| av免费在线看不卡| 亚洲成人久久爱视频| 国产毛片a区久久久久| 国产国拍精品亚洲av在线观看| 99热这里只有是精品在线观看| 国产有黄有色有爽视频| 高清av免费在线| 久久久久久久国产电影| 久久久久久九九精品二区国产| 久久综合国产亚洲精品| 又黄又爽又刺激的免费视频.| 国产精品国产三级专区第一集| 成人无遮挡网站| av专区在线播放| 国内揄拍国产精品人妻在线| 中文乱码字字幕精品一区二区三区 | av福利片在线观看| 亚洲欧洲日产国产| 99热这里只有精品一区| av.在线天堂| 69av精品久久久久久| 国产色爽女视频免费观看| 日本与韩国留学比较| 久久久色成人| 搞女人的毛片| 久久久久久国产a免费观看| 久热久热在线精品观看| 男女边摸边吃奶| 卡戴珊不雅视频在线播放| 欧美潮喷喷水| 菩萨蛮人人尽说江南好唐韦庄| 欧美xxxx性猛交bbbb| 九九在线视频观看精品| 国产高清国产精品国产三级 | 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 成人综合一区亚洲| 久久精品夜夜夜夜夜久久蜜豆| 三级经典国产精品| 乱人视频在线观看| 国产在线一区二区三区精| 午夜福利在线在线| 久久久久网色| h日本视频在线播放| 欧美性猛交╳xxx乱大交人| 欧美一级a爱片免费观看看| 最新中文字幕久久久久| 国产人妻一区二区三区在| 国产熟女欧美一区二区| 中国国产av一级| 真实男女啪啪啪动态图| 黄片wwwwww| 亚洲成人久久爱视频| 大话2 男鬼变身卡| 一二三四中文在线观看免费高清| 欧美激情国产日韩精品一区| 国产精品人妻久久久影院| 尤物成人国产欧美一区二区三区| 神马国产精品三级电影在线观看| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 国产成人福利小说| 丝袜喷水一区| 91av网一区二区| 女人被狂操c到高潮| 亚洲电影在线观看av| 亚洲国产最新在线播放| 伦精品一区二区三区| 久久久精品欧美日韩精品| 老女人水多毛片| 国产 亚洲一区二区三区 | 国产精品嫩草影院av在线观看| 网址你懂的国产日韩在线| 欧美性猛交╳xxx乱大交人| 国产色婷婷99| 天堂√8在线中文| 狠狠精品人妻久久久久久综合| 蜜臀久久99精品久久宅男| 大香蕉97超碰在线| 欧美xxⅹ黑人| 嫩草影院入口| 免费大片黄手机在线观看| 成人特级av手机在线观看| 午夜免费观看性视频| 我的女老师完整版在线观看| 观看免费一级毛片| 97人妻精品一区二区三区麻豆| 韩国高清视频一区二区三区| 午夜福利网站1000一区二区三区| 国产精品久久视频播放| 国产毛片a区久久久久| 免费看美女性在线毛片视频| 免费看不卡的av| 天堂√8在线中文| 可以在线观看毛片的网站| 最近最新中文字幕免费大全7| 国产亚洲av片在线观看秒播厂 | 国产精品精品国产色婷婷| 日韩视频在线欧美| 亚洲一级一片aⅴ在线观看| 狠狠精品人妻久久久久久综合| 久久99蜜桃精品久久| 午夜激情福利司机影院| 1000部很黄的大片| 精品欧美国产一区二区三| 国产老妇女一区| 婷婷色综合www| 国产精品99久久久久久久久| 2018国产大陆天天弄谢| 午夜免费观看性视频| 嫩草影院精品99| av在线亚洲专区| 18禁在线播放成人免费| 精品少妇黑人巨大在线播放| 两个人视频免费观看高清| 国产亚洲精品av在线| 美女xxoo啪啪120秒动态图| 国产在线男女| 国产亚洲午夜精品一区二区久久 | 国产精品久久久久久久久免| 91aial.com中文字幕在线观看| 99热6这里只有精品| 2022亚洲国产成人精品| 成人av在线播放网站| 永久免费av网站大全| 美女cb高潮喷水在线观看|