Maximum determinants of FIM related to different number of beacons.
In order to achieve more scientific returns for Mars, future Mars landers will be required to land at certain landing point with special scientific interest. Therefore, autonomous navigation is indispensable during the Mars approach, entry, and landing phase. However, the number of beacons or the Mars orbiters which can provide the navigation service is so limited and the line-of-sight visibility cannot be guaranteed during the landing period. So the navigation scheme especially the beacon configuration has to be optimized in order to efficiently use the limited navigation information. This chapter aims to analyze the feasibility and optimize the performance of the Mars Networks-based navigation scheme for the Mars pinpoint landing. The observability of navigation system is used as an index describing the navigation capability. Focusing on the relationship between the configuration of radio beacons and observability, the Fisher information matrix is introduced to analytically derive the degree of observability, which gives valuable conclusions for navigation system design. In order to improve the navigation performance, the navigation scheme is optimized by beacon configuration optimization, which gives the best locations of beacons (or the best orbit of navigation orbiters). This is the main approach to improve the navigation capability.
- Mars networks
As the most similar planet to the Earth in the Solar system, Mars is considered as an ideal target for planetary exploration [1, 2]. Since the 1960s, humans have investigated the Mars exploration missions in the near distance. With the development of aerospace science and technology, the manner of Mars exploration has shifted from flyby/orbiting to landing and roving explorations. Considering scientific returns and exploration capabilities, Mars landing exploration is also essential and is one of the most popular tasks of human deep space exploration in the near future. The representative Mars landing missions including NASA’s Viking 1 and 2, Mars Pathfinder (MPF), Mars Exploration Rovers (MER, including the Spirit and Opportunity rovers), Phoenix, Mars Science Laboratory (MSL, including the Curiosity rover), and ESA’s Mars Express/Beagle 2 mission. All of these greatly inspire the development of advanced guidance, navigation, and control (GNC) technologies.
During the past 50 years of Mars exploration, 46 Mars exploration spacecraft have been launched. The overall success rate is only 41.3% though. Furthermore, among the 20 Mars landing attempts, only 7 robotic rovers were successful. The success rate for Mars landing missions is only 35%. Among the failed landing missions, most failures occur during the landing phase. The pinpoint landing has to be based on the precise autonomous navigation technology.
In the entry phase of a Mars landing, the lander is covered by a heat shield which blocks the optical sensor measurement, causing that all landers relied on the Inertial Measurement Unit (IMU) recursion. The initial errors of the lander cannot be corrected by IMU data. Even worse, the recursion errors using IMU are accumulated due to the sensor bias and noise. To overcome the incapability of IMU, the Mars Network-based Mars entry navigation is developed based on high frequency radio communication between the lander with ground or orbiting radio beacons [3, 4, 5]. Involving the radio measurement date into a navigation filter, the position and velocity of the lander can be optimally estimated.
The Mars Network-based Mars entry navigation is faced with two challenging. One is that the geometric configuration of the radio beacons affects the navigation performance. The other is that the available beacons at present are very limited. Considering these two factors, effort should be devoted to optimizing the configuration of radio beacons to maximize the function of the limited beacons. In , the navigation accuracy from the Extend Kalman Filter (EKF) by processing the radio measurements is analyzed, and the optimal configuration of ground beacons is selected among potential beacon position. Yu focused on the navigation observability and take it as a performance index to optimize the configuration of radio beacons . The research on ground beacons, to some extent, inspired the future Mars landing navigation. However, the practice application of ground beacon-based navigation is hardly applied in practice. The first concern is that no ground beacon is available. Even if several beacons are distributed on Mars surface, it’s still a tough job to place them exactly at the optimal locations. Moreover, the accurate positions of the beacons are hardly obtained accurately. Considering the immovability of ground beacons, the potential location areas are constrained by the line-of-sight visibility, resulting in an unsatisfactory beacon configuration during the entry phase.
As a substitution of ground radio beacons, the Mars orbiters which can also serve as beacons for Mars Network-Based Navigation are of more practice value. Currently, the operational orbiter around Mars includes 2001 Mars Odyssey and 2005 Mars Reconnaissance Orbiter. With another forthcoming spacecraft Mars Atmosphere and Volatile Evolution (MAVEN) , the capability of Mars network can be further increased. Focusing on how to fulfill the function of a Mars network, Ely firstly established the basic principle to design a constellation for navigation . Then, taking the Mean of the Position Accuracy Response Time (MPART) as the performance index, the constellation configuration was optimized . In , the number of orbiters and the coverage was considered to design the Martian navigation constellations envisaged in the ESA’s Martian Constellation for Precise Object Location program. The optimization method of the above researches is inherited from the Global Positioning System (GPS). The global navigation performance was emphasized. For the limited amount of Mars orbiters, global coverage is difficult to realized, and local navigation performance should be investigated thoroughly for specific missions. Moreover, the effect of geometric configuration of the Mars network on the navigation performance should be revealed clearly. Inspired by these requirements, Yu et al. optimized the orbits of Mars orbiters in the observability point of view, and tried to explain the relationship between the configuration of beacons and orbiters and the navigation capability [13, 14].
To optimize the configuration of the radio beacons, a performance index should be firstly setup. The observability of the navigation system is selected as the performance index since it reflects the navigation capability directly. A lot of work has investigated the observability of linear and nonlinear dynamic systems [6, 15, 16, 17]. However, the analytic relationship between geometric configuration and observability has never been revealed. According to Cramér-Rao inequality , the inverse of the Fisher Information Matrix (FIM) estimates the lower bound of the estimation error. Therefore, FIM can be used to quantify the observability of the navigation system [19, 20, 21]. In this circumstance, some valuable analytic conclusions about the navigation design can thus be obtained.
Based on the requirement of the navigation optimization for Mars pinpoint landing, this chapter discusses the design and optimization of the Mars Networks-based navigation during Mars entry phase. Firstly, the Mars Networks-based navigation scheme is introduced, and the dynamic model and the observation model are given. Based on the navigation system, the observability of the Mars entry navigation analysis, and the analysis methods based on the quadratic approximation and Fisher information matrix are proposed. The relationship between the observability and the beacon configuration is derived, and the theoretically optimal configuration is given. Considering the constraints of Mars entry scenario, the ground beacons and the orbit of Mars orbiters are optimized based on observability based on an entry trajectory. The simulations also indicate the improved navigation performance.
2. Mars networks-based navigation scheme
2.1. Dynamic model of Mars entry phase
In the dynamical model with respect to a stationary atmosphere of a rotating planet, the 6 dimensional states x of the entry vehicle include
In the equation, is the banking angle, which is fixed at 0 in the following analysis. refers to the rotation rate of Mars. For simplicity, the second order terms of are neglected, which is feasible because the value of is quite small. Then the gravity acceleration g, lift and drag accelerations l and d are given by
where is the Martian gravitational constant. S and m denotes the reference area and mass of the entry vehicle, and and are the lift and drag coefficients respectively. Furthermore, the Mars atmospheric density is approximated by the conventional exponential model
where kg/m3 is the reference density, km is the reference radial position, and m refers to the atmospheric scale height. The dynamical model of the entry vehicle is abbreviated as .
2.2. Observation model
The radio ranging and velocity data between the lander and the radio beacon can be measured through radio communication, given by
where is the real range between the lander and the
The relative velocity model is given by
where is the real line-of-sight relative velocity between the lander and the
With different radio beacons come different navigation scenarios. Without losing the generality, the observation model can be summarized as . Obviously, both radio measurements in Eqs. (6) and (7) are nonlinear. Moreover, the navigation performance is closely related to the geometric configuration of radio beacons. Therefore, the beacon configuration needs to be optimized based on the observability analysis.
3. Observability of the navigation system
3.1. Observability analysis based on the quadratic approximation
Consider the following nonlinear system:
where is the n-dimensional state vector and is the m-dimensional observation vector. Define as the nonlinear measurement operator.
The Lie algebra is an efficient tool for observability analysis. For the
The differential of is defined as
Regarding the zero-order Lie derivative of the jth measurement function as itself, the matrixis given as
It is proven that the dynamical system at state is locally observable if the observability matrix given below has the rank of n.
It’s a heavy burden to calculate the observability matrix in Eq. (12) due to the existence of high order differential, especially for the 6-dimensional dynamics of Mars entry phase which requires the calculation of 5th order Lie derivatives. Next, a quadratic approximation method is developed to simplify the computation of the observability matrix.
First of all, the quadratic approximation of the kth order Lie derivative is given as
where is the value of at , and and refer to, respectively, the Jacobian and Hessian matrix of at . The linearized state equation is given by
in which refers to the value of at , and is the Jacobi matrix of at .
This leads to
The observability matrix can be computed as
Obtaining and , the observability matrix can be iteratively calculated. Only 2nd order differential of is needed here to compute the Jacobian and Hessian matrices, reducing largely the computation cost.
Linearize the dynamical and observation model by first-order approximation
Construct the observability matrix according to the linear system theory
The Hessian matrix is involved in the quadratic approximation, improving the accuracy of observability analysis compared with the linearized observability analysis. However, the higher order terms of may appear when computing in Eq.(9) if the state equation is approximated to a higher order. In this case, the predetermined presentation form in Eq.(13) is no longer valid. One way to defeat this case is to increase the approximation order of Lie derivatives. Note that tensor calculus can be involved and the computation complexity is increased. Thus, the trade between accuracy and computation cost is balanced by the quadratic approximation of Lie derivatives and the linearization of state equation.
In the optimization of observability, the condition number of observability matrix is selected as the performance index, given by
where and are, respectively, the maximum and minimum singular value of the matrix. The condition number measures the singularity of the matrix. A larger condition number means a more singular matrix. Here we take the inverse of the condition number to quantify the system observability.
Obviously, the observability degree is in the interval [0, 1]. When , the observability matrix is rank defect, and the navigation system is locally unobservable. When , the observability is full rank, indicating an observable navigation system.
3.2. Observability analysis based on the fisher information matrix
Without loss of generality, we will consider the nonlinear observation models
This equation may describe the measurement of relative range and range-rate according to Eq. (6) and (7). Meanwhile, in order to investigate the impact of different measurement methods on the observability of position and velocity of the entry vehicle separately, the 3-dimensional state x may be r or v of the entry vehicle. The likelihood function of x is defined as the joint probability density function of multiple measurements given by
Then, take the negative of the natural log of Eq. (23) and omitting the terms not related to x, and the loss function can be derived as
Find a state vector to minimize and the state vector is the optimal estimation of the lander’s states. The FIM of the state is given by
The estimate error covariance and FIM satisfy the following equation
where is the estimate error covariance, and “” means that () is positive semidefinite. According to Eq. (26), the FIM can be used to evaluate the lower bound of the estimation error covariance, and further the system observability. Give the trace of in Eq. (27).
where are the eigenvalues of F. It’s illustrated from Eq. (27) that larger eigenvalues of the FIM leads to smaller trace of estimation error covariance and stronger system observability. Quantify the observability by the determinant of FIM . The following relationship can be obtained.
Eq. (28) means that the trace of FIM measures the lower bound of estimation errors.
4. Observability analysis of Mars networks-based navigation
4.1. Observability analysis using only range measurements
In this subsection, the system observability using only range measurements between the lander and ground beacons is analyzed. Since no velocity information is included in Eq. (6), only the observability of the position vector is studied. The cases with different amount of beacons are studied.
4.2. One-beacon case
In this case, the FIM is given by
The rank of the matrix is only one. Solving the following equation
Clearly, the eigenvalues of are given by twice repeated 0 and . Therefore, the eigenvalues of are given by .
Next, we have the eigenvector corresponding to
The vector corresponds to the observable state combination, and means that only the state component along the vector can be observable.
According to Eq. (28), the lower bound of estimation errors can be obtained as
Eq. (32) means the lower bound of estimation errors is higher than the estimation accuracy. In another word, the estimation accuracy cannot be higher than the measurement accuracy. Note that, even if multiple beacons are involved in the navigation system, the observability is still deteriorated if the beacons are located in similar direction.
4.3. Two-beacon case
Assume two non-collinear beacons, the FIM in Eq. (25) is derived by
Involving one more measurement, the rank of is increased to two. The observable state combinations can be obtained by solving the eigenvalue and eigenvector. In this case, the eigenvalues of are given by , . The eigenvector corresponding to the zero eigenvalue is obtained as
The vector gives the unobservable state component which is in the direction perpendicular to the plane constructed by and . From an opposite view, all state components in plane are observable.
Since, in this case, the observability matrix is still zero, the navigation system is unobservable. According to Eq.(28), the lower bound of the estimation errors can be obtained as
4.4. More-than-two-beacon case
In this case, the FIM is given by
The matrix has a full rank, indicating an observable system. The determinant of is given in Eq. (37).
The detailed derivation can be found in Ref. . From Eq.(7), we can know that more radio beacons, no matter where they are, increase the determinant of the FIM, thus increase the system observability. To analyze the maximum value of , Eq. (37) is reorganized as
where is the minimum value among . The selection of the direction of radio beacons to maximize the observability can be described by the following optimization problem
Note that the locations of radio beacons are not constrained. In cases with three beacons, the determinant of is maximized if and only if
|Number of beacons||Maximum determinant of FIM|
According to the results in Table 1, the relationship between the maximum determinant and the number of beacons can be induced by an exponential formulation, given by
The lower bound of estimation errors is derived as
The change of lower bound of estimation errors with number of beacons is shown in Figure 1. It’s shown that with more beacons comes more accurate estimation. However, the increasing rate of accuracy is slowed down, indicating that the navigation accuracy cannot be improved endlessly by only increasing the number of beacons.
4.5. Observability analysis of the navigation using range-rate measurements
4.5.1. Observability analysis of vehicle’s velocity
The FIM of vehicle’s velocity using range-rate data is given by
Eq. (42) has a similar form with Eq. (36) which describes the FIM of position. The only difference lies in the measurement deviation. Hence the same conclusion of the observability of velocity can be obtained as that in Section 4.1. The detailed analysis is omitted here.
4.6. Observability analysis of vehicle’s position
Using the range-rate measurements, the FIM of the lander’s position is derived as
where is given by
The FIM here is much more complicated than that in Section 4.1 due to the involvement of both range and velocity information in FIM. Define , the following equation can be obtained.
When there are one or two beacons, the FIM is rank defect, and the navigation system is also unobservable. With three or more beacons comes the full-rank FIM. In this section, only the observable cases are focused on.
It is also concluded that the determinant of FIM will be zero if only one or two beacons is used, which indicates that the position of entry vehicle will be observable if more than two beacons are used. Furthermore, we focus on three-beacon and more-than-three beacon cases. Compare Eq. (46) with Eq. (36), we can find that the determinant of FIM for range-rate measurement cases has a similar format as Eq. (37)
It’s shown that Eq.(46) has the similar format with Eq. (37). Thus, the change of the observability with the number of radio beacons is similar with the results in Table 1. However, due to involving relative range and velocity information, the optimal geometric configuration is different with the cases using only range measurements.
The lower bound of estimation errors in this case is evaluated by
where is the lander’s velocity value. Obviously, more radio beacons lead to more accurate estimation. Since the value of relative range is much bigger than relative velocity, the lower bound of estimation errors using range-rate date is larger than that using range data. Besides, it’s concluded that more accurate range-rate measurement, closer relative range, and slower velocity can realizer more accurate position estimation.
5. Orbit optimization based on observability analysis
5.1. Optimization of navigation using ground beacons
The configuration radio beacons is expressed by the following set
where is the position of the
To realize the Mars network-based navigation, the visibility of the beacons to the lander should be guaranteed. Define two unit vectors as follows
where is the relative position vector from the lander to the radio beacon, obtained as
To guarantee the visibility, the two vectors in Eq. (50) should satisfy
The schematic of visibility is shown in Figure 2.
The optimization problem of beacon configuration is given as
where is the set of the areas of radio beacons that satisfy the visibility during the whole entry phase. In this optimization problem, the global optimization algorithm is selected to obtain the optimal beacon configuration.
The initial states of the lander are listed in Table 2.
It is assumed that the Mars entry phase lasts approximately 240 seconds. The entry trajectory and the corresponding visible area are shown in Figure 3.
Three radio ranging measurements at a certain time can geometrically determine the position of the lander. Thus the navigation scenario with three beacons is first analyzed with respect to the observability. The optimal locations of beacons are displayed in Figure 4.
The optimal three beacons are located close to the edge of both sides of the visible area. The beacon on the east side is almost along the entry trajectory, while the west two beacons are separated on the north and south side of the entry trajectory. The observability degree in this situation calculated by different methods is illustrated in Figure 5, and the computation time for each method is listed in Table 3.
|Analysis approach||Computation time, s|
|Method with Lie algebra||>10,000|
|Method based on quadratic approximation||2.1558|
Figure 5 shows a huge undulation in observability degree during the Marts entry phase. The maximum and minimum value are and respectively. Considering the machine precision, the navigation system is observable only if the observability exceeds . The observability degree during the entire entry phase passes through the threshold, and thus, the navigation system is observable. The minimum degree of observability occurs at the beginning of the entry phase when the entry vehicle is at its greatest distance from radio beacons, while the maximum degree of observability occurs when the entry vehicle approaches two beacons on the west side. In order to explain the evolution of the degree of observability. An index related to the geometric configuration of the lander and radio beacons is given in Eq. (54) to explain the evolution of the degree of observability.
where , , and are the unit vectors from the beacon to the lander,
The observability degree obtained from the three methods is quite close to each other. However, the method based on Lie algebra consumes the most time. The linearization method provides the largest deviations, especially at the peak time, indicating a relatively low accuracy. The proposed quadratic approximation method achieves a performance balance in accuracy and complexity. To analyze the navigation accuracy, the Extended Kalman Filter (EKF) is used to estimate the lander’s states. The range measurement error is assumed to be Gaussian white noise with a standard deviation of 100 m. The initial errors are listed in Table 4. The estimation errors and the 1-sigma uncertainty bounds are depicted in Figure 7.
It’s illustrated that
5.2. Optimization of navigation using Mars orbiters
Compared with ground beacons, the Mars obiters are constrained by the orbital dynamics, which is considered to be two-body dynamics here. In this subsection, the initial states of the Mars orbiters are considered as the optimized variables. Furthermore, assuming that the Mars orbiters moves in a circular orbit, the variables to be optimized are simplified as inclination i, longitude of ascending node Ω, and the true anomaly
Given the initial states of the Mars obiter, the subsequent states can be obtained by propagating the two-body dynamics. Likewise, the trajectory of the lander can be also obtained by propagating the entry dynamics. To evaluate the overall performance of the observability of the entry phase, the integration of the observability is taken as the performance index, given by
The angle between the position vectors
The visibility requires that
The schematic of the visibility is illustrated in Figure 8. The gray part represents the area in which the Mars orbiter is invisible to the lander.
Then the orbit optimization problem is given by
In the optimization problem, the performance index cannot be expressed explicitly by the optimization variables, and the gradient cannot be obtained. Thus, the heuristic global optimization algorithm is chosen to solve the optimization problem. The lander’s initial states are listed in Table 5 with the assumption of a ballistic entry having a banking angle of zero. The duration of entry phase is setup as 240 seconds.
The navigation scenario with three Mars orbiters is analyzed. The nominal orbit altitude of the three orbiters is 725 km. The observability is quantified by
At a certain epoch, the maximum value of
|Orbit element||Obiter 1||Orbiter 2||Orbiter 3|
It’s shown that the three orbiters keep a relatively stable configuration, and stays orthogonal approximately to each other. The value of maximized performance index is 237.963. The observability almost reaches the maximum value all the time during the Mars entry phase. The comparison of Mars obiters-based navigation and ground beacon-based navigation is performed. The observability degree of these two scenarios is shown in Figure 10.
The fixed ground beacons have limited locations due to the visibility constrain and the geometric configuration cannot remain optimal during the entry phase. Thus, the observability is undulated to a large extent. The Mars orbiters overcome this defect with its moving property. To show straightforward the geometric configuration., the observability degree is close to maximum value at each epoch during the Mars entry phase. The angles between the vectors
It’s shown that, using the ground beacons, the angles between the three vectors change dramatically in the entry phase. The optimal configuration can be met only at the epoch of 75 s. However, for the orbiter-based navigation scheme,
Next, 500-time Monte Carlo simulations of navigation systems based on EKF are carried out. The initial position and velocity have standard deviations of 1 km and 0.5 m/s respectively. The measurement error is set to be 50 m, and considered as Gaussian white noise. The simulation results are shown in Figure 12.
Since no information of entry vehicle’s velocity is provided from range measurements, the convergence of velocity estimation is not as quick as position estimation. A much better navigation performance can be achieved by the Mars orbiter-based navigation. It can be concluded that the configuration of orbiters is a main contributor to the navigation performance. The Mars orbiter-based navigation, which can achieve a better configuration, is more practical for Mars entry navigation.
This chapter introduced the Mars Networks-based navigation for the Mars entry phase. Based on the navigation scheme, the observability of the navigation system was analyzed using the proposed two novel observability analysis methods. Furthermore, the beacon configuration was optimized based on observability considering the line-of-sight constraints were concluded that the beacon configuration is a main contributor to the Mars Networks-based navigation. The observability analysis showed that an improved behavior of observability and more flexibility of beacon configuration determination can be achieved using more beacons. Navigation also demonstrated this conclusion. Meanwhile, compared with the ground beacons, Mars orbiters may be a better choice as Mars Network which gives a more accurate navigation result.