小编Jur*_*iën的帖子

我有2个数组:一个有x坐标,另一个有y坐标.由于蒙特卡罗模拟,两者都是正态分布.我知道如何找到两个阵列的sigma和mu,并获得95%的置信区间:

[mu,sigma]=normfit(x_array);
hist(x_array);
x=norminv([0.025 0.975],mu,sigma)

但是,两个阵列都是相互关联的.为了绘制组合阵列的概率分布,我使用多元正态分布.在MATLAB中,这给了我:

[MuX,SigmaX]=normfit(x_array);
[MuY,SigmaY]=normfit(y_array);
mu = [MuX MuY];
Sigma=cov(x_array,y_array);
x1 = MuX-4*SigmaX:5:MuX+4*SigmaX; x2 = MuY-4*SigmaY:5:MuY+4*SigmaY;
[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(:))]);
set(gca,'Ydir','reverse')
xlabel('x0-as'); ylabel('y0-as'); zlabel('Probability Density');

到现在为止还挺好.现在我想计算95%的概率区域.我正在寻找一个功能mndinv,就像norminv.然而,MATLAB中不存在这样的功能,这是有道理的,因为有无穷无尽的可能性......有人有关于如何获得95%概率区域的提示吗？提前致谢.