8matlab在概率统计中的应用课件

在概率统计中的应用许海洋...1在概率统计中的应用许海洋...18.1分布率和概率密度函数8.1.1通用函数计算概率密度函数值8.1.2专用函数计算概率密度函数值8.1.3常见分布的密度函数作图8.2随机变量的累积概率值(分布函数值)8.2.1通用函数计算累积概率值8.2.2专用函数计算累积概率值8.3样本描述8.3.1集中趋势8.3.2离中趋势...28.1分布率和概率密度函数...28.4随机变量的数字特征8.4.1期望8.4.2方差8.4.3常见分布的期望和方差8.4.4协方差与相关系数8.5参数估计8.5.1点估计8.5.2区间估计8.5.3常见分布的参数估计...38.4随机变量的数字特征...38.6假设检验8.6.12已知，单个正态总体的均值μ的假设检验(U检验法）8.6.22未知，单个正态总体的均值μ的假设检验(t检验法)8.6.3两个正态总体均值差的检验（t检验）8.6.4两个总体一致性的检验——秩和检验8.6.5两个总体中位数相等的假设检验——符号秩检验8.6.6两个总体中位数相等的假设检验——符号检验8.7方差分析8.7.1单因素方差分析8.7.2双因素方差分析

...48.6假设检验...4利用MATLAB统计工具箱，可以进行基本概论和数理统计分析，以及进行比较复杂的多元统计分析。概率论...5利用MATLAB统计工具箱，可以进行基本概论和8.1分布率和概率密度函数8.1.1通用函数计算概率密度函数值格式Y=pdf(name，K，A)Y=pdf(name，K，A，B)Y=pdf(name，K，A，B，C)返回在X=K处、参数为A、B、C的概率密度值，对于不同的分布，参数个数是不同；name为分布函数名，其取值如表。...68.1分布率和概率密度函数8.1.1通用函数计算概率密度name的取值函数说明'beta'或'Beta'Beta分布'bino'或'Binomial'二项分布'chi2'或'Chisquare'卡方分布'exp'或'Exponential'指数分布'f'或'F'F分布'gam'或'Gamma'GAMMA分布'geo'或'Geometric'几何分布'norm'或'Normal'正态分布'poiss'或'Poisson'泊松分布't'或'T'T分布'unif'或'Uniform'均匀分布常见分布函数表...7name的取值函数说明'beta'或'Beta'Beta分布例如二项分布：设一次试验，事件A发生的概率为p，那么，在n次独立重复试验中，事件A恰好发生K次的概率P_K为：P_K=P{X=K}=pdf('bino'，K，n，p)例计算正态分布N（0，1）的随机变量X在点0.6578的密度函数值。

pdf('norm',0.6578,0,1)ans=0.3213例自由度为8的卡方分布，在点2.18处的密度函数值。

pdf('chi2',2.18,8)ans=0.0363...8例如二项分布：设一次试验，事件A发生的概率为p，8.1.2专用函数计算概率密度函数值（1）二项分布的概率值binopdf(k,n,p)p—每次试验事件A发生的概率；K—事件A发生K次；n—试验总次数（2）泊松分布的概率值poisspdf(k,Lambda)

（3）正态分布的概率值normpdf(K,mu,sigma)%计算参数为μ=mu，σ=sigma的正态分布密度函数在K处的值...98.1.2专用函数计算概率密度函数值（1）二项分布的概函数名调用形式注

释Unifpdfunifpdf(x,a,b)[a,b]上均匀分布(连续)概率密度在X=x处的函数值unidpdfUnidpdf(x,n)均匀分布（离散）概率密度函数值Exppdfexppdf(x,Lambda)参数为Lambda的指数分布概率密度函数值normpdfnormpdf(x,mu,sigma)参数为mu，sigma的正态分布概率密度函数值chi2pdfchi2pdf(x,n)自由度为n的卡方分布概率密度函数值Tpdftpdf(x,n)自由度为n的t分布概率密度函数值Fpdffpdf(x,n1,n2)第一自由度为n1,第二自由度为n2的F分布概率密度函数值binopdfbinopdf(x,n,p)参数为n,p的二项分布的概率密度函数值geopdfgeopdf(x,p)参数为p的几何分布的概率密度函数值专用函数计算概率密度函数表...10函数名调用形式注释Unifpdfunifp例绘制卡方分布密度函数在自由度分别为1、5、15的图形x=0:0.1:30;y1=chi2pdf(x,1);plot(x,y1,':')holdony2=chi2pdf(x,5);plot(x,y2,'+')y3=chi2pdf(x,15);plot(x,y3,'o')axis([0,30,0,0.2])%指定显示的图形区域...11例绘制卡方分布密度函数在自由度分别为例计算参数为mu和1的正态分布概率密度函数在1.5处的值，其中mu为1到2之间以0.2为间隔的小数。mu=[0:0.2:2];y=normpdf(1.5,mu,1)y=0.12950.17140.21790.26610.31230.35210.38140.39700.39700.38140.3521...12例计算参数为mu和1的正态分布概率密度函数在1.8.1.3常见分布的密度函数作图1．二项分布x=0:10;y=binopdf(x,10,0.5);plot(x,y,'+')2．卡方分布x=0:0.2:15;y=chi2pdf(x,4);plot(x,y)

...138.1.3常见分布的密度函数作图1．二项分布2．卡方分3．非中心卡方分布x=(0:0.1:10)';p1=ncx2pdf(x,4,2);p=chi2pdf(x,4);plot(x,p,'--',x,p1,'-')4．指数分布x=0:0.1:10;y=exppdf(x,2);plot(x,y)

...143．非中心卡方分布4．指数分布...148.2随机变量的累积概率值(分布函数值)为X的分布函数。如果知道X的分布函数，就可以知道落在任一区间(x1,x2)上的概率。若X为随机变量，x为任意实数，则函数...158.2随机变量的累积概率值(分布函数值)为X8.2.1通用函数计算累积概率值通用函数cdf用来计算随机变量X≤K的概率之和（累积概率值）

cdf('name',K,A)cdf('name',K,A,B)cdf('name',K,A,B,C)返回以name为分布、随机变量X≤K的概率之和的累积概率值，name的取值见’常见分布函数表’...168.2.1通用函数计算累积概率值通用函数c例求标准正态分布随机变量X落在区间(-∞，0.4)内的概率。cdf('norm',0.4,0,1)ans=0.6554例求自由度为16的卡方分布随机变量落在[0，6.91]内的概率cdf('chi2',6.91,16)ans=0.0250...17例求标准正态分布随机变量X落在区间(-∞，0.48.2.2专用函数计算累积概率值（随机变量X≤K的概率之和）...188.2.2专用函数计算累积概率值...18则有：p1=normcdf(5,3,2)-normcdf(2,3,2)p1=0.5328p2=normcdf(10,3,2)-normcdf(-4,3,2)p2=0.9995p3=1-normcdf(2,3,2)-normcdf(-2,3,2)p3=0.6853p4=1-normcdf(3,3,2)p4=0.5000

...19则有：...19专用函数的累积概率值函数表函数名调用形式注

释unifcdfunifcdf(x,a,b)[a,b]上均匀分布(连续)累积分布函数值F(x)=P{X≤x}unidcdfunidcdf(x,n)均匀分布（离散）累积分布函数值F(x)=P{X≤x}expcdfexpcdf(x,Lambda)参数为Lambda的指数分布累积分布函数值F(x)=P{X≤x}normcdfnormcdf(x,mu,sigma)参数为mu，sigma的正态分布累积分布函数值F(x)=P{X≤x}chi2cdfchi2cdf(x,n)自由度为n的卡方分布累积分布函数值F(x)=P{X≤x}tcdftcdf(x,n)自由度为n的t分布累积分布函数值F(x)=P{X≤x}fcdffcdf(x,n1,n2)第一自由度为n1,第二自由度为n2的F分布累积分布函数值gamcdfgamcdf(x,a,b)分布累积分布函数值F(x)=P{X≤x}betacdfbetacdf(x,a,b)分布累积分布函数值F(x)=P{X≤x}...20专用函数的累积概率值函数表函数名调用形式注描述样本数据集中趋势的统计量有算术平均值、中位数、众数、几何均值、调和均值和截尾均值等。8.3样本描述8.3.1集中趋势...21描述样本数据集中趋势的统计量有算术平均值、中位数（1）利用mean求算术平均值

mean(X)%X为向量，返回X中各元素的平均值

mean(A)%A为矩阵，返回A中各列元素的平均值构成的向量

mean(A,dim)%在给出的维数内的平均值...22（1）利用mean求算术平均值...22例A=[1345;2346;1315]A=134523461315mean(A)ans=1.33333.00003.00005.3333mean(A,1)ans=1.33333.00003.00005.3333...23例...23（2）利用geomean计算几何平均数

M=geomean(X)

X为向量，返回X中各元素的几何平均数。M=geomean(A)

A为矩阵，返回A中各列元素的几何平均数构成的向量。...24（2）利用geomean计算几何平均数...24例B=[1345]B=1345M=geomean(B)M=2.7832A=[1345;2346;1315]A=134523461315M=geomean(A)M=1.25993.00002.51985.3133...25例...25（3）利用median计算中值（中位数）

median(X)%X为向量，返回X中各元素的中位数。

median(A)%A为矩阵，返回A中各列元素的中位数构成的向量。

median(A,dim)%求给出的维数内的中位数...26（3）利用median计算中值（中位数）...26例A=[1345;2346;1315]

A=134523461315median(A)

ans=1345...27例...27（4）利用harmmean求调和平均值

M=harmmean(X)%X为向量，返回X中各元素的调和平均值。

M=harmmean(A)%A为矩阵，返回A中各列元素的调和平均值构成的向量。调和平均值的数学含义是其中：样本数据非0，主要用于严重偏斜分布。...28（4）利用harmmean求调和平均值...28例B=[1345]B=1345M=harmmean(B)M=2.2430A=[1345;2346;1315]A=134523461315M=harmmean(A)M=1.20003.00002.00005.2941...29例...298.3.2离中趋势描述离散趋势的统计量包括均值绝对差、极差、方差和标准差。（1）均值绝对差y=mad(X)若X为矢量，则y用mead(abs(X-mean(X)))计算；若X为矩阵，则y为包含X中每列数据均值绝对差的行矢量。mad(X,0):与mad(X)相同，使用均值。mad(X,1):基于中值计算y,即:median(abs(X-median(X))).（2）极差y=range(X)返回X中数据的最小值与最大值之间的差值。...308.3.2离中趋势描述离散趋势的统计量包括8.4随机变量的数字特征8.4.1期望计算样本均值函数meanX-2-1012P0.30.10.20.10.3例设随机变量X的分布律为：求E(X)E(X2-1)解：在Matlab编辑器中建立M文件如下：X=[-2-1012];p=[0.30.10.20.10.3];EX=sum(X.*p)Y=X.^2-1EY=sum(Y.*p)运行后结果如下：EX=0Y=30-103EY=1.6000

...318.4随机变量的数字特征8.4.1期望计算样本均值函8.4.2方差求样本方差...328.4.2方差求样本方差...32求标准差...33求标准差...33例求下列样本的样本方差和样本标准差，方差和标准差14.7015.2114.9015.3215.32解：X=[14.715.2114.914.9115.3215.32];DX=var(X,1)%方差

DX=0.0559sigma=std(X,1)%标准差sigma=0.2364DX1=var(X)%样本方差DX1=0.0671sigma1=std(X)%样本标准差sigma1=0.2590...34例求下列样本的样本方差和样本标准差，方差8.4.3常见分布的期望和方差（1）均匀分布（连续）的期望和方差

[M,V]=unifstat(A,B)

%A、B为标量时，就是区间上均匀分布的期望和方差，A、B也可为向量或矩阵，则M、V也是向量或矩阵。例a=1:6;b=2.*a;[M,V]=unifstat(a,b)M=1.50003.00004.50006.00007.50009.0000V=0.08330.33330.75001.33332.08333.0000...358.4.3常见分布的期望和方差（1）均匀分布（连续）的期（2）正态分布的期望和方差[M,V]=normstat(MU,SIGMA)

%MU、SIGMA可为标量也可为向量或矩阵，则M=MU，V=SIGMA2。例n=1:4;[M,V]=normstat(n'*n,n'*n)M=1234246836912481216V=149164163664936811441664144256...36（2）正态分布的期望和方差...36（3）二项分布的均值和方差[M,V]=binostat(N,P)

%N，P为二项分布的两个参数，可为标量也可为向量或矩阵。例n=logspace(1,5,5)n=10100100010000100000[M,V]=binostat(n,1./n)M=11111V=0.90000.99000.99900.99991.0000[m,v]=binostat(n,1/2)m=550500500050000v=1.0e+04*0.00030.00250.02500.25002.5000...37（3）二项分布的均值和方差...37常见分布的均值和方差函数名调用形式注

释unifstat[M,V]=unifstat(a,b)均匀分布(连续)的期望和方差，M为期望，V为方差unidstat[M,V]=unidstat(n)均匀分布（离散）的期望和方差expstat[M,V]=expstat(p,Lambda)指数分布的期望和方差normstat[M,V]=normstat(mu,sigma)正态分布的期望和方差chi2stat[M,V]=chi2stat(x,n)卡方分布的期望和方差tstat[M,V]=tstat(n)t分布的期望和方差fstat[M,V]=fstat(n1,n2)F分布的期望和方差gamstat[M,V]=gamstat(a,b)分布的期望和方差betastat[M,V]=betastat(a,b)分布的期望和方差...38常见分布的均值和方差函数名调用形式注释unifstat8.4.4协方差与相关系数

（1）协方差

cov(X)%求向量X的协方差

cov(A)%求矩阵A的协方差矩阵，该协方差矩阵的对角线元素是A的各列的方差，即：var(A)=diag(cov(A))。

cov(X,Y)%X,Y为等长列向量，等同于cov([XY])。...398.4.4协方差与相关系数（1）协方差...39例X=[0-11]';Y=[122]‘;C1=cov(X)%X的协方差C1=1C2=cov(X,Y)%列向量X、Y的协方差矩阵，对角线元素为各列向量的方差C2=1.0000000.3333A=[123;40-1;173]A=12340-1173...40例...40C1=cov(A)%求矩阵A的协方差矩阵C1=3.0000-4.5000-4.0000-4.500013.00006.0000-4.00006.00005.3333C2=var(A(:,1))%求A的第1列向量的方差C2=3C3=var(A(:,2))%求A的第2列向量的方差C3=13C4=var(A(:,3))C4=5.3333...41C1=cov(A)%求矩阵A的协方差矩阵...41（2）相关系数

corrcoef(X,Y)%返回列向量X,Y的相关系数，等同于corrcoef([XY])。

corrcoef(A)%返回矩阵A的列向量的相关系数矩阵...42（2）相关系数...42例A=[123;40-1;139]A=12340-1139C1=corrcoef(A)%求矩阵A的相关系数矩阵C1=1.0000-0.9449-0.8030-0.94491.00000.9538-0.80300.95381.0000C1=corrcoef(A(:,2),A(:,3))%求A的第2列与第3列列向量的相关系数矩阵C1=1.00000.95380.95381.0000...43例...438.5参数估计8.5.1点估计：用单个数值作为参数的估计（1）矩法：用总体的样本矩来估计总体的同阶矩。例随机取8个活塞环，测得它们的直径为(以mm计):74.00174.00574.00374.00174.00073.99874.00674.002,设环直径的测量值服从正态分布，现估计总体的方差。解：因为样本的2阶中心矩是总体方差的矩估计量，所以可以用moment函数进行估计。X=[74.00174.00574.00374.00174.00073.99874.00674.002];moment(X,2)ans=6.0000e-006...448.5参数估计8.5.1点估计：用单个数值作（2）最大似然法p=mle(‘dist’,data)使用data矢量中的样本数据，返回dist指定的分布的最大似然估计。例用最大似然估计法解上例。X=[74.00174.00574.00374.00174.00073.99874.00674.002];p=mle('norm',X);p(2)*p(2)ans=6.0000e-006...45（2）最大似然法p=mle(‘dist’,data)使8.5.2区间估计区间估计不仅仅给出了参数的近似取值，还给出了取该值的误差范围。...468.5.2区间估计区间估计不仅仅给出了参

phat=mle(‘dist’,X)

%返回用dist指定分布的最大似然估计值[phat,pci]=mle(‘dist’,X)

%置信度为95%的区间[phat,pci]=mle(‘dist’,X,alpha)

%置信度由alpha确定[phat,pci]=mle(‘dist’,X,alpha,pl)

%仅用于二项分布，pl为试验次数。dist为分布函数名，如：beta(分布)、bino（二项分布）等，X为数据样本，alpha为显著水平α，(1-α)×100%为置信度。...47phat=mle(‘dist’,X)例X=binornd(20,0.75)%产生二项分布的随机数X=11[p,pci]=mle('bino',X,0.05,20)%求概率的估计值和置信区间，置信度为95%p=0.5500pci=0.31530.7694...48例...48例从一批灯泡中随机地取5只作寿命试验,测得寿命(以小时计)为:10501100112012501280,设灯泡寿命服从正态分布,求灯泡寿命平均值的95%置信区间.X=[10501100112012501280];[p,ci]=mle('norm',X,0.05)p=1.0e+003*1.16000.0892ci=1.0e+003*1.0751-0.04991.24490.2283...49例从一批灯泡中随机地取5只作寿命试验,测得寿8.5.3常见分布的参数估计（1）β分布的参数a和b的最大似然估计值和置信区间PHAT=betafit(X)[PHAT,PCI]=betafit(X,ALPHA)PHAT为样本X的β分布的参数a和b的估计量PCI为样本X的β分布参数a和b的置信区间，是一个2×2矩阵，其第1例为参数a的置信下界和上界，第2例为b的置信下界和上界，ALPHA为显著水平，（1-α）×100%为置信度。...508.5.3常见分布的参数估计（1）β分布的参数a和b的最例随机产生100个β分布数据，相应的分布参数真值为4和3。则4和3的最大似然估计值和置信度为99%的置信区间为：解：X=betarnd(4,3,100,1);%产生100个β分布的随机数[PHAT,PCI]=betafit(X,0.01)%求置信度为99%的置信区间和参数a、b的估计值...51例随机产生100个β分布数据，相应的分结果显示PHAT=4.03173.1046PCI=2.72492.16445.33864.0448说明：估计值3.9010的置信区间是[2.52445.2776]，估计值2.6193的置信区间是[1.74883.4898]。...52结果显示...52（2）正态分布的参数估计

[muhat,sigmahat,muci,sigmaci]=normfit(X)[muhat,sigmahat,muci,sigmaci]=normfit(X,alpha)muhat,sigmahat分别为正态分布的参数μ和σ的估计值，muci,sigmaci分别为置信区间，其置信度为；alpha给出显著水平α，缺省时默认为0.05，即置信度为95%...53（2）正态分布的参数估计...53例有两组(每组100个元素)正态随机数据，其均值为10，均方差为2，求95%的置信区间和参数估计值。解：r=normrnd(10,2,100,2);%产生两列正态随机数据[mu,sigma,muci,sigmaci]=normfit(r)则结果为mu=10.09599.7460%各列的均值的估计值sigma=1.73701.8894%各列的均方差的估计值muci=9.75129.371110.440510.1209sigmaci=1.52511.65892.01782.1949说明muci，sigmaci中各列分别为原随机数据各列估计值的置信区间，置信度为95%...54例有两组(每组100个元素)正态随机例分别使用金球和铂球测定引力常数（1）用金球测定观察值为：6.6836.6816.6766.6786.6796.672（2）用铂球测定观察值为：6.6616.6616.6676.6676.664设测定值总体为，μ和σ为未知。对(1)、(2)两种情况分别求μ和σ的置信度为0.9的置信区间。解：建立M文件：X=[6.6836.6816.6766.6786.6796.672];Y=[6.6616.6616.6676.6676.664];[mu,sigma,muci,sigmaci]=normfit(X,0.1)%金球测定的估计[MU,SIGMA,MUCI,SIGMACI]=normfit(Y,0.1)%铂球测定的估计...55例分别使用金球和铂球测定引力常数...55运行后结果显示如下：mu=6.6782sigma=0.0039muci=6.67506.6813sigmaci=0.00260.0081MU=6.6640SIGMA=0.0030MUCI=6.66116.6669SIGMACI=0.00190.0071...56运行后结果显示如下：...56由上可知，金球测定的μ估计值为6.6782，置信区间为[6.6750，6.6813]；σ的估计值为0.0039，置信区间为[0.0026，0.0081]。泊球测定的μ估计值为6.6640，置信区间为[6.6611，6.6669]；σ的估计值为0.0030，置信区间为[0.0019，0.0071]...57由上可知，金球测定的μ估计值为6.6782，常用分布的参数估计函数函数名调

用

形

式函

数

说

明binofitPHAT=binofit(X,N)[PHAT,PCI]=binofit(X,N)[PHAT,PCI]=binofit(X,N,ALPHA)二项分布的概率的最大似然估计置信度为95%的参数估计和置信区间返回水平α的参数估计和置信区间poissfitLambdahat=poissfit(X)[Lambdahat,Lambdaci]=poissfit(X)[Lambdahat,Lambdaci]=poissfit(X,ALPHA)泊松分布的参数的最大似然估计置信度为95%的参数估计和置信区间返回水平α的λ参数和置信区间normfit[muhat,sigmahat,muci,sigmaci]=normfit(X)[muhat,sigmahat,muci,sigmaci]=normfit(X,ALPHA)正态分布的最大似然估计，置信度为95%返回水平α的期望、方差值和置信区间betafitPHAT=betafit(X)[PHAT,PCI]=betafit(X,ALPHA)返回β分布参数a和b的最大似然估计返回最大似然估计值和水平α的置信区间unifit[ahat,bhat]=unifit(X)[ahat,bhat,ACI,BCI]=unifit(X)[ahat,bhat,ACI,BCI]=unifit(X,ALPHA)均匀分布参数的最大似然估计置信度为95%的参数估计和置信区间返回水平α的参数估计和置信区间expfitmuhat=expfit(X)[muhat,muci]=expfit(X)[muhat,muci]=expfit(X,alpha)指数分布参数的最大似然估计置信度为95%的参数估计和置信区间返回水平α的参数估计和置信区间...58常用分布的参数估计函数函数名调用形式函数各函数返回已给数据向量X的参数最大似然估计值和置信度为（1-α）×100%的置信区间。α的默认值为0.05，即置信度为95%gamfitphat=gamfit(X)[phat,pci]=gamfit(X)[phat,pci]=gamfit(X,alpha)γ分布参数的最大似然估计置信度为95%的参数估计和置信区间返回最大似然估计值和水平α的置信区间weibfitphat=weibfit(X)[phat,pci]=weibfit(X)[phat,pci]=weibfit(X,alpha)韦伯分布参数的最大似然估计置信度为95%的参数估计和置信区间返回水平α的参数估计及其区间估计Mlephat=mle('dist',data)[phat,pci]=mle('dist',data)[phat,pci]=mle('dist',data,alpha)[phat,pci]=mle('dist',data,alpha,p1)分布函数名为dist的最大似然估计置信度为95%的参数估计和置信区间返回水平α的最大似然估计值和置信区间仅用于二项分布，pl为试验总次数...59各函数返回已给数据向量X的参数最大似然8.6假设检验检验关于分布或参数未知的总体的假设是否合理。...608.6假设检验检验关于分布或参数未知8.6.12已知，单个正态总体的均值μ的假设检验（U检验法）h=ztest(x,m,sigma)%x为正态总体的样本，m为均值μ0，sigma为标准差，显著性水平为0.05(默认值)h=ztest(x,m,sigma,alpha)%显著性水平为alpha[h,sig,ci,zval]=ztest(x,m,sigma,alpha,tail)%sig为观察值的概率，当sig为小概率时则对原假设提出质疑，ci为真正均值μ的1-alpha置信区间，zval为统计量的值。...618.6.12已知，单个正态总体的均值μ的假设检验h=...62...62例某车间用一台包装机包装葡萄糖，包得的袋装糖重是一个随机变量，它服从正态分布。当机器正常时，其均值为0.5公斤，标准差为0.015。某日开工后检验包装机是否正常，随机地抽取所包装的糖9袋，称得净重为（公斤）0.497,0.506,0.518,0.524,0.498,0.511,0.52,0.515,0.512问机器是否正常？...63例某车间用一台包装机包装葡萄糖，包得的袋装糖...64...648.6.22未知，单个正态总体的均值μ的假设检验(t检验法)h=ttest(x,m)%x为正态总体的样本，m为均值μ0，显著性水平为0.05h=ttest(x,m,alpha)%alpha为给定显著性水平[h,sig,ci]=ttest(x,m,alpha,tail)%sig为观察值的概率，当sig为小概率时则对原假设提出质疑，ci为真正均值μ的1-alpha置信区间。...658.6.22未知，单个正态总体的均值μ的假设检验h=...66...66例某种电子元件的寿命X（以小时计）服从正态分布，μ

、σ2均未知。现测得16只元件的寿命如下159280101212224379179264222362168250149260485170问是否有理由认为元件的平均寿命大于225（小时）？...67例某种电子元件的寿命X（以小时计）服从正态分布...68...688.6.3两个正态总体均值差的检验（t检验）两个正态总体方差未知但等方差时，比较两正态总体样本均值的假设检验[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置信区间。...698.6.3两个正态总体均值差的检验（t检验）两...70...70...71...71解：两个总体方差不变时，在水平下检验假设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)结果显示为：h=1sig=2.1759e-004%说明两个总体均值相等的概率很小ci=-Inf-1.9083结果表明：H=1表示在水平下，应该拒绝原假设，即认为建议的新操作方法提高了产率，因此，比原方法好。...72解：两个总体方差不变时，在水平下检验假设...728.6.4两个总体一致性的检验——秩和检验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，则不一致较明显。...738.6.4两个总体一致性的检验——秩和检验p=ra例某商店为了确定向公司A或公司B购买某种商品，将A和B公司以往的各次进货的次品率进行比较，数据如下所示，设两样本独立。问两公司的商品的质量有无显著差异。设两公司的商品的次品的密度最多只差一个平移，取α=0.05。A：7.03.59.68.16.25.110.44.02.010.5B：5.73.24.111.09.76.93.64.85.68.410.15.512.3...74例某商店为了确定向公司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];[p,h,stats]=ranksum(A,B,0.05)结果为：p=0.8041h=0stats=zval:-0.2481ranksum:116结果表明:1)两样本总体均值相等的概率为0.8041，不接近于0；2)，H=0也说明可以接受原假设，即认为两个公司的商品的质量无明显差异。...75解：设、分别为A、B两个公司的商品次品率总体的均值。则该问题8.6.5两个总体中位数相等的假设检验——符号秩检验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的正态统计量的值。...768.6.5两个总体中位数相等的假设检验——符号秩检验p例两个正态随机样本的中位数相等的假设检验x=normrnd(0,1,20,1);y=normrnd(0,2,20,1);[p,h,stats]=signrank(x,y,0.05)p=0.0251h=1stats=zval:-2.2400signedrank:45

