# X-Ray Pulsar Profile Recovery Based on Tracking-Differentiator

Authors: Dapeng Zhang; Wei Zheng; Yidi Wang; Lu Zhang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

The profile recovery is an important work in X-ray pulsar-based navigation. It is a key step for the analysis on the pulsar signal’s characteristic and the computing of time of arrival (TOA). This paper makes an argument for an algorithm based on the tracking-differentiator (TD) to recover the profile from the low Signal-to-Noise Ratio (SNR) signals. In the method, a TD filter with cascade structure is designed which has very low phase delay and amplitude distortion. In the simulation experiment, two typical pulsars (PSR B0531+21 and PSR B1937+21) are used to verify the algorithm’s performance. The simulation results show that the method satisfies the application requirements in the aspects of SNR and profile fidelity. By processing the data collected by the Rossi X-Ray Timing Explorer (RXTE) satellite in space, similar results can also be achieved.


---

## Body

## 1. Introduction

A pulsar is a kind of neutron star rotating rapidly in space. An intersection angle exists between the pulsars’ rotation axis and magnetic axis. The magnetic poles of the pulsars radiate X-ray beam. When the X-ray beam sweeps on the measurement device, the device can receive the pulsar’s signal [1]. Since the discovery of the first pulsar, scientists gradually recognized that the pulsar can be used as a natural clock because of its long-term stability. The X-ray pulsar-based navigation was initiated in the 1970s, and since then a complete navigation framework has already been developed. Since 1993, the United States carried out a series of research projects on X-ray pulsar-based navigation [2–5]. The European Space Agency (ESA) studied on the feasibility of X-ray pulsar-based navigation in 2004 [6]. Chinese scientists also began their research in 2004. In the aspects of signal processing and navigation algorithms, many valuable results have been achieved [7–9].According to the distinctions of the energy bands of the pulsars, the pulsars can be divided into X-ray pulsars, radio pulsars, gamma-ray pulsars, and so forth. In order to reduce the mass and power dissipation of sensors, X-ray pulsars are mostly selected as navigation pulsars. Objectively, the signal of the X-ray pulsars is extremely weak so that only discrete photons can be recorded [10, 11]. The background noise of the pulsar signals includes space background noise, high-energy particle noise, nebula noise, and detector noise. These noises can influence the precision of the pulsar profiles and increase navigation errors. Particularly, for the millisecond pulsar (e.g., PSR B1937+21), its signal is almost flooded in the noise. Figure 1 showed the influences caused by noise. Literature indicates that these noises can be modeled as Gauss white noise, and the variance is proportional to the number of the bins of the folding signal [12]. However, to get higher-precision profile needs to select more bins. One alternative method is to execute long-time pulsar observation. Nonetheless, owing to the limitation of the space environments and the task scheme of the aircraft, long-time observation is hard to implement. Therefore, it is quite necessary to research the profile recovery method in a shorter observation time.Figure 1
(a) The profile of PSR B0531+21 when observing 360 s. (b) The profile of PSR B1937+21 when observing 720 s.
(a) (b)The discrete tracking-differentiator (TD) is based on the second-order minimal time system. Since being proposed by Professor Jingqing [13, 14], the method has been applied in active disturbance rejection controllers (ADRC). Based on it, Xie and Long proposed a tracking-differentiator with a higher precision [15]. It can solve the disturbance problem of the signals effectively, track the input signal without overshooting, and obtain high-quality differential signal [16]. The structures of these two kinds of tracking-differentiators are different. Han’s method is based on the transform of TD’s nonlinear border and reachable region. Xie and Long’s method is based on the switching curve and the time of reaching original point. In theory, the criterion of Xie and Long’s method is stricter and faster.


## 2. X-Ray Pulsar Profile Recovery

### 2.1. The Epoch Folding of Discrete Photon Signal of Pulsar

In order to analyze the characteristics of pulsar signal and compute the TOA, the signal profile should be recovered first. At present, the method we usually use is direct epoch folding. The method includes 4 steps. (1) Perform barycenter correction. (2) Uniform the collected TOA photon series to a single pulse period. (3) Divide the period into some equal bins, and put the uniformed photons into their bins. (4) Compute the number of the photons in every bin, and plot the profile of the pulse. The process of direct epoch folding is shown in Figure2.Figure 2
The principle of direct epoch folding.


