Discussion :

[el nino]

Bonjour,

Voilà j'ai juste une question toute simple mais je ne sais pas comment faire :

J'ai un vecteur [0,360]. Je veux lui assigner un pas 0:pas:360 puis pour un nombre réel donné, je veux trouver l'index correspondant à sa valeur la plus proche dans le vecteur.

Exemple :

index 1 2 3 4 5 ... 37
vec 0 10 20 30 40 ... 360

Par exemple pour 23.6, index -> 3.

J'ai essayé de trouver une formule simple entre autre, pour un réel d :

Code :

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

1
2
ratio = (nvec-1) /(vec(nvec)-vec(1))
index = round(ratio*(d-vec(1)+1))

Mais ça ne marche pas ma formule doit être fausse (par exemple d=23,6).
En attendant pour une valeur d pour trouver son index je fais find(min(abs(vec-d))) mais une formule mathématique serait plus rapide dans mes boucles.


Je vous remercie.

[Invité]

Bonjour,

Ta ligne find(min(abs(vec-d))) est fausse : min(...) te renvoie la valeur minimum, et appliquer find ne sert à rien.
[~ index] = min(abs(vec-d)).
Sinon index = round(d/pas+1)

[el nino]

Merci pour les conseils. Desole j avais mal recopié la précédente formule du find, j'avais dans mon script :

find(abs(vec-d)==min(abs(vec-d)))

Merci pour la formule, il semble que :

index = round(d/pas+1)

marche pour des vecteurs commencant a vec(1) = 0 et pour generaliser :

index = round((d-vec(1))/pas+1)

il me semble.


Merci beaucoup pour votre aide !

[le fab]

sinon en continuant d'utiliser find tu peux ecrire:

Code :

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

index= find(vec>=d,1,'first')

[Invité]

@le fab : Petit bémol sur ta solution, elle va retourner l'indice de la première valeur supérieure égale à d et non la plus proche.
Par exemple pour vec = 0:10:360; :

Code :

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

1
2
>> find(vec>=17,1,'first') % retourne 3 => OK
>> find(vec>=12,1,'first') % retourne 3 et non 2

[le fab]

pas faux, j'avais lu un peu vite