|
|
| Line 1: |
Line 1: |
| The '''invariant extended Kalman filter (IEKF)'''<ref name="Bonnabel_cdc09">S. Bonnabel, Ph. Martin and E. Salaün, "Invariant Extended Kalman Filter: theory and application to a velocity-aided attitude estimation problem", 48th IEEE Conference on Decision and Control, pp. 1297–1304, 2009.</ref> is a new version of the [[extended Kalman filter]] (EKF) for nonlinear systems possessing symmetries (or ''invariances''). It combines the advantages of both the EKF and the recently introduced [[symmetry-preserving filter]]s. Indeed, instead of using a linear correction term based on a linear output error, it uses a geometrically adapted correction term based on an invariant output error; in the same way the gain matrix is not updated from a linear state error, but from an invariant state error. The main benefit is that the gain and covariance equations converge to constant values on a much bigger set of trajectories than equilibrium points than is the case for the EKF, which results in a better convergence of the estimation.
| | Hi, everybody! <br>I'm Dutch male ;=). <br>I really like Amateur geology!<br><br>my web site: [http://www.massivewallpapers.com/profile/sudaniels Popular mountain bike sizing.] |
| | |
| == Motivation ==
| |
| | |
| Most physical systems possess natural symmetries (or invariance), i.e. there exist transformations (e.g. rotations, translations, scalings) that leave the system unchanged. From mathematical and engineering viewpoint, it makes sense that a filter well-designed for the considered system should preserve the same invariance properties. The idea for the IEKF is a modification of the EKF equations to take advantage of the symmetries of the system.
| |
| | |
| == Definition ==
| |
| | |
| Consider the system
| |
| | |
| :<math>
| |
| \begin{align}
| |
| \dot x&=f(x,u)+M(x) w \\
| |
| y&=h(x,u) +N(x)v
| |
| \end{align}
| |
| </math>
| |
| | |
| where <math> w,v </math> are independent [[white Gaussian noise]]s.
| |
| Consider <math>G</math> a [[Lie group]] with identity <math>e</math>, and
| |
| (local) [[transformation group]]s <math>\varphi_g, \psi_g, \rho_g</math> (<math>g \in G</math>) such that <math>(X,U,Y)=(\varphi_g(x),\psi_g(u),\rho_g(y))</math>. The previous system with noise is said to be ''[[Invariant (mathematics)|invariant]]'' if it is left unchanged by the action the transformations groups <math>\varphi_g, \psi_g, \rho_g</math>; that is, if
| |
| | |
| :<math> | |
| \begin{align}
| |
| \dot X&=f(X,U)+M(X) w \\
| |
| Y&=h(X,U) +N(X)v.
| |
| \end{align}
| |
| </math>
| |
| | |
| == Filter equations and main result ==
| |
| | |
| Since it is a [[symmetry-preserving filter]], the general form of an IEKF reads <ref name="Bonnabel_tac">S. Bonnabel, Ph. Martin, and P. Rouchon, “Symmetry-preserving observers,”
| |
| ''IEEE Transaction on Automatic and Control'', vol. 53, no. 11, pp. 2514–2526, 2008.</ref>
| |
| | |
| :<math>\dot{\hat x}=f(\hat x,u)
| |
| +W(\hat x)L\Bigl(I(\hat x,u),E(\hat x,u,y)\Bigr)E(\hat x,u,y)</math>
| |
| | |
| where
| |
| | |
| * <math>E(\hat x,u,y)</math> is an invariant output error, which is different to the usual output error <math>\hat y-y</math>
| |
| | |
| *<math>W(\hat x)=\bigl(w_1(\hat x),..,w_n(\hat x)\bigr)</math> is an invariant frame
| |
| | |
| *<math>I(\hat x,u)</math> is an invariant vector
| |
| | |
| *<math>L(I,E)</math> is a freely chosen gain matrix.
| |
| | |
| To analyze the error convergence, an invariant state error <math> \eta(\hat x,x)</math> is defined, which is different from the standard output error <math> \hat x-x </math>, since the standard output error usually does not preserve the symmetries of the system.
| |
| | |
| Given the considered system and associated transformation group, there exists a constructive method to determine <math> E(\hat x,u,y),W(\hat x),I(\hat x,u),\eta(\hat x,x)</math>, based on the moving frame method.
| |
| | |
| Similarly to the EKF, the gain matrix <math>L(I,E)</math> is determined from the equations<ref name="Bonnabel_cdc09" />
| |
| | |
| :<math> | |
| \begin{align}
| |
| L&=PC^T\bigl(N(e)N^T(e)\bigr)^{-1}\\
| |
| \dot P&=AP+PA^T+M(e)M^T(e)-PC^T\bigl(N(e)N^T(e)\bigr)^{-1}CP,
| |
| \end{align}
| |
| </math>
| |
| | |
| where the matrices <math>A,C</math> depend here only on the known invariant vector <math>I(\hat x,u)</math>, rather than on <math>(\hat x,u)</math> as in the standard EKF. This much simpler dependence and its consequences are the main interests of the IEKF. Indeed, the matices <math>A,C</math> are then constant on a much bigger set of trajectories (so-called ''permanent trajectories'') than equilibrium points as it is the case for the EKF. Near such trajectories, we are back to the "true", i.e. linear, Kalman filter where convergence is guaranteed. Informally, this means the IEKF converges in general at least around any slowly varying permanent trajectory, rather than just around any slowly varying equilibrium point for the EKF.
| |
| | |
| == Application example in aerospace engineering==
| |
| | |
| Invariant extended Kaman filters are for instance used in [[Attitude and Heading Reference Systems|attitude and heading reference systems]]. In such systems the orientation, velocity and/or position of a moving rigid body,
| |
| e.g. an aircraft, are estimated from different embedded sensors, such as inertial sensors, magnetometers, GPS or sonars. The use of an IEKF naturally leads<ref name="Bonnabel_cdc09" /> to consider the [[quaternion]] error <math>\hat qq^{-1}</math>, which is often used as an ''ad hoc'' trick to preserve the constraints of the quaternion group. The benefits of the IEKF compared to the EKF are experimentally shown for a large set of trajectories.<ref name="Martin_gnc10">Ph. Martin and E. Salaün, "Generalized Multiplicative Extended Kalman Filter for Aided Attitude and Heading Reference System", AIAA Guidance, Navigation and Control Conference, 2010</ref>
| |
| | |
| == References ==
| |
| | |
| <references/>
| |
| | |
| {{DEFAULTSORT:Invariant Extended Kalman Filter}}
| |
| [[Category:Nonlinear filters]]
| |
| [[Category:Estimation theory]]
| |
Hi, everybody!
I'm Dutch male ;=).
I really like Amateur geology!
my web site: Popular mountain bike sizing.