Le cours
Résumé en cinq points
- Toutes les simulations partent de la loi uniforme sur . Bernoulli s'obtient par le test , la loi géométrique en comptant les épreuves jusqu'au premier succès.
- Si est continue et strictement croissante, suit la loi de fonction de répartition . Pour l'exponentielle, . Pour la normale, on inverse numériquement.
- Dans le cas discret, on renvoie la première valeur dont la probabilité cumulée dépasse . La loi de Poisson s'obtient ainsi, ou comme limite de lois binomiales de paramètres et .
- La moyenne empirique estime l'espérance avec une erreur en . La fonction de répartition empirique estime la fonction de répartition.
- Le théorème central limite se voit sur l'histogramme de , qui se rapproche de la densité gaussienne quand grandit. Une chaîne de Markov se simule pas à pas avec la ligne de la matrice de transition.
Pourquoi ce chapitre
Le hasard est partout en biologie :
- mutations ;
- transmission des allèles ;
- rencontres entre individus ;
- erreurs de mesure.
Beaucoup de modèles probabilistes n'ont pas de résultat exact calculable, ou le calcul est long. La simulation offre une autre voie. On fait tirer au hasard par l'ordinateur un grand nombre de réalisations du modèle, puis on observe les fréquences, les moyennes, les histogrammes.
Ce chapitre apprend à :
- fabriquer toutes les lois usuelles à partir d'une seule, la loi uniforme ;
- estimer une espérance ou une loi à partir de simulations ;
- voir de ses yeux des résultats théoriques comme la loi des grands nombres ou le théorème central limite.
Ces compétences servent ensuite pour les tests statistiques et dans de nombreux sujets de modélisation en sciences de la vie.
Classiques et HP
Méthodes et approfondissements du chapitre · le statut est précisé pour BCPST2.
Classiques et HP réservés
Prolongements et démonstrations de Informatique.
Fin de l'aperçu gratuit
Encore 15 parties de cours et 10 corrigés.
Méthodes et exercices
Exercices classiques corrigés
Les exercices sont rédigés pour ce document. Les mentions de niveau indiquent une difficulté, sans attribution à une épreuve. Les fonctions des modules utilisées ci-dessous sont supposées documentées.
Rangés par thème, du plus simple au plus exigeant. Difficulté et durée indicatives.
Simulation de lois discrètes 3
Exercice 1. Attendre une germination ★
Application du cours · 10 min · Niveau khôlle
Chaque graine germe indépendamment avec probabilité , où .
Programmer le rang de la première germination et étudier le coût.
Voir l'indice
Technique compter des essais jusqu’au premier succès.
Voir le corrigé offert
- Après j échecs, k vaut j + 1, le rang du prochain essai. Un succès fait renvoyer k, un échec l’incrémente.
- On utilise une nouvelle uniforme à chaque essai.
import random
import math
def premier_succes(p):
k = 1
while random.random() >= p:
k = k + 1
return k
- La probabilité de dépasser m essais vaut (1 − p)m, qui tend vers zéro. L’arrêt est presque sûr.
- Le coût moyen est proportionnel à 1/p, la mémoire est constante. Pour 0 < p < 1, aucun coût maximal fini n’existe. Pour p = 1, un seul tirage suffit.
piège : Invoquer un variant entier décroissant. Ici, la terminaison est probabiliste.
Exercice 2. Trois issues possibles ★
Application du cours · 10 min · Niveau khôlle
Une observation produit 0, 2 ou 5 colonies, avec probabilités respectives 0,2, 0,5 et 0,3. Programmer un tirage et proposer deux contrôles.
Voir le début du corrigé
- Les intervalles sont [0, 0,2[, [0,2, 0,7[ et [0,7, 1[. …
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Exercice 4. Poisson par cumulées ★★
Classiques de khôlle et concours · 20 min · Niveau écrit
Pour , programmer l’inversion d’une loi de Poisson sans calculer de factorielles.
Prévoir un signal d’échec si l’accumulation numérique cesse de progresser.
Voir le début du corrigé
- On initialise p à e−λ. Au rang k, p représente P(X = k) et c représente P(X ≤ k).
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Simulation de lois continues 3
Exercice 3. Première défaillance ★★
Classiques de khôlle et concours · 15 min · Niveau khôlle
Deux composants ont des durées indépendantes, exponentielles de paramètres 2 et 3 par heure. Simuler la première défaillance et déterminer sa durée moyenne.
Voir le début du corrigé
- Pour t ≥ 0, les deux composants fonctionnent encore avec probabilité e−2te−3t = e−5t.
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Exercice 7. Estimer une durée plafonnée ★★
Classiques de khôlle et concours · 20 min · Niveau écrit
suit une loi exponentielle de paramètre .
Pour , estimer par simulations indépendantes, puis obtenir une valeur exacte de contrôle.
Voir le début du corrigé
- On transforme chaque durée en sa valeur plafonnée. …
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Exercice 8. Inverser la normale ★★★
Approfondissement · 25 min · Niveau écrit, approfondissement
Une fonction phi évalue la répartition normale centrée réduite.
Pour , on fournit avec Approcher le quantile par étapes de dichotomie.
Voir le début du corrigé
- L’unique solution reste dans [a, b]. …
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Approximation de lois 2
Exercice 5. Approcher les événements rares ★★
Classiques de khôlle et concours · 15 min · Niveau écrit
Approcher une Poisson de paramètre par une binomiale.
Comparer les choix et , puis programmer le tirage.
Voir le début du corrigé
- On impose n entier positif et n ≥ λ pour que λ/n soit une probabilité. …
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Exercice 9. Voir le théorème central limite ★★
Classiques de khôlle et concours · 20 min · Niveau écrit
Avec des exponentielles indépendantes de paramètre 1, construire une somme centrée réduite. Décrire une expérience distinguant les rôles de n et N.
Voir le début du corrigé
- Chaque variable a moyenne 1 et variance 1. …
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Répartition empirique 1
Exercice 6. Répartition avec répétitions ★
Application du cours · 15 min · Niveau khôlle
Pour l’échantillon [1, 1, 3, 5], déterminer la répartition empirique, puis programmer son évaluation sur une liste non vide.
Voir le début du corrigé
- La fonction vaut 0 avant 1, 1/2 sur [1, 3[, 3/4 sur [3, 5[, puis 1 à partir de 5.
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Chaîne à deux états 1
Exercice 10. Canal à deux états ★★★
Approfondissement · 25 min · Niveau écrit, approfondissement
Un canal fermé s’ouvre avec probabilité par pas.
Un canal ouvert se ferme avec probabilité .
Il est initialement fermé.
Programmer son état après pas et calculer sa probabilité d’ouverture.
Voir le début du corrigé
- L’état courant vaut 0 ou 1. Chaque tour réalise exactement une transition avec une nouvelle uniforme.
La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.
Approfondissement hors programme 1
Exercice HP. Estimer e−1 par tirages antithétiques ★★★
Approfondissement hors programme
Approfondissement facultatif · hors programme en BCPST2
Soit de loi uniforme sur et . On veut estimer . Soit indépendantes de même loi que . On pose et .
- Calculer et vérifier que et sont sans biais.
- Calculer et en déduire .
- Calculer . Sans la calculer, justifier le signe de .
- En déduire , puis le rapport (valeur approchée).
- On veut un écart-type d'au plus . Combien d'évaluations de faut-il pour chaque méthode ? On donne .
Corrigé détaillé dans Sup et Spé+ : essai gratuit 7 jours.
Fiche de cours rédigée à partir du programme officiel d'informatique de BCPST (2021), programme des semestres 3 et 4, partie 3.1 « Méthodes numériques et statistiques ». Elle couvre les thèmes suivants du tableau :