### 2.2. Profile Recovery Method Based on the Track-Differentiator

#### 2.2.1. The Algorithm of the Track-Differentiator

The discrete form of the track-differentiator can be expressed as follows:(1) x 1 k + 1 x 2 k + 1 = 1 h 0 1 x 1 k x 2 k + h 2 2  h u k - h 3 3  h 2 2 b k , u k ≤ r , k = 0,1 , 2 , … ,where h is the step of sampling, r is the factor of tracking, and b(k) is the correction term. h is selected according to time interval of data collection. r decides the tracking speed and bk can be calculated as (9).If arbitrary pointM(x1,x2) in the phase plane Σ=R2 is not on the switching curve Γ(x1,x2)=x1+x2x2/2r, the time it takes to reach the switching curve is(2)TA=x2sr+x1sr+12x22r2.If M(x1,x2) is on the switching curve, the time it takes to reach the original point is(3)TB=x2r.When M(x1,x2) is not on the switching curve and cannot reach the switching curve in one step, the control law is(4)u=-r·sign⁡Γx1,x2.Otherwise, the states x1 and x2 will be controlled to switching point in one step:(5)u=-r·s·ua,where(6)ua=-12+x2srh+121+4rhx2+2x1hs,s=sign⁡Γx1,x2.If M(x1,x2) is on the switching curve Γ and cannot reach the original point in one step, the control law is(7)u=r·sign⁡x1.Otherwise,(8)u=2hx2+3x1h.The correction term is given as follows:(9)b=6h2x2+2x1h.


## 2.1. The Epoch Folding of Discrete Photon Signal of Pulsar

In order to analyze the characteristics of pulsar signal and compute the TOA, the signal profile should be recovered first. At present, the method we usually use is direct epoch folding. The method includes 4 steps. (1) Perform barycenter correction. (2) Uniform the collected TOA photon series to a single pulse period. (3) Divide the period into some equal bins, and put the uniformed photons into their bins. (4) Compute the number of the photons in every bin, and plot the profile of the pulse. The process of direct epoch folding is shown in Figure2.Figure 2
The principle of direct epoch folding.


## 2.2. Profile Recovery Method Based on the Track-Differentiator

### 2.2.1. The Algorithm of the Track-Differentiator

The discrete form of the track-differentiator can be expressed as follows:(1) x 1 k + 1 x 2 k + 1 = 1 h 0 1 x 1 k x 2 k + h 2 2  h u k - h 3 3  h 2 2 b k , u k ≤ r , k = 0,1 , 2 , … ,where h is the step of sampling, r is the factor of tracking, and b(k) is the correction term. h is selected according to time interval of data collection. r decides the tracking speed and bk can be calculated as (9).If arbitrary pointM(x1,x2) in the phase plane Σ=R2 is not on the switching curve Γ(x1,x2)=x1+x2x2/2r, the time it takes to reach the switching curve is(2)TA=x2sr+x1sr+12x22r2.If M(x1,x2) is on the switching curve, the time it takes to reach the original point is(3)TB=x2r.When M(x1,x2) is not on the switching curve and cannot reach the switching curve in one step, the control law is(4)u=-r·sign⁡Γx1,x2.Otherwise, the states x1 and x2 will be controlled to switching point in one step:(5)u=-r·s·ua,where(6)ua=-12+x2srh+121+4rhx2+2x1hs,s=sign⁡Γx1,x2.If M(x1,x2) is on the switching curve Γ and cannot reach the original point in one step, the control law is(7)u=r·sign⁡x1.Otherwise,(8)u=2hx2+3x1h.The correction term is given as follows:(9)b=6h2x2+2x1h.


## 2.2.1. The Algorithm of the Track-Differentiator

The discrete form of the track-differentiator can be expressed as follows:(1) x 1 k + 1 x 2 k + 1 = 1 h 0 1 x 1 k x 2 k + h 2 2  h u k - h 3 3  h 2 2 b k , u k ≤ r , k = 0,1 , 2 , … ,where h is the step of sampling, r is the factor of tracking, and b(k) is the correction term. h is selected according to time interval of data collection. r decides the tracking speed and bk can be calculated as (9).If arbitrary pointM(x1,x2) in the phase plane Σ=R2 is not on the switching curve Γ(x1,x2)=x1+x2x2/2r, the time it takes to reach the switching curve is(2)TA=x2sr+x1sr+12x22r2.If M(x1,x2) is on the switching curve, the time it takes to reach the original point is(3)TB=x2r.When M(x1,x2) is not on the switching curve and cannot reach the switching curve in one step, the control law is(4)u=-r·sign⁡Γx1,x2.Otherwise, the states x1 and x2 will be controlled to switching point in one step:(5)u=-r·s·ua,where(6)ua=-12+x2srh+121+4rhx2+2x1hs,s=sign⁡Γx1,x2.If M(x1,x2) is on the switching curve Γ and cannot reach the original point in one step, the control law is(7)u=r·sign⁡x1.Otherwise,(8)u=2hx2+3x1h.The correction term is given as follows:(9)b=6h2x2+2x1h.


## 3. The Cascaded Filter Group

The track-differentiator mentioned above can track low-frequency signal quickly. In other words, it has a characteristic of low pass filter (TD filter), as Figure3 shows. However, the signal processed by the single filter may lead to phase delay which cannot be accepted in the profile recovery application. Therefore, a filter group with cascaded structure is designed below. The basic idea is to compensate the delay by multiple outputs of filters which has been applied in several engineering fields [15, 16].Figure 3
The IO model of the TD filter.The relation between the output and input signals is(10)W~t=Wt-τ+ξt,where W~t is output signal, Wt is input signal, and ξt is residual.In order to compensate the delayτ in single filter, a filter group that consists of n TD filters is established as Figure 4 shows.Figure 4
The cascaded filter group.The filter group can be expressed as(11)W1t=Wt-τ+ξ1t,W2t=Wt-2τ+ξ2t,W3t=Wt-3τ+ξ3t,W4t=Wt-4τ+ξ4t,⋮Wnt=Wt-nτ+ξnt.The Taylor expansion of Wt-nτ is(12)Wt-iτ=Wt+dWtdt·-iτ+d2Wtdt2·-iτ22!+d3Wtdt3·-iτ33!+Rt,τ.If we make i=1,2,3,4 (four-order Taylor expansion), (13) W 1 t W 2 t W 3 t W 4 t = 1 - τ τ 2 2 - τ 3 6  1 - 2 τ 2 τ 2 2 - 2 τ 3 6  1 - 3 τ 3 τ 2 2 - 3 τ 3 6  1 - 4 τ 4 τ 2 2 - 4 τ 3 6 W t W ′ t W ′ ′ t W ′ ′ ′ t + ξ 1 t ξ 2 t ξ 3 t ξ 4 t .Furthermore, separate the delay parameter τ,(14) W 1 t W 2 t W 3 t W 4 t = 1 - 1 1 2 - 1 6  1 - 2 2 - 4 3  1 - 3 9 2 - 9 2  1 - 4 8 - 32 3 W t τ W ′ t τ 2 W ′ ′ t τ 3 W ′ ′ ′ t + ξ 1 t ξ 2 t ξ 3 t ξ 4 t .Therefore,(15)WtτW′tτ2W′′tτ3W′′′t=4-64-14.3-9.57-1.83-87-21-33-1W1tW2tW3tW4t-4-64-14.3-9.57-1.83-87-21-33-1ξ1tξ2tξ3tξ4t.The input signal Wt can be denoted as(16)Wt=4W1t-6W2t+4W3t-W4t-Γt,where Γt=4ξ1t-6ξ2t+4ξ3t-ξ4t is residual noise of input signal.Therefore, when the output signal is constructed to be(17)W~t=4W1t-6W2t+4W3t-W4tthe output signal’s phase is the same as the input signal. In other words, the phase delay is compensated. The structure of the filter group is as Figure 5 shows.Figure 5
The structure of the filter group.


## 4. Experiments

### 4.1. Performance Analysis of the Filter Group

