在 MATLAB 中绘制给定 MVN 的多元正态 PDF 的轮廓？

Question

我有一个二元高斯，我定义如下：

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.

现在我想绘制上述矩阵的 pdf 的轮廓。我做了什么：

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

现在我想绘制二元高斯 beta 的等高线。我知道我应该使用命令轮廓，但这个要求 Z 是一个方阵。我该如何解决？我很困惑，不知道如何绘制双变量高斯的轮廓！！非常感谢任何帮助..

谢谢

Answer 1

您需要定义x、y 轴并使用meshgrid（或ndgrid）以两个矩阵X 和Y 的形式生成x、y 值的所有组合.然后，您计算 X 和 Y 的 Z 值（您的高斯 pdf），并使用 contour（等高线图）绘制 Z 作为 X、Y 的函数surf（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