Répondre à la discussion
Affichage des résultats 1 à 4 sur 4

[Maple] Marche au hasard/ Marche de l'ivrogne.



  1. #1
    IexGIfate.

    [Maple] Marche au hasard/ Marche de l'ivrogne.


    ------

    Bonjour à tous,

    Je solicite votre aide pour mon tipe, pour la partie expérimentale, j'aimerai concevoir un programme maple qui peut simuler une marche de l'ivrogne (en gros on part d'un point sur l'origine, dans le plan XY, dont son ordonnée sera toujours nulle, et on ajoute aléatoirement 1 ou en enlève 1 à son abscisse) et une marche aléatoire (même principe dans le plan, y différent de 0).

    Bon j'ai déja du mal avec ca, mais je voudrais ploter tout ca pour retrouver quelques résultats théoriques remarquables (loi poissonienne, gaussienne, et que toutes les lois tendent vers la loi normale)

    Le principe serait je pense de représenter sur un grand nombre de tirages la répartition des points, et de montrer que la probabilité de se trouver à une distance de plus en plus grande de l'origine décroit exponentiellement... retrouver également une expression de la moyenne, variance, écart type, "expérimentalement" que je pourrais comparer avec les résultats théoriques...

    Voila, merci encore!

    {Je suis en Mpsi, il s'agit d'un tipe de physique sur la physique statistique}


    IexGIfate.

    -----

  2. #2
    HAL 9000

    Re : [Maple] Marche au hasard/ Marche de l'ivrogne.

    Quelle est ta question du coup ?

    "et on ajoute aléatoirement 1 ou en enlève 1 à son abscisse"

    Avec quelle probabilité du coup ?

  3. #3
    IexGIfate.

    Re : [Maple] Marche au hasard/ Marche de l'ivrogne.

    Huum, il s'agit de deux programmes différents, pour celui de la marche de l'ivrogne, la probabilité d'avoir x(n+1)=x(n)+1 est de 0,5 et celle d'avoir x(n+1)=x(n)-1 est aussi de 0,5...

  4. #4
    Seirios

    Re : [Maple] Marche au hasard/ Marche de l'ivrogne.

    Bonjour,

    Est-ce que ça correspond à ce que tu cherches ? :

    Code:
    marche:=proc(n)
    > local j,p,l;
    > l:=[[0,0]];
    > for j from 1 to n do
    > p:=rand(0..1)();
    > if p=0 then
    > l:=[op(l),[j,l[j][2]+1]];
    > fi;
    > if p=1 then 
    > l:=[op(l),[j,l[j][2]-1]];
    > fi;
    > od;
    > plot(l);
    > end;
    ;
    If your method does not solve the problem, change the problem.

Sur le même thème :

Discussions similaires

  1. [Thermique] Chauffe eau marche que sur marche forcée
    Par rosette28 dans le forum Dépannage
    Réponses: 15
    Dernier message: 21/07/2009, 15h45
  2. Windows7 : quoi qui marche, quoi qui marche pas?
    Par archeos dans le forum Logiciel - Software - Open Source
    Réponses: 11
    Dernier message: 18/01/2009, 16h04
  3. Maple qui ne marche pas
    Par ABN84 dans le forum Logiciel - Software - Open Source
    Réponses: 6
    Dernier message: 17/03/2008, 22h26
  4. tv qui marche pas et qui marche!!
    Par bobluc dans le forum Dépannage
    Réponses: 12
    Dernier message: 19/09/2006, 18h12
  5. Réponses: 5
    Dernier message: 11/08/2006, 21h29