C++

Problème de multiplication de doubles

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : Problème de multiplication de doubles

inblack

Bonjour tout lemonde

J'ai un probème que je n'arrive pas à résoudre depuis quelques jours et ce n'est pas faute d'avoir chérché...
Voici mon code

Code :

double x=exp(-((_Touverture+Tarc)-_Tfermeture)/_T);
double y=_ModuleI*sin(0+(_alpha-_phi)*pi/180);
double z=_ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000+0+(_alpha-_phi)*pi/180);
double ph1tc = z-y*x;

J'ai ph1tc=-9.2559631349317831e+061 alors que z=0.81462211321630484, x=0.25394129165067114 et y=-0.99026806534391931

Le problème doit venir de la reconversion en décimale aprés la multiplication.
Est ce que quelqu'un peut m'aider??
J'avais penser à limiter le nombre de chiffres après la virgule dans tout les calculs (opérations arithmétiques mais aussi dans les calculs de sinus et d'exponentiels) mais ca je ne sais pas faire et je dois avoir une bonne précision donc pas mal de chiffres...
Merci baucoup pour votre aide..

Publicité

Taz

bisounours-codeur

euh 0.8 - (-0.99 * 0.25) ça fait bien du 1.06.

Code :

cat mult.cpp
#include <iostream>
int main()
{
double z=0.81462211321630484, x=0.25394129165067114, y=-0.99026806534391931;
double ph1tc = z - y * x;
std::cout << ph1tc << '\n';
}
[l0247308@x-pau-101624 tmp]$ g++ -Wall mult.cpp && ./a.out
1.06609

inblack

Je suis sous visual c++ et ca ne marche pas pareil...

Message cité 1 fois

Elmoricq

Modérateur

Je ne vois pas bien ce que le compilateur vient faire là-dedans.

Y a pas un "e+??" après z, x ou y ?

Taz

bisounours-codeur

t'as compilé le même bout de code que moi et tu n'as pas de résultat similaire ?

kyntriad

inblack a écrit :

Je suis sous visual c++ et ca ne marche pas pareil...

moi aussi et ça marche pareil (encore heureux d'ailleurs)

Message édité par kyntriad le 13-06-2008 à 13:30:53

---------------
You can't start a fire with moonlight

inblack

Taz quand j'essaie ton code ca marche mais en réalité les valeurs reelle de x y et z(celles calculés par le compilateur) doivent contenir bcp plus de chiffres apres la virgule (plus que les 64 bits autorisés par un double donc ils sont tranqués et stockés dans x,y et z) et donc puisque les valeurs de x y et z correspondent déja au maximum que peut contenir un double, le fait de les multiplier donne des valeurs qui devrait etre stockés sur 128 bits par exemple...ce qui conduit a des valeurs aberrantes si on les stocke sur 64 bits....
C'est mon intérprétation.... est ce que je me trompe...
c'est pour ca que j'ai pensé a limiter le nombre des chiffres apres la virgule dans tous les calculs que je fait...
en fait j'ai décomposé ca en x y z pour voir d'ou vient le probleme mais j'ai pleins de calculs comme ca a faire et il me faut une méthode qui marche a tout les coups...

Elmoricq

Modérateur

Tu te trompes. T'as beau multiplier 0.33 ou 0.330000000000000012468416574687416354124254, ça ne donnera jamais 9e+61.

Message édité par Elmoricq le 13-06-2008 à 14:35:33

kyntriad

Est-ce que toutes tes variables intermédiaires sont bien des doubles aussi ?

Envoie un peu plus de code sinon qu'on puisse voir le problème. Ca m'étonnerai fort que tes souçis soient dus à un dépassement des capacité.

---------------
You can't start a fire with moonlight

inblack

Code :

int A=0;
int & a=A;
double & I1ret=I1exact;
double Equations_CO::I1(double t)
{
double ph1=_ModuleI*sin(2*pi*_f*t/1000+0+(_alpha-_phi)*pi/180)-_ModuleI*sin(0+(_alpha-_phi)*pi/180)*exp(-(t-_Tfermeture)/_T);
double x = exp(-((_Touverture+Tarc)-_Tfermeture)/_T);
double y = _ModuleI*sin(0+(_alpha-_phi)*pi/180);
double z = _ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000+0+(_alpha-_phi)*pi/180);
double ph1tc = z-y*x;
double d = 1.732/2*_ModuleI*sin(2*pi*_f*t/1000-pi/6+(_alpha-_phi)*pi/180);
double e = (_Touverture+Tarc)/_T;
double f = (ph1tc-1.732/2*_ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000-pi/6+(_alpha-_phi)*pi/180)*exp(-(t-_Tfermeture)/_T));
double ph1bi = d-e*f;

Message cité 1 fois

Publicité

inblack

au niveau du debogage je m'arrete la donc c pas la peine de vous donner la suite...sinon pour les valeurs des variable utilisées:

Code :

Equations_CO::Equations_CO()
{
this->_f=50;
this->_ModuleI=1;
this->_T=23;
this->_phi=82.12l;
this->_alpha=0;
this->_t=0l;
this->_Tarcmin=5.82l;
this->_Tfermeture=20.00l;
this->_Touverture=100.00l;
this->_Tessai=120.00l;
this->_Tarcmin=5.82l;
this->_Parc=40l;
this->_Marc=5*_Parc;
this->_Pas=0.015l;
this->_offset=0l;
this->_bruit=0l;
}

Message cité 1 fois

inblack

j'ai mis des "l" à la suite des valeurs parce que j'avais lu un truc dans un bouquin que j'ai pas tres bien compris d'ailleurs mais en fait ca ne change rien à mon problème...

bjone

Insert booze to continue

inblack a écrit :

Code :

int A=0;
int & a=A;
double & I1ret=I1exact;
double Equations_CO::I1(double t)
{
double ph1=_ModuleI*sin(2*pi*_f*t/1000+0+(_alpha-_phi)*pi/180)-_ModuleI*sin(0+(_alpha-_phi)*pi/180)*exp(-(t-_Tfermeture)/_T);
double x = exp(-((_Touverture+Tarc)-_Tfermeture)/_T);
double y = _ModuleI*sin(0+(_alpha-_phi)*pi/180);
double z = _ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000+0+(_alpha-_phi)*pi/180);
double ph1tc = z-y*x;
double d = 1.732/2*_ModuleI*sin(2*pi*_f*t/1000-pi/6+(_alpha-_phi)*pi/180);
double e = (_Touverture+Tarc)/_T;
double f = (ph1tc-1.732/2*_ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000-pi/6+(_alpha-_phi)*pi/180)*exp(-(t-_Tfermeture)/_T));
double ph1bi = d-e*f;

ça te dit pas de sortir quelques constantes ?

Message édité par bjone le 13-06-2008 à 15:23:15

inblack

deja fait..regarde plus haut

Message cité 1 fois

kyntriad

Il manque les déclarations de _ModuleI, _phi, _alpha, _T, _t...

---------------
You can't start a fire with moonlight

inblack

inblack a écrit :

au niveau du debogage je m'arrete la donc c pas la peine de vous donner la suite...sinon pour les valeurs des variable utilisées:

Code :

Equations_CO::Equations_CO()
{
this->_f=50;
this->_ModuleI=1;
this->_T=23;
this->_phi=82.12;
this->_alpha=0;
this->_t=0;
this->_Tarcmin=5.82;
this->_Tfermeture=20.00;
this->_Touverture=100.00;
this->_Tessai=120.00;
this->_Tarcmin=5.82;
this->_Parc=40;
this->_Marc=5*_Parc;
this->_Pas=0.015;
this->_offset=0;
this->_bruit=0;
}

c'est ce que tu veux?
je l'avais posté plus haut mais en
fait je retire les "l" parce que ca ne sert a rien

inblack

Sinon t'as ca aussi:

Code :

class Equations_CO
{
/*
f: Frequence du signal a générer
ModuleI: Module du courant(Courant maximun).
ModuleV: Module de tension (Tension maximum).
T: Constante de temps.
phi: Déphasage entre le courant et la tension.
alpha: Etablissement de I1 par rapport au zero en V1 (en °).
bruit: Bruit blanc à ajouter au signal.
Offset:
Tarcmin: Temps d'arc (caractéristique des appareils a tester).
Tfermeture: Instant de fermeture des contatcts.
Touverture: Instant de séparation des contacts.
Parc: Pied d'arc caractéristique des appareils a tester (permet de determiner la tension d'arc).
Marc: Maximum de l'arc (permetant aussi la détermination de la tension d'arc.
t: temps
*/
private:
int _f,_ModuleI,_ModuleV,_T;
double _phi,_alpha,_bruit,_offset,_Tarcmin,_Tfermeture,_Touverture,_Parc,_Marc,_t;
public:
double _Tessai,_Pas;
Equations_CO();//Constructeur par defaut
Equations_CO(int,int,int,int,double,double,double,double,
double,double,double,double,double,double,double,double);
~Equations_CO(void);
double I1(double);
double I2(double);
double I3(double);
double V1(double);
double V2(double);
double V3(double);
};

bjone

Insert booze to continue

inblack a écrit :

deja fait..regarde plus haut

non:

Code :

double x = exp(-((_Touverture+Tarc)-_Tfermeture)/_T);
double y = _ModuleI*sin(0+(_alpha-_phi)*pi/180);
double z = _ModuleI*sin(2*pi*_f*(_Touverture+Tarc)/1000+0+(_alpha-_phi)*pi/180);

ça peut faire:

Code :

const double Toa = _Touverture+Tarc;
const double alphi = ( _alpha - _phi )*M_PI/180.;
double x = exp( ( Tfermeture -Toa ) / _T );
double y = _ModuleI * sin( alphi );
double z = _ModuleI * sin( 2*M_PI * _f * Toa / 1000. + alphi );

enfin simplifier, réduire, réutiliser des calculs, aérer, tout ça....
et les variables avec _ ça pue.
pi est déjà défini dans maths.h (cmaths) comme M_PI, y'a un define a mettre stout.

Message édité par bjone le 13-06-2008 à 16:04:26

Taz

bisounours-codeur

utilise un truc de vérif mémoire. tu dois bousiller ta pile à tous les coups. Ou bien fourni nous un code minimal et compilable.

FORUM HardWare.fr

Programmation

C++

Problème de multiplication de doubles

Sujets relatifs
[Flash/actionscript]probleme clip qui reste	Problème Flash -> PHP
Probleme de requete a la base de donnees	probleme de socket c/java
Probleme : Executer un fichier bat sans le sauvegarder	Problème d'affichage de tableau : vertical Firefox / Horizontal IE
Probleme d'installation JAVA	Problème fonction excel
Probleme d'affichage perl/tk	[OPENGL SFML C++] Probléme executable.
Plus de sujets relatifs à : Problème de multiplication de doubles

Page générée en 0.094 secondes