与预报中应用.PDFVIP

第卷第期 天 文 学 报 年月 经验模态分解方法在脉冲星计时残差分析 及预报中的应用 1,2,3,4 1,2 1,2 1,2 1,2 高 峰 高玉平 童明雷 杨廷高 赵成仕 中国科学院国家授时中心西安 中国科学院时间频率基准重点实验室西安 中国科学院大学北京 西安科技大学应用物理系西安 摘要 高精度脉冲星计时分析在脉冲星诸多应用领域中占重要地位 而多个未能消除的 误差源的存在 限制了脉冲星计时精度的进一步提高 分析其观测特征 这种影响将导致 脉冲星计时残差波动的多样性 利用从北美计时阵观测的颗毫秒脉冲星数 据 采用经验模态分解方法对计时残差的波动特征进行了研究 通过分析发现 计 时残差的振荡可用几个内在本征模态函数和一个剩余趋势揭示其多时间尺度的振荡结构 和非平稳性 进一步比对脉冲星间多振荡结构对计时残差波动贡献的差异 发现经验模态 分解非常有利于寻找脉冲星计时残差数据间的共性及个性 同时 采用自回归滑动平均模 型进行残差预报 发现应用经验模态分解方法将残差序列分解后再预报 在短期预报时具 有明显的优势 上述结果有利于拓展脉冲星计时残差在时域的分析 关键词 天体测量学 脉冲星 计时残差 方法 数据分析 中图分类号 文献标识码 引言 脉冲星计时分析是整个脉冲星领域研究的基石, 在诸如强引力场中的广义相对论 [14] 效应、检验引力波、探测太阳系外行星以及提供计时标准等研究中 , 都需要高精度 的脉冲星计时分析的结果. 当前, 由于计时模型不完善及一些未能消除的误差源的存在, 限制了脉冲星计时分析精度进一步提高, 其影响主要体现在脉冲星计时残差的波动性 上. 脉冲星计时残差是脉冲星时与原子时比对的产物, 其波动性主要来源于脉冲星计时 观测误差、脉冲到达时间(Time of Arrival, TOA)转换误差、脉冲星本身内部噪声等不 [57] 确定因素的影响. 不同的误差源对计时残差的影响各不相同 , 分析计时残差诸多波 收到原稿 收到修改稿 国家自然科学基金项目、、、、、、 中国科学院西部之光青年学者项目、陕西省青年科技新星项目和 西安科技大学校培育基金项目资助 卷 天 文 学 报 期 动特征, 将有利于解读上述各误差的来源. 参考文献[8 10]已报道了利用指数模型、频 域方法、功率谱方法等实现计时残差中低频噪声的度量及分类. 频域分析的最大优点 是通过傅立叶或小波变换能将波动信号的某些特征清楚地显现出来. 然而, 上述变换所 用的基底函数系和“窗口”分析都具有明显的主观选择性. Huang等人于1998年提出了经 验模态分解(Empirical Mode Decomposition, EMD)方法用于分析非线性非平稳时间序 [11] 列 . EMD方法具有良好的功能: 对数据具有自适应性, 能用几个内在的本征模态函 数(Intrinsic Mode Function, IMF)和一个剩余