besoin de confirmation

**salseropom** · 10/10/2006, 10h57

Bonjour, je désire tirer un nombre aléatoire suivant une loi normale. Voici mon code

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
16
17
18
 
/* 
   RAND is a macro which returns a pseudo-random numbers from a uniform
   distribution on the interval [0 1]
*/
#define RAND ((double) rand())/((double) RAND_MAX)
 
/*
   TWOPI = 2*pi
*/
#define TWOPI 6.28318530717959
 
/* 
   RANDN is a macro which returns a pseudo-random number 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)

Est-ce que les deux RAND de ma macro RANDN sont bien indépendants ? Je pense que oui car le #define remplace "bêtement" le code par son expression. Est-ce que je me trompe ?
Merci.

**zooro** · 10/10/2006, 11h04

Envoyé par salseropom

Est-ce que les deux RAND de ma macro RANDN sont bien indépendants ? Je pense que oui car le #define remplace "bêtement" le code par son expression. Est-ce que je me trompe ?
Merci.

Je dirais qu'ils sont indépendants, puisque remplacés par le préprocesseur. Mais si j'étais toi, je rajouterais des parenthèses autour des deux expressions dans les #define:

Code :

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

1
2
3
#define RAND ( ((double) rand())/((double) RAND_MAX) )
#define RANDN ( sqrt(-2.0*log(RAND))*cos(TWOPI*RAND) )

Pour éviter des effets de bord possibles lors du remplacement.

**salseropom** · 10/10/2006, 11h06

OK je rajoutes les parenthèses.
Merci.

**homeostasie** · 10/10/2006, 11h09

Est-ce que les deux RAND de ma macro RANDN sont bien indépendants ?

Je dirais que oui même si ce n'est pas vraiment le cas avec la fonction rand(), qui spécifie qu'elle retourne des valeurs pseudo aléatoires.
C'est à dire qu'à chaque lancement du programme, ces nombres pseudo-aléatoires seront toujours les mêmes. Ce qui peut être contraignant. Une solution pour parer à cela est d'initialiser le générateur de nombres pseudo aléatoires avec l'heure courante de ton système (Utilisation de srand).