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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 5 / Missile State Observer<br />

_ _<br />

Propagating the state covariance matrix over the time interval (∆t),<br />

−<br />

k+<br />

1<br />

k+<br />

1<br />

k<br />

4-11<br />

−<br />

k<br />

k+<br />

1 T<br />

T<br />

[ Φ ] + Γ ⋅ Q ⋅<br />

C : = Φ ⋅ C ⋅<br />

Γ<br />

k<br />

k<br />

k<br />

k<br />

Equation 4.4-8<br />

The symmetric, semi-definite, covariance matrix (Q) is computed between<br />

measurement updates from,<br />

t<br />

∫ + k 1<br />

t<br />

k<br />

Γ<br />

⋅ Q<br />

⋅ Γ<br />

( τ ) ⋅ diag q ( τ )<br />

k<br />

: =<br />

( ) ⋅ G ( τ )<br />

T k+<br />

1 T<br />

[ ] ⋅ [ ] ⎟<br />

⎞ dτ<br />

⎜<br />

⎛ k+<br />

1<br />

Φτ<br />

⋅ G<br />

Φ<br />

⎝<br />

k<br />

Equation 4.4-9<br />

If ∆t is small, the transition matrix collapses to the identity matrix and,<br />

Q<br />

k<br />

: =<br />

G<br />

k<br />

⋅ diag<br />

T<br />

k<br />

T<br />

( q ( k ) ) ⋅ G ⋅ ∆t<br />

k<br />

τ<br />

⎠<br />

Equation 4.4-10<br />

Assuming that the sensor measurements ( Z ~ ) are non-linear functions of the<br />

reference state vector at time (tk+1) plus zero mean Gaussian noise,<br />

E<br />

Z ~<br />

: = h X<br />

( ) + υ<br />

T<br />

( υ ) : = 0 ⇒ E ( υ ⋅ υ ) : = R<br />

Equation 4.4-11<br />

Equation 4.4-12<br />

[R] is the symmetric, positive definite measurement covariance matrix. The<br />

expected measurement is computed using the state approximation<br />

propagated to the measurement reference time (tk+1),<br />

Z ˆ<br />

: =<br />

h<br />

⎜<br />

⎛<br />

⎝<br />

−<br />

Xk+1 ˆ<br />

⎟<br />

⎞<br />

⎠<br />

Equation 4.4-13<br />

Expanding the measurement functions in a series about the estimated state,<br />

h<br />

( ) ⎛ − ⎞<br />

−1<br />

T<br />

= X + H ⋅ ∆X<br />

+ 2 ⋅ ∆X<br />

⋅ G ⋅ ∆X<br />

+ L<br />

ˆ<br />

X : h<br />

k+<br />

1<br />

⎜<br />

⎝<br />

k+<br />

1<br />

⎟<br />

⎠<br />

Equation 4.4-14

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

Saved successfully!

Ooh no, something went wrong!