结果表明：h=0表示X与Y的中位数之差不显著...77例两个正态随机样本的中位数相等的假设检验...778.6.6两个总体中位数相等的假设检验——符号检验p=signtest(X,Y,alpha)%X、Y为两个总体的样本，长度必须相同，alpha为显著性水平，P两个样本X和Y的中位数相等的概率，p接近于0则可对原假设质疑。[p,h]=signtest(X,Y,alpha)%h为检验结果：h=0表示X与Y的中位数之差不显著，h=1表示X与Y的中位数之差显著。[p,h,stats]=signtest(X,Y,alpha)%stats中sign为符号统计量的值...788.6.6两个总体中位数相等的假设检验——符号检验p=s例两个正态随机样本的中位数相等的假设检验X=normrnd(0,1,20,1);Y=normrnd(0,2,20,1);[p,h,stats]=signtest(X,Y,0.05)p=0.8238h=0stats=sign:9结果表明：h=0表示X与Y的中位数之差不显著...79例两个正态随机样本的中位数相等的假设检验...798.7方差分析8.7.1单因素方差分析单因素方差分析是比较两组或多组数据的均值，它返回原假设——均值相等的概率...808.7方差分析8.7.1单因素方差分析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为分析结果的构造...81p=anova1(X)%X的各列为彼此独立的样本观说明anova1函数产生两个图：标准的方差分析表图和盒图。方差分析表中有6列：第1列(source)显示：X中数据可变性的来源；第2列(SS)显示：用于每一列的平方和；第3列(df)显示：与每一种可变性来源有关的自由度；第4列(MS)显示：是SS/df的比值；第5列(F)显示：F统计量数值，它是MS的比率；第6列显示：从F累积分布中得到的概率，当F增加时，p值减少。...82说明anova1函数产生两个图：标准的方例设有3台机器，用来生产规格相同的铝合金薄板。取样测量薄板的厚度，精确至‰厘米。得结果如下：机器1：0.2360.2380.2480.2450.243机器2：0.2570.2530.2550.2540.261机器3：0.2580.2640.2590.2670.262检验各台机器所生产的薄板的厚度有无显著的差异？...83例设有3台机器，用来生产规格相同的铝合金薄解：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

...84解：...84例建筑横梁强度的研究：3000磅力量作用在一英寸的横梁上来测量横梁的挠度，钢筋横梁的测试强度是：8286798384858687；其余两种更贵的合金横梁强度测试为合金1：748278757677；合金2：797977788279]。检验这些合金强度有无明显差异？解：strength=[82867