Matlab概率统计Word格式文档下载.docx

Matlab概率统计Word格式文档下载.docx

《Matlab概率统计Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《Matlab概率统计Word格式文档下载.docx（9页珍藏版）》请在冰豆网上搜索。

解:

总体μ和σ已知,该问题是当为已知时,在水平下,根据样本值判断μ=0.5还是.为此提出假设:

备择假设:

>

X=[0.497,0.506,0.518,0.524,0.498,0.511,0.52,0.515,0.512];

[h,sig,ci,zval]=ztest（X,0.5,0.015,0.05,0）

结果显示为

h=

1

sig=

0.0248%样本观察值的概率

ci=

0.50140.5210%置信区间,均值0.5在此区间之外

zval=

2.2444%统计量的值

结果表明:

h=1,说明在水平下,可拒绝原假设,即认为包装机工作不正常.

4.8.2未知,单个正态总体的均值μ的假设检验（t检验法）

函数ttest

格式h=ttest（x,m）%x为正态总体的样本,m为均值μ0,显著性水平为0.05

h=ttest（x,m,alpha）%alpha为给定显著性水平

[h,sig,ci]=ttest（x,m,alpha,tail）%sig为观察值的概率,当sig为小概率时则对原假设提出质疑,ci为真正均值μ的1-alpha置信区间.

若tail=0,表示备择假设:

例4-75某种电子元件的寿命X（以小时计）服从正态分布,,σ2均未知.现测得16只元件的寿命如下

159280101212224379179264222362168250

149260485170

问是否有理由认为元件的平均寿命大于225（小时）

未知,在水平下检验假设:

:

:

X=[159280101212224379179264222362168250149260485170];

[h,sig,ci]=ttest（X,225,0.05,1）

结果显示为:

0.2570

198.2321Inf%均值225在该置信区间内

H=0表示在水平下应该接受原假设,即认为元件的平均寿命不大于225小时.

4.8.3两个正态总体均值差的检验（t检验）

两个正态总体方差未知但等方差时,比较两正态总体样本均值的假设检验

函数ttest2

格式[h,sig,ci]=ttest2（X,Y）%X,Y为两个正态总体的样本,显著性水平为0.05

[h,sig,ci]=ttest2（X,Y,alpha）%alpha为显著性水平

[h,sig,ci]=ttest2（X,Y,alpha,tail）%sig为当原假设为真时得到观察值的概率,当sig为小概率时则对原假设提出质疑,ci为真正均值μ的1-alpha置信区间.

（为X为期望值,为Y的期望值）

例4-76在平炉上进行一项试验以确定改变操作方法的建议是否会增加钢的产率,试验是在同一只平炉上进行的.每炼一炉钢时除操作方法外,其他条件都尽可能做到相同.先用标准方法炼一炉,然后用建议的新方法炼一炉,以后交替进行,各炼10炉,其产率分别为

（1）标准方法:

78.172.476.274.377.478.476.075.576.777.3

（2）新方法:

79.181.077.379.180.079.179.177.380.282.1

设这两个样本相互独立,且分别来自正态总体和,,,均未知.问建议的新操作方法能否提高产率（取α=0.05）

两个总体方差不变时,在水平下检验假设:

X=[78.172.476.274.377.478.476.075.576.777.3];

Y=[79.181.077.379.180.079.179.177.380.282.1];

[h,sig,ci]=ttest2（X,Y,0.05,-1）

2.1759e-004%说明两个总体均值相等的概率很小

-Inf-1.9083

H=1表示在水平下,应该拒绝原假设,即认为建议的新操作方法提高了产率,因此,比原方法好.

4.8.4两个总体一致性的检验——秩和检验

函数ranksum

格式p=ranksum（x,y,alpha）%x,y为两个总体的样本,可以不等长,alpha为显著性水平

