Je réserve ici
Maths expertes (option Tle) Terminale

Ressources · Terminale · Maths expertes (option Tle)

Graphes probabilistes et chaînes de Markov

Rappels sur l’arbre pondéré, les probabilités totales et les suites arithmético-géométriques ; chaîne de Markov à deux ou trois états, graphe pondéré et matrice de transition ; distribution πₙ₊₁ = πₙP et πₙ = π₀Pⁿ, coefficients de Pⁿ ; distribution invariante ; convergence à long terme et seuil.

À propos de cette page
Cinq problèmes corrigés : la fidélisation des abonnés d’une plateforme avec un seuil en Python, les parts de marché de deux enseignes avec le calcul de Pⁿ par récurrence, la répartition d’une population entre trois zones, la disponibilité d’un serveur et le modèle de diffusion d’Ehrenfest, dont la distribution invariante existe sans que la chaîne converge.

5 problèmes, du plus simple au plus difficile. Écris tes réponses, puis « Vérifier ».

Résoudre un problème : la méthode

  1. Je traduis l’énoncé : états et ordre des états, matrice de transition P (chaque ligne somme à 1), distribution initiale π₀ en ligne.
  2. Je calcule les premières distributions par πₙ₊₁ = πₙP pour me faire une idée, et je contrôle que leurs coefficients somment à 1.
  3. Je cherche la distribution invariante : πP = π avec la condition « somme égale à 1 ».
  4. Je prouve le comportement à long terme : à deux états, suite arithmético-géométrique et limite de qⁿ ; à trois états, résultat admis dont je vérifie l’hypothèse.
  5. J’interprète dans la situation : proportion limite, seuil trouvé à la calculatrice ou avec une boucle « tant que », en vérifiant les deux rangs qui encadrent la réponse.
Problème 1 · Fidélisation des abonnés d’une plateforme de lecture★★★Moyen
Une plateforme de lecture en ligne suit ses clients mois après mois. Chaque client est soit abonné (état A), soit non abonné (état N). D’un mois au suivant, un abonné garde son abonnement avec la probabilité 0,88 ; un non-abonné souscrit un abonnement avec la probabilité 0,05. Au mois 0, 20% des clients sont abonnés. On note an la proportion d’abonnés au bout de n mois et πn=(an1−an), les états étant rangés dans l’ordre (A;N).
  1. 1. Écrire la matrice de transition P et la distribution initiale π0. Coefficient p22 :
    Un abonné reste abonné avec la probabilité 0,88, donc se désabonne avec la probabilité 1−0,88=0,12. Un non-abonné souscrit avec la probabilité 0,05 et reste non abonné avec la probabilité 0,95. D’où P=(0,880,120,050,95) et π0=(0,20,8).
  2. 2. Calculer π1 et π2. Proportion d’abonnés au bout de deux mois, arrondie à 10−4 :
    π1=π0P : 0,2×0,88+0,8×0,05 =0,176+0,04 =0,216, donc π1=(0,2160,784). Puis 0,216×0,88+0,784×0,05 =0,19008+0,0392 =0,22928 : π2=(0,229280,77072), soit environ 22,93% d’abonnés.
  3. 3. À l’aide d’un arbre pondéré entre le mois n et le mois n+1, justifier que an+1=0,83an+0,05. Contrôle : a3 arrondi à 10−4 =
    Rédaction attendue. Au premier niveau de l’arbre, le client est abonné au mois n (probabilité an) ou non (probabilité 1−an) ; au second niveau, il est abonné au mois n+1 avec la probabilité 0,88 dans le premier cas, 0,05 dans le second. Ces deux cas forment une partition ; d’après la formule des probabilités totales, an+1=0,88an+0,05(1−an) =0,83an+0,05. Alors a3=0,83×0,22928+0,05 =0,2403024, soit 0,2403.
  4. 4. Valeur exacte du premier coefficient de la distribution invariante : (fraction irréductible)
    π=(x1−x) est invariante si et seulement si 0,83x+0,05=x, soit 0,17x=0,05 et x=517. La distribution invariante est (5171217) ≈(0,29410,7059).
  5. 5. On pose un=an−517. Démontrer que (un) est géométrique, puis exprimer an en fonction de n et donner sa limite. Premier terme : u0= (fraction irréductible)
    Rédaction attendue. Comme 517 est le point fixe, 0,05=517−0,83×517. Pour tout n, un+1=0,83an+0,05−517 =0,83an−0,83×517 =0,83un. La suite (un) est géométrique de raison 0,83, de premier terme u0=0,2−517 =3,4−517 =−885. Donc an=517−885×0,83n. Comme −1<0,83<1, liman=517≈0,294 : la proportion d’abonnés se stabilise vers 29,4%.
  6. 6. Recopier et compléter la fonction ci-dessous pour qu’elle renvoie le nombre de mois nécessaires pour que la proportion d’abonnés atteigne au moins s.
    def seuil(s):
        a = ...
        n = 0
        while ... :
            a = ...
            n = n + 1
        return n
    Valeur renvoyée par seuil(0.28) :
    def seuil(s):
        a = 0.2
        n = 0
        while a < s:
            a = 0.83 * a + 0.05
            n = n + 1
        return n
    Par le calcul : an⩾0,28 ⟺885×0,83n⩽517−0,28 =1,285 ⟺0,83n⩽0,15. Or 0,8310≈0,1552>0,15 et 0,8311≈0,1288<0,15 ; la suite (0,83n) est décroissante, donc le premier rang est 11 (a10≈0,2795 et a11≈0,2820). La fonction renvoie 11.
  7. 7. La direction vise 35% d’abonnés. Justifier que c’est impossible avec ces probabilités, puis donner le taux de souscription mensuel (à la place de 0,05) qui donnerait une proportion limite de 35%, le taux de désabonnement 0,12 restant inchangé, arrondi à 10−4 :
    Rédaction attendue. u0<0 et 0<0,83<1, donc un=u0×0,83n est négatif et croissant : (an) est croissante et an<517 ≈0,294 <0,35 pour tout n. L’objectif est hors d’atteinte. Avec un taux de souscription β, la proportion limite est β0,12+β (cours, deux états). β0,12+β =0,35⟺β =0,042+0,35β ⟺0,65β =0,042⟺β =21325 ≈0,0646 : il faudrait près de 6,5% de souscriptions par mois.