In general, the characteristics of amplitude-frequency and phase-frequency can better present the filters performance. Below, the designed filter group can be stimulated by the sinusoidal waves with different frequencies and the Bode diagrams were plotted.Because the data processing here is off-line processing, the step of samplingh can be arbitrarily set as 0.001 and the filtering factor C0 can be set as 10, 12, or 15. According to the experience equation r=2πf·A/h·C0, the tracking factor r can be roughly set as 500 where the cut-off frequency f is 4 and the gain A is 1. The order of Taylor expansion is 4, 5, or 6. The results are shown in Figure 6.Figure 6
(a) Amplitude-frequency curve of the filter group with 4-order Taylor expansion. (b) Phase-frequency curve of the filter group with 4-order Taylor expansion. (c) Amplitude-frequency curve of the filter group with 5-order Taylor expansion. (d) Phase-frequency curve of the filter group with 5-order Taylor expansion. (e) Amplitude-frequency curve of the filter group with 6-order Taylor expansion. (f) Phase-frequency curve of the filter group with 6-order Taylor expansion.
(a) (b) (c) (d) (e) (f)From the results, it can be seen that the larger the filtering factor selected, the higher the cut-off frequency and the overshoot of the gain. Meanwhile, the variation trend of the signal phase is steeper. When the filter group with higher order is used, the overshoot is lower but the ability to suppress the high-frequency component is weaker. As shown in the amplitude-frequency and phase-frequency curve, the advantage of the TD filter group is the fact that, in the low-frequency band, the processed signal has very low amplitude distortion and phase delay. The performance can almost get an ideal value.


### 4.2. Simulation Experiments of Pulsar Profile Recovery

In order to research the application range of the method the paper mentioned above, this section simulated two typical pulsars’ data. The pulsars are PSR B0531+21 and PSR B1937+21. PSR B0531+21 is a young pulsar with very high radiation flux. PSR B1937+21 is an old pulsar with very low radiation flux. The parameters are shown in Table1. In this paper, two aspects are considered to measure the algorithm’s performance. They are SNR (Signal-to-Noise Ratio) and fidelity. The fidelity is measured by two indices.Table 1
Simulation parameters.
Parameters PSR B0531+21 PSR B1937+21 Period (ms) 33.4 [17] 1.56 [4] Flux (ph/s) 660 [18] 0.5 [4] Flux of background noise (ph/s) 13860 [18] 50a Observation time (s) 3600 7200 Bins 334000 15600 aEstimated value. Cosmetic noise and detector noise are taken into consideration.


#### 4.2.1. Relative Error of the Recovered Profileer

Consider(18)er=∑i=1binλ~ϕi-λϕi2∑i=1binλ2ϕi×100%,where λ~· is the recovered profile, λ· is the standard profile, and bin is the number of the bins.


#### 4.2.2. The Phase Bias of the Recovered ProfileΔϕ

This value is computed by the method of cross-correlation between the recovered profile and standard profile. In this paper, the results are obtained by 1000 Monte Carlo simulations.The parameters of the filter group arer=800, h=0.001, C0=50, and the order of the Taylor expansion is 5.The performance of the recovered profile is shown in Table2 and the result pictures are shown in Figure 7.Table 2
Performance of the recovered profile.
PSR B0531+21 PSR B1937+21 SNR (dB) 20.9 17.06 e r (%) 0.82 1.99 Δ ϕ (cycle) 4.9 × 10−5 1.34 × 10−3Figure 7
(a) The standard profile of PSR B0531+21. (b) The standard profile of PSR B1937+21. (c) The original wave of PSR B0531+21. (d) The original wave of PSR B1937+21. (e) The recovered profile of PSR B0531+21. (f) The recovered profile of PSR B1937+21.
(a) (b) (c) (d) (e) (f)


### 4.3. RXTE Data Experiments of Pulsar Profile Recovery

The RXTE (Rossi X-Ray Timing Explorer) satellite was launched by the USA in 1995. It has operated over 16 years and obtains a large amount of pulsar photon data [19, 20]. The detector on RXTE has high time resolution and high sensitivity. The data observed by RXTE has supported a large number of scientific researches [21–24].In order to verify the feasibility of the method, the data detected by RXTE is used in this section. We choose three pulsars to verify the method. They are PSR B0531+21, PSR B1509-58, and PSR B0540-69. The basic parameters are shown in Table3. In this experiment we used PCA Event Encoded Mode data. The PSR B0521+21 data packages we used are “20804-01-01-00,” “40805-01-05-01,” and “50098-01-01-00.” The PSR B1509-58 data package is “94803-01-11-06” and PSR B0540-69 data package is “10218-01-01-00.” Heasoft (v6.17) was used to process the observational data as follows. (1) Remove clock events with the “sefilter” command for PSR B0531+21 or combine FITS files with the “se_make” command for PSR B1509-58 and PSR B0540-69. (2) Create filter file with the “xtefilt” command. (3) Generate the Good Time Interval (GTI) file with the “maketime” command. (4) Create good events file with the “grosstimefilt” command. (5) Correct the barycenter with the “faxbary” command. Therefore the original experiment data is obtained. The data information is shown in Table 4. Taking advantage of the method of χ2 searching, the signal’s period is computed and the algorithm is as the following equation shows:(19)χ2=∑i=1mni-N/m2N/m,P^=arg⁡max⁡χ2,where m is the bin number of the profile and total photon number N=∑ni, where ni is the photon number in the ith bin. And then, by epoch folding, the original profile waveform is obtained.Table 3
Pulsar parameters [4].
Parameters Period (ms) Fluxa (ph/cm2/s) Distance (kpc) Galactic longitude (deg.) Galactic latitude (deg.) PSR B0531+21 33.4 1.54 2 184.56 −5.78 PSR B1509-58 150.23 1.62E − 02 4.3 320.32 −1.16 PSR B0540-69 50.37 5.15E − 03 47.3 279.72 −31.52 aThe energy range is 2–10 keV.Table 4
Experiment data information.
Parameters Observational ID Effective data length Discriminant conditions of effective data elva Offsetb num_pcu_onc PSR B0531+21 20804-01-01-00 9921 s >10 <0.01 ≥1 40805-01-05-01 8789 s >10 <0.01 ≥1 50098-01-01-00 1829 s >10 <0.01 ≥1  PSR B1509-58 94803-01-11-06 10298 s >10 <0.01 ≥1  PSR B0540-69 10218-01-01-00 21816 s >10 <0.01 ≥1 aelv means elevation angle. bOffset means pointing position. cnum_pcu_on means the number of the working sensors.Use the profile recovery method mentioned above. The data processing results of PSR B0531+21 are shown in Figure8.Figure 8
The experiment results of PSR B0531+21.
(a) (b) (c) (d) (e) (f)The data processing results of PSR B1509-58 are shown in Figure9 and of PSR B0540-69 in Figure 10.Figure 9
The experiment results of PSR B1509-58.
(a) (b)Figure 10
The experiment results of PSR B0540-69.
(a) (b)Figures8–10 show the recovered profile based on the tracking-differentiator. By comparing with the simulation results in Section 4.2, the results in this section are consistent with them. The method is proved to be effective.


## 4.1. Performance Analysis of the Filter Group

In general, the characteristics of amplitude-frequency and phase-frequency can better present the filters performance. Below, the designed filter group can be stimulated by the sinusoidal waves with different frequencies and the Bode diagrams were plotted.Because the data processing here is off-line processing, the step of samplingh can be arbitrarily set as 0.001 and the filtering factor C0 can be set as 10, 12, or 15. According to the experience equation r=2πf·A/h·C0, the tracking factor r can be roughly set as 500 where the cut-off frequency f is 4 and the gain A is 1. The order of Taylor expansion is 4, 5, or 6. The results are shown in Figure 6.Figure 6
(a) Amplitude-frequency curve of the filter group with 4-order Taylor expansion. (b) Phase-frequency curve of the filter group with 4-order Taylor expansion. (c) Amplitude-frequency curve of the filter group with 5-order Taylor expansion. (d) Phase-frequency curve of the filter group with 5-order Taylor expansion. (e) Amplitude-frequency curve of the filter group with 6-order Taylor expansion. (f) Phase-frequency curve of the filter group with 6-order Taylor expansion.
(a) (b) (c) (d) (e) (f)From the results, it can be seen that the larger the filtering factor selected, the higher the cut-off frequency and the overshoot of the gain. Meanwhile, the variation trend of the signal phase is steeper. When the filter group with higher order is used, the overshoot is lower but the ability to suppress the high-frequency component is weaker. As shown in the amplitude-frequency and phase-frequency curve, the advantage of the TD filter group is the fact that, in the low-frequency band, the processed signal has very low amplitude distortion and phase delay. The performance can almost get an ideal value.


