estx matlab,概率统计的matlab求解.ppt

《概率统计的matlab求解.ppt》由会员分享，可在线阅读，更多相关《概率统计的matlab求解.ppt(76页珍藏版)》请在装配图网上搜索。

1、2020/7/31,1,第八讲 概率统计的MATLAB求解,刘北战,2020/7/31,2,随机试验、样本空间与随机事件的关系,每一个随机试验相应地有一个样本空间, 样 本空间的子集就是随机事件.,随机试验样本空间,随机事件,随机事件,基本事件,必然事件,不可能事件,两个特殊事件,1.1 随机变量及其分布,2020/7/31,3,随,量,机,变,简记为 r.v.,定义 对于定义在样本空间 上的每一个样本点 有一个实数 与之对应， 则称实值函数,在样本空间 上定义一种实值单值函数.,2020/7/31,4,离散型,定义 如果随机变量所取的可能值是有限多个或 无限可列个,并以确定的概率取这些不同的。

2、 值， 则为离散型随机变量.,随机变量,连续型,非离散型,其它,2020/7/31,5,分布律也可用表格形式表示,设 xk k1,2, 是离散型随机变量 所取 的一切可能值，称,为离散型随机变量 的概率函数或分布律.,分布律的定义,2020/7/31,6,分布函数的定义,2020/7/31,7,连续型随机变量及其概率密度的定义,连续型随机变量的分布函数在 上连续,则称 为连续型随机变量, 称 为 的概率密度 函数，简称为概率密度或密度函数 .,2020/7/31,8,超几何分布Hn,M,N,命令1Fxhygecdfx,M,N,K 功能计算超几何分布的累积概率，总共M件产品，其中次品N 件，抽。

3、取K件检查，计算发现次品不多于x件的概率FxP次品数XxFx 命令2xhygeinvp,M, N,K 功能在已知参数M、N 、 K和p的情况下计算随机量x，使得pP0次品数Xx 命令3XhygerndM,N,K,m,n 功能在已知参数M,N ,K的情况下产生m*n维符合超几何分布的随机数矩阵X,2020/7/31,9,命令4Pxhygepdfx,M, N, K 功能总共M件产品，其中次品N 件，抽取K件检查，计算发现恰好x件次品的概率PxPXx 注以后碰到命令末尾为 rnd产生随机数X； cdf产生分布函数Fx pdf产生密度函数px或分布律PxPXx inv计算xF-1p pF x,2020。

4、/7/31,10,二项分布Bn,p,命令1Fxbinocdfx,n,p 功能计算二项分布的累积概率FxPXxFx 命令2xbinoinvy, n,p 功能计算随机量x，使得yPXx 命令3Xbinorndn,p,M,N 功能产生M*N维符合二项分布的随机数矩阵X 命令4Pxbinopdfx,n, p 功能计算试验中事件恰好发生x次的概率PxPXx,2020/7/31,11,泊松分布XP,命令1Fxpoisscdfx,lambda 功能计算累积概率FxPXxFx 命令2xpoissinvp, lambda 功能计算随机量x，使得pPXx 命令3Xpoissrndlambda,M,N 功能产生M*。

5、N维随机数矩阵X 命令4Pxpoisspdfx,lambda 功能计算概率PxPXx,2020/7/31,12,正态分布XN,2,命令1Fxnormcdfx, mu,sigma 功能计算累积概率FxPXxFx 命令2xnorminvp, mu,sigma 功能计算随机量x，使得pPXx 命令3Xnormrndmu,sigma,M,N 功能产生M*N维随机数矩阵X 命令4Pxnormpdfx, mu,sigma 功能计算分布密度px在x的值,2020/7/31,13,指数分布Xexp,命令1Fxexpcdfx, lambda 功能计算累积概率FxPXxFx 命令2xexpinvp, lambda。

6、 功能计算随机量x，使得pPXx 命令3Xexprndlambda,M,N 功能产生M*N维随机数矩阵X 命令4Pxexppdfx, lambda 功能计算分布密度px在x的值,2020/7/31,14,均匀分布XUa,b,命令1Fxunifcdfx, a,b 功能计算累积概率FxPXxFx 命令2xunifinvp, a,b 功能计算随机量x，使得pPXx 命令3Xunifrnda,b,M,N 功能产生M*N维随机数矩阵X 命令4Pxunifpdfx, a,b 功能计算分布密度px在x的值,2020/7/31,15,分布,命令gamcdfx, a, lambda, gaminvp, a, l。

7、ambda gampdfx, a,lambda, gamrnda, lambda,m,n,2020/7/31,16,卡方分布,命令chi2cdfx, k, chi2invp, k,chi2pdfx, k chi2rndk,m,n,2020/7/31,17,t分布,命令tcdfx, k, tinvp, k,tpdfx, k trndk,m,n,2020/7/31,18,F分布,命令fcdfx, p,q, finvF,p,q,fpdfx, p,q frndp,q,m,n,2020/7/31,19,例1.1某人向空中抛硬币100次，落下为正面的概率为0.5。这100次中正面向上的次数记为X 1试计算。

8、x45的概率和x45的概率； 2绘制分布函数图象和分布列图象。,程序clear; pxbinopdf45,100,0.5 计算x45的概率 px 0.0485 fxbinocdf45,100,0.5 计算x45的概率 fx 0.1841 x1100;p1binocdfx,100,0.5;plotx,p1,; title分布函数图,2020/7/31,20,p2binopdfx,100,0.5;plotx,p2,*r；title概率分布图),2020/7/31,21,例1.2设XN2,0.25 1 求概率P1X2.5； 2绘制分布函数图象和分布密度图象; 3画出区间1.5,1.9上的分布密。

9、度曲线下方区域。,程序1pnormcdf2.5,2,0.5- normcdf1,2,0.5 p 0.8186 2 x00.14;pxnormpdfx,2,0.5; fx normcdfx,2,0.5; plotx,px,b;hold on; plotx,fx,*r;legend正态分布函数,正态分布密度; 3 specs1.5,1.9; ppnormspecspecs,2,0.5,2020/7/31,222020/7/31,23,1.2 随机变量的数字特征,随机变量的数学期望,1.统计数据的平均值-YmeanX 功能当X为向量时，输出一个平均数；当X为矩阵时，输出为行向量，对应于矩阵每列的平。

10、均值；因此计算矩阵所有数的平均值，应用嵌套meanmeanX或mmeanX 与此类似的有求和sum,最大max,最小min等 2.离散型随机变量的期望EXsumX.*P 功能计算随机值向量X与对应概率向量P的乘积之和 3.连续型随机变量的期望EXintx*fx,x,a,b 功能用积分计算期望,2020/7/31,24,例1.4设随机变量X的分布列，求期望。,程序clear; x-1,0,2,3; p1/8,1/4,3/8,1/4; EXsumx.*p 1.3750,2020/7/31,25,例1.5设随机变量X的分布密度为 且EX3/5，求常数a,b的值。,程序clear;syms a b x。

11、;fxab*x2; EXintx*fx,x,0,1 EX1/4*b1/2*a Fintfx,x,0,1 Fa1/3*b f1EX-3/5;f2F-1; a,bsolvef1,f2 a3/5,b6/5,2020/7/31,26,例1.6设随机变量X的分布密度为 求随机变量Y|X|的期望。,程序clear;syms x; fx10.5*expx; fx20.5*exp-x; EYint-x*fx1,x,-inf,0 intx*fx2,x,0, inf EY 1,2020/7/31,27,随机变量的方差,1.统计数据的方差-DvarX,flag 功能当X为向量时，输出一个标量；当X为矩阵时，输出为。