Problème 2 · Parts de marché de deux enseignes★★★Moyen
Deux enseignes de vente de pièces détachées, X et Y, se partagent tout un marché. D’une semaine à la suivante, un client de X passe chez Y avec la probabilité 0,15, et un client de Y passe chez X avec la probabilité 0,25. La semaine 0, chaque enseigne détient la moitié du marché. On note xn la part de marché de X au bout de n semaines et πn=(xn1−xn), les états étant rangés dans l’ordre (X;Y).
  1. 1. Écrire P, puis calculer π1 et π2 : x2=
    Un client de X y reste avec la probabilité 0,85, un client de Y y reste avec la probabilité 0,75 : P=(0,850,150,250,75) et π0=(0,50,5). π1=(0,550,45) (car 0,5×0,85+0,5×0,25 =0,55), puis 0,55×0,85+0,45×0,25 =0,4675+0,1125 =0,58 : π2=(0,580,42).
  2. 2. Premier coefficient de la distribution invariante :
    0,85x+0,25(1−x) =x⟺0,6x+0,25 =x⟺0,4x =0,25⟺x =0,625. Distribution invariante : (0,6250,375).
  3. 3. Établir la relation entre xn+1 et xn, puis exprimer xn en fonction de n à l’aide d’une suite auxiliaire. Contrôle : x3=
    Rédaction attendue. Le premier coefficient de πn+1=πnP donne xn+1=0,85xn+0,25(1−xn) =0,6xn+0,25. On pose vn=xn−0,625 : vn+1=0,6xn+0,25−0,625 =0,6xn−0,375 =0,6(xn−0,625) =0,6vn. La suite (vn) est géométrique de raison 0,6 et de premier terme v0=0,5−0,625=−0,125, donc xn=0,625−0,125×0,6n. Contrôle : x3=0,625−0,125×0,216 =0,625−0,027 =0,598.
  4. 4. Au bout de combien de semaines la part de marché de X dépasse-t-elle 62% ?
    xn>0,62⟺0,125×0,6n <0,005⟺0,6n <0,04. Or 0,66≈0,0467>0,04 et 0,67≈0,0280<0,04 ; la suite (0,6n) est décroissante, donc il faut 7 semaines (x6≈0,6192 et x7≈0,6215).
  5. 5. On note L la matrice carrée dont les deux lignes sont égales à la distribution invariante, et B=I2−L. Vérifier que LP=L et BP=0,6B, puis démontrer par récurrence que Pn=L+0,6nB pour tout n∈ℕ. En déduire la limite de Pn et l’interpréter. Contrôle : (P2)21=
    Rédaction attendue. L=(0,6250,3750,6250,375) et B=I2−L=(0,375−0,375−0,6250,625). Chaque ligne de L est la distribution invariante, donc LP=L. Première ligne de BP : 0,375×0,85−0,375×0,25 =0,225 et 0,375×0,15−0,375×0,75 =−0,225, soit 0,6×(0,375−0,375) ; de même la seconde ligne vaut (−0,3750,375) =0,6×(−0,6250,625). Donc BP=0,6B. Initialisation. L+B=I2=P0. Hérédité. Si Pn=L+0,6nB, alors Pn+1=PnP =LP+0,6nBP =L+0,6n+1B. Conclusion. Pn=L+0,6nB pour tout n. Comme −1<0,6<1, Pn→L : les deux lignes de la limite valent la distribution invariante. À long terme, un client achète chez X avec la probabilité 0,625, quelle que soit l’enseigne de départ. Contrôle : (P2)21=0,625−0,36×0,625 =0,4, et directement 0,25×0,85+0,75×0,25 =0,4.
