Discussion :

[HAL-9000]

Salut à tous les Matlabeurs

Je travaille sur des grandes quantités de données, et je cherche à récupérer, parmis un vecteur de n lignes, les numéros de lignes correspondant aux m premières valeurs maximales. J'utilise actuellement la méthode du

Code :

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

find(ligne =  min(vecteur))

dans une boucle à m répétitions, en supprimant à chaque fois la ligne obtenue (me laissant ainsi par élimination les m-1 valeurs max).
Seulement sur un vecteur de taille m = 1 000 000 cela me prend énormément de temps, car je dois à chaque fois balayer les m valeurs afin de trouver le min et le supprimer...

Auriez une piste pour une solution moins coûteuse en temps ?

Merci à vous.

[duf42]

Bonjour,

Peut-être une idée (avec data ton vecteur de données):

Code :

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

1
2
 
find(data<inf,m,'last')

A priori ca doit donner les m plus grandes valeurs du vecteur data (à tester plus avant, je n'ai fais qu'un test rudimentaire )

Duf

[HAL-9000]

Merci duf42, ça à l'air de marcher oui, bien plus efficace

[HAL-9000]

Erf, en fait non... ta technique me donne les x valeurs inférieures à...
Or dans mon problème, je recherche les x valeurs les plus grandes.

Une autre idée ?

Merci

[duf42]

Ah oui autant pour moi, les tests n'étaient pas suffisants

Une autre solution plus simple serait peut-être de trier ton vecteur avec SORT.

[HAL-9000]

J'ai pensé à sort en effet, sauf que cela devient vite une usine à gaz dans mon cas, puisque je recherche les m valeurs les plus grandes à partir d'un vecteur P qui fait environ 20 000 lignes, et qu'à partir de ces m valeurs je retiens derrière X structures de m*n valeurs de simulations, (ou n est proche de 50, ce qui me fait des matrices de 20 000*50 à trier, et X = 1000...), et ensuite je vais refaire une seconde sélection de l valeurs sur un autre vecteur Q, etc... etc.

D'autres idées ?

[Jerome Briot]

Le tri me semble pourtant à moi aussi la solution la plus efficace... ou alors je n'ai pas bien compris le contexte du problème

[HAL-9000]

Je vais re-tester le SORT alors... Mais est-il possible de récupérer l'algo de tri qui a été effectué ? Exemple, je tri le vecteur [1,5,3,9] par ordre croissant, donc [1,3,5,9]. Est-il possible de récupérer l'algo afin que le vecteur [1,2,3,4] devienne [1,3,2,4] ?

Merci pour vos réponses.

[Jerome Briot]

Code :

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

[idx,idx] = sort([1,5,3,9])

C'est ça ?

[HAL-9000]

Oui, mais est-ce possible par la suite de faire une sorte de

Code :

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

sort([a,b,c,d], idx)

qui me donnerai [a,c,b,d] ?

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
>> A = [1,5,3,9];
 
>> [B,idx] = sort([1,5,3,9])
 
B =
 
     1     3     5     9
 
 
idx =
 
     1     3     2     4
 
>> B(idx)
 
ans =
 
     1     5     3     9

?

[HAL-9000]

C'est top

Merci Dut

[ol9245]

Si Dut a résolu ton problème, merci de cocher la case "résolu".
Moi, j'dis rien, mais quand Dut a résolu un problème, c'est en général difficile d'ajouter quelque chose d'intelligent.

DUT = thread killer (pour notre plaisir à tous )