[p,h]=ranksum（x,y,alpha）%h为检验结果,h=0表示X与Y的总体差别不显著h=1表示X与Y的总体差别显著

[p,h,stats]=ranksum（x,y,alpha）%stats中包括:

ranksum为秩和统计量的值以及zval为过去计算p的正态统计量的值

说明P为两个总体样本X和Y为一致的显著性概率,若P接近于0,则不一致较明显.

例4-77某商店为了确定向公司A或公司B购买某种商品,将A和B公司以往的各次进货的次品率进行比较,数据如下所示,设两样本独立.问两公司的商品的质量有无显著差异.设两公司的商品的次品的密度最多只差一个平移,取α=0.05.

A:

7.03.59.68.16.25.110.44.02.010.5

B:

5.73.24.111.09.76.93.64.85.68.410.15.512.3

设,分别为A,B两个公司的商品次品率总体的均值.则该问题为在水平α=0.05下检验假设:

A=[7.03.59.68.16.25.110.44.02.010.5];

B=[5.73.24.111.09.76.93.64.85.68.410.15.512.3];

[p,h,stats]=ranksum（A,B,0.05）

结果为:

p=

0.8041

stats=

zval:

-0.2481

ranksum:

116

一方面,两样本总体均值相等的概率为0.8041,不接近于0;

另一方面,H=0也说明可以接受原假设,即认为两个公司的商品的质量无明显差异.

4.8.5两个总体中位数相等的假设检验——符号秩检验

函数signrank

格式p=signrank（X,Y,alpha）%X,Y为两个总体的样本,长度必须相同,alpha为显著性水平,P两个样本X和Y的中位数相等的概率,p接近于0则可对原假设质疑.

[p,h]=signrank（X,Y,alpha）%h为检验结果:

h=0表示X与Y的中位数之差不显著,h=1表示X与Y的中位数之差显著.

[p,h,stats]=signrank（x,y,alpha）%stats中包括:

signrank为符号秩统计量的值以及zval为过去计算p的正态统计量的值.

例4-78两个正态随机样本的中位数相等的假设检验

x=normrnd（0,1,20,1）;

y=normrnd（0,2,20,1）;

[p,h,stats]=signrank（x,y,0.05）

0.3703

-0.8960

signedrank:

81

h=0表示X与Y的中位数之差不显著

4.8.6两个总体中位数相等的假设检验——符号检验

函数signtest

格式p=signtest（X,Y,alpha）%X,Y为两个总体的样本,长度必须相同,alpha为显著性水平,P两个样本X和Y的中位数相等的概率,p接近于0则可对原假设质疑.

[p,h]=signtest（X,Y,alpha）%h为检验结果:

[p,h,stats]=signtest（X,Y,alpha）%stats中sign为符号统计量的值

例4-79两个正态随机样本的中位数相等的假设检验

X=normrnd（0,1,20,1）;

Y=normrnd（0,2,20,1）;

[p,h,stats]=signtest（X,Y,0.05）

0.2632

sign:

7

4.8.7正态分布的拟合优度测试

函数jbtest

格式H=jbtest（X）%对输入向量X进行Jarque-Bera测试,显著性水平为0.05.

H=jbtest（X,alpha）%在水平alpha而非5%下施行Jarque-Bera测试,alpha在0和1之间.

[H,P,JBSTAT,CV]=jbtest（X,alpha）%P为接受假设的概率值,P越接近于0,则可以拒绝是正态分布的原假设;

JBSTAT为测试统计量的值,CV为是否拒绝原假设的临界值.

说明H为测试结果,若H=0,则可以认为X是服从正态分布的;

若X=1,则可以否定X服从正态分布.X为大样本,对于小样本用lillietest函数.

例4-80调用MATLAB中关于汽车重量的数据,测试该数据是否服从正态分布

loadcarsmall

[h,p,j,cv]=jbtest（Weight）

0.0267

j=

7.2448

cv=

5.9915

