Comparación de gausianas. Distribución normal

http://www.mathworks.es/es/help/stats/multivariate-normal-distribution.html
close all;
mu = [0 0];
Sigma = [1 0; 0 100];
x1 = -10:.2:10; x2 = -10:.2:10;
[X1,X2] = meshgrid(x1,x2);
F = mvnpdf([X1(:) X2(:)],mu,Sigma);
F = reshape(F,length(x2),length(x1));
surf(x1,x2,F);
caxis([min(F(:))-.5*range(F(:)),max(F(:))]);
%axis([-3 3 -3 3 0 .4])
xlabel('x1' ); ylabel('x2'); zlabel('Probability Density');




 Distribución normal bivariante.

Parámetros:

mu = [0 0]; La media me indica donde está centrada
Con una sigma de por ejemplo: 
 
Sigma = [2 1
         1 2];

Obtenemos

 
 

Como se puede observar, se encuentra rotada al ser la autocorrelación en ambos ejes mayor que la correlación xy


Otro ejemplo:


Sigma = [1 0; 0 10]; 
(mayor dispersión en x2 ya que la varianza es mayor que la de x1)
   
 

Caso contrario al anterior

Sigma = [10 0; 0 1]; 


 



Sigma = [2 1; 1 2];
 

Sigma = [10 5; 5 10];