一种脉冲星周期的时域估计新方法

脉冲星信号的周期估计在脉冲星搜索、观测以及定时研究中具有重要作用.在分析脉冲星周期频域识别方法的基础上,基于脉冲星观测信号的二阶循环平稳模型,提出了一种周期估计的时域新方法-----最大相关方差搜索法.通过对观测信号按假定周期(可变)分组,计算分组后不同时刻的相关系数,进一步得到系数的方差,然后寻找不同周期下的方差最大值,即可得到脉冲星的辐射周期估计.分别对单峰脉冲星(PSRJ0437-4715)的实测信号、双峰脉冲星(PSRB1821-24)的仿真信号以及ATNFEPN数据库中的多颗脉冲星数据进行了周期估计,均得到了令人满意的结果.对众多脉冲星的计算结果表明,与傅里叶频谱法相比,该时域方法对噪声不敏感,尤其在低信噪比情况下仍表现出较优越的性能;与双谱相干统计量方法^[6]相比,运算复杂度更低;尽管其在处理非连续观测数据时能力不足,但仍然可为微弱脉冲星的实时周期估计提供有力参考.     

基于TOA信息的脉冲星信号周期估计方法

针对脉冲星信号周期估计需逐步搜索、计算量大等问题,提出脉冲星信号周期直接估计算法.根据周期不确定误差对时域上脉冲到达时间(Pulse Time of Arrival,TOA)估计方法的影响,推导了基于TOA信息的脉冲星信号周期估计方法的数学模型.该方法将一组观测到的脉冲光子到达时间序列(Photon Time of Arrivals,TOAs)进行等时间间隔分段,在时域上对每段TOAs进行相位估计获得相应的TOA信息,并根据建立的TOA信息与周期误差的关系采用最小二乘原理估计脉冲星信号周期.对周期准确性的评价不同于传统的依赖折叠轮廓的好坏,而是以TOA与时间图像的斜率作为依据,这一依据更直观、易于理解.理论分析与物理、计算机仿真数据结果表明,本文所提的周期估计方法能够根据短时间脉冲光子到达时间序列获得高精度和高分辨率的脉冲信号周期,有利于X射线脉冲星导航的工程应用.     

基于Bayesian估计的X射线脉冲星微弱信号检测

累积轮廓、流量和周期是X射线脉冲星辐射信号的三个重要特征,将其应用于X射线脉冲星信号检测中,提出了一种基于Bayesian估计的X射线脉冲星周期辐射信号时域检测方法.该方法以非脉冲区噪声观测为先验知识,利用X射线脉冲星辐射信号的泊松分布模型推导了信号概率密度分布函数,以该函数的累积分布函数为判据,对X射线脉冲星微弱信号进行检测,并提取位相偏移量.利用仿真数据和RXTE卫星的实测数据进行实验验证,结果表明:本文方法性能优于同类的基于高斯分布模型的检测方法,在检测信号的同时能在一定精度下给出信号位相偏移值. 关键词： 脉冲星 Bayesian估计 位相测量 时域检测     

空间X射线观测确定脉冲星星历表参数精度分析

脉冲星星历表维持着脉冲星导航所需的时空基准, 其精度直接影响着航天器导航定位结果, 是脉冲星导航系统的基本要素. 本文分析了脉冲星空间观测精度的估计方法, 探索性地研究了基于空间X射线观测获取星历表参数的可行性. 通过建立星历表参数拟合模型, 采用大样本重复事件仿真分析了空间X射线观测精度, 研究了星历表参数确定精度与观测精度、观测时间及观测频次的关系. 研究结果表明, 空间X射线观测可以确定脉冲星星历表参数, 但受限于脉冲星信号特征及探测器技术水平, 当前高精度导航用脉冲星星历表难以通过空间X射线观测手段获得, 可通过地面射电观测技术较好地建立与维持, 提出了推进我国大口径射电望远镜建设的建议.     

基于TOA信息的脉冲星信号周期估计方法（英文）

针对脉冲星信号周期估计需逐步搜索、计算量大等问题,提出脉冲星信号周期直接估计算法.根据周期不确定误差对时域上脉冲到达时间(Pulse Time of Arrival,TOA)估计方法的影响,推导了基于TOA信息的脉冲星信号周期估计方法的数学模型.该方法将一组观测到的脉冲光子到达时间序列(Photon Time of Arrivals,TOAs)进行等时间间隔分段,在时域上对每段TOAs进行相位估计获得相应的TOA信息,并根据建立的TOA信息与周期误差的关系采用最小二乘原理估计脉冲星信号周期.对周期准确性的评价不同于传统的依赖折叠轮廓的好坏,而是以TOA与时间图像的斜率作为依据,这一依据更直观、易于理解.理论分析与物理、计算机仿真数据结果表明,本文所提的周期估计方法能够根据短时间脉冲光子到达时间序列获得高精度和高分辨率的脉冲信号周期,有利于X射线脉冲星导航的工程应用.     