说明p=2.67%表示应该拒绝服从正态分布的假设;

h=1也可否定服从正态分布;

统计量的值j=7.2448大于接受假设的临界值cv=5.9915,因而拒绝假设（测试水平为5%）.

4.8.8正态分布的拟合优度测试

函数lillietest

格式H=lillietest（X）%对输入向量X进行Lilliefors测试,显著性水平为0.05.

H=lillietest（X,alpha）%在水平alpha而非5%下施行Lilliefors测试,alpha在0.01和0.2之间.

[H,P,LSTAT,CV]=lillietest（X,alpha）%P为接受假设的概率值,P越接近于0,则可以拒绝是正态分布的原假设;

LSTAT为测试统计量的值,CV为是否拒绝原假设的临界值.

若X=1,则可以否定X服从正态分布.

例4-81

Y=chi2rnd（10,100,1）;

[h,p,l,cv]=lillietest（Y）

0.0175

l=

0.1062

0.0886

说明h=1表示拒绝正态分布的假设;

p=0.0175表示服从正态分布的概率很小;

统计量的值l=0.1062大于接受假设的临界值cv=0.0886,因而拒绝假设（测试水平为5%）.

hist（Y）

从图中看出,数据Y不服从正态分布.

4.8.9单个样本分布的Kolmogorov-Smirnov测试

函数kstest

格式H=kstest（X）%测试向量X是否服从标准正态分布,测试水平为5%.

H=kstest（X,cdf）%指定累积分布函数为cdf的测试（cdf=[]时表示标准正态分布）,测试水平为5%

H=kstest（X,cdf,alpha）%alpha为指定测试水平

[H,P,KSSTAT,CV]=kstest（X,cdf,alpha）%P为原假设成立的概率,KSSTAT为测试统计量的值,CV为是否接受假设的临界值.

说明原假设为X服从标准正态分布.若H=0则不能拒绝原假设,H=1则可以拒绝原假设.

例4-82产生100个威布尔随机数,测试该随机数服从的分布

x=weibrnd（1,2,100,1）;

[H,p,ksstat,cv]=kstest（x,[xweibcdf（x,1,2）],0.05）%测试是否服从威布尔分布

H=

0.3022

ksstat=

0.0959

0.1340

说明H=0表示接受原假设,统计量ksstat小于临界值表示接受原假设.

[H,p,ksstat,cv]=kstest（x,[xexpcdf（x,1）],0.05）%测试是否服从指数分布

0.0073

0.1653

说明H=1表明拒绝服从指数分布的假设.

[H,p,ksstat,cv]=kstest（x,[],0.05）%测试是否服从标准正态分布

3.1285e-026

0.5380

说明H=1表明不服从标准正态分布.

4.8.10两个样本具有相同的连续分布的假设检验

函数kstest2

格式H=kstest2（X1,X2）%测试向量X1与X2是具有相同的连续分布,测试水平为5%.

H=kstest2（X1,X2,alpha）%alpha为测试水平

[H,P,KSSTAT]=kstest（X,cdf,alpha）%与指定累积分布cdf相同的连续分布,P为假设成立的概率,KSSTAT为测试统计量的值.

说明原假设为具有相同连续分布.测试结果为H,若H=0,表示应接受原假设;

若H=1,表示可以拒绝原假设.这是Kolmogorov-Smirnov测试方法.

例4-83

x=-1:

1:

5;

y=randn（20,1）;

[h,p,k]=kstest2（x,y）

0.0444

k=

0.5643

说明h=1表示可以认为向量x与y的分布不相同,相同的概率只有4.4%.

4.9方差分析

4.9.1单因素方差分析

单因素方差分析是比较两组或多组数据的均值,它返回原假设——均值相等的概率

函数anova1

格式p=anova1（X）%X的各列为彼此独立的样本观察值,其元素个数相同,p为各列均值相等的概率值,若p值接近于0,则原假设受到怀疑,说明至少有一列均值与其余列均值有明显不同.

p=anova1（X,group）%X和group为向量且group要与X对应

p=anova1（X,group,'

displayopt'

）%displayopt=on/off表示显示与隐藏方差分析表图和盒图

[p,table]=anova1（…）%table为方差分析表

[p,table,stats]=anova1（…）%stats为分析结果的构造

说明anova1函数产生两个图:

标准的方差分析表图和盒图.

方差分析表中有6列:

第1列（source）显示:

X中数据可变性的来源;

第2列（SS）显示:

用于每一列的平方和;

第3列（df）显示:

与每一种可变性来源有关的自由度;

第4列（MS）显示:

是SS/df的比值;

第5列（F）显示:

F统计量数值,它是MS的比率;

第6列显示:

从F累积分布中得到的概率,当F增加时,p值减少.

例4-84设有3台机器,用来生产规格相同的铝合金薄板.取样测量薄板的厚度,精确至‰厘米.得结果如下:

机器1:

0.2360.2380.2480.2450.243

机器2:

0.2570.2530.2550.2540.261

机器3:

0.2580.2640.2590.2670.262

检验各台机器所生产的薄板的厚度有无显著的差异

X=[0.2360.2380.2480.2450.243;

0.2570.2530.2550.2540.261;

…

0.2580.2640.2590.2670.262];

P=anova1（X'

）

P=

1.3431e-005

还有两个图,即图4-22和图4-23.

图4-22图4-23

例4-85建筑横梁强度的研究:

3000磅力量作用在一英寸的横梁上来测量横梁的挠度,钢筋横梁的测试强度是:

8286798384858687;

其余两种更贵的合金横梁强度测试为合金1:

748278757677;

合金2:

797977788279].

检验这些合金强度有无明显差异

strength=[8286798384858687748278757677797977788279];

alloy={'

st'

'

'

al1'

…

'

al2'

};

[p,table,stats]=anova1（strength,alloy,'

on'

结果为

1.5264e-004

table=

Source'

'

SS'

df'

MS'

F'

Prob>

Groups'

[184.8000][2][92.4000][15.4000][1.5264e-004]

Error'

[102.0000][17][6.0000][][]

Total'

[286.8000][19][][][]

gnames:

{3x1cell}

n:

[866]

source:

anova1'

means:

[847779]

df:

17

s:

2.4495

图4-24图4-25

说明p值显示,3种合金是明显不同的,盒图显示钢横梁的挠度大于另两种合金横梁的挠度.

4.9.2双因素方差分析

函数anova2

格式p=anova2（X,reps）

p=anova2（X,reps,'

[p,table]=anova2（…）

[p,table,stats]=anova2（…）

说明执行平衡的双因素试验的方差分析来比较X中两个或多个列（行）的均值,不同列的数据表示因素A的差异,不同行的数据表示另一因素B的差异.如果行列对有多于一个的观察点,则变量reps指出每一单元观察点的数目,每一单元包含reps行,如:

reps=2

其余参数与单因素方差分析参数相似.

例4-86一火箭使用了4种燃料,3种推进器作射程试验,每种燃料与每种推进器的组合各发射火箭2次,得到结果如下:

推进器（B）B1B2B3

A158.200056.200065.3000

52.600041.200060.8000

A249.100054.100051.6000

燃料A42.800050.500048.4000

A360.100070.900039.2000

58.300073.200040.7000

A475.800058.200048.7000

71.500051.000041.4000

考察推进器和燃料这两个因素对射程是否有显著的影响

建立M文件

X=[58.200056.200065.3000

49.100054.100051.6000

42.800050.500048.4000

60.100070.900039.2000

75.800058.200048.7000

71.500051.000041.4000];

P=anova2（X,2）

0.00350.02600.0001

显示方差分析图为图4-26.

图4-26

Matlab6.0数学手册