正态总体参数的假设检验matlab处理

正态总体参数的检验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.05Tail：尾部类型变量,‘both’双侧检验（默认），u不等于uo；‘right’右侧检验，u>u0; ‘left’左侧检验,u<u0；返回值：h：假设的结果（0,1），h=0时，接受假设H0；h=1,拒绝假设H0p：检验的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; %原假设中的mu0sigma=2; %总体标准差Alpha=0.05; %显著性水平%调用ztest函数做总体均值的双侧检验（默认）,%返回变量h，检验的p值，均值的置信区间muci，检验统计量的观测值zval[h,p,muci,zval]=ztest(x,mu0,sigma,Alpha)h =1p =0.0282muci =100.1212 102.1455zval =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; %原假设中的mu0sigma=2; %总体标准差Alpha=0.05; %显著性水平tail='right';%调用ztest函数做总体均值的双侧检验（默认）,%返回变量h，检验的p值，均值的置信区间muci，检验统计量的观测值zval[h,p,muci,zval]=ztest(x,mu0,sigma,Alpha,tail)h =1p =0.0141muci =100.2839 Infzval =2.1947返回值中，h=1且p=0.0141<0.05，在显著性水平下=0.05下拒绝了原假设H0：u<=u0=100;如果：H0：u>uo=100; H1=u<u0那么tail应为‘left’；返回值将会h=0，p>0.05，接受原假设H0，即认为u>100;2 总体标准差未知时的单个正态总体均值的t检验例：化肥厂用自动包装机包装化肥，某日测得9包化肥的质量如下：49.4 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9设每包化肥的质量服从正态分布，是否可以认为每包化肥的平均质量为50？取显著性水平a=0.05。

分析：这是总体标准差未知时的单个正态总体均值的检验，根据题目要求可以写出如下假设：H0：u=u0=50，H1：u/=u0(u不等于u0) MATLAB统计工具箱中提供了ttest函数用来做总体标准差未知时的正态总体均值的检验，调用格式和ztest类似，返回值有点不同[h,p,muci,stats]=ttest(x,mu0,Alpha,Tail)输入参数中没有标准差，其它都一样返回值stats是一个结构体变量，包括t检验统计量的观测值，自由度，和样本的标准差；其它都一样%定义样本观测值向量x=[49.9 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9];mu0=50; %原假设中的均值u0=50Alpha=0.05; %显著性水平alpha%调用ttest函数做总体均值的双侧检验%返回变量h，检验值p，均值的置信空间muci，结构体变量stats[h,p,muci,stats]=ttest(x,mu0,Alpha)h =p =1.0000muci =49.0625 50.9375stats =tstat: -1.7478e-14df: 8sd: 1.2196由于返回值h=0，p=1>0.05,所以在显著性水平=0.05下接受原假设H0：u=u0=50，认为每包化肥的平均质量为50，并且总体均值u的置信水平为95%的置信区间为[49.0625 50.9375]3 总体标准差未知时的两个正态总体均值的比较t 检验（1）两独立样本的t 检验例：甲、乙两台机床加工同一种产品，从这两台机床加工的产品中随机抽取若干件，测得产品直径为：甲机床：20.1,20.0,19.3,20.6,20.2,19.9,20.0,19.9,19.1,19.9乙机床：18.6, 19.1,20.0,20.0,20.0,19.7,19.9,19.6,20.2设甲、乙两机床加工的产品的直径分布服从正态分布N（u1，a1^2）和N(u2,a2^2)，试比较甲、乙两台机床加工的产品的直径均值是否有显著性差异，取显著性水平a=0.05 分析：这是总体标准差未知，并且两样本是相互独立的，对这两样本均值做比较检验，根据题目要求，可写出如下假设H0：u1=u2 H1:u1/=u2(u1不等于u2)MATLAB统计工具箱中的ttest2函数可以用来做总体标准差未知时的两个正态总体均值的比较检验；调用格式：[h,p,muci,stats]=ttest(x,y,Alpha,Tail,vartype)x，y为输入的两个样本观测值Alpha为显著性水平tail为尾部类型cartype：为方差类型，用来指定两总体方差是否相等，‘equal’表示等方差，‘unequal’表示异方差返回值与ttest函数一致，muci是指均值差的置信空间%定义甲机床的样本观测值向量x=[20.1,20.0,19.3,20.6,20.2,19.9,20.0,19.9,19.1,19.9];%定义乙机床的样本观测值向量y=[18.6, 19.1,20.0,20.0,20.0,19.7,19.9,19.6,20.2];Alpha=0.05; %显著性水平tail='both'; %尾部类型为双侧vartype='equal'; %方差类型为等方差%调用ttest2函数作两个正态总体均值的比较检验%返回变量h，检验的p值，均值差的置信区间，结构体变量stats[h,p,muci,stats]=ttest2(x,y,Alpha,tail,vartype)h =p =0.3191muci =-0.2346 0.6791stats =tstat: 1.0263df: 17sd: 0.4713返回的检验值p>0.05,所以在显著性水平=0.05下，接受原假设H0：u1=u2，认为甲、乙两台机床加工的产品的直径没有显著差异。

此时，u1-u2的置信水平为95%的置信区间为[-0.2346 0.6791]（2）配对样本的t 检验（两样本不是独立的）例：两组（各10名）有资质的评酒员分别对12种不同的酒进行品评，每个评酒员在品尝后进行评分，然后对每组的每个样品计算其平均分，评分结果如下样品1 样本2 样品3 样品4 样品5 样品6 样品7 样品8 样品9 样品10 样品11 样品12第一组80.3 68.6 72.2 71.5 72.3 70.174.6 73.0 58.7 78.6 85.6 78.0第二组74.0 71.2 66.3 65.3 66.0 61.668.8 72.6 65.7 72.6 77.1 71.5设两组评酒员的评分分布服从正态分布N(u1,a1^2)和N(u2,a2^2)，试比较两组评酒员的评分是否有显著差异，取显著性水平a=0.05分析：由于每个红酒样本都对应两个评分，显然样本等长，并且两样本不独立，这是配对样本的比较问题，根据题目要求可写出如下的假设：H0：u1=u2，H1：u1/=u2(u1不等于u2) 由于两个样本不独立，通常的做法是将两个样本对应数据最差，把两个正态总体均值的比较检验转化为单个正态总体均值的检验，然后就可用ttest函数进行检验上面的假设改写为如下假设H0：u=u1-u2=0，H1：u/=0(u不等于0) 然后调用ttest函作配对样本的比较t 检验调用格式[h,p,muci,stats]=ttest(x,y,Alpha,Tail)x，y为输入的观测样本观测值向量，其它参数与ttest一致%样本1x=[80.3,68.6,72.2,71.5,72.3,70.1,74.6,73.0,58.7,78.6,85.6, 78.0];%样本2。