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

    基于雙譜濾波的綜合脈沖星時(shí)算法研究?

    2021-03-29 12:33:00周慶勇魏子卿劉思偉孫鵬飛奮姜劉曉剛
    天文學(xué)報(bào) 2021年2期
    關(guān)鍵詞:雙譜脈沖星穩(wěn)定度

    周慶勇 魏子卿 張 華 劉思偉 孫鵬飛 張 奮姜 坤 劉曉剛 李 奎

    (1 地理信息工程國家重點(diǎn)實(shí)驗(yàn)室西安710054)(2 西安測(cè)繪研究所西安710054)(3 信息工程大學(xué)地理空間信息學(xué)院鄭州450051)(4 西安電子科技大學(xué)空間科學(xué)與技術(shù)學(xué)院西安710126)(5 中國科學(xué)院國家授時(shí)中心西安710600)(6 北京通信與跟蹤技術(shù)研究所北京100090)

    1 引言

    脈沖星是一類高速自轉(zhuǎn)的中子星, 具有穩(wěn)定度極高的自轉(zhuǎn)頻率, 部分毫秒脈沖星自轉(zhuǎn)周期長期穩(wěn)定性與原子鐘相當(dāng)甚至更優(yōu)[1], 這表明脈沖星能夠提供一種基于遙遠(yuǎn)自然天體并持續(xù)數(shù)百萬乃至數(shù)十億年的獨(dú)立時(shí)間頻率, 并且長期穩(wěn)定性優(yōu)于原子鐘. 基于脈沖星高穩(wěn)定度自轉(zhuǎn)頻率可構(gòu)建一種新的時(shí)間尺度, 稱為脈沖星時(shí). 脈沖星時(shí)具有以下優(yōu)勢(shì): (1)能夠?qū)Φ厍驎r(shí)(Terrestrial time, TT)提供獨(dú)立外部檢核, 檢測(cè)現(xiàn)有原子時(shí)的長期穩(wěn)定性; (2)與原子時(shí)的工作原理不同, 脈沖星時(shí)是建立在恒星質(zhì)量天體的運(yùn)動(dòng)過程上,不易受干擾; (3)原子鐘工作壽命有限,隨著原子鐘器件的老化,其守時(shí)性能逐漸變差,而脈沖星時(shí)可穩(wěn)定持續(xù)工作數(shù)百萬乃至數(shù)十億年, 且作用范圍廣. 脈沖星時(shí)具有高穩(wěn)定性、全自主性和全宇宙性的特點(diǎn), 在時(shí)空基準(zhǔn)建設(shè)方面有良好的應(yīng)用前景. 利用脈沖星計(jì)時(shí)數(shù)據(jù)可改進(jìn)原子時(shí)的長期穩(wěn)定性, 同時(shí)結(jié)合原子鐘的短期穩(wěn)定性和脈沖星的長期可用性、穩(wěn)定性, 可構(gòu)建一個(gè)新的綜合時(shí)間尺度.

    國外學(xué)者在脈沖星發(fā)現(xiàn)不久就意識(shí)到脈沖星自轉(zhuǎn)頻率在時(shí)間基準(zhǔn)方面的應(yīng)用價(jià)值. 1982年, Backer等人發(fā)現(xiàn)極其穩(wěn)定的毫秒脈沖星[2]. 1984年II’in等人提出了脈沖時(shí)間尺度的概念, 并對(duì)脈沖星時(shí)間建立原理進(jìn)行了初步討論[3]. 1988年, Guinot指出脈沖星計(jì)時(shí)有利于對(duì)地球時(shí)的認(rèn)識(shí)[4]. 1991年, Guniot等人分析發(fā)現(xiàn)毫秒脈沖星可改善原子時(shí)間的長期穩(wěn)定性義,可用于建立新的時(shí)間基準(zhǔn)[5]. 1996年P(guān)etit等人提出了綜合脈沖星時(shí)的概念, 研究了綜合脈沖星時(shí)經(jīng)典加權(quán)方法[6]. 2008年以來, Rodin等人研究了基于維納濾波的綜合脈沖星時(shí)加權(quán)平均算法, 并利用帕克斯脈沖星計(jì)時(shí)陣(Parkes Pulsar Timing Array, PPTA)數(shù)據(jù)進(jìn)行了分析[7–8]. 2012年, Hobbs等人提出將鐘信號(hào)包含在計(jì)時(shí)模型中, 采用全局最小二乘擬合來構(gòu)建綜合脈沖星時(shí), 并利用PPTA的19顆毫秒脈沖星觀測(cè)數(shù)據(jù)構(gòu)建了第1個(gè)具有與國際原子時(shí)精度相當(dāng)?shù)拿}沖時(shí)間尺度TT(PPTA11)[9], 數(shù)字11表示采用截至2011年P(guān)PTA數(shù)據(jù). 2019年, Hobbs等人采用譜估計(jì)和貝葉斯估計(jì)方法分別對(duì)國際脈沖星計(jì)時(shí)陣(International Pulsar Timing Array,IPTA)觀測(cè)數(shù)據(jù)噪聲進(jìn)行建模, 并構(gòu)建了更加穩(wěn)定的脈沖星時(shí)間尺度TT (IPTA16),發(fā)現(xiàn)TT (IPTA16)與TT (BIPM17)具有較好的一致性, 并將TT (IPTA16)的功率譜與TT (BIPM2017)功率譜進(jìn)行比較, 發(fā)現(xiàn)脈沖星時(shí)可檢查原子時(shí)間尺度的不穩(wěn)定性[10].BIPM是國際計(jì)量局(International Bureau of Weights and Measures)的簡(jiǎn)稱, 負(fù)責(zé)國際原子時(shí)(International Atomic Time, TAI)的統(tǒng)一和發(fā)布. 近年來, 我國學(xué)者跟蹤國際上脈沖星時(shí)相關(guān)技術(shù)的發(fā)展動(dòng)態(tài), 在綜合脈沖星時(shí)算法方面取得較好的進(jìn)展[11–14]. 2007年仲崇霞[15]將小波分解方法引入綜合脈沖星時(shí)算法構(gòu)建中, 提高了綜合脈沖星時(shí)的穩(wěn)定性. 2011年, 陳鼎等人介紹了脈沖星時(shí)的應(yīng)用, 并分析了基于PPTA項(xiàng)目4顆毫秒脈沖星的綜合脈沖星時(shí)[16]. 2016年, 尹東山等人使用Vondrak濾波處理了北美納米赫茲天文臺(tái)(North American Nanohertz Observatory for Gravitational Waves, NANOGrav) 9 yr的觀測(cè)數(shù)據(jù), 構(gòu)建的脈沖星時(shí)長期穩(wěn)定性為3.4×10?15[17–18]. 2017年, 張彩紅等人利用PPTA數(shù)據(jù)構(gòu)建綜合脈沖星時(shí), 比較了脈沖星時(shí)與原子時(shí)的長期穩(wěn)定性[19].

    脈沖星時(shí)最新發(fā)展方向之一是改善衛(wèi)星導(dǎo)航系統(tǒng)時(shí)間基準(zhǔn)的穩(wěn)定性. 2018年12月,歐空局宣布運(yùn)行了一個(gè)基于脈沖星的時(shí)鐘項(xiàng)目(Pulsar Chronos, PulChron), 其目的是通過觀測(cè)毫秒脈沖星監(jiān)測(cè)和改善伽利略衛(wèi)星導(dǎo)航系統(tǒng)時(shí)間的穩(wěn)定性[20]. 報(bào)道稱PulChron已經(jīng)運(yùn)行了較長時(shí)間, 且其結(jié)果令人鼓舞, 但未見更多詳細(xì)文獻(xiàn), 也可預(yù)見將來脈沖星時(shí)在我國北斗導(dǎo)航系統(tǒng)中的應(yīng)用.

    單顆毫秒脈沖星的自轉(zhuǎn)頻率可定義一個(gè)脈沖星時(shí)(Pulsar Time, PT), 然而毫秒脈沖星觀測(cè)及數(shù)據(jù)處理中受到各種誤差的影響, 如原子鐘誤差、行星歷表精度、星際介質(zhì)不確定性、脈沖星自轉(zhuǎn)不穩(wěn)性(周期躍變及計(jì)時(shí)噪聲)、設(shè)備觀測(cè)誤差等, 這些誤差會(huì)使脈沖星觀測(cè)數(shù)據(jù)中包含高頻噪聲、高斯白噪聲、紅噪聲等, 往往這些噪聲影響是耦合的.除脈沖星自轉(zhuǎn)不穩(wěn)定性和太陽系行星歷表誤差外, 其他誤差源對(duì)不同脈沖星觀測(cè)的影響是獨(dú)立的, 可通過多顆毫秒脈沖星定義的綜合脈沖星時(shí), 削弱這些噪聲源的影響, 提高脈沖星時(shí)的穩(wěn)定性. 此外, 也可使用譜分析方法, 進(jìn)一步降低脈沖星計(jì)時(shí)殘差中的非高斯噪聲影響, 提高綜合脈沖星時(shí)的長期穩(wěn)定度. 本文研究了一種基于雙譜濾波的綜合脈沖星時(shí)構(gòu)建算法, 并利用IPTA 4顆觀測(cè)時(shí)間超過10 yr的毫秒脈沖星數(shù)據(jù), 來驗(yàn)證該方法的有效性.

    2 綜合脈沖星時(shí)計(jì)時(shí)及雙譜濾波基本原理

    毫秒脈沖星自轉(zhuǎn)極其穩(wěn)定, 其周期變化率在10?20s/s量級(jí), 1 d周期變化累積量為10?14s. 脈沖星自轉(zhuǎn)具有穩(wěn)定性、可觀測(cè)性, 故毫秒脈沖星計(jì)時(shí)觀測(cè)可構(gòu)建基于脈沖星自轉(zhuǎn)頻率的時(shí)間尺度, 稱為脈沖星時(shí). 脈沖星時(shí)是一種相對(duì)的時(shí)間尺度, 主要通過脈沖星時(shí)與參考原子時(shí)間的鐘差測(cè)量而實(shí)現(xiàn). 對(duì)毫秒脈沖星進(jìn)行高精度計(jì)時(shí)觀測(cè), 可得到脈沖星的計(jì)時(shí)殘差, 計(jì)時(shí)殘差是太陽系質(zhì)心處脈沖到達(dá)時(shí)間預(yù)報(bào)值和觀測(cè)值之差. 預(yù)報(bào)值是基于反映脈沖星自轉(zhuǎn)規(guī)律的脈沖相位模型(也稱鐘模型), 代表了脈沖星時(shí), 而通過射電觀測(cè)得到脈沖到達(dá)時(shí)間是以原子時(shí)(AT)為參考的, 故計(jì)時(shí)殘差含有脈沖星時(shí)與參考原子時(shí)之間的鐘差. 如果脈沖星計(jì)時(shí)處理中, 各種效應(yīng)得到完全修正, 那么計(jì)時(shí)殘差即為PT與AT之差. 脈沖信號(hào)到達(dá)測(cè)站的時(shí)刻是通過測(cè)站氫鐘記錄的, 之后通過鐘差修正溯源到TAI. 當(dāng)前, BIPM提供了TAI的兩種實(shí)現(xiàn)方式—準(zhǔn)實(shí)時(shí)TT (TAI)和事后TT(BIPM).

    隨著大量高精度原子鐘的加入、新型基準(zhǔn)鐘的涌現(xiàn)以及原子時(shí)算法的改進(jìn), TAI的準(zhǔn)確度和穩(wěn)定度不斷提升.然而考慮實(shí)時(shí)性的需求,BIPM每月計(jì)算并發(fā)表一個(gè)版本TAI,稱為TT (TAI). 由于兼顧修正模型優(yōu)化和計(jì)算效率, TT (TAI)精度會(huì)受到一定程度上的限制, 故其不是一個(gè)最高精度的時(shí)間尺度. 同時(shí), BIPM提供了另一種更高精度的時(shí)間尺度, 即TT (BIPM). TT (BIPM)每年計(jì)算一個(gè)版本, 利用了全球所有基準(zhǔn)鐘數(shù)據(jù)計(jì)算重新得到TT的精確實(shí)現(xiàn), 并在計(jì)算時(shí)段內(nèi)對(duì)所采用的基準(zhǔn)鐘加入了黑體輻射、相對(duì)論效應(yīng)等更多改正, 且對(duì)頻率進(jìn)行平滑和插值后得到, 使得TT (BIPM)比TT (TAI)更準(zhǔn)確和更穩(wěn)定. 本文選擇了TT (BIPM2015)作為脈沖星觀測(cè)數(shù)據(jù)處理的時(shí)間基準(zhǔn).

    綜合脈沖星時(shí)定義為多個(gè)脈沖星計(jì)時(shí)殘差的加權(quán)平均:

    式中PTi為第i個(gè)脈沖星時(shí),ωi為權(quán)重, 以脈沖星觀測(cè)誤差σi作為計(jì)算依據(jù).

    式中n為脈沖星的數(shù)量. 對(duì)于脈沖星時(shí)間穩(wěn)定度, Matsakis等人提出了σz估計(jì)方法進(jìn)行評(píng)估[1]. 假設(shè)脈沖星的天體測(cè)量參數(shù)、脈沖相位模型參數(shù)及雙星模型參數(shù)都得到精確測(cè)定, 暫不考慮星際介質(zhì)延遲、計(jì)時(shí)噪聲、行星歷表、參考原子時(shí)等誤差, 那么扣除二次多項(xiàng)式后的殘差就是脈沖星時(shí)與參考原子時(shí)剩余的最低階偏差, 因此σz適用于脈沖星時(shí)間穩(wěn)定度的估計(jì)[1].σz估計(jì)方法計(jì)算基本思路如下, 其將時(shí)間殘差測(cè)量序列分成m個(gè)子序列(m=1,2,4,8,···), 利用每一個(gè)子序列的殘差測(cè)量數(shù)據(jù)進(jìn)行最小二乘擬合, 使

    式中χ(tj)為觀測(cè)時(shí)間tj的計(jì)時(shí)殘差,ψj為其誤差,j為觀測(cè)的子序列號(hào),m為子序列個(gè)數(shù).X(tj)的表達(dá)式為:

    式中ck(k=0,1,2,3)為其系數(shù),t0為參考?xì)v元.σz定義如下:

    式中τ為子序列連續(xù)采樣時(shí)間,??表示對(duì)所有子序列的c3進(jìn)行加權(quán)平均, 權(quán)重類似于(2)式的計(jì)算.σz估計(jì)的上下限采用統(tǒng)計(jì)方法估計(jì), 詳見參考文獻(xiàn)[1].

    雙譜濾波是通過計(jì)算計(jì)時(shí)殘差的雙譜, 利用雙譜頻域?yàn)V波器對(duì)雙譜幅值和相位進(jìn)行濾波處理, 消除信號(hào)噪聲, 再對(duì)信號(hào)進(jìn)行重構(gòu)獲取濾波后的信號(hào). 非參數(shù)直接法雙譜估計(jì)的方差較小, 方法簡(jiǎn)單且容易實(shí)現(xiàn), 所以在信號(hào)處理中得到了較好的應(yīng)用, 直接法雙譜估計(jì)的原理請(qǐng)參考文獻(xiàn)[21]. 雙譜對(duì)高斯過程是“盲”的, 因此通過雙譜濾波可以很好地抑制信號(hào)中的高斯噪聲, 同時(shí)在雙譜域結(jié)合閾值濾波器, 也可有效地抑制加性噪聲、乘性噪聲等一些非高斯噪聲成分.

    3 算例分析

    射電望遠(yuǎn)鏡都是多科學(xué)任務(wù)的實(shí)體, 脈沖星計(jì)時(shí)觀測(cè)無法做到等間隔觀測(cè). 從IPTA發(fā)布的最新數(shù)據(jù)來看, 不同脈沖星觀測(cè)數(shù)據(jù)起始觀測(cè)時(shí)間不一樣, 觀測(cè)間隔非常不均勻,有些脈沖星間隔1 yr才有觀測(cè), 特別是早期數(shù)據(jù), 這種現(xiàn)象尤為明顯. 隨著技術(shù)的發(fā)展,脈沖星觀測(cè)精度提高, 計(jì)時(shí)殘差誤差減小. 綜合脈沖星時(shí)的計(jì)算及估計(jì)流程如下:

    (1)計(jì)算原始數(shù)據(jù)的綜合脈沖星時(shí). 對(duì)所有觀測(cè)序列去重, 并按照時(shí)間排序. 以天為單位, 對(duì)位于該天所有觀測(cè)系列進(jìn)行加權(quán)平均, 權(quán)重計(jì)算見(2)式. 如果該天沒有觀測(cè), 不均勻化處理, 設(shè)置為無數(shù)據(jù), 避免引入未知插值誤差. 得到新的計(jì)時(shí)殘差系列, 即綜合脈沖星時(shí)與參考原子時(shí)之差;

    (2)計(jì)算脈沖星時(shí)的σz. 按照(3)–(5)式計(jì)算綜合脈沖星時(shí)的穩(wěn)定度, 同時(shí)考慮兩個(gè)約束條件: 每個(gè)計(jì)時(shí)殘差子序列中觀測(cè)量數(shù)大于4; 子序列中時(shí)間間隔要大于子序列時(shí)間長度的同時(shí)對(duì)于m=1的情況, 采用最新觀測(cè)序列的計(jì)時(shí)殘差;

    (3)使用雙譜濾波對(duì)每個(gè)脈沖星的觀測(cè)序列進(jìn)行處理, 消除高斯噪聲的影響, 然后重復(fù)(1)–(2)步, 得到濾波后綜合脈沖星時(shí)的穩(wěn)定度.

    2019年IPTA發(fā)布了第2批共65顆脈沖星的觀測(cè)數(shù)據(jù)[22], 包括兩種類型的數(shù)據(jù).第1種類型數(shù)據(jù)繼承了第1批IPTA數(shù)據(jù)的形式, 擬合了色散量(Dispersion Measure, DM)和白噪聲參數(shù)化. 第2種數(shù)據(jù)在第1類基礎(chǔ)上進(jìn)行紅噪聲擬合. 本文選擇了比較穩(wěn)定的PSR J0437?4715、J0613?0200、J1713+0747和J1909?37444顆脈沖星的第1類數(shù)據(jù), 沒有處理紅噪聲. 4顆脈沖星的基本情況見表1, 擬合后殘差見圖1, 計(jì)時(shí)殘差的時(shí)間參考為TT (BIPM2015), 歷表采用DE436行星歷表. 表1中PSR為脈沖星名稱, Period為毫秒脈沖星的自轉(zhuǎn)頻率, RMS為計(jì)時(shí)殘差(residual, Res)的均方根值, Span為觀測(cè)周期長度, NTOA為參與解算的觀測(cè)量數(shù)量, Start和Finish為觀測(cè)數(shù)據(jù)的起始、結(jié)束時(shí)刻, 用簡(jiǎn)化儒略日(Modified Julian Date, MJD). 由圖1可知, 由于觀測(cè)設(shè)備更新和時(shí)間系統(tǒng)提升, 脈沖星觀測(cè)精度也相應(yīng)地提高.

    表1 IPTA 4顆脈沖星的基本信息Table 1 The basic information of four IPTA pulsars

    圖1 IPTA 4顆脈沖星的計(jì)時(shí)殘差(時(shí)間基準(zhǔn)TT (BIPM2015))Fig.1 Timing residuals of four IPTA Pulsars (time reference is TT (BIPM2015))

    脈沖星觀測(cè)及數(shù)據(jù)處理中, 不可避免地受到各種誤差源的影響, 如各種殘余的微小射電干擾、氫鐘的不穩(wěn)定性及設(shè)備的熱噪聲、數(shù)據(jù)處理中DM偏差、行星歷表原點(diǎn)誤差、計(jì)時(shí)噪聲難以白化處理等. 這些影響會(huì)在計(jì)時(shí)殘差中產(chǎn)生高斯及非高斯噪聲, 而雙譜濾波可抑制其中部分噪聲. 每顆脈沖星的計(jì)時(shí)殘差包含的噪聲成分不盡相同, 因此表現(xiàn)出不同殘差形態(tài), 雙譜濾波方法可能對(duì)不同脈沖星的處理效果不盡相同. 作為對(duì)比,按照上節(jié)計(jì)算過程, 采用原始和濾波后殘差數(shù)據(jù)分別構(gòu)建4顆脈沖星綜合時(shí), 見圖2.

    脈沖星觀測(cè)具有獨(dú)立性, 不具有時(shí)間關(guān)聯(lián)性, 故在數(shù)據(jù)處理中不建議均勻插值采樣,那樣做盡管可填補(bǔ)空白時(shí)間段的觀測(cè)數(shù)據(jù), 但插值沒有真實(shí)觀測(cè)量對(duì)應(yīng), 容易白化和平滑各種噪聲, 導(dǎo)致綜合脈沖星時(shí)評(píng)估精度過高. 由圖1和圖2可知, 基于原始數(shù)據(jù)得到綜合脈沖星時(shí), 實(shí)際上是4顆脈沖星時(shí)的加權(quán)平均, 綜合脈沖星時(shí)精度主要受觀測(cè)精度略差的PSR J0613?0200、J1713+0747的影響. 基于雙譜濾波數(shù)據(jù)的綜合脈沖星時(shí)精度得到明顯的提升, 計(jì)時(shí)殘差的平均值從5.69×10?7s減小到1.12×10?7s, 方差從1.25×10?6s降低到9.58×10?8s, 精度有大幅度提高. 使用σz評(píng)估兩種方法的綜合脈沖星時(shí)穩(wěn)定性, 分別見圖3、圖4.

    圖2 IPTA 4顆脈沖星構(gòu)建的綜合脈沖星時(shí)(時(shí)間基準(zhǔn)參考于TT (BIPM2015), 上圖基于原始數(shù)據(jù), 下圖基于濾波后數(shù)據(jù)).Fig.2 Ensemble pulsar time of four IPTA Pulsars (time reference is TT (BIPM2015), the top panel is based on the original data, and the bottom panel is based on the filtered data).

    圖3 基于IPTA 4顆脈沖星原始數(shù)據(jù)構(gòu)建的綜合脈沖星時(shí)的穩(wěn)定度Fig.3 The stability of ensemble pulsar time established by the original data of four IPTA Pulsars

    由圖3可知, 4顆脈沖星J0437?4715、J0613?0200、J1713+0747、J1909?3744和綜合脈沖星時(shí)的月穩(wěn)定度分別為2.26×10?13、2.37×10?11、1.03×10?9、2.55×10?13、3.77×10?13; 年穩(wěn)定度分別為3.30×10?14、1.94×10?13、7.32×10?14、1.65×10?14、7.77×10?14; 10 yr的穩(wěn)定度為1.23×10?15、2.25×10?14、2.86×10?15、8.75×10?16、8.56×10?16. 基于所有濾波數(shù)據(jù)實(shí)現(xiàn)的綜合脈沖星時(shí)穩(wěn)定度為2.03×10?16, 可見4顆脈沖星時(shí)與綜合脈沖星時(shí)的穩(wěn)定度隨著觀測(cè)時(shí)間增加而更加穩(wěn)定, 穩(wěn)定度較好的兩顆脈沖星時(shí)是觀測(cè)精度高的J0437?4715和J1909?3744. 當(dāng)觀測(cè)時(shí)間超過10 yr,綜合脈沖星時(shí)的穩(wěn)定性更優(yōu), 原因在于觀測(cè)數(shù)據(jù)充分, 且有大批高質(zhì)量的新觀測(cè)數(shù)據(jù)加入. 由圖4可知, 基于雙譜濾波算法構(gòu)建4顆脈沖星時(shí)和綜合脈沖星時(shí), 其穩(wěn)定度的基本規(guī)律與基于原始數(shù)據(jù)構(gòu)建的脈沖星時(shí)一樣, 但穩(wěn)定度更高. 4顆脈沖星J0437?4715、J0613?0200、J1713+0747、J1909?3744和綜合脈沖星時(shí)的月穩(wěn)定度分別為3.77×10?14、1.60×10?12、6.75×10?11、4.76×10?14、7.40×10?14; 年穩(wěn)定度分別為6.94×10?15、4.39×10?14、2.03×10?14、3.19×10?15、1.50×10?14; 10 yr的穩(wěn)定度為2.91×10?15、1.74×10?14、6.53×10?16、7.36×10?17、3.50×10?16. 基于所有濾波數(shù)據(jù)實(shí)現(xiàn)的綜合脈沖星時(shí)穩(wěn)定度為8.21×10?17. 相比于經(jīng)典加權(quán)算法, 基于雙譜濾波算法的綜合脈沖星時(shí)構(gòu)建算法, 能夠較好地提升單脈沖星時(shí)及綜合脈沖星時(shí)的穩(wěn)定性. 需要說明的是, 當(dāng)時(shí)間尺度大于數(shù)據(jù)觀測(cè)周期的一半時(shí), 即m=1, 會(huì)優(yōu)先使用最新的觀測(cè)數(shù)據(jù), 這樣的數(shù)據(jù)處理策略會(huì)導(dǎo)致脈沖星時(shí)的穩(wěn)定度較好.

    圖4 基于IPTA 4顆脈沖星濾波數(shù)據(jù)構(gòu)建的綜合脈沖星時(shí)的穩(wěn)定度Fig.4 The stability of the ensemble pulsar time established by the filtering data of four IPTA Pulsars

    BIPM定期發(fā)布各實(shí)驗(yàn)室保持的地方原子時(shí)TAI (XX)與國際原子時(shí)TAI的偏差文件, XX代表具體實(shí)驗(yàn)室, 用英文名稱簡(jiǎn)寫表示. 同樣使用σz方法估算了中科院國家授時(shí)中心(National Time Service Center, NTSC)、德國技術(shù)物理研究所(Physikalisch-Technische Bundesanstalt, PTB)、俄羅斯國家技術(shù)物理及無線電工程研究院(VNIIFTRI, 簡(jiǎn)寫SU)、美國國家標(biāo)準(zhǔn)技術(shù)研究院(National Institute of Standards and Technology, NIST) 4家授時(shí)單位原子時(shí)的穩(wěn)定性, 時(shí)間從MJD50814—MJD58904. PTB是目前國際上所有授時(shí)單位中原子鐘精度最高的單位之一, 在國際綜合原子時(shí)建立中的權(quán)重較大, TAI (PTB)的數(shù)據(jù)由于基準(zhǔn)鐘升級(jí)存在明顯的2次跳變, 且第2段數(shù)據(jù)較短, 故本文分別處理了第1段和第3段數(shù)據(jù), 第3段數(shù)據(jù)精度明顯優(yōu)于第1段, 下面分別標(biāo)記為TAI (PTB1)、TAI (PTB2), 第1段數(shù)據(jù)周期為MJD50814—MJD55439, 第2段數(shù)據(jù)周期為MJD56079—MJD58904. 4家單位原子鐘穩(wěn)定性見圖5.

    原子時(shí)的穩(wěn)定度一般使用Allan方差表示, 由于方便與綜合脈沖星(Ensemble Pulsar Time, EPT)的比較, 也使用σz方法評(píng)估其穩(wěn)定度. 由圖5可知,σz方法表示的原子時(shí)穩(wěn)定度基本上保持穩(wěn)定, 首先隨著觀測(cè)時(shí)間原子時(shí)穩(wěn)定度緩慢提高, 約6 yr后穩(wěn)定性下降.TAI (NIST)、TAI (SU)、TAI (NTSC)、TAI (PTB1)和TAI (PTB2)的年穩(wěn)定度分別為1.19×10?12、1.38×10?13、3.53×10?14、2.31×10?13、1.46×10?15; 5 yr穩(wěn)定度分別為1.05×10?13、1.12×10?14、2.51×10?15、5.51×10?14、2.22×10?16. 綜合脈沖星時(shí)的年穩(wěn)定度優(yōu)于TA(NIST)、TA(SU)、TAI(NTSC)、TAI(PTB1), 遜色于TAI(PTB2),TAI (PTB2)代表當(dāng)前原子鐘的最高水平. 下面以TAI (PTB2)為例進(jìn)行分析, 隨著觀測(cè)精度的提高, 綜合脈沖星時(shí)穩(wěn)定性提升, 而原子時(shí)穩(wěn)定性恒定, 5 yr時(shí)間兩者穩(wěn)定性相當(dāng),然而5 yr后綜合脈沖星時(shí)優(yōu)于TAI (PTB2), 可用于改善原子時(shí)穩(wěn)定性. 脈沖星時(shí)短穩(wěn)主要受頻率閃爍噪聲和隨機(jī)游走頻率噪聲等系統(tǒng)噪聲影響, 主要來源包括脈沖星自轉(zhuǎn)的不規(guī)則性、觀測(cè)系統(tǒng)噪聲或是星際閃爍噪聲. 脈沖星時(shí)長穩(wěn)主要受脈沖星自轉(zhuǎn)不規(guī)律性的影響, 主要包括計(jì)時(shí)噪聲, 故精確的計(jì)時(shí)噪聲模型、優(yōu)秀的觀測(cè)設(shè)備及更合理的觀測(cè)方案有利于脈沖星時(shí)穩(wěn)定度的提升.

    圖5 國際上4家授時(shí)單位原子時(shí)的穩(wěn)定度Fig.5 Stability of atomic time of four international time service units

    4 結(jié)論與分析

    本文研究了一種基于雙譜濾波的綜合脈沖星時(shí)算法, 并利用IPTA的4顆毫秒脈沖星觀測(cè)數(shù)據(jù)展示其良好性能, 比較分析了脈沖星時(shí)與原子時(shí)的穩(wěn)定性. 由于脈沖星的動(dòng)能會(huì)因自身的輻射和周圍吸積的旋轉(zhuǎn)物質(zhì)同磁層間的相互作用這兩個(gè)因素逐漸消耗, 使脈沖星旋轉(zhuǎn)角動(dòng)量減少, 所以脈沖星自轉(zhuǎn)周期都是隨著時(shí)間逐漸變長, 然而也存在一類自轉(zhuǎn)不穩(wěn)定性現(xiàn)象, 如計(jì)時(shí)噪聲. 計(jì)時(shí)噪聲是指脈沖星自轉(zhuǎn)參數(shù)發(fā)生的連續(xù)、時(shí)標(biāo)較長(通常為幾個(gè)月或幾年)的擾動(dòng). 一般來說, 計(jì)時(shí)噪聲表現(xiàn)為無規(guī)律的隨機(jī)信號(hào), 是一種均值為零的高斯分布, 該類噪聲可使用雙譜濾波較好地消除. 然而不少脈沖星的計(jì)時(shí)噪聲表現(xiàn)為準(zhǔn)周期性、紅噪聲, 其功率在低頻端更強(qiáng). 當(dāng)前計(jì)時(shí)噪聲產(chǎn)生原因不明確, 普遍認(rèn)為其與中子星內(nèi)部的超流過程、內(nèi)部溫度的變化以及磁層中的物理過程有關(guān). 下一步可對(duì)紅噪聲進(jìn)行建模擬合, 如采用Cholesky方法和貝葉斯估計(jì)方法對(duì)計(jì)時(shí)噪聲進(jìn)行估計(jì),減弱紅噪聲的影響.

    本文的研究也可為我國脈沖星時(shí)間系統(tǒng)的建設(shè)及應(yīng)用提供參考. 我國新疆天文臺(tái)南山射電望遠(yuǎn)鏡于1996年率先在國內(nèi)開展脈沖星觀測(cè), 對(duì)脈沖星已經(jīng)連續(xù)觀測(cè)20多年, 同時(shí)我國于2016年建成的全球最大的射電望遠(yuǎn)鏡500 m口徑球面射電望遠(yuǎn)鏡(Fivehundred-meter Aperture Spherical radio Telescope, FAST)能實(shí)現(xiàn)當(dāng)前毫秒脈沖星最高精度的觀測(cè), 它們都將為我國脈沖星時(shí)間尺度的研究提供寶貴的資料.

    猜你喜歡
    雙譜脈沖星穩(wěn)定度
    “中國天眼”已發(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
    高穩(wěn)晶振短期頻率穩(wěn)定度的仿真分析
    雙譜圖在語音分析中的應(yīng)用
    基于雙譜特征融合的通信輻射源識(shí)別算法
    基于虛擬觀測(cè)值的X射線單脈沖星星光組合導(dǎo)航
    多MOSFET并聯(lián)均流的高穩(wěn)定度恒流源研究
    工藝參數(shù)對(duì)橡膠球鉸徑向剛度穩(wěn)定度的影響
    基于小波包域雙譜的風(fēng)力機(jī)振動(dòng)信號(hào)監(jiān)測(cè)方法研究
    旱澇不穩(wěn)定度的定量化研究及應(yīng)用
    国产男女超爽视频在线观看| 午夜福利在线观看吧| 精品久久久久久久久亚洲| 国产av码专区亚洲av| 亚洲国产成人一精品久久久| or卡值多少钱| 3wmmmm亚洲av在线观看| 欧美潮喷喷水| 午夜激情福利司机影院| 激情 狠狠 欧美| 免费观看a级毛片全部| 老司机影院毛片| 黄色欧美视频在线观看| 亚洲成人一二三区av| 日本一二三区视频观看| 成人漫画全彩无遮挡| 婷婷色麻豆天堂久久| 国产av码专区亚洲av| 精品国产一区二区三区久久久樱花 | 亚洲欧美一区二区三区黑人 | 久久国产乱子免费精品| 麻豆成人av视频| 亚洲一区高清亚洲精品| 卡戴珊不雅视频在线播放| 草草在线视频免费看| 在线天堂最新版资源| 国产色爽女视频免费观看| 国产精品美女特级片免费视频播放器| 18禁在线播放成人免费| 热99在线观看视频| 麻豆久久精品国产亚洲av| av在线亚洲专区| 久久久久久久久久久丰满| 亚洲欧美中文字幕日韩二区| 一级毛片久久久久久久久女| 亚洲国产精品专区欧美| 日韩欧美 国产精品| 亚洲在久久综合| 极品少妇高潮喷水抽搐| 亚洲人与动物交配视频| av在线播放精品| av又黄又爽大尺度在线免费看| 五月伊人婷婷丁香| 国产免费福利视频在线观看| 色5月婷婷丁香| 日韩电影二区| 日韩欧美一区视频在线观看 | 晚上一个人看的免费电影| 成年人午夜在线观看视频 | 国产精品熟女久久久久浪| 欧美成人a在线观看| 九九久久精品国产亚洲av麻豆| 一级二级三级毛片免费看| 大香蕉97超碰在线| 亚洲国产成人一精品久久久| 日日摸夜夜添夜夜添av毛片| 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 久久综合国产亚洲精品| 在线免费十八禁| 精品久久久久久久久久久久久| 国产伦在线观看视频一区| 我的老师免费观看完整版| 纵有疾风起免费观看全集完整版 | 亚洲第一区二区三区不卡| 精品少妇黑人巨大在线播放| 欧美一级a爱片免费观看看| 18禁动态无遮挡网站| 在线免费观看不下载黄p国产| 性色avwww在线观看| 校园人妻丝袜中文字幕| 亚洲av.av天堂| 黄片无遮挡物在线观看| 少妇丰满av| 亚洲激情五月婷婷啪啪| www.av在线官网国产| 色5月婷婷丁香| 赤兔流量卡办理| 日韩欧美国产在线观看| 欧美3d第一页| 亚洲av免费高清在线观看| 精品久久久久久久末码| 偷拍熟女少妇极品色| 中国国产av一级| 少妇的逼好多水| 美女大奶头视频| 日韩人妻高清精品专区| 国产高清不卡午夜福利| 国产精品久久久久久久电影| 一区二区三区高清视频在线| a级一级毛片免费在线观看| 成人亚洲精品av一区二区| 中文乱码字字幕精品一区二区三区 | 有码 亚洲区| 国产精品国产三级国产av玫瑰| 亚洲成人中文字幕在线播放| 亚洲精品日韩在线中文字幕| 免费观看在线日韩| 九九在线视频观看精品| 看免费成人av毛片| av.在线天堂| 熟女电影av网| 精品国产露脸久久av麻豆 | 国内精品一区二区在线观看| 日韩电影二区| videos熟女内射| 女的被弄到高潮叫床怎么办| 久久精品夜色国产| 亚洲最大成人手机在线| 免费观看无遮挡的男女| 青春草亚洲视频在线观看| 国产亚洲91精品色在线| www.色视频.com| 亚洲美女视频黄频| 久久精品国产自在天天线| 久久精品久久久久久噜噜老黄| 婷婷六月久久综合丁香| 国产伦精品一区二区三区视频9| 亚洲一级一片aⅴ在线观看| 丰满人妻一区二区三区视频av| 激情五月婷婷亚洲| 噜噜噜噜噜久久久久久91| 亚洲内射少妇av| 九色成人免费人妻av| 激情 狠狠 欧美| 一个人观看的视频www高清免费观看| 国产综合精华液| 少妇的逼好多水| 色吧在线观看| 国产精品久久久久久精品电影| 国内揄拍国产精品人妻在线| 日韩制服骚丝袜av| 亚洲欧洲国产日韩| 日韩av不卡免费在线播放| 永久网站在线| 插逼视频在线观看| 寂寞人妻少妇视频99o| 麻豆久久精品国产亚洲av| 全区人妻精品视频| freevideosex欧美| 婷婷色综合大香蕉| 国产黄频视频在线观看| 蜜臀久久99精品久久宅男| 汤姆久久久久久久影院中文字幕 | 国产一区二区三区av在线| 色综合色国产| kizo精华| 18+在线观看网站| 日韩在线高清观看一区二区三区| 国产综合懂色| 欧美一级a爱片免费观看看| 欧美潮喷喷水| 国内精品一区二区在线观看| 人人妻人人看人人澡| 免费黄频网站在线观看国产| 日日啪夜夜爽| 精品久久久久久成人av| 99热这里只有精品一区| 99久久人妻综合| 成人综合一区亚洲| 久久久久久久久中文| 天天一区二区日本电影三级| 看黄色毛片网站| 99久久精品一区二区三区| 日韩电影二区| 久久这里有精品视频免费| 国产伦精品一区二区三区四那| 欧美潮喷喷水| 成人欧美大片| 国产激情偷乱视频一区二区| 一级毛片aaaaaa免费看小| 大又大粗又爽又黄少妇毛片口| 成人美女网站在线观看视频| 观看免费一级毛片| 真实男女啪啪啪动态图| 最近最新中文字幕大全电影3| 亚洲天堂国产精品一区在线| 在线观看人妻少妇| 最近最新中文字幕大全电影3| 免费黄频网站在线观看国产| 日本一本二区三区精品| 亚洲人成网站高清观看| 久久久久久久亚洲中文字幕| 日本一本二区三区精品| 国产成人a区在线观看| 精品国内亚洲2022精品成人| av专区在线播放| 亚洲欧美日韩东京热| 国产成人午夜福利电影在线观看| 国产精品一区二区三区四区免费观看| 在线天堂最新版资源| 免费看a级黄色片| 一级毛片 在线播放| xxx大片免费视频| 嫩草影院新地址| 国产在视频线精品| 亚洲va在线va天堂va国产| 日韩 亚洲 欧美在线| 成人亚洲精品av一区二区| 欧美变态另类bdsm刘玥| 国产精品一区二区在线观看99 | 麻豆久久精品国产亚洲av| 国产爱豆传媒在线观看| 自拍偷自拍亚洲精品老妇| 精品久久久久久久久亚洲| 国产成人福利小说| 偷拍熟女少妇极品色| 国模一区二区三区四区视频| 国产乱来视频区| av线在线观看网站| 极品少妇高潮喷水抽搐| 国产高清三级在线| 成人午夜精彩视频在线观看| 亚洲国产高清在线一区二区三| 亚洲av中文av极速乱| 亚洲乱码一区二区免费版| 日韩三级伦理在线观看| 日韩成人av中文字幕在线观看| 亚洲精品影视一区二区三区av| 亚洲国产精品成人久久小说| 99热6这里只有精品| 18禁动态无遮挡网站| 精品久久久久久久久久久久久| 中文字幕av在线有码专区| 黄色日韩在线| 免费看a级黄色片| 亚洲欧美成人综合另类久久久| 男人舔女人下体高潮全视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 又爽又黄a免费视频| 赤兔流量卡办理| 麻豆成人av视频| 欧美不卡视频在线免费观看| 肉色欧美久久久久久久蜜桃 | 国产男女超爽视频在线观看| 永久免费av网站大全| 美女主播在线视频| 欧美精品国产亚洲| 69人妻影院| 国产成人a∨麻豆精品| 久久久久网色| 99热这里只有是精品在线观看| 精品久久久久久成人av| 国产一区有黄有色的免费视频 | 欧美日韩一区二区视频在线观看视频在线 | 欧美日韩国产mv在线观看视频 | 国产白丝娇喘喷水9色精品| av.在线天堂| 午夜福利在线在线| 国产精品av视频在线免费观看| 777米奇影视久久| 欧美性猛交╳xxx乱大交人| 大片免费播放器 马上看| 亚洲成人久久爱视频| 18禁动态无遮挡网站| 亚洲天堂国产精品一区在线| 国产伦一二天堂av在线观看| 国产免费福利视频在线观看| 久久精品国产自在天天线| 日韩强制内射视频| 日本欧美国产在线视频| 精品久久久久久久久久久久久| 少妇熟女aⅴ在线视频| 日韩电影二区| 三级男女做爰猛烈吃奶摸视频| 91aial.com中文字幕在线观看| 亚洲欧美成人综合另类久久久| 亚洲av福利一区| 日日干狠狠操夜夜爽| 水蜜桃什么品种好| 特大巨黑吊av在线直播| 免费看不卡的av| 寂寞人妻少妇视频99o| 男人舔奶头视频| 日本午夜av视频| 国产 亚洲一区二区三区 | 亚洲自偷自拍三级| 日日摸夜夜添夜夜爱| 日韩一区二区三区影片| 又爽又黄a免费视频| 看免费成人av毛片| 深夜a级毛片| 国产成年人精品一区二区| 两个人视频免费观看高清| 亚洲国产精品成人久久小说| 免费黄频网站在线观看国产| 搞女人的毛片| 两个人的视频大全免费| 国内揄拍国产精品人妻在线| 国产成人a∨麻豆精品| 亚洲av电影不卡..在线观看| 一个人看视频在线观看www免费| 国产成人午夜福利电影在线观看| 99热这里只有是精品50| 99久久人妻综合| 国产伦精品一区二区三区四那| 欧美三级亚洲精品| 日本wwww免费看| 亚洲内射少妇av| 一区二区三区免费毛片| 欧美3d第一页| 高清视频免费观看一区二区 | 精品一区在线观看国产| 国模一区二区三区四区视频| 一级av片app| 看黄色毛片网站| 日本黄大片高清| 婷婷六月久久综合丁香| 两个人视频免费观看高清| 99久国产av精品| 国产伦理片在线播放av一区| 免费观看在线日韩| 七月丁香在线播放| 天美传媒精品一区二区| 乱系列少妇在线播放| videos熟女内射| www.av在线官网国产| 日日啪夜夜撸| 又粗又硬又长又爽又黄的视频| 性插视频无遮挡在线免费观看| 麻豆乱淫一区二区| 99热这里只有是精品50| 黄色日韩在线| 精华霜和精华液先用哪个| 久久精品综合一区二区三区| 成人无遮挡网站| 极品教师在线视频| 日日摸夜夜添夜夜添av毛片| 一区二区三区高清视频在线| 高清日韩中文字幕在线| 日日摸夜夜添夜夜爱| 麻豆av噜噜一区二区三区| 成人午夜精彩视频在线观看| 免费观看在线日韩| 亚洲怡红院男人天堂| 日本爱情动作片www.在线观看| 免费黄色在线免费观看| 日韩欧美一区视频在线观看 | 色综合色国产| 久久国产乱子免费精品| 欧美一区二区亚洲| 超碰97精品在线观看| 熟妇人妻久久中文字幕3abv| 97人妻精品一区二区三区麻豆| av在线播放精品| 亚洲一区高清亚洲精品| 80岁老熟妇乱子伦牲交| 亚洲欧美精品自产自拍| 精品午夜福利在线看| 亚洲四区av| 成年女人在线观看亚洲视频 | 97超碰精品成人国产| 精品一区二区三区人妻视频| xxx大片免费视频| 毛片女人毛片| 日本一本二区三区精品| 免费观看在线日韩| 69人妻影院| 日韩 亚洲 欧美在线| 日日啪夜夜爽| 黄色日韩在线| 免费av毛片视频| 欧美xxxx黑人xx丫x性爽| 国产 一区 欧美 日韩| 亚洲欧美一区二区三区黑人 | 成年免费大片在线观看| 国产日韩欧美在线精品| 国产女主播在线喷水免费视频网站 | 亚洲av电影在线观看一区二区三区 | 街头女战士在线观看网站| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品电影小说 | 精品久久久噜噜| 天堂影院成人在线观看| 三级经典国产精品| 菩萨蛮人人尽说江南好唐韦庄| 最近手机中文字幕大全| 伦精品一区二区三区| 国产精品国产三级国产av玫瑰| 国产精品一区www在线观看| 成人亚洲欧美一区二区av| 免费看美女性在线毛片视频| 成人二区视频| 日韩国内少妇激情av| 国产一区二区在线观看日韩| 大片免费播放器 马上看| 老女人水多毛片| 欧美精品一区二区大全| 国产国拍精品亚洲av在线观看| 男人和女人高潮做爰伦理| 中文字幕久久专区| 久久久久久九九精品二区国产| 日本猛色少妇xxxxx猛交久久| 熟女人妻精品中文字幕| 99久久精品一区二区三区| 中文字幕av成人在线电影| av一本久久久久| 秋霞在线观看毛片| 亚洲精品一二三| 国产精品无大码| 久久97久久精品| 丝瓜视频免费看黄片| 在线a可以看的网站| 最近最新中文字幕免费大全7| 亚洲精品国产av成人精品| 亚洲人成网站高清观看| 欧美精品一区二区大全| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av涩爱| 国产精品国产三级国产av玫瑰| 久久精品久久久久久久性| av福利片在线观看| xxx大片免费视频| 成年人午夜在线观看视频 | 亚洲性久久影院| 欧美bdsm另类| 18禁在线无遮挡免费观看视频| 又爽又黄a免费视频| 亚洲电影在线观看av| 九色成人免费人妻av| 免费在线观看成人毛片| 国产伦一二天堂av在线观看| 日韩中字成人| 日韩欧美精品免费久久| 欧美97在线视频| 亚洲精品日韩av片在线观看| 亚洲国产最新在线播放| 国产伦精品一区二区三区四那| 黄片wwwwww| 美女cb高潮喷水在线观看| 日日摸夜夜添夜夜添av毛片| 美女国产视频在线观看| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 久久久国产一区二区| 免费看日本二区| 1000部很黄的大片| 亚洲精品久久午夜乱码| 日日啪夜夜爽| 日韩欧美 国产精品| 女人十人毛片免费观看3o分钟| 精品人妻一区二区三区麻豆| 在线观看av片永久免费下载| 一区二区三区四区激情视频| 天堂网av新在线| 日韩精品青青久久久久久| 视频中文字幕在线观看| 日本三级黄在线观看| 纵有疾风起免费观看全集完整版 | 18禁在线无遮挡免费观看视频| 成人一区二区视频在线观看| 久久国产乱子免费精品| 在现免费观看毛片| 国产一区二区三区综合在线观看 | 久久精品久久久久久噜噜老黄| 国产伦理片在线播放av一区| 日日摸夜夜添夜夜爱| 国产精品麻豆人妻色哟哟久久 | 能在线免费看毛片的网站| 国产又色又爽无遮挡免| 国产一区亚洲一区在线观看| 国产一级毛片在线| 肉色欧美久久久久久久蜜桃 | 国产伦精品一区二区三区视频9| 99久国产av精品国产电影| 色5月婷婷丁香| 日日撸夜夜添| 国产视频内射| 欧美日韩一区二区视频在线观看视频在线 | 欧美高清性xxxxhd video| 国产免费福利视频在线观看| 一级毛片久久久久久久久女| 观看免费一级毛片| 免费观看a级毛片全部| 少妇人妻精品综合一区二区| 波野结衣二区三区在线| 久久久精品94久久精品| 亚洲天堂国产精品一区在线| 男人舔奶头视频| 国模一区二区三区四区视频| 欧美变态另类bdsm刘玥| 久久久久精品久久久久真实原创| 国产淫片久久久久久久久| 天天一区二区日本电影三级| 亚洲自偷自拍三级| 久久草成人影院| 人妻夜夜爽99麻豆av| 91在线精品国自产拍蜜月| 日韩视频在线欧美| 精品久久久久久久久亚洲| 在线 av 中文字幕| 国产成人freesex在线| 亚洲18禁久久av| 久久草成人影院| 1000部很黄的大片| 91久久精品国产一区二区成人| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 老司机影院成人| 欧美激情久久久久久爽电影| 国产成人福利小说| 国产久久久一区二区三区| 日韩精品青青久久久久久| 日韩欧美精品免费久久| av在线亚洲专区| 国产精品人妻久久久影院| 欧美成人精品欧美一级黄| 国产免费一级a男人的天堂| 91久久精品国产一区二区成人| 国产伦精品一区二区三区四那| 久久这里只有精品中国| 99热全是精品| 成年女人看的毛片在线观看| 69人妻影院| 欧美变态另类bdsm刘玥| 高清日韩中文字幕在线| 亚洲精品456在线播放app| 一区二区三区高清视频在线| 亚洲高清免费不卡视频| 97超碰精品成人国产| 日韩欧美三级三区| 中文字幕亚洲精品专区| 日韩av在线免费看完整版不卡| 性色avwww在线观看| 国产一区二区三区综合在线观看 | 久久人人爽人人爽人人片va| 亚洲人成网站在线播| 午夜福利在线在线| 亚洲欧美一区二区三区国产| 亚洲欧美精品自产自拍| 精品久久久久久久末码| 亚洲av二区三区四区| 亚洲高清免费不卡视频| 色网站视频免费| 成人无遮挡网站| 欧美不卡视频在线免费观看| 国产精品女同一区二区软件| 中文字幕av在线有码专区| 国产亚洲91精品色在线| 久久久a久久爽久久v久久| 亚洲国产av新网站| 高清毛片免费看| 欧美bdsm另类| 国产精品伦人一区二区| 亚洲精品亚洲一区二区| 国产精品一区www在线观看| 亚洲电影在线观看av| 成人无遮挡网站| h日本视频在线播放| 一夜夜www| 国产精品99久久久久久久久| 听说在线观看完整版免费高清| 99久久精品一区二区三区| 亚洲精品国产成人久久av| 中国美白少妇内射xxxbb| 欧美日韩视频高清一区二区三区二| 国产一区二区三区综合在线观看 | 午夜福利视频1000在线观看| 亚洲色图av天堂| 一区二区三区免费毛片| 欧美日韩综合久久久久久| 大香蕉97超碰在线| 男的添女的下面高潮视频| 日韩 亚洲 欧美在线| 久久99精品国语久久久| 欧美日韩亚洲高清精品| 99热网站在线观看| 国产在线一区二区三区精| 26uuu在线亚洲综合色| 中文字幕制服av| 男人舔女人下体高潮全视频| 蜜桃亚洲精品一区二区三区| 97人妻精品一区二区三区麻豆| 亚洲成人一二三区av| 国产午夜精品久久久久久一区二区三区| 99热6这里只有精品| 亚洲精品亚洲一区二区| 乱系列少妇在线播放| 日韩人妻高清精品专区| 国产成人aa在线观看| 日韩av不卡免费在线播放| 非洲黑人性xxxx精品又粗又长| 亚洲激情五月婷婷啪啪| 日韩 亚洲 欧美在线| freevideosex欧美| 国产精品一二三区在线看| 蜜臀久久99精品久久宅男| a级毛色黄片| 亚洲欧美日韩无卡精品| 超碰97精品在线观看| 日本色播在线视频| 99久国产av精品| 亚洲精品成人久久久久久| 狠狠精品人妻久久久久久综合| 少妇人妻一区二区三区视频| 韩国高清视频一区二区三区| av播播在线观看一区| 免费看av在线观看网站| 哪个播放器可以免费观看大片| 国产一级毛片在线| 哪个播放器可以免费观看大片| 婷婷色av中文字幕| 免费观看a级毛片全部| 欧美3d第一页| 深夜a级毛片| 精品久久久久久成人av| 婷婷色综合大香蕉| 女人被狂操c到高潮| 国产午夜精品久久久久久一区二区三区| 少妇人妻精品综合一区二区| 午夜日本视频在线| 韩国av在线不卡| 日韩欧美精品v在线| 国产精品久久久久久久久免|