Library: | kalman |
See also: | rlsfil calibrLS |
Quantlet: | rlsbnorm | |
Description: | Auxiliary routine for rlsfil: solves E [ |X-MYw_b(MY)|^2]=(1+e)E [ |X-MY|^2] - if possible - by MC-integration for X ~ N_n(0,Sigt), v ~ N_m(0,Q) indep. M = Sigt H'(Q+HSigt H')^{-1} Y = HX+v, w_b(x)=min(1,b/|x|) |
Usage: | {b,ctrl}=rlsbnorm(e,Sigt,H,Q,N,eps,itmax) | |
Input: | ||
e | numeric; efficiency loss to attain | |
Sigt | n x n covariance matrix of the initial state | |
H | m x n matrix | |
Q | m x m variance-covariance matrix | |
N | integer; MC-sample size | |
eps | numeric; exactitude | |
itmax | integer; maximal number of iterations | |
Output: | ||
b | numeric clipping height achieving the relative effiency loss | |
ctrl | integer; 1 if error occurs, 0 else |
to be looked up in rlsfil
n/a