## 4.2. Simulation Experiments of Pulsar Profile Recovery

In order to research the application range of the method the paper mentioned above, this section simulated two typical pulsars’ data. The pulsars are PSR B0531+21 and PSR B1937+21. PSR B0531+21 is a young pulsar with very high radiation flux. PSR B1937+21 is an old pulsar with very low radiation flux. The parameters are shown in Table1. In this paper, two aspects are considered to measure the algorithm’s performance. They are SNR (Signal-to-Noise Ratio) and fidelity. The fidelity is measured by two indices.Table 1
Simulation parameters.
Parameters PSR B0531+21 PSR B1937+21 Period (ms) 33.4 [17] 1.56 [4] Flux (ph/s) 660 [18] 0.5 [4] Flux of background noise (ph/s) 13860 [18] 50a Observation time (s) 3600 7200 Bins 334000 15600 aEstimated value. Cosmetic noise and detector noise are taken into consideration.


### 4.2.1. Relative Error of the Recovered Profileer

Consider(18)er=∑i=1binλ~ϕi-λϕi2∑i=1binλ2ϕi×100%,where λ~· is the recovered profile, λ· is the standard profile, and bin is the number of the bins.


### 4.2.2. The Phase Bias of the Recovered ProfileΔϕ

This value is computed by the method of cross-correlation between the recovered profile and standard profile. In this paper, the results are obtained by 1000 Monte Carlo simulations.The parameters of the filter group arer=800, h=0.001, C0=50, and the order of the Taylor expansion is 5.The performance of the recovered profile is shown in Table2 and the result pictures are shown in Figure 7.Table 2
Performance of the recovered profile.
PSR B0531+21 PSR B1937+21 SNR (dB) 20.9 17.06 e r (%) 0.82 1.99 Δ ϕ (cycle) 4.9 × 10−5 1.34 × 10−3Figure 7
(a) The standard profile of PSR B0531+21. (b) The standard profile of PSR B1937+21. (c) The original wave of PSR B0531+21. (d) The original wave of PSR B1937+21. (e) The recovered profile of PSR B0531+21. (f) The recovered profile of PSR B1937+21.
(a) (b) (c) (d) (e) (f)


## 4.2.1. Relative Error of the Recovered Profileer

Consider(18)er=∑i=1binλ~ϕi-λϕi2∑i=1binλ2ϕi×100%,where λ~· is the recovered profile, λ· is the standard profile, and bin is the number of the bins.


## 4.2.2. The Phase Bias of the Recovered ProfileΔϕ

This value is computed by the method of cross-correlation between the recovered profile and standard profile. In this paper, the results are obtained by 1000 Monte Carlo simulations.The parameters of the filter group arer=800, h=0.001, C0=50, and the order of the Taylor expansion is 5.The performance of the recovered profile is shown in Table2 and the result pictures are shown in Figure 7.Table 2
Performance of the recovered profile.
PSR B0531+21 PSR B1937+21 SNR (dB) 20.9 17.06 e r (%) 0.82 1.99 Δ ϕ (cycle) 4.9 × 10−5 1.34 × 10−3Figure 7
(a) The standard profile of PSR B0531+21. (b) The standard profile of PSR B1937+21. (c) The original wave of PSR B0531+21. (d) The original wave of PSR B1937+21. (e) The recovered profile of PSR B0531+21. (f) The recovered profile of PSR B1937+21.
(a) (b) (c) (d) (e) (f)


## 4.3. RXTE Data Experiments of Pulsar Profile Recovery

