基于两级压缩感知的脉冲星时延估计方法

Pulsar time delay estimation method based on two-level compressed sensing

为了快速获得高精度的脉冲星累积脉冲轮廓时延估计,提出了一种基于两级压缩感知的时延估计方法.压缩感知主要包括三个部分:字典、测量矩阵、恢复算法,其中字典尺寸是影响压缩感知估计精度的重要因素.针对压缩感知中字典的原子数增加虽能提高估计精度但又带来计算量大的问题,该方法采用粗估计与精估计两级字典相结合,先利用粗估计字典原子间隔大的特点进行累积脉冲轮廓全相位估计,得到预估时延值,再利用精估计字典的原子间隔小且个数少适合局部估计的特点对累积脉冲轮廓进行精确时延估计.理论分析与实验结果表明:两级字典数据量比传统字典小两个数量级,在相同的时延估计精度下,该方法比传统压缩感知方法计算量大幅度减少,是一种能保持高估计精度并有效降低计算量的脉冲星时延估计方法.

In the traditional compressed sensing Mgorithms, the precision of the time delay estimation is closely related to the number of atoms in the dictionary. The bigger the atom number, the smaller the atomic interval becomes, thus the higher the accuracy of the time delay estimation will be. However, the bigger atom number leads to a higher calculation load. Considering the limited calculation capacity of on-board computer, in order to fast obtain high-accuracy time delay estimation value of the integrated pulsar profile of pulsar in the X-ray pulsar-based navigation, we propose a time delay estimation method based on two-level compression sensing. Compressed sensing mainly includes three parts： the dictionary, the measurement matrix, and the recovery algorithm. Among them, the dictionary size is one of the most important factors that affect the estimation accuracy of the compressed sensing. Aiming to solve the problem of the greater computational load with the increase of the atom number in the dictionary of compressed sensing while improving the accuracy of estimation, we combine the rough estimation with the precision estimation as a two-level dictionary, in the first level, the global phase estimation of the low-dimensional integrated pulsar profile is carried out by making use of the feature of the large atomic interval and the small atomic amount of the rough estimation dictionary. Specifically, first, construct a coarse estimation dictionary according to the low-dimensional standard pulsar profile. Then make dimension reduction sampling on the low-dimensional integrated pulsar profile by the rough estimation measurement matrix based on low-dimensional Hadamard matrix. Finally, use an orthogonal matching pursuit method to obtain the predictive estimation of delay value. In the second level, by taking advantage of the small atomic intervals and numbers of the precise estimation dictionary which are suitable for local estimation, the exact time delay estimation of the high dimensional integrated pulsar profile is performed. Specifically, the original position is first corrected by using the predictive estimation of time delay value, that is, shifting the initial high-dimensional integrated pulsar profile as the input signal of the second level. Then the precise estimation dictionary is constructed according to the partial signal of the length of the high dimension standard pulse profile, using the precise estimation measurement matrix sampling on high-dimensional integrated pulsar profile to obtain measurement value. Finally, the optimal matching position is obtained through the recovery algorithm, which is then combined with the predictive estimation of delay value to calculate the precis time delay estimation value. Theoretical analysis and experimental results show that the quantity of data in the two level dictionary is two orders of magnitude smaller than in the traditional dictionary. The proposed method reduces the computational complexity greatly compared with traditional compression sensing method in the same time delay estimation accuracy. Therefore, this method has the advantages of high precision and small calculation load.