mahalan.m

来自「patten regnization source从1-14章能运行」· M 代码 · 共 28 行

function dist=mahalan(X,mu,C)
% MAHALAN 付窍扼畴厚胶 芭府甫 拌魂茄促. 
%
% 侩过:
%  dist=mahalan(X,mu,C)
%
% 汲疙:
%  青纺 X客 氦磐 mu 弊府绊 青纺 C牢 痢甸 埃狼 付窍扼畴厚胶 芭府甫 拌魂茄促. 
%    dist = (x-mu)'*inv(C)*(x-mu),
%  咯扁辑 x绰 X狼 凯氦磐捞促. 
%
% 涝仿:
%  X [dim x num_data] 涝仿 单捞磐.
%  mu [dim x 1] 氦磐 mu.
%  C [dim x dim] 青纺 C.
%
% 免仿:
%  dist [1 x num_data] 付窍扼畴厚胶 芭府.
%


[dim, num_data] = size( X );

XC = X - repmat(mu,1,num_data);
dist= sum((XC'*inv( C ).*XC')',1);

return;