12、行向量，对应于矩阵每列的方差值；因此计算矩阵所有数的方差值，应用嵌套varvarX或varX flag0时，计算 X的样本方差 flag1时，计算 X的简单方差 2.统计数据的标准差-SstdX,flag 功能flag的解释同上 3. 一般随机变量的方差DXEX2-EX2 功能用积分计算,2020/7/31,28,例1.7设随机变量X的分布密度为 求随机变量X的期望和方差。,程序clear;syms x;fx2/pi*cosx2; EXintx*fx,x,-pi/2,pi/2 E2Xintx2*fx,x,-pi/2,pi/2 DXE2X-EX2,2020/7/31,29,常见分布的期望和方差,。

13、1.二项分布-E,Dbinostatn,p 说明n,p可以是标量,向量,矩阵,则E,D是对应的标量,向量,矩阵 2.超几何分布-E,DhygestatM,N,K 3.泊松分布-E,Dpoissstatlambda 4.均匀分布-E,Dunifstata,b 5.指数分布-E,Dexpstatlambda 6.正态分布-E,Dnormstatmu,sigma 其他gamstat,tstat,fstat,chi2stat等等,2020/7/31,30,协方差与相关系数的计算,1.随机变量的协方差-covX,YEX-EXY-EY 2.随机变量的相关系数-covX,Y/sqrtDX*DY 3.统计数据。

14、的协方差 covX-当X为向量时,covXvarX;当X为矩阵时,结果为X的协方差矩阵.对角线是X每列的方差,Xij为X的第i列和第j列的协方差值。 covX,Y-计算向量X和Y的协方差值 4.统计数据的相关系数 corrcoefX,corrcoefX,Y-说明与用法与cov相同,2020/7/31,31,矩的计算,1.随机变量的k阶原点矩-EXk 2.随机变量的k阶中心矩-EX-EXk 3.统计数据的k阶原点矩-AksumX.k/lengthX 4.统计数据的k阶中心矩-BkmomentX,k2020/7/31,32,1.3 参数估计,常用分布的参数估计,1.正态分布的参数估计 格式muha。

15、t,sigmahat,muci,sigmacinormfitX,alpha 功能数组X服从正态分布,给定显著水平alpha,缺省时为0.05,前二项给出点估计，后二项给出区间估计。X为矩阵时,针对列进行计算。 2.二项分布的参数估计n已知,p未知 格式phat,pcibinofitX,n,alpha 3.泊松分布的参数估计 格式lbdhat,lbdcipoissfitX, alpha 4.均匀分布的参数估计 格式ahat,bhat,aci,bciunifitX,alpha,2020/7/31,33,5.指数分布的参数估计 格式lbdhat, lbdciexpfitX,alpha 6.通用命令m。

16、le 格式输出参数项mle分布函数名,X,alpha ,N 说明分布函数名有bino二项,geo几何,hyge超几何 poiss泊松,uinf均匀,unid离散均匀,exp指数 norm正态,tT分布,fF分布,beta贝塔,gam伽吗 N仅当二项分布时需要，其他没有(因为二项分布的参数估计需要n已知)。,2020/7/31,34,例1.8设生成一组均值为15,方差为2.52的正态分布的随机数据，然后对这组数据进行置信度97的参数估计。,程序clear; wnormrnd15,2.5,50,1; 或w152.5*randn50,1; alpha0.03; mh,sh,mc,scnormfitw。

17、,alpha 运行一次mh15.1076 sh2.4038 mc14.347815.8674 sc1.97093.0703,2020/7/31,35,例1.9设从一大批产品中抽取100个产品，经检验知有60个一级品，求这批产品的一级品率置信度95。,程序clear; alpha0.05;N100;X60; Ph,Pcmlebino,X,alpha,N 运行一次Ph0.6000 Pc0.49720.69672020/7/31,36,1.参数检验如果观测的分布函数类型已知，这时构造出的 统计量依赖于总体的分布函数，这种检验称为参数检验. 参数检验的目的往往是对总体的参数及其有关性质作出明 确的判断。

