matlab数理统计综合实验--excel

Excel Excel基本操ExcelA1A1(a1)选取单元格范围(矩形区域鼠标指针移到终点（右下角）位置，先按下Shift键不放，而后点击鼠标左键．入公式，公式输入完毕后按Enter键即可确认．.如果单击了“编辑公式”按钮或“粘贴函数”按钮，Excel将自动一个等号．提示：(1)通过先选定一个区域，再键入公式，然后按CTRL+ENTER组合键，可乘法等运算．公式可以同一工作表中的其它单元格、同一工作簿不同工作表中的单B4中的数值加上25，再除以单元格D5、E5和F5中数值的和．在默认状态下，Excel从等号（＝）开始，从左到右计算公式．可以通过修改公式语法来控制计算的顺序．例如，公式=5+2*31123（6），然5Excel先计算乘法再计算加法；可以使用圆括号来改变语法，圆括号内的内容将首先被计算．公式=(5+2)*32152，再用其结果乘以3．称为从属单元格，它的值依赖于被单元格的值．只要被单元格做了修改，包含公式的单元格也就随之修改．例如，公式“=B15*5”将单元格B15中的数值乘以在默认状态下，Excel使用A1类型．这种类型用字母标志列（从A到IV，共256列，用数字标志行（从1到65536如果要单元格，请顺序输入列字母和行数字．例如，D50了列D和行50交叉处的单元格．如果要单元格区域， A10A102015BEA10E205510HHJExcel包含许多预定义的，或称内置的公式，它们被叫做函数．函数可以进行简单的ENTERs２=(xins２=(xin均值（μ）：分布的均值；σ得到的结果为0.158655，即P{X≤500}=0.158655．ExcelNORMSDIST(x)，及标准正态分布的反函数NORMSINV(概率)．范例：已知Ｘ～Ｎ(0,1),计算(2=P{X<2}0.97725，即(2)

)若求临界值uα，则使用公式＝NORMSINV(1-tt（T>1.711得到0.05，即P（T>1.711）=0.05．t2,得到1.812462，即t005(10)=1.812462．F

范例：设X服从度n1=5，n2=15的F分布，求P(X>2.9)的值．输入公F范例：已知随量X服从F(9,9)分布，临界值α=0.05，求其上侧0.05分位025（9,9F0975（9,9,0.95，即1F(5.226)=0.95F(5.226)0得到值5.2262(1202若求临界值χα(n)，则使用公式＝CHIINV(α,2范例：设Ｘ服从参数为4的泊松分布，计算P{X=6}及P{X≤6}．输入公式在下面的实验中，还将碰到一些其它函数，例如：COUNT，开平方函数SQRT，和函数SUM，等等．关于这些函数的具体用法，可以查看Excel的关 区间估计置信区间为x

x

nn nn 例1随机从一批苗木中抽取16株，测得其高度（单位：m）为 =average(b3:g5)-normsinv(1-0.5*)*=average(b3:g5)+normsinv(1-0.5*)*标准差是具体的数值而不是符号．本例中，

0.01计算结果为（1.148225,1158025置信区间为x

(n

n xn

(n Sn． Sn．例2 同例1，但未1.172555未知时2 置信区间为(n1)s2 (n1) ．(n (n1) 例 试求总体方差20.9的置信区间(设总体为正计算置信区间下限：在单元格C14中输入公式=(C9-1)*C11/C12 11

xyt(nn 例4 甲 乙 （α＝0.05）在单元格B11输入公式＝AVERAGE(A2:A6)，在单元格B12中输入公式和μ,μ未知时方差比σ２/σ２ 22

1F(n1,n1) 12 (n1,n1)ss ss =0.5419，S=0.6065．设σ2和σ2分别是A、 所测量的数据总体（设为正态分布）的方差．求方差比σ2σ20.95 和60株林木组成样本．样22.3,21.2,19.2,16.6,23.1,23.9,24.8,26.4,26.6,24.8,23.9,23.2,23.3,21.4,19.8,20.0,21.5,18.7,22.4,26.6,23.9,24.8,18.8,27.1,20.6,25.0,22.5,23.5,23.9,25.3,22.6,21.5,20.6,25.8,24.0,23.5,22.6,21.8,20.8,19.5,20.9,22.1,22.7,23.6,24.5,21.0,21.3,22.4,18.7,21.3,15.4,22.9,17.8,21.7,19.1,20.3,已知某种玉米的产量服从正态分布，现有种植该玉米的两个实验区，各分为10第一实验区： 第二实验区： 假设检验Excel进行假设检验的方法和操作步骤；(3)能够利用试验结果的信息，对所关心的σ2已知时μU检800，8，780，90，82（k后1000m2产量Ｘ也服从正态分布，试问：σ2未知时的t检2某一引擎制造商新生产某一种引擎，将生产的引擎装入汽车内进行速度测试， 计算样本均值：在单元格D8输入公式计算标准差：在单元格D9输入公式在单元格D11输入t检验值计算公式=(D8-250)/(D9/(SQRT(D10))已知t统计量的度为(n-1)=20-1=19，域为t>t0 =2.093．由上面计算得到 例 第1组： 第2组： 1的区域”输入(7)在“输出区域”输入D1；(8)在显著水平“α”框，输入0.05；(9)在“假设平均差”窗口输入0；0.210（ 1的区域”输入入0． H：σ２=σ２；H 1的区域”输入(6在“变量2B2:B9．(7)在显著水平“α”框，输入0.025．(8)在“输出区域”框输F2.35049小于“F5.118579

k(nnp2

0其中２为上述统计量(1的值．在应用中，可根据临界概率000k(nnp

p0，利用函数CHIINV(p0,k1)inp 6设总体Ｘ120x209，样本方s=42.77及下表．试检验Ｘ是否服从正态分布（α＝0051234567867986∑建立如下图所示工作表，输入区间(A2:A10)，端点值（B2:B10），实测频数的值是＋∞时，为了便于处理，可输入一个很大的数（本例取10000）代替＋∞．则P{a<X≤b}=F(b)-F(a),因此，事件{a<X≤b}发生的理论频数为n(F(b)-F(a)).将计算的理论频数值放入D在D3输入=120*(NORMDIST(B3,209,6.539877675,TRUE)类似地，可算出D4至D10应用小技巧：D4D10D3单元格，单击鼠标右键弹出中选择＂粘贴＂，即可得到D4到D10的值．2个(,2)=

0E13输入＝CHIINV(E12,8),20值1.104413小于临界值12.5918，故接受原假设，认为X服从正态分布．某春小麦良种千粒重μ=34克，方差σ2=1.968个小为防止某种害虫而将一种施入土中，规定经三年后土壤中5ppm以上浓 甲种：24.3，17.4，23.7，20.8，21.3 乙种：20.2，16.9，16.7，18.2 （123456频数 6210

方差分析13个剂量进行试验，加上对照（处重1（对照2341234A7结果分析：Fcrit=3.4903α＝0.05F统计量临界值，F=1.305047F统计量的例2将土质基本相同的一块耕地分成均等的五个地块，每块又分成均等的四个小：kg．品(6)在“输出区域”输入A7指定显著水平“α”为“0.05”或“0.01例3 (6)在“输出区域”输入A11．行．(8)指定显著水平“α”为“005”．本例假设：HA：因素A对试验结果无显著影HBB对试验结果无显著影MSAB=294.7821,误差=236.95,MSe=19.74583,St=2638.293123435675身高128.1，134.1，133.1150.3，147.9，136.8140.6，143.1，144.51234操作甲乙丙§5回归分析实利用Excel11950～1975年的工农业总产值（X）与货运周转量（Y）的历X与Y间的关本Y数据存放在B1:B27，其中B1存标记Ｙ．B27．(6)在“输入Ｘ区域”框输入A1：95%(10)在“输出区域”框，确定单元格利用EXCEL2 69 解ExcelX数据输入在A2：D16Y数据输入在E2：E16．在“输入Ｙ区域”框输E2：E16；(9)在“输出区域”框，确定单元格G1；关于Ｘ2，由于P{t>4.369}=0.0014，故Ｘ2是显著关于Ｘ3P{t>0.413}=0.68896，故Ｘ3是不显著的；4气体在容器中被吸引的比率YX1和吸收液体的蒸汽压力X2有关，其数学模型为Y=A＋B1X1＋B2X2，测得试验数据为：YYY关于X1、X2§6数据分析综12个苹果园下表）前三列为为单株产量（kg），最后一列为土壤有机质含量%因变量操作步骤：于本实为数据析综实验，不出具体操步骤．分发挥自己的科研潜力，充分利用自己的数学知识和农业科学知识，分析农业科研问．α=0.05B1、B2B3三个水平下的产量观测数据分别来自正态总体N（μ1，σ12、2（μ2，σ2）N（μ3，σ32，分别进行任意两个总体均值差的t检验和方差齐性的212个实验作为因素A12B3个水平，分别以三种施肥方式下的单株产量作为自变量X1，X2，X3，以土壤有机质含量作为因变量Y，即选择XC2:D13，YE2:E13进行多元线性value=0.938941048X3仍不显著（tStat=0.531040343,P-value=0.608244525E2:E13回归分析，结果表明，回归方程和方程中包含的自变量（仅有X2）均极显著．某地区部门对玉米的4个新品种B1、B2、B3、B4进行产量和环境适应性实验，每个品种分别在该地区的沙壤土地（A1、粘壤土地（A2）和中壤土地（A3）上进X（kg查数据如表2．试对蔬菜需求量的影响因素进行分析．表 品种品种品种品种沙壤粘壤中壤表 销量9§ 数理7.1基命令窗口；左边的两个窗口分别为启动平台（LaunchPad）和命令历史（Commandhistory）窗口，也可以切换到工作空间（Workspace）和当前浏览器窗口（CurrentDirectory按F5键执行程序．a=[123;456;78或a=[1,2,3;4,5,6;a123456789a=2345678可以先将一个矩阵输入到Excel工作表，即先建立数据文件．然后打开“file”—“ImportDataExcel工作表，将数据导入到的工作空间中去．a(i,j)．x=[4/3format1.3333formatshort1.3333e+000formatshort1.3333formatformatlongformatlongformat4/3 的缺省格式为format窗口依次执行a=1/3，formatlong，a这三个命令，结果如下：>>a>>formata=若最大的元素大于1000或小于0.001，则显示short或long格式时会加上一个比例．在命令产生100*100的魔方矩阵，但并不在屏幕上显示．A=s=1–1/2+1/3–1/4+1/5–1/6+1/7...–1/8+1/9–1/10+1/11– savejune10.mat．也可只保存工作间中的部分变量值，例如要保存变量x,y,z到文件june10.mat，使用命令savejune10xyload要想将变量ab，可使用命令clear' +加 -减 *乘 ^幂 /左除 \右除对于两个标量a,b来说，a/b=b\a=a÷b，而对于两个矩阵a,b来说，用符号“x,yxyx的元素分别除以y的对应元素得到的数组．z=x./z(ij)=x(ij)/y(ij)．xyyxxy，但运算是在对应元素间进矩阵的乘方用“^”符号表示，a^p的意思是ap则z=x.^y=[x1^y1,x2^y2,…,z=x.^c=[x1^c,x2^c,…,y=poisscdf(x,)y=expcdf(x,正态分布N(μσ2)y=normcdf(x,μ,χ2(n)分布的分布函数)x=expinv(1-α,正态分布N(μσ2)αF(n1,n2)分布的上侧αx是单个向量（可以是行向量，也可以是列向量mx的均值，x是矩阵，则返回结果m是行向量，它包含x的每列数据的均值．即若 x

x1kx2kxx

xn

xnkn x1n n

,k)则 [x1,x2 ,xk]，其中

ij( 准化处理，其中n为数据的长度．若要运用n进行标准化，可使用var(x,1)格式．Cxx是矩阵，则返回各列数据的协方差构成的协方差矩阵．cov(x)运用n-1进行标准化处理．,素C(i,j)的关系为C(i,C(iC(i,C(i,i)C(j, 语法：hztestx,mh=ztest(x,m,sigma,[h,sig,ci,zval]=ztest(x,m,sigma,alpha,hztestx,msigma)0.05Usigma态总体的均值是否等于m．即总体X~N(,sigma2)，样本x(x,x ,x)来自 X

H0:mhztestx,msigmaalpha)[h,sig,ci,zval]=ztest(x,m,sigma,alpha,tail)可以通过制定tail的值来控择假设H1．tail的取值及意义为：sigU

x

tail0

sigPuUtail1sigPuUtail1

sigPuU．函数：ttest语法：httestx,h=ttest(x,m,[h,sig,ci]=ttest(x,m,alpha,)hztestx,malpha)给出了显著性水平控制参数[h,sig,ci]=ttest(x,m,alpha,tail)可以通过制定tail的值来控择假设H1．tail的取值及意义与ztest函数一致．xs 返回值sig是t(n1)分布以统计量T s tail0时sigPtTtail