Discussion :

[khoeds]

Bonjour,
J'ai deux matrices, A(n lignes,3 colones) et B(m lignes,3 colones). Je souhaite connaitre la distance minimum de chaque point de A avec l'ensemble des points de B. Pour l'instant j'ai une boucle qui ralenti enormément mon programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
for i=1:n
    dist=min(sqrt((B(:,1)-A(i,1)).^2+(B(:,2)-A(i,2)).^2+(B(:,3)-A(i,3)).^2))
end

Je recherche un moyen d'éliminer cette boucle, si quelqu'un à une idée, merci d'avance.

[Jerome Briot]

Tu trouveras quelques idées ici : [Défis][2] Pourquoi les matrices et les vecteurs ne s'entendent pas mieux ?

[paradize3]

Quelques pistes...

Pour une solution 'brute force' sous matlab, tu peux utiliser 'repmat' puis calculer toutes les distances d'un coup et ensuite faire un 'sort'.
Si c'est toujours trop lent tu peux le faire dans une fonction MEX..
Si c'est toujours trop lent, essaye une structure de donnes plus appropriee a ton probleme: par ex: kd-tree, hashed-grid, etc..

Salutations,

Greg