正交试验结果的方差分析技巧

4.1.8正交试验结果的方差 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 计算公式和项目试验指标的加和值i-i，试验指标的平均值与 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 4-13一样,Ij”水平所对应的试验指标的数值之和IIj――“2”水平所对应的试验指标的数值之和⑶……kj――同一水平出现的次数。等于试验的次数除以第j列的水平数.1j/kj――“水平所对应的试验指标的平均”11j/kj――“2”水平所对应的试验指标的平均值……以上各项的计算方法，与“极差法”同，见4.1.7节偏差平方和■+AA(4-1)(4-2)(4-3)fj――自由度.fj第j列的水平数一1.V――方差.vj=S/fjve—误差列的方差Fj方差之比(4-4)**W查F分布数值表(见附录6),做显著性检验。显著性检验结果的具体表示方法与第3章相同。总的偏差平方和(4-5)应-y(Ti-y?2-1总的偏差平方和等于各列的偏差平方和之和。即(4-6)式中，m为正交表的列数若误差列由5个单列组成，则误差列的偏差平方和Se等于5个单列的偏差平方和之和，即：S=Si+S2+S3+S4+Se5；也可用Se=S总一S'来计算，其中：S为安排有因素或交互作用的各列的偏差平方和之和应引出的结论。与极差法相比，方差分析方法可以多引出一个结论：各列对试验指标的影响是否显著，在什么水平上显著。在数理统计上，这是一个很重要的问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 。显著性检验强调试验误差在分析每列对指标影响中所起的作用。如果某列对指标的影响不显著，那么，讨论试验指标随它的变化趋势是毫无意义的。因为在某列对指标的影响不显著时，即使从表中的数据可以看出该列水平变化时，对应的试验指标的数值也在以某种“规律”发生变化，但那很可能是由于实验误差所致，将它作为客观规律是不可靠的。有了各列的显著性检验之后，最后应将影响不显著的交互作用列与原来的“误差列”合并起来，组成新的“误差列”，重新检验各列的显著性。方差分析方法应用举例例4-6为了提高猪发酵饲料的营养和猪爱吃的程度，选择了四个因素进行正交试验，其因素水平见表4-18。表4-18例4-6的因素水平表因素发酵温度/C发酵时间/h初始的PH值投曲量/%符号X1X2X3X4水平110127522024610330485450724试验指标(y)为成品的总酸度。要求写出应用正交试验 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 方法的全过程，用方差分析方法分析正交试验的结果。解：试验的目的：为改善猪发酵饲料的品质，寻找适宜的发酵条件。试验指标(y)：成品的总酸度因素水平表：见表4-18。理论和经验都不知道有应该考虑的交互作用。四个因素的水平数不完全相同，所以应选择混合水平正交表。因为3个因素是4水平，1个因素是2水平，所以选Li6(43X26)正交表，见表4-19(a)表头设计：见表4-19(a)表中数据的计算举例：(以第3列为例)13=y1+y6+yn+y16=6.36+5.39+8.03+16.54=36.323=y2+y5+y12+y15=7.43+8.66+12.45+9.80=38.343=ya+y8+y9+y14=10.36+19.53+12.08+10.77=52.743=y4+y7+y1o+y13=11.56+15.50+13.13+13.49=53.68k3=43ka=36.32/4=9.083/ka=38.34/4=9.593/ka=52.74/4=13.19IVa/ka=53.68/4=13.42极差D3=13.42-9.08=4.34181.0811.32表4-19(a)使用正交表L16(43X26)的正交试验数据表列号123456789总酸度/%试验号X2X3X4X1eeeeey11111(10)(12)(7)11111⑸6.3612222(10)(24)(6)11222(10)7.4313323(10)(48)(5)22112(10)10.3614414(10)(72)(4)22221⑸11.5621225(20)(12)(6)22121(10)8.6622116(20)(24)(7)22212⑸5.3923417(20)(48)(4)11122⑸15.5024328(20)(72)(5)11211(10)19.5331319(30)(12)(5)12222⑸12.08324210(30)(24)(4)12111(10)13.13331211(30)(48)(7)21221(10)8.031234221112112.45皆羽-y)l=i-l218.35(30)(7'2)⑹⑸134(50)1(12)4(4)212122(10)13.49144(50)2(24)3(5)211211⑸10.77154(50)3(48)2(6)122111⑸9.80164(50)4(72)1(7)121222(10)16.54表4-19(b)[上接表4-19(a)]列号名称1X12X23X34e5e6e7e8e9X4Ij35.7140.5936.3283.91-LSIJOSIIj49.0836.7238.3497.17y-1132IIIj45.6943.69二52.74M50.6060.0853.68kj4n44888888Ij/kj3.9310.159.0810.49IIj/kj12.279.189.5912.15IIIj/kj11.4210.9213.19Vj/kj12.6515.0213.42极差D3.72③5.84①4.34②1.66④偏差平方和S33.57③79.19①63.67②se=30.911.02④u218.35自由度f.1j333fe=51方差V11.1926.4021.22U=6.1811.02方差比F1.814.273.431.78F0。251.881.881.881.69F0.103.623.623.624.06F0.055.41F0.01显著性0*(0.25)2*(0.10)1*(0.25)1*(0.25)偏差平方和I飼曜T唱=4(9.08—11.32)2+4(9.59—11.32)2+4(13.19—11.32)2+4(13.42—11.32)=63.67自由度f3=4—1=3方差Se=S总一(S+S+S3+S9)=218.35—(33.57+79.19+63.67+11.02)=30.9fe=(16-1)—(3+3+3+1)=5618a=0.01,f1=3,f2=5)=12.06〉F3a=0.05,f1=3,f2=5)=5.41>F3a=0.10,f1=3,f2=5)=3.62>F3a=0.25,f1=3,f2=5)=1.88<F3查F分布数据值表得:F(F(F(F(所以，第3列对试验指标的影响在a=0.25水平上显著。其它列的计算结果见表4-19(b)用方差分析方法分析正交试验结果，应该引出如下几点结论:(1)关于显著性的结论1015段曲壺眄(驸图4-6指标随各因素的变化趋势发酵时间(X2)对指标的影响在a=0.10水平上显著；初始的PH值(X3)和投曲量(X4)在a=0.25水平上显著；发酵温度(Xi)在a=0.25水平上仍不显著。⑵试验指标随各因素的变化趋势：见图4-6图4-6是用表4-18及表4-19(b)中的Ij/kj,nj/kj,Mj/kj,IVj/kj值来标绘的。(3)适宜的操作条件在确定适宜操作条件时，对于F检验中a=0.25不显著的因素，如本例中的因素X1,一方面因为图4-6(a)所示的“规律”不可靠，不能作为确定X1适宜水平的依据。另一方面，F检验不显著，Fj太小，可能是因为Ve太大，误差太大；也可能是因为V太小，该因素对指标影响太小。所以，对于F检验不显著的因素,适宜的水平可以是任意的。如本例，可认为X1=(20~50)C即可，不必非50C不可。所以在本例中为提高总酸度，适宜的操作条件为：Xi=(20~50)C,X2=72h,X3=4,X4=10%(4)对所得结论及进一步研究方向的讨论。由图4-6(d)可见，投曲量X4这个水平为试验范围的边上(最大值或最小值)所以X4增大，成品的总酸度也增大的结论尚需作进一步的研究。应研究投曲量大于10%寸试验指标随投曲量的变化规律。从图4-6(c)可见，初始PH值等于5时的总酸度与初始PH值等于4时的总酸度差不多。但与令PH=4相比较，令PH=5比较容易实现。所以进一步研究的方向之一，是研究令PH=5的好处和问题。从图4-6(b)可见，发酵时间愈长，成品的总酸度愈大，所以进一步研究的方向之一，是研究为提高总酸度而增长发酵时间的优缺点。例4-7为了提高某种产品的产量，寻求较好的工艺条件。考虑三个因素：反应温度、反应压力和溶液浓度。它们都取三个水平［见表4-20(a)］。表4-20(a)例4-7的因素水平表因素温度/C压力/kPa浓度/%符号(A)(B)(C)水平160200.5265251.0370302.0为考察3个因素间所有的两因素交互作用的影响，选正交表L27(313)，依该表的表头设计表得到的表头设计如表4-20(b)所示。表4-20(b)例4-7正交表表头设计列号12345]6789|10111213因素符号AB(AXB)1(AXB)2C(AXC)1(AXC)2(BXC)1e1e2(BXC)2e3e4可见，3水平两因素的交互作用占两列试验结果见表4-20(c)试验结果的方差分析计算见表4-20(d)17V=100.64y=3.72742?17]27H昭心壬卫七钿电=2^谕cc总的偏差平万和^161.02②两个三水平因素的交互作用占两列，它的S、f、V如何计算?以交互作用BXC为例。BXC占第8和第11列。偏差平方和Sbxc=S+S1=0.09187465+0.08907423=0.18095表4-20(c)正交表L27(313)的试验设计计算表列号12345678910111213产量/T因素符号AB(AXB)1(AXB)2C(AX6(AXC)2(BXC>1e1e2(BxC)2es8y试验号11111111111.3021111222224.6331111333337.2341222111230.5051222222313.6761222333126.2371333111321.3781333222134.7391333333217.07102123123110.47112123231223.47122123312336.13132231123230.33142231231313.40152231312125.80162312123320.63172312231133.97182312312216.50193132132110.03203132213223.40213132321336.80223213132230.57233213213313.97243213321126.83253321132321.0726332121313二3.97273321321216.57表4-20(d)方差分析计算表列号12345678119,10,12,13因素符号AB(AXB)1(AXB)2C(AXC)1(AxC)2(BXC)1(BxC)2e项目Ij36.7333.4635.6334.306.2732.9434.2133.3332.98IIj310.7031.3032.0831.7335.2134.6633.1333.0433.43IIIj33.2135.8832.9334.6159.1633.0433.3034.2734.23kj999999999j/kj4.0811113.7177783.9588893.8111110.69666673.6600003.8011113.7033333.664444IIj/3.4113.4773.56443.52553.91223.85113.68113.67113.7144kj11177844562211111144IIIj3.6903.9863.65883.84556.57333.67113.70003.80773.80331/kj000667895633110078330.76350.55530.20710.07490.09180.08902.0389411.166608201843155.8695406406974657423SeSSAxb=S3+S4=1.3189Saxc=Ss+Sf=0.282081SBXC=S3+Sh=0.18095=0.3445fj22fAXB=f3+f4=2+2=42fAXC=f6+f7=2+2=4fBXC=f8+f11=2+2=4fe=8VAxB=VAXC=Vj1.01950.58330SjWZk®=0.329726077.935=0.070520Vbxc=^Atc//AC=0.0453238Ve=0.043062Fj23.713.5Faxb=7.661.81X103Faxc=1.64Fbxc=1.051.0F0.018.658.657.018.657.01Fo.053.84F0.102.81Fo.251.661.66显著性4*4*(AXB):4*4*(AXC):0*(0.25)(BXC):0*(0.25)Fj/20.211.66.53F/°.°i6.236.234.77（显著性厂4*4*4*