next up previous contents
Next: Calculating the a priori Up: Filtering Stochastic Control Systems Previous: and are not

Shadows in Infinite Dimensional Space

I imagine the systems's behavior in time as an unfolding in the dark recesses of tex2html_wrap_inline782 . Each time step brings a new element to the basis, a new axis. The State Sequence goes twisting and writhing along, while we, in some hyperdimensional Allegory of the Cave, peer at it's flickering images in the space of our observations. While the Probability Space spanned by the system builds over time into some baroque geometry of causality, we poke and peer at its shadow, trying to glean some hint of the Truth which underlies surface appearances, rife as they are with noise and error.

We begin this process by performing a Gramm-Schmidt orthogonalization of both our State and Observation sequences. This is a recursive decomposition of the Random Variables. By Theorem 1 we may decompose tex2html_wrap_inline915 , the tex2html_wrap_inline977 observation, into its projections on tex2html_wrap_inline979 and its complement tex2html_wrap_inline979 for any n.

eqnarray497

for notational convenience, and more or less in keeping with Kalman's original notation, tex2html_wrap_inline985

eqnarray507

To generate an orthogonal basis for tex2html_wrap_inline969 we merely note that if we assume some knowledge of initial conditions we may set

eqnarray517

tex2html_wrap_inline989 are pairwise independent, and thus tex2html_wrap_inline991 . We will thus use the tex2html_wrap_inline993 's for our ``reference basis''. Our next step is to perform a similar decomposition on the state variable sequence.

tex2html_wrap_inline995

eqnarray524

Since tex2html_wrap_inline997 we call tex2html_wrap_inline999 the estimate of tex2html_wrap_inline913 based on the (orthogonalized) observations tex2html_wrap_inline1003 , and tex2html_wrap_inline1005 the error.

Now, by orthogonality of the tex2html_wrap_inline1007 we have

equation186

so we may decompose the projections:

eqnarray545

That is, we may split the process of calculating this turn's estimate into a sum of an estimate based on last turn's estimate and one based entirely on the current observation. This has obvious computational as well as aesthetic merit, and the fact that the two spaces are orthogonal complements allows us to treat them seperately; the two estimates are independent random variables gif . If we consider tex2html_wrap_inline915 the ``current'' observation, we may term the projection of tex2html_wrap_inline913 on tex2html_wrap_inline1013 the a priori estimate, that on tex2html_wrap_inline1015 the update, and on tex2html_wrap_inline1017 the posteriori estimate, because it is our best estimate based on all the information availible at that time.




next up previous contents
Next: Calculating the a priori Up: Filtering Stochastic Control Systems Previous: and are not

Scot Free Kennedy
Sat Sep 13 00:27:51 PDT 1997