6 MATLAB参数估计与假设-正态总体参数的检验

正态总体参数的检验

更多MATLAB数据分析视频请点击，或者在网易云课堂上搜索《MATLAB数据分析与统计》 http://study.163.com/course/courseMain.htm?courseId=1003615016

  1 总体标准差已知时的单个正态总体均值的U检验

  例：

 某切割机正常工作时，切割的金属棒的长度服从正态分布N(100,4)。从该切割机切割的一批金属棒中随机抽取15根，测得长度为：

 97 102 105 112 99 103 102 94 100 95 105 98 102 100 103

假设总体的方差不变，试检验该切割机工作是否正常，即检验总体均值是否等于100？，取显著性水平a=0.05。

分析：

这是总体标准差已知时的单个正态总体均值的检验，根据题目要求可写出如下假设：

  H0：u=u0=100，                H1=u /=u0(u不等于u0)

 H0称为原假设，H1称为被择假设（或对立假设）

MATLAB统计工具箱中的ztest函数用来做总体标准差已知时的单个正态总体均值的检验

调用格式ztest

[h,p,muci,zval]=ztest(x,mu0,Sigma,Alpha,Tail)

x:是输入的观测向量

mu0：假设的均值

Sigma：总体标准差

Alpha：显著性水平，默认0.05

Tail：尾部类型变量,‘both’双侧检验（默认），u不等于uo；‘right’右侧检验，u>u0; ‘left’左侧检验,u<u0；

返回值：

h：假设的结果（0,1），h=0时，接受假设H0；h=1,拒绝假设H0

p：检验的p值，p>Alpha时，接受原假设H0；p<=Alpha时，拒绝原假设H0.

muci：总体均值u的置信水平为1-Alpha的置信区间

zval：检验统计量的观测值

%定义样本观测值向量

x=[97 102 105 112 99 103 102 94 100 95 105 98 102 100 103];

mu0=100;       %原假设中的mu0

sigma=2;       %总体标准差

Alpha=0.05;    %显著性水平

%调用ztest函数做总体均值的双侧检验（默认）,

%返回变量h，检验的p值，均值的置信区间muci，检验统计量的观测值zval

[h,p,muci,zval]=ztest(x,mu0,sigma,Alpha)

h =

     1

p =

    0.0282

muci =

  100.1212  102.1455

zval =

    2.1947

由ztest函数返回值可以看到，h=1，且p=0.0282<0.05，所以在显著性水平=0.05下拒绝的原假设H0:u=u0=100，因此认为该切割机不能正常工作，同时还返回了总体均值的置信水平为95%（1-0.05）的置信区间为[100.1212  102.1455]。

现在我们已经知道u不等于u0（100），那么接下来还需要做如下的检验

H0：u<=u0=100;               H1:u>u0

这里就需要设置一下尾部类型变量了，Tail用来指定备择假设H1的形式，它可能取的字符串‘both’，‘right’和'left'，对应的备择假设分别为H1：u/=u0(不等于)，（双侧检验）、

H1：u>u0,(右尾检验)和H1：u<u0(左尾检验)

%定义样本观测值向量

x=[97 102 105 112 99 103 102 94 100 95 105 98 102 100 103];

mu0=100;       %原假设中的mu0

sigma=2;       %总体标准差

Alpha=0.05;    %显著性水平

tail='right';

%调用ztest函数做总体均值的双侧检验（默认）,

%返回变量h，检验的p值，均值的置信区间muci，检验统计量的观测值zval

[h,p,muci,zval]=ztest(x,mu0,sigma,Alpha,tail)

h =

     1

p =

    0.0141

muci =

  100.2839       Inf

zval =

    2.1947

返回值中，h=1且p=0.0141<0.05，在显著性水平下=0.05下拒绝了原假设H0：u<=u0=100;

如果：

H0：u>uo=100;      H1=u<u0

那么tail应为‘lef