1. Introduction
The strapdown inertial navigation system (SINS) has emerged as the most attractive choice for autonomous navigations, particularly for military applications due to its immunity to external interference. However, the inertial instruments errors and initial misalignment errors of SINS lead to accumulated errors over time [
1]. Global Positioning System (GPS) is the most popular system used for the integration with SINS to enhance the performance and reliability of the navigation system [
2]. In practice, the anti-jamming problems and degraded accuracy in hostile environments make the SINS/GPS integration system a poor choice for autonomous and precise navigation. In comparison to GPS, the celestial navigation system (CNS) is accurate, autonomous, reliable, inexpensive and practically independent of all external inputs, but requires low dynamic or no rotation motion [
3]. As a result, the SINS/CNS integrated navigation can take advantage of these methods and achieve high precision in near-Earth flight vehicles [
4,
5].
Advancement in Optoelectronics and image processing techniques has enabled the development of Charged Coupled Device (CCD) star sensor which can provide the arc-seconds level accuracy of attitude in the inertial frame [
6]. The accuracy of a SINS/CNS integration system degrades to the level of arc-minutes without calibration of the installation error angles between SINS and the star sensor [
7]. Besides, the installation error calibrated on the ground would undergo changes after the launch of the flight vehicle and large acceleration [
8]. The star sensor measurement error induced by its installation error is considerably larger than its own error in the SINS/CNS integration system. The installation error angles, as one of the main factors, degrade the measurement accuracy of the star sensor in the practical application [
9]. At the same time, the observability analysis of the installation error is more difficult due to its definition in a body frame, rather than a navigation frame.
Currently, there is not enough literature data to design estimating methods of star sensor installation error, and its observability is rarely compared to other related inertial device errors in SINS/CNS integration for a near-Earth flight vehicle. Gerald and Shuster put forward an algorithm that can make full use of Satellite redundant attitude information to calibrate the star-sensor installation error on an orbit for a satellite. The method applied to satellites is not completely available in an SINS/CNS system for near-Earth flight vehicles [
10]. As shown in [
11], the misalignments of the star tracker are input as the system states, and the federated filter algorithm is used to estimate the state errors in the SINS/CNS/GNSS (Global Navigation Satellite System) integrated navigation system. An integration approach of SINS/CNS by Uscend Kalman Filter (UKF), Extend Kalman Filter (EKF) and neural network respectively was proposed by Fang [
12,
13,
14], and these methods require large computation efforts when applied in practical cases. Khan Badshah gave a method for fusing the data of SINS and CNS with a modified measurement model based on Linearized Kalman Filter (LKF) without considering the star sensor installation error [
15]. Considering the complicated circumstances of near-Earth flight, in order to achieve optimal attitude accuracy, it is necessary to estimate the star sensor installation error, and put forward its observable condition in the case of SINS/CNS. The concept of local observability and locally observable assertion is introduced in [
16,
17]. For near-Earth flight vehicles, star sensors could be frequently allowed to use shorter measurement intervals [
18]. Moreover, a locally observable system is sure to be globally observable, while a globally observable system is not sure to be locally observable [
2]. So the observability of this paper stands for local observability. The rank of observability matrix is also used to determine local observability [
19]. Consequently, an improved observability analysis of the installation error based on the eigenvalues of the observability matrix for a near-Earth flight vehicle is presented in this study. The simplest measurement schemes of the star sensor for a near-Earth flight vehicle are designed in order to estimate on-line and compensate for the installation error between SINS and a star sensor. Simulations show that the estimated value of errors are in good agreement with the true value, and the platform error angles are less than 7 arc-seconds over the measurement schemes which include rolling and heading adjustment for LFOV. In this work, an integration algorithm for a conventional star sensor is proposed in a local level mechanization. A linearized error model of SINS is used for analysis of the stand-alone SINS and integrated system. KF is used for errors estimation and analysis of the integrated system. The measurement model of KF is developed with the installation error between SINS and the star sensor.
This paper is categorized into five different sections. In
Section 2, the integration algorithm is deduced, and Kalman Filtering for the integrated system is designed. The local observability analysis method is described and an observable condition is proposed in detail in
Section 3. Simulation results are produced and discussed in
Section 4. Finally, the paper is concluded in
Section 5.
3. Local Observability Analysis of the Installation Errors
Formally, a system is said to be locally observable, for any possible sequence of state vectors, if the current state can be determined in finite time only by the measurements. It can be evaluated by forming the system observability matrix and checking the rank of the matrix [
16]. A full rank observability matrix means that all the states can be estimated when the noise characteristic of the system is known. It is noticeable that the observability of the star sensor installation errors is only related to the platform's error angles. Therefore, we can choose the installation errors and platform's error angles as the states. The star sensor and the near-Earth flight vehicle undergo strap down installation. The vehicle cannot perform pose adjustment when star sensor conducts measurement for the accuracy of star sensor. Hence, we can suppose that platform's error angles and installation errors remain constant in the ultra-short time period for simplicity, i.e.,:
It also should be noted that local observability is related to the number of stars or the field of view. Accordingly, two categories are included in the analysis below. One is a large field of view (LFOV) star sensor which measures two or more star vectors and the other is a narrow field of view (NFOV) star sensor measuring only one star vector.
3.1. Local Observability of LOFV Star Sensor
According to the principle of SINS/CNS, the following equations can be obtained in the three different-time measurements of two star vectors guaranteeing that the rank of observability matrix is equal to 6. It is determined by the process of measuring a star. Proof is given below with properties of a matrix rank and orthogonal matrix.
The observability matrix of an equation can be expressed as follows,
where
The rank of the observability matrix meets the following inequality Equation (12).
Apparently, the rank of
A is 9, and the rank of
B is less than or equal to 6. Therefore, the rank of the observation matrix is only determined by the rank of
B.
Row elementary operations on
B matrix yield:
Based on the reduced order theorem of the matrix rank, the rank of
B can be given below
We rewrite the above formula with the theorem of the block matrix rank, then we have:
If T is an attitude matrix and if it is orthogonal, then the difference can be expressed as:
Note that the attitude matrix is the direction cosine matrix and its rank is equal to 3. Hence, the rank of the difference of attitude matrices is presented as follows
The matrix
of the above formula is an order 3 real orthogonal matrix the determinant of which is equal to 1, and the number of eigenvalues is 1 or 3. The rank of
can be calculated after diagonalization of the matrix
,
where
eig is the function of matrix eigenvalue.
and
are mutual complex conjugates . The results are summarized by the following conclusion in the observation of two star vectors.
Case LOFV-1: The rank of the observation matrix is 3, when the attitude matrix remains constant in three different time measurements. Accordingly, not all of the installation errors can be observable. The observability is mainly constrained by the rank of the
B matrix. The lack of posture adjustment during three different time measurements determines that the rank of the B matrix is equal to 3.
Case LOFV-2: One of the installation error angles is not observable, in which case just one of the attitude angles changes during three different time measurements. The one which is unable to be estimated is predictable. Obviously, the direction of the unobservable installation error angle is parallel to the rotation of the virtual axis of the star sensor.
is defined as the transition matrix that rotates at the angle of around the r axis.
Case LOFV-3: Two of the attitude angles change during the three different time measurements, all of the installation errors are clearly observable.
This conclusion gives a sufficient and necessary condition to confirm the installation errors’ observability for the LFOV star sensor in a SINS/CNS integrated system. The above results are suitable for more than 2 star vectors. The increasing number of star vectors does not change observability, but can improve the degree of observability and the capacity of the resisting disturbance.
3.2. The Observability of the NFOV Star Sensor
The observability matrix of Equation (11) can be obtained for the three different time measurements of one star vector, and it can be expressed as follows,
where
The rank of observability matrix satisfies the following inequalities.
Apparently, the rank of
A is 6, and the rank of
B is less than or equal to 6. Therefore, the rank of the observation matrix is still determined by the rank of
B.
The results are summarized by the following conclusion.
Case NOFV-1: When the attitude matrix remains constant during three different time measurements, none of the installation errors are observable.
Case NOFV-2: At least one of the installation error angles are not observable, in which case one of the attitude angles changes during the three different time measurements.
Case NOFV-3: Two of the attitude angles change during three different time measurements, all of the installation errors between SINS and CNS may be clearly observable.
The above-mentioned conclusion gives a necessary condition observable for the NFOV star sensor installation error in a SINS/CNS integrated system. The installation error around the optical axis is not observable in any case for NFOV.
4. Simulation and Discussion
In this section, a simulation is carried out to check the observability analysis of the installation errors. In order to make the simulation result more credible, most of the error terms are included in a simulation error source. The simulating condition is set as follows:
A Fiber Optic Gyroscope (FOG) is chosen as a test gyro, the constant and random drifts of which are 0.03°/h, 0.006°/h, respectively. The scale factor error is 20 ppm, and the installation error is 5 arc-seconds. The constant and random biases of accelerometers are 50 g and 15 g. The scale factor error is 50 ppm, and the installation error is 5 arc-seconds. The SINS measurements are generated with a sample rate of 100 Hz. The precision of the star sensor is 5 arc-seconds. NFOV is about 2° and LFOV is about 20°. The error of the star catalog is 1arc-seconds. The misalignment error between INS and star sensor is about 200 arc-seconds. The star sensor measurement update is provided at a frequency of 10 Hz. The initial position error is set as 100 m. The initial attitude error is 1 arc-minute. One measurement operation lasts 5 s.
According to typical maneuver of near-Earth flight vehicles, the movement track of a vehicle in the simulation is designed as
Figure 2, ground initial alignment for 260 s and accelerating to near-Earth space. The vehicle then moved with a horizontal velocity that changed with small discontinuous accelerating and decelerating vertical velocity and pose adjustment. We measured the star vector three times until the end of the experiment. Simulations were performed in three cases based on the Monte Carlo method. Case 1 represents the situation that the attitude matrix remains constant during three different time measurements. Case 2 represents the situation that the roll attitude angles changes after the first measurements and the two remaining measurements are performed with the same attitude. Case 3 represents the situation that the roll attitude angles changes after the first measurements and the heading attitude angles changes after the second measurements. The estimated values of the installation errors are illustrated in
Table 1,
Table 2 and
Table 3 respectively, which indicate the observability in the SINS/CNS.
In Case 1, the near-Earth flight vehicle just finished the first starlight measurement. Three installation errors of the star sensor would not be observed. The observability is unable to be improved whether for LFOV or NFOV. The simulation results shown in
Table 1 are in agreement with the analysis and observability is mainly influenced by posture adjustment times instead of the increasing number of stars.
In Case 2, the near-Earth vehicle had two different observation attitudes, the heading of which changed. Compared with Case 1, the observability is improved greatly. Two installation error angles of the star sensor can be observed for the LFOV star sensor. The unobservable error
is just parallel to heading axis. So the optimal observability can be based on more than two star vectors. However, only one of the installation error angles can be estimated due to
being parallel to the optical axis for NFOV. Thus, the simulation results shown in
Table 2 are in agreement with the analysis and it is still mainly influenced by the viewing posture for LFOV.
In Case 3, the near-Earth vehicle had three different observation attitudes, the heading and rolling of which are changed. The results are shown in
Table 3. Compared with Case 1 and Case 2, the observability is improved. All the installation error angles can be observed for the LFOV star sensor due to the times of posture adjustment being equal to 2. It must be noted, however, that the installation error is still not fully observable owing to the unobservable installation error being parallel to the optical axis in the case of NFOV. Consequently, a simple position adjustment around the optical axis can make the observability optimal for NFOV. With more than one star vector, the installation errors are completely observable. In summary, the number of star vectors should be greater than or equal to 2, and the times of posture should be greater than or equal to 2.
Simulation results verify the correctness of observability analysis above. The installation error convergence curve is given in
Figure 3 measuring the starlight three times for LFOV. At the end of the first rolling posture adjustment, the estimated value of installation error
parallel to the rolling axis is unchanged, the same as
parallel to the heading axis at the end of second heading posture adjustment. Other installation error angles basically converged to the true value. At the end of the second heading posture adjustment, the estimated installation error
is close to its true value, and still a little weaker than other installation error angles parallel to the optical axis. In all cases, the installation error is observable by means of rolling adjustment and heading adjustment for LFOV.
The previous conclusion showed a remarkable efficiency in the estimation of the installation errors. In order to verify that the accuracy of SINS/CNS can be clearly increased, it is necessary to undertake simulations under three conditions, namely: SINS, SINS/CNS integration method without considering star sensor installation errors (SINS/CNS-1), SINS/CNS integration method considering star sensor installation errors (SINS/CNS-2).
The simulation results are only intended to emphasize the fact that it is essential to precisely estimate the installation error for the SINS/CNS integration system. In this specific case, it is not only a filtering algorithm issue, and previous cited work could help only in part. The real need is for a better modeling, i.e., for the proposed algorithm that includes installation errors. As can be seen in
Figure 4d–f, attitude errors and azimuth error are completely estimated along with vehicle maneuvering in the case of SINS/CNS-2, increasing the observability of star sensor installation errors in order to improve the accuracy of the CNS/SINS integrated system. According the principle of SINS/CNS integration, observability of the velocity and position errors is poor. From the simulation results in
Figure 4a–c, position accuracy was also improved significantly with increasing attitude accuracy. The main reason for this is that the velocity as well as the position errors due to platform misalignment are compensated for. However, the slight increase in longitude error owing to the accelerometer cannot be compensated for after the third calibration measurement in the case of SINS/CNS-2. For the case of SINS/CNS-1, the integration method is not convergent as a result of inaccurate measurement, and navigation accuracy is inferior to the case of SINS. It is essential to precisely estimate the installation errors for the SINS/CNS integration system.
The detailed navigation error results of simulations are illustrated in
Table 4. As in the case of SINS/CNS-1, the heading error is close to the installation errors of the star sensor, which may restrict the actual accuracy of SINS/CNS. The position error is greater than in the case of SINS due to inaccurate estimation of inertial device errors caused by the star sensor installation errors. Compared with SINS, the SINS/CNS-2 integration method presented in this work can achieve the same precision as the star sensor. The attitude errors are accurately compensated for and reduced from 559.48 arc-seconds to 2.18 arc-seconds. The position errors are limited to less than 103 m in about a 220-s flight. Consequently, the result is indicative of effective estimation about the star sensor installation errors. Meanwhile, the CNS/SINS integration method and observability conclusion are validated.