非等间隔计时数据的X射线脉冲星周期快速搜索算法

利用X射线脉冲星信号的计时信息实现航天器深空导航具有重要意义,精确的时间模型(周期及周期导数)是获得高精度自主导航的基础.针对X射线光子并非等间隔地到达探测器的特点,在对比分析两种非等间隔数据处理方法(x2量评估方法和Lomb算法)的基础上,提出了用Lomb算法估计脉冲星周期初值,进一步使用x2量评估方法精化周期值的方法.同时运用快速傅里叶变换算法的思想对Lomb算法进行改进,较大地提高了程序的运算效率.最后利用X射线源仿真系统的计时数据对算法进行试验分析,搜索出准确的脉冲星周期和折叠出正确的脉冲轮廓.     

多谐波脉冲星信号时延估计方法

针对脉冲星导航技术中延时估计这一关键问题, 提出了频域上直接使用脉冲星信号测量到达时间集合进行时延估计的方法——多谐波脉冲星信号时延估计(MHSPE)方法. 该方法建立在频域上相位时延的极大似然估计的基础上, 通过高次谐波对脉冲星观测信号提取出各谐波相位的极大似然估计, 然后取频谱上各谐波的幅值进行归一化作为各谐波相位的权值, 最后取各谐波相位的加权平均作为该时刻的相位估计. 理论上证得MHSPE算法对相位的估计是无偏、一致的, 相比于频域上一次谐波的极大似然估计, MHSPE方法的信噪比随谐波数m的增加而增加, 当各谐波幅值相同时, 信噪比可提高m_1/2倍; 与脉冲星信号时延的克拉美罗界比较, 脉冲星信号时域的导数在频域上的反映就是各谐波分量的数量, 因此随着谐波次数的增加脉冲星信号时延估计可极大趋近克拉美罗界. 采用RXTE航天器对Crab脉冲星的实测数据检验MHSPE方法的性能, 实验结果表明, 针对低信噪比的脉冲星信号, MHSPE可获得高精度的相位估计, 随观测时间增加, 估计精度快速收敛于克拉美罗界.     

脉冲星脉冲轮廓累积的最小熵方法

获取脉冲星辐射信号累积轮廓的传统方法是先识别信号周期,再按周期对信号进行折叠累加,这种方法依赖于周期识别的准确性,对累积轮廓质量也无评价标准.本文分析了脉冲星累积轮廓和周期的强相关性,提出了周期未知时直接进行轮廓累积的最小熵方法,并对其进行证明,在此基础上,讨论了该方法在周期识别中的作用.利用仿真数据和RXTE卫星的实测数据实验验证了该方法的有效性. 关键词： 脉冲星 最小熵 累积轮廓 周期识别     

基于小波Fisz变换的X射线脉冲星信号降噪研究

针对X射线脉冲星信号信噪比低,难以获取精确脉冲星信号相位的问题,提出了基于WaveletFisz变换的X射线脉冲星信号的估计方法.根据X射线脉冲星信号特性建立了信号模型,推导了Wavelet-Fisz变换理论,利用该变换将X射线脉冲星信号逼近于高斯分布,对变换后的信号数据进行估计.对大量RXTE的实测数据进行分析,结果表明在Wavelet-Fisz变换下能获得较好的估计性能,且耗时小,有利于X射线脉冲星导航的工程应用.     

一种新的X射线脉冲星信号模拟方法

由于无法在地面直接观测到X射线脉冲星信号, 且空间飞行探测耗时长, 成本高, 模拟具有真实物理特征的X射线脉冲星信号对X射线脉冲星信号处理方法及导航方案的验证具有重要意义. 本文提出了一种利用太阳系质心处脉冲星信号模型和航天器轨道信息, 建立航天器处实时光子到达速率函数, 再利用尺度变换法产生光子到达时间序列的脉冲星信号模拟新方法. 该方法真实还原了脉冲星信号的频率缓变特性, 考虑了动态探测环境下的相对论效应, 且避免了现有模拟航天器处光子序列方法中的迭代过程, 产生非齐次泊松光子到达时间序列的运算量也低于常用的齐次泊松过程筛选法和反函数递推法. 数值仿真从频率特性、流量特性、轮廓相似度及与实测数据的接近程度四个方面验证了该模拟方法的正确性.     

一种基于最优频段的X射线脉冲星累积轮廓时延估计方法

