(a) Suppose that
(XY)∼N((μXμY),(VXXVYXVXYVYY))
Prove that conditional on Y=y, the distribution of X is again multivariate normal, with mean μX+VXYVYY−1(y−μY) and covariance VXX−VXYVYY−1VYX.
(b) The Rd-valued process X evolves in discrete time according to the dynamics
Xt+1=AXt+εt+1
where A is a constant d×d matrix, and εt are independent, with common N(0,Σε) distribution. The process X is not observed directly; instead, all that is seen is the process Y defined as
Yt=CXt+ηt,
where ηt are independent of each other and of the εt, with common N(0,Ση) distribution.
If the observer has the prior distribution X0∼N(X^0,V0) for X0, prove that at all later times the distribution of Xt conditional on Yt≡(Y1,…,Yt) is again normally distributed, with mean X^t and covariance Vt which evolve as
X^t+1Vt+1=AX^t+MtCT(Ση+CMtCT)−1(Yt+1−CAX^t)=Mt−MtCT(Ση+CMtCT)−1CMt
where
Mt=AVtAT+Σε
(c) In the special case where both X and Y are one-dimensional, and A=C=1, Σε=0, find the form of the updating recursion. Show in particular that
Vt+11=Vt1+Ση1
and that
Vt+1X^t+1=VtX^t+ΣηYt+1
Hence deduce that, with probability one,
t→∞limX^t=t→∞limt−1j=1∑tYj