Next Article in Journal
Approach for Cost Functions for the Use in Trade-Off Investigations Assessing the Environmental Impact of a Future Energy Efficient European Aviation
Next Article in Special Issue
A Beam Search-Based Channel Allocation Method for Interference Mitigation of NGSO Satellites with Multi-Beam Antennas
Previous Article in Journal
Multi-Objective Design of a Distributed Ducted Fan System
Previous Article in Special Issue
Constrained Motion Planning of 7-DOF Space Manipulator via Deep Reinforcement Learning Combined with Artificial Potential Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination

1
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China
2
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(3), 166; https://doi.org/10.3390/aerospace9030166
Submission received: 11 February 2022 / Revised: 11 March 2022 / Accepted: 12 March 2022 / Published: 18 March 2022
(This article belongs to the Special Issue Recent Advances in Spacecraft Dynamics and Control)

Abstract

:
Cooperative orbit determination (OD) using inter-spacecraft optical measurements is an important technology for space constellation missions. In this paper, the observability of a two-spacecraft cooperative OD system is investigated. The influence of geometric configuration on the observability is analyzed, and two special unobservable configurations are identified. Then, an approach to improve the observability by involving an additional spacecraft is proposed. Comparative analysis of system observability shows that an extra spacecraft in the system could change the coplanar and symmetric configuration and improve the observability of the cooperative OD system. Monte-Carlo simulations are carried out, and results verify the observability improvement conclusion.

1. Introduction

Autonomous orbit determination (OD) is one of the basic technologies to ensure the completion of space missions. With the dramatic increase of the number of spacecraft and the development of space cooperative behaviors, such as space constellation and formation missions [1,2], OD based on inter-spacecraft measurement becomes an important feature [3,4,5,6]. The techniques for estimating the absolute orbits of both the observer and the target using the inter-spacecraft measurements, which is called the cooperative OD, are increasingly desired.
The work on cooperative OD can be traced back to the 1980s, when Markley and Psiaki [7] proposed the OD method for two spacecraft using inter-spacecraft relative position measurements. The relative position measurement includes both distance and direction information, which requires several types of sensors, such as a laser sensor measuring relative distance and an optical camera measuring the relative direction [8]. With the miniaturization development of spacecraft, the onboard sensors are limited [9,10,11]. Another two widely used methods for cooperative OD are based on range-only measurements [12,13,14,15,16] and line-of-sight (LOS) measurements [17,18,19]. Ranging information can be obtained through radio, radar, or laser, and thus range measurements are hardly affected by the complex space environment and low cost [20]. However, the lack of angle information inevitably leads to overall rotation of the constellation, and absolute OD is unavailable [21,22]. The LOS measurement can be easily obtained using optical cameras with the advantages of long observation arc, high measurement accuracy and long detection range [23,24,25]. Moreover, the optical measurement, including the right ascension and the declination [26,27], can provide more information than the relative range. Therefore, angle-only measurements are of great value in the cooperative OD problem.
Regarding the cooperative OD problem based on inter-spacecraft optical measurements, observability analysis is particularly important [28,29,30,31,32]. It has been demonstrated that an OD system including one observer and one target is observable using inertial LOS measurements [33,34,35]. The inertial LOS measurements can be obtained by combining the relative LOS information and attitude of the observer in the inertial coordinate. The cooperative OD system using inertial LOS measurement is first analyzed by Yim et al. [33]. The conclusion that the two spacecrafts are observable was obtained through numerical simulations. Hu further investigated the observability of an inertial optical cooperative OD system with three spacecraft using the Lie-derivative criterion [34,35]. It was declared in Hu’s work that the three-spacecraft system is third-order locally, weakly observable. However, the previous conclusion does not hold in some cases. Specifically, the inertial cooperative OD system is unobservable under certain orbital geometric configurations, and observability of the OD system is poor when configuration is close to those special cases. Thus, it is necessary to further investigate the influence of orbital configuration on observability of a cooperative OD system.
This paper investigates the observability of the cooperative optical OD system. Differently from the previous studies of this problem focusing on the measurements and dynamics, this paper handles the problem in view of the orbit configuration. The contributions of this paper are summarized as follows. Firstly, the observability of a two-spacecraft optical cooperative OD system with different configurations is analyzed based on the observability matrix (OM). Several unobservable configurations are dug out and the corresponding unobservable elements are discussed. Secondly, an observability improvement strategy for the unobservable configurations is proposed by adding an additional observer. Simulations are performed to verify the enhancement effect on observability, and OD accuracy is analyzed and compared in conditions with and without the additional observer.
This paper is organized as follows: In Section 2, the observability of the inertial optical cooperative OD is analyzed. In Section 3, the observability improvement method is proposed and studied. Numerical simulations are implemented for comparisons through the Monte-Carlo method for accuracy analysis.

2. Observability Analysis of the Two-Spacecraft OD System

2.1. System Description

2.1.1. State Model

