Dans ce TP on se propose de simuler une loi exponentielle avec Python.

La loi exponentielle est utilisée dans le cadre de durée de vie. On fabrique une ampoule et on se demande combien de temps elle va tenir. La durée de vie de l'ampoule est aléatoire, c'est donc une variable aléatoire et on peut s'interroger sur sa densité de probabilité.

Pour bien comprendre le fonctionnement de la loi exponentielle, on va passer par un jeu.

Dans ce jeu, à chaque tour, on jette un dé (ou une pièce, ou un osselet, ou ce qu'on veut…)

Par exemple on pourrait dire qu'à chaque tour on jette le dé et que si on fait 6, on a perdu. On dirait alors qu'on la probabilité $p = \frac{1}{6}$ de perdre à chaque tour.

De façon plus générale, on supposera qu'on a la probabilité $p$ de perdre à chaque tour – c'est à dire $p$ pas forcément égal à $\frac{1}{6}$.

On s'intéresse au nombre de tours où on a survécu.

On réalise d'abord une fonction qui simule le dé (ou tirage de cartes, ou autre chose…)

On veut écrire une fonction perdu(p) qui renvoie vrai avec la probabilité p.

from random import random

def perdu(p):
    '''
    renvoie True avec la probabilité p
    '''
    if p > random():
        return True
    else:
        return False

Entrez ce code et vérifiez que perdu(0.5) renvoie True à peu près une fois sur deux.

Une partie est une succession de lancers de dés. Chaque fois qu'on gagne, on compte +1. Quand c'est perdu, c'est fini.

Ajoutez à la suite :

def simulation(p):
    '''
    simule une partie et renvoie le nombre de lancer gagnants
    '''
    n = 0
    while ...
        n = ...
    return n

On sait qu'à chaque tour, on a la probabilité $p$ de perdre.

en voir le fonctionnement, nous allons