Tracer le contour de la normale multivariée PDF de un MVN dans MATLAB?

J'ai un bivariée gaussien je définis comme suit:

I=[1 0;0 1];
mu=[0,0];
sigma=0.5*I;
beta = mvnrnd(mu,sigma,100);  %100x2 matrix where each column vector is a variable.

maintenant, je veux tracer un contour de la pdf de la matrice. Ce que j'ai fait:

Z = mvnpdf(beta,mu,sigma);  %100x1 pdf matrix

Maintenant, je veux tracer un contour de la gaussienne bidimensionnelle bêta. Je sais que je dois utiliser la commande contour mais celui-ci nécessite Z une matrice carrée. comment puis-je résoudre ce problème? Je suis très confus et ne sais pas comment tracer le contour de la gaussienne bidimensionnelle!! TOUTE AIDE EST GRANDEMENT APPRÉCIÉE..

Merci

OriginalL'auteur Evan | 2013-11-24

  1. 15

    Vous devez définir votre x, y axes et l'utilisation meshgrid (ou ndgrid) pour générer toutes les combinaisons de valeurs x, y, sous la forme de deux matrices X et Y. Ensuite, vous calculez la Z valeurs (votre Gaussien pdf) pour ceux X et Y, et l'intrigue Z en fonction de X , Y à l'aide de contour (tracé de contour), ou peut-être surf (graphique 3D).

    mu = [0,0]; %//data
sigma = [.5 0; 0 .5]; %//data
x = -5:.1:5; %//x axis
y = -4:.1:4; %//y axis

[X Y] = meshgrid(x,y); %//all combinations of x, y
Z = mvnpdf([X(:) Y(:)],mu,sigma); %//compute Gaussian pdf
Z = reshape(Z,size(X)); %//put into same size as X, Y
%//contour(X,Y,Z), axis equal  %//contour plot; set same scale for x and y...
surf(X,Y,Z) %//... or 3D plot

    Tracer le contour de la normale multivariée PDF de un MVN dans MATLAB?
    Tracer le contour de la normale multivariée PDF de un MVN dans MATLAB?

    Merci beaucoup a parfaitement fonctionné
    Ce code fonctionne aussi parfaitement en Octave, à condition que le paquet de statistiques est chargé pkg load statistics 🙂

    OriginalL'auteur Luis Mendo