matlab中exsolu函数,正态总体参数区间估计的MATLAB实现

2010年2月

吉林师范大学学报(自然科学版)l .1第1期Journal of Jilin Normal University (Natural Science Edition)Feb.2010

收稿日期:2010-01-05

作者简介:陈少云(1969-),男,重庆合川人,现为四川建筑职业技术学院副教授.研究方向:高等数学教学.

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

陈少云

(四川建筑职业技术学院计算机系,四川德阳618000)

摘 要:本文介绍了MATLAB 软件的normfit()函数在求解正态总体参数的区间估计中的长处和短处,结合实例编写了MATLAB 程序求解标准差R 已知时均值L 的置信区间和均值L 已知时标准差R 的置信区间,弥补了normfit()函数在该方面的不足.

关键词:正态总体;均值;标准差;置信区间

中图分类号:O212.2 文献标识码:A 文章编号:1674-3873-(2010)01-0076-03

0 引言

总体参数的点估计作为待估参数的近似值给出了明确的数量描述,在统计分析中有多方面的应用.但点估计没有给出这种近似的精确程度和可信程度,使其在实际应用中受到很大的限制,区间估计却可以弥补这一不足.

在工程技术中广泛使用的数学软件MATLAB 提供了现成的函数normfit()可以很方便地求出正态总体标准差R 未知时均值L 的置信区间和均值L 未知时标准差R 的置信区间,与数理统计公式和查相关的临界值表计算的结果完全吻合.

其函数调用格式为:[muhat,sigmahat,muci,sigmaci]=normfit(x,alpha).其中muhat,sigmahat 分别为正态分布的参数L 和R 的估计值,muci,sigmaci 分别为它们的置信区间,置信度为(1-alpha)@100%,alpha 为显著性水平.

例1 商店用机器包装某种商品,每包重量X 服从正态分布,为检查包装的质量,对机器包装的商品抽测8包,其重量为

5.08 4.97 5.12 5.05 4.95 4.90 5.00 5.01

试估计机器包装的商品重量的均值L 和标准差R 的置信度为0.95的置信区间

解 在MATLAB 的命令窗口输入如下语句

>>x =[5.08 4.97 5.12 5.05 4.95 4.90 5.00 5.01];

>>[muhat,sigmahat,muci,sigmaci]=nor mfit(x,0.05)

结果显示为

muhat=5.0100

sigmahat=0.0717

muci=4.9500 5.0700

sigmaci=0.0474 0.1460

结果表明均值L 的估计值为5.0100,其置信度为95%的置信区间为(4.9500,5.0700),这与标准差R 未知时运用数理统计公式和查T 分布临界值表计算的结果一致;标准差R 的估计值为0.0717,其置信度为95%的置信区间为(0.0474,0.1460),这与均值L 未知时运用数理统计公式和查卡方分布临界值表计算的结果一致.

但在计算标准差R 已知时均值L 的置信区间和均值L 已知时标准差R 的置信区间方面有较大的误#

76#