Matlab : Problème de format dans boucle

[invite2d9939ba]

Bonjour,

J'ai un souci dans ma boucle :

cpt=0
While ..
...
cpt=cpt+0.001;
if cpt==1
 display('test1');
elseif cpt==2
 display('test2');
else
end


end

Le souci c'est que je suis sur que cpt passe par 1 et 2 mais pas moyen de rentrer dans mes conditions if et elseif??? peut être un problème de format de données???

Merci
-----

[lou_ibmix_xi]

Est tu sûr qu'aprés ta boucle "while", cpt est égale à 0.999 ? De manière générale, c'est une très mauvaise idée de comparer des nombres à virgules flottantes, il vaut mieux fair un truc du genre:

if (abs (cpt - reference) < max_diff)
    ...
end

[danyvio]

Mets une instruction juste derrière ou juste avant
cpt=cpt+0.001;

pour afficher le cpt et vérifier ta certitude de passer par les valeurs 1 ou 2, quitte à mettre un pas de +.01 pour ne pas encombrer par des milliers d'affichage

[invite2d9939ba]

j'ai tout essayé, et cpt passe bien par 1 et et 2... mais j'ai trouvé le problème, en faisant "format long" je me suis apercu qu'il ne m'ajoutait pas toujours 0.001, mais des fois 0.00100000001 ou 0.000999999999 et du coup lors de ma différence le rsultat n'est pas nul... c'est un bug de matlab je pense. Pour résoudre mon problème, j'utilise une boucle for combiner avec un if..break..

[A voir en vidéo sur Futura]

[lou_ibmix_xi]

Ce n'est absolument pas un bug de MATLAB, et c'est la démonstration de ce que je te disais:

c'est une très mauvaise idée de comparer des nombres à virgules flottantes

Regarde comment sont codés les chiffres flottants dans un ordinateur!

[ketchupi]

effectivement, ce n'est pas un bug, le codage des nombres à virgules (sur 32 bits je crois), et surtout des petits nombres, n'est pas précis. Il existe des incertitudes sur le codage des bits. Comme l'a dit Lou, mieux vaut comparer des nombres en relatif, par exemple avec l'une des inégalités dont il (elle) a parlées. Généralement, on utilise la variable eps.
++

[lou_ibmix_xi]

Matlab utilise du 64 bits pour coder les nombres, mais ca ne change rien au problème...

[HS]
C'est bien "il", mais malheureusement on ne choisit pas son prénom sur les forums...
[/HS]

[danyvio]

Je ne connais pas spécifiquement matlab, mais, m'étant frotté au cours de ma longue carrière à pas mal d'em*ts dans différents langages, je suggère ceci : au lieu de comparer sous la syntaxe :
if cpt==1

mettre

if cpt==1.0

C'est peut-être inutile, mais allez savoir...