正态总体参数区间估计的MATLAB实现

正态总体参数区间估计的MATLAB实现陈少云【摘要】本文介绍了MATLAB软件的normfit()函数在求解正态总体参数的区间估计中的长处和短处，结合实例编写了MATLAB程序求解标准差o已知时均值p的置信区间和均值p已知时标准差o的置信区间,弥补了normfit()函数在该方面的不足.【期刊名称】《吉林师范大学学报(自然科学版)》【年(卷),期】2010(031)001【总页»】3页(P76-78)【关键词】正态总体;均值;标准差;置信区间【作者】陈少云【作者单位】四川建筑职业技术学院,计算机系,四川,德阳,618000【正文语种】中文【中图分类】O212.2总体参数的点估计作为待估参数的近似值给出了明确的数量描述,在统计分析中有多方面的应用.但点估计没有给出这种近似的精确程度和可信程度,使其在实际应用中受到很大的限制,区间估计却可以弥补这一不足.在工程技术中广泛使用的数学软件MATLAB提供了现成的函数normfit()可以很方便地求出正态总体标准差o未知时均值p的置信区间和均值p未知时标准差o的置信区间,与数理统计公式和查相关的临界值表计算的结果完全吻合.其函数调用格式为:[muhat,sigmahat,muci,sigmaci]二normfit(x,alpha).其中muhat,sigmahat分别为正态分布的参数p和。的估计值,muci,sigmaci分别为它们的置信区间，置信度为(1-alpha)x100%,alpha为显著性水平.例1商店用机器包装某种商品，每包重量X服从正态分布，为检查包装的质量，对机器包装的商品抽测8包,其重量为试估计机器包装的商品重量的均值p和标准差o的置信度为0.95的置信区间解在MATLAB的命令窗口输入如下语句结果显示为结果表明均值p的估计值为5.0100,其置信度为95%的置信区间为(4.9500,5.0700),这与标准差o未知时运用数理统计公式和查T分布临界值表计算的结果一致;标准差o的估计值为0.0717,其置信度为95%的置信区间为(0.0474,0.1460),这与均值p未知时运用数理统计公式和查卡方分布临界值表计算的结果一致.但在计算标准差o已知时均值p的置信区间和均值p已知时标准差o的置信区间方面有较大的误差•下面是结合数理统计知识设计的MATLAB程序求解标准差o已知时均值p的置信区间和均值p已知时标准差o的置信区间的两个例子.由数理统计知识,标准差o已知时正态总体均值的区间估计应采用U统计量，置信度为1-a的p置信区间为•其中.例2某课程命题初衷,其成绩Z~N(p,13.52),考毕抽查其中10份试卷的成绩为：试求该课程平均成绩p的置信区间.(置信度1-a=0.95)解在MATLAB的编辑窗口建立如下的M-文件(并保存为myfun1.m),以便以后套用.运行后显示结果为即置信度为0.95时均值p的置信区间为(62.8328,79.5672).这与运用数理统计公式和查标准正态分布函数数值表计算的结果完全一致.运用normfit()函数计算该问题的结果为在相同置信度0.95时均值p的置信区间为(59.9912,82.4088),误差较大.根据数理统计知识，均值p已知时正态总体标准差o的区间估计采用自由度为n的卡方统计量，置信度为1-a的o置信区间为其中入1=x2(1-a/2;n),入2=x2(a/2;n).例3设总体Z~N(p,o2),o为待估参数•样本的一组观察值为(14.6,15.1,14.9,14.8,15.2,15.1),置信度为95%,求p=14.5时o的置信区间.解建立如下M-文件(并保存为myfun2.m)运行后结果显示为即置信度为95%时所求的置信区间为(0.3190,1.0900),这与运用数理统计公式和查卡方分布临界值表计算的结果完全一致.运用normfit()函数计算该问题的结果为在相同置信度95%下的置信区间为(0.1410,0.5539),误差极大.综上所述,在计算标准差o已知时均值p的置信区间和均值p已知时标准差o的置信区间不能再套用MATLAB所提供的现成函数normfit(),而必须重新编写程序•笔者在本文中写出的程序较好地解决了这两个问题的MATLAB实现并且有较强的实用性,有兴趣的读者只需调整显著性水平和更改样本数据便可求出实际问题在给定置