14.02.2013 Views

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 4 / Target Tracking<br />

_ _<br />

The innovation (∆Z) is the difference between the actual and predicted<br />

measurements,<br />

∆<br />

− Zˆ Z ~<br />

⎛<br />

⎜<br />

⎝<br />

∂ h<br />

⋅ ∆X<br />

+ 2<br />

∂ X<br />

T<br />

Z : = : =<br />

⋅ ∆X<br />

2<br />

4-12<br />

−1<br />

2<br />

∂ h ⎞<br />

⋅ ⋅ ∆X<br />

+ L ⎟<br />

+ υ<br />

∂ X ⎠<br />

Equation 4.4-15<br />

The measurement Jacobian (H), and Hessian [G], are computed using the<br />

state at time (tk+1) - in the EKF the Hessian is ignored. The measurement<br />

innovation is then apportioned between the states by the Kalman gains so as<br />

to minimise the trace of the covariance matrix, Gelb [G.12] (p186),<br />

K<br />

k+<br />

1<br />

: =<br />

X ˆ<br />

H<br />

+<br />

k+<br />

1<br />

k+<br />

1<br />

: =<br />

C<br />

⋅ C<br />

−<br />

k+<br />

1<br />

−<br />

k+<br />

1<br />

X ˆ<br />

−<br />

k+<br />

1<br />

⋅ H<br />

⋅ H<br />

T<br />

k+<br />

1<br />

T<br />

k+<br />

1<br />

+ K<br />

+ R<br />

k+<br />

1<br />

k+<br />

1<br />

⋅ ∆Z<br />

k+<br />

1<br />

≡<br />

C<br />

−<br />

k+<br />

1<br />

⋅ H<br />

A<br />

The covariance matrix is update using the simple expression,<br />

C<br />

+<br />

k+<br />

1<br />

−<br />

( I − K ⋅ H ) ⋅ C<br />

= 3 k+<br />

1 k+<br />

1 k+<br />

1<br />

Equation 4.4-16<br />

k+<br />

1<br />

T<br />

k+<br />

1<br />

Equation 4.4-17<br />

Equation 4.4-18<br />

Nishimara [N.3] showed that Joseph's form of this equation yields the correct<br />

covariance whether the Kalman gains are optimal or not. This form should<br />

always be used in the IEKF, Kerr [K.4] ,<br />

−<br />

T<br />

T<br />

( I − K ⋅ H ) ⋅ C ⋅ ( I − K ⋅ H ) + K ⋅ R ⋅ K<br />

+<br />

C k+<br />

1 : =<br />

k+<br />

1 k+<br />

1 k+<br />

1<br />

k+<br />

1 k+<br />

1 k+<br />

1 k+<br />

1 k+<br />

1<br />

Equation 4.4-19<br />

Measurements should be processed serially to avoid matrix inversion, and to<br />

accommodate reduced order and asynchronous measurement sets. The<br />

effect of process model and linearisation errors on [A], [H], [Φ] and (∆Z) is<br />

quantified by Hanlon [H.8] .<br />

4.4.2 Iterated EKF<br />

EKF measurement linearisation can be performed about any trajectory. If<br />

the filter is to be well conditioned the trajectory about which linearisation is<br />

performed should ideally be close to the reference trajectory. If the state<br />

estimates become too inaccurate, as might be the case after a long period of<br />

prediction, transformed measurements are often substituted. The IEKF<br />

iterates process and measurement linearisation, a process that if taken to

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!