Problème 3 · Répartition d’une population entre trois zones★★★Difficile
Une agence d’urbanisme étudie, année après année, la répartition des habitants d’une agglomération entre le centre (C), la première couronne (K) et la zone rurale (R). Le tableau donne les probabilités de changement de zone d’une année sur l’autre (en ligne : zone de départ ; en colonne : zone d’arrivée).
vers Cvers Kvers R
depuis C0,700,150,15
depuis K0,050,900,05
depuis R0,150,150,70
Au début de l’étude, 50% des habitants vivent au centre, 20% en première couronne et 30% en zone rurale. On note πn=(cnknrn) la distribution au bout de n années. On admet le résultat du cours à trois états.
  1. 1. Écrire P et π0, puis calculer π1 : c1=
    P=(0,700,150,150,050,900,050,150,150,70) et π0=(0,50,20,3). c1=0,5×0,70+0,2×0,05+0,3×0,15 =0,35+0,01+0,045 =0,405 ; k1=0,075+0,18+0,045=0,3 ; r1=0,075+0,01+0,21=0,295.
  2. 2. Montrer que kn+1=0,75kn+0,15 pour tout n. Contrôle : k2=
    Rédaction attendue. Les événements « être en C », « être en K », « être en R » l’année n forment une partition ; d’après la formule des probabilités totales, kn+1=0,15cn+0,90kn+0,15rn =0,15(cn+rn)+0,90kn. Or cn+kn+rn=1, donc cn+rn=1−kn et kn+1=0,15(1−kn)+0,90kn =0,75kn+0,15. Contrôle : k2=0,75×0,3+0,15=0,375.
  3. 3. Exprimer kn en fonction de n, puis donner limkn=
    Point fixe : ℓ=0,75ℓ+0,15, soit ℓ=0,6. La suite (kn−0,6) est géométrique de raison 0,75 et de premier terme 0,2−0,6=−0,4, donc kn=0,6−0,4×0,75n. Comme −1<0,75<1, limkn=0,6.
  4. 4. Déterminer la distribution invariante (xyz). Part du centre : x=
    D’après la question 2 appliquée à une distribution invariante, y=0,75y+0,15, donc y=0,6 et x+z=0,4. La première équation donne 0,70x+0,05×0,6+0,15z=x, soit 0,30x=0,03+0,15(0,4−x), donc 0,45x=0,09, x=0,2 et z=0,2. Vérification sur le troisième coefficient : 0,2×0,15+0,6×0,05+0,2×0,70 =0,03+0,03+0,14 =0,2.
  5. 5. Justifier que (πn) converge, puis donner le nombre d’années au bout duquel la part de la première couronne dépasse 59% :
    Rédaction attendue. Tous les coefficients de P sont strictement positifs : d’après le résultat admis, (πn) converge vers l’unique distribution invariante (0,20,60,2), quelle que soit π0 (la question 3 le confirme sur kn). Puis kn>0,59⟺0,4×0,75n <0,01⟺0,75n <0,025. Or 0,7512≈0,0317 et 0,7513≈0,0238 ; la suite (0,75n) décroît, donc il faut 13 années (k12≈0,5873 et k13≈0,5905).
  6. 6. L’agence affirme : « le centre et la zone rurale finiront par se vider ». Discuter cette affirmation. Pour une agglomération de 450000 habitants, nombre d’habitants de la zone rurale à long terme :
    Rédaction attendue. L’affirmation est fausse : (cn) et (rn) convergent vers 0,2. À long terme, le centre et la zone rurale gardent chacun 20% de la population : le modèle prévoit une stabilisation, où les départs de chaque zone sont compensés par les arrivées. Zone rurale : 0,2×450000=90000 habitants.
