mahalanobis.m

来自「用EKF」· M 代码 · 共 10 行

function d2 = mahalanobis (z, P)
global calculations;
%-------------------------------------------------------
% 
%-------------------------------------------------------

Cp = chol(P);
y = Cp'\z;
d2 = full(y'*y);