Code covered by the BSD License  

Highlights from
UD Factorization & Kalman Filtering

UD Factorization & Kalman Filtering


Gerard Van Willigenburg


15 Aug 2011 (Updated )

UD and LD factorization of nonnegative matrices and associated Kalman filter implementations.

% LDT2SYM: Recover P=L'DL, LD an LDL factorization of P
%          L upper triangular, D diagonal, P nonnegative symmetric
%          [p]=udt2sym(ld)
%          p=l'*d*l
%          Used to compute : inv(p)=inv(l)'*inv(d)*inv(l)
%          [ldi]=ldinv(ld) computes inv(l),inv(l)
%          Then: inv(p)=ldt2sym(ldinv(sym2ud(p)))
%          See also ldinv, ld2sym, sym2ld
%          References: Factorization methods for discrete sequential estimation
%                      1977, Gerald J. Bierman
% L.G. van Willigenburg, W.L. de Koning, Update August 2011

  function [p]=ldt2sym(ld)
  if n~=m; error('  ld must be square'); end;
  if n==0; error('  Compatible but empty inputs'); end;


Contact us