hnormalize.m

数学优化工具箱

function [T] = hnormalize(inPoints)
  
  dims = size(inPoints,1)-1;
  nbrpoints = size(inPoints,2);
  
  for i = 1:dims
    m(i) = mean(inPoints(i,:));
  end
  
  T1 = eye(dims+1);
  for i = 1:dims
      T1(i,dims+1) = -m(i);
  end
    
  %for i = 1:dims
  %  s(i) = std(inPoints(i,:));
  %end
  
  inPoints = T1*inPoints;
  for i = 1:nbrpoints
      a(i) = sqrt(sum((inPoints(1:dims,i)).^2));
  end
  s = sqrt(dims)/mean(a);
  T2 = eye(dims+1); T2(1:dims,1:dims) = s*eye(dims);
  
  T = T2 * T1;
    
  %T = eye(dims+1);
  %for i = 1:dims
  %  T(i,i) = 1/s(i);
  %  T(i,dims+1) = -m(i)/s(i);
  %end