The cooperative OD system includes two spacecraft orbiting around Earth, known as the observer spacecraft (SO) and the target spacecraft (ST), with the Keplerian orbit elements given as:
E O = [ a O ,   e O ,   i O ,   ω O ,   Ω O ,   n O ] E T = [ a T ,   e T ,   i T ,   ω T ,   Ω T ,   n T ]
where the subscripts O and T represent the two spacecraft, SO and ST, respectively. The observation spacecraft can actively observe the target with onboard camera, while the target spacecraft is passively observed and is usually a non-cooperative target. For each spacecraft, the orbit elements denote semi-major axis, eccentricity, inclination, argument of the periapsis, longitude of the ascending node, and true anomaly, respectively. The goal of a cooperative OD system is to determine the orbit elements of both SO and ST. The system state vector x is given by:
x = [ a O ,   e O ,   i O ,   ω O ,   Ω O ,   n O , a T ,   e T ,   i T ,   ω T ,   Ω T ,   n T ] T
Considering point-mass two-body dynamics without any non-gravitational perturbation [36], the size, shape, and spatial position of each orbit described by the first five parameters of E O and E T are time-invariant [35,36,37,38]. The true anomaly, which represents the position of the spacecraft in orbit, is time varying. Then, the state model of the cooperative OD system is given by:
x ˙ = f ( x ) = [   0 ,   0 ,   0 ,   0 ,   0 ,   n ˙ O ,   0 ,   0 ,   0 ,   0 ,   0 ,   n ˙ T ] T
n ˙ O = μ a O ( 1 e O 2 ) ( 1 + e O cos ( n O ) ) 2 a O ( 1 e O 2 ) n ˙ T = μ a T ( 1 e T 2 ) ( 1 + e T cos ( n T ) ) 2 a T ( 1 e T 2 )

2.1.2. Measurement Model

Space-based optical measurements use optical cameras to obtain images of the surrounding environment, and the measurement information is then passed through a signal processor to extract the target. For a distant space target, its image is often seen as a light spot on the image plane, so that we can obtain only the target’s line-of-sight (LOS) information relative to the observation spacecraft. Combining the attitude of SO received via star sensors, the inertial LOS measurement can be obtained.
The inertial LOS measurement model is given as:
h ( x ) = r T r O r T r O + υ
where r O and r T represent the inertial position vectors of SO and ST at the time of measurement regardless of the time it takes for the light to travel, and υ is the measurement noise. The LOS measurement is defined in the Earth centered inertial frame ( O E X E Y E Z E ), and a schematic diagram of the measurement model is shown in Figure 1.
To analyze the effect of orbital configurations, the LOS measurement should be expressed as a function of the system state variables, and we should describe the inertial position vectors by the orbit elements. The coordinate transformation is the same for each spacecraft, so the subscripts are omitted during the derivation of the inertial position vector.
The near-focal coordinate system ( O x ¯ y ¯ z ¯ ) describes the natural properties of an orbit [39], which takes the focus of the orbit (center of the earth for earth-orbiting spacecraft) as the origin, and the orbital plane is the reference plane. The x ¯ -axis points to the perigee direction from the center of the earth, the y ¯ -axis is obtained by rotating the x ¯ -axis by 90° along the motion direction in the orbital plane, and the O x ¯ y ¯ z ¯ system is a right hand coordinate system.
In the O x ¯ y ¯ z ¯ system, the position vector of a spacecraft can be expressed by the orbit elements as follows:
r x ¯ = a ( 1 e 2 ) 1 + e cos n [ cos n sin n 0 ]
where the subscript x ¯ represents the O x ¯ y ¯ z ¯ system, and r x ¯ is the position vector in O x ¯ y ¯ z ¯ system. The orthogonal transformation matrix from O E X E Y E Z E to O x ¯ y ¯ z ¯ is:
Q X x ¯ = T z ( ω ) T x ( i ) T z ( Ω )
where T z and T x are Euler orientation cosine matrixes.
Then the transformation matrix from O x ¯ y ¯ z ¯ to O E X E Y E Z E is the transpose of Q X x ¯ , that is:
Q x ¯ X = Q X x ¯ T = [ cos ( Ω ) sin ( Ω ) 0 sin ( Ω ) cos ( Ω ) 0 0 0 1 ] [ 1 0 0 0 cos ( i ) sin ( i ) 0 sin ( i ) cos ( i ) ] [ cos ( ω ) sin ( ω ) 0 sin ( ω ) cos ( ω ) 0 0 0 1 ]
The inertial position vector r can be described by the orbit elements as:
r = Q x ¯ X r x ¯ = Q x ¯ X a ( 1 e 2 ) 1 + e cos n [ cos n sin n 0 ]

2.1.3. Observability Matrix

