Discussion :

[ERCO503]

Bonjour,

Je suis à la recherche d'une astuce pour trouver une liste d'index des valeurs les plus proches.

Soit un vecteur ligne A contenant un très grand nombre de valeurs (par exemple 10 000).

Soit un vecteur ligne B contenant un très grand nombre de valeurs (par exemple 50 000).

Comment trouver les 10 000 index du vecteur B correspondant chacun à une des 50 000 valeurs du vecteur B les plus proches pour chacune des 10 000 valeurs du vecteur A.

L'astuce la plus simple et évidente consiste à effectuer ceci :

Code :

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

[~,imin] = min(abs(A-B'));

Mais cela implique la création d'une matrice 50 000 x 10 000. Cette astuce n'est donc pas appropriée à ces grands vecteurs.

J'aimerais minimiser les itérations par boucles pour une meilleure rapidité d'exécution.

Merci,

Éric

[abel413]

Bonjour,

Je suis loin de pouvoir réellement t'aider mais bon je peux toutefois partager quelques idées : tu pourrais par exemple essayer de faire un code où A et B ont même taille. On peut imaginer que tu rajoutes des 0 à A, en l’occurrence 40000 zéros après les 10000 premières valeurs :

Code :

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

1
2
A1=zeros(size(B,1),1);
A1(1:size(A))=A;

Ainsi A1 et B ont même taille : 50000. Tu peux ensuite récupérer les positions indicielles non nulles de A1 qui sont en fait les valeurs de A :

Code :

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

A3=find(A1~=0);

Ensuite il faudrait essayer de réaliser une boucle sur i variant de 1 jusqu'à la taille de A3 et essayer de déterminer les positions que tu veux. Malheureusement, quelque chose comme ce qui suit ne marche pas :

Code :

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

1
2
3
4
5
 
A4=[];
for i=1:size(A3,1)
    A4=[A4;find(min(abs(B-A1(A3(i)))))];
end

Dans A4 on doit retrouver les indices qui t'intéresse mais cela ne marche pas avec la fonction find. Enfin bon, c'est à adapter en écrivant éventuellement une ou d'autres fonctions.

J'ignore si cela t'aide un peu.

EDIT:

Ceci ne passe pas non plus chez moi :

Code :

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

1
2
3
4
5
A4=[];
imin=[];
for i=1:size(A3,1)
    [A4,imin]=min(abs(B-A1(A3(i)));
end

[Jerome Briot]

Le code suivant ne prend "que" 1 seconde sur mon portable équipé d'un simple Intel Core i5 :

Code :

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

1
2
3
4
5
6
7
A = rand(1,10000);
B = rand(1,50000);
 
imin = zeros(1,10000);
for n = 1:10000    
    [~,imin(n)] = min(abs(A(n)-B));    
end