Aller au contenu
Télécharger l'application

Accueil › BCPST › Informatique

BCPST2 · Informatique · semaine 5

Simulation de lois et convergences

Simulation d'une loi géométrique, d'une loi de Poisson, exponentielle et normale, méthode de la réciproque de la fonction de répartition. Fonction de répartition empirique, illustration du théorème central limite, chaînes de Markov.

Sommaire

Le cours

Résumé en cinq points

  1. Toutes les simulations partent de la loi uniforme sur [0,1][0, 1]. Bernoulli s'obtient par le test U<pU < p, la loi géométrique en comptant les épreuves jusqu'au premier succès.
  2. Si FF est continue et strictement croissante, F−1(U)F^{-1}(U) suit la loi de fonction de répartition FF. Pour l'exponentielle, X=−ln⁡(1−U)/λX = -\ln(1 - U)/\lambda. Pour la normale, on inverse Φ\Phi numériquement.
  3. Dans le cas discret, on renvoie la première valeur dont la probabilité cumulée dépasse UU. La loi de Poisson s'obtient ainsi, ou comme limite de lois binomiales de paramètres nn et λ/n\lambda/n.
  4. La moyenne empirique estime l'espérance avec une erreur en σ/N\sigma/\sqrt{N}. La fonction de répartition empirique estime la fonction de répartition.
  5. Le théorème central limite se voit sur l'histogramme de (Sn−nμ)/(σn)(S_n - n\mu)/(\sigma\sqrt{n}), qui se rapproche de la densité gaussienne quand nn 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.

Essayer 7 jours gratuitement Sans engagement · cours, exercices et fiches complets

Fin de l'aperçu gratuit

Encore 15 parties de cours et 10 corrigés.

Essayer 7 jours gratuitement Sans engagement · cours, exercices et fiches complets

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é pp, où 0<p≤10 < p \leq 1.

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
  1. Après j échecs, k vaut j + 1, le rang du prochain essai. Un succès fait renvoyer k, un échec l’incrémente.
  2. 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
  1. La probabilité de dépasser m essais vaut (1 − p)m, qui tend vers zéro. L’arrêt est presque sûr.
  2. 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é
  1. 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 λ>0\lambda > 0, 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é
  1. 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é
  1. 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

XX suit une loi exponentielle de paramètre λ>0\lambda > 0.

Pour T>0T > 0, estimer E(min⁡(X,T))E(\min(X, T)) par NN simulations indépendantes, puis obtenir une valeur exacte de contrôle.

Voir le début du corrigé
  1. 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 0<u<10 < u < 1, on fournit a<ba < b avec phi(a)≤u≤phi(b).\mathrm{phi}(a) \leq u \leq \mathrm{phi}(b). Approcher le quantile par mm étapes de dichotomie.

Voir le début du corrigé
  1. 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 44 par une binomiale.

Comparer les choix n=20n = 20 et n=400n = 400, puis programmer le tirage.

Voir le début du corrigé
  1. 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é
  1. 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é
  1. 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é a=0,2a = 0{,}2 par pas.

Un canal ouvert se ferme avec probabilité b=0,3b = 0{,}3.

Il est initialement fermé.

Programmer son état après nn pas et calculer sa probabilité d’ouverture.

Voir le début du corrigé
  1. 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 UU de loi uniforme sur ]0,1[]0,1[ et f(x)=exf(x)=e^x. On veut estimer I=∫01ex dxI=\int_0^1 e^x\,\dd x. Soit U1,…,U2nU_1,\dots,U_{2n} indépendantes de même loi que UU. On pose N2n=12n∑k=12neUkN_{2n}=\frac1{2n}\sum_{k=1}^{2n}e^{U_k} et An=1n∑k=1neUk+e1−Uk2A_n=\frac1n\sum_{k=1}^n\frac{e^{U_k}+e^{1-U_k}}2.

  1. Calculer II et vérifier que N2nN_{2n} et AnA_n sont sans biais.
  2. Calculer σ2=Var⁡(eU)\sigma^2=\operatorname{Var}(e^U) et en déduire Var⁡(N2n)\operatorname{Var}(N_{2n}).
  3. Calculer c=Cov⁡(eU,e1−U)c=\operatorname{Cov}(e^U,e^{1-U}). Sans la calculer, justifier le signe de cc.
  4. En déduire Var⁡(An)\operatorname{Var}(A_n), puis le rapport Var⁡(N2n)/Var⁡(An)\operatorname{Var}(N_{2n})/\operatorname{Var}(A_n) (valeur approchée).
  5. On veut un écart-type d'au plus 10−310^{-3}. Combien d'évaluations de ff faut-il pour chaque méthode ? On donne e≈2,718282e\approx2{,}718282.

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 :