- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Abstract and Applied Analysis
Volume 2014 (2014), Article ID 494060, 11 pages
Finite-Horizon Robust Kalman Filter for Uncertain Attitude Estimation System with Star Sensor Measurement Delays
1College of Automation, Harbin Engineering University, Harbin 150001, China
2Dongguan Techtotop Microelectronics Co., Ltd., Chengdu R & D Branch, Chengdu 610041, China
Received 10 October 2013; Accepted 6 January 2014; Published 23 February 2014
Academic Editor: Xiaojie Su
Copyright © 2014 Hua-Ming Qian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper addresses the robust Kalman filtering problem for uncertain attitude estimation system with star sensor measurement delays. Combined with the misalignment errors and scale factor errors of gyros in the process model and the misalignment errors of star sensors in the measurement model, the uncertain attitude estimation model can be established, which indicates that uncertainties not only appear in the state and output matrices but also affect the statistic of the process noise. Meanwhile, the phenomenon of star sensor measurement delays is described by introducing Bernoulli random variables with different delay characteristics. The aim of the addressed attitude estimation problem is to design a filter such that, in the presence of model uncertainties and star sensors delays for the attitude estimation system, the optimized filter parameters can be obtained to minimize the upper bound on the estimation error covariance. Therefore, a finite-horizon robust Kalman filter is proposed to cope with this question. Compared with traditional attitude estimation algorithms, the designed robust filter takes into account the effects of star sensor measurement delays and model uncertainties. Simulation results illustrate the effectiveness of the developed robust filter.
Attitude estimation has played an important role in many actual applications, such as aerospace, satellites, marine, and robots. For attitude estimation system, due to the high measurement precision of star sensor, the rate gyro and star sensor are often integrated to determinate the spacecraft attitude. Furthermore, the filter design is one of the key technologies in attitude estimation. As is well known, Kalman filter has been employed to solve the attitude estimation filtering problem [1–3]. Although these attitude estimation filtering algorithms are available for handling attitude estimation problem, they need to know the accurate model with Gaussian noises and assume exact alignment of gyro and star sensor. However, in practical problems, the measurement misalignment errors of these sensors are inevitable, which will severely degrade the filtering performance. To overcome the sensor misalignment problem in attitude estimation, many researches have been reported in some recent notes [4–8]. For example, Shuster et al.  utilize the batch estimation technique to calibrate the misalignment of the sensors. Pittelkau [5, 6] develops the Kalman filtering technique to estimate the calibration parameters of gyro and star sensor. Lai and Crassidis  derived a new spacecraft sensor alignment estimation approach based on the unscented filter. Vandersteen  presents the real-time moving horizon estimation of a spacecraft’s attitude and sensor calibration parameters. Unfortunately, even though misalignment calibration is accomplished, the measurement misalignment error of gyro and star sensor cannot be removed completely, which lead to model uncertainty. Therefore, in the case that an exact uncertain model is established, the robust filtering technique can be used to deal with the filtering problem with model uncertainty. For the purpose, in the past few decades, many researchers’ attentions have been drawn to the robust filtering problem with model uncertainties, including the filter , new energy-to-peak FIR filter , fuzzy filter design , and robust Kalman filter [12, 13]. Among them, the robust Kalman filter design based on the minimum variance theory has been approved to be an effective methodology. Based on this, Wang et al.  proposed a regularized robust filter for attitude determination system to deal with the installation error of star trackers. In this work, the installation error of star trackers is expressed as model uncertainty in the measurement model, but the measurement misalignment error of gyros is not taken into account.
In this paper, all the misalignment errors and scale factor errors of gyros introduced into the process model and the misalignment errors of star sensors in the measurement model are described as model uncertainties, so that uncertain attitude estimation model is established. From the uncertain model, we can find that the attitude estimation filtering problem suffers from uncertainties in the state and output matrices and uncertainty in the process noise matrix. A typical way is to represent the model uncertainties as norm-bounded uncertainties. Recently, the finite-horizon robust Kalman filter design has been investigated to be available for handling the filtering problem with model uncertainties in the state, output, and noise matrices by getting an optimized upper bound on the estimation error covariance . Souto and Ishihara  extend this work by considering correlated noises with unknown mean and variance.
However, the above works have been based on this assumption that all the observations should be available at the time of estimation. In many situations, the sensor measurements are disturbed by complicated signal processing circuit, leading to the sensor measurement delays or measurement failures [17–20]. Therefore, the filtering problems with sensor measurement delays have stirred considerable research attention, such as [21–26]. In attitude estimation system, due to optics imaging, star recognition, and attitude determination, the attitude information output of star sensors has the random delay characteristic. Up to now, the attitude estimation filtering problem with model uncertainties in the state, output, and process noise matrices and star sensor delays has not been reported. So, there is great desire to present a robust Kalman filter for uncertain attitude estimation system with star sensor delays.
Based on the above discussion, a finite-horizon robust Kalman filter is proposed for uncertain attitude estimation system with star sensor delays. The star sensor measurement is assumed as one-step randomly delayed measurement with different delay characteristics. The main contributions of the paper are as follows. The uncertain attitude estimation model is established to take into consideration measurement errors of gyros and star sensors, which indicates that the norm-bounded uncertainties appear in the state, output, and process noise matrices. Combined with star sensor delays, a new finite-horizon robust Kalman filter design is derived for the uncertain attitude estimation system. The Hadamard product is employed to help the robust Kalman filter development. The presented robust filter is recursive, which is suitable for online applications.
This paper is organized as follows. In Section 2, the uncertainty attitude estimation model with star sensor delays is set up. In Section 3, a finite-horizon robust Kalman filter for uncertainty attitude estimation system with star sensor delays is developed. In Section 4, the simulation results and analysis are given. In Section 5, some conclusions are drawn.
2. Uncertainty Attitude Estimation Model with Star Sensor Delays
2.1. Gyro Error Model
As Figure 1 shows, , , and are three axes of the gyro, respectively. The gyro model with misalignment errors and scale factor errors is given as follows: where is the gyro measured output, is the actual gyro angular rate, is the gyro bias, and are independent Gaussian white-noise processes with zero means and covariance and covariance , and is an unknown matrix with misalignment errors and scale factor errors, which is defined by where is the unknown scale factor error vector and is the projection of the -gyro axis on the body-axis, which is assumed to be a small and unknown misalignment angle.
2.2. Uncertainty Process Model
The quaternion is employed to express the attitude for the attitude estimation system consisting of the gyro and star tracker. So, the quaternion orientation equation is described as where is the attitude quaternion, is the quaternion vector, is the quaternion scalar part, is the quaternion product, and can be defined as follows: where is a cross-product matrix defined by
Using quaternion multiplication, the quaternion error is expressed as where is the true quaternion, is the estimated quaternion, is the inverse quaternion, which is given by , and is the quaternion error vector part.
The gyro error angular rate is assumed as the difference between the estimated and actual angular rate: . According to (1), using a small angle approximation, we have where is the gyro bias error vector.
Differentiating (6) with respect to time and combining the quaternion multiplication, we obtain
In order to avoid the quaternion normalization constraint, only the vector component of the quaternion error is considered in the states. Inserting (7) into (8), by neglecting the second-order terms, we have
The quaternion error vector part and the gyro bias error vector are constructed as the error state vector: . An error state process model with unknown misalignment errors and scale factor errors can be expressed as
According to (10), the discrete-time process equation can be developed as where is the zero mean Gaussian noise with covariance ,
From (11), it can be seen that the unknown error matrix not only appears in the state matrix but also affects the statistic of the process noise. In order to realize the robust filtering design, the unknown error matrix can be rewritten as where ; ,
The parameters and are positive constants, which can be chosen by the priori information of the gyro installation errors. If the and are set to be large enough, the inequalities and can be fulfilled, so that the inequality is satisfied. According to (13), the model error matrices and can be described as where ,
2.3. Uncertainty Measurement Model with Star Sensor Delays
To obtain the attitude information, three star sensors are chosen. Considering the misalignment error of star sensors, the measurement model with model errors is expressed as where , , and denote different star sensors, are the measured quaternion error vector parts of star sensors, is the unknown misalignment error vector, and are the zero mean Gaussian white noises with covariance matrix . The unknown cross-product matrix can be written as where
The parameters are positive constants. If are large enough, the inequalities and can be satisfied. For convenience, the uncertain measurement with misalignment errors can be rewritten as where is the zero mean Gaussian white-noise process with covariance ,
Since the attitude information of star sensors is given by the complicated data processing and transmission, the output measurement is delayed. The delayed star sensor measurement is assumed as where is the true measurement output vector, accounts for the different delay rates, and are independent random variables taking the values of 1 or 0 with where is a known scalar.
Remark 1. The measurement errors of sensors are inevitable in real applications, and the gyro and star sensor are no exception. Though the literature  takes into consideration the misalignment error of star sensors, no attention is paid to the misalignment errors and scale factor errors of the gyro. However, they can lead to the uncertainty process model. As shown in (11), the uncertainty error matrix exists in the state matrix and noise matrix, which influences the design of robust filter. Besides, the sensor measurement signal transmission is susceptible to interference from the external environment and limited bandwidth of network, which makes the sensor measurement delay occur. In (22), the delayed star sensor measurement model is established. As discussed in the work , different delay rates are taken into account by introducing the diagonal matrix .
3. Finite-Horizon Robust Kalman Filter for Attitude Estimation
3.1. Problem Description
Considering the uncertain discrete-time linear stochastic system with sensor delays where is the state vector, is the measurement vector, is the true measurement output vector, and are uncorrelated process and measurement Gaussian noises with zero means and covariance and covariance , and , , and are known matrices with appropriate dimensions. The matrices , , and represent uncertainties in the state, process noise, and output matrices, which have the following form: where , , , , and are known matrices with appropriate dimensions and and are the norm-bounded uncertainties satisfying and .
Due to the delayed measurement model, we need to obtain a concise model for convenience. By defining
we have the following form: where it is known that
According to the definition, can be expressed as where .
For the uncertain system (27), a required filter form is assumed as where is the state estimation value with and and are the filter parameters to be determined. According to the above analysis, the robust filtering problem for delayed uncertain system (24) can be converted to the robust filter design problem for uncertain system (27). Therefore, our aim is to find an upper bound on the estimation error covariance and design a finite-horizon robust filter for (30) to minimize the upper bound.
Remark 2. Compared with the literature [15, 16], it is obvious that the designed robust Kalman filter does not apply to the case that the measurement delay appears in the system. Meanwhile, the definition of the uncertainty matrices and is different from the definition of the corresponding matrices in [15, 16]. So, in order to facilitate the robust filter design, we need to utilize the state augmentation method to obtain a new uncertain system in (27).
3.2. Upper Bound of the Estimation Error Covariance
Because there are uncertain and delay rate terms for the system (27), it is difficult to obtain the true estimation error covariance. Our objective is to find the upper bound , where
Considering the system (27) and the filter structure (30), we define an augmented state . Then, the augmented state-space model is expressed as where The state covariance matrix of in augmented system (32) is denoted as
So, the evolution equation can be expressed as where
In order to obtain the upper bound of the error covariance in (35), the following two lemmas are employed.
Lemma 3 (see ). Given matrices , , , and with compatible dimensions such that , let be a symmetric positive definite matrix and let be an arbitrary positive constant such that Then, the following matrix inequality holds:
Lemma 4 (see ). Let be a real matrix and let be a diagonal random matrix. Then, where is the Hadamard product.
Then, the following conclusion can be given by making use of the two lemmas.
Theorem 5. If there exist three positive scalars , , and , such that and there exists a symmetric positive-definite matrix , such that where with initial value , then for .
Proof. According to Lemma 4, we have
Since the assumptions in (40) hold, using Lemma 3, we can get
Then, when , we have .
When , suppose that .
When , comparing (41) with (44), it is easy to obtain that . According to the mathematical induction, the proof is complete.
Assume that ; then
Therefore, for the upper bound , we need to choose the filter parameters and to obtain an optimized upper bound. In the next part, we will find the value of and design the finite-horizon robust Kalman filter to minimize the upper bound.
3.3. Finite-Horizon Robust Kalman Filter Design
In order to obtain a solution to the above question and design the robust filter, the main result is presented in the following theorem.
Theorem 6. Assume that the positive scalars , , and fulfill the assumptions in (40); then the upper bound can be expressed as If the filter parameters and can be written as where then is minimized. So, the state covariance can be obtained as and the estimation error covariance can be given as
Proof. when , we have . When , assume that (47) is valid. When , we will prove that (47) is still valid. From (47), suppose that the upper bound can be partitioned as . According to the definitions of (50), inserting into (41), we have
Due to the fact that , using (54)–(56), the required upper bound can be calculated as
Computing the first-order variation of (58) with respect to and and making them equal to zero, we have
According to (59) and (60), the optimal parameters and to minimize the required upper bound can be calculated by
Then, substituting (61) and (62) into (55) and (56), we can obtain
Thus, when , (47) is still valid. We can deduce that (50) is valid, for . From (47), we can know that . Utilizing (57), we have
Substituting (64) into (61), (48) can be given. Furthermore, using (59), the required upper bound can be rewritten as (56). Therefore, the theorem has been proved.
Based on the above theorems, the finite-horizon robust Kalman filter can be summarized as follows.
Step 1. The initial values can be given as and .
Remark 7. The finite-horizon robust Kalman filter design is accomplished by using Theorems 5 and 6 for uncertain attitude estimation system with star sensor delays. Different from the most existing attitude estimation filtering algorithms, the finite-horizon robust Kalman filter presented in this paper has the advantage to consider the misalignment errors and scale factor errors of gyros and measurement delays of star sensors for attitude estimation system. Note that these phenomena of the misalignment errors of gyro and star sensor and star sensor measurement delays are often encountered in real attitude estimation systems. To compensate the misalignment errors of sensor and star sensor measurement delays, we have designed a finite-horizon robust Kalman filter by finding the upper bound of the estimation error covariance and minimizing the upper bound. It is worth mentioning that, due to the presence of the star sensor measurement delays, the upper bound of the estimation error covariance in (35) and the filter parameters and in (48) and (49) distinguished our work from the counterpart in [15, 16]. In addition, this paper talks about the filter problem with only one type of noise disturbance. The filter problem with multiple disturbances can be considered to achieve more practical oriented results, as discussed in [28, 29], which will be one of our future research topics.
4. Simulations and Analysis
4.1. Simulation Conditions
The simulation utilizes the data from a satellite. The initial orbit elements of the satellite are set as follows: the semimajor axis , the eccentricity , the inclination , the ascending node longitude , and the argument of perigee . Numerical simulation gives the measurement data of gyro. The standard deviation of gyros’ measurement noise is ; the standard deviation of gyros’ drift noise is ; the gyro sampling interval is ; the components of the gyro scale factor error vector are chosen randomly at the interval ; the misalignment error of gyro is chosen randomly at the interval ; three star sensors are used and the standard deviation of star sensors’ measurement noise is all . The misalignment error of star sensor is set to . Three star sensors have different delay rates, so the random variables satisfy the Bernoulli distribution with
The initial attitude quaternion of the system is taken as ; the gyros’ initial bias is set as ; all the filters are initialized with no attitude errors and zero bias estimate; the initial attitude error covariance is set at for the quaternion components and for the bias components. In order to complete the robust filter design and ensure the estimation precision, the parameters , , and are set to satisfy the condition (40).
4.2. Simulation Results
To validate the effectiveness of the proposed robust filter for controlling the measurement errors of gyros and star sensors and star sensor delays, the proposed approach (FRKF) is compared with the traditional Kalman filter (KF) and the robust Kalman filter (RKF) in the literature . For a fair comparison, the root-mean square error (RMSE) and accumulative RMSE (ARMSE) [30, 31] of the attitude are employed to describe the quality of the attitude estimation. Monte-Carlo simulation runs are set as , and the RMSE of attitude angles can be defined by where expresses the attitude estimation error vector at the th Monte-Carlo run. Then, the ARMSE of the attitude is defined by where denotes the simulation time. The simulation results are shown in Figures 2–4.
From Figures 2 and 3, it is obvious to be seen that the FRKF performs much better than the RKF and KF, and the KF performs the worst. This is because the traditional KF is not suitable for handling the model uncertainties in the system model and star sensor delays. However, the RKF compensates the measurement errors of gyros and star sensors, whose precisions are better than the KF. But the RKF cannot control the effect of the star sensor delays. Compared with the RKF, the FRKF has higher estimation precision that the RKF, which indicates that the FRKF not only deals with the measurement errors of gyros and star sensors but also cope with the star sensor delays.
For the sake of accounting for the effect of star sensor delays, three star sensors are assumed to have the same delay rate, such as . Figure 4 shows the ARMSE of the attitude angles from three filters when . From Figure 4, it can be seen that if there are no star sensor delays, that is, , the FRKF is equal to the RKF. Meanwhile, it is apparent that the KF and RKF increase faster than the FRKF as the delay probability is greater, which illustrates the efficiency of the FRKF to control the star sensor delays.
By the fact that the misalignment errors and scale factor errors of gyros and the misalignment errors of star sensors are difficult to be removed entirely by the attitude estimation filter calibration, these measurement errors of sensors are assumed as the norm-bounded uncertainties in the attitude estimation model. At the same time, due to the complicated signal processing of star sensors, the star sensor delay is one of the most important problems for attitude estimation system. Therefore, a finite-horizon robust Kalman filter for the uncertain attitude estimation system with star sensor delays is proposed in this paper. The uncertain attitude estimation model with star sensor delays is constructed, and the finite-horizon robust Kalman filter design is developed. Finally, the applicability and effectiveness of the proposed filter have been demonstrated by the simulation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
- E. J. Lefferts, F. L. Markley, and M. D. Shuster, “Kalman filtering for spacecraft attitude estimation,” Journal of Guidance, Control, and Dynamics, vol. 5, no. 5, pp. 417–429, 1982.
- F. L. Markley, “Attitude error representations for Kalman filtering,” Journal of Guidance, Control, and Dynamics, vol. 26, no. 2, pp. 311–317, 2003.
- J. L. Crassidis, F. L. Markley, and Y. Cheng, “Survey of nonlinear attitude estimation methods,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 1, pp. 12–28, 2007.
- M. D. Shuster, D. S. Pitone, and G. J. Bierman, “Batch estimation of spacecraft sensor alignments—relative alignment estimation,” Journal of the Astronautical Sciences, vol. 39, no. 4, pp. 519–546, 1991.
- M. Pittelkau, “Kalman filtering for spacecraft system alignment calibration,” Journal of Guidance, Control, and Dynamics, vol. 24, no. 6, pp. 1187–1195, 2001.
- M. Pittelkau, “Survey of calibration algorithms for spacecraft attitude sensors and gyros,” Advances in the Astronautical Sciences, vol. 129, paper no. AAS 07-295, pp. 651–706, 2008.
- K. L. Lai and J. L. Crassidis, “In-space spacecraft alignment calibration using the unscented filter,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, pp. 1–11, 2003.
- J. Vandersteen, M. Diehl, C. Aerts, and J. Swevers, “Spacecraft attitude estimation and sensor calibration using moving horizon estimation,” Journal of Guidance, Control, and Dynamics, vol. 36, no. 3, pp. 734–742, 2013.
- A. Hmamed, C. E. Kasri, E. H. Tissir, T. Alvarez, and F. Tadeo, “Robust H∞ filtering for uncertain 2-D continuous systems with delays,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 5, pp. 2167–2183, 2013.
- C. K. Ahn and P. S. Kim, “New energy-to-peak FIR filter design for systems with disturbance: a convex optimization approach,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 5, pp. 1987–1993, 2013.
- L. G. Wu and D. W. C. Ho, “Fuzzy filter design for Itô stochastic systems with application to sensor fault detection,” IEEE Transactions on Fuzzy Systems, vol. 17, no. 1, pp. 233–242, 2009.
- J. Hu, Z. D. Wang, H. J. Gao, and L. K. Stergioulas, “Extended Kalman filtering with stochastic nonlinearities and multiple missing measurements,” Automatica, vol. 48, no. 9, pp. 2007–2015, 2012.
- J. Hu, Z. D. Wang, B. Shen, and H. J. Gao, “Quantised recursive filtering for a class of nonlinear systems with multiplicative noises and missing measurements,” International Journal of Control, vol. 86, no. 4, pp. 650–663, 2013.
- J. Q. Wang, Z. M. He, H. Y. Zhou, and Y. Y. Jiao, “Regularized robust filter for attitude determination system with relative installation error of star trackers,” Acta Astronautica, vol. 87, pp. 88–95, 2013.
- Z. Dong and Z. You, “Finite-horizon robust Kalman filtering for uncertain discrete time-varying systems with uncertain-covariance white noises,” IEEE Signal Processing Letters, vol. 13, no. 8, pp. 493–496, 2006.
- R. F. Souto and J. Y. Ishihara, “Comments on ‘finite: horizon robust Kalman filtering for uncertain discrete time: varying systems with uncertain: covariance white noises’,” IEEE Signal Processing Letters, vol. 17, no. 2, pp. 213–216, 2010.
- X. Su, L. Wu, and P. Shi, “Sensor networks with random link failures: distributed filtering for T-S fuzzy systems,” IEEE Transactions on Industrial Informatics, vol. 9, no. 3, pp. 1739–1750, 2013.
- X. Su, P. Shi, L. Wu, and S. K. Nguang, “Induced l2 filtering of fuzzy stochastic systems with time-varying delays,” IEEE Transactions on Cybernetics, vol. 43, no. 4, pp. 1251–1264, 2013.
- R. Yang, P. Shi, and G. Liu, “Filtering for discrete-time networked nonlinear systems with mixed random delays and packet dropouts,” IEEE Transactions on Automatic Control, vol. 56, no. 11, pp. 2655–2660, 2011.
- R. Caballero-Aguila, A. Hermoso-Carazo, J. D. Jimenez-Lopez, J. Linares-Perez, and S. Nakamori, “Recursive estimation of discrete-time signals from nonlinear randomly delayed observations,” Computers & Mathematics with Applications, vol. 58, no. 6, pp. 1160–1168, 2009.
- R. Caballero-Aguilaa, A. Hermoso-Carazo, J. D. Jimenez-Lopez, J. Linares-Perez, and S. Nakamori, “Signal estimation with multiple delayed sensors using covariance information,” Digital Signal Processing, vol. 20, no. 2, pp. 528–540, 2010.
- H. H. Dam, “Variable fractional delay filter with sub-expression coefficients,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 7, pp. 2995–3003, 2013.
- F. O. Hounkpevi and E. Yaz, “Minimum variance generalized state estimators for multiple sensors with different delay rates,” Signal Processing, vol. 87, no. 4, pp. 602–613, 2007.
- X. X. Wang, Y. Liang, Q. Pan, and C. H. Zhao, “Gaussian filter for nonlinear systems with one-step randomly delayed measurements,” Automatica, vol. 49, no. 4, pp. 976–986, 2013.
- J. Hu, Z. D. Wang, B. Shen, and H. J. Gao, “Gain-constrained recursive filtering with stochastic nonlinearities and probabilistic sensor delays,” IEEE Transactions on Signal Processing, vol. 61, no. 5, pp. 1230–1238, 2013.
- L. Xie, Y. C. Soh, and C. E. de Souza, “Robust Kalman filtering for uncertain discrete-time systems,” IEEE Transactions on Automatic Control, vol. 39, no. 6, pp. 1310–1314, 1994.
- R. A. Horn and C. R. Johnson, Topics in Matrix Analysis, Cambridge University Press, New York, NY, USA, 1991.
- X. M. Yao and L. Guo, “Composite anti-disturbance control for Markovian jump nonlinear systems via disturbance observer,” Automatica, vol. 49, no. 8, pp. 2538–2545, 2013.
- X. M. Yao, Z. Dong, and D. F. Wang, “Full-order disturbance-observer-based control for singular hybrid system,” Mathematical Problems in Engineering, Article ID 352198, 7 pages, 2013.
- X. J. Tang, Z. B. Liu, and J. S. Zhang, “Square-root quaternion cubature Kalman filtering for spacecraft attitude estimation,” Acta Astronautica, vol. 76, pp. 84–94, 2012.
- I. Arasaratnam, S. Haykin, and T. R. Hurd, “Cubature Kalman filtering for continuous-discrete systems: theory and simulations,” IEEE Transactions on Signal Processing, vol. 58, no. 10, pp. 4977–4993, 2010.