Research Article  Open Access
Hao Liang, Yafeng Zhan, "A Fast Detection Algorithm for the XRay Pulsar Signal", Mathematical Problems in Engineering, vol. 2017, Article ID 9607821, 5 pages, 2017. https://doi.org/10.1155/2017/9607821
A Fast Detection Algorithm for the XRay Pulsar Signal
Abstract
The detection of the Xray pulsar signal is important for the autonomous navigation system using Xray pulsars. In the condition of short observation time and limited number of photons for detection, the noise does not obey the Gaussian distribution. This fact has been little considered extant. In this paper, the model of the Xray pulsar signal is rebuilt as the nonhomogeneous Poisson distribution and, in the condition of a fixed false alarm rate, a fast detection algorithm based on maximizing the detection probability is proposed. Simulation results show the effectiveness of the proposed detection algorithm.
1. Introduction
Xray pulsars are rapidly rotating neutron stars which could emit Xray signals periodically, stably, and uniquely [1]. Since Downs first proposed interplanetary navigation using pulsars in 1974 [2], Xray Pulsar Navigation (XPNAV) system has been considered in lots of scenarios for autonomous navigation [3–7]. XPNAV gets the spacecraft’s position vector from measuring Time Difference of Arrival (TDOA) of the Xray pulsar signals between the spacecraft and the Solar System Barycenter (SSB). Because of the large distance (thousands of light years) between the Xray pulsars and the solar system, the received Xray pulsar signal is always weak and easy to be disturbed by kinds of noise [8]. It is difficult to recover the Xray pulsar profile [9] and get high precision TDOA, and each observation time for the Xray pulsars may be as long as hundreds of seconds [3, 10]. However, before using the Xray pulsar signals for navigation, we should firstly determine whether we have received the desired Xray pulsar signals, which would help the detector’s altitude controller to confirm the right direction of the Xray pulsar. Unlike measuring TDOA, the signal detection process should be both accurate and fast. Therefore, it is a challenge to detect the weak Xray pulsar signal with short time and limited photons.
The first widely used Xray pulsar signal detection algorithm is based on the Fast Fourier Transform (FFT), which has been proved to be invalid when the signal’s profile is nonsinusoidal or the noise is nongaussian [11]. The GregoryLoredo method parameterizes the number of bins for epoch folding and calculates the average ratio of the probability for the periodic model to the constant model [12]. The GregoryLoredo method is effective on discovering new pulsars when there is no prior knowledge about the signal’s shape; however, it does not use the prior information and is less efficient. The other detection algorithms, like transform [13, 14] and Match Filter (MF) [15], consider the noise obeying the Gaussian distribution after epoch folding. However, that happens when the number of photons in the bins of epoch folding is larger than 20 [8]. For the detection problem, the observation time is not long enough to support approximating the Poisson noise by the Gaussian [16], and the algorithms above have this deficiency. Zhang et al. built the Poisson distribution model of Xray pulsar and proposed a detection method by using the Bayesian theory [16]. However, it simply built the signal and noise by homogeneous Poisson distribution, which is imprecision. Besides that, it did not consider the influences of the unknown signal phase and the Doppler shift, which make it hard to be used in the practical scenario.
For XPNAV, the accurate prior knowledge about the Xray pulsar signal is essential. The prior knowledge includes the Xray background noise’s rate, the signal’s flux, period, and profile. By using the prior information, we may build the Probability Density Function (PDF) [10] and get a fast detection algorithm based on NeymanPearson (NP) criterion [17]. Even in the condition of few photons for detection, the proposed algorithm’s performance could meet detection requirements of XPNAV.
The rest of this paper is organized as follows. Section 2 introduces the Xray pulsar signal model. Section 3 describes the fast detection algorithm in detail. Computer simulations are shown in Section 4 and some conclusions are given in the final part.
2. XRay Pulsar Signal Model
Unlike the radio signal, the Xray pulsar signal shows the particle property, which means the output of Xray detector on the spacecraft is not an analog signal but a discrete arriving time sequence of individual photons. Assuming the observation time is , the number of received photons is , and the PDF of can be described by nonhomogeneous Poisson process [10]:where is the arriving rate of the photons asHere is the Xray background noise’s rate which is timeinvariant in the cosmos environment, is the Xray pulsar signal’s flux, is the Xray pulsar signal’s normalization standard profile template, is the Xray pulsar signal’s phase, is the received frequency, and is the initial phase. Since the Xray pulsar signal is periodic, , where is the period.
3. The Fast Detection Algorithm
3.1. The Optimal NP Algorithm
We divide the observation timeline into equal bins, which is shown in Figure 1. The duration time of each of the bins is . When is small enough, the arriving rate in the th bin could be treated as a determined variable:where ; the supper index represents the transpose of a vector. The number of photons in th bin is a random variable with the PDF:
Since the random variables in the different bins are independent, the Joint Probability Density Function (JPDF) could be expressed aswhere .
If the Xray pulsar signal exists in the received photons, , then the JPDF of could be expressed as
If the Xray pulsar signal does not exist in the received photons, , then the JPDF of could be expressed as
In the condition of a fixed probability of false alarm , NP criterion maximizes the detection probability [17]. Based on NP criterion, we could do the Likelihood Ratio Test (LRT):
If is larger than the threshold , we can determine that we have received the Xray pulsar signal. Otherwise, the Xray pulsar signal does not exist in the received photons. Equation (8) could be simplified as follows:
Since is independent of the measurement variable , the left side of (11) could be added to the threshold, and then
If , there is only one or none photon in each bin, and then (12) could be converted aswhere is the total number of the photons, , is the phase of th photons, which could be calculated by
From (13), the optimal detection variable is . Thus, the steps of the optimal NP detection algorithm for the Xray pulsar signal are as below:(I)The Xray detector onboard receives the photons and notes each photons’ arriving time .(II)Convert the transferred arriving time into using (14).(III)Search the arriving rate of each arriving time from the XPNAV database (since is periodic, the XPNAV database is a table which stores one period of the function ), and calculate the logarithm .(IV)Add up all , and compare the summation with .
The prior knowledge before the detection process includes the XPNAV database about the arriving rate functions and . The XPNAV database could be built by observing the Xray pulsars for months and some explorers have tried to do that, such as the Rossi XRay Timing Explorer (RXTE) from the National Aeronautics and Space Administration (NASA) and the Röntgen Satellite (ROSAT) from the German Aerospace Center (DLR). It is difficult to get an analytical expression for . Monte Carlo simulations would be used to calculate . The calculation process is as below:(I)Set a fixed probability of false alarm .(II)Generate the arriving time of the photons based on the Poisson distribution with the arriving rate .(III)Calculate based on the above steps of the optimal NP detection algorithm.(IV)Repeat step (II) and step (III) for times; is the times of the Monte Carlo simulations.(V)Rank each simulation result in descending order, and set a threshold to make the proportion of greater than equal . is the threshold we pursue.
3.2. The Suboptimal NP Algorithm
Besides the prior knowledge, how to convert the arriving time into the phase is another important issue. is different due to the spacecraft’s position, and could be affected by the Doppler shift due to the spacecraft’s velocity:where is the fixed frequency of the Xray pulsar, is the Doppler frequency, is the spacecraft’s velocity vector, is the Xray pulsar’s direction vector, and is the light velocity. Supposing there are no prior distributions of and , we could use the twodimension search based on the Maximum Likelihood (ML) estimation [10] to deal with the issue.
Assuming ’s searching range is , dividing the range into segments, the th searching phase is . Dividing the ’s searching range into segments, the th searching frequency is . Then, take and into step (II) of the optimal NP detection algorithm and calculate . and in have the maximum probability equaling the real ones, where is the maximum function. Therefore, we could use as the detection variable to be compared with . We may call the algorithm resolving the unknown and as the suboptimal NP detection algorithm.
For the practical navigation application, and are always unknown, so the suboptimal NP detection algorithm is more useful than the optimal one.
4. Computer Simulations
Three weak Xray pulsars (B054069, B174424A, and B182313) are selected to simulate the proposed detection algorithm’s performance. The parameters of the Xray pulsars are shown in Table 1, and the normalization standard profile templates [18] are shown in Figure 2. The Xray background noise’s rate is 0.005 photons/s/cm^{2}. The Xray pulsar signal’s flux could be obtained from Table 1 [1]. The Xray detector’s area is 1 m^{2}, and the fixed probability of false alarm is set as 0.001.