18、.,对总体X的分布或分布参数作某种假设，根据抽取的样本观察值，运用数理统计的分析方法，检验这种假设是否正确，从而决定接受假设或拒绝假设.,2.非参数检验如果观测的分布函数类型未知，所检验的假设 不是对某个参数作出明确的判断，此时必须要求构造出的 检验统计量的分布函数不依赖于观测值的分布函数类型， 这种检验叫非参数检验. 如要求判断总体分布类型的检验就是非参数检验.,1.4 假设检验,2020/7/31,37,假设检验的一般步骤是,2020/7/31,38,(一)单个正态总体均值检验,1.4.1参数检验,2020/7/31,39,2020/7/31,40,(二)单个正态总体方差检验,2020/7。

19、/31,41,(三)两个正态总体均值检验,2020/7/31,42,(四)两个正态总体方差检验,2020/7/31,43,单正态总体均值的假设检验,1.方差已知u检验或z检验 格式H,P,ci,ZvalztestX,Mu,sigma,alpha,tail 功能对正态分布总体的采样X进行Z检验，判断采样的均值在已知的标准差sigma下是否等于假设值Mu；给定显著水平alpha,缺省时为0.05；tail是假设的备选项即备择假设 ，有三个值tail0是默认值，可省略,说明备选项为均值不等于 M； tail1,说明备选项为均值大于 M； tail-1,说明备选项为均值小于 M。H0说明接受原假设，H。

20、1拒绝原假设；P为假设成立的概率，P值非常小时对假设置疑；ci给出均值的置信；Zval给出统计量的值。,2020/7/31,44,例1.10某面粉厂的包装车间包装面粉，每袋面粉的重量服从正态分布，机器正常运转时每袋面粉重量的均值为50kg，标准差1。某日随机的抽取了9袋，重量分别为 49.7，50.6，51.8，52.4，49.8，51.1，52，51.5，51.2 机器运转是否正常,程序clear;x49.7,50.6,51.8,52.4,49.8,51.1,52,51.5,51.2; sigma1;mu50; h,p,ci,zztestx,mu,sigma 结果h1 拒绝原假设即认为机器不。

21、正常 p 7.6083e-004 p0.00076083很小，对假设置疑 ci 50.4689 51.7755 均值偏高 z 3.3667,2020/7/31,45,单正态总体均值的假设检验,2.方差未知t检验 格式H,P,ci,statsttestX,Mu, alpha,tail 功能对正态分布总体的采样X进行t检验， 对H,Mu,alpha,tail,P,ci的解释同上；stats是个结构，包含三个元素tstat统计值、df自由度和sd样本标准差。,例1.11某灯泡厂出厂的标准是寿命不少于2000小时，现随机的从该厂生产的一批灯泡中抽取了20只，寿命分别为 1558，1627，2101，。

22、1786，1921，1843，1655，1675 1935，1573，2023，1968，1606，1751，1511，1247 2076，1685，1905，1881 假设灯泡的寿命服从正态分布问这批灯泡是否达到了出厂标准 a0.01,2020/7/31,46,原假设H0mu2000 备择假设H1mu2000 程序clear; x1558,1627,2101,1786,1921,1843,1655,1675,1935,1573,2023,1968,1606,1751,1511,1247,2076,1685,1905,1881; alpha0.01;mu2000; h,p,ci,statstt。

23、estx,mu,alpha,-1 结果h1 拒绝原假设即认为不符合出厂标准 p 5.9824e-005 p很小，对假设置疑 ci 1.0e003 * -Inf 1.8895 均值偏低 stats tstat - 4.8176 df 19 sd 216.8973,2020/7/31,47,双正态总体均值的假设检验,比较两个方差相等的正态总体的均值是否相等T检验 格式H,P,ci,statsttest2X,Y, alpha,tail 功能对两个正态分布总体的采样X、Y进行T检验，对H, P,alpha的解释同上； tail是假设的备选项即备择假设 ，有三个值tail0是默认值，可省略,说明备选。

24、项为均值不相等； tail1,说明备选项为X的均值大于Y的均值； tail-1,说明备选项为X的均值小于Y的均值。ci给出均值差的置信区间； stats是个结构，包含三个元素tstat统计值、df自由度和sd标准差Sw。,2020/7/31,48,原假设H0mu2mu1 备择假设H1mu2mu1 程序clear; x2461,2404,2407,2439,2394,2401,2543,2463,2392,2458; y2496,2485,2538,2596,2556,2582,2494,2528,2537,2492; alpha0.01;h,p,ci,stttest2x,y,alpha,-1 。

