Robust LCEKF for Mismatched Nonlinear Systems with Non-Additive Noise/Inputs and Its Application to Robust Vehicle Navigation

It is well known that the standard state estimation technique performance is particularly sensitive to perfect system knowledge, where the underlying assumptions are: (i) Process and measurement functions and parameters are known, (ii) inputs are known, and (iii) noise statistics are known. These are rather strong assumptions in real-life applications; therefore, a robust filtering solution must be designed to cope with model misspecifications. A possible way to design robust filters is to exploit linear constraints (LCs) within the filter formulation. In this contribution we further explore the use of LCs, derive a linearly constrained extended Kalman filter (LCEKF) for systems affected by non-additive noise and system inputs, and discuss its use for model mismatch mitigation. Numerical results for a robust tracking and navigation problem are provided to show the performance improvement of the proposed LCEKF, with respect to state-of-the-art techniques, that is, a benchmark EKF without mismatch and a misspecified EKF not accounting for the mismatch.


Introduction
The design and use of state estimation techniques is fundamental in a plethora of applications, such as robotics, tracking, guidance and navigation systems [1][2][3][4]. For a linear dynamic system, the Kalman filter (KF) is the the best linear minimum mean square error (MSE) estimator. The most widespread solution for nonlinear systems is to resort to system linearisations, leading to the so-called linearised or extended KF (EKF) [3]. In both cases, as well as for more advanced techniques such as sigma-point filters [5], the main assumption is perfect system knowledge: (i) Known process and measurement functions, and their parameters, (ii) known inputs, and (iii) known noise statistics (i.e., first and second order moments for the KF and EKF). However, these are rather strong assumptions in real-life applications, where the noise statistics' parameters may be unknown to a certain extent, inputs may be uncertain and system parameters may be misspecified. The performance degradation of minimum MSE estimators under model mismatch has been widely reported in the literature [6][7][8][9][10], which is a reason why there exists a need to develop robust filtering techniques able to cope with mismatched systems.
Notice that in the robust filtering literature, a lot of effort has been devoted to the mitigation of possible noise statistics' parameters misspecification, that is, to coping with non-nominal or unexpected noise behaviours. In that perspective, we may cite the robust KFs able to counteract the presence of outliers [11], the ones that estimate Gaussian noise covariances [12,13], or alternative formulations to deal with heavy-tailed distributions [14,15]. In contrast, few contributions explored how to counteract a mismatch on system matrices/functions or the filter initialisation. Within the KF framework, the latter can be solved by either the information filter form of the KF ([3] §6.2) or the so-called Fisher initialisation [16], which can be generalised by imposing initial distortionless constraints [17], i.e., the so-called minimum variance distortionless response (MVDR) estimators. In order to further generalise these MVDR results, how to incorporate non-stationary constraints within the KF has been recently proposed by Villà-Valls et al. [18], leading to a general linearly constrained KF (LCKF) formulation, which has also been shown concerning linear systems to generalise the results proposed by Teixeira et al. [19]. Indeed, the use of linear constraints (LCs) has been shown to be a promising robust filter design solution in order to mitigate the impact of misspecified system matrices [18], which complements the existing results cited above for noise parameter mismatch mitigation. Notice that other alternatives exist to estimate unknown system parameters/inputs, such as augmenting the filter state, which are not discussed in this contribution.
Once the general LCKF has been established, it is of interest to extend its use to more general nonlinear settings. A first attempt in the context of additive nonlinear systems was recently proposed by Hrustic et al. [20], where a linearly constrained EKF (LCEKF) was introduced, together with its use to mitigate parametric misspecifications on both system functions. Notice that this approach is fundamentally different from state constrained solutions [21,22], where LCs are imposed on the state and not on the filter. Overall, the use of non-stationary linear constraints for robust filtering in general nonlinear systems is still an open issue. In this contribution we further explore the use of LCs to mitigate the impact of model misspecifications in nonlinear dynamic systems with both non-additive noise and system inputs. (Notations: Italic, lower case boldface and upper case boldface indicate scalar, vector and matrix quantities; E[.] denotes the expectation operator; the filter estimates based on measurements up to time k are denoted by x k|k ; N (µ, R) is the normal distribution of mean µ and covariance R; m z and C z denote the mean and covariance of a given variable z.) The main contributions of this article are: • The derivation of a LCEKF for general non-additive nonlinear systems. • A demonstration of how to exploit LCs for robust filter design towards the mitigation of parametric modelling errors in both system functions.

•
The LCEKF performance improvement with respect to state-of-the-art EKF solutions is validated and discussed for an illustrative tracking and navigation problem.
Notice that the proposed methodology, even if not directly stated throughout the article, can also be used to mitigate a mismatch on both process and measurement noise means (constant or time-varying), which can be regarded as inputs.

Background
Consider a linear discrete state-space model, where the state vector x k ∈ C P k must be estimated from the available measurements y k ∈ C N k (for k ≥ 1), with F k−1 ∈ C P k ×P k−1 and H k ∈ C N k ×P k are known system model matrices, w k ∈ C P k and v k ∈ C N k are the process and measurement noise with zero mean and known covariance. If a minimum set of uncorrelation conditions holds [17], the recursive linear estimator of x k , which minimises the MSE (for k ≥ 1) is the KF (the superscript (·) b stands for the best solution in the MSE sense), where the optimal gain K b k ∈ C P k ×N k is the one that minimises the MSE, A possible way to robustify the KF is to incorporate LCs. How to use non-stationary constraints within the KF framework, leading to the so-called LCKF, has been recently proposed in [18]. In this case, the filter is computed from a "constrained" KF recursion [18], with L b k ∈ C P k ×N k . Such LCKF is fully adaptive and allows to incorporatenew LCs at every time k (∆ k ∈ C N k ×Q k and T k ∈ C P k ×Q k ). The use of such LCKF to robustify the filter under a mismatched model was also discussed in [18].
If we consider now an additive noise nonlinear discrete state-space model (NLDSSM), where f k−1 (·) and h k (·) are the known system model (process and measurement) functions, and E[dw k−1 ] = 0, E[dv k ] = 0; a standard approach to derive a nonlinear filter of x k is to assume that Equations (6) and (7) can be linearised at the vicinity of a nominal trajectory [3] yielding the standard EKF, where K b k is computed as in the unconstrained KF with F k−1 As has been recently proposed in [20], the corresponding LCEKF for additive noise systems is , where the constrained gain L b k is computed as in the original LCKF but using the linearised matrices, as done in the EKF. Notice that the LCKF in [18] is not applicable to nonlinear systems, and the LCEKF in [20] is only valid for additive systems. The goal of this article is to further analyse these methods and extend their use to more general settings.

An LCEKF with Non-Additive Noise and System Inputs, and Its Use in Robust Filtering
We consider now a more general NLDSSM, represented by the following state and measurement equations (k ≥ 1) with f k−1 (·) and h k (·) being known system model functions, and both noises w k and v k with known mean and covariance. First, we want to obtain the EKF-type linear MMSE filter of x k that can be written as in the previous cases. If tr(C x 0 ) 1, tr C w k−1 1, tr C x k−1 1 and tr C v k 1, we can resort to a first order Taylor expansion of f k−1 (x k−1 , w k−1 ) at the vicinity of m x k−1 , m w k−1 , and of h k (x k , v k ) at the vicinity of m x k , m v k . In such a case, if we restrict ourselves, for legibility, to the usual uncorrelation conditions, i.e., C x 0 ,w k = 0, C x 0 ,v k = 0, C w l ,w k = C w k δ l k , C v l ,v k = C v k δ l k , C w l ,v k = 0, the general form of the EKF becomes Notice that we do not explicitly exhibit the dependence on the system inputs in the previous EKF, which can be embedded in the noise means m w k−1 and m v k .

