Salut je ne suis pas sûr si mon algorithme est correct, je suis en train de reproduire la fonction mvnrnd de matlab mais en gsl. J'ai trouvé un algorithme dans quelques articles de journal qui produit un vecteur de normale multivariée, mais j'ai besoin d'une matrice de nombres aléatoires normaux multivariésmatlab mvnrnd en gsl
disons que la distribution est Z ~ (mu, sigma);
en supposant que sigma est une matrice déjà définie positive.
un algorithme je l'ai trouvé sur le web dit
1. cholskey(sigma) = A
2. generate uniform gaussian vector r
3. matrix vector triangular product with gsl_blas_dtrmv A * r
4. add mu to Ar and that will be a vector of multivariate normal random numbers
ma méthode ci-dessous
sont les modifications suivantes belowcorrect au produit une matrice de variables aléatoires
1. cholskey(sigma) = A
2. generate uniform gaussian matrix R
3. matrix matrix scalar product AR
4. add mu to AR and that will be a matrix of multivariate normal random numbers
merci! Je ne pensais pas que j'aurais quelqu'un célèbre comme toi pour répondre à ma question – pyCthon