Discussion :

[Madalen]

Bonjour à tous!

Je cherche à modéliser la mutation génétique d'une population d'individus. Pour cela, je sélectionne 10% des individus et leur applique une mutation tirée selon une gaussienne (à savoir que les petites mutations interviennent souvent et les grosses mutation rarement...) Cependant, je ne sais pas comment effectuer ce tirage aléatoire. Je pensais utiliser rand(), mais est-il possible de l'initialiser sur une gaussienne (au lieu du temps généralement utilisé.)
Merci de vos réponses.

[Médinoc]

Ce n'est pas la même chose.

• Rien dans rand() ne règle les distributions de probabilités, tu dois faire ça toi-même dans tes fonctions qui utilisent rand().
• Initialiser avec le temps sert juste à rendre les données pseudo-aléatoires moins prévisibles: Sans ça, tu auras toujours les mêmes nombres dans le même ordre à chaque lancement du programme.

[Matthieu Brucher]

Je te conseille de voir mon tutoriel et d'utiliser Boost : http://matthieu-brucher.developpez.c.../boost/random/ (et pas utiliser rand())

[uriotcea]

Voic une fonction tirant gaussian si sigma est > 0

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
double randg(double sigma)
{
  double x,y,r2;
  static double m=2./RAND_MAX;
 
  if (sigma==0) return(rand()*m-1.);
  do {
	x=rand()*m-1.;
    y=rand()*m-1.;
	r2=x*x+y*y;
  } while (r2>1.0 || r2==0);
  return (sigma*y*Psqrt(-2.*log(r2)/r2));
}

[Médinoc]

m devrait être const en plus de static, non?

[uriotcea]

Oui ca pourrait etre CONST, mais cela change t-il quelque chose au résulat ou du moins à la compilation (optimisation vitesse)?

[Médinoc]

Ça change surtout le risque de se retrouver avec une vraie variable statique (donc, durée de vie supérieure à celle de la fonction, initialisée seulement une fois, etc.) là où ce n'est pas nécessaire...

En fait, elle devrait soit être const, soit static const, soit rien du tout, mais surtout pas static tout court...

[uriotcea]

Ok, l'intérêt d'un static ici, est de ne faire ce calcul qu'une fois, sachant que moi, j'y passe des millard de fois

[Médinoc]

Sauf que dans tous les cas, la division 2./RAND_MAX sera faite par le compilo et son résultat hard-codé dans le programme...

[Madalen]

Matthieu Brucher, les fonctions que tu présentes dans ce tuto sont-elles rapides en vitesse d'exécution ? j'ai besoin de quelque chose qui tourne le plus vite possible !
En tout cas elles sont bien pratiques toutes ces méthodes.
Merci à vous.

[Matthieu Brucher]

Oui, c'est très rapide. La génération uniforme est plus robuste que le rand() et pour la transformation en bruit gaussien, c'est en gros le code proposé par uriotcea, mais en potentiellement plus raide.

[Madalen]

Ok, merci beaucoup à tous! Grâce à vous je vais pouvoir faire muter mes souris! ^^