181187_es-ES

liliconcom

Medir la distancia entre todos los puntos

Estoy tratando de determinar una medida de cómo la distancia desde un punto afecta las mediciones para encontrar un punto mínimo/máximo dentro del rango de puntos.

 

Traté de hacer esto en los archivos adjuntos, pero no puedo encontrar una manera de hacerlo programáticamente.Traté de hacer esto concatenando una tabla de N filas, N veces, copiando la columna XY, ordenando Y luego X y luego pegando esta columna ordenada en la tabla NxN.

 

Traté de hacer esto con EVAL / SUSTITUCIÓN pero no puedo lograr que se ejecute con éxito.

Re: medir la distancia entre todos los puntos

¿Viste la función Distancia()? Aquí hay un ejemplo del índice de secuencias de comandos de una secuencia de comandos que lo usa:

 

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 );