The accuracy of autonomous orbit determination of Lagrangian navigation constellation will affect the navigation accuracy for the deep space probes. Because of the special dynamical characteristics of Lagrangian navigation satellite, the error caused by different estimation algorithm will cause totally different autonomous orbit determination accuracy. We apply the extended Kalman filter and the fading–memory filter to determinate the orbits of Lagrangian navigation satellites. The autonomous orbit determination errors are compared. The accuracy of autonomous orbit determination using fading-memory filter can improve 50% compared to the autonomous orbit determination accuracy using extended Kalman filter. We proposed an integrated Kalman fading filter to smooth the process of autonomous orbit determination and improve the accuracy of autonomous orbit determination. The square root extended Kalman filter is introduced to deal with the case of inaccurate initial error variance matrix. The simulations proved that the estimation method can affect the accuracy of autonomous orbit determination greatly.
basic scientific research fund of national defense2016110C019Natural Science Foundation of Jiangsu ProvinceBK20160811Shanghai Deep Space Detection Technology Key Laboratory Open FundingDS2016-01specialized research fund for the doctoral program of higher education of china201332181200371. Introduction
Deep space exploration has become a hot spot of aerospace. Several deep space probes have been launched. The autonomous navigation is important for deep space probes to deal with communication delay as well as reducing the dependency on ground stations. As early as 1968, the sextant had been used for autonomous navigation in “Apollo program” [1]. In 1999, “Deep Space 1” achieved autonomous orbit determination by tracking small celestial bodies with an optical sensor, which is the first successful on-orbit application of the deep space autonomous navigation technology [2]. The comet probe “Deep Impact” which was launched in 2005 also carried out its navigation and control automatically based on an optical navigation system with a high resolution imager [3]. In [4], Downs proposed to use X-ray pulsar radiation signal to navigate spacecraft. The rotation period of X-ray pulsar is extremely stable; therefore, time and the location of spacecraft can be determined by tracking several X-ray pulsars with given and fixed frequency [5]. The satellite navigation constellation can also provide navigation information for deep space probes. GPS can navigate the deep space probes when they are running in low-Earth orbits and medium Earth orbits. For deep space transfer orbits and deep space target orbits, the GPS is not good enough. Several researchers investigated weak GNSS signal navigation for the deep space probes [6–8]. Witternigg et al. introduced how GPS and Galileo could be used for orbit determination in future missions to the Moon [8]. Farquhar introduced a concept of using Earth-Moon libration point satellites for lunar navigation [9, 10]. In 2005, Hill suggested placing navigation constellation on the periodic orbits in the vicinity of libration points of the Earth-Moon system to support deep space navigation [11]. Zhang and Xu analyzed the architecture and navigation performance of the Lagrangian point satellite navigation system [12–14]. The Lagrangian navigation constellation is introduced to navigate the deep space probes autonomously. Hence, the navigation constellation itself should have the ability of autonomous orbit determination (AOD). The methods introduced in [1–5] can be considered as absolute navigation (or absolute autonomous orbit determination) because the estimated orbit refers to an inertial or quasi-inertial frame. Methods introduced in [6–14] can be classified as relative navigation. Relative navigation seeks optimal estimates for the position and velocity of one satellite relative to the other one.
Relative navigation usually is applied to satellites in a formation or constellation involved using GPS which restricts the spacecraft formation to near-Earth applications, such as Deep Space Mission 3 [15] and Grace project [16]. Relative navigation is primarily proposed for formation configuration control and formation reconfiguration. However, Lagrangian navigation constellation should provide absolute navigation information to deep space probes to achieve absolute navigation. Therefore, autonomous orbit determination of the Lagrangian navigation constellation is actually using the relative measurement to achieve absolute navigation. The laser interferometer space antenna (LISA) mission is an example which uses the relative range to assist the absolute orbit determination. The LISA mission which consists of three spacecraft separated by 5 million kilometers forming an equilateral triangle is a huge Michelson interferometer in space for gravitational wave detection [17]. The deep space-network provides a raw estimation of the absolute positioning for the three satellites. An accurate relative positioning will be provided with a laser-based ranging measurement in order to obtain an accuracy of positioning of tens of meters [18, 19]. Psiaki [20] and Markley [21] suggested using crosslink range, attitude information, and an optical tracker to determine the orbits autonomously. Yim et al. proposed using optical tracking and attitude information to find the direction vector between the two spacecraft and determine both orbits [22]. But these methods require extensive hardware development. In order to reduce the operational cost, size, and weight of spacecraft for formation missions, crosslink range can be used as the only measurement for orbit determination of a constellation. However, for the Earth navigation satellite constellation, there is a rank deficiency problem when only crosslink range is used to determine the orbit [23, 24]. Hill’s study illustrated that the rank deficiency problem does not exist for the Lagrangian navigation satellites because of the special dynamics near the libration points [11]. Thus, the Lagrangian navigation satellites can autonomously determine their orbits using only crosslink range. In [11], Hill discussed the unique distribution of Lagrangian orbit from the view of dynamics, which theoretically proved the autonomy of the Lagrangian navigation constellation. From the perspective of identifiability of epoch state, Qian et al. verified the feasibility of AOD for satellites in quasiperiodic orbits about the Earth-Moon libration point [25]. Based on circular restricted three-body problem (CR3BP), Du et al. researched the autonomous orbit determination method of satellites in halo orbits, and only crosslink range was used as observation [26]. In [23], Gao et al. discussed the feasibility of autonomous orbit determination using only the crosslink range measurement for a combined Lagrangian navigation constellation and GNSS. The most widely used algorithm for Lagrangian navigation satellite is extended Kalman filter (EKF) method. However, EKF is based on linearization of the system dynamics and the assumption of Gaussian process/measurement noise. These can seriously affect the performance of the state estimation and even lead to divergence. Unscented Kalman filter (UKF) which is based on the unscented transform can achieve higher accuracy than EKF while the added computational cost is not significant. More importantly, UKF is robust with respect to the initial conditions. Therefore, Sun et al. introduced UKF to relative navigation for multiple spacecraft formation flying [27]. Giannitrapani et al. analyze the performance of EKF and UKF for the localization of a spacecraft [28]. In order to improve the robustness and stability accuracy, Wang and Gu applied fault tolerant UKF in autonomous determination of relative orbit for satellite formation flying [29]. Rigatos introduced the technical analysis and implementation cost assessment of sigma-point Kalman filtering and particle filtering in autonomous navigation systems [30]. Reali and Palmerini provided a preliminary comparison of different estimation techniques to be used in formation flying navigation [31].
Since the CR3BP is sensitive to the state error and calculation error, the AOD of Lagrangian navigation satellite may refer to the accuracy of the estimation algorithms. One factor which must be considered in the AOD of Lagrangian navigation constellation is to prevent the divergence of AOD error. Therefore, we introduce four estimation methods to achieve the AOD of Lagrangian navigation constellation to analyze the effect on AOD by estimation method.
2. Dynamical Model of Lagrangian Navigation Satellite
For satellites in Lagrangian point orbits, the equation of CR3BP should be an appropriate model to describe the satellites’ dynamical characteristics. Consider two massive bodies m1 and m2 moving under the action of just their mutual gravitation, and let their orbit around each other be a circle of radius r12. As shown in Figure 1, a noninertial, comoving frame of reference o-xyz is defined. The origin of frame o-xyz lies at the center of mass of the two-body system. The positive x direction goes from m1 to m2. The positive y-axis is parallel to the velocity vector. The z-axis is perpendicular to the orbital plane. Now the third body of mass m which is vanishingly small compared to the primary masses m1 and m2 is introduced. We assume that the mass m is so small that it has no effect on the motion of the primary bodies. This is called the restricted three-body problem.
Diagram of the circular restricted three-body problem in o-xyz frame.
The nondimensional equations of motion for the CR3BP are shown as follows [32]:(1)x¨-2y˙=x-1-μx+μr13-μx+μ-1r23,y¨+2x˙=1-1-μr13-μr23y,z¨=μ-1r13-μr23z,where m1 and m2 are the masses of the two primaries. r1 and r2 are the distances of the spacecraft from the massive and secondary primaries, respectively. r1, r2, and μ are given by (2)r1=x+μ2+y2+z2,r2=x+μ-12+y2+z2,μ=m2m1+m2.
The nondimensional scales are defined as(3)M=m1+m2,L=r12,T=r123Gm1+m21/2,where r12 is the distance between the two primaries. G is the Universal Constant of Gravitation.
Equations (1) have five equilibrium points which are called Lagrange points (or libration point), as shown in Figure 2.
Schematic location of the five Lagrange points in the CR3BP.
Three Lagrange points on the x-axis are unstable, and two Lagrange points which form an equilateral triangle with the two primary bodies in the x-y plane are stable. Many interesting periodic orbits exist in the vicinity of the Lagrange points no matter whether they are stable. The Lagrange navigation satellite constellation discussed in this paper is distributed in these periodic orbits.
3. AOD of Lagrangian Navigation Constellation with Different Estimation Methods
EKF method is a commonly used method for AOD of satellite. EKF is summarized as follows [33]:(4)X^k=Φk,k-1X^k-1+Kkρk-HkΦk,k-1X^k-1,Kk=PkHkTRk-1,Pk=I-KkHkPk,k-1,Pk,k-1=Φk,k-1Pk-1Φk,k-1T+Γk,k-1Qk-1Γk,k-1T,where X^k-1 is the estimated state at tk-1. Φtk,tk-1 is the transition matrix. Pk is the covariance matrix. Hk is the mapping matrix that relates the observation deviation vector to the state deviation vector. Kk is the Kalman gain matrix. In this paper, only crosslink range between two satellites is used as observation. The crosslink range is described as(5)ρi,j=xi-xj2+yi-yj2+zi-zj2.So (6)Hk=∂ρxk,k∂xkTxk,k-1.Two periodic orbits around L1 and L2, respectively, are chosen to achieve the simulation. The initial states of the two orbits are shown in Table 1. The values are nondimensional. 10 m initial state errors and 10 m measurement noise are added. In Figures 3 and 4, the orbits based on different initial state is shown. As we can see, the trajectory will diverge within a short time if the initial state error is 10 m. Therefore, the initial state error and the measurement noise are chosen based on the navigation requirement, not the ability of the existing technology.
Initial states of the periodic orbits around Lagrange points L1 and L2.
Satellites
L1
L2
x(L)
0.8440212401521470
1.1726789595745000
y(L)
0.0000000000000000
0.0000000000000000
z(L)
0.0592695845762629
0.0834298981288340
x˙(L/T)
0.00000000000000000
0.0000000000000000
y˙(L/T)
-0.0053009560909076
-0.1864489128036080
z˙(L/T)
0.0000000000000000
0.0000000000000000
The periodic orbit based on accurate initial state.
The orbit based on the initial state with 10 m error.
The initial covariance matrix is defined as(7)Pt0=σx02000000σy02000000σz02000000σx˙02000000σy˙02000000σz˙02,where σx0, σy0, and σz0 are the initial position error in three coordinate axes and σx˙0, σy˙0, and σz˙0 are the initial velocity error in three coordinate axes.
Figures 5 and 6 show the results of AOD of Lagrangian navigation satellite using EKF. As we can see, the maximum errors in three axes of L1 are 541 m, 254 m and 168 m. The maximum errors in three axes of L2 are 555 m, 195 m, and 263 m. And the estimation process is divergent.
The AOD position error of L1 using EKF.
The AOD position error of L2 using EKF.
Since CR3BP is sensitive to calculation error, the state estimation error covariance matrix of EKF may lose its properties of positive definite when the computations are carried out with finite digit arithmetic of the computer. Nonpositive definite covariance matrix will cause the gain matrix to lose its function gradually and lead to divergence of AOD error. In order to suppress the divergence of estimation process, the fading-memory filter (FMF) is introduced to improve the AOD characteristics. The FMF is essentially the old estimate plus a gain times a residual (difference between current measurement and previous estimate). The difference between the EKF and FMF is that the FMF is recursive and weights new measurements more heavily than older measurements [34].
The basic equations of FMF are listed as follows [33]:(8)X^k=Φk,k-1X^k-1+Kkρk-HkΦk,k-1X^k-1,Kk=PkHkTRk-1,Pk=I-KkHkPk,k-1,Pk,k-1=Φk,k-1Pk-1Φk,k-1TeCk-1+Γk,k-1Qk-1Γk,k-1T.FMF raises the role of new observation data and reduces the influence of old observation to the filter in order to suppress the filtering divergence. Figures 7 and 8 show the AOD position errors of L1 and L2 using FMF. It can be seen that the divergence of errors is suppressed by using FMF method. The accuracy of the AOD is improved as well.
The AOD position error of L1 using FMF.
The AOD position error of L2 using FMF.
One disadvantage of FMF is that the oscillation amplitude of filtering error at early stage is larger. Comparing Figure 5 with Figure 7, we can see that estimation process of EKF is relatively ideal during the first 100 s. Therefore, we can use EKF to determine the orbit of Lagrangian navigation satellite at the early stage of the AOD. When the AOD error is bigger than a threshold value, the estimation method will be switching to FMF. We call this method integrated Kalman fading filter (IKFF) which is illustrated as (9)X^k=Φk,k-1X^k-1+Kkρk-HkΦk,k-1X^k-1,Kk=PkHkTRk-1,Pk=I-KkHkPk,k-1,Pk,k-1=Φk,k-1Pk-1Φk,k-1T+Γk,k-1Qk-1Γk,k-1T,if e≥es,Φk,k-1Pk-1Φk,k-1TeCk-1+Γk,k-1Qk-1Γk,k-1T,if e<es,where es is the threshold value of the AOD error. Since the estimation error cannot be got in practical applications, the residual can be used to approximate the threshold value of estimation error. The linearized relationship between the residual and the AOD error can be expressed as(10)y=Ce,where y is the residual and C is denoted as(11)C=x^-x^2ρy^-y^2ρz^-z^2ρ000.x^y^z^T and x^2y^2z^2T are the estimated position variable of Lagrangian satellite 1 and Lagrangian satellite 2, respectively. ρ is the crosslink range between these two satellites.
The AOD results of IKFF are shown in Figures 9 and 10. We list the maximum error of the above three methods in Tables 2 and 3. As shown in Table 2, the maximum error of L1 satellite in three axes using IKFF reduced to 204 m, 104 m, and 62 m. Meanwhile, oscillation characteristic at the early stage of AOD is improved in Figure 9. The similar improvement of AOD can be seen about L2 satellite. Therefore, as long as the threshold value is selected appropriately, IKFF will produce better AOD results than FMF or EKF.
The maximum AOD error of L1 using different methods.
Maximum error in x-axis/m
Maximum error in y-axis/m
Maximum error in z-axis/m
EKF
541
253
168
FMF
273
129
92
IKFF
204
104
62
The maximum AOD error of L2 using different methods.
Maximum error in x-axis/m
Maximum error in y-axis/m
Maximum error in z-axis/m
EKF
554
195
262
FMF
296
226
163
IKFF
215
81
128
The AOD position error of L1 using IKFF.
The AOD position error of L2 using IKFF.
To overcome the divergence of AOD error caused by the state estimation error covariance matrix losing its properties of positive definite symmetric, we also introduce the square root extended Kalman filter (SR-EKF) to estimate the state of the Lagrangian satellite using only crosslink range. The SR-EKF is described as follows [33]:(12)X^k=Φk,k-1X^k-1+Kkρk-HkΦk,k-1X^k-1,X^0=EX0,Kk=akSk,k-1Fk,Fk=Sk,k-1THkT,ak=FkTFk+Rk-1,Sk,k-1=Φk,k-1Sk-1,P0=S0S0T,Sk=Sk,k-1I-akrkFkFkT,rk=11±akRk.The state estimation error covariance matrix P is replaced by SST where S is the square root of P [35]. Here we define the order of σx02, σy02, and σz02 as P0. This method can guarantee the P is symmetric nonnegative definite matrix at any time. In Tables 4 and 5, we give the maximum error of AOD of L1 satellite and L2 satellite with different initial state estimation error covariance matrix. As we can see, with a small initial value of error variance matrix, SR-EKF algorithm cannot estimate the state very well. The estimation error using SR-EKF is bigger than the EKF results. However, the accuracy of AOD is improved with the increase of initial error covariance matrix. Thus, the effectiveness of SR-EKF method is affected by the initial value of state estimation error covariance matrix P. SR-EKF is useful when the initial error variance matrix is not accurate. However, the AOD error will increase when the initial error variance matrix is bigger than a boundary. In our simulation, the best AOD result coincides with the P0 = 10−10. Figures 11 and 12 exhibit the AOD error using SR-EKF under the same initial condition as the above three methods with P0 = 10−10. Even when the initial error variance matrix is not accurate, the precision of AOD of two Lagrangian satellites is better than EKF.
The maximum AOD error of L1 using SR-EKF with different P0.
Maximum error in x-axis/m
Maximum error in y-axis/m
Maximum error in z-axis/m
P0=10-16
17936
8507
4378
P0=10-14
853
393
184
P0=10-12
404
186
88
P0=10-10
377
178
104
The maximum AOD error of L2 using SR-EKF with different P0.
Maximum error in x-axis/m
Maximum error in y-axis/m
Maximum error in z-axis/m
P0=10-16
17290
5216
6532
P0=10-14
790
283
301
P0=10-12
374
134
154
P0=10-10
365
107
134
The AOD position error of L1 using SR-EKF with P0=10-10.
The AOD position error of L2 using SR-EKF with P0=10-10.
4. Conclusion
Since the special dynamical characteristics of CR3BP, the AOD accuracy of the Lagrangian navigation constellation is sensitive to all the errors introduced by the AOD process. Different estimation methods will cause great different AOD accuracy.
Four estimation methods are used to determinate the orbit of Lagrangian navigation constellation. EKF is a commonly used method, but it is not the best choice for the AOD of Lagrangian navigation constellation. The accuracy of AOD using FMF can improve 50% more than using EKF. And the divergence of the AOD error is suppressed.
We proposed a new method, IKFF, to estimate the state of Lagrangian navigation satellite. The accuracy of AOD using IKFF can improve 20% more than using FMF. Furthermore, IKFF can smooth the process of the AOD.
For the case of not accurate initial error variance matrix, SR-EKF can be chosen as the AOD estimation method. The accuracy of SR-EKF is close to the FMF.
In order to improve the AOD accuracy, all the factors which will cause errors should be analyzed. As the estimation methods, we will find other advanced technologies to get a much more accurate AOD of Lagrangian navigation satellite.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This study was supported by the basic scientific research fund of national defense (no. 2016110C019), the Natural Science Foundation of Jiangsu Province (no. BK20160811), Shanghai Deep Space Detection Technology Key Laboratory Open Funding (no. DS2016-01), and the specialized research fund for the doctoral program of higher education of china (no. 20133218120037).
BattinR. H.1999Reston, Va, USAAmerican Institute of Aeronautics and Astronautics (AIAA)AIAA Education Series10.2514/4.861543MR1755760BhaskaranS.RiedelJ. E.SynnottS. P.WangT. C.The deep space 1 autonomous navigation system: a post-flight analysisProceedings of the AIAA/AAS Astrodynamics Specialist ConferenceAugust 2000Denver, Colo, USA4252AIAA-2000-39352-s2.0-62949087398MastrodemosN.KubitschekD. G.SynnottS. P.Autonomous navigation for the Deep Impact mission encounter with comet Tempel 120051171-29512110.1007/s11214-005-3394-42-s2.0-22044458422DownsG. S.Interplanetary navigation using pulsating radio sources1974N74-34150/4SheikhS. I.PinesD. J.RayP. S.WoodK. S.LovelletteM. N.WolffM. T.Spacecraft navigation using X-ray pulsars2006291496310.2514/1.133312-s2.0-32544432047SilvaP. F.LopesH. D.PeresT. R.SilvaJ. S.OspinaJ.CichockiF.DovisF.MusumeciL.SerantD.CalmettesT.PessinaI.PerellóJ. V.Weak GNSS signal navigation to the moonProceedings of the 26th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS '13)September 2013Nashville, Tenn, USA335733672-s2.0-84899026165Manzano-JuradoM.Alegre-RubioJ.PellacaniA.Seco-GranadosG.López-SalcedoJ. A.GuerreroE.García-RodríguezA.Use of weak GNSS signals in a mission to the moonProceedings of the 7th ESA Workshop on Satellite Navigation Technologies and European Workshop on GNSS Signals and Signal Processing (NAVITEC '14)2014WitterniggN.ObertaxerG.SchönhuberM.PalmeriniG. B.RodriguezF.CapponiL.SoualleF.FlochJ.Weak GNSS signal navigation for Lunar exploration missionsProceedings of the 28th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS '15)September 2015Tampa, Fla, USA39283944FarquharR. W.Lunar communications with libration-point satellites19674101383138410.2514/3.29095FarquharR. W.The control and use of libration point satellites1970TR R-346HillK. A.2007Boulder, Colo, USAUniversity of ColoradoZhangL.XuB.A universe light house—candidate architectures of the libration point satellite navigation system201467573775210.1017/s03734633140001372-s2.0-84905513664ZhangL.XuB.Navigation performance of the libration point satellite navigation system in cislunar space201568236738210.1017/S03734633140006172-s2.0-84922578598ZhangL.XuB.Navigation performance of the libration point satellite navigation system for future Mars exploration2016691415610.1017/s03734633150004782-s2.0-84949533104LauK.LichtenS.YoungL.HainesB.An innovative deep space application of GPS technology for formation flying spacecraftProceedings of the Guidance, Navigation, and Control ConferenceJuly 1996San Diego, Calif, USAAmerican Institute of Aeronautics and Astronautics10.2514/6.1996-3819LongY.XurongD.DongmeiX.Precise navigation for GRACE formation flyingProceedings of the 1st International Symposium on Systems and Control in Aerospace and AstronauticsJanuary 2006Harbin, China7927972-s2.0-33750940622The LISA Study Team1998Max-Planck-Institute for Quantum OpticsEsteban DelgadoJ. J.Garcia MarínA. F.BykovI.HeinzelG.DanzmannK.Free-space laser ranging and data communicationProceedings of the 6th Workshop on Positioning, Navigation and Communication (WPNC '09)March 2009Hannover, GermanyIEEE27528110.1109/wpnc.2009.49078382-s2.0-70349137947EstebanJ. J.GarcíaA. F.EichholzJ.PeinadoA. M.BykovI.KullmannJ.HeinzelG.DanzmannK.Optical ranging and data communication in space-based applicationsProceedings of the 7th Workshop on Positioning, Navigation and Communication (WPNC '10)March 2010Dresden, Germany192210.1109/wpnc.2010.56537832-s2.0-78651417426PsiakiM. L.Autonomous orbit determination for two spacecraft from relative position measurements199922230531210.2514/2.43792-s2.0-0033098945MarkleyF. L.Autonomous navigation using landmark and intersatellite dataProceedings of the AIAA/AAS Astrodynamics Conference1984Seattle, Wash, USAAIAA Paper 1984–19872-s2.0-0021304837YimJ. R.CrassidisJ. L.JunkinsJ. L.Autonomous orbit navigation of two spacecraft system using relative line of sight measurementsProceedings of the AAS/AIAA Astrodynamics Specialist Conference2004Maui, Hawaii, USAAAS Paper 04-257GaoY. T.XuB.ZhangL.Feasibility study of autonomous orbit determination using only the crosslink range measurement for a combined navigation constellation20142751199121010.1016/j.cja.2014.09.0052-s2.0-84922559743HillK.BornG. H.Autonomous interplanetary orbit determination using satellite-to-satellite tracking200730367968610.2514/1.245742-s2.0-34249822819QianY.-J.JingW.-X.GaoC.-S.WeiW.-S.Autonomous orbit determination for quasi-periodic orbit about the translunar libration point201334562563310.3873/j.issn.1000-1328.2013.05.0052-s2.0-84879174376DuL.ZhangZ.YuL.ChenS.SST orbit determination of halo-LMO constellation in CRTBP20134221841902-s2.0-84877913630SunD.ZhouF. Q.ZhouJ.Relative navigation based on UKF for multiple spacecraft formation flyingProceedings of the AIAA Guidance, Navigation, and Control Conference and ExhibitAugust 2004Providence, RI, USAGiannitrapaniA.CeccarelliN.ScortecciF.GarulliA.Comparison of EKF and UKF for spacecraft localization via angle measurements2011471758410.1109/TAES.2011.57056602-s2.0-79551615434WangH.GuX. M.Autonomous determination of relative orbit for satellite formation flying using fault tolerant UKFProceedings of the IEEE International Conference on Mechatronics and Automation (ICMA '09)August 2009Changchun, China3905390910.1109/icma.2009.52448882-s2.0-77449132397RigatosG. G.Technical analysis and implementation cost assessment of sigma-point Kalman filtering and particle filtering in autonomous navigation systemsProceedings of the 2010 IEEE 71st Vehicular Technology Conference (VTC '10-Spring)May 2010Taipei, TaiwanIEEE10.1109/vetecs.2010.54936392-s2.0-77954939451RealiF.PalmeriniG.Estimate problems for satellite clustersProceedings of the IEEE Aerospace Conference (AC '08)March 2008Big Sky, Mont, USA10.1109/aero.2008.45262722-s2.0-49349095811LiuL.HouX. Y.2012Beijing, ChinaPublishing House of Electronics IndustryFuM. Y.DengZ. H.YanL. P.2010Beijing, ChinaScience PressZarchanP.MusoffH.20154thAmerican Institute of Aeronautics and Astronautics(AIAA)Progress in Astronautics and AeronauticsVan Der MerweR.WanE. A.The square-root unscented Kalman filter for state and parameter-estimationProceedings of the International Conference on Acoustics, Speech, and Signal ProcessingMay 2001Salt Lake City, Utah, USA