Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Method for Extracting Joints on Mountain Tunnel Faces Based on Mask R-CNN Image Segmentation Algorithm
Previous Article in Journal
State-of-the-Art Research on Loess Microstructure Based on X-ray Computer Tomography
Previous Article in Special Issue
Three-Dimensional Guidance Laws for Spacecraft Propelled by a SWIFT Propulsion System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

X-ray Pulsar-Based Navigation Using Pulse Phase Delay between Spacecraft and Verification with Real Data

1
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2
Beijing Institute of Tracking and Telecommunications Technology, Beijing 100094, China
3
College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China
4
China Astronaut Research and Training Center, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(15), 6401; https://doi.org/10.3390/app14156401
Submission received: 19 May 2024 / Revised: 5 July 2024 / Accepted: 11 July 2024 / Published: 23 July 2024
(This article belongs to the Special Issue Advances in Deep Space Probe Navigation)

Abstract

:
Pulsars are neutron stars with high rotation speeds and have extraordinary long-term rotational stability. X-ray pulsar-based navigation (XNAV) is a navigation method that estimates the position and velocity of a spacecraft using the X-ray radiation from pulsars. Flight experiments on Insight-Hard X-ray Modulation Telescope (Insight-HXMT) and Neutron Star Interior Composition Explorer (NICER) have successfully verified the feasibility of using XNAV for a single spacecraft. For spacecraft in formation, a pulsar-based navigation method that uses the pulse phase delay between spacecraft is derived. Moreover, a direct estimation method for pulse phase delay, which is independent from the pulsar template, is proposed. The proposed method is verified with simulation data of the Crab pulsar and real data of the same pulsar obtained from Insight-HXMT and NICER.

1. Introduction

With the development of aerospace technology, on-orbit spacecraft missions have become more complex [1]. Formation spacecraft have the advantages of low system cost and high flexibility [2,3]. Compared to a single spacecraft, spacecraft formation flying is more flexible for complex space missions [4]. Thus, spacecraft formation flying has attracted a lot of attention [5,6]. Furthermore, navigation services are important for spacecraft formation flying missions [7].
XNAV is a developing autonomous navigation method providing navigation services for Earth-orbit and deep-space spacecraft [8,9]. The basic idea of the XNAV for a single spacecraft is similar with that of the Global Navigation Satellite System. A spacecraft estimates its position and velocity using received photons from different pulsars [10]. Since its introduction in 1981, a significant growth in the navigation algorithm for XNAV has occurred in the past four decades [11,12]. In addition, some flight experiments, such as NICER and Insight-HXMT, were performed [13,14,15,16].
Aside from providing navigation service for a single spacecraft, XNAV can also be used for spacecraft in formation [17,18,19]. In [17], Xiong proposed the use of time difference of arrival (TDOA) measurements from pulsars and pseudo-range observations from inter-satellite links to determine the position of a satellite in the constellation. In [18], Zhang proposed using TDOA, phase increment, the baseline of double satellites, and the plane formed by three satellites for spacecraft formation navigation. In [19], Emadzadeh proposed the use of pulse delay between two spacecraft to estimate the relative position and studied the effect of velocity errors on position estimation. However, in [19], the pulse delay is calculated after pulse phases of spacecraft are estimated, respectively. This can lead to unnecessary redundant calculations. Moreover, in [19], the pulse phase of spacecraft is estimated using the template of pulsar and empirical profiles recovered form pulsar photons. However, this can affect the estimation performance of the pulse phase based on the accuracy of the template. In [17,18,19], the proposed methods were only verified through simulation.
In this paper, we first described a pulse phase estimation method and the measurement model of XNAV for a single spacecraft. Then, taking the formation of two spacecraft as an example, a pulsar-based navigation method is used to derive the pulse phase delay between spacecraft. To improve computational efficiency and eliminate dependence on the template of pulsars, a direct estimation method for pulse delay without using a template is proposed. The proposed method is verified with simulation data of the Crab pulsar and real data of the Crab pulsar derived from NICER and Insight-HXMT. For convenience, we only consider the Romer delay and Shapiro delay when determining the transform photon arrival time from spacecraft to solar system barycenter (SSB). To provide an accurate navigation service for spacecraft near gravitational sources, an accurate transform model should be considered in further studies [20].
The remainder of the paper is organized as follows: The related method of XNAV for a single spacecraft is shown in Section 2; Section 3 derives a pulsar-based navigation method using the pulse phase delay between spacecraft, and a direct estimation method is proposed for pulse phase delay. In Section 4, the simulation data of the Crab pulsar and the real data of the Crab pulsar derived from NICER and Insight-HXMT are used to verify the proposed method. The conclusions are presented in Section 5.

2. Related Method of XNAV for Single Spacecraft

2.1. Measurement Model

