181187_fr-FR

liliconcom

Mesurer la distance entre tous les points

J'essaie de déterminer une mesure de la façon dont la distance d'un point affecte les mesures pour trouver un point minimum/maximum dans la plage de points.

 

J'ai essayé de le faire dans les fichiers joints mais je ne trouve pas le moyen de le faire par programme.J'ai essayé de le faire en concaténant un tableau de N lignes, N fois, en copiant le coulmn XY, en triant Y puis X puis en collant cette colonne triée dans le tableau NxN.

 

J'ai essayé de le faire avec EVAL / SUBSTITUTION mais je n'arrive pas à l'exécuter avec succès.

Re : Mesurer la distance entre tous les points

Avez-vous vu la fonction Distance() ? Voici un exemple tiré du Scripting Index d'un script l'utilisant :

 

Names Default To Here( 1 );
/*1-D example*/
exX1 = [1, 2, 3, 4];
exX2 = [2, 4, 6, 8]; 
/*Compute squared Euclidean distance*/
exD = Distance( exX1, exX2 ); 
/*Verify result*/
exDm = J( 4, 4, . );
For( exi = 1, exi <= 4, exi++,
 For( exj = 1, exj <= 4, exj++,
  exDm[exi, exj] =
  Sum( (exX1[exi, 0] - exX2[exj, 0]) ^ 2 )
 )
);
Show( exDm == exD ); 

/*2-D example*/
exX1 = [1 1, 2 2, 3 3, 4 4];
exX2 = [2 1, 4 2, 6 0, 8 7]; 
/*Compute squared Euclidean distance*/
exD = Distance( exX1, exX2 ); 
/*Verify result*/
exDm = J( 4, 4, . );
For( exi = 1, exi <= 4, exi++,
 For( exj = 1, exj <= 4, exj++,
  exDm[exi, exj] =
  Sum( (exX1[exi, 0] - exX2[exj, 0]) ^ 2 )
 )
);
Show( exDm == exD ); 

/*2-D example*/
exX1 = [1 1, 2 2, 3 3, 4 4];
exX2 = [2 1, 4 2, 6 0, 8 7]; 
/*Compute squared Euclidean distance, with a scaler [0.5 2.0]*/
exD = Distance( exX1, exX2, [0.5 2.0] ); 
/*Verify result*/
exDm = J( 4, 4, . );
For( exi = 1, exi <= 4, exi++,
 For( exj = 1, exj <= 4, exj++,
  exDm[exi, exj] =
  Sum(
   [0.5 2.0] :*</img> (
   Abs( exX1[exi, 0] - exX2[exj, 0] ) ^ 2)
  )
 )
);
Show( exDm == exD ); 

/*2-D example*/
exX1 = [1 1, 2 2, 3 3, 4 4];
exX2 = [2 1, 4 2, 6 0, 8 7]; 
/*Compute squared Euclidean distance, with scalers [0.5 2.0] and powers [1.5 2.0]*/
exD = Distance( exX1, exX2, [0.5 2.0], [1.5 2.0] ); 
/*Verify result*/
exDm = J( 4, 4, . );
For( exi = 1, exi <= 4, exi++,
 For( exj = 1, exj <= 4, exj++,
  exDm[exi, exj] =
  Sum(
   [0.5 2.0] :*</img> (
   Abs( exX1[exi, 0] - exX2[exj, 0] ) :^ [
   1.5 2.0])
  )
 )
);
Show( exDm == exD );