next up previous contents
Next: Filtering Stochastic Control Systems Up: Spaces of Random Variables Previous: Random Variables

Expectations and Products

Like any other vector space, tex2html_wrap_inline782 has a dual space of functionals (ie, mappings from a vector space to its set of scalars) gif. In fact, since it is an infinite-dimensional Hilbert Space, it is self dual. One very important functional is the expectation. This can be thought of as the center of mass (or moment) of a function on the state space. That is,

displaymath798

Where P(x) is the measure used in the definition of x as a Random variable. In the usual terminology P(x) is called the Probability Density Function, and tex2html_wrap_inline826 .

We can define an inner product in tex2html_wrap_inline782 by taking the expected value of the standard inner product in tex2html_wrap_inline715 ,

displaymath799

Note that this is called the Covariance by statisticians,and that the Norm in tex2html_wrap_inline782 induced by this inner product is simply the standard deviation.

displaymath800

We will generally just work with the variance, tex2html_wrap_inline834 . Also of great interest is the induced distance function,

displaymath801

though again, we generally encounter tex2html_wrap_inline836 in this context. This is because, if we consider x to be in some sense an estimate of y, then a natural measure of error is the distance function. We use its square for convenience, since it works just as well. That is, minimizing tex2html_wrap_inline836 gif (called the ``Sum Of Squares'' error) minimizes d(x,y) as well.

Another construction of great utility here is the Outer Product, which is a function of two random variables yielding their covariance matrix. This is defined as

displaymath802

Note that this is simply the covariance matrix of x with y, though we often also write, for notational convenience,

displaymath803

which is, of course, the variance matrix as it is defined in standard statistics texts gif.

It is worth making explicit some properties of this operation which follow directly from the definition. Given tex2html_wrap_inline848 and A,B matrices of appropriate dimension to pre-multiply them with X and Y, respectively:

eqnarray423

These properties gif will be used again and again in what follows.

In fact, it is this outer product we will use here to define orthogonality. We call two vectors tex2html_wrap_inline856 orthogonal if [x,y]=0

Why are we revising our definition? Because it covers more general situations, and frees us from worrying unneccessarily about dimensionality. We really aren't changing anything, just generalizing. In the case where the two vectors corresspond to vectors of equal dimension in the state space, we note that the Trace of the Outer Product equals the inner product,

displaymath804

So that if the Outer Product is the Zero Matrix, the Inner Product is the Zero Scalar. So the definition coincides exactly with the usual definition of orthogonality in terms of the inner product.

If we cannot directly calculate the inner product gif due to seemingly conflicting dimensionality, this more general definition still makes sense. Consider our original construction of tex2html_wrap_inline782 as the Cartesian Product of tex2html_wrap_inline874 , or tex2html_wrap_inline737 . The tex2html_wrap_inline878 element of [x,y] gives the (scalar) covariance between the tex2html_wrap_inline882 component of x and the tex2html_wrap_inline886 of y. So if [x,y] = 0 it means that every component is probabalistically orthogonal to each component of the other vector. That is, looking at everything as happenning in tex2html_wrap_inline790 , with one-dimensional random variables, the vectors are component-wise orthogonal. Note that independence implies orthogonality. If both RV's are Gaussian, the converse can be shown to be true as well.

exercise Characterize the unit hypersphere in tex2html_wrap_inline782 , tex2html_wrap_inline896 .

displaymath803

That is, the set of Random Variables with Covariance one, or Trace([x])=1.

Of slightly more interest is a Sphere around a specific point in tex2html_wrap_inline782 ,

eqnarray152

So we have some (infinite dimensional) manifold in this wild space of random variables. What dynamics can we imagine sliding across the surface of this bubble of causality, shimmering in the piercing light of reason? What atlas maps its surface?

The orthogonal projection of tex2html_wrap_inline898 on tex2html_wrap_inline900 (assuming, for simplicity, both are zero-mean) is the random variable

eqnarray452

proof By the uniqueness of the orthogonal projection on a subspace, we need only demonstrate an element of tex2html_wrap_inline904 such that

displaymath806

which requires that

displaymath807

a simple case of the celebrated Wiener-Hopf equation. Now, since tex2html_wrap_inline906 , we can express it as Ay where A is some dim(x) by dim(y) matrix.

eqnarray460

So we can solve for A, possibly using psuedo-inverses, which yields:

displaymath808


next up previous contents
Next: Filtering Stochastic Control Systems Up: Spaces of Random Variables Previous: Random Variables

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