Levenberg-Marquardt Method Based Iterative Square Root Cubature Kalman Filter and its Applications to Maneuvering Re-entry Target Tracking

Publications

Share / Export Citation / Email / Print / Text size:

International Journal of Advanced Network, Monitoring and Controls

Xi'an Technological University

Subject: Computer Science , Software Engineering

GET ALERTS

eISSN: 2470-8038

DESCRIPTION

8
Reader(s)
29
Visit(s)
0
Comment(s)
0
Share(s)

SEARCH WITHIN CONTENT

FIND ARTICLE

Volume / Issue / page

Related articles

VOLUME 3 , ISSUE 2 (May 2018) > List of articles

Levenberg-Marquardt Method Based Iterative Square Root Cubature Kalman Filter and its Applications to Maneuvering Re-entry Target Tracking

Mu Jing * / Wang Changyuan *

Keywords : Nonlinear Filtering, Square Root Cubature Kalman Filter, Levenberg-Marquardt Method, Maneuvering Reentry Targets Tracking

Citation Information : International Journal of Advanced Network, Monitoring and Controls. Volume 3, Issue 2, Pages 10-13, DOI: https://doi.org/10.21307/ijanmc-2018-027

License : (CC-BY-NC-ND-4.0)

Published Online: 07-May-2018

ARTICLE

ABSTRACT

Levenberg-Marquardt (abbr.L-M) method based iterative square root cubature Kalman filter (abbr. ISRCKFLM) inherits the numerical stability of square root Cubature Kalman filter and effectively suppresses the influence of the larger initial estimation error and the nonlinearity of the measurement equation on the state estimation in the nonlinear state estimation due to obtaining the optimal state and variance estimates using the latest measurement through L-M method. We apply the ISRCKFLM algorithm to the state estimation of maneuvering re-entry target tracking, the simulation results demonstrate that the ISRCKFLM algorithm has better accuracy of state estimation, comparable to Unscented Kalman filter and square root Cubature Kalman filter, according to estimation error analysis of the position, velocity, drag coefficient, turn coefficient and climbing force coefficient, and has fast convergence rate.

Graphical ABSTRACT

I. INTRODUCTION

The state estimate of maneuver re-entry target is highly non-linear filtering problem and has been much attention in the academic and engineering domain [1, 2]. Up to now, the commonly used non-linear filtering is the extended Kalman filter (EKF) [3]. The EKF is based on first-order Taylor approximations of state transition and observation equation about the estimated state trajectory under Gaussian assumption, so EKF may introduce significant bias, or even convergence problems due to the overly crude approximation [4]. Especially, for highly nonlinear problems such as state estimation of maneuver re-entry target, EKF may produce large filtering errors and even divergence. Based on the Unscented Transformation (UT), the Unscented Kalman filter (UKF) use a set of deterministic sampling points to approximate the posterior probability distribution of the system state, that is, the Sigma points are used to capture the mean and variance information of the state [5].

Recently, as a new way to solve the nonlinear estimation problem, cubature rules based cubature Kalman filter (CKF) proposed in [6] uses numerical multi-dimensional integral to approximate the recursive Bayesian estimation integrals under the Gaussian assumption. The CKF can solve high-dimensional nonlinear filtering problems with minimal computational effort. Then the square root cubature Kalman filter algorithm (SRCKF) was proposed in order to improve the numerical stability [6].

On the other hand, in order to decrease the effect of initial estimation error and nonlinearity of measurement equation, Levenberg-Marquardt method based iterative square root cubature Kalman filter (ISRCKFLM) was developed on the basis of the SRCKF in [7]. In this paper, we apply the ISRCKFLM algorithm to state estimation of maneuver re-entry target. Simulations demonstrate that the ISRCKFLM algorithm can greatly improve the tracking accuracy of maneuver re-entry target and obtain fast convergence, compared with the UKF and SRCKF algorithms.

The rest of the paper is organized as follows. We begin with a description of the ISRCKFLM algorithm in Section 2. Then we apply the ISRCKFLM algorithm to track re-entry ballistic target (RBT) with unknown ballistic coefficient and discuss the simulation results in Section 3. Finally, we draw conclusion in Section 4.

II. L-M BASED ITERATIVE SQUARE ROOT CUBATURE KALMAN FILTER

Consider the following non-linear dynamics system:

(1)
xk=f(xk1)+wk1,
(2)
zk=h(xk)+vk.

where f and h are some known nonlinear functions; xk ∈ ℝnx and zk ∈ ℝnz is state and the measurement vector, respectively; wk−1 and vk are process and measurement Gaussian noise sequences with zero means and covariance Qk−1 and Rk, respectively, and {wk−1} and {vk} are mutually uncorrelated.

Suppose that the state distribution at k-1 time is xk1(x^k1,Sk1Sk1T), Levenberg-Marquardt based Iterative square root cubature Kalman filter (ISRCKFLM) is described as follows.

