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.

% UTT2SYM : Recovers P=UT'*UT, UT upper triangular, P nonnegative symmetric
%           [p]=utt2sym(ut)
%           p=ut'*ut
%           Used to compute : inv(p)=inv(u)'*inv(u)
%           [uti]=utinv(ut) computes inv(u)
%           Then: inv(p)=utt2sym(utinv(sym2ud(p)))
%           See also utinv, ut2sym, sym2ut
%           References: Factorization methods for discrete sequential estimation
%                       1977, Gerald J. Bierman
% L.G. van Willigenburg, W.L. de Koning, Update August 2011

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


Contact us