Discussion :

[pierolefou]

Salut à vous tous! Et enchanté, je suis tout nouveau sur le forum!
Ce qui m'amène sur le forum, probablement plus que pour cet échange, c'est un petit problème qui semble rien comme ça mais qui me bloque complètement dans mon travail:

Code :

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

1
2
3
4
5
6
 
>> 1 - 0.8 - 0.2
 
ans =
 
  -5.5511e-17

Au lieu de me sortir 0 comme résultat, cette opération très simple me renvoie ce petit chiffre qui m'ennuie beaucoup.
Savez vous d'où vient ce chiffre et comment je pourrais retrouver un beau "0"?
Je précise que pour toutes les opérations avec des nombres décimaux (à virgule, ou plutôt avec un point dans Matlab) me sortent des résultats avec une quantité infinitésimale en plus ou en moins.
Merci pour votre aide!
Piero le fou

[Jerome Briot]

Il n'y a aucun problème ici. Le résultat est conforme à la norme IEEE 754 à partir de la laquelle MATLAB gère les valeurs réelles.

=> : Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?

[ol9245]

Ca a beau être dans la norme IEEE 754, je confirme que c'est ausssi dans la norme "emmerdement maximum" mais c'est pas trop difficile à gérer.

Il y a deux options :
* soit tu arrondis les valeurs pour retrouver la précision que tu connais sur tes données
* soit quand tu fais des tests, tu prends en compte que tu ne peux pas tester des égalités parfaites mais des intervalles.

option 1 :

Code :

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

1
2
precision = 10^12 ;
X = round(X .* precision) ./ precision ;

option 2 :

Code :

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

1
2
3
4
5
6
7
8
9
10
% au lieu de :
if X==k
...
end
 
% préférer :
epsilon = 10^-12 ;
if abs(X - k) <= epsilon 
...
end

[FR119492]

Salut!
C'est clair que si on t'avait offert un ordinateur sur lequel le processeur travaille avec des mots comportant une infinité de bits, le problème ne se poserait pas.
Jean-Marc Blanc

[ol9245]

Salut!
C'est clair que si on t'avait offert un ordinateur sur lequel le processeur travaille avec des mots comportant une infinité de bits, le problème ne se poserait pas.
Jean-Marc Blanc

ou qui fait du calcul formel
(mais là on digresse hein )