The initial phase obeys the uniform distribution with range . The unknown velocity obeys the uniform distribution with range km/s, in which 7.9 km/s is the first cosmic velocity. The Doppler frequency could be calculated by (16). The method to generate the photons refers to reference [10].
The MF algorithm calculates the correlation peak for the Xray pulsar signal’s profile with the standard signal template [15]. Since the MF algorithm is optimal when the noise is obeying the Gaussian distribution, we choose the MF algorithm for comparison.
We divide the simulation algorithm into three groups: (I) the optimal NP detection algorithm; (II) the suboptimal NP detection algorithm, and the number of segments for searching both and is set as 10; (III) the MF algorithm with unknown information about and , and the number of bins for epoch folding is 64. We perform the Monte Carlo simulations for 1000 times to calculate the detection probability, and the results are shown from Figures 3–5.
From Figures 3–5, the detection probability increases with the increasing observation time. Meanwhile, in the same simulation conditions, the suboptimal NP detection algorithm’s detection probability is better than the MF’s. The MF algorithm is effective when the noise is obeying the Gaussian distribution. However, in the condition of short observation time and limited photons, the noise obeys Poisson distribution, which makes the MF algorithm less effective than the proposed detection algorithm.
Since the limited number of photons, the estimation for and based on ML may be imprecise, and the detection probability of the suboptimal NP detection algorithm is lower than the optimal NP detection algorithm. Therefore, the optimal NP detection algorithm is a performance bound.
5. Conclusion
The detection for the Xray pulsar signal is important for XPNAV and some other navigation scenarios using Xray pulsars. In this paper, we propose a fast detection algorithm for the Xray pulsar signal based on NP criteria. Computer simulations show that, in the condition of short observation time and limited number of photons, the detection probability of the proposed algorithm is still high. Besides that, the fast detection algorithm may be extended to other detection issues when the signals obey the Poisson distribution.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by National Natural Science Foundation of China (61271265 and 61671263) and Tsinghua University Initiative Scientific Research Program (2013089244 and 20161080057).
References
 S. I. Sheikh, The use of variable celestial Xray sources for spacecraft navigation [Ph.D. thesis], University of Maryland, 2005.
 G. S. Downs, Interplanetary Navigation Using Pulsating Radio Sources, NASA, 1974.
 S. I. Sheikh, D. J. Pines, P. S. Ray, K. S. Wood, M. N. Lovellette, and M. T. Wolff, “Spacecraft navigation using Xray pulsars,” Journal of Guidance, Control, and Dynamics, vol. 29, no. 1, pp. 49–63, 2006. View at: Publisher Site  Google Scholar
 A. A. Emadzadeh and J. L. Speyer, “Relative navigation between two spacecraft using Xray pulsars,” IEEE Transactions on Control Systems Technology, vol. 19, no. 5, pp. 1021–1035, 2011. View at: Publisher Site  Google Scholar
 E. Wei, S. Jin, Q. Zhang, J. Liu, X. Li, and W. Yan, “Autonomous navigation of Mars probe using Xray pulsars: modeling and results,” Advances in Space Research, vol. 51, no. 5, pp. 849–857, 2013. View at: Publisher Site  Google Scholar
 X. Kai, W. Chunling, and L. Liangdong, “The use of Xray pulsars for aiding navigation of satellites in constellations,” Acta Astronautica, vol. 64, no. 4, pp. 427–436, 2009. View at: Publisher Site  Google Scholar
 Y. Zhou, P. Wu, and X. Li, “Autonomous navigation for lunar satellite using Xray pulsars with measurement faults,” IET Science, Measurement & Technology, vol. 10, no. 3, pp. 239–246, 2016. View at: Publisher Site  Google Scholar
 J. Hanson, S. Sheikh, P. Graven, and J. Collins, “Noise analysis for Xray navigation systems,” in Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS '08), pp. 704–713, Monterey, Calif, USA, May 2008. View at: Publisher Site  Google Scholar
 D. Zhang, W. Zheng, Y. Wang, and L. Zhang, “Xray pulsar profile recovery based on trackingdifferentiator,” Mathematical Problems in Engineering, vol. 2016, Article ID 4238165, 10 pages, 2016. View at: Publisher Site  Google Scholar
 A. A. Emadzadeh and J. L. Speyer, “On modeling and pulse phase estimation of Xray pulsars,” IEEE Transactions on Signal Processing, vol. 58, no. 9, pp. 4484–4495, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 P. C. Gregory and T. J. Loredo, “Bayesian periodic signal detection: analysis of ROSAT observations of PSR 0540693,” The Astrophysical Journal, vol. 473, no. 2, pp. 1059–1066, 1996. View at: Publisher Site  Google Scholar
 P. C. Gregory and T. J. Loredo, “A new method for the detection of a periodic signal of unknown shape and period,” The Astrophysical Journal, vol. 398, no. 1, pp. 146–168, 1992. View at: Publisher Site  Google Scholar
 L. Wang, L.P. Xu, H. Zhang, and N. Luo, “Pulsar signal detection based on Stransform,” Acta Physica Sinica, vol. 62, no. 13, pp. 596–605, 2013. View at: Publisher Site  Google Scholar
 L. Wang, L.P. Xu, and H. Zhang, “A constant false alarm rate detection algorithm for Xray pulsar signal based on Stransform,” Journal of Astronautics, vol. 35, no. 8, pp. 931–937, 2014. View at: Publisher Site  Google Scholar
 C. Kabakchiev, V. Behar, P. Buist et al., “Detection and estimation of pulsar signals for navigation,” in Proceedings of the IEEE 16th International Radar Symposium, pp. 688–693, Dresden, Germany, June 2015. View at: Google Scholar
 H. Zhang, L.P. Xu, Q. Xie, and N. Luo, “Xray pulsar weak signal detection based on Bayesian estimation,” Acta Physica Sinica, vol. 60, no. 4, pp. 829–835, 2011. View at: Google Scholar
 S. M. Key, Fundamentals of Statistical Signal Processing, Vol. II, Detection Theory, Prentice Hall, 1998.
 HEASARC (High Energy Astrophysics Science Archive Research Center) Browse: Main Interface, http://heasarc.gsfc.nasa.gov/cgibin/W3Browse/w3browse.pl.
Copyright
Copyright © 2017 Hao Liang and Yafeng Zhan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.