Simulation de la planche de Galton

Simulation d'une variable aléatoire suivant une loi de Bernoulli de paramètre p

Lorsqu'une bille lancée sur une planche de Galton rencontre un clou, celle-ci a une probabilité 

p

 de passer à droite de ce clou et une probabilité 

1 - p

 de passer à gauche.
Cette expérience à deux issues est donc une expérience de Bernoulli. Pour la modéliser, on  utilise une fonctionbernoulliqui prendra en argument un réel

p

 entre 0 et 1 et qui renverra 0 si la bille passe à gauche et 1 si la bille passe à droite.

Compléter la fonction suivante qui permet de simuler une variable aléatoire suivant une loi de Bernoulli de paramètre p.

from random import random

def bernoulli(p) :
    # On génère un nombre au hasard entre 0 et 1. On le stocke dans une variable a.
    a = random()
    # Si ce nombre est inférieur ou égal à p, on renvoie 1. Sinon, on renvoie 0.
    if ... :
        return ...
    else :
        return ...

Trajectoire d'une bille

On note

h

 la hauteur de la planche de Galton, c'est-à-dire le nombre de rangées de clous que les billes rencontreront avant d'atteindre les compartiments en bas de la planche.
Pour simuler la trajectoire de la bille, on construit alors une liste de

h

 nombres générés aléatoirement selon une loi de Bernoulli de paramètre 

p

.

Exercice 1

Compléter la fonction suivante qui permet de simuler la trajectoire d'une bille lancée sur la planche de Galton. La liste renvoyée pourra avoir été générée par ajouts successifs avec la méthode append ou par compréhension.

def trajectoire(h, p):
    """
    Renvoie une liste de 0 et de 1 correspondant à la trajectoire d'une bille lancée sur la planche de Galton de hauteur h.
    Les 0 et 1 sont générés en utilisant une loi de Bernoulli de paramètre 

p

."""
    # A compléter

trajectoire(10, 0.5)

On suppose par la suite que les compartiments d'arrivée en bas de la planche de Galton sont numérotés de 0 à

h

 (il y en a donc 

h+1

 au total !).

Exercice 2

On suppose que la trajectoire d'une bille est donnée par [0, 1, 0, 1, 0, 0, 1].

Quelle est la hauteur de la planche de Galton ?
Quel est le numéro du compartiment dans lequel la bille est arrivée ?
Reprendre les questions précédentes si la trajectoire de la bille est donnée par la liste [1, 0, 1, 1, 1, 0].
D'une manière générale, si l'on se donne une trajectoire d'une bille, comment peut-on déterminer le compartiment dans lequel elle finit sa course ?

Simulation du lancer de n billes

On note

h

 la hauteur de la planche de Galton, 

p

la probabilité qu'une bille a de passer à droite d'un clou à chaque rangée et 

n

le nombre de billes lancées sur la planche.

Compléter la fonction suivante qui renverra une liste detaille \(h+1\) correspondant aux compartiments situés en bas de la planche de Galton. Le nombre à la position 

i

de cette liste correspondra alors au nombre de billes qui sont arrivées dans le compartiment numérotée

i

.

Remarque

Ne pas hésiter à tester le programme en changeant les valeurs de 

h

,

p

et

n

.

def somme(L):
    """renvoie la somme des éléments de la liste L
    Cette fonction pourra être utilisée dans la fonction suivante"""
    total = 0
    for x in L :
        total = total + x
    return total

def galton(h = 10, p = 0.5, n = 1000):
    # A compléter

galton(h = 10, p = 0.5, n = 1000)

On peut alors visualiser le résultat d'une simulation à l'aide du programme suivant (modifier les paramètres pour observer différents cas de figure).
Les ordonnées indiquent les fréquences d'apparition des billes dans les compartiments en abscisse.

from matplotlib.pyplot import *

def visualisation(h, p, n):
    bar(range(h+1),[x/n for x in galton(h, p, n)])
    show()
    close()
visualisation(10, 0.5, 1000)

Comparaison avec la théorie

On considère une planche de Galton de hauteur 

h

. On suppose qu'à chaque rangée, la probabilité pour qu'une bille passe à droite du clou vaut

p

. On note alors \(X\) la variable aléatoire qui donne le compartiment d'arrivée d'une bille à l'issue d'un lancer sur cette planche.

1.Quelle est la loi de la variable aléatoire 

X

? Quelle est alors la probabilité que la billearrive dans le compartiment 

k

à l'issue de sa traversée de la planche ?

2.Compléter la fonction suivante qui permet de calculer

P(X=k)

(utiliser la formule du cours).

Pour calculer la factorielle d'un entier, on utilise la fonction factorial du module math de Python.

from math import factorial

def proba(k, h, p):
    # A compléter

Le programme suivant permet alors de comparer visuellement la distribution empirique et la distribution théorique des billes.

from matplotlib.pyplot import *

def visualisation(h, p, n):
    subplot(1,2,1)
    bar(range(h+1),[x/n for x in galton(h, p, n)])
    title('Distribution empirique')
    subplot(1,2,2)
    bar(range(h+1),[proba(k, h, p) for k in range(h+1)])
    title('Distribution théorique')
    show()
    close()
visualisation(20, 0.7, 1000)

Simulation de la planche de Galton

Présentation et activité complète

Présentation de l'activité

Illustration de la planche de Galton

Simulation de la planche de Galton

Activité CAPYTALE : simulation de la planche de Galton

Rappels nécessaires pour l'activité

Quelques rappels : librairie random

Quelques rappels sur les listes

Exercice d'application

Simulation de la planche de Galton

Simulation d'une variable aléatoire suivant une loi de Bernoulli de paramètre p

Trajectoire d'une bille

Simulation du lancer de n billes

Comparaison avec la théorie