Problème 4 · Disponibilité d’un serveur★★★Difficile
À chaque heure, un serveur est soit disponible (état D), soit indisponible (état I). S’il est disponible, il tombe en panne au cours de l’heure suivante avec une probabilité a, où a∈]0;1[. S’il est indisponible, il est réparé au cours de l’heure suivante avec la probabilité 0,6. Les états sont rangés dans l’ordre (D;I) et dn désigne la probabilité que le serveur soit disponible au bout de n heures. À l’instant initial, il est disponible.
  1. 1. Écrire P en fonction de a, puis π0. Pour a=0,08, coefficient p11 :
    P=(1−aa0,60,4) et π0=(10) ; p21=0,6 est la probabilité d’être réparé en une heure. Pour a=0,08 : p11=1−0,08=0,92, probabilité de rester disponible une heure de plus.
  2. 2. Montrer que l’unique distribution invariante est (0,6a+0,6aa+0,6). Contrôle pour a=0,1 : disponibilité limite = (fraction irréductible)
    Rédaction attendue. Soit π=(d1−d). πP=π⟺(1−a)d+0,6(1−d) =d⟺0,6 =(a+0,6)d. Comme a+0,6>0, d=0,6a+0,6 et 1−d=aa+0,6 : la distribution invariante existe et elle est unique. Pour a=0,1 : 0,60,7=67.
  3. 3. L’hébergeur s’engage sur une disponibilité limite de 96%. Valeur de a correspondante :
    0,6a+0,6 =0,96⟺0,6 =0,96a+0,576 ⟺0,96a =0,024⟺a =0,025. Alors P=(0,9750,0250,60,4).
  4. 4. Désormais, a prend la valeur trouvée à la question précédente. Montrer que dn+1=0,375dn+0,6. Contrôle : d2=
    dn+1=0,975dn+0,6(1−dn) =0,375dn+0,6. Avec d0=1 : d1=0,975 et d2=0,375×0,975+0,6 =0,365625+0,6 =0,965625.
  5. 5. Exprimer dn, étudier le sens de variation de (dn) et donner le plus petit entier n tel que dn<0,9601 :
    Rédaction attendue. Point fixe : ℓ=0,375ℓ+0,6, soit 0,625ℓ=0,6 et ℓ=0,96. La suite (dn−0,96) est géométrique de raison 0,375 et de premier terme 0,04, donc dn=0,96+0,04×0,375n. Le terme 0,04×0,375n est strictement positif et décroît : (dn) est strictement décroissante, de limite 0,96, et reste au-dessus de 0,96. Puis dn<0,9601 ⟺0,375n <0,0025 ; or 0,3756≈0,00278 et 0,3757≈0,00104 : n=7. Au bout de 7 heures, la probabilité de disponibilité est à moins de 10−4 de sa limite : le départ sur un serveur disponible n’a presque plus d’influence.
  6. 6. En régime stationnaire, durée moyenne d’indisponibilité sur une année de 8760 heures, en heures :
    La probabilité d’indisponibilité vaut 1−0,96=0,04 ; sur 8760 heures, la durée moyenne est 0,04×8760=350,4 heures, soit un peu plus de 14 jours par an.
Problème 5 · Le modèle de diffusion d’Ehrenfest avec deux boules★★★Difficile
Les physiciens Paul et Tatiana Ehrenfest ont proposé en 1907 un modèle de diffusion d’un gaz entre deux récipients. On l’étudie avec deux boules numérotées 1 et 2, réparties entre deux urnes A et B. À chaque étape, on choisit au hasard un numéro, 1 ou 2 avec la même probabilité, et l’on change d’urne la boule qui porte ce numéro. On note Xn le nombre de boules dans l’urne A après n étapes : les états sont 0, 1 et 2, rangés dans cet ordre. Au départ, les deux boules sont dans l’urne A : X0=2. On admet le résultat du cours à trois états.
  1. 1. Justifier que (Xn) est une chaîne de Markov et écrire P. Probabilité de passer de l’état 1 à l’état 2 :
    Rédaction attendue. Le nombre de boules dans A à l’étape n+1 ne dépend que de ce nombre à l’étape n et du numéro tiré, indépendant du passé : c’est une chaîne de Markov. Depuis 0 (les deux boules en B), la boule tirée passe en A : on va en 1 avec la probabilité 1. Depuis 1, on tire la boule de A (on va en 0) ou celle de B (on va en 2), chacune avec la probabilité 0,5. Depuis 2, on va en 1. Donc P=(0100,500,5010).
  2. 2. Calculer π1, π2, π3 et π4. p(X4=2)=
    π0=(001), π1=π0P =(010), π2=(0,500,5), π3=(010) et π4=(0,500,5). Donc p(X4=2)=0,5.
  3. 3. Calculer P2, puis démontrer que P3=P et que P2k+1=P pour tout k∈ℕ. Coefficient de P2 en ligne de l’état 1, colonne de l’état 1 :
    Rédaction attendue. P2=(0,500,50100,500,5) : partant de 1, on est sûr d’être revenu en 1 deux étapes plus tard. Puis P3=P2P : la ligne (0,500,5) multipliée par P donne (010) et la ligne (010) donne (0,500,5) : P3=P. Récurrence : P1=P ; si P2k+1=P, alors P2k+3=P2k+1P2=P3=P. Donc P2k+1=P pour tout k, et P2k=P2k−1P=P2 pour k⩾1.
  4. 4. Déterminer la distribution invariante. Probabilité de l’état 0 :
    π=(xyz) : πP=π s’écrit 0,5y=x, x+z=y et 0,5y=z. Avec x+y+z=1 : 2y=1, donc y=0,5 et x=z=0,25. C’est la loi binomiale de paramètres 2 et 0,5 : chaque boule est dans A avec la probabilité 0,5.
  5. 5. La suite (πn) converge-t-elle ? Justifier. p(X100=0)=
    Rédaction attendue. Pour k⩾1, π2k=π0P2k =π0P2=(0,500,5) et π2k+1=π0P =(010). La suite (πn) alterne entre ces deux distributions : elle ne converge pas, bien qu’une distribution invariante existe. Le résultat du cours ne s’applique pas : aucune puissance de P n’a tous ses coefficients strictement positifs. 100 est pair : p(X100=0)=0,5.
  6. 6. Variante : à chaque étape, on lance d’abord une pièce équilibrée ; sur « pile », rien ne bouge ; sur « face », on applique la règle précédente. La matrice devient Q=12(P+I3). Vérifier que la distribution invariante trouvée plus haut l’est encore pour Q, et que Q2 a tous ses coefficients strictement positifs. Coefficient de Q2 en ligne de l’état 0, colonne de l’état 2 :
    Rédaction attendue. Q=(0,50,500,250,50,2500,50,5). Pour π=(0,250,50,25) : πQ=12(πP+π) =12(π+π) =π. Puis Q2=(0,3750,50,1250,250,50,250,1250,50,375), par exemple (Q2)02=0,5×0+0,5×0,25+0×0,5 =0,125. Tous ses coefficients sont strictement positifs : d’après le résultat admis, (πn) converge vers (0,250,50,25), quelle que soit la répartition de départ.
Bloqué sur ce chapitre ?

Moi c'est Ben, j'ai créé ce site pour aider le plus d'élèves possible, partout en France. Je peux aussi reprendre avec toi ce qui coince en cours particulier : en visio, ou à domicile si tu es à Marseille.

Prof de maths à Marseille · Cours particuliers au lycée · Aide aux devoirs