数据的统计分析与描述-matlab.ppt

Matlab相关命令介绍 频数直方图或频数表 频数直方图或频数表 频数直方图或频数表 频数直方图或频数表 参数估计 参数估计：点估计 点估计举例 参数估计：区间估计 区间估计举例 假设检验 假设检验 在Matlab中，假设检验问题都提出两种假设：即原假设和备择假设。对于正态总体均值的假设检验给出了检验函数： ztest 已知 ，检验正态总体均值 ； ttest 未知 ，检验正态总体均值 ； ttest2 两个正态总体均值比较。 1 单个正态总体N（ ）的假设检验 已知，对期望 的假设检验—Z检验法 调用函数 H=ztest（X，m,sigma） H=ztest（X，m,sigma,alpha） [H,sig,ci]=ztest（X，m,sigma,alpha,tail） 说明：X：样本；m:期望值；sigma:正态总体标准差；alpha:经验水平 ； tail:备择假设的选项，若tail=0(缺省)，则 ； 若tail=1，则 ；若tail=-1，则 。 即tail=0(缺省)为双边检验，其余为单边检验问题。 H：检验结果，分两种情况：若H=0，则在水平 下，接受原假设；若H=1，则在水平 下，拒绝原假设。 sig为当原假设为真时（即 成立），得到观察值的概率，当sig为小概率时，则对原假设提出质疑。Ci:均值 的1-alpha置信区间。 应用举例 例7．1 某车间用一台包装机包装葡萄糖，包得的袋装糖重是一个随机变量，它服从正态分布。当机器正常时，其均值为0.5公斤，标准差为0.015。某日开工后检验包装机是否正常，随机地抽取所包装的糖9袋，称得净重为：（公斤） 0.497 0.518 0.524 0.498 0.511 0.52 0.515 0.512 问机器是否正常？ 解： 已知，在水平 =0.05下检验假设： 原假设： 备择假设： 程序： X=[0.497,0.506,0.518,0.524,0.498,0.511,0.52,0.515,0.512]; [H,SIG]=ztest(X,0.5,0.015,0.05,0) 运行后显示结果如下： H=1 SIG=0.0248 结果表明：H=1，说明在水平=0.05下，可拒绝原假设，即认为包装机工作不正常。 1 单个正态总体N（ ）的假设检验 未知，对期望 的假设检验—t检验法 调用函数 H=ttest(X,m,sigma) %在水平 =sigma下检验是否成立。 说明：X：样本；m:期望值；alpha:经验水平 ； tail:备择假设的选项，若tail=0(缺省)，则备择假设为 ；若tail=1，则 ；若tail=-1，则 。即tail=0(缺省)为双边检验，其余为单边检验问题。 H：检验结果，分两种情况：若H=0，则在水平 下，接受原假设；若H=1，则在水平 下，拒绝原假设。 sig为当原假设为真时（即 成立），得到观察值的概率，当sig为小概率时，则对原假设提出质疑。Ci:均值 的1-alpha置信区间。 应用举例 例7.2 某种电子元件的寿命X（以小时计）服从正态分布， 均未知，现测得16只元件寿命如下： 159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170 问是否有理由认为元件的平均寿命大于225（小时）？ 解： 未知，在水平 =0.05下检验假设： 程序： X=[159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170]; [H,SIG]=ttest(X,225,0.05,1) 运行后显示结果如下： H=0 SIG=0.2570 结果表明：H=0，说明在水平=0.05下，应接受原假设，即认为元件的平均寿命不大于225小时。 2、两个正态总体均值差的检验（t检验） 调用函数 [h,sig,ci]=ttest(X,Y) [h,sig,ci]=ttest2(X,Y,alpha) [h,sig,ci]=ttest2(X,Y,alpha,tail) 说明：原假设为：当tail=0时，表示 （缺省）；当tail=1时，表示