Mismatched and True System Models
As previously stated, it is unlikely that the practitioner has a full knowledge of the system, and therefore we want to cope with the situation where there is a true M T and a mismatched/assumed M A NLDSSM as follows, M T : Since the EKF of x k is based on the measurements and our knowledge of the model dynamics, any mismatch between the true model dynamics M T and the assumed one M A leads to a suboptimal filter, and possibly to a filter with bad performance, as the discrepancy between the two models increases. The existence of uncertainty on system nonlinear functions, either because of a parametric model error or a mismatch on system inputs, is taken into account as f k−1 (.) f k−1 (., !) and h k (.) h k (.,`), where ! and`are deterministic vector values or system inputs, and the possible parametric model mismatch is given by ! = ! + d ! and `=`+ d `. If errors d ! and d `are small, then the true state and measurement nonlinear functions differ from the assumed ones via first order Taylor series as follows, At time k ≥ 1, the EKF of x k is obtained from the Kalman-like recursion Equation (12) computed with the assumed model M A , with the gain L k obtained from the MSE minimisation, and also computed with M A .

Impact of System Parametric Modelling Errors
The first step towards the system model mismatch mitigation is to compute the estimation error induced by the use of the assumed model M A . Among the possible estimation error breakdowns, we look for the one that makes the terms in (17) appear, and which reduces to the error analysed in [18] when the system becomes linear (as shown in Section 3.5). An iterative approach was taken to integrate these two requirements and led to breaking down the error aŝ with the additional error term, Indeed, if we assume that x b k−1|k−1 is a good enough unbiased estimate of x k−1 , that is, then the first order approximations of Equations (21) and (22) arê where the linearised matrices for the assumed model are now

Mitigation of Parametric Modelling Errors through LCs
Once we have established the additional error induced by the model mismatch, the practitioner needs an efficient way to mitigate it. Notice that therefore, we can conclude that which defines a sensible set of constraints in order to mitigate at first order, which mainly consists of the bias, the error introduced by parametric modelling errors in the nonlinear system functions. However, expectations are probably not computable in several applications, and These LCs provide a non degenerate solution [18]

Exploiting LCs for Robust Filtering under Mismatch
Notice that by imposing LCs in Equation (30), the estimate obtained with the mismatched NLDSSM Equation (15) is matched to the true observation in Equation (16). Indeed, the LCKF minimises the MSE associated to the true state x k , matching the true observations to the assumed model. We detail the robust LCEKF methodology in the sequel. At every time k, Linear Constraints k|k−1 instead of at the unknown state means.

Special Cases: Linear and Additive Systems
If we assume that the NLDSSMs Equations (15) and (16) become linear, we have that, then, Equation (26) reduces to the linear system error analysed in [18], which confirms the relevance of Equations (25) and (26) as a first order approximation of Equations (21) and (22) and the choice of the form of Equation (21). If we assume the intermediate case where the NLDSSMs Equations (15) and (16) are still nonlinear but with additive noise, then, and m x k f k−1 m x k−1 , ! . Because the noise terms do not depend on the system parameters or inputs, we simply have that, Finally, if we consider nonlinear additive inputs of the form, and we recover the results in [20].

Illustrative Example: Robust Vehicle Tracking and Navigation
To illustrate the validity of the proposed LCEKF and its performance improvement with respect to standard (not robust) EKF solutions, we analyse a two vehicle tracking and navigation problem. We consider for both vehicles a tricycle dynamic model (i.e., two rear wheels provide the speed, and the front one controls the direction), also known as the Ackerman model [23], which is used to describe the kinematic behavior (2D position and orientation) of most parts of vehicles with three and four wheels (see Figure 1). In the sequel subscript (·) T refers to the tracking vehicle, and (·) B to the tracked one. Refer to Appendix A for details.
As shown in Figure 2, we consider a tracking vehicle at (x T , y T ) with linear velocity V T , orientation α and steering angle ψ, and a tracked vehicle at (x B , y B ) with linear velocity V B , orientation θ and steering angle φ. The state at discrete time k is

State Model with Both Additive and Non-Additive Noise
We consider a nonlinear state model of the form The additive noise w k−1 ∼ N (m w , C w ) models deviations on the vehicles' position and orientation, due to wind perturbations or unmodeled accelerations, with C w = diag(σ 2 x B , σ 2 y B , σ 2 θ , σ 2 x T , σ 2 y T , σ 2 α ). The non-additive noise w k−1 = [w v , w ψ ] T k−1 ∼ N (m w , C w ) models perturbations on the input parameters V T and ψ (i.e., the odometer may be affected by tire pressure and road conditions, and the steering information is subject to noise due to friction and wind), with C w = diag(σ 2 v , σ 2 ψ ), with V T = V T + w v and ψ = ψ + w ψ .

Nonlinear Measurement Model and Mismatch
The measurement principle is illustrated in Figure 2, where the tracker locates itself at point O and the tracked vehicle at point B in its coordinate frame (x'O'y'), and h k (x k ) is given by We consider the case where we may have an imperfect knowledge on the system inputs. In practice, the distance D B and D T are not accurate, i.e., there exists a deviation between the true values, D B and D T , and the assumed ones, D B and D T , dD B = D B − D B and dD T = D T − D T . In order to mitigate the possible impact of such error, we resort to the LCEKF in Section 3.4, incorporating the constraint

Scenarios
We consider the following setup: dt = 1 ms, C

Results
Numerical results for the three scenarios are summarised in Figure 3, where we show the position MSE vs time for the five filters obtained over 100 Monte Carlo runs. First, notice that in each scenario (S1, S2 and S3) and regardless of the mismatch value dD the standard EKF always diverges and is not able to correctly estimate the states of interest. That is, even a minor mismatch on the distance from the front wheel to the rear axle induces a performance breakdown on the standard misspecified EKF; therefore, such a mismatch must be accounted for. In addition, the larger the mismatch the faster the performance breakdown. Second, the LCEKF DB accounting only for a mismatched input on the tracked vehicle (D B ) does not improve the standard EKF performance (results are superimposed in the plots). Third, the LCEKF DT , which mitigates the tracker input mismatch (D T ), slightly improves the performance w.r.t. the EKF and LCEKF dB but does not prevent a performance breakdown. Then, we can conclude that with a mismatch on both vehicles' distance D, and regardless of whether the system is affected by either additive or non-additive noise, either we correctly deal with both mismatched inputs or the filters are not useful anymore.
In contrast to the previous filters, the LCEKF DB,DT , which properly accounts for both input mismatched values, is able to correctly estimate the states and avoid the performance breakdown exhibited by the mismatched filters, regardless of the mismatched value dD, and therefore it is a powerful solution to coping with a possible system model mismatch. Notice that the LCEKF DB,DT performance is almost superimposed with the optimal EKF re f in the plots, which again proves the good behaviour of the LCEKF for both additive and non-additive noises. It is worth pointing out that the non-additive noise has a larger impact than the additive one on the system dynamics, which can be seen from the optimal MSE results (i.e., MSE EKF re f ≈ 10 −4 for the additive noise case S1 and MSE EKF re f ≈ 0.5 × 10 −2 for both non-additive noise cases S2 and S3); this drives the estimation error (MSE EKF re f S2 ≈ MSE EKF re f S3). This further supports the need for filters able to cope with non-additive noise and system inputs, and the interest of the robust LCEKF formulation proposed in this contribution.

Conclusions
In this contribution we explored the use of linear constraints to design robust nonlinear KF-type filtering strategies. This was shown in previous contributions to be a promising solution to mitigate parametric modelling errors in both system functions, a problem of interest in real-life applications. In that perspective, we derived a linearly constrained EKF for systems affected by non-additive noise and system inputs, which generalises previous solutions only able to cope with either linear systems or additive noises. Within this framework, it was shown how to exploit linear constraints to mitigate possible nonlinear system model parametric misspecifications on both process and measurement functions. This approach encompasses both errors on system function parameters and system inputs, as well as on noise mean values. A two vehicle robust tracking and navigation problem was used to show the validity and performance improvement provided by the robust LCEKF with respect to state-of-the-art solutions. In addition, it was shown that the non-additive noise has a larger impact than the additive one on the system dynamics, and drives the estimation error, which further supports the need for filters able to cope with non-additive noise and system inputs, and the interest of the robust LCEKF formulation.
As future work, it would be of interest to analyse the use of linear constraints within the context of sigma-point filtering, to better cope with the system nonlinearity and avoid the EKF linearisation errors. Moreover, the proposed solution should be compared to alternative system identification techniques, which, rather than mitigating the impact of mismatched parameters, estimate them together with the states of the system. Another research line would be the generalisation to more general mismatched systems, for instance, considering a possible mismatch on covariance matrices, as well as the comparison or combination with noise statistics' parameter estimation techniques.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Details on the Robust Vehicle Tracking and Navigation Example
Appendix A. 1

. Tricycle Dynamic Model
The tricycle model (also known as Ackerman model) can be used to model the kinematic behavior of most part of vehicles with three and four wheels. As it is shown in Figure 1, it consists of two main wheels which provide the vehicle speed, and a third wheel which controls the vehicle direction. The turn radius ρ is given by where D denotes the distance from the front wheel to the rear axle. We consider the following equations where, and then, Considering the state at discrete time k, x T k = [(x B ) k , (y B ) k , θ k , (x T ) k , (y T ) k , α k ] (in the global frame (xOy)), and the previous equations, we can obtain the state equation for the general case with additive and non-additive noise terms, x k = f k−1 (x k−1 , u, w k−1 ) + w k−1 , detailed in Equations (41)-(47). For the three scenarios of interest, we give in the sequel the corresponding approximations (i.e., we omit the dependency with discrete time k for the sake of simplicity): • (S1) Only additive noise (w k−1 = 0, N k−1 = 0): Let then • (S2) Only non-additive noise (w k−1 = 0): Let then where ∂f 4 (x, u, w) ∂w v = dt sinc(q T (x))cos(g T (x)) − dt q T (x)sinc(q T (x))sin(g T (x)) + dt (cos(q T (x)) − sinc(q T (x)))cos(g T (x)), )sin(g T (x)) + dt q T (x) sinc(q T (x))cos(g T (x)) + dt (cos(q T (x)) − sinc(q T (x)))sin(g T (x)), ∂f 5 (x, u, w) ∂w ψ = (V T + w v ) 2 dt 2 2D T 1 + tan 2 (ψ + w ψ ) sinc(q T (x))cos(g T (x)) + (V T + w v ) dt 1 + tan 2 (ψ + w ψ ) tan(ψ) sin(g T (x))cos(q T (x)) − (V T + w v ) dt 1 + tan 2 (ψ + w ψ ) tan(ψ + w ψ ) sin(g T (x))sinc(q T (x)), • (S3) Both additive and non-additive noise: in this case F k−1 and N k−1 are computed as in S2)