A. Time Update

  • 1) Calculate the cubature points and propagate the cubature points through the state equation

    (3)
    Xi,k1=Sk1ξi+x^k1.
    (4)
    Xi,k*=f(Xi,k1).

    where ξi=m/2[1]i, ωi = 1/m, i = 1, … m = 2nx, the [1] xi is a dimensional vector and is generated according to the way described in [2].

  • 2) Evaluate the predicted state and square root of the predicted covariance

    (5)
    x¯k=i=1mωiXi,k*.
    (6)
    S¯k=Tria([χk*SQ,k1]).

    here, SQ, k−1 denotes a square-root factor of Qk −1 and Tria() is denoted as a general triagularization algorithm. The matrix χk* is defined as:

    (7)
    χk*=1m[X1,k*x¯k X2,k*x¯k,,Xm,k*x¯k].

  • 3) Evaluate the modified covariance:

    (8)
    P˜k=[IS¯kS¯kT(S¯kS¯kT+1μiI)1]S¯kS¯kT.

    where μi is adjusting parameter.

B. Measurement update

  • 1) Set the initial value as: x^k(0)=x¯k.

  • 2) Assuming the i-th iterate x^k(i), calculate the matrix

    (9)
    Lk(i)=P˜kJhT(x^k(i))[Jh(x^k(i))P˜kJhT(x^k(i))+Rk]1.

  • 3) Calculate the i-th iterate

    (10)
    x^k(i+1)=x¯k+Lk(i){zkh(xk(i))J(x^k(i))(x¯kx^k(i))}.μi{ILk(i)J(x^k(i))}P˜k(x¯kx^k(i))

  • 4) Calculate the iteration termination condition

    (11)
    x^k(i+1)x^k(i)ε or i=Nmax.

    ε and Nmax are predetermined threshold and maximum iterate number, respectively. If the termination condition meets, the iterate return to 5); otherwise, set x^k(i)=x^k(i+1), continue to 2).

  • 5) Calculate the state estimation at time instant

    (12)
    x^k=x^k(N).

  • 6) Evaluate the cross-covariance and square root of innovation covariance at k time

    (13)
    Pxz=S¯kS¯kTJhT(x^k(N)).
    (14)
    Szz=Chol([Jh(x^k(N))S¯k SR,k]).

  • 7) Calculate the square root of covariance at k time

    (15)
    Kk=Pxz/SzzT/Szz.
    (16)
    Sk=Chol([S¯kKkJh(x^k(N))S¯k KkSR,k]).

    where symbol “/” represents the matrix right divide operator.

III. APPLICATIONS TO MANEUVERING REENTRY TARGET TRACKING

In the simulation, the trajectory of the maneuvering reentry generated in [8] is used with the same parameters, initial state and covariance estimation. To compare the performance of the ISRCKFLM with the UKF and SRCKF algorithms, we obtain the root mean square errors (RMSEs) in the position, velocity, drag coefficient, turn coefficient and climb coefficient showed in the Figure 1-5. All performance curves were obtained by averaging over 100 independent Monte Carlo runs. Table. 1 lists the accumulated root mean square errors (ARMSEs).

Figure 1.

RMSE in Position

10.21307_ijanmc-2018-027-f001.jpg
Figure 2.

RMSE in Velocity

10.21307_ijanmc-2018-027-f002.jpg
Figure 3.

RMSE in Drag Coefficient

10.21307_ijanmc-2018-027-f003.jpg
Figure 4.

RMSE in Turn Coefficient

10.21307_ijanmc-2018-027-f004.jpg
Figure 5.

RMSE in Climbing Coefficient

10.21307_ijanmc-2018-027-f005.jpg
TABLE I.

ARMSES IN THREE FILTERS

10.21307_ijanmc-2018-027-tbl1.jpg

From Figur. 1-2, we can see the ISRCKFLM’s RMSEs in positon and velocity are lower than those of UKF and SRCKF algorithms, especially, the RMSEs of ISRCKFLM algorithm are effectively suppressed, and those of UKF and SRCKF algorithms have a big jump when the target experiences the high maneuver at altitude 23km and 10km.

As for the estimate on drag coefficient, turn coefficient and climbing coefficient of maneuver re-entry target, from Figure. 3-5 we can see the RMSEs of three algorithms decrease between 35km-30km at altitude, and the RMSEs of three algorithms reach the smallest at about 30km at altitude. Then the RMSEs of three algorithms increases when a dive and turn maneuver occurs on re-entry target, however, the ISRCKFLM’s RMSEs are lower than those of the other two algorithms. The RMSEs of all filters gradually increase when the target experiences an increased climb force at altitude 16km-10km, but the ISRCKFLM’s RMSE increases slightly slower. After the maneuver was withdrawn, the RMSE of the three algorithms began to decline.

Moreover, the ARMSEs of the ISRCKFLM in the positon, velocity, drag coefficient, turn coefficient and climbing coefficient are lower than those of UKF and SRCKF algorithms from Table. 1.

