【问题标题】:Standard deviation ellipse for a bivariate normal distribution双变量正态分布的标准差椭圆
【发布时间】:2016-03-03 04:38:32
【问题描述】:

如何绘制仅显示从 1-σ 到 7-σ 的七个等值线的二元正态分布的标准差等值线(使用 MATLAB)?

mu=[2 3];     %mean
Sigma=[1 1];  %standard deviation
x1 = -10:.2:10;   
x2 = -10:.2:10;
[X1,X2] = meshgrid(x1,x2);
F=mvnpdf([X1(:) X2(:)],mu,Sigma);   %compute Gaussian pdf
F=reshape(F,length(x2),length(x1));
contour(x1,x2,F,'ShowText','on');  %Equally spaced contour

我不需要等间距的轮廓,而是只显示从 1σ 到 7σ 的七个轮廓

https://en.wikipedia.org/wiki/Multivariate_normal_distribution 此链接中图 1 概率密度函数中的绿线(3-sigmas 椭圆)显示了我想要做什么,但我不知道他们是如何做到的。如何为不同的 sigma 绘制绿线 (3-σ) 和类似的线。

【问题讨论】:

标签: matlab contour normal-distribution


【解决方案1】:

我想这段代码可以完成这项工作。它基于tolerance intervals

mu = [2 3];
Sigma = [1 1];
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));

%Values drawn from a normal distribution are 
%5-σ,4-σ and 3-σ awayfrom the mean
TL=1-erf([5 4 3]/sqrt(2));             
[C,h]=contour(x1,x2,F,TL,'ShowText','on');

但是,我不确定是否可以将此轮廓的标签更改为 5-σ、4-σ 和 3-σ 而不是数字?

【讨论】:

    最近更新 更多