The measurement model only reflects the relationship between the observation and the spacecraft state at a certain observation moment. The spacecraft state changes with time during the whole observation period, and its changing rule is described by the dynamic model. The observability matrix (OM) combines the measurement equation and the dynamics equation, and thus characterizes the relationship between observations and states over the entire observation period [18].
Using the method in [14], the OM of the two-spacecraft OD system, signed as M , is constructed based on the state equation in Equation (3) and the observation equation in Equation (5). The partial differential matrix of the measurement model at the observation time t i is given by:
H i ( x ) = [ h x ] i
For each observation time t i , H i ( x ) is a 3 by 12 matrix, and it should be mapped to the initial epoch t 0 through the state transformation matrix (STM), given as:
H ˜ i ( x ) = H i ( x ) Φ ( t i , t 0 )
where Φ ( t i , t 0 ) is the STM from t 0 to t i , and given as:
{ Φ ˙ ( t i , t 0 ) = [ f ( x 1 ) x 1 ] i Φ ( t i , t 0 ) Φ ( t 0 , t 0 ) = I 12 × 12
In the cooperative OD system, the STM has the following format:
Φ ( t i , t 0 ) = [ A O 0 6 × 6 0 6 × 6 A T ]
where A k ( k = O , T ) is derived as:
A k = [ 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 Φ a k n k ( t i , t 0 ) Φ e k n k ( t i , t 0 ) 0 0 0 Φ n k n k ( t i , t 0 ) ]
The elements in matrix A k shows the state transformations for specific orbit elements. For the time-invariant elements, the state transformations caused by themselves are equal to 1, showing as the diagonal elements in the first five rows. The state transformations caused by the irrelevant orbit elements are equal to 0. It is known from Equation (4) that n ˙ is relevant to a , e and n for each spacecraft, so the state transformations of n k are expressed in the differential equations as:
Φ ˙ a k n k ( t i , t 0 ) = n ˙ k a k + n ˙ k n k Φ a k n k ( t i , t 0 ) ,           Φ a k n k ( t 0 , t 0 ) = 0 Φ ˙ e k n k ( t i , t 0 ) = n ˙ k e k + n ˙ k n k Φ e k n k ( t i , t 0 ) ,           Φ e k n k ( t 0 , t 0 ) = 0 Φ ˙ n k n k ( t i , t 0 ) = n ˙ k n k     Φ n k n k ( t i , t 0 ) ,             Φ n k n k ( t 0 , t 0 ) = 1  
where Φ a k n k ( t i , t 0 ) describes the state transformation of n k (marked by the superscript) caused by a k (marked by the subscript). Similarly, Φ e k n k ( t i , t 0 ) and Φ n k n k ( t i , t 0 ) represent the state transformation of n k caused by e k and n k .
Then, the OM is updated with observation in the form of:
M = [ H ˜ 1 ( x ) H ˜ 2 ( x ) ]
For each observation time t i , H ˜ i ( x ) is also a 3 by 12 matrix. It can be seen that H ˜ i ( x ) makes up M over time, and each column of M corresponds to the orbital element.
Further mathematical analysis based on OM is implemented to investigate the observability of an OD system. In the first place, checking the rank condition of OM is necessary for qualitatively analyzing the observability of the system. A full rank OM means the OD system is observable, while when OM is not full rank, we can draw the conclusion that the system is not observable. The rank of OM represents the number of observable states or state combinations. Secondly, the quantitative analysis of system observability is important for checking the degree of system observability, which is generally based on the computation of the condition number (CN). CN is the ratio of the singular values of OM, given by:
C N = max ( σ ( O M ) ) min ( σ ( O M ) )
A larger CN indicates a less observable problem, and the estimation problem is ill-conditioned when its value exceeds a certain limit. An OD system is regarded as unobservable if the condition number is greater than 10 16 [32].

2.2. Observability of the Angle-Only Cooperative OD System

In this subsection, the influence of geometric configuration on the observability of the two-spacecraft system is discussed, and two unobservable configurations, including the symmetric case and the same circular orbit case, are analyzed. Numerical simulations are implemented to verify the observability results for different orbit configurations.
The obstruction of Earth is considered to ensure continuous visibility between the two spacecraft, and more complex environmental factors such as light conditions are not considered in this paper. The target spacecraft with an orbital altitude of 5000 km and the observation spacecraft with the orbital altitude greater than 4000 km is adopted to obtain continuous observation. The nominal orbital elements for the following three cases are shown in Table 1.

2.2.1. General Case

Assume that two spacecraft move on orbits with totally different elements that satisfy Equations (18)–(20):
Ω 1 Ω 2 2 k π   or   i 1 i 2
Ω 1 Ω 2 ( 2 k 1 ) π   or   i 1 i 2
a T a O ,   e T e O ,   e T 0 ,   e O 0
where k is an integer. Equations (18) and (19) mean that the two orbits are not coplanar, and Equation (20) reflects that they have different sizes and shapes. The orbits are shown in Figure 2.
In this case, the rank of the OM is calculated as 12, illustrating that all states in Equation (3) are observable, and the OD system is observable. The simulations of the OD problem are implemented by using the unscented Kalman filter (UKF). For each spacecraft in the OD system, the initial deviations of the triaxial position and triaxial velocity are set as 10 km and 1 m/s, respectively. Then, the initial covariance of state is defined as a diagonal matrix whose elements are the squares of the initial deviations. Then, the initial deviations of the orbital elements are obtained by converting state vectors to the six orbital elements, and the numerical values are shown in Table 2.
The standard deviation of measurement noise is set to 0.01 deg for LOS measurements according to the payload capabilities [40,41]. The process noise of the UKF is set to 10−12, and the simulation step is 60 s, depending on the measurement frequency. The initial simulation parameters are shown in Table 3.
Taking the case in which SO1 observes ST1 as a numerical example, the errors of orbit elements are shown in Figure 3. The estimation errors converge within four hours, and the sub-windows show the local magnifications of the estimation errors within 8–12 h. The estimation accuracies of semi-major axes are within 0.1 km, the accuracies of eccentricities are within 3 × 10−5, the accuracy of i reaches the level of 10−5, the estimated accuracies of Ω , ω and n reach the level of 10−4.
The orbit determination precisions are obtained through Monte-Carlo simulations [42]. The final standard deviations (STDs) and convergency ratio (CR) of the ST are shown in Table 4.
It can be seen that the states of target spacecraft can be determined in a general case. The STDs for triaxial position errors are within 0.1 km, and the convergency ratios are higher than 99%. The STDs for triaxial velocity errors are within 0.1 m/s, and the convergency ratios are higher than 96%.

2.2.2. Symmetric Case

In this case, the positions of the two spacecraft always remain symmetric with a stationary inertial plane. The restrictions in Equations (21) and (22) should be satisfied:
a T = a O ,   e T = e O ,   i T i O ,   n T = n O
{ ω T = ω O Ω T = Ω O   or   { | ω T ω O | = π | Ω T Ω O | = π
Equation (21) shows that the two orbits have the same semi-major axis and eccentricity, but different orbital inclinations. The restrictions of ω and Ω in Equation (22) indicate that the two orbits are mirror image symmetric in inertial space, and the same n means that the two spacecraft are symmetric with the same inertial plane.
Let SO2 observe ST1, all orbit elements of the two orbits are different except for the inclinations. Numerical analysis shows that the rank of OM decreases to 6, meaning that there are 6 observable elements or element combinations. The OM can be simplified to the reduced row echelon form through Gauss–Jordan elimination method. Only six rows of the simplified matrix are valid, and the elements in the remaining rows are zero. The first six rows are shown in Table 5, and the numerical values that are not listed in the table are all equal to 0.
Each row of the simplified matrix has more than one number, indicating that no independent element is observable, and there are six pairs of observable element combinations. Taking the first row as an example, there are two numbers in the columns that are corresponding to a O and a T , showing that the combination ( a O a T ) is observable. Similarly, the 2nd, 3rd and 6th rows indicate that the combinations ( e O e T ) , ( i O + i T ) and ( n O n T ) are observable. The numbers in the 4th and 5th row show that the remaining angles, Ω O , Ω T , ω O and ω T , make up two observable combinations in an irregular form.
The estimation errors of orbit elements over time are shown in Figure 4. It can be known that none of the error’s convergence, which verifies the observability conclusion that no independent element is observable. The estimation errors of the combinations ( a O a T ) , ( e O e T ) , ( i O + i T ) and ( n O n T ) in Figure 5 are convergent, which proves the observability results for orbit element combinations.
From the point of view of OD, this symmetrical configuration is unfavorable. The LOS direction is always perpendicular to the symmetry plane of the two spacecraft, and therefore doesn’t contain the information related to the semi-major axis, eccentricity and the true anomaly. Moreover, the LOS vectors before and after the cross epoch are opposed to each other. In this case, the drastic changes of the measurements cause the estimations to diverge (as shown in Figure 4). Although this symmetric case rarely appears in practice, it is worth investigating because the observability and convergence will also be poor in configurations close to the symmetric case.

2.2.3. Same Circular Orbit Case

For two spacecraft on the coplanar orbits, the elements should satisfy one of the following equations:
i 1 = i 2   and   Ω 1 Ω 2 = 2 k π
i 1 = i 2   and   Ω 1 Ω 2 = ( 2 k 1 ) π
In this case, the two spacecraft move on the same circular orbit with only different true anomalies, and the remaining elements should satisfy:
a T = a O ,   e T = e O = 0 ,   n T n O
Let SO3 observe ST2. The rank of OM becomes 6, and the reduced echelon form is shown in Table 6. For the 2nd and 6th row, there is only one number in each row, indicating that the corresponding elements, e O and e T , are independently observable. The numbers corresponding to the column of a T are in the 1st and 4th row. There is another number corresponding to a O in the 1st row, and there are 7 numbers in the 4th row. Meanwhile, the numbers corresponding to i T and Ω T also appear in the 3rd, 4th and 5th rows. Therefore, the other 4 observable variables are composed of the semi-major axes and the eight angles in a combinatorial form.
The estimation errors of orbit elements over time are shown in Figure 6. It can be seen that the errors of a ,   i ,   ω ,   Ω and n are not convergent, proving that the 10 elements are unobservable. For the errors of e , there are tendencies of convergence, but the accuracies are not high due to the influence of other orbit elements.

3. Improvement Approach for Cooperative Optical Orbit Determination

It can be drawn from the previous section that the observability of a two-spacecraft cooperative OD system is limited by the relative geometric configuration. To solve this problem, an observability improvement approach is proposed by adding an additional observer into the cooperative OD system to obtain more measurement information. In this section, the optical cooperative OD system with an additional observer is modeled, and the system observability is analyzed.

3.1. System Description for the Cooperative OD System with an Additional Observer

In this cooperative OD scenario, each observation spacecraft measures the LOS from itself to the target spacecraft. The measurement model is given as follows:
h = [ h 1 T h 2 T ] T
where h 1 and h 2 are the LOS measurements taken by the two-observation spacecraft. The measurement model is shown in Figure 7, S O 1 and S O 2 are observation spacecraft.
The elements of the target and the two observers should be estimated, and there are six orbital elements for each spacecraft. Then the state variables become:
x = [ E O 1 ,   E O 2 ,   E T ] T
The OM of the cooperative OD system with an additional observer, defined as M , is constructed in the same way as in Equations (10)–(16), with the observation matrix H i ( x ) and the state transformation matrix Φ described as:
H i ( x ) = [ h x ] i
Φ ( t i , t 0 ) = [ A O 1 0 6 × 6 0 6 × 6 0 6 × 6 A O 2 0 6 × 6 0 6 × 6 0 6 × 6 A T ]
where the matrix A k ( k = O 1 ,   O 2 ,   T ) has the same form shown in Equations (14)–(16).
Then, we can get the intermediate matrix H ˜ i ( x ) , and M is given as:
M = [ H ˜ 1 ( x ) H ˜ 2 ( x ) ]
Similarly, the 18 columns of M correspond to the 18 orbital elements.

3.2. Observability Analysis for the Cooperative OD System with an Additional Observer

With an additional spacecraft involved in the OD system, not only the measurement information, but also the number of states, are increased. The observability of the system requires further analysis. In this part, three cases corresponding to the configurations in Sec. II are analyzed to show the observability improvement. The orbital elements of the additional spacecraft, signed as SO4, are shown in Table 7.

3.2.1. General Case

The first scenario considers two general non-coplanar observers. According to the observability result of the general case in a two-spacecraft system, either of the observation spacecraft can make up an observable OD system with the target spacecraft, so the cooperative OD system with both observers is completely observable. Therefore, it is feasible to estimate all orbital elements, although the dimension of the state variables increases.
Let SO1 and SO4 observe ST1, and the orbits are shown in Figure 8. The rank of OM is calculated to be 18, as expected, and the observability of cooperative OD system is enhanced, reflected in the fact that the numerical value of CN decreases to 6.7037 × 103. The errors in Figure 9 show that the accuracies of semi-major axes are around 20 m, the estimated errors of eccentricities reach the level of 10−6, the estimated accuracies of i and Ω reach the level of 10−6, and the estimated accuracies of ω and n reach the level of 10−4.
To validate the advantage of using an additional observer in the inertial LOS cooperative OD system, the OD accuracies of the three-spacecraft system proposed are analyzed. For the general case with two general observation spacecraft, the STD and CR results of 100 case Monte-Carlo simulations are shown in Table 8. Comparing with the Monte-Carlo results for two-spacecraft system in Table 4, the values of STDs are much smaller. The STDs for triaxial position errors are within 0.015 km, and the STDs for triaxial velocity errors are within 0.01 m/s. The convergency ratios are higher than 99%.
To compare the OD accuracies with and without the additional observer, the root-mean-square error (RMSE) of Monte-Carlo simulations are shown in Figure 10. The top two images are the RMSEs of the case without additional observer (i.e., two-spacecraft case), and the lower subgraphs show the RMSEs of the case with additional observer (i.e., three-spacecraft case). For the two-spacecraft case, the RMSEs of the position are smaller than 1.5 km, and the RMSEs of the velocity are smaller than 8 × 10−4 km/s. For the three-spacecraft case, the RMSEs of position are within 0.1 km, and the RMSEs of velocity are within 5 × 10−5 km/s. The OD accuracies of the three-spacecraft system are better. Moreover, the estimated errors of the three-spacecraft OD system converge after about 2 h, which are much faster than in the two-spacecraft system.

3.2.2. Symmetric Case

It has been proven that the cooperative OD system is unobservable with only an observer on the symmetric orbit with the target, and neither of the two orbits can be determined.
Let SO2 and SO4 observe ST1. The OM is calculated to be full rank, meaning that the system is observable and even the orbit of the symmetric observer is available. The cooperative OD system determines the orbits of SO1 and ST1 firstly, as in the general OD case, and then determines the orbit of SO3 based on this reference. The effect of symmetric configuration has been diminished but not eliminated, and the value of CN in this case is an order of magnitude larger than that in the general case, indicating that the system is still observable, but with weaker observability. The errors of orbital elements in Figure 11 show that the element accuracies can reach the same level, as in general case with an observation time longer than 3 h. It should be noted that the additional observer cannot be in special configurations, because the combining of singular cases is not able to break the previous unobservable configurations.
The OD accuracy results for the system with a symmetric observer and a general observer are shown in Table 8. It can be seen that the accuracies reach the same level as the case with two general observers.

3.2.3. Same Circular Observer Case

Let SO3 and SO4 observe ST. The value of CN increases sharply and reaches the level of 1016. The rank of OM is calculated to be 16, while the unobservable situation is caused by the circular orbits, which leads to the combination of ( ω O 3 + n O 3 ) and ( ω T 2 + n T 2 ) observable instead of the four individual elements. Therefore, it is proven that adding a general observer can make the same circular orbits observable.
The advantage of including two observers is that the observability of the OD system improves. The orbit configuration between each observation spacecraft and the target is different, so that the system can obtain more effective information. The analysis above shows that the combination of various orbital configurations can enrich measurement information to change the unobservable configurations, and thus improve the system observability. Same as in symmetric case, an additional observer on the same circular orbit or symmetric orbit with the target is not able to make the unobservable cooperative OD system observable.

4. Conclusions

In this paper, the observability of cooperative orbit determination system was analyzed based on numerical calculation of OM and analysis of the convergence of the orbit elements. The general case where the two spacecraft are on orbits with different elements is proven to be observable, while two unobservable orbit configurations are identified: (1) in the symmetric case, no independent element but six pairs of element combinations are observable; (2) in the same circular orbit case, there exist two observable individual elements and four element combinations. Involving an additional spacecraft has a remarkable effect on changing unobservable orbital configurations, and thus makes the OD system observable. For general cases, the numerical value of condition number decreases from the level of 106 in two-spacecraft system to the level of 103 in a three-spacecraft system. In terms of OD accuracy, the final STDs for triaxial position errors decreased from 0.1 km in a two-spacecraft system to within 0.015 km in a three-spacecraft system. The simulations show that an additional observer can improve the observability and OD accuracy of the cooperative OD system.

Author Contributions

Conceptualization, T.Q.; methodology, Y.L.; software, X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Natural Science Foundation of China, grant number 51827806, and The National Defense Basic Scientific Research Project, grant number JCKY2020903B002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, X.; Cheng, Y.; Qiao, D. An adaptive surrogate model-based fast planning for swarm safe migration along halo orbit. Acta Astronaut. 2022, 194, 309–322. [Google Scholar] [CrossRef]
  2. Wu, W.R.; Tang, Y.H.; Zhang, L.H.; Qiao, D. Design of communication relay mission for supporting lunar-farside soft landing. Sci. China Inf. Sci. 2018, 61, 040305. [Google Scholar] [CrossRef]
  3. Markley, F. Autonomous navigation using landmark and intersatellite data. In Proceedings of the Astrodynamics Conference, Seattle, WA, USA, 20–22 August 1984; pp. 1987–1990. [Google Scholar] [CrossRef]
  4. Psiaki, M.L. Autonomous orbit determination for two spacecraft from relative position measurements. J. Guid. Control Dyn. 1999, 22, 305–312. [Google Scholar] [CrossRef] [Green Version]
  5. Tang, Y.H.; Wu, W.R.; Qiao, D.; Li, X.Y. Effect of orbital shadow at an Earth--Moon Lagrange point on relay communication mission. Sci. China Inf. Sci. 2017, 60, 112301. [Google Scholar] [CrossRef] [Green Version]
  6. Li, X.Y.; Qiao, D.; Barucci, M.A. Analysis of equilibria in the doubly synchronous binary asteroid systems concerned with non-spherical shape. Astrodynamics 2018, 2, 133–146. [Google Scholar] [CrossRef]
  7. Psiaki, M.L. Absolute orbit and gravity determination using relative position measurements between two satellites. J. Guid. Control Dyn. 2011, 34, 1285–1297. [Google Scholar] [CrossRef] [Green Version]
  8. Li, X.Y.; Qiao, D.; Li, P. Frozen orbit design and maintenance with an application to small body exploration. Aerosp. Sci. Technol. 2019, 92, 170–180. [Google Scholar] [CrossRef]
  9. Li, X.Y.; Qiao, D.; Li, P. Bounded Trajectory design and self-adaptive maintenance control near non-synchronized binary systems comprised of small irregular bodies. Acta Astronaut. 2018, 152, 768–781. [Google Scholar] [CrossRef]
  10. Wang, Z.; Li, Z.; Wang, N.; Hoque, M.; Wang, L.; Li, R.; Zhang, Y.; Yuan, H. Real-Time Precise Orbit Determination for LEO between Kinematic and Reduced-Dynamic with Ambiguity Resolution. Aerospace 2022, 9, 25. [Google Scholar] [CrossRef]
  11. Colagrossi, A.; Lavagna, M. Fault Tolerant Attitude and Orbit Determination System for Small Satellite Platforms. Aerospace 2022, 9, 46. [Google Scholar] [CrossRef]
  12. Wang, J.; Butcher, E.A.; Lovell, T.A. Ambiguous relative orbits in sequential relative orbit estimation with range-only measurements. Acta Astronaut. 2018, 151, 626–644. [Google Scholar] [CrossRef]
  13. Christian, J.A. Relative navigation using only intersatellite range measurements. J. Spacecr. Rocket 2016, 54, 13–18. [Google Scholar] [CrossRef]
  14. Qin, T.; Qiao, D.; Macdonald, M. Relative orbit determination using only intersatellite range measurements. J. Guid. Control Dyn. 2019, 42, 703–710. [Google Scholar] [CrossRef]
  15. Qin, T.; Qiao, D.; Macdonald, M. Relative orbit determination for unconnected spacecraft within a constellation. J. Guid. Control Dyn. 2021, 44, 614–621. [Google Scholar] [CrossRef]
  16. Qin, T.; Macdonald, M.; Qiao, D. Fully Decentralized Cooperative Navigation for Spacecraft Constellations. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2383–2394. [Google Scholar] [CrossRef]
  17. Woffinden, D.C. Angle-only navigation for autonomous orbital rendezvous. J. Guid. Control Dyn. 2007, 30, 1455–1469. [Google Scholar] [CrossRef]
  18. Kaufman, E.; Lovell, T.A.; Lee, T. Nonlinear observability measure for relative orbit determination with angle-only measurements. J. Astronaut. Sci. 2016, 63, 60–80. [Google Scholar] [CrossRef]
  19. Ou, Y.; Zhang, H.B.; Li, B. Absolute orbit determination using line-of-sight vector measurements between formation flying spacecraft. Astrophys. Space Sci. 2018, 363, 76. [Google Scholar] [CrossRef]
  20. Cappuccio, P.; Notaro, V.; Ruscio, A.D. Report on First Inflight Data of BepiColombo’s Mercury Orbiter Radio Science Experiment. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 4984–4988. [Google Scholar] [CrossRef]
  21. Wang, Y.; Zhang, Y.; Qiao, D. Transfer to near-Earth asteroids from a lunar orbit via Earth flyby and direct escaping trajectories. Acta Astronaut. 2017, 133, 177–184. [Google Scholar] [CrossRef]
  22. Zhang, L. Development and Prospect of Chinese Lunar Relay Communication Satellite. Space Sci. Technol. 2021, 2021, 3471608. [Google Scholar] [CrossRef]
  23. Jia, B.; Pham, K.D.; Blasch, E. Cooperative space object tracking using space-based optical sensors via consensus-based filters. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 1908–1936. [Google Scholar] [CrossRef]
  24. Li, X.Y.; Warier, R.R.; Sanyal, A.K.; Qiao, D. Trajectory Tracking Near Small Bodies Using Only Attitude Control. J. Guid. Control Dyn. 2019, 42, 109–122. [Google Scholar] [CrossRef]
  25. Huang, J.; Lei, X.; Zhao, G.; Liu, L.; Li, Z.; Luo, H.; Sang, J. Short-Arc Association and Orbit Determination for New GEO Objects with Space-Based Optical Surveillance. Aerospace 2021, 8, 298. [Google Scholar] [CrossRef]
  26. Stokes, G.; Vo, C.; Sridharan, R. The space-based visible program. In Proceedings of the Space 2000 Conference and Exposition, Long Beach, CA, USA, 19–21 September 2000. [Google Scholar]
  27. Dong, L.; Yi, D. Initial Orbit Determination of Space Object Using Space-Based Optical Observations. In Proceedings of the 2012 International Conference on Computer Science and Electronics Engineering, Hangzhou, China, 23–25 March 2012. [Google Scholar] [CrossRef]
  28. Hermann, R.; Krener, A. Nonlinear controllability and observability. IEEE Trans. Autom. Control 1977, 22, 728–740. [Google Scholar] [CrossRef] [Green Version]
  29. Ham, F.M.; Brown, R.G. Observability, eigenvalues, and Kalman filtering. IEEE Trans. Aerosp. Electron. Syst. 1983, AES-19, 269–273. [Google Scholar] [CrossRef]
  30. Wang, Y.; Qiao, D.; Cui, P.Y. Design of optimal impulse transfers from the Sun–Earth libration point to asteroid. Adv. Space Res. 2015, 56, 176–186. [Google Scholar] [CrossRef]
  31. Molli, S.; Durante, D.; Gascioli, G.; Proietti, S. Performance analysis of a Martian polar navigation system. In Proceedings of the 72nd International Astronautical Congress, Dubai, United Arab Emirates, 25–29 October 2021. [Google Scholar]
  32. Woffinden, D.C.; Geller, D.K. Observability criteria for angle-only navigation. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 1194–1208. [Google Scholar] [CrossRef]
  33. Yim, J.; Crassidis, J.; Junkins, J. Autonomous orbit navigation of two spacecraft system using relative line of sight vector measurements. Adv. Astronaut. Sci. 2005, 1, 257. [Google Scholar]
  34. Hu, Y.; Sharf, I.; Chen, L. Three-spacecraft autonomous orbit determination and observability analysis with inertial angle-only measurements. Acta Astronaut. 2020, 170, 106–121. [Google Scholar] [CrossRef]
  35. Hu, Y.; Sharf, I.; Chen, L. Distributed orbit determination and observability analysis for satellite constellations with angle-only measurements. Automatica 2021, 129, 109626. [Google Scholar] [CrossRef]
  36. Stefano, I.D.; Cappuccio, P.; Iess, L. The BepiColombo solar conjunction experiments revisited. Class. Quantum Gravity 2022, 38, 0550002. [Google Scholar] [CrossRef]
  37. Cui, P.Y.; Qiao, D.; Cui, H.T. Target selection and transfer trajectories design for exploring asteroid mission. Sci. China Technol. Sci. 2010, 53, 1150–1158. [Google Scholar] [CrossRef]
  38. Li, X.Y.; Sanyal, A.K.; Warier, R.R.; Qiao, D. Landing of Hopping Rovers on Irregularly-shaped Small Bodies Using Attitude Control. Adv. Space Res. 2020, 65, 2674–2691. [Google Scholar] [CrossRef]
  39. Howard, D.C. Orbital Mechanics for Engineering Students, 4th ed.; Butterworth-Heinemann: Oxford, UK, 2021; pp. 59–61. [Google Scholar]
  40. Maskell, C.P. Sapphire: Canada’s Answer to Space-Based Surveillance of Orbital Objects. In Proceedings of the Advanced Maui Optical and Space Surveillance Technologies Conference, Wailea Beach, Maui, HI, USA, 16–19 September 2008. [Google Scholar]
  41. Li, X.Y.; Qiao, D.; Jia, F.D. Investigation of Stable Regions of Spacecraft Motion in Binary Asteroid Systems by Terminal Condition Maps. J. Astronaut. Sci. 2021, 68, 891–915. [Google Scholar] [CrossRef]
  42. Merel, M.H.; Mullin, F.J. Analytic Monte Carlo error analysis. J. Spacecr. Rocket. 1968, 5, 1304–1308. [Google Scholar] [CrossRef]
Figure 1. Orbits and inertial LOS Measurement.
Figure 1. Orbits and inertial LOS Measurement.
Aerospace 09 00166 g001
Figure 2. Orbits in general orbit configuration case.
Figure 2. Orbits in general orbit configuration case.
Aerospace 09 00166 g002
Figure 3. Estimation errors of orbit elements in general case.
Figure 3. Estimation errors of orbit elements in general case.
Aerospace 09 00166 g003
Figure 4. Estimation errors of orbit elements in the symmetric case.
Figure 4. Estimation errors of orbit elements in the symmetric case.
Aerospace 09 00166 g004
Figure 5. Estimation errors of orbit element combinations in the symmetric case.
Figure 5. Estimation errors of orbit element combinations in the symmetric case.
Aerospace 09 00166 g005
Figure 6. Estimation errors of orbit elements in the same circular orbit case.
Figure 6. Estimation errors of orbit elements in the same circular orbit case.
Aerospace 09 00166 g006
Figure 7. LOS measurement model for a cooperative OD system with an additional observer.
Figure 7. LOS measurement model for a cooperative OD system with an additional observer.
Aerospace 09 00166 g007
Figure 8. Orbits in general orbit configuration case with an additional observer.
Figure 8. Orbits in general orbit configuration case with an additional observer.
Aerospace 09 00166 g008
Figure 9. Errors of orbit elements in general case.
Figure 9. Errors of orbit elements in general case.
Aerospace 09 00166 g009
Figure 10. RMSE results of Monte-Carlo simulations.
Figure 10. RMSE results of Monte-Carlo simulations.
Aerospace 09 00166 g010
Figure 11. Errors of orbit elements in symmetric case.
Figure 11. Errors of orbit elements in symmetric case.
Aerospace 09 00166 g011
Table 1. Nominal orbital elements for observability analysis.
Table 1. Nominal orbital elements for observability analysis.
Spacecrafta/kmei/degΩ/degω/degn/degAnnotation
ST111,378.1370.014594.8199.0−54.13Elliptical orbit
ST20Circular orbit
SO110,378.1370.0545.0529.93132.9−17.74General case
SO211,378.1370.013094.8199.0−54.13Symmetric case
SO311,378.13704594.8199.0−24.13Same circular case
Table 2. Initial deviations of orbital elements.
Table 2. Initial deviations of orbital elements.
Spacecrafta/kmei/degΩ/degω/degn/deg
SO133.31.3 × 10−38.5 × 10−51.2 × 10−37.7 × 10−30.039
ST118.62.7 × 10−41.1 × 10−41.6 × 10−40.0150.129
Table 3. Simulation parameters used in the UKF.
Table 3. Simulation parameters used in the UKF.
NameValue
Initial state deviation of each spacecraft (km, km/s)[10, 10, 10, 1 × 10−3, 1 × 10−3, 1 × 10−3]T
Initial covariance of each spacecraft (km2, km2/s2)diag ([100, 100, 100, 1 × 10−6, 1 × 10−6, 1 × 10−6])
Standard deviation of measurement noise0.01 deg (approximately 40 arcsec)
Process noisediag (10−12)
Simulation step60 s
Table 4. Final STD and CR results for the two-spacecraft OD system in general cases.
Table 4. Final STD and CR results for the two-spacecraft OD system in general cases.
Indexx/kmy/kmz/kmvx/(km/s)vy/(km/s)vz/(km/s)
STD0.03060.02000.09933.074 × 10−52.9191 × 10−52.4185 × 10−5
CR99.69%99.80%99.01%96.93%97.08%97.85%
Table 5. The reduced row echelon in symmetric case.
Table 5. The reduced row echelon in symmetric case.
RowaOeOiOωOΩOnOaTeTiTωTΩTnT
11 −1
2 1 −1
3 1 1
4 1 −1−1.3
5 1 0.92
6 1 −1
Table 6. The reduced row echelon in same circular orbit case.
Table 6. The reduced row echelon in same circular orbit case.
RowaOeOiOωOΩOnOaTeTiTωTΩTnT
11 1
2 1
3 1 −0.866 −0.354
4 1 −1−6.5 × 10−4 −0.511.319−1
5 1 0.7071 −0.866
6 1
Table 7. Nominal orbital elements for observability analysis.
Table 7. Nominal orbital elements for observability analysis.
Spacecrafta/kmei/degΩ/degω/degn/deg
SO411,878.1370.02300010
Table 8. Final STD and CR results for the cooperative OD with an additional observer.
Table 8. Final STD and CR results for the cooperative OD with an additional observer.
ConfigurationIndexx/kmy/kmz/kmvx/(km/s)vy/(km/s)vz/(km/s)
General orbitSTD0.00910.00100.01171.1744 × 10−66.0360 × 10−69.2720 × 10−7
CR99.99%99.99%99.99%99.88%99.40%99.99%
Symmetric orbitSTD0.00180.00760.1135.2919 × 10−65.0495 × 10−67.3012 × 10−6
CR99.99%99.99%99.99%99.47%99.50%99.27%
Same circular orbitSTD0.00360.00600.01055.4042 × 10−61.3093 × 10−64.3750 × 10−6
CR99.99%99.99%99.89%99.46%99.87%99.56%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, Y.; Qin, T.; Zhou, X. Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination. Aerospace 2022, 9, 166. https://doi.org/10.3390/aerospace9030166

AMA Style

Luo Y, Qin T, Zhou X. Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination. Aerospace. 2022; 9(3):166. https://doi.org/10.3390/aerospace9030166

Chicago/Turabian Style

Luo, Yan, Tong Qin, and Xingyu Zhou. 2022. "Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination" Aerospace 9, no. 3: 166. https://doi.org/10.3390/aerospace9030166

APA Style

Luo, Y., Qin, T., & Zhou, X. (2022). Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination. Aerospace, 9(3), 166. https://doi.org/10.3390/aerospace9030166

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