Discussion :

[vetchang]

Bonjour à tous,
j'aimerais savoir s'il existe une fonction équivalente de la fonction randn de matalab sous en langage c . car je cherche à générer des variables aléatoire de moyenne 0 et d'écart type =1 par exemple. merci

[souviron34]

voir la fonction rand() et srand() pour intialiser.

Liens ici pour les man pages : rand et srand.

Lien ici pour une explication de l'utilisation :
http://c.developpez.com/faq/?page=no..._random_bornes

[salseropom]

Salut, oui c'est possible. Il faut que tu utilises l'algo de Box-Muller

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
#define TWOPI (6.2831853071795864769252867665590057683943387987502) /* 2 * pi */
 
/* 
   RAND is a macro which returns a pseudo-random numbers from a uniform
   distribution on the interval [0 1]
*/
#define RAND (rand())/((double) RAND_MAX)
 
/* 
   RANDN is a macro which returns a pseudo-random numbers from a normal
   distribution with mean zero and standard deviation one. This macro uses Box
   Muller's algorithm
*/
#define RANDN (sqrt(-2.0*log(RAND))*cos(TWOPI*RAND))

[DaZumba]

Code :

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

1
2
 
#define RANDN (sqrt(-2.0*log(RAND))*cos(TWOPI*RAND))

RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?

[droggo]

Fia,

RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?

+1

Et en plus, on en prend le log, alors qu'il peut renvoyer 0.

Et la valeur pour TWOPI dépasse largement le nombre de digits valides pour un double. Je sais, le compilateur va tronquer, mais ça peut donner une fausse impression de précision.

[Jean-Marc.Bourguet]

RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?

Oui. Si on veut, pour chaque paire de nombres uniformement distribues retourner une paire de nombres distribues normalement en utilisant aussi sin a la place de cos.

[salseropom]

RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?

oui c'est normal. Dans l'algorithme de Box Muller, il faut tirer 2 nombres aléatoires suivant la loi uniforme sur [0, 1]. Ces 2 nombres doivent être INDEPENDANTS !! C'est pourquoi il y a 2 fois RAND.

[salseropom]

Fia,

+1

Et en plus, on en prend le log, alors qu'il peut renvoyer 0.

Et la valeur pour TWOPI dépasse largement le nombre de digits valides pour un double. Je sais, le compilateur va tronquer, mais ça peut donner une fausse impression de précision.

Certes le log peut renvoyer 0 mais si U suit une loi uniforme sur [0, 1], la probabilité que U soit nul est égale à 0...

TWOPI dépasse largement le nombre de digits. Je le sais, mais je n'avais pas envie de me prendre la tête entre un compilo 64 bits, 32 bits, PC, LINUX, MAC. Donc sur Maple, j'ai fait un Digits=50 pour une question de "sécurité"

[salseropom]

Oui. Si on veut, pour chaque paire de nombres uniformement distribues retourner une paire de nombres distribues normalement en utilisant aussi sin a la place de cos.

.

Tout à fait, c'est l'algo de Box Muller. J'ai opté arbitrairement pour le sinus