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

Multiple harmonic X-ray pulsar signal phase estimation method

作　　者：宋佳凝 徐国栋[1] 李鹏飞[2] 

机构地区：[1]哈尔滨工业大学,卫星技术研究所,哈尔滨150001 [2]哈尔滨理工大学,电气与电子工程学院,哈尔滨150001

出　　处：《物理学报》2015年第21期490-499,共10页Acta Physica Sinica

基　　金：国家高技术研究发展计划(863计划)(批准号:2008AA8051602)资助的课题~~

摘　　要：针对脉冲星导航技术中延时估计这一关键问题,提出了频域上直接使用脉冲星信号测量到达时间集合进行时延估计的方法——多谐波脉冲星信号时延估计(MHSPE)方法.该方法建立在频域上相位时延的极大似然估计的基础上,通过高次谐波对脉冲星观测信号提取出各谐波相位的极大似然估计,然后取频谱上各谐波的幅值进行归一化作为各谐波相位的权值,最后取各谐波相位的加权平均作为该时刻的相位估计.理论上证得MHSPE算法对相位的估计是无偏、一致的,相比于频域上一次谐波的极大似然估计,MHSPE方法的信噪比随谐波数m的增加而增加,当各谐波幅值相同时,信噪比可提高m1/2倍;与脉冲星信号时延的克拉美罗界比较,脉冲星信号时域的导数在频域上的反映就是各谐波分量的数量,因此随着谐波次数的增加脉冲星信号时延估计可极大趋近克拉美罗界.采用RXTE航天器对Crab脉冲星的实测数据检验MHSPE方法的性能,实验结果表明,针对低信噪比的脉冲星信号,MHSPE可获得高精度的相位估计,随观测时间增加,估计精度快速收敛于克拉美罗界.Pulsars, a small portion of celestial sources that emit radiation at varying intensity, provide new possible navigation algorithms which are different from steady point sources. Time-delay estimation is one of the key aspects of pulsar-based navigation technology. Previous work for pulse phase estimation uses a maximum likelihood estimator（MLE） for the phase-in time domain, which is seen as one of the most useful phase estimators. However, the analytic solution for phase cannot be found using MLE. As a result, a brute-force method involving nested, iterative grid-searches is needed to solve this MLE issue, which leads to lots of computations. In order to solve this problem, a multiple harmonic X-ray pulsar signal phase estimation（MHSPE） method is proposed. This method uses the times of arrivals（TOAs） measured pulsar signal to estimate the time-delay in the frequency domain. In this paper, firstly we use the arrival time to derive the maximum-likehood（ML） estimation of phase-delay by fundamental frequency, then an analytic expression for the fundamental frequency phase is obtained. The MHSPE method based on the fundamental frequency phase equation,calculates different harmonic phases by generalizing the analytic expression of fundamental frequency phase, and the normalized amplitude of each harmonic in the spectrum is used as the weight of each harmonic phase. Finally, the weighted average of harmonic phases, which is given by the final analytic expression, is used as the estimation of the moment. To evaluate the MHSPE method, the error and variance equations are calculated and the MHSPE method is demonstrated to be unbiased and consistent. Moreover, by comparing with the ML estimation of the first harmonic, if the amplitudes of harmonic in the spectrum are almost the same, the signal-to-noise ratio（SNR） of MHSPE improves m1/2times when the number of harmonic waves is m. Compared with the Cramer-Rao bound of pulsar time-delay estimation, the derivative of pulsar signal in the time domain reflec

关 键 词：脉冲星导航 克拉美罗界 极大似然估计 谐波相位 

分 类 号：P128.4[天文地球—天体测量]