- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验项目8: 数据的统计分析
一、实验目的
掌握假设检验和因素分析方法以及如何用MATLAB实现拟合。
二、实验要求
掌握用Matlab进行假设检验和因素分析的方法。
三、实验内容
1、主要命令和注意事项:
(1)在Matlab 中Z 检验法由函数ztest 来实现,命令为
[h,p,ci]=ztest(x,mu,sigma,alpha,tail)
其中输入参数x 是样本,mu 是H0 中的μ0 ,sigma 是总体标准差σ ,alpha 是显著性水平α (alpha 缺省时设定为0.05),tail 是对备选假设H1 的选择:H1 为μ0 -μ ≠ 0时用tail=0(可缺省);H1 为μ-μ0 0 时用tail=1; H1 为μ-μ0 0时用tail=-1。输出参数h=0 表示接受H0 ,h=1 表示拒绝H0,p 表示在假设H0 下样本均值出现的概率,p越小H0 越值得怀疑,ci 是μ0 的置信区间。
例 某车间用一台包装机包装糖果。包得的袋装糖重是一个随机变量,它服从正
态分布。当机器正常时,其均值为0.5 公斤,标准差为0.015 公斤。某日开工后为检验
包装机是否正常,随机地抽取它所包装的糖9 袋,称得净重为(公斤):
0.497 0.506 0.518 0.524 0.498 0.511 0.520 0.515 0.512
问机器是否正常?
解 总体σ 已知,x ~ N(μ,0.0152 ),μ 未知。于是提出假设: H0: μ = μ0 =0.5 和H1 μ ≠ 0.5。
Matlab 实现如下:
x=[0.497 0.506 0.518 0.524 0.498 0.511 0.520 0.515 0.512];
[h,p,ci]=ztest(x,0.5,0.015)
求得h=1,p=0.0248,说明在0.05 的水平下,可拒绝原假设,即认为这天包装机工作不正常。
(2)在Matlab 中t 检验法由函数ttest 来实现,命令为
[h,p,ci]=ttest(x,mu,alpha,tail)
例 某种电子元件的寿命x (以小时计)服从正态分布,μ,σ 2 均未知.现得16 只元件的寿命如下:
159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170
问是否有理由认为元件的平均寿命大于225(小时)?
解 按题意需检验
H0 μ ≤225 H1: μ 225 ,
取α = 0.05。Matlab实现如下:
x=[159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170];
[h,p,ci]=ttest(x,225,0.05,1)
求得 h=0,p=0.2570,说明在显著水平为0.05 的情况下,不能拒绝原假设,认为元件的平均寿命不大于225 小时。
(3)两个正态总体均值差的检验( t 检验)
在Matlab 中由函数ttest2 实现,命令为:
[h,p,ci]=ttest2(x,y,alpha,tail)
与上面的ttest 相比,不同处只在于输入的是两个样本x,y(长度不一定相同),而不是一个样本和它的总体均值;tail 的用法与ttest 相似,可参看帮助系统。
例 在平炉上进行一项试验以确定改变操作方法的建议是否会增加钢的得率,试验是在同一平炉上进行的。每炼一炉钢时除操作方法外,其它条件都可能做到相同。先用标准方法炼一炉,然后用建议的新方法炼一炉,以后交换进行,各炼了10 炉,其得率分别为
1°标准方法 78.1 72.4 76.2 74.3 77.4 78.4 76.0 75.6 76.7 77.3
2°新方法 79.1 81.0 77.3 79.1 80.0 79.1 79.1 77.3 80.2 82.1
设这两个样本相互独立且分别来自正态总体, , 均未知,问建议的新方法能否提高得率?(取α = 0.05。)
解 (i)需要检验假设
(ii)Matlab 实现
x=[78.1 72.4 76.2 74.3 77.4 78.4 76.0 75.6 76.7 77.3];
y=[79.1 81.0 77.3 79.1 80.0 79.1 79.1 77.3 80.2 82.1];
[h,p,ci]=ttest2(x,y,0.05,-1)
求得h=1,p=2.2126×10-4。表明在α = 0.05的显著水平下,可以拒绝原假设,即认
为建议的新操作方法较原方法优。
(4)Matlab 统计工具箱中单因素方差分析的命令是anoval。
若各组数据个数相等,称为均衡数据。若各组数据个数不等,称非均衡数据。
1)均衡数据
处理均衡数据的用
文档评论(0)