Tag Archives: aléatoire

Marches aléatoires en dimension 1, 2 ou 3, ou histoire éthylique

Revenons à nos histoires d’alcoolique un instant (les précédentes sont en ligne iciet ). L’avant dernière fois (ici), nous avions parlé de la marche aléatoire, et du temps mis avant de trouver la sortie de la rue de la soif.

Et j’avais évoqué la probabilité la probabilité qu’il ne récupère pas les clés (i.e. qu’il ne repasse pas par l’origine) s’il peut sortir un jour de la rue. En probabilité, on appelle cela la probabilité de repasser par l’origine et cela vérifie une propriété étrange dès lors que l’on change de dimension… On va suivre notre alcoolique et noter https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy01.png sa position à la date https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy02.png, puis poser

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy03.png

autrement dit la probabilité de se retrouver à sa position initiale à la date n (et de retrouver ainsi ses clés), le https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png pouvait être en n’importe quelle dimension, on posera

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy05.png

correspondant à la probabilité de premier retour à la position initiale, et enfin

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy06.png

qui est la probabilité de revenir un jour au point de départ, c’est à dire https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png… C’est cette dernière grandeur qui nous intéresse….

  • De la probabilité de récupérer ses clés dans la rue

La rue est ici infinie (contrairement à l’autre fois, où on pouvait en sortir, soit vers la place de Lices, soit vers la place Sainte Anne). Autrement dit, on s’intéresse à la probabilité qu’une marche aléatoire sur https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy08.png issue de https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png repasse un jour par https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png. S’il va à droite ou à gauche avec probabilités https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl99.png et https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl98.png respectivement, notons que https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy11.png, car il faut un nombre de pas pairs pour revenir à l’origine, et

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy12.png

Il est possible de montrer que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy13.png

et en particulier

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy14.png

On peut relier les https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy15.png et les https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy16.png en notant que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy17.png

Or https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy16.png constitue une distribution de probabilité, dont on peut définir la fonction génératrice

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl18.png

On notera que que l’on peut écrire

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy19.png

soit, en utilisant l’expression donnée auparavant

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl19.png

Soit finalement

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl20.png

ce qui se traduit par https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy21.png si la série https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy22.png converge. Aussi,

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl25.png

donc si https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy26.png, notre alcoolique repassera presque sûrement par le point d’origine, alors que si https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl97.png, il ne repasse pas par le bar où il a oublié ses clés avec probabilité https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy27.png|.

  • De la probabilité de récupérer ses clés sur une place

La place est ici infinie de part et d’autre. Autrement dit, on s’intéresse à la probabilité qu’une marche aléatoire sur https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy09.png issue de https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png repasse un jour par https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png. Il a alors 1 chance sur 4 d’aller devant lui, autant d’aller en arrière, d’aller à gauche ou d’aller à droite. On peut noter que là encore
https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy11.png, car il faut un nombre de pas pairs pour revenir à l’origine, et

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl26.png

Or on peut réécrire cette somme en notant que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy28.png

et que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl29.png

aussi,

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl30.png

que l’on peut approcher (via la formule de Stirling) sous le forme

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy31.png

La série https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy22.png est alors divergente, et donc https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl31.png. Aussi, même sur une place, notre alcoolique reviendra presque sûrement à sa place initiale (et ce une infinité de fois, car une fois qu’il y est revenu et qu’il en repart, il repassera presque sûrement à cette position, etc).

  • De la probabilité de récupérer ses clés dans l’espace

Bon, cette fois notre alcoolique passe en dimension 3…. j’admets, il faut de l’imagination. On s’intéresse à la probabilité qu’une marche aléatoire sur https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy10.png issue de https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png repasse un jour par https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethy04.png. Il a alors 1 chance sur 6 d’aller devant lui, autant d’aller en arrière, d’aller à gauche ou d’aller à droite, d’aller en haut ou d’aller en bas.

Il y a désormais 6 voisins, et en raisonnant comme avant (le maillage est un peu plus complexe), on obtient que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl32.png

Je laisse les amateurs de combinatoires faire des calculs ou des approximations, mais on peut montrer que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl33.png

La série \sum u_n est alors convergente, et

https://perso.univ-rennes1.fr/arthur.charpentier/latex/ethyl34.png

qui vaut un peu plus de 1/3. Autrement dit, on a seulement 1 chance sur 3 de revenir à l’origine et de retrouver ses clés…
Bref, en dimension 1 ou 2, on repasse presque sûrement une infinité de fois à l’origine, mais en dimension 3, on a 2 chances sur 3 de ne jamais y repasser…. J’ai toujours trouvé ça étonnant….

Histoire éthylique, et de clés

En tant que statisticien, j’aime beaucoup les alcooliques… ils nous permettent de voir certains problèmes de modélisation sous un autre angle.
Considérons une personne ayant un trousseau de 10 clés, ou disons n pour généraliser. Il souhaite ouvrir la porte de chez lui, mais ne sait plus laquelle est la bonne. Combien d’essais – en moyenne – faudra-t-il avant de tomber sur la bonne clé ?
Pour une personne à jeun (et saine d’esprit), le plus simple est de tenter une clé au hasard, puis, si ça ne marche pas, de prendre celle à sa droite (ou à sa gauche, peut importe, mais il faut une stratégie et s’y tenir, à moins de pouvoir noter quelle clé a été utilisée). Bref, au bout de n essais, on est certain de tomber sur la bonne clé. Si le fait de manière un peu formelle (histoire de justifier une intuition), et que l’on note https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-01.png l’évènement “la bonne clé est tirée au bout de https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-02.png essais“, alors

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-03.png

Raisonnons alors simplement,

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-04.png

puis

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-05.png

puis

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-06.png

etc. Aussi, on peut légitimement penser que

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-07.png

Moralité (il s’agit de la moyenne d’une loi uniforme)

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-08.png

Pour une personne ivre (ou un peu moins à jeun), à chaque essai, elle prend une clé au hasard dans le trousseau (on pourrait appeler ça la stratégie du poisson rouge). Le nombre d’essais peut être potentiellement infini. Et en reprenant les notations précédentes, on note que si

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-09.png

en revanche

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-10.png

puisqu’à l’étape 2, on n’a pas écarté la clé qui ne marche pas… Aussi,

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-11.png

qui correspond à une loi géométrique, ou loi de Pascal, de paramètre https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-12.png. Pour ceux qui ne savent plus comment calculer l’espérance d’une telle loi, on peut prendre trois lignes pour refaire le calcul

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-13.png
soit tout simplement

https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-14.png

i.e.
https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-15.png

Aussi https://perso.univ-rennes1.fr/arthur.charpentier/latex/etyl-16.png.
Grosso modo, il faut deux fois plus de temps pour trouver une clé quand on est rond comme une queue d’pelle, que quand on est à jeun (sans compter le temps où l’on parle à son trousseau de clé, en espérant qu’il trouvera un moyen de désigner la bonne clé). La prochaine partie du billet parlera de “marche aléatoire“, et de difficultés pour quelqu’un d’ivre mort à quitter a rue Saint Michel….