[Scilab] Simulation de la loi géométrique

Tout ce qui concerne les différents logiciels mathématiques (Matlab, Maxima, Scilab, R, Maple, Mathematica, ...) et les calculatrices scientifiques.
Règles du forum
Merci d'éviter le style SMS dans vos messages et de penser à utiliser la fonction Recherche avant de poster un message. Pour joindre des fichiers à vos messages, consulter ce sujet.
> Penser à utiliser les balises Code pour poster du code.
> Préciser le nom du logiciel concerné dans le titre de votre message.

[Scilab] Simulation de la loi géométrique

Messagepar melbachir » Vendredi 30 Juillet 2010, 21:54

Bonjour à tous,

Je veux simuler la loi géométrique à l'aide de Scilab. Donc j'écris :

Code: Tout sélectionner
//- Loi géométrique de paramètre p
function X = geo(p)
  X = 1;
  while %t
  //- Si échec +1, si Succès sort
    if(rand() < p) then break; else, X = X + 1;  end
  end 
endfunction


Cette fonction si je l'ai bien codée, nous donne le rang du premier succès. Mon problème est celui là :
Comment puis-je représenter P(X=k) en fonction de k.

Merci d'avance,
melbachir
Hecto-utilisateur
 
Messages: 88
Inscription: Mercredi 17 Septembre 2008, 23:03
Statut actuel: Post-bac | Licence

Publicité

Re: [Scilab] Simulation de la loi géométrique

Messagepar guiguiche » Samedi 31 Juillet 2010, 08:10

Il faut répéter ta fonction un grand nombre de fois et faire compter la fréquence d'apparition de k (plus le nombre est grand plus sa probabilité d'être proche de la P(X=k) est importante).
Pas d'aide par MP : les questions sont publiques, les réponses aussi.
Tu as apprécié l'aide qui t'a été fournie ? Alors n'hésite pas à rendre la pareille à quelqu'un d'autre.
Un peu d'autopromotion.
guiguiche
Modérateur
 
Messages: 8063
Inscription: Vendredi 06 Janvier 2006, 15:32
Localisation: Le Mans
Statut actuel: Actif et salarié | Enseignant

Re: [Scilab] Simulation de la loi géométrique

Messagepar melbachir » Samedi 31 Juillet 2010, 10:30

Merci pour ta réponse,

Alors ce qui me vient à l'esprit :

Code: Tout sélectionner
//- Taille de l'échantillon
  N = 10000;
//- Le paramètre p
  p = 0.3;
//- Créer un vecteur colonne de composantes nulles
  Z = zeros(N,1);
//- Remplir ce vecteur
  for k=1:N
     Z(k) = geo(p);
  end


Mais après je ne sais pas quoi faire pour obtenir P(X=k) en fonction de k :(
melbachir
Hecto-utilisateur
 
Messages: 88
Inscription: Mercredi 17 Septembre 2008, 23:03
Statut actuel: Post-bac | Licence

Re: [Scilab] Simulation de la loi géométrique

Messagepar guiguiche » Samedi 31 Juillet 2010, 13:06

J'aurai mis un truc du genre (je ne connais pas la syntaxe de scilab) :
Code: Tout sélectionner
k=geo(p) ;
Z(k)=Z(k)+1 ;
Pas d'aide par MP : les questions sont publiques, les réponses aussi.
Tu as apprécié l'aide qui t'a été fournie ? Alors n'hésite pas à rendre la pareille à quelqu'un d'autre.
Un peu d'autopromotion.
guiguiche
Modérateur
 
Messages: 8063
Inscription: Vendredi 06 Janvier 2006, 15:32
Localisation: Le Mans
Statut actuel: Actif et salarié | Enseignant

Re: [Scilab] Simulation de la loi géométrique

Messagepar melbachir » Samedi 31 Juillet 2010, 16:13

guiguiche a écrit:J'aurai mis un truc du genre (je ne connais pas la syntaxe de scilab) :
Code: Tout sélectionner
k=geo(p) ;
Z(k)=Z(k)+1 ;



Mon gros problème, c'est que j'obtiens sur l'axe des ordonnées des nombres > 1 alors que je m'attends à voir une probabilité...
melbachir
Hecto-utilisateur
 
Messages: 88
Inscription: Mercredi 17 Septembre 2008, 23:03
Statut actuel: Post-bac | Licence

Re: [Scilab] Simulation de la loi géométrique

Messagepar guiguiche » Samedi 31 Juillet 2010, 16:34

Quelle est la définition d'une fréquence ?
Pas d'aide par MP : les questions sont publiques, les réponses aussi.
Tu as apprécié l'aide qui t'a été fournie ? Alors n'hésite pas à rendre la pareille à quelqu'un d'autre.
Un peu d'autopromotion.
guiguiche
Modérateur
 
Messages: 8063
Inscription: Vendredi 06 Janvier 2006, 15:32
Localisation: Le Mans
Statut actuel: Actif et salarié | Enseignant


Retourner vers Logiciels mathématiques

 


  • Articles en relation
    Réponses
    Vus
    Dernier message

Qui est en ligne

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