A voir également:
- Simulation d'une variable de pareto
- Vba range avec variable ✓ - Forum VB / VBA
- Simulation place accor hotel arena - Forum Loisirs / Divertissements
- Simulateur de navigation maritime gratuit - Télécharger - Loisirs créatifs
- Logiciel simulation ombre solaire - Télécharger - Architecture & Déco
3 réponses
Salut
Je me pose exactement la même question mais en Matlab. Enfin le principe est le même. Tu peux me le dire si tu as la réponse. Merci!
Je me pose exactement la même question mais en Matlab. Enfin le principe est le même. Tu peux me le dire si tu as la réponse. Merci!
/**
* Return a real number with a Pareto distribution with parameter alpha.
*/
public static double pareto(double alpha) {
return Math.pow(1 - uniform(), -1.0/alpha) - 1.0;
}
* Return a real number with a Pareto distribution with parameter alpha.
*/
public static double pareto(double alpha) {
return Math.pow(1 - uniform(), -1.0/alpha) - 1.0;
}
##inverse de la fonction de répartition d'une pareto:
invfdrpareto=function(y,l,a)
{x=l*(1/(1-y)^(1/a)-1)
x}
##génération de nombres suivant la loi de Pareto :
rpareto=function(l,a)
{f=function(n=1)
{U=runif(n)
Y=sapply(U,invfdrpareto,l,a)
Y}
f}
##Exemple d'exécution:
V=rpareto(1.4,1.5e6)(1000)
plot(V)
hist(V,prob=TRUE)
invfdrpareto=function(y,l,a)
{x=l*(1/(1-y)^(1/a)-1)
x}
##génération de nombres suivant la loi de Pareto :
rpareto=function(l,a)
{f=function(n=1)
{U=runif(n)
Y=sapply(U,invfdrpareto,l,a)
Y}
f}
##Exemple d'exécution:
V=rpareto(1.4,1.5e6)(1000)
plot(V)
hist(V,prob=TRUE)