Discussion :

[gregnorg]

Bonjour,

Voici mon problème:

J'ai le vecteur suivant C = [0 150 180 480 510 810 840 ... 150+m*330 180+m*330 ...]
J'ai un reel ro, tiré au hasard, tel que ro >0 et ro<max(C). Je cherche à trouver pour chaque ro, si 150+ m*330<ro<180+m*330 et trouver la valeur de m. Le vecteur C pouvant être trés long, j aimerais faire cela de la facon la plus rapide possible.
Est ce qu'il existe une fonction matlab permettant de faire cela? dois je programmer quelque chose? si oui quelle est la manière la plus efficace de traiter le problème?
D'avance merci.

Greg

[Jean Dumoncel]

Bonjour,

L'intervalle que tu cherches est obligatoirement 150+ m*330<ro<180+m*330, il ne peut pas être 180+ m*330<ro<150+(m+1)*330 ?

[gregnorg]

Bonjour,

L'intervalle que tu cherches est obligatoirement 150+ m*330<ro<180+m*330, il ne peut pas être 180+ m*330<ro<150+(m+1)*330 ?

Salut Magelan,

Si, ma question était mal posée. Disons que j´ai envie de déterminer à quel intervalle ro appartiens, 150+ m*330<ro<180+m*330, ou bien 180+ m*330<ro<150+(m+1)*330, le tout sans être obligé à faire des boucles.

Grégory

[Jean Dumoncel]

Il doit y avoir plusieurs façon de le faire, en voici une :

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
19
20
21
22
23
>> C = [0 150 180 480 510 810 840]
 
C =
 
     0   150   180   480   510   810   840
 
>> r0 = rand(1)*850
 
r0 =
 
  107.9388
 
>> res = C - r0
 
res =
 
 -107.9388   42.0612   72.0612  372.0612  402.0612  702.0612  732.0612
 
>> find(res>0,1,'first')
 
ans =
 
     2

Avec ce petit code, je cherche à quel moment les éléments du vecteur change de signe. Ensuite avec la parité de l'indice, tu peux déterminer si tu te trouves dans la situation 150+ m*330<ro<180+m*330 ou 180+ m*330<ro<150+(m+1)*330.

[gregnorg]

Il doit y avoir plusieurs façon de le faire, en voici une :

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
19
20
21
22
23
>> C = [0 150 180 480 510 810 840]
 
C =
 
     0   150   180   480   510   810   840
 
>> r0 = rand(1)*850
 
r0 =
 
  107.9388
 
>> res = C - r0
 
res =
 
 -107.9388   42.0612   72.0612  372.0612  402.0612  702.0612  732.0612
 
>> find(res>0,1,'first')
 
ans =
 
     2

Avec ce petit code, je cherche à quel moment les éléments du vecteur change de signe. Ensuite avec la parité de l'indice, tu peux déterminer si tu te trouves dans la situation 150+ m*330<ro<180+m*330 ou 180+ m*330<ro<150+(m+1)*330.

Ok merci beaucoup, je vais tester ca.

Gregory