NIRK-based Cholesky-factorized square-root accurate continuous-discrete unscented Kalman filters for state estimation in nonlinear continuous-time stochastic models with discrete measurements

被引:22
作者
Kulikov, G. Yu. [1 ]
Kulikova, M. V. [1 ]
机构
[1] Univ Lisbon, Inst Super Tecn, CEMAT, Ave Rovisco Pais, P-1049001 Lisbon, Portugal
关键词
Continuous-discrete stochastic system; Square-root unscented Kalman filter; Sigma point differential equations; Gauss-type nested implicit Runge-Kutta formula with local and global error controls; Stochastic radar tracking model; GLOBAL ERROR ESTIMATION; ORDINARY DIFFERENTIAL-EQUATIONS; KUTTA COLLOCATION METHODS; 2-STEP W-METHODS; TOLERANCE PROPORTIONALITY; EFFICIENT FAMILY; STIFF SYSTEMS; PEER METHODS; CUBATURE; COVARIANCES;
D O I
10.1016/j.apnum.2019.08.021
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
This paper further advances the idea of accurate Gaussian filtering towards efficient unscented-type Kalman methods for estimating continuous-time nonlinear stochastic systems with discrete measurements. It implies that the differential equations evolving sigma points utilized in computations of the predicted mean and covariance in time-propagations of the Gaussian distribution are solved accurately, i.e. with negligible error. The latter allows the total error of the unscented Kalman filtering technique to be reduced significantly and gives rise to the novel accurate continuous-discrete unscented Kalman filtering algorithm. At the same time, this algorithm is rather vulnerable to round-off and numerical integration errors committed in each state estimation run because of the need for the Cholesky decomposition of covariance matrices involved. Such a factorization will always fail when the covariance's positivity is lost. This positivity lost issue is commonly resolved with square-root filtering implementations, which propagate not the full covariance matrix but its square root (Cholesky factor) instead. Unfortunately, negative weights encountered in applications of the accurate continuous-discrete unscented Kalman filter to high-dimensional stochastic systems preclude from designing conventional square-root methods. In this paper, we address this problem with low-rank Cholesky factor update procedures or with hyperbolic Q R transforms used for yielding J-orthogonal square roots. Our novel square-root algorithms are justified theoretically and examined numerically in an air traffic control scenario. (C) 2019 IMACS. Published by Elsevier B.V. All rights reserved.
引用
收藏
页码:196 / 221
页数:26
相关论文
共 102 条
[1]  
Aggoun L., 2005, MEASURE THEORY FILTE
[2]   Numerical investigations on global error estimation for ordinary differential equations [J].
Aid, R ;
Levacher, L .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 1997, 82 (1-2) :21-39
[3]   A SQUARE ROOT FORMULATION OF KALMAN COVARIANCE EQUATIONS [J].
ANDREWS, A .
AIAA JOURNAL, 1968, 6 (06) :1165-&
[4]  
[Anonymous], P 2018 22 INT C SYST
[5]  
[Anonymous], VARIABLE STEPS UNPUB
[6]  
[Anonymous], 2013, Model predictive control: Theory and design
[7]  
[Anonymous], 1994, ANN NUMER MATH
[8]  
[Anonymous], 2008, Numerical Methods for Ordinary Differential Equations
[9]  
[Anonymous], P 35 INT C DEC CONTR
[10]  
[Anonymous], GENERAL LINEAR METHO