ddistm.m

模式识别工具箱,希望对大家有用!

function d = ddistm(a,b,dtype,par1,par2)
% D = ddistm(A,B,dtype,p1,p2)
%
% Distance computation between two sets A and B. The following
% distances are defined:
%
%  dtype:       meaning:
%  'lp'         (sum (|x_i - y_i|.^p1))^(p2/p1)
%  'gauss'      1-exp( -((x_i - y_i).^2))/(p1*p1) )

% Copyright: D. Tax, R.P.W. Duin, [e-mail]
% Faculty of Applied Physics, Delft University of Technology
% P.O. Box 5046, 2600 GA Delft, The Netherlands

if nargin<5
  par2 = 1;
end
if nargin<4
  par1 = 2;
end
if nargin<3
  dtype = 'lp';
end

switch dtype
  case 'lp'
    d = lpdistm(+a,+b,par1).^par2;
  case 'gauss'
    d = 1-exp(-(lpdistm(+a,+b,2)/par1).^2);
  otherwise
    error('Unknown distance requested');
end

return