Estoy tratando de mejorar el rendimiento del algoritmo de agrupación OPTICS. La implementación que he encontrado en código abierto utiliza un bucle for para cada muestra y puede ejecutarse durante horas ...

Creo que el uso de la función repmat () puede ayudar a mejorar su rendimiento cuando el sistema tiene suficiente cantidad de RAM. Le invitamos a sugerir otras formas de mejorar la implementación.

Aquí está el código:

x son los datos: una matriz [mxn] donde m es el tamaño de la muestra yn es la dimensionalidad de la característica, que es la mayoría de las veces significativamente mayor que uno.

[m,n] = size(x);

for i = 1:m
    D(i,:) = sum(((repmat(x(i,:),m,1)-x).^2),2).';
end

muchas gracias.