Filtering problem (stochastic processes)
In the theory of
The problem of optimal non-linear filtering (even for the non-stationary case) was solved by
In general, if the
The mathematical formalism
Consider a
of the formwhere B denotes standard p-dimensional Brownian motion, b : [0, +∞) × Rn → Rn is the drift field, and σ : [0, +∞) × Rn → Rn×p is the diffusion field. It is assumed that observations Ht in Rm (note that m and n may, in general, be unequal) are taken for each time t according to
Adopting the Itō interpretation of the stochastic differential and setting
this gives the following stochastic integral representation for the observations Zt:
where W denotes standard r-dimensional Brownian motion, independent of B and the initial condition Y0, and c : [0, +∞) × Rn → Rn and γ : [0, +∞) × Rn → Rn×r satisfy
for all t and x and some constant C.
The filtering problem is the following: given observations Zs for 0 ≤ s ≤ t, what is the best estimate Ŷt of the true state Yt of the system based on those observations?
By "based on those observations" it is meant that Ŷt is
By "best estimate", it is meant that Ŷt minimizes the mean-square distance between Yt and all candidates in K:
Basic result: orthogonal projection
The space K(Z, t) of candidates is a Hilbert space, and the general theory of Hilbert spaces implies that the solution Ŷt of the minimization problem (M) is given by
where PK(Z,t) denotes the
is exactly the conditional expectation operator E[·|F], i.e.,
Hence,
This elementary result is the basis for the general Fujisaki-Kallianpur-Kunita equation of filtering theory.
More advanced result: nonlinear filtering SPDE
The complete knowledge of the filter at a time t would be given by the probability law of the signal Yt conditional on the sigma-field Gt generated by observations Z up to time t. If this probability law admits a density, informally
then under some regularity assumptions the density satisfies a non-linear stochastic partial differential equation (SPDE) driven by and called Kushner-Stratonovich equation,[10] or a unnormalized version of the density satisfies a linear SPDE called Zakai equation.[10] These equations can be formulated for the above system, but to simplify the exposition one can assume that the unobserved signal Y and the partially observed noisy signal Z satisfy the equations
In other terms, the system is simplified by assuming that the observation noise W is not state dependent.
One might keep a deterministic time dependent in front of but we assume this has been taken out by re-scaling.
For this particular system, the Kushner-Stratonovich SPDE for the density reads
where T denotes transposition, denotes the expectation with respect to the density p, and the forward diffusion operator is
where . If we choose the unnormalized density , the Zakai SPDE for the same system reads
These SPDEs for p and q are written in Ito calculus form. It is possible to write them in Stratonovich calculus form, which turns out to be helpful when deriving filtering approximations based on differential geometry, as in the projection filters. For example, the Kushner-Stratonovich equation written in Stratonovich calculus reads
From any of the densities p and q one can calculate all statistics of the signal Yt conditional on the sigma-field generated by observations Z up to time t, so that the densities give complete knowledge of the filter. Under the particular linear-constant assumptions with respect to Y, where the systems coefficients b and c are linear functions of Y and where and do not depend on Y, with the initial condition for the signal Y being Gaussian or deterministic, the density is Gaussian and it can be characterized by its mean and variance-covariance matrix, whose evolution is described by the Kalman-Bucy filter, which is finite dimensional.[10] More generally, the evolution of the filter density occurs in an infinite-dimensional function space,[5] and it has to be approximated via a finite dimensional approximation, as hinted above.
See also
- The smoothing problem, closely related to the filtering problem
- Filter (signal processing)
- Kalman filter, a well-known filtering algorithm for linear systems, related both to the filtering problem and the smoothing problem
- Extended Kalman filter, an extension of the Kalman filter to nonlinear systems
- Smoothing
- Projection filters
- Particle filters
References
- ^ Stratonovich, R. L. (1959). Optimum nonlinear systems which bring about a separation of a signal with constant parameters from noise. Radiofizika, 2:6, pp. 892-901.
- ^ Stratonovich, R.L. (1960). Application of the Markov processes theory to optimal filtering. Radio Engineering and Electronic Physics, 5:11, pp.1-19.
- ^ Kushner, Harold. (1967). Nonlinear filtering: The exact dynamical equations satisfied by the conditional mode. Automatic Control, IEEE Transactions on Volume 12, Issue 3, Jun 1967 Page(s): 262 - 267
- ^ a b Mireille Chaleyat-Maurel and Dominique Michel. Des resultats de non existence de filtre de dimension finie. Stochastics, 13(1+2):83-102, 1984.
- ^ Maybeck, Peter S., Stochastic models, estimation, and control, Volume 141, Series Mathematics in Science and Engineering, 1979, Academic Press
- ^ Damiano Brigo, Bernard Hanzon and François LeGland, A Differential Geometric approach to nonlinear filtering: the Projection Filter, I.E.E.E. Transactions on Automatic Control Vol. 43, 2 (1998), pp 247--252.
- ^ Damiano Brigo, Bernard Hanzon and François Le Gland, Approximate Nonlinear Filtering by Projection on Exponential Manifolds of Densities, Bernoulli, Vol. 5, N. 3 (1999), pp. 495--534
- .
- ^ a b c Bain, A., and Crisan, D. (2009). Fundamentals of Stochastic Filtering. Springer-Verlag, New York, https://doi.org/10.1007/978-0-387-76896-0
Further reading
- Jazwinski, Andrew H. (1970). Stochastic Processes and Filtering Theory. New York: Academic Press. ISBN 0-12-381550-9.
- ISBN 3-540-04758-1. (See Section 6.1)