If() et les nombres à virgules

Question

Bonjour,
J'ai un prog qui tourne comme ceci : il élève un nombre x à une puissance n et la compare avec une équation contraire (n puissance x). Les X sont increnmenté nomalement (x=x+interval;).
Tout fonctionne très bien avec un interval entier ou qui est dans la suite 0.5/2, 0.5/4 ....
Mais avec un interval de 0.1 : aucune solution alors que j'ai entré exactment les meme données !!!!!

Comment faut-il faire pour résoudre ce problème ?
Mes variables sont des float (j'ai essayé avec double, ca change rien...)

Jolfulorc,
Qui vous remercie d'avance  :)

hilairenicolas · Answer

La je pense qu'il serait temps de nous montrer un bout de code

Jolfulorc · Answer

le voilà en entier (il n'est pas long, ne vous en fait pas :D) #include #include #include int main() { float x, max, resultat, resultat2, intervalle, equation, nbres=0; cout<<"|------------------------| | Math2 : Equations supp | |------------------------|"; getchar(); cout<<" Recherche des X pour l'equation n^x = x^n ... Veuillez taper le nombre N (N+) : "; cin>>equation; cout<<"Veuillez taper le nombre max de la serie (N+) : "; cin>>max; cout<<"Veuillez taper le nombre min de la serie (N-) : "; cin>>x; cout<<"Veuillez taper le nombre d'intervalle dans la serie (R+) : "; cin>>intervalle; cout<<"--------------------------------------------------------"<

ymca2003 · Answer

le test
if (resultat == resultat2)

est assez dangeureux avec des nombres à virgules. En effet, la représentation en mémoire est une troncature des nombres rééls. De plus, les nombres sont représentés en base 2 et pas 10. C'est pour ça que ça marche avec les nombres entiers et les intervales du type 0.5/2^n car dans ce cas les nombres ne sont pas tronqués. Par contre 0.1 (1/10) est arrondi en base 2 de la même façon que 1/3 en base 10.

Au lieu de faire un test d'égalite, fait un test avec une marge d'erreur (0.0001 ou 0.00001 par exemple).

Jolfulorc · Answer

Ok, merci bcp
Jolfulorc

Jolfulorc · Answer

lol, mais tu ne pourrais pas me montrer un exemple de test de marge d'erreur ? J'arrive pas à en trouver en plus c'est assez urgent  :big) 

Merci d'avance ...
Jolfulorc

ymca2003 · Answer

float val1;
float val2
float err = 0.0001f;

marge absolue :
if(fabs(val1-val2)<err)
{
//val1 et val2 égaux à 0.0001 près 
}

marge relative (tester si égaux à 0 avant de faire la division):
if(fabs(val1-val2)/val1 < err)
{
// val1 et val2 égaux à un pourcentage près
}

à toi de faire des tests pour voir ce qui est le mieux.

Jolfulorc · Answer

Merci bcp pour la technique, mais ca change rien au résultat, toujours aucun résultat !
Jolfulorc

ymca2003 · Answer

Peux-tu me mettre un exemple qui est censé donner un résultat et que le prog ne trouve pas ?

Jolfulorc · Answer

Vois par toi-même :
correct :

incorrect :

Jolfulorc

ymca2003 · Answer

On faisant le test :

float err = 0.00001f;
if (fabs(resultat-resultat2) < err)
{
}
on obtient les résultats mais il y a trop d'approximation (l'ajout de 0.1 à chaque itération et les calculs de pow accumulent les erreurs) du coup pour x 2 on obtient N^x 3.99999 (érreur de 1e-6 par raport à la valeur attendue).

en passant avec des double, l'erreur n'est plus que de 1e-15 mais les valeurs ne sont toujours pas strictement égales. Il faudra donc effectuer le test avec écart absolu.
double err = 1e-9;
if (fabs(resultat-resultat2) < err)
{
}

Jolfulorc · Answer

Ca fonctionne :D  :-p  Merci merci bcp !
Jolfulorc

If() et les nombres à virgules

11 réponses

Votre réponse

Discussions similaires