25、结果h1 拒绝原假设即认为寿命未提高 p 6.3361e-005 p很小，对假设置疑 ci -Inf -44.6944 st tstat -4.8567 df 18 sd 43.3705,例1.12某灯泡厂在采用一项新工艺前后，分别抽取了10只进行寿命试验，寿命分别为 旧灯泡2461,2404,2407,2439,2394,2401,2543,2463,2392,2458 新灯泡2496,2485,2538,2596,2556,2582,2494,2528,2537,2492 假设灯泡的寿命服从正态分布，采用新工艺前后方差相等，能否认为采用新工艺后，灯泡的寿命提高了a0.01,2020/7。

26、/31,49,1.两个总体一致性的假设检验,比较两个不知道确切分布的总体均值是否相等 格式P, H,statsranksumX,Y, alpha 功能对两个总体的采样X、Y进行检验，对H, P,alpha的解释同上； stats是个结构，包含二个元素zval均值差的正态统计值和ranksum统计的秩和值。 秩和检验详见word文档,1.4.2 非参数检验,2020/7/31,50,程序clear;x33.592,33.862,33.751,33.673,33.847, 33.778,33.631,33.911,33.785,33.928; y34.221,33.947,33.856,34.03。

27、9,34.000,33.924,34.125, 34.273,33.968,33.923; alpha0.05;p, h,stranksumx,y,alpha 结果 p 7.6854e-004 p很小，对假设置疑 h1 拒绝原假设即认为直径没有显著不同 st zval -3.3639 ranksum 60,例1.13两台机床加工同一种轴，抽样测量产品的直径mm 机床甲33.592,33.862,33.751,33.673,33.847,33.778, 33.631,33.911,33.785,33.928 机床乙34.221,33.947,33.856,34.039,34.000,33.9。

28、24, 34.125,34.273,33.968,33.923 在a0.05下能否认为两台机床加工的直径没有显著不同,2020/7/31,51,2.两个样本具有相同类型连续分布的假设检验,检验两个样本是否具有相同类型的连续分布 格式 H ,P, ksstatkstest2X,Y, alpha,tail 功能对两个总体的采样X、Y进行检验，对H, P,alpha的解释同上； tail是假设的备选项即备择假设 ，有三个值tail0是默认值，可省略,说明备选项为不相等； tail1,说明备选项为大于； tail-1,说明备选项为小于。 ksstat表示测试统计量的值。,2020/7/31,52,程序。

29、clear;xrandn1,10;yrandn1,10x; h, p,ksstatkstest2x,y 结果 h0 接受原假设即认为两样本具有相同类型的连续分布 p 0.6751 表示假设成立的概率为0.6751 Ksstat 0.3000,例1.14,2020/7/31,53,3.单个样本分布的 Kolmogorov-Smirnov 测试,指定累积分布函数为cdf的测试,格式 H ,P, ksstat,cvkstest X,cdf, alpha,tail 功能对采样X进行检验是否服从名为cdf类型的连续累积概率 分布，cdf缺省为 ,默认为标准正态分布,声明格式为两 个相同长度的列向量采。

30、样和采样对应的分布函数； H0，接受原假设， H1拒绝原假设，接受备择假设； P为原假设成立的概率； ksstat为测试统计量的值； cv为是否接受假设的临界值。,2020/7/31,54,产生100个二项分布随机数，测试该随机数服从的分布 程序clear; xbinornd100,0.2,100,1; H,p,ksstat,cvkstestx,x binocdfx,100,0.2,0.05 测试是否服从二项分布 结果 H 0 接受原假设即认为随机数服从二项分布 p 0.1680 表示假设成立的概率为0.1680 ksstat 0.1098 cv 0.1340,例1.15,2020/7。