为提高X射线脉冲星导航中累积脉冲轮廓的时间延迟估计精度, 分析了X射线脉冲星累积脉冲轮廓的频谱特性和现有Taylor快速傅立叶变换时延估计算法的缺陷, 提出了一种基于最优频段的累积轮廓时延估计算法, 并通过建立不同信噪比下时延估计误差与所采用频段之间的关系以确定最优频段. 数值及实测数据实验结果表明: 在短时观测或光子流量较小时, 该算法优于常用的近似最大似然 (FAML)、相关 (CC)、最小二乘 (NLS) 及加权最小二乘 (WNLS) 方法; 在观测时间较长或光子流量较大时, 该算法的估计精度与CC及NLS方法相当, 但其运算量低于NLS, FAML 及WNLS方法. 本文所提算法适用于短时观测脉冲轮廓或低流量脉冲星的高精度时延估计.     

First results of millisecond pulsar 1937+21 timing by TNA-1500 radio telescope of Special Research Bureau of Moscow Power Engineering Institute

Regular timing observations of the millisecond pulsar PSR 1937+21 at a frequency of 610 MHz by means of the TNA-1500 radio telescope with a fully steerable 64-m antenna have begun at the Bear Lakes Radio Observatory of the Special Research Bureau of the Moscow Power Engineering Institute. A 40-channel spectrum analyzer and a system for high-speed recording and acquisition of observation data are connected to a reference clock by means of CAMAC interfaces and an 80286 control computer. The reception band is 1.6 MHz (40 channels of 40 kHz each). The total noise temperature of the system is on the order of 160 K. The signal-to-noise ratio is on the order of 50 for a 30-min observation of the pulsar. The mean-square error of measurement of the pulse arrival time is 400 nsec. Data on pulsar 1937+21 along with timing results for normal reference pulsars PSR 0834+06 and PSR 1919+21 obtained at the Pushino Radio Observatory of the Lebedev Physics Institute over a period of 13 years will be used to construct a pulsar time scale.Scientific-Research Radio-Physics Institute. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 37, No. 11, pp. 1388–1398, November, 1994.     

高稳定度X射线脉冲星信号模拟

X射线脉冲星导航是一种完全自主的导航方式, 在深空乃至行星际空间具有潜在的工程应用价值.由于空间飞行试验系统复杂, 成本巨大, 在实验室环境下高精度地模拟X射线脉冲星信号对数据处理方法和导航方案的验证具有重要意义. 针对当前机械转盘式模拟系统中时间稳定度和轮廓精度的不足, 提出了一种通过产生的轮廓电压信号直接控制可见光光源, 再利用衰减获得光子流, 最后经单光子探测和处理电路输出光子到达时间序列的模拟新方法.该方法实现成本低, 支持任意X射线脉冲星信号的模拟, 且具有高时间稳定度和轮廓精度. 详细地讨论了该方法的原理和涉及的关键技术, 搭建了X射线脉冲星信号模拟系统, 并进行了实验. 实验结果表明: 该系统大幅提高了X射线脉冲星信号的模拟效果, 将模拟脉冲星自转周期的稳定度从现有的10^-4提高到10^-9; 当探测器面积为1 m², 探测能谱范围为2–10 keV, 积分时间为1200 s时, 模拟的PSR B1509-58 观测脉冲轮廓与标准脉冲轮廓的相关系数达到了0.993. 关键词： X 射线光子信号 时间稳定度 观测脉冲轮廓 导航算法     

The effect of radiometric noise on pulsar timing precision

The effect of radio-telescope characteristics on the precision of pulsar timing is studied. We found a definite correlation between the accuracy of long-term pulsar time-series data, instrumental uncertainties of arrival-time timing, and pulsar-flux density. This correlation suggests that modern radio telescopes are still not sensitive enough to investigate the noise that is intrinsic to pulsars or the noise that is due to the propagation of the pulsar signal toward the solar system in the interstellar medium. The next-generation radio telescope, the Square Kilometer Array (SKA), is estimated to have an instrumental uncertainty that will be low enough to be insignificant and not affect the quality of pulsar-timing data.     

一种脉冲星信号模拟新方法

在脉冲星导航地面实验中,为解决地面无法接收脉冲星辐射的X射线信号的问题,必须在实验系统中模拟脉冲星辐射信号. 本文提出两种全新脉冲星信号模拟方法. 第一种方法基于统计物理模型,该方法克服了常用泊松模型模拟方法中仿真适应性窄、仿真速度慢和时间分辨率低等不足,实现了对任意随机过程的脉冲信号进行快速高时间分辨率的仿真;第二种方法基于优化统计模型,该方法针对非齐次泊松过程信号的特性进行优化,仿真速度比常用泊松模型模拟方法提高至少30 倍,对于低流量脉冲星,甚至提高4个数量级,可实现纳秒级时间分辨率仿真. 关键词： 脉冲星导航 非齐次泊松过程 X射线光子信号