Therefore, on the basis of the simulation results presented in Figure. 1-5, we can draw a conclusion that the ISRCKFLM yields on the superior performance over the UKF and SRCKF on the state estimation of maneuvering re-entry target.

IV. CONCLUSION

In this paper we apply the ISRCKFLM algorithm to maneuvering re-entry target tracking. The latest measurement are fully used in the ISRCKFLM algorithm, and innovation covariance and cross-covariance are improved in the iterative process, so we can obtain the optimal state estimation and covariance estimation. Simulation results demonstrate that the performance of ISRCKFLM algorithm is superior to UKF and CKF algorithms by analysis of errors in position, velocity, drag coefficient, turn coefficient and climbing coefficient, and has the faster convergence rate.

ACKNOWLEDGMENT

The authors would like to thank the support of fund of the State and Local Joint Laboratory of Advanced Network and Monitoring Control Engineering (No. GSYSJ201603); Natural Science Basic Research Program of Department of science and technology of Shaanxi Province - Youth Foundation (No. 2017JQ6056); The Key Project of Department of Education of Shaanxi Province (No. 16JF012).

References


  1. Li XR, Jilkov VP. A survey of maneuvering target tracking - part II: ballistic target models. Signal and Data Processing of Small Targets, 2001, 4473:559-581.
  2. Farina A, Ristic B, Benvenuti D. Tracking a ballistic target: comparison of several nonlinear filters. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(3):854-867.
    [CROSSREF]
  3. Bar-Shalom Y, Li XR, Kirubarajan T. Estimation with Applications to Tracking and Navigation. New York: John Wiley &Son, Inc, 2001.
    [CROSSREF]
  4. Athans M, Wishner RP, Bertolini A. Suboptimal state estimation for continuous-time nonlinear system from discrete noisy measurements[J]. IEEE Transactions on Automatic Control, 1968, 13(5):504-514.
    [CROSSREF]
  5. Julier SJ, Uhlmann JK. Unscented filtering and nonlinear estimation. Proceedings of the IEEE, 2004, 92(12):1958-1958.
    [CROSSREF]
  6. Arasaratnam I, Haykin S. Cubature Kalman filters. IEEE Transactions on Automatic Control, 2009, 54(6):1254-1269.
    [CROSSREF]
  7. Mu. J.; Cai. Y.; Wang. C.Y. L-M Method Based Iteration Cubature Kalman Filter and its Applications, Journal of Xi’an Technological University, 2013, 33(1):1-6.
  8. Jing Mu, Yuanli Cai, Changyuan Wang. State estimation of Maneuvering Reentry Target Using Likelihood Based Iterated Divided Difference Filter, Journal of Xi’an Technological University, 32(5): 349-354, 2012.5.
XML PDF Share

FIGURES & TABLES

Figure 1.

RMSE in Position

Full Size   |   Slide (.pptx)

Figure 2.

RMSE in Velocity

Full Size   |   Slide (.pptx)

Figure 3.

RMSE in Drag Coefficient

Full Size   |   Slide (.pptx)

Figure 4.

RMSE in Turn Coefficient

Full Size   |   Slide (.pptx)

Figure 5.

RMSE in Climbing Coefficient

Full Size   |   Slide (.pptx)

TABLE I.

ARMSES IN THREE FILTERS

Full Size   |   Slide (.pptx)

REFERENCES

  1. Li XR, Jilkov VP. A survey of maneuvering target tracking - part II: ballistic target models. Signal and Data Processing of Small Targets, 2001, 4473:559-581.
  2. Farina A, Ristic B, Benvenuti D. Tracking a ballistic target: comparison of several nonlinear filters. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(3):854-867.
    [CROSSREF]
  3. Bar-Shalom Y, Li XR, Kirubarajan T. Estimation with Applications to Tracking and Navigation. New York: John Wiley &Son, Inc, 2001.
    [CROSSREF]
  4. Athans M, Wishner RP, Bertolini A. Suboptimal state estimation for continuous-time nonlinear system from discrete noisy measurements[J]. IEEE Transactions on Automatic Control, 1968, 13(5):504-514.
    [CROSSREF]
  5. Julier SJ, Uhlmann JK. Unscented filtering and nonlinear estimation. Proceedings of the IEEE, 2004, 92(12):1958-1958.
    [CROSSREF]
  6. Arasaratnam I, Haykin S. Cubature Kalman filters. IEEE Transactions on Automatic Control, 2009, 54(6):1254-1269.
    [CROSSREF]
  7. Mu. J.; Cai. Y.; Wang. C.Y. L-M Method Based Iteration Cubature Kalman Filter and its Applications, Journal of Xi’an Technological University, 2013, 33(1):1-6.
  8. Jing Mu, Yuanli Cai, Changyuan Wang. State estimation of Maneuvering Reentry Target Using Likelihood Based Iterated Divided Difference Filter, Journal of Xi’an Technological University, 32(5): 349-354, 2012.5.

EXTRA FILES

COMMENTS