besoin de confirmation

Version imprimable

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

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.

10/10/2006, 11h04
zooro
Citation:

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:

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.
10/10/2006, 11h06
salseropom

OK je rajoutes les parenthèses.
Merci.
10/10/2006, 11h09
homeostasie

Citation:

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).

;)
10/10/2006, 11h15
DaZumba

Citation:

Envoyé par salseropom

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

Je te conseille d'utiliser les fonctions des Numerical Recipes (sections 7.1 et 7.2) pour ce genre de chose. C'est plus robuste que rand().
10/10/2006, 12h17
salseropom

Merci pour ton lien !
10/10/2006, 12h46
jowo

Un autre lien randomnumbers si tu as besoin d'une bonne source aléatoire de nombre.