Discussion :

[einrish]

Bonjour,

Après de vaines recherches je sollicite votre aide.
Je n'arrive pas à générer des nombres aléatoires suivant une loi de puissance (exemple : p(x) = 1/x^3).
Quelle est la technique la plus efficace ?

Merci pour votre aide

ps : je travaille en fortran

[Nebulix]

L'intégrale de ton p(x) entre 0 et l'infini est infinie ...

[einrish]

oui pardon, prenons p(x) = 1/(1+x)^a

[Aleph69]

Bonjour,

à quel domaine appartient x?

[benDelphic]

cherches sur Google avec les mots " générateur aléatoire" ou "simulation de lois "
. pour générer des nombres aléatoires suivant une loi quelconque il suffit de l'écrire en fonction d'une loi uniforme.

[Nebulix]

oui pardon, prenons p(x) = 1/(1+x)^a

si y a une distribution uniforme sur 0..1 alors
x=(y)^(1/(1-a))-1, sauf erreur

[Aleph69]

Bonjour,

je trouve x=p(x)^(-1/a)-1 mais tant qu'on ne sait pas à quel domaine appartient x rien n'est sîr.