IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

C Discussion :

Equivalent de randn (de matlab) en c


Sujet :

C

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre averti
    Inscrit en
    Janvier 2007
    Messages
    29
    Détails du profil
    Informations forums :
    Inscription : Janvier 2007
    Messages : 29
    Par défaut Equivalent de randn (de matlab) en c
    Bonjour à tous,
    j'aimerais savoir s'il existe une fonction équivalente de la fonction randn de matalab sous en langage c . car je cherche à générer des variables aléatoire de moyenne 0 et d'écart type =1 par exemple. merci

  2. #2
    Expert confirmé

    Profil pro
    Inscrit en
    Janvier 2007
    Messages
    10 610
    Détails du profil
    Informations personnelles :
    Âge : 67
    Localisation : France

    Informations forums :
    Inscription : Janvier 2007
    Messages : 10 610
    Billets dans le blog
    2
    Par défaut
    voir la fonction rand() et srand() pour intialiser.

    Liens ici pour les man pages : rand et srand.

    Lien ici pour une explication de l'utilisation :
    http://c.developpez.com/faq/?page=no..._random_bornes

  3. #3
    Membre éprouvé
    Profil pro
    Inscrit en
    Décembre 2004
    Messages
    1 299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 299
    Par défaut
    Salut, oui c'est possible. Il faut que tu utilises l'algo de Box-Muller

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
     
    #define TWOPI (6.2831853071795864769252867665590057683943387987502) /* 2 * pi */
     
    /* 
       RAND is a macro which returns a pseudo-random numbers from a uniform
       distribution on the interval [0 1]
    */
    #define RAND (rand())/((double) RAND_MAX)
     
    /* 
       RANDN is a macro which returns a pseudo-random numbers from a normal
       distribution with mean zero and standard deviation one. This macro uses Box
       Muller's algorithm
    */
    #define RANDN (sqrt(-2.0*log(RAND))*cos(TWOPI*RAND))

  4. #4
    Membre Expert
    Inscrit en
    Décembre 2004
    Messages
    1 478
    Détails du profil
    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 478
    Par défaut
    Citation Envoyé par salseropom Voir le message
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
     
    #define RANDN (sqrt(-2.0*log(RAND))*cos(TWOPI*RAND))
    RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?

  5. #5
    Expert confirmé

    Inscrit en
    Août 2006
    Messages
    3 967
    Détails du profil
    Informations forums :
    Inscription : Août 2006
    Messages : 3 967
    Par défaut
    Fia,
    Citation Envoyé par DaZumba Voir le message
    RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?
    +1

    Et en plus, on en prend le log, alors qu'il peut renvoyer 0.

    Et la valeur pour TWOPI dépasse largement le nombre de digits valides pour un double. Je sais, le compilateur va tronquer, mais ça peut donner une fausse impression de précision.

  6. #6
    Membre éprouvé
    Profil pro
    Inscrit en
    Décembre 2004
    Messages
    1 299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 299
    Par défaut
    Citation Envoyé par droggo Voir le message
    Fia,

    +1

    Et en plus, on en prend le log, alors qu'il peut renvoyer 0.

    Et la valeur pour TWOPI dépasse largement le nombre de digits valides pour un double. Je sais, le compilateur va tronquer, mais ça peut donner une fausse impression de précision.
    Certes le log peut renvoyer 0 mais si U suit une loi uniforme sur [0, 1], la probabilité que U soit nul est égale à 0...

    TWOPI dépasse largement le nombre de digits. Je le sais, mais je n'avais pas envie de me prendre la tête entre un compilo 64 bits, 32 bits, PC, LINUX, MAC. Donc sur Maple, j'ai fait un Digits=50 pour une question de "sécurité"

  7. #7
    Expert confirmé

    Inscrit en
    Novembre 2005
    Messages
    5 145
    Détails du profil
    Informations forums :
    Inscription : Novembre 2005
    Messages : 5 145
    Par défaut
    Citation Envoyé par DaZumba Voir le message
    RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?
    Oui. Si on veut, pour chaque paire de nombres uniformement distribues retourner une paire de nombres distribues normalement en utilisant aussi sin a la place de cos.

  8. #8
    Membre éprouvé
    Profil pro
    Inscrit en
    Décembre 2004
    Messages
    1 299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 299
    Par défaut
    Citation Envoyé par Jean-Marc.Bourguet Voir le message
    Oui. Si on veut, pour chaque paire de nombres uniformement distribues retourner une paire de nombres distribues normalement en utilisant aussi sin a la place de cos.
    .

    Tout à fait, c'est l'algo de Box Muller. J'ai opté arbitrairement pour le sinus

  9. #9
    Membre éprouvé
    Profil pro
    Inscrit en
    Décembre 2004
    Messages
    1 299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 299
    Par défaut
    Citation Envoyé par DaZumba Voir le message
    RAND est evalue deux fois et va sortir deux valeurs differentes. C'est normal ?
    oui c'est normal. Dans l'algorithme de Box Muller, il faut tirer 2 nombres aléatoires suivant la loi uniforme sur [0, 1]. Ces 2 nombres doivent être INDEPENDANTS !! C'est pourquoi il y a 2 fois RAND.

Discussions similaires

  1. La fonction randn de matlab en C
    Par magol1 dans le forum C
    Réponses: 2
    Dernier message: 07/06/2012, 11h05
  2. Equivalent Java de nchoosek() MatLab
    Par Tehko dans le forum Général Java
    Réponses: 5
    Dernier message: 28/05/2012, 17h40
  3. Réponses: 6
    Dernier message: 02/03/2008, 02h47
  4. Réponses: 2
    Dernier message: 02/09/2005, 15h23
  5. Equivalent de Maple ou Matlab en open source ou gratuit
    Par wikers dans le forum Autres Logiciels
    Réponses: 2
    Dernier message: 02/09/2005, 15h23

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo