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

Tout ce qui concerne par exemple les outils de calcul numérique, de calcul formel ou de géométrie.
[participation réservée aux membres inscrits]
Règles du forum
Merci de soigner la rédaction de vos messages et de consulter ce sujet avant de poster. Pensez également à utiliser la fonction recherche du forum.
melbachir

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

Message non lu par melbachir »

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,
guiguiche
Modérateur général
Modérateur général
Messages : 8129
Inscription : vendredi 06 janvier 2006, 15:32
Statut actuel : Enseignant
Localisation : Le Mans

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

Message non lu par guiguiche »

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

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

Message non lu par melbachir »

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 :(
guiguiche
Modérateur général
Modérateur général
Messages : 8129
Inscription : vendredi 06 janvier 2006, 15:32
Statut actuel : Enseignant
Localisation : Le Mans

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

Message non lu par guiguiche »

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

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

Message non lu par melbachir »

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é...
guiguiche
Modérateur général
Modérateur général
Messages : 8129
Inscription : vendredi 06 janvier 2006, 15:32
Statut actuel : Enseignant
Localisation : Le Mans

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

Message non lu par guiguiche »

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.