Proba : loi uniforme

Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
djocin
Messages: 4
Enregistré le: 02 Juil 2009, 15:05

Proba : loi uniforme

par djocin » 02 Juil 2009, 15:09

Bonjour

Voici mon problème : dans un programme info, après avoir tiré des nombres aléatoires x suivant une loi uniforme sur [0,1], je fais l'opération b=min+x*err.

Pour les x, le graphe de la densité de probabilité est bien un échelon entre 0 et 1, mais par contre le graphe de la densité de proba de b ne donne plus un échelon... Est-ce normal ?

Merci de votre aide.



kazeriahm
Membre Irrationnel
Messages: 1608
Enregistré le: 04 Juin 2006, 09:49

par kazeriahm » 02 Juil 2009, 18:08

Salut

si min et err sont des constantes, tu dois avoir un échelon entre min et min+err

parcontre si min et err sont elles memes des variables aléatoires alors...

djocin
Messages: 4
Enregistré le: 02 Juil 2009, 15:05

par djocin » 02 Juil 2009, 18:44

non non min et err sont des constantes

pourtant je n'ai pas l'échelon attendu... :marteau:

kazeriahm
Membre Irrationnel
Messages: 1608
Enregistré le: 04 Juin 2006, 09:49

par kazeriahm » 02 Juil 2009, 19:32

t'as quoi ? tu codes en quoi ?

djocin
Messages: 4
Enregistré le: 02 Juil 2009, 15:05

par djocin » 03 Juil 2009, 08:32

Je code en fortran.

Mais en fait j'ai refait des graphiques ce matin avec les idées au clair et j'ai bien l'échelon comme prévu. Le vrai problème vient d'ailleurs.

Je fais ces deux boucles :

Code: Tout sélectionner


do i=1,n
x(i)=-log(uniform())
tmp(i)=min(i)+x(i)*err(i)
S=S+tmp(i)
enddo

do i=1,n
RB(i)=tmp(i)/S
enddo


où RB et tmp sont des vecteurs de dimension n.

Le but de ces boucles est d'obtenir les RB, tels que leur somme fasse 1.
Uniform() est un tirage aléatoire uniforme entre 0 et 1.

Le problème est que la répartition des RB(i) est censée être uniforme également, alors que quand je trace le graphe j'ai plutôt une répartition normale, voire log-normale.

J'attends votre aide ! :briques:

 

Retourner vers ✯✎ Supérieur

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 31 invités

Tu pars déja ?



Fais toi aider gratuitement sur Maths-forum !

Créé un compte en 1 minute et pose ta question dans le forum ;-)
Inscription gratuite

Identification

Pas encore inscrit ?

Ou identifiez-vous :

Inscription gratuite