31、/31,55,程序clear; mu1;sigma 2; xnormrndmu,sigma,20,1;alpha0.01;a0;b1; h, p,ksstat,cvkstest x,x,unifcdfx,a,b,alpha 结果 h 1 拒绝原假设即认为随机数不服从U0,1 p 1.4115e-006 表示假设成立的概率很小，拒绝 ksstat 0.5764 cv 0.3524,例1.16,2020/7/31,56,4.正态分布拟合优度测试,命令1 H ,P, jbstat,cvjbtest X, alpha 功能对采样X进行检验是否服从正态分布，对H, P,alpha的解释同上； jbst。

32、at表示测试统计量的值;cv为是否拒绝假设的临界值。适合大样本 命令2 H ,P, lstat,cvlillietest X, alpha 功能对采样X进行检验是否服从正态分布，对H, P,alpha的解释同上； lstat表示测试统计量的值;cv为是否拒绝假设的临界值。适合小样本,2020/7/31,57,程序clear; m1ones1,11*2.55; m2ones1,12*2.65; m3ones1,17*2.75; m4ones1,19*2.85; m5ones1,26*2.95; m6ones1,24*3.05; m7ones1,22*3.15; m8ones1,19*3.25; 。

33、m9ones1,13*3.35; Mm1,m2,m3,m4,m5,m6,m7,m8,m9; h, p,lst,cvlillietest M histM,例1.17从一批零件中随机抽取一组样品，下面是零件样品直径的统计表。在显著水平a0.05下能否认为这批零件的直径服从正态分布绘出统计数据的直方图。,2020/7/31,58,结果 h1 拒绝原假设 即认为直径不服从正态分布 p Nan 表示假设成立 的概率很小 lst 0.1062 cv 0.0694 测试统计值 大于临界值也表明应拒绝 histM,n -绘制向量M的频数直方图,n定义条方的数目,默认为10,h, p,jbst,cvjbte。

34、st M 结果h 1；p 0.0368；jbst 6.6043；cv 5.9915,2020/7/31,59,5. 概率纸检验法,概率纸是一种判断总体分布的简便工具.使用它们，可以很快地判断总体分布的类型.概率纸的种类很多.,2020/7/31,60,Matlab工具箱提供了两种分布的概率纸检验命令,(1)h normplotx,(2)h weibplotx,此命令显示数据矩阵x的正态概率图.如果数据来自于正态分布，则图形显示出直线性形态.而其它概率分布函数显示出曲线形态.,此命令显示数据矩阵x的Weibull概率图.如果数据来自于Weibull分布，则图形将显示出直线性形态.而其它。

35、概率分布函数将显示出曲线形态.,2020/7/31,61,程序 Xnormrnd0,1,50,1; normplotX 结果,例1.18,2020/7/31,62,程序 r weibrnd1.2,1.5,100,1; weibplotr 结果,例1.191,2020/7/31,63,程序 r weibrnd1.2,1.5,1000,1; weibplotr 结果,例1.192,2020/7/31,64,例1.20 一道工序用自动化车床连续加工某种零件，由于刀具损坏等会出现故障.故障是完全随机的，并假定生产任一零件时出现故障机会均相同.工作人员是通过检查零件来确定工序是否出现故障的.现积累有。

36、100次故障纪录，故障出现时该刀具完成的零件数如下 459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659 775 859 755 49 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 48。

37、4 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851 试观察该刀具出现故障时完成的零件数属于哪种分布.,2020/7/31,65,解 1、数据输入,MATLAB(liti201),2、作频数直方图 histx,10,3、分布的正态性检验 normplotx,4、参数估计 muhat,sigmahat,muci,sigmaci normfitx,(看起来刀具寿命服从正态分布),(刀具寿命近似服从正态。

38、分布),估计出该刀具的均值为594，方差204，均值的0.95置信区间为 553.4962，634.5038，方差的0.95置信区间为 179.2276，237.1329.,MATLAB(liti204),MATLAB(liti202),MATLAB(liti203),2020/7/31,66,5、假设检验,MATLAB(liti205),已知刀具的寿命服从正态分布，现在方差未知的情况下，检验其均值 m 是否等于594.,结果h 0，p 1，ci 553.4962，634.5038.,检验结果 1. h0, 表示不拒绝零假设. 说明提出 的假设寿命均值594是合理的. 2. 95的置信区间。