The RXTE (Rossi X-Ray Timing Explorer) satellite was launched by the USA in 1995. It has operated over 16 years and obtains a large amount of pulsar photon data [19, 20]. The detector on RXTE has high time resolution and high sensitivity. The data observed by RXTE has supported a large number of scientific researches [21–24].In order to verify the feasibility of the method, the data detected by RXTE is used in this section. We choose three pulsars to verify the method. They are PSR B0531+21, PSR B1509-58, and PSR B0540-69. The basic parameters are shown in Table3. In this experiment we used PCA Event Encoded Mode data. The PSR B0521+21 data packages we used are “20804-01-01-00,” “40805-01-05-01,” and “50098-01-01-00.” The PSR B1509-58 data package is “94803-01-11-06” and PSR B0540-69 data package is “10218-01-01-00.” Heasoft (v6.17) was used to process the observational data as follows. (1) Remove clock events with the “sefilter” command for PSR B0531+21 or combine FITS files with the “se_make” command for PSR B1509-58 and PSR B0540-69. (2) Create filter file with the “xtefilt” command. (3) Generate the Good Time Interval (GTI) file with the “maketime” command. (4) Create good events file with the “grosstimefilt” command. (5) Correct the barycenter with the “faxbary” command. Therefore the original experiment data is obtained. The data information is shown in Table 4. Taking advantage of the method of χ2 searching, the signal’s period is computed and the algorithm is as the following equation shows:(19)χ2=∑i=1mni-N/m2N/m,P^=arg⁡max⁡χ2,where m is the bin number of the profile and total photon number N=∑ni, where ni is the photon number in the ith bin. And then, by epoch folding, the original profile waveform is obtained.Table 3
Pulsar parameters [4].
Parameters Period (ms) Fluxa (ph/cm2/s) Distance (kpc) Galactic longitude (deg.) Galactic latitude (deg.) PSR B0531+21 33.4 1.54 2 184.56 −5.78 PSR B1509-58 150.23 1.62E − 02 4.3 320.32 −1.16 PSR B0540-69 50.37 5.15E − 03 47.3 279.72 −31.52 aThe energy range is 2–10 keV.Table 4
Experiment data information.
Parameters Observational ID Effective data length Discriminant conditions of effective data elva Offsetb num_pcu_onc PSR B0531+21 20804-01-01-00 9921 s >10 <0.01 ≥1 40805-01-05-01 8789 s >10 <0.01 ≥1 50098-01-01-00 1829 s >10 <0.01 ≥1  PSR B1509-58 94803-01-11-06 10298 s >10 <0.01 ≥1  PSR B0540-69 10218-01-01-00 21816 s >10 <0.01 ≥1 aelv means elevation angle. bOffset means pointing position. cnum_pcu_on means the number of the working sensors.Use the profile recovery method mentioned above. The data processing results of PSR B0531+21 are shown in Figure8.Figure 8
The experiment results of PSR B0531+21.
(a) (b) (c) (d) (e) (f)The data processing results of PSR B1509-58 are shown in Figure9 and of PSR B0540-69 in Figure 10.Figure 9
The experiment results of PSR B1509-58.
(a) (b)Figure 10
The experiment results of PSR B0540-69.
(a) (b)Figures8–10 show the recovered profile based on the tracking-differentiator. By comparing with the simulation results in Section 4.2, the results in this section are consistent with them. The method is proved to be effective.


## 5. Conclusions

In order to enhance the SNR of the recovered profiles of the pulsars, the paper introduces a profile recovery method based on tracking-differentiator. Applying the method, a cascaded structure filter group is constructed. The filter group reduces the phase distortion and enhances the performance of the recovered profile. By simulating the photon data of the two typical pulsars (PSR B0531+21 and PSR B1937+21), the paper analyzes the performance of the method. The results show that the method can meet requirements in the aspects of the SNR and fidelity and is appropriate for both young pulsars and old pulsars. By processing the measured data from RXTE satellite, the simulation results are verified. The method does not need the a priori information of the pulsar standard profile, and it can be applied in the X-ray pulsar-based navigation fields.



---

*Source: 4238165-2016-05-31.xml*