数学建模常用各种检验方法及常用方法

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学建模各种检验方法

1.单个总体2

Nμσ的均值μ的检验:

(,)

2

σ已知,关于均值的检验用ztest命令来实现.

[h,p,ci]=ztest(x,mu,sigma,alpha,tail)

2

σ已知,关于均值的检验用ttest命令来实现.

[h,p,ci]=ttest(x,mu,alpha,tail)

2.两个正态总体均值差的检验(t 检验)

还可以用t 检验法检验具有相同方差的2 个正态总体均值差的假设。在Matlab 中

由函数ttest2 实现,命令为:

[h,p,ci]=ttest2(x,y,alpha,tail)

3.分布拟合检验

在实际问题中,有时不能预知总体服从什么类型的分布,这时就需要根据样本来检

验关于分布的假设。下面介绍2χ检验法和专用于检验分布是否为正态的“偏峰、峰度

检验法”。

2

χ检验法

0 H :总体x的分布函数为F(x) ,

1 H : 总体x的分布函数不是F(x).

在用下述χ 2检验法检验假设0 H 时,若在假设0 H 下F(x)的形式已知,但其参数

值未知,这时需要先用极大似然估计法估计参数,然后作检验。

偏度、峰度检验

4.其它非参数检验

Wilcoxon秩和检验

在Matlab中,秩和检验由函数ranksum实现。命令为:

[p,h]=ranksum(x,y,alpha)

其中x,y可为不等长向量,alpha为给定的显著水平,它必须为0和1之间的数量。p返回

产生两独立样本的总体是否相同的显著性概率,h返回假设检验的结果。如果x和y的总

体差别不显著,则h为零;如果x和y的总体差别显著,则h为1。如果p 接近于零,则可对

原假设质疑。

5.中位数检验

在假设检验中还有一种检验方法为中位数检验,在一般的教学中不一定介绍,但在

实际中也是被广泛应用到的。在Matlab中提供了这种检验的函数。函数的使用方法简单,

下面只给出函数介绍。

signrank函数

signrank Wilcoxon符号秩检验

[p,h]=signrank(x,y,alpha)

其中p给出两个配对样本x和y的中位数相等的假设的显著性概率。向量x,y的长度必须

相同,alpha为给出的显著性水平,取值为0和1之间的数。h返回假设检验的结果。如果

这两个样本的中位数之差几乎为0,则h=0;若有显著差异,则h=1。signtest函数

signtest 符号检验

[p,h]= signtest(x,y,alpha)

其中p给出两个配对样本x和y的中位数相等的假设的显著性概率。x

和y若为向量,二者

的长度必须相同;y亦可为标量,在此情况下,计算x的中位数与常数y之间的差异。alpha

和h同上。

matlab 判断正态分布

总体分布正态性检验

进行参数估计和假设检验时,通常总是假定总体服从正态分布,虽然在许多情况下这个假定是合理的,但是当要以此为前提进行重要的参数估计或假设检验,或者人们对它有较大怀疑的时候,就确有必要对这个假设进行检验,

进行总体正态性检验的方法有很多种,以下针对MATLAB统计工具箱中提供的程序,简单介绍几种方法。

1)Jarque-Bera检验

利用正态分布的偏度g1和峰度g2,构造一个包含g1,g2的分布统计量(自由度n=2),对于显著性水平,当分布统计量小于分布的分位数时,接受H 0:总体服从正态分布;否则拒绝H0,即总体不服从正态分布。这个检验适用于大样本,当样本容量n较小时需慎用。Matlab命令:h =jbtest(x),[h,p,jbstat,cv] =jbtest(x,alpha)。

2)Kolmogorov-Smirnov检验

通过样本的经验分布函数与给定分布函数的比较,推断该样本是否来自给定分布函数的总体。容量n的样本的经验分布函数记为F n(x),可由样本中小于x的数据所占的比例得到,给定分布函数记为G(x),构造的统计量为,即两个分布函数之差的最大值,对于假设H0:总体服从给定的分布G(x),及给定的,根据D n的极限分布(n®¥时的分布)确定统计量关于是否接受H0的数量界限。

因为这个检验需要给定G(x),所以当用于正态性检验时只能做标准

正态检验,即H0:总体服从标准正态分布。Matlab命令:h =kstest(x)。

3)Lilliefors检验

它将Kolmogorov-Smirnov检验改进用于一般的正态性检验,即H0:总体服从正态分布,其中由样本均值和方差估计。Matlab 命令:

h =lillietest(x),[h,p,lstat,cv]=lillietest(x,alpha)。

4)另外还有一种方法:首先对于数据进行标准化:Z = ZSCORE(X),然后在进行2)的Kolmogorov-Smirnov检验,检验是否为标准正态分布,类似于对于方法2)的改进

相关文档
最新文档