J0613-0200驾驭铯钟的方法_李变.pdf
《J0613-0200驾驭铯钟的方法_李变.pdf》由会员分享,可在线阅读,更多相关《J0613-0200驾驭铯钟的方法_李变.pdf(8页珍藏版)》请在咨信网上搜索。
1、526500-1航空学报Acta Aeronautica et Astronautica SinicaFeb.15 2023 Vol.44 No.3ISSN 1000-6893 CN 11-1929/VJ0613-0200驾驭铯钟的方法李变1,2,*,屈俐俐1,2,高玉平1,2,31中国科学院 国家授时中心,西安 7106002中国科学院 时间频率基准重点实验室,西安 7106003中国科学院大学,北京 100049摘 要:近年来随着脉冲星观测设备性能提升,尤其是中国 FAST 的投入运行,毫秒脉冲星的计时观测精度不断提高,观测资料的数量、质量及时效性也显著改善,为脉冲星时的建立及应用提供有力
2、的数据及技术支持。脉冲星时与原子时具有优势互补的特点,利用中国科学院国家授时中心的铯钟和相关时频设备,结合 IPTA dr2数据资料,建立脉冲星驾驭铯钟的实验系统。并以 J0613-0200为参考,开展单颗毫秒脉冲星驾驭铯钟的研究。截止目前,实验结果符合单颗毫秒脉冲星驾驭铯钟的预期,证明了脉冲星驾驭原子钟的可行性,并且实验系统建立的综合时间尺度保持了铯钟稳定度的优势。关键词:毫秒脉冲星;铯钟;驾驭;计时噪声;稳定度中图分类号:V556;P16 文献标识码:A 文章编号:1000-6893(2023)03-526500-08目 前,精 度 最 高 的 时 间 尺 度 国 际 原 子 时(TAI)
3、、协调世界时(UTC)和地球时的另一个实现 TT(BIPMXX)都是由国际权度局(BIPM)基于原子钟建立并保持的。BIPM 首先收集分布在全球 80 余个守时实验室约 500 台原子钟的数据资料,采用 ALGOS 算法加权平均算法得到自由原 子 时(EAL),然 后 利 用 10 余 台 基 准 频 标 对EAL进行频率校准后获得 TAI。UTC是 TAI与世界时相互协调的产物,与 TAI 只差整数秒,截止目前二者相差 37 s1。参与 EAL 计算的原子钟主要有铯钟和氢钟2 种类型,分别约占总钟数的 59%和 40%,TT(BIPMXX)是 BIPM 每年年初利用所有可用的基准频标和原子钟
4、的数据资料,采用事后处理的方式获得的滞后一年的时间尺度,也是目前全球最 准 确、最 稳 定 的 原 子 时 尺 度2。然 而 TT(BIPMXX)、TAI和 UTC 是滞后的纸面时间,无法满足时间/频率作为参考实时性的要求。因此,各个国家或地区的守时实验室都建立与保持着具有实时、连续、稳定的物理信号 UTC(k),作为 UTC 的物理实现用于实时测量比对,k为实验室名称缩写。毫秒脉冲星通常被认为是由双星系统演化而来的,双星系统中的正常脉冲星通过吸积伴星质量来增加角动量,将脉冲星自转周期加速到毫秒量级。与年轻正常脉冲星相比,毫秒脉冲星具有辐射流量弱、自转稳定度高、脉冲形状陡、计时观测精度高、且很
5、少发生自转频率的突变(glitch)http:/ 引用格式:李变,屈俐俐,高玉平.J0613-0200 驾驭铯钟的方法 J.航空学报,2023,44(3):526500.LI B,QU L L,GAO Y P.Methods for J0613-0200 steering a cesium clock J.Acta Aeronautica et Astronautica Sinica,2023,44(3):526500(in Chinese).doi:10.7527/S1000-6893.2021.26500收稿日期:2021-10-11;退修日期:2021-11-05;录用日期:2021-1
6、1-17;网络出版时间:2021-11-23 18:31网络出版地址:https:/ SKA专项(2020SKA0120103)*通信作者E-mail:航空学报526500-2等 特 征,有 些 自 转 周 期 变 化 率 可 达 到 10-1910-21,被誉为自然界最稳定的“时钟”。1982 年Backer 等3发 现 第 一 颗 毫 秒 脉 冲 星 PSR B1937+21 之后,人们便开始了毫秒脉冲星的计时观测,关于脉冲星计时资料分析及脉冲星时算法 的 研 究 也 随 之 兴 起4-7,2012 年 澳 大 利 亚Parkes 天文台的 Hobbs 等得到了可以与地球时TT(BIPMX
7、X)相媲美的脉冲星时8-9。受观测设备、观测技术、星际介质传播延迟等因素的影响,目前毫秒脉冲星最好的计时精度仅为 100 ns,远低于原子钟测量比对精度。但原子时的长期频率稳定度受原子钟性能、季节性变化、周年效应、以及原子时算法等因素的影响,长期稳定度不如脉冲星时。由于脉冲星时没有独立定义时刻起点,因此脉冲星只能被视为长期稳定度极高的频率源,用于其他时频设备(例如:原子钟)输出频率长期性能的验证或校准。近年来随着脉冲星观测设备性能提升,观测到的毫秒脉冲星的数量有所增加、观测精度不断提高。尤其是我国 FAST 的投入运行,将积累大量的、优质的毫秒脉冲星观测资料10,同时中国SKA先导项目已顺利开
8、展,这些为改变目前时频体系单一依赖原子钟的格局,进而提高当前原子时体系的可靠性和稳定性提供有力的数据及技术支持。脉冲星时与原子时具有优势互补的特点,将脉冲星时用于原子时的检测、校准与驾驭,利用二者的优势,可以建立一种能够在保持原子时短期稳定度优势的基础上,又具有脉冲星时长期稳定度的优势综合时间尺度。1驾驭铯钟的参考毫秒脉冲星选取国际脉冲星计时阵(IPTA)2019 年发布了第二批脉冲星计时观测数据(IPTA dr2)。这些观测数据分别来自欧洲脉冲星计时阵(EPTA)、北 美 脉 冲 星 计 时 阵(NANOGrav)和 澳 大 利 亚Parkes脉冲星计时阵(PPTA)。与 IPTA 发布的第
9、一批数据(IPTA dr1)相比,IPTA dr2 中包含65颗毫秒脉冲星,比 IPTA dr1多了 16颗11-12。1.1IPTA dr2概况IPTA dr2 中 共 有 65 颗 毫 秒 脉 冲 星,其 中EPTA 包含 42颗,观测时间跨度 718 a,NANOGrav包含 37颗,观测时间跨度 0.69.2 a,PPTA包含 20颗,观测时间跨度 5.629.4 a。与 IPTA dr1不同,IPTA dr2分别采用脉冲星计时软件 TEMPO2 和 TEMPOnest 拟合计时模型得到 A、B 两个版本的脉冲星计时残差。IPTA dr2 具有以下特点:对于同一颗脉冲星多个观测带宽的情
10、况,根据脉冲到达时间(TOA)的不确定度计算权重,并将权重最大的数据序列作为计时模型的参考,进而得到其他时间序列相对于该参考的系统差;处理方法更全面,计时残差精度更高。例如:对低偏心率脉冲双星,采用 ELL1模型拟合拉普拉斯-拉格朗日参数。对低偏心率中-高倾角的脉冲双星,采用 DDH 模型对轨道周期的 3次谐波振幅和连续谐波的振幅比进行拟合等;计时观测的时间跨度相差大。毫秒脉冲星J1939+2134 的观测间隔最长,约 29.4 a,J0931-1902的观测间隔最短,只有 0.6 a;计时残差数据质量差异大,不仅表现在不同的脉冲星之间,而且同一颗脉冲星不同观测带宽之间的数据质量差异也非常明显
11、。1.2参考毫秒脉冲星选取用于铯钟驾驭的脉冲星参考必须满足测量精度高、观测间隔尽可能长、长期稳定度高,并且计时观测数据的缺失尽可能少等条件。为此,首先 在 IPTA dr2 中 选 取 满 足 观 测 间 隔 15 a、TOA 测量精度TOA 1.5s条件,并且计时观测数据缺失尽可能少的 6颗毫秒脉冲星作为候选参考(图 1)。图 1 是 对 初 步 选 取 的 编 号 分 别 为 J0437-4715、J1640+2224、J1713+0747、J0613-0200、J1744-1134和 J1939+2134的 6颗毫秒脉冲星不同观测波段计时观测数据,采用 TEMPO2 软件处理后的计时残差
12、。图中横坐标是约化儒略日(MJD)表示的时间序列,纵坐标计时残差。利用毫秒脉冲星驾驭铯钟的宗旨是在保持原子钟短期稳定度的基础上,发挥脉冲星长期稳定度的优势。因此,稳定度是必须考虑的关键因素。脉冲星数据采样具有不规则性,其稳定度估计采用与 Allan方差类似的z()估计方法13:航空学报526500-3z()=225c21/2(1)式中:为采样间隔;c为间隔上多项式拟合的三次项系数;表示统计平均。对选取的 6颗候选毫秒脉冲星,采用z()分别计算稳定度,结果如图 2所示。计时噪声也是决定毫秒脉冲星参考选取的一个重要考虑因素。计时噪声由低频信号组成,是计时残差中不可预报的长期变化趋势14-17。表
13、1是 6 颗毫秒脉冲星的观测间隔、TOA 测量精度、计时残差精度及是否含有计时噪声等情况。根据表 1、图 1 及图 2 中 6 颗候选毫秒脉冲星的稳定度,选取 J0613-0200 作为铯钟驾驭的参考,开展毫秒脉冲星驾驭铯钟的实验研究。2J0613-0200驾驭铯钟实验2.1实验系统利用中国科学院国家授时中心(NTSC)的一台 5071A 铯钟(编号:Cs2928)、相位微调仪、计图 16颗毫秒脉冲星的计时残差Fig.1Timing residuals of 6 millisecond pulsars图 26颗毫秒脉冲星的稳定度Fig.2Stabilities of 6 millisecond
14、 pulsars表 16颗毫秒脉冲星比较Table 1Comparison of 6 millisecond pulsars毫秒脉冲星J0437-4715J1713+0747J1744-1134J1640+2224J0613-0200J1939+2134观测间隔/a18.622.519.917.216.029.4TOA测量精度/s0.110.210.730.771.140.24计时残差精度/s1.121.240.810.851.2150.63计时噪声有无有无无有航空学报526500-4数器和相位频率分配放大器等相关时频设备以及 IPTA dr2 的计时观测资料,建立单颗毫秒脉冲星驾驭铯钟实验系
15、统(图 3)。图 3 中:f为 10 MHz 频率信号;p 为 1 PPS 信号。时间基准为 NTSC 建立和保持的中国标准时间 UTC(NTSC),作为计数器的开门信号。计数器测量 Cs2928 及其被驾驭后输出的实验系统时间(APTCs)与 UTC(NTSC)的相位差,用于实验数据的分析处理。TEMPO2 软件用于分析处理 IPTA dr2 中的 J0613-0200 计时观测资料,并基 于 J0613-0200 的 参 数 文 件,定 期 模 拟 用 于Cs2928驾驭的计时观测资料。2.2J0613-0200模拟数据IPTA dr2 发 布 的 J0613-0200 计 时 观 测 数
16、据截止到 2014 年 5 月 18 日(MJD=56 795),为了开展毫秒脉冲星驾驭铯钟的实验研究,根据J0613-0200 的自转参数、天体测量参数及计时残差的统计特性,利用 TEMPO2 软件中的 fake插件,每月定期模拟用于 Cs2928 驾驭的 J0613-0200 计时观测资料。图 4 是利用 J0613-0200 已有观测资料和模拟数据计算得到的计时残差。图中横坐标是 MJD 表示的时间序列,纵坐标是J0613-0200 的计时残差。图中黑色曲线(Obs)是根据 IPTA dr2 中 J0613-0200(MJD=50 93156 795)的计时观测数据,得到的计时残差,灰色
17、曲线(Sim)是根据 J0613-0200 模拟数据得到的计时残差。J0613-0200 模拟数据中计时红噪声的功率谱密度P(f)为18P(f)=P01+()ffc2-2(2)式中:P0为振幅;fc为截止频率;为谱指数。2.3实验系统铯钟性能分析守时原子钟组通常由铯钟和氢钟组成,铯钟拥有优良的长期稳定度(铯钟的闪烁底噪 1101411015,典型值 51015),氢钟拥有出色 的 短 期 稳 定 度(典 型 氢 钟 的 闪 烁 底 噪 1图 3驾驭实验系统Fig.3Experimental system for steering图 4J0613-0200计时残差Fig.4Timing resi
18、duals of J0613-0200航空学报526500-5101511016,典型值 51016),但因存在频率漂移,其长期稳定度取决于频率漂移的变化与大小。本实验的目的是证明毫秒脉冲星驾驭原子钟的可行性,因此,选择没有频率漂移或频率漂移很小的铯钟作为被驾驭的原子钟。Cs2928 是 NTSC 一台连续工作近 8 年的铯钟,该钟的频率稳定性略高于同类铯钟平均水平,并且拥有完整的历史测量资料,噪声特性长期持续稳定。图 5 是 Cs2928 的稳定度曲线。由图 5可知:Cs2928的噪声主要包括白色调频噪声(WFM)、频率随机游走调频噪声(RWFM)和闪烁调频噪声(FFM)。Cs2928 的频
19、率漂移量值很小(11017),被淹没在随机游走频率噪声中。因此,Cs2928 的驾驭实验不考虑频率漂移的影响,只需根据历史测量资料,结合其噪声特点,准确预报其频率变化。2.4驾驭算法实验系统旨在产生和保持既具有 Cs2928 短期稳定度,同时又具有 J0613-0200长期稳定度特性的时间/频率信号。因此,铯钟驾驭实验系统采用双重驾驭的算法。首先,基于 Cs2928 的历史 数 据 特 征 及 噪 声 特 点,以 地 方 原 子 时 TA(NTSC)为参考,计算用于每天驾驭的日驾驭量(offset1)。由 2.3 节内容可知,Cs2928 的频率漂移量值很小,被随机游走噪声淹没。因此,每天的日
20、驾驭量仅考虑 Cs2928 的频率估计;其次,每月根据 fake 插件模拟的 J0613-0200 计时观测资料,分析计算用于每月驾驭的月驾驭量(offset2)。对于月驾驭,需要评估 J0613-0200和实验系统产生的时间 APTCs之间的偏差。以 TAI 为共同参考,评估 J0613-0200与 APTCs之间的偏差:Z=J0613-0200-APTCs=J0613-0200-TAI-()APTCs-TAI(3)采用卡尔曼滤波方法,估计滤波后的相位和频率19-20,即J0613-0200-APTCs=H B0B1(4)式中:B0和B1分别为滤波后的相位和频率;H为测量矩阵。H=10 设X
21、为状态向量;为状态转移矩阵;Zk为测量输入数据(J0613-0200)-APTCs;P为状态向量的协方差矩阵;Q为滤波器的过程噪声矩阵;R为测量噪声矩阵,则 Xk|k-1=k-1Xk-1|k-1Pk|k-1=k-1Pk-1|k-1Tk-1+Qk-1Kk=Pk|k-1HTk()HkPk|k-1HTk+Rk-1Xk|k=Xk|k-1+Kk()Zk-HkXk|k-1Pk|k=()I-KkHkPk|k-1(5)式中:Qk=(Xk|k-Xk|k-1)(Xk|k-Xk|k-1)T滤 波 器 初 始 值 由 前 一 计 算 间 隔 的(J0613-0200)-APTCs计算得到,卡尔曼滤波估计的频率即为用于
22、每月驾驭的月驾驭量 offset2:offset2=B1(6)这样,每天的驾驭量是这 2项之和:总驾驭量=offset1+offset2(7)3实验结果J0613-0200 驾驭 Cs2928 的日驾驭量每天定时送入相位微调仪对 Cs2928 进行初次驾驭;月驾驭量每月根据 J0613-0200的模拟数据计算,对Cs2928 进 行 二 次 驾 驭,输 出 实 验 系 统 时 间APTCs。利用时间间隔计数器将 APTCs时间信号与 NTSC 保持的我国标准时间 UTC(NTSC)进行测量比对,用于实验系统时间 APTCs的评估。作为国际标准时间的 UTC 由 TAI闰秒后得到,通常BIPM滞
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- J0613 0200 驾驭 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。