04.01.2013 Views

Springer - Read

Springer - Read

Springer - Read

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.

8.8 Generalized State-Space Models 293<br />

and that Xt+1 given Xt, X (t−1) , Y (t) is independent of X (t−1) , Y (t) with conditional<br />

density function<br />

p(xt+1|xt) : p xt+1|xt, x (t−1) , y (t)<br />

t 1, 2,.... (8.8.2)<br />

We shall also assume that the initial state X1 has probability density p1. The joint<br />

density of the observation and state variables can be computed directly from (8.8.1)–<br />

(8.8.2) as<br />

p(y1,...,yn,x1,...,xn) p yn|xn, x (n−1) , y (n−1) p xn, x (n−1) , y (n−1)<br />

p(yn|xn)p xn|x (n−1) , y (n−1) p y (n−1) , x (n−1)<br />

p(yn|xn)p(xn|xn−1)p y (n−1) , x (n−1)<br />

···<br />

<br />

n<br />

n<br />

p(yj|xj) p(xj|xj−1) p1(x1),<br />

j1<br />

and since (8.8.2) implies that {Xt} is Markov (see Problem 8.22),<br />

<br />

n<br />

p(y1,...,yn|x1,...,xn) p(yj|xj) . (8.8.3)<br />

j1<br />

We conclude that Y1,...,Yn are conditionally independent given the state variables<br />

X1,...,Xn, so that the dependence structure of {Yt} is inherited from that of the state<br />

process {Xt}. The sequence of state variables {Xt} is often referred to as the hidden<br />

or latent generating process associated with the observed process.<br />

In order to solve the filtering and prediction problems in this setting, we shall<br />

determine the conditional densities p xt|y (t) of Xt given Y (t) , and p xt|y (t−1) of Xt<br />

given Y (t−1) , respectively. The minimum mean squared error estimates of Xt based<br />

on Y (t) and Y (t−1) can then be computed as the conditional expectations, E Xt|Y (t)<br />

and E Xt|Y (t−1) .<br />

An application of Bayes’s theorem, using the assumption that the distribution of<br />

Yt given Xt, X (t−1) , Y (t−1) does not depend on X (t−1) , Y (t−1) , yields<br />

and<br />

p xt|y (t) p(yt|xt)p xt|y (t−1) /p yt|y (t−1)<br />

p xt+1|y (t) <br />

<br />

j2<br />

(8.8.4)<br />

p xt|y (t) p(xt+1|xt)dµ(xt). (8.8.5)<br />

(The integral relative to dµ(xt) in (8.8.4) is interpreted as the integral relative to dxt<br />

in the continuous case and as the sum over all values of xt in the discrete case.) The<br />

initial condition needed to solve these recursions is<br />

p x1|y (0) : p1(x1). (8.8.6)

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

Saved successfully!

Ooh no, something went wrong!