The measurement of pulsar-based navigation is the pulse phase. The pulse phase at t i can be expressed as [16] follows:
ϕ ( t i ) = F 0 [ g ( t i ) T 0 ] + F 1 2 [ g ( t i ) T 0 ] 2 ,
where F 0 and F 1 are the frequency of the pulsar and the time derivative of the frequency at T 0 , respectively, and F 0 and F 1 can be found in the public pulsar ephemeris [21]. In addition, g ( t i ) can be described by [20,22] the following:
g ( t i ) = t i + n ( r i + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i + r k , i | ,
where c is the speed of light, r i is position vector of the spacecraft relative to the Earth at t i , n is the direction vector of the pulsar, r E , i is the position vector of Earth relative to the SSB at t i , r k , i is the position of the spacecraft relative to the kth celestial body(unit m), μ k is the gravitational coefficient of the kth celestial body, and ‖ ‖ is the magnitude of the vector. The second and third terms on the right side of Equation (2) are the Romer delay and the Shapiro delay, respectively. The Romer delay represents a vacuum delay between the spacecraft and the SSB, while the Shapiro delay represents a time delay caused by the pulse passing through curved space–time.
Assuming that x i = x ( t i ) = [ r i , v i ] , y i = ϕ ( t i ) = ϕ i , and h ( x i ) = F 0 [ g ( t i ) T 0 ] + F 1 2 [ g ( t i ) T 0 ] 2 , v i is the velocity of the spacecraft relative to the Earth at t i . Thus, Equation (1) can be rewritten as [16] follows:
y i = h ( x i ) .

2.2. Pulse Phase Estimation Method

The pulsar is very weak such that a spacecraft can only receive a series of events rather than a continuous pulsed signal. Thus, a spacecraft must estimate the pulse phase from received events [23].
Let us assume a spacecraft observes a pulsar from t i to t i + 1 , i = 1, 2, …I and receives N events { τ j } j = 1 N . At first, { τ j } j = 1 N should be corrected to SSB by (2); thus, we have 16:
τ j , S S B = g ( τ j ) .
Therefore, the pulse phase at τ j , S S B is expressed as
ϕ ( τ j , S S B ) = ϕ i + F i [ g ( τ j ) g ( t i ) ] + 1 2 F ˙ i [ g τ j g t i ] 2 ,
where ϕ i = ϕ ( t i ) is the pulse phase at g ( t i ) , F i is the frequency at g ( t i ) , and F ˙ i is the time derivative of the frequency at g ( t i ) , respectively.
However, the true position of the spacecraft, r ( τ j ) in (2), is inaccessible in an autonomous navigation task. r ( τ j ) is expressed as r ( τ j ) = r ~ ( τ j ) + δ r ( τ j ) to estimate the pulse phase, where r ~ ( τ j ) is the predicted position of the spacecraft and δ r ( τ j ) is the error of r ~ ( τ j ) . Moreover, δ r ( τ j ) can be derived as functions of δ v ( τ i ) and δ r ( τ i ) [24,25,26]:
δ r ( τ j ) = Φ r r ( τ j , τ i ) δ r ( τ i ) + Φ r v ( τ j , τ i ) δ v ( τ i ) ,
where δ v ( τ i ) is the error of the predicted velocity of the spacecraft v ~ τ j .   I n a d d i t i o n , Φ r r ( τ j , τ i ) and Φ r v τ j , τ i can be expressed as [24,25] follows:
Φ r r ( τ j , τ i ) = I 3 × 3 + m = 1 1 m ! φ m τ j τ i m Φ r v ( τ j , τ i ) = m = 1 1 m ! γ m τ j τ i m ,
where φ m and γ m are constant matrixes that can be calculated using analytical, numerical, or finite difference techniques [26].
Therefore, (5) can be rewritten as [27]:
ϕ τ j , S S B = ϕ i + F i D o p p l e r τ j t i + ϕ ~ τ j , S S B ,
where
ϕ ~ ( τ j , S S B ) = F i [ g ~ ( τ j ) g ~ ( t i ) ] + 1 2 F ˙ i [ g ~ τ j g ~ t i ] 2 ,
g ~ ( t i ) = t i + n r ~ t i + r E t i c + k = 1 K 2 μ k c 3 ln | n r ~ k ( t i ) + r ~ k t i | ,
F i o r b i t = F i [ G i φ 1 δ r ( t i ) + G i γ 1 δ v ( t i ) ] ,
G i = g ( t i ) r ( t i ) | r ( t i ) =   r ~ ( t i ) = 1 c n + 2 k μ k c 3 n   r ~ k ( t i ) +   r ~ k ( t i ) | n ·   r ~ k ( t i ) +   r ~ k ( t i ) |   r ~ k ( t i ) .
In the equations above, ϕ ~ τ j , S S B and g ~ ( t i ) are the pulse phase and photon arrival time calculated by r ~ ( t i ) , respectively. Applying (8) to every event { τ j } j = 1 N , we have { ϕ ( τ j , S S B ) } j = 1 N . Thus, we have to estimate ϕ 0 i and F i o r b i t from { ϕ ( τ j , S S B ) } j = 1 N .
Here, F i o r b i t and ϕ 0 i can be estimated using the epoch folding method. The observation period is assumed to consist of M pulsar periods in the epoch folding method. The period can be divided into M b equal-length bins. The length of each bin is assumed as t b , and c l m is the number of detected photons in the l t h bin of the m t h period. Therefore, the empirical profile can be recovered. Here, we define the empirical profile function at φ l as [19] the following:
λ E F φ l ; F i o r b i t = 1 M t b m = 1 M c l m φ l ; F i o r b i t ,
where φ l is the pulse phase at the middle of l t h bin.
We adopt the Chi-squared test to calculate the statistical characteristics of the empirical profile as follows [28]:
χ 2 = l = 1 M b ( C l C ¯ ) 2 C ¯ ,
where
C l = m = 1 M c l m ( φ l ; F i o r b i t ) ,
C ¯ = m = 1 M l = 1 M b c l m ( φ l ; F i o r b i t ) M b ,
Using Equations (14)–(16), υ ¯ 1 can be estimated with
F ^ i o r b i t = arg m a x F i D o p p l e r χ 2 .
I n a d d i t i o n ,   ϕ i can be estimated by solving a nonlinear least squares estimation problem [19]:
J ϕ i = l = 0 M b λ E F ( φ l ; F ^ i o r b i t ) λ T e ( φ l ) 2 ,
where λ T e is the template of the pulsar.
At the same time, ϕ i can be estimated with
ϕ ^ i = arg min J ( ϕ i ) ϕ i [ 0 , 1 ) .

3. XNAV for Spacecraft in Formation

3.1. Measurement Model

The diagram of XNAV for spacecraft in formation is shown in Figure 1. Assume that spacecraft A and spacecraft B observe the same pulsar; thus, according to (1), the pulse phases are given by
ϕ i A = ϕ A ( t i ) = F 0 [ g A ( t i ) T 0 ] + F 1 2 [ g A ( t i ) T 0 ] 2 ,
ϕ i B = ϕ B ( t i ) = F 0 [ g B ( t i ) T 0 ] + F 1 2 [ g B ( t i ) T 0 ] 2 ,
where
g A ( t i ) = t i + n ( r i A + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i A + r k , i A | ,
g B ( t i ) = t i + n ( r i B + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i B + r k , i B | ,
where r i A is the position of spacecraft A relative to the Earth at t i   and r i B is the position of spacecraft B relative to the Earth at t i . Given that the distance between the two spacecraft is far less than the distances between the celestial bodies and the spacecraft, we assume that r k , i A r k , i B , and we denote it as r k , i .
We use pulse phase delay between spacecraft A and spacecraft B as a measurement:
Δ ϕ i A , B = ϕ i A ϕ i B
Substituting (20) and (21) into (24), we have
Δ ϕ i A , B = n ( r i A r i B ) c · [ F 0 + F 1 · ( t i + n ( r i A + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i + r k , i | T 0 ) F 1 2 · n ( r i A r i B ) c ) ] .
Since c is about 3 × 108 m/s, F 1 is smaller than 10−10 HZ/s for most pulsars. Assume that the distance between spacecraft A and B is shorter than 100 km. Therefore, F 1 2 n ( r i A r i B ) c is close to 0, and we have
Δ ϕ i A , B = n ( r i A r i B ) c · [ F 0 + F 1 · ( t i + n ( r i A + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i + r k , i | T 0 ) ] .
We also define F i A , B as follows:
F i A , B = F 0 + F 1 · ( t i + n ( r i A + r E , i ) c + k = 1 K 2 μ k c 3 ln | n r k , i + r k , i | T 0 ) .
Substituting (27) into (25), we have
Δ ϕ i A , B = F i A , B n r i A r i B . c
Based on Equation (28), the difference in pulse phases of spacecraft A and spacecraft B reflects their distance in the direction of the pulsar. Assuming that Δ x i A , B = x A ( t i ) x B ( t i ) , y i , A B = Δ ϕ i A , B and h A , B ( Δ x i A , B ) = F i A , B n ( r i A r i B ) c , then Equation (28) can be rewritten as a standard measurement model of
y i , A B = h A , B ( Δ x i A , B ) .

3.2. Pulse Phase Delay Estimation

In this section, we explain how the pulse phase delay can be abstained. Assume that spacecraft A and spacecraft B observe the same pulsar from t i to t i + 1 , i = 1, 2,…I and receives { τ j A } j = 1 N A and { τ j B } j = 1 N B , respectively. Applying (8) to { τ j A } j = 1 N A and { τ j B } j = 1 N B , we have { ϕ ( τ j , S S B A ) } j = 1 N and { ϕ ( τ j , S S B B ) } j = 1 N , where
ϕ ( τ j , S S B A ) = ϕ i A + F i o r b i t , A ( τ j t i ) + ϕ ~ ( τ j , S S B A ) ,
ϕ τ j , S S B B = ϕ i B + F i o r b i t , B τ j t i + ϕ ~ τ j , S S B B .
As shown in Section 2.2, F i o r b i t , A and F i o r b i t , B can be estimated by (17). In Equation (19), ϕ i A and ϕ i B can be estimated by (18), and Δ ϕ i A , B can be calculated by (24). However, when estimating ϕ i A and ϕ i B , respectively, we have to solve the optimization problem shown in (18) twice, and the template is indispensable. To improve computational efficiency and eliminate dependence on the template, we proposed a direct estimation method for Δ ϕ i A , B .
Similar to (19), we define the cost function for estimating Δ ϕ i A , B as follows:
J Δ ϕ i A , B = l = 0 M b λ E F A ( φ l ; F ^ I orbit , A ) λ E F B ( φ l ; F ^ I orbit , B ) 2 .
Therefore, Δ ϕ i A , B can be estimated by
Δ ϕ ^ i A , B = arg min J ( Δ ϕ i A , B ) Δ ϕ i A , B [ 0 , 1 ) .
As shown in (32) and (33), Δ ϕ i A , B can be estimated by solving (33) once and is independent from the template of pulsars.

3.3. Procedure of Pulsar-Based Navigation Using Pulse Phase Delay

In this section, we discuss the procedure of pulsar-based navigation using pulse phase delay between two spacecraft. We assume that the true position and velocity of spacecraft A is known, and the position and velocity of spacecraft B can be estimated by measuring pulse phase delay between spacecraft A and B. By measuring Δ ϕ i A , B , we can obtain the distance of the direction of the pulsar between spacecraft A and spacecraft B. Using one measurement, we can determine a ring of coordinate triples for spacecraft B. Moreover, given that the dynamic model of spacecraft B can be modeled in advance, we use an unscented Kalman filter (UKF) to estimate the state of spacecraft B using the measurement and dynamic model. The procedure of pulsar-based navigation using pulse phase delay is shown in Algorithm 1. At first, we apply (30) and (31) to { τ j A } j = 1 N A and { τ j B } j = 1 N B and obtain { ϕ ( τ j , S S B A ) } j = 1 N and { ϕ ( τ j , S S B B ) } j = 1 N , respectively. Then, we estimate F ^ i o r b i t , A , F ^ i o r b i t , B by (17) and estimate Δ ϕ i A , B by (33). The state of spacecraft B is then estimated using UKF.
Algorithm 1. Procedure of pulsar-based navigation using pulse phase delay
Input:  x A ( t i ) , { τ j A } j = 1 N A , { τ j B } j = 1 N B  i = 1,2,…I
1. Initialization: Set initial state guess x ~ B ( t 1 ) .
2. for i = 1, 2,…I do
3.           Δ x ~ i A , B = x A ( t i ) x ~ B ( t i )
4.        Predict the state of spacecraft B at { τ j B } j = 1 N B by propagating x ~ B t i .
5.        Apply (30) and (31) to { τ j A } j = 1 N A and { τ j B } j = 1 N B .
6.        Estimate F ^ i o r b i t , A and F ^ i o r b i t , B by (17).
7.        Recover empirical profile λ E F A and λ E F B by F ^ i o r b i t , A and F ^ i o r b i t , B , respectively.
8.        Calculate y i , A B by (28) and (33).
9.        Estimate x ^ B ( t i ) by UKF.
10.      Predict x ~ B ( t i + 1 ) by propagating x ^ B ( t i ) .
11. end for
Output:  x ^ B ( t i )

4. Experiments

In this section, we used the simulation data of Crab pulsar and real data of Crab pulsar derived from NICER and Insight-HXMT to verify the proposed pulsar-based navigation method using pulse phase delay.

4.1. Experiments with Simulation Data

Here, the Crab pulsar with simulation parameters shown in Table 1 was taken as the observed pulsar. Figure 2 shows the template of Crab pulsar. The initial orbit elements of spacecraft A and spacecraft B are shown in Table 2. Both spacecraft were assumed to observe Crab pulsar with observation t o b s = 100   s . In this section, we assumed that the true position and velocity of spacecraft A was known, and the observation data of Crab pulsar were used to estimate state of spacecraft B with Algorithm 1. The initial state errors of spacecraft B were set as [5 km, 5 km, 5 km] and [5 m/s, 5 m/s, 5 m/s].
The initial state of spacecraft A and B was used to recover empirical profiles with { ϕ ( τ j , S S B H X M T ) } and { ϕ ( τ j , S S B N I C E R ) } via (13). The recovered empirical profiles are shown in Figure 3. Due to the distance between spacecraft A and B, pulse phase of Spacecraft A and B in the same exposure were different. Using (33) to estimate Δ ϕ ^ i A , B , the value of estimated pulse phase delay was 0.2650, and Δ ϕ i A , B , calculated by (29), was 0.2652. It can be seen that Δ ϕ ^ i A , B was close to Δ ϕ i A , B , computed via (29). Therefore, the pulse phase delay can reflect the distance between spacecraft A and B.
Figure 4 and Figure 5 show the navigation result of spacecraft B. As can be seen, the position and velocity gradually reduced to below 1 km and 1 m/s, respectively, thus verifying the feasibility of proposed navigation method using pulse phase delay.

4.2. Experiments with Real Data

We use real data of Crab pulsar from LE (Low Energy X-ray Telescope) on Insight-HXMT and NICER on 13 March 2018 and 14 March 2018 (Obs_ID: P0111605054 (Insight-HXMT), 1013010125 (NICER), and 1013010126 (NICER)). The start and finish time of real data are shown in Table 3.
To verify the proposed pulsar-based navigation method using pulse phase delay, we chose the real data when Insight-HXMT and NICER simultaneously observe the Crab pulsar. The exposures on Crab pulsar from NICER and Insight-HXMT are shown in Figure 6. As shown in the figure, 13 March 2018 (MJD:58190) Insight-HXMT and NICER did not observe Crab pulsar simultaneously. On 14 March 2018 (MJD:58191), there are nine exposures that meet the criteria. We select these nine exposures and align the start time of each exposure. The start times of the selected exposures are shown in Table 4.
Given that there are only nine exposures in which Insight-HXMT and NICER simultaneously observe Crab pulsar, we verified the proposed pulsar-based navigation method by comparing pulse phase delay calculated by (28) and that estimated by (33). With the pulse phases of events received by Insight-HXMT and NICER, { ϕ ( τ j , S S B H X M T ) } and { ϕ ( τ j , S S B N I C E R ) } were calculated by (20) and (21) using the real position of NICER and Insight-HXMT. The empirical profiles of NICER and Insight-HXMT in each exposure were then recovered using { ϕ ( τ j , S S B H X M T ) } and { ϕ ( τ j , S S B N I C E R ) } via (13). The empirical profiles in each exposure are shown in Figure 7. Given the distance between NICER and Insight-HXMT in the direction of pulsar, the pulse phases of NICER and Insight-HXMT in the same exposure were different. In accordance with (24), the pulse phase delay Δ ϕ ^ i H X M T , N I C E R reflects the distant between NICER and Insight-HXMT in the direction of the pulsar.
Next, we used the estimation error of Δ ϕ ^ i H X M T , N I C E R   to reflect the measurement error of the distance between Insight-HXMT and NICER. The pulse phase delays Δ ϕ ^ i H X M T , N I C E R , i = 1 , , 9 were estimated by (33) using the empirical profiles of NICER and Insight-HXMT shown in Figure 6. The real pulse phase delays were calculated by (28) using the real positions of NICER and Insight-HXMT. Furthermore, the estimation error of pulse phase delay was defined as the difference between estimated pulse phase delay and real pulse phase delay. Table 5 shows the estimation errors of the pulse phase delay of nine exposures. The results revealed that the estimated pulse phase delay can accurately reflect the distance between Insight-HXMT and NICER, thus verifying the feasibility of the proposed pulsar-based navigation using the pulse phase delay.
Finally, we compared the performance of the separate estimation method for pulse phase delay used in (19) and the proposed direct estimation method. In the proposed direct estimation method, pulse phase delays Δ ϕ ^ i H X M T , N I C E R ,   i = 1 , , 9 were directly estimated by (33) using the empirical profiles of Insight-HXMT and NICER in each exposure. In the respective estimation method, Δ ϕ ^ i H X M T and Δ ϕ ^ i N I C E R were estimated, respectively, by (19) using empirical profiles and templates, after which Δ ϕ ^ i H X M T , N I C E R ,   i = 1 , , 9 were then calculated by (24). The evaluation criterion of the estimation performance is defined as RMS (root mean square) error of the pulse phase delay. We defined CPU time cost as evaluation criterion of computational cost. The computation environment contained Intel Core [email protected] GHZ (Intel, Santa Clara, CA, USA) and python 3.8. Table 6 shows the RMS error and mean CPU time cost of the proposed direct estimation method and the separate estimation method for pulse phase delay. The results showed that the estimation performance of the proposed direct estimation method was close to that of respective estimation method. Furthermore, the CPU time cost of the proposed direct estimation method was about 36.99% lower than that of respective estimation method.

5. Conclusions

This study derived a pulsar-based navigation method using pulse phase delay between spacecraft. A direct estimation method for pulse phase delay is proposed to eliminate reliance on templates of pulsars for pulse phase delay estimation. The proposed method is verified with the simulation data of the Crab pulsar and the real data of the Crab pulsar form Insight-HXMT and NICER. Real data experiment results show that the pulse phase delay estimated by the real data of the Crab pulsar can accurately reflect the distance between Insight-HXMT and NICER. Moreover, compared with the estimation method that estimates the pulse phase of two spacecraft, the proposed direct estimation method can estimate pulse phase delay between two spacecraft by running the pulse phase estimation algorithm just once. The experimental results show that the estimation performance of the proposed direct estimation method is close to that of the respective estimation method and that the CPU time cost of the proposed method is about 36.99% lower. Notably, the transform model shown as (2) only considers the Romer delay and Shapiro delay. An accurate transform model that considers the Einstein delay and high-order terms of relativistic effects should be considered in further studies to provide an accurate navigation service for spacecraft near gravitational sources.

Author Contributions

Conceptualization, K.J. and H.Y. (Hong Yuan); methodology, K.J. and Y.W.; software, K.J. and Y.W.; validation, K.J., Y.W. and H.Y. (Hui Yang); formal analysis, K.J.; investigation, K.J.; resources, Y.W.; data curation, Y.W.; writing—original draft preparation, K.J.; writing—review and editing, K.J.; visualization, Y.W.; supervision, H.Y. (Hong Yuan); project administration, K.J. and H.Y. (Hong Yuan); funding acquisition, H.Y. (Hong Yuan). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available at [http://archive.hxmt.cn/iapi/opendownload/10014825/1L/zI4YsMOiXkMI3EO4jUpfM3ugAQ] (19 May 2024) [P0111605054] [NICER] [https://heasarc.gsfc.nasa.gov/FTP/nicer/data/obs/2018_03//1013010125/] (19 May 2024) [1013010125] [NICER] [https://heasarc.gsfc.nasa.gov/FTP/nicer/data/obs/2018_03//1013010126/] (19 May 2024) [1013010126].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fraser, C.T.; Ulrich, S. Adaptive extended Kalman filtering strategies for spacecraft formation relative navigation. Acta Astronaut. 2021, 178, 700–721. [Google Scholar] [CrossRef]
  2. Vasile, M.; Torre, F.; Serra, R.; Grey, S. Autonomous orbit determination for formations of Cubesats beyond LEO. Acta Astronaut. 2018, 153, 327–336. [Google Scholar] [CrossRef]
  3. Scharnagl, J.; Haber, R.; Dombrovski, V.; Schilling, K. NetSat—Challenges and lessons learned of a formation of 4 nano-satellites. Acta Astronaut. 2022, 201, 580–591. [Google Scholar] [CrossRef]
  4. Yang, C.; Zhang, H.; Gao, Y. Analysis of a neural-network-based adaptive controller for deep-space formation flying. Adv. Space Res. 2021, 68, 54–70. [Google Scholar] [CrossRef]
  5. Wang, X.; Gong, D.; Jiang, Y.; Mo, Q.; Kang, Z.; Shen, Q.; Wu, S.; Wang, D. A Submillimeter-Level Relative Navigation Technology for Spacecraft Formation Flying in Highly Elliptical Orbit. Sensors 2020, 20, 6524. [Google Scholar] [CrossRef]
  6. Zink, M.; Moreira, A.; Hajnsek, I.; Rizzoli, P.; Bachmann, M.; Kahle, R.; Fritz, T.; Huber, M.; Krieger, G.; Lachaise, M.; et al. TanDEM-X: 10 Years of Formation Flying Bistatic SAR Interferometry. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 3546–3565. [Google Scholar] [CrossRef]
  7. Wang, S.; He, C.; Gong, B.; Ding, X.; Yuan, Y. Cooperative Angles-Only Relative Navigation Algorithm for Multi-Spacecraft Formation in Close-Range. Comput. Model. Eng. Sci. 2023, 134, 121–134. [Google Scholar] [CrossRef]
  8. Liu, J.; Fang, J.; Yang, Z.; Kang, Z.-W.; Wu, J. X-ray pulsar/Doppler difference integrated navigation for deep space exploration with unstable solar spectrum. Aerosp. Sci. Technol. 2015, 41, 144–150. [Google Scholar] [CrossRef]
  9. Wang, Y.; Zheng, W.; Sun, S.; Li, L. X-ray pulsar-based navigation using time-differenced measurement. Aerosp. Sci. Technol. 2014, 36, 27–35. [Google Scholar] [CrossRef]
  10. Wang, Y.; Zheng, W.; Zhang, S.; Ge, M.; Li, L.; Jiang, K.; Chen, X.; Zhang, X.; Zheng, S.; Lu, F. Review of X-ray pulsar spacecraft autonomous navigation. Chin. J. Aeronaut. 2023, 36, 44–63. [Google Scholar] [CrossRef]
  11. Chen, P.-T.; Bayard, D.S.; Sharrow, R.F.; Majid, W.A.; Dunst, B.A.; Speyer, J.L. A Gaussian Sum Filter for Pulsar Navigation: Processing Single Photon Arrival Time Measurements. IEEE Trans. Contr. Syst. Technol. 2023, 31, 2499–2514. [Google Scholar] [CrossRef]
  12. Sheikh, S.I.; Pines, D.J.; Ray, P.S.; Wood, K.S.; Lovellette, M.N.; Wolff, M.T. Spacecraft Navigation Using X-Ray Pulsars. J. Guid. Control Dyn. 2006, 29, 49–63. [Google Scholar] [CrossRef]
  13. Winternitz, L.B.; Hassouneh, M.A.; Mitchell, J.W.; Price, S.R.; Yu, W.H.; Semper, S.R.; Ray, P.S.; Wood, K.S.; Arzoumanian, Z.; Gendreau, K.C. SEXTANT X-ray Pulsar Navigation Demonstration: Additional On-Orbit Results. In Proceedings of the 2018 SpaceOps Conference, Marseille, France, 28 May–1 June 2018. [Google Scholar]
  14. Zheng, S.J.; Zhang, S.N.; Lu, F.J.; Zheng, W.; Chen, X.; Zheng, S.; Lu, F. In-orbit demonstration of X-ray pulsar navigation with the Insight-HXMT satellite. Astrophys. J. Suppl. Ser. 2019, 244, 1. [Google Scholar] [CrossRef]
  15. Zheng, S.J.; Zhang, S.N.; Lu, F.J.; Wang, W.B.; Gao, Y.; Li, T.P.; Song, L.M.; Ge, M.Y.; Han, D.W.; Chen, Y.; et al. Fast On-Orbit Pulse Phase Estimation of X-Ray Crab Pulsar for XNAV Flight Experiments. IEEE Trans. Aerosp. Electron. Syst. 2022, 59, 3395–3404. [Google Scholar]
  16. Wang, Y.; Zheng, W.; Ge, M.; Zheng, S.; Zhang, S.; Lu, F. Use of Statistical Linearization for Nonlinear Least-Squares Problems in Pulsar Navigation. J. Guid. Control Dyn. 2023, 46, 1850–1856. [Google Scholar] [CrossRef]
  17. Xiong, K.; Wei, C.L.; Liu, L.D. The use of X-ray pulsars for aiding navigation of satellites in constellations. Acta Astronaut. 2009, 64, 427–436. [Google Scholar]
  18. Zhang, H.; Jiao, R.; Xu, L.P.; Xu, C.; Mi, P. Formation of a Satellite Navigation System Using X-Ray Pulsars. Publ. Astron. Soc. Pac. 2019, 131, 045002. [Google Scholar]
  19. Emadzadeh, A.A.; Speyer, J.L. X-Ray Pulsar-Based Relative Navigation using Epoch Folding. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2317–2328. [Google Scholar] [CrossRef]
  20. Zhou, Q.Y.; Ji, J.F.; Ren, H.F. Timing equation in X-ray pulsar autonomous navigation. Acta Phys. Sin. 2013, 62, 139701. [Google Scholar] [CrossRef]
  21. Lyne, A.G.; Pritchard, R.S.; Graham Smith, F. 23 Years of Crab Pulsar Rotational History. Mon. Not. R. Astron. Soc. 1993, 265, 1003–1012. [Google Scholar] [CrossRef]
  22. Edwards, R.T.; Hobbs, G.B.; Manchester, R.N. TEMPO2, a new pulsar timing package—II. The timing model and precision estimates. Mon. Not. Roy. Astron. Soc. 2006, 372, 1549–1574. [Google Scholar] [CrossRef]
  23. Emadzadeh, A.A.; Speyer, J.L. On Modeling and Pulse Phase Estimation of X-ray Pulsars. IEEE Trans. Signal Process 2010, 58, 4484–4495. [Google Scholar] [CrossRef]
  24. Wang, Y.; Zheng, W. Pulse Phase Estimation of X-ray Pulsar with the Aid of Vehicle Orbital Dynamics. J. Navig. 2016, 69, 414–432. [Google Scholar] [CrossRef]
  25. Wang, Y.; Zhang, W. Pulsar phase and Doppler frequency estimation for XNAV using on-orbit epoch folding. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 2210–2219. [Google Scholar] [CrossRef]
  26. Robert, G.M. Fundamentals of Astrodynamics and Applications. J. Guid. Control Dyn. 1998, 21, 803–812. [Google Scholar]
  27. Wang, Y.S.; Wang, Y.D.; Jiang, K.; Zheng, W.; Song, M. Adaptive Grid Search Based Pulse Phase and Doppler Frequency Estimation for XNAV. IEEE Trans. Aerosp. Electron. Syst. 2024, 60, 3707–3717. [Google Scholar] [CrossRef]
  28. Sun, H.; Su, J.; Deng, Z.; Shen, L.; Bao, W.; Li, X.P.; Li, L.; Su, Z.; Wang, W. Grouping bi-chi-squared method for pulsar navigation experiment using observations of Rossi X-ray timing explorer. Chin. J. Aeronaut. 2022, 36, 386–395. [Google Scholar] [CrossRef]
  29. Winternitz, L.M.B.; Hassouneh, M.A.; Mitchell, J.W.; Valdez, J.E.; Price, S.R.; Semper, S.R.; Yu, W.H.; Ray, P.S.; Wood, K.S.; Arzoumanian, Z.; et al. X-ray pulsar navigation algorithms and testbed for SEXTANT. In Proceedings of the 2015 IEEE Aerospace Conference, Big Sky, MT, USA, 7–14 March 2015. [Google Scholar]
Figure 1. Diagram of XNAV for spacecraft in formation.
Figure 1. Diagram of XNAV for spacecraft in formation.
Applsci 14 06401 g001
Figure 2. Template of the Crab pulsar.
Figure 2. Template of the Crab pulsar.
Applsci 14 06401 g002
Figure 3. Empirical profiles of spacecraft A and B.
Figure 3. Empirical profiles of spacecraft A and B.
Applsci 14 06401 g003
Figure 4. Position estimation error of spacecraft B.
Figure 4. Position estimation error of spacecraft B.
Applsci 14 06401 g004
Figure 5. Velocity estimation error of spacecraft B.
Figure 5. Velocity estimation error of spacecraft B.
Applsci 14 06401 g005
Figure 6. Exposure on Crab pulsar from Insight-HXMT and NICER.
Figure 6. Exposure on Crab pulsar from Insight-HXMT and NICER.
Applsci 14 06401 g006
Figure 7. Empirical profiles of each exposure.
Figure 7. Empirical profiles of each exposure.
Applsci 14 06401 g007aApplsci 14 06401 g007b
Table 1. Simulation parameters of Crab pulsar [29].
Table 1. Simulation parameters of Crab pulsar [29].
ParametersValue
Period [ms]33.4
α [ph/s]660
β [ph/s]13,860.2
Table 2. Orbital elements of spacecraft A and B.
Table 2. Orbital elements of spacecraft A and B.
Orbital ElementsSpacecraft ASpacecraft B
Inclination/°51.699563.2
Semi-major axis/km6781.58886922.0854
Eccentricity8.63 × 10−44.76 × 10−4
Right ascension of the ascending node/°137.2927252.4538
Argument of perigee/°124.2185121.7897
True anomaly/°183.5224279.4325
Table 3. Start and finish time of real data.
Table 3. Start and finish time of real data.
Obs_IDStart Time [UTC]Finish Time [UTC]
P011160505413 March 2018 03:38:2414 March 2018 18:02:10
101301012513 March 2018 02:58:1513 March 2018 23:02:12
101301012613 March 2018 00:32:4914 March 2018 14:58:16
Table 4. The start times of the selected nine exposures.
Table 4. The start times of the selected nine exposures.
NumberStart Time [UTC]
114 March 2018 02:06:27
214 March 2018 03:40:51
314 March 2018 05:14:53
414 March 2018 06:43:16
514 March 2018 08:16:06
614 March 2018 09:50:46
714 March 2018 11:22:42
814 March 2018 12:58:12
914 March 2018 14:34:19
Table 5. Estimation errors of the pulse phase delay.
Table 5. Estimation errors of the pulse phase delay.
Exposure NumberEstimation Error [Cycle]
12.63 × 10−3
2−2.82 × 10−3
32.80 × 10−3
43.89 × 10−4
5−9.16 × 10−4
63.39 × 10−5
7−2.51 × 10−3
8−2.44 × 10−4
9−1.42 × 10−4
Table 6. RMS error and mean CPU time cost for different pulse phase delay estimation methods.
Table 6. RMS error and mean CPU time cost for different pulse phase delay estimation methods.
MethodRMS Error [cycle]Mean CPU Time Cost [s]
Proposed direct estimation method1.83 × 10−30.0092
Respective estimation method1.71 × 10−30.0146
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiang, K.; Wang, Y.; Yang, H.; Yuan, H. X-ray Pulsar-Based Navigation Using Pulse Phase Delay between Spacecraft and Verification with Real Data. Appl. Sci. 2024, 14, 6401. https://doi.org/10.3390/app14156401

AMA Style

Jiang K, Wang Y, Yang H, Yuan H. X-ray Pulsar-Based Navigation Using Pulse Phase Delay between Spacecraft and Verification with Real Data. Applied Sciences. 2024; 14(15):6401. https://doi.org/10.3390/app14156401

Chicago/Turabian Style

Jiang, Kun, Yusong Wang, Hui Yang, and Hong Yuan. 2024. "X-ray Pulsar-Based Navigation Using Pulse Phase Delay between Spacecraft and Verification with Real Data" Applied Sciences 14, no. 15: 6401. https://doi.org/10.3390/app14156401

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop