利用MATLAB绘制置信区域

最新推荐文章于 2024-08-14 20:15:48 发布

weixin_30433075

最新推荐文章于 2024-08-14 20:15:48 发布

阅读量9.8k

点赞数 5

文章标签： matlab 人工智能

原文链接：http://www.cnblogs.com/caizhao/p/8081604.html

版权

本文介绍了如何使用MATLAB的自编函数ConfidenceRegion来绘制一维、二维和三维的置信区域，包括置信区间、置信椭圆和置信椭球。通过实例展示了在正态分布和经验分布下绘制95%置信水平的区域，并提供了详细的代码解释。

摘要由CSDN通过智能技术生成

<MATLAB小技巧之二十四：利用MATLAB绘制置信区域>

***************************************

统计中经常会遇到求置信区间、置信区域（如置信椭圆、置信椭球）等，有时候需要把置信区域画出来，这样看起看更为直观，下面结合具体案例介绍调用自编函数ConfidenceRegion绘制置信区域。

【例1】绘制置信区间。产生一元正态分布随机数向量，绘制样本数据的95%置信区间。

x = normrnd(10,4,100,1);
subplot(1,2,1)
ConfidenceRegion(x)
subplot(1,2,2)
ConfidenceRegion(x,'exp')

效果图如下图所示：

【例2】绘制置信椭圆。产生二元正态分布随机数矩阵，绘制样本数据的95%置信椭圆区域。

x = mvnrnd([1;2],[1 4;4 25],100);
subplot(1,2,1)
ConfidenceRegion(x)
subplot(1,2,2)
ConfidenceRegion(x,'exp')

果图如下图所示：

【例3】绘制置信椭球。产生三元正态分布随机数矩阵，绘制样本数据的95%置信椭球区域。
x = mvnrnd([1;2;3],[3 0 0;0 5 -1;0 -1 1],100);
subplot(1,2,1)
ConfidenceRegion(x)
subplot(1,2,2)
ConfidenceRegion(x,'exp')

效果图如下图所示：

自编函数ConfidenceRegion程序代码如下：

function ConfidenceRegion(xdat,alpha,distribution)
%   绘制置信区域（区间、椭圆区域或椭球区域）
%   ConfidenceRegion(xdat,alpha,distribution)
%   xdat：样本观测值矩阵,p*N 或 N*p的矩阵，p = 1,2 或 3
%   alpha：显著性水平，标量，取值介于[0,1]，默认值为0.05
%   distribution：字符串（'norm'或'experie