Appendix A.3. Nonlinear Measurement Model
The nonlinear measurement model is given in Equation (48). To use the EKF and LCEKF we need to compute the approximation H k = ∂h k (x k ,v k , `) which is given by Notice that in this case M k = 0.

Appendix A.4. Linear Constraints
As already stated, the goal is to mitigate the possible mismatch impact by exploiting the constraints I − L k H k dF k−1 = 0, dF k = ∂f k (x k , u, w k ) ∂a T with u = [V B , φ, D B , V T , ψ, D T ] T and a = [D B , D T ] T . Omitting the dependency with discrete time k for the sake of simplicity, such constraint is computed as follows. Let q(x) = V B dt tan(φ)/(2D B ), g(x) = θ + q(x), q T (x) = (V T + w v ) dt tan(ψ + w ψ )/(2D T ), g T (x) = α + q T (x), sinc(q(x))cos(g(x)) − V B dt D B cos(q(x))cos(g(x)) + V B dt D B q(x)sinc(q(x))sin(g(x)), sinc(q(x))cos(g(x)) − V B dt D B cos(q(x))sin(g(x)) − V B dt D B q(x)sinc(q(x))cos(g(x)), ∂f 4 (x, u) ∂D T = V T dt D T sinc(q T (x))cos(g T (x)) − V T dt D T cos(q T (x))cos(g T (x)) + V T dt D T q T (x)sinc(q T (x))sin(g T (x)), ∂f 5 (x, u) ∂D T = − V T dt D T sinc(q T (x))sin(g T (x)) − V T dt D T cos(q T (x))sin(g T (x)) − V T dt D T q T (x)sinc(q T (x))cos(g T (x)),