kalmanTools
所属分类:matlab编程
开发工具:matlab
文件大小:12KB
下载次数:14
上传日期:2018-03-03 15:31:13
上 传 者:
冬日里的影子
说明: Functions
kalman_filter
kalman_smoother - implements the RTS equations
learn_kalman - finds maximum likelihood estimates of the parameters using EM
sample_lds - generate random samples
AR_to_SS - convert Auto Regressive model of order k to State Space form
SS_to_AR
learn_AR - finds maximum likelihood estimates of the parameters using least squares
(This toolbox supports filtering, smoothing and parameter estimation (using EM) for Linear Dynamical Systems.)
文件列表:
COPYING.LESSER (7639, 2007-11-25)
kalmanFilter.m (2293, 2007-12-13)
kalmanLL.m (2108, 2007-12-06)
kalmanMLE.m (3226, 2007-12-09)
kalmanMStep0.m (1639, 2007-12-13)
kalmanMStep1.m (1796, 2007-12-13)
kalmanMStep2.m (1900, 2007-12-13)
kalmanMStep3.m (2045, 2007-12-13)
kalmanSmooth.m (2843, 2007-12-13)
The following notation is used throughout these packages:
x(t+1) = F*x(t) + B*u(t) + e(t), e(t) ~ N(0,S)
z(t) = a + A*x(t) + w(t), w(t)~ N(0,W)
x(t) is the k*1 state vector at time t
z(t) is the n*1 observation vector at time t
u(t) is the m*1 vector of known inputs at time t
F is the k*k state evolution matrix
B is the k*m input response matrix
S is the k*k covariance matrix for the Gaussian state innovations
A is the n*n observation matrix
a is an optional n*1 affine term in the measurement equation
W is the n*n covariance matrix for the Gaussian observation noise
Z is the T*n matrix of observations
U is the T*m matrix of known inputs
Vfilt & Vmfilt are n*n*T arrays containing the filtered prior and posterior
error covariance matrices for each time.
Vsmooth and VVsmmoth are n*n*T arrays containing smoothed error covariance
matrices for each time. Vsmooth is the estimated covariance between
contemporaneous state estimation errors. VVsmooth is the estimated covariance
between state estimation errors at times t & t-1.
All inputs suffixed with 0 (A0, W0, F0, S0, etc.) are the starting values
for the EM algorithm. These should be chosen quite carefully, as the EM algorithm
is only capable of finding local maximum, not the global one.
For more information on the EM algorithm & its application to linear state-space
systems, I recommend:
Ghahramani and Hinton, "Parameter Estimation for LDS", U. Toronto tech. report, 1996
Digalakis, Rohlicek and Ostendorf, "ML Estimation of a stochastic linear system with the EM
algorithm and its application to speech recognition",
IEEE Trans. Speech and Audio Proc., 1(4):431--442, 1993.
Borman, Sean, "The Expectation Maximization Algorithm, A Short Tutorial",
Tutorial Notes, University of Notre Dame EE Department,
http://www.seanborman.com/publications/EM_algorithm.pdf, July 2004.
近期下载者:
相关文件:
收藏者: