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

Probabilités Discussion :

[Proba] Générer un échantillon gaussien "oblique" ~> factorisation de Cholesky


Sujet :

Probabilités

Mode arborescent

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Invité
    Invité(e)
    Par défaut [Proba] Générer un échantillon gaussien "oblique" ~> factorisation de Cholesky
    Bonjour,

    Pour m'exercer, j'essaye de générer un échantillon gaussien 2D de cette forme :



    Ce que j'ai à disposition, c'est une fonction prédéfinie randn() qui me donne un nombre aléatoire selon la loi gaussienne d'espérance zéro et de variance un.

    Pour l'instant, j'arrive à créer un échantillon gaussien "vertical" ou "horizontal" de cette façon :

    Code pseudocode : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    coordonnée_x = randn() * écart_type_x + moyenne_x
    coordonnée_y = randn() * écart_type_y + moyenne_y

    Si "écart_type_x" est plus grand que "écart_type_y" alors l'échantillon a une forme horizontale :



    Mais je ne sais pas comment introduire une covariance différente de zéro entre l'axe des x et l'axe des y.

    Merci d'avance pour vos indications.
    Dernière modification par Invité ; 16/05/2011 à 12h45.

+ Répondre à la discussion
Cette discussion est résolue.

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