Discussion :

[gegon]

Bonjour,

j'aimerai definir un point qui a une probabilité de sortir selon une gausienne.
Y-a-t-il sous matlab, une fonction connue pour definir ceci facilement?

Merci d'avance

[shaiHulud]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
doc randn
doc normrnd
doc mvnrnd

[gegon]

Bonjour,

en fait, je cherche à réaliser un tirage aléatoire dans une loi normale sur 2dimensions.

J'y arrive avec une dimension :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

norme= normrnd (0,sigma)

mais quelqu'un aurait une idée sur 2 dimensions?

merci d'avance

[FLB]

Salut,
L'utilisation des paramètres 3 et 4 de la fonction normrnd ne te conviennent pas?

R = normrnd(mu,sigma,m,n,...) or R = normrnd(mu,sigma,[m,n,...]) generates an m-by-n-by-... array. The mu, sigma parameters can each be scalars or arrays of the same size as R

[gegon]

peut-etre mais je ne dois pas bien les comprendre, je peux creer une matrice de 2 colonnes où chacune des colonnes se referent à un sigma particulier?

[FLB]

une loi normale sur 2dimensions

Tu veux plutôt un tableau à 2 dimensions de taille [m,n], plutot que un tableau de taille [2xn], non?
Si tu veux le même mu et sigma :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
mu = 1;
sigma = 0.5;
R = normrnd(mu,sigma,2,3);

Si tu veux des mus différents :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

normrnd([1 1 1;2 2 2],0.5,2,3);

Tu obtiens ainsi un tirage selon une distribution normale, avec mu = 1 sur la première ligne, mu = 2 sur la seconde.

[gegon]

oui, il est vrai.

Je réflechi en meme temps et en fait il me faut une valeur et une seule qui soit définie par une loi normale bidimensionelle de meme moyenne mais avec 2 sigmas differents.

En essayant de visualiser, cela donne un point qui peut prendre les valeurs d'une ellipse.

j'espere être comprehénsible..

[FLB]

Ah ok!
Dans ce cas il faut utiliser la fonction mvnrndque VV33D te conseille.
En reprenant l'exemple de la doc :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
mu = [1 1];
SIGMA = [1 0; 0 3]; % Décorrélation entre les 2 dimensions : ellipse selon les axes X Y. Applique une rotation matricielle pour faire tourner l'ellipse.
r = mvnrnd(mu,SIGMA,100);%100 pour visualiser comment sont dispersés tes points, tu remplaces 100 par 1 pour ne demander qu'une valeur
plot(r(:,1),r(:,2),'+');