39、为553.5，634.5, 它 完全包括594, 且精度很高. 3. p1, 远超过0.5, 不能拒绝零假设,2020/7/31,67,1.5 方差分析ANOVA-analysis of variance 详见word文档,单因素方差分析,命令 P, anovatab,statsanova1 X,group, displayopt 功能比较多组数据的均值，返回这些均值相等的概率，从而判断因素对结果是否有显著影响。X为输入数据，列向量表示相互独立的样本观测值，具有相同长度；P为X的各列均值相等的概率，P越小，则质疑原假设(即均值不相等)，表示因素的影响显著；group是与X对应的字符或字符串数组。

40、，用来声明X每一列中数据的名字或意义，可以省略；displayopt表示参数on表示显示图，off表示隐藏图；anovatab返回方差分析表；stats返回一个附加的统计数据结构。,2020/7/31,68,1.21将同一批同种牌号丝袜在不同温度下作弹力试验，得到数据表试检验温度对弹力有无显著影响。0.05,2020/7/31,69,程序clear;X4.3,6.1,10.0,6.5,9.3,9.5; 7.8,7.3,4.8,8.3,8.7,8.8; 3.2,4.2,5.4,8.6,7.2,11.4; 6.5,4.1,9.6,8.2,10.1,7.8; p,tab,statsanova1X, 。

41、,on 结果p 0.0214 p很小，拒绝原假设， 即温度对弹力有影响 tab Source SS df MS F ProbF Columns 55.5471 5 11.1094 3.5254 0.0214 Error 56.7225 18 3.1512 Total 112.2696 23 stats gnames 6x1 char n 4 4 4 4 4 4 source anova1 means 5.4500 5.4250 7.4500 7.9000 8.8250 9.3750 df 18 s 1.7752 总体标准差的无偏估计,2020/7/31,70,2020/7/31,71,1.。

42、22将四种工艺下生产的灯泡进行寿命测试，得到数据表试检验工艺对寿命有无显著影响。0.05,2020/7/31,72,程序clear;X1620,1670,1700,1750,1800,1580,1600,1640,1720,1460,1540,1620,1500,1550,1610,1680; group1,1,1,1,1,2,2,2,2,3,3,3,4,4,4,4; p,tab,statsanova1X,group,on 结果p 0.0331 p很小，拒绝原假设， 即工艺对寿命有影响 tab Source SS df MS F ProbF Groups 62820 3 20940 4.0。

43、608 0.0331 Error 61880 12 5.1567e003 Total 124700 15 stats gnames 4x1 cell n 5 4 3 4 source anova1 means 1708 1635 1540 1585 df 12 s 71.8099,2020/7/31,73,双因素方差分析,命令 P, anovatabanova2 X,reps, displayopt 功能判断因素对结果是否有显著影响。X为输入数据，列向量表示因素1的差异，行向量表示2的差异；P是概率向量，P越小，则质疑原假设，表示因素的影响显著；reps声明每一状态下的试验次数；displa。

44、yopt表示参数on表示显示图，off表示隐藏图；anovatab返回方差分析表。,2020/7/31,74,1.23有3个工人分别在4台机器上加工某种零件，工作的3天中日产量列表如下试检验操作工人的技术水平有无显著差异机器性能有无显著差异交互作用的影响是否显著0.05,2020/7/31,75,程序clear; A115,15,17,19,19,16,16,18,21; A217,17,17,15,15,15,19,22,22; A315,17,16,18,17,16,18,18,18; A418,20,22,15,16,17,17,17,17; XA1,A2,A3,A4; reps3; p,tabanova2X,reps,on,2020/7/31,76,解释 操作工人的技术水平无显著差异p0.6645 机器性能有显著差异p0.0023 交互作用的影响显著p0.0002。