正交试验结果的极差分析法

正交试验结果的极差分析法 正交试验方法能得到科技工作者的重视，在实践中得到广泛的应用，原因之一是不仅试验的次数减少，而且用相应的方法对试验结果进行分析可以引出许多有价值的结论。因此，在正交试验中，如果不对试验结果进行认真的分析，并明确地引出应该引出的结论，那就失去用正交试验法的意义和价值。 3下面以L(2)为例讨论正交试验结果的极差分析法。 4 3表4-13 L(2 )正交试验计算表 4 列 号 1 2 3 试验指标y i 试验号 1 1 1 1 y 1 2 1 2 2 y 2 3 2 1 2 y 3 n=4 2 2 1 y 4 I I = y+y I = y+y I = y+y j1j2213314 II II = y+y II = y+y II = y+y j134224323 k k =2 k =2 k =2 j123 I/k I/k I/k I/k jj112233 II/k II/k II/k II/k jj112233 极差(D) max{}-min{} max{}-min{} max{}-min{} j 在表4-13中: I——第j列“1”水平所对应的试验指标的数值之和。 j ?——第j列“2”水平所对应的试验指标的数值之和。 j (第j列有“3”，“4”水平时) k——第j列同一水平出现的次数。等于试验的次数(n)除以第j列的水平数。 j ——第j列“1”水平所对应的试验指标的平均值。 ——第j列“2”水平所对应的试验指标的平均值。 D——第j列的极差。等于第j列各水平对应的试验指标平均值中的最大值减最小值，即 j 用极差法分析正交试验结果应引出以下几个结论: ?在试验范围内，各列对试验指标的影响从大到小的排队。 某列的极差最大，表示该列的数值在试验范围内变化时，使试验指标数值的变化最大。 所以各列对试验指标的影响从大到小的排队，就是各列极差D的数值从大到小的排队。 试验指标随各因素的变化趋势。 ? ?使试验指标最好的适宜的操作条件(适宜的因素水平搭配)。 ?对所得结论和进一步研究方向的讨论。 例4-5 要求对例4-4的试验问题，写出应用正交试验设计方法的全过程，用极差法分析正 交实验的结果。 解:试验目的:提高磺化反应的乙酰胺苯的收率。 试验指标:乙酰胺苯的收率 表4-14 因素水平表 因素 反应温度/? 反应时间/h 硫酸浓度/% 操作方法 符号 A B C D 水平 1 A=50 B=1 C=17 D=搅拌 11112 =70 =2 =27 =不搅拌 ABCD2222应考虑的交互作用:A×B，A×C 7选择的正交表:L(2)，表头设计及计算结果均见表4-15。 8 应该引出的四个结论如下: 1( 各列对试验指标影响大小的排队问题。 因为极差 D=D=4.75最大，依次是D=2.75，D=D=2.25，D=1.75，最小是D=0.75。 3412765 所以，在本试验范围内，各因素和交互作用对试验指标的影响是交互作用A×B及因素C最 大，其次是因素A，再次要的是因素B、D，交互作用A×C的影响比“其他”因素的影响还小， 可以认为交互作用A×C不必考虑，让A×C参与排队毫无意义。 由此可见，正交试验法的优点是，即使无法直接利用表中几个试验来考察某列单独变化的 影响，它仍然能够求出每列或每个因素对指标影响的大小，求出对应于每个水平的指标平均值。 7表4-15 L(2)正交表应用计算表 8 列号 1 2 3 4 5 6 7 反应温度 因素 反应时间/h 硫酸浓度/% 收率y/% i /? 符号 A B A×B C A×C 其他 D 试验号 水平 1 1(50) 1(1) 1 1(17) 1 1 1(搅拌) y=65 12 1 1 1 2(27) 2 2 2(不搅拌) y=74 23 1 2(2) 2 1 1 2 2 y=71 34 1 2 2 2 2 1 1 y=73 45 2(70) 1 2 1 2 1 2 y=70 56 2 1 2 2 1 2 1 y=73 67 2 2 1 1 2 2 1 Y=62 78 2 2 1 2 1 1 2 y=67 8I 283 282 268 268 276 275 273 j II 272 273 287 287 279 282 282 j k 4 4 4 4 4 4 4 j I/k 70.75 70.50 67.00 67.00 69.00 68.75 68.25 jj II/k 68.00 68.25 71.75 71.75 69.75 70.50 70.50 jj 极差D 2.75 2.25 4.75 4.75 0.75 1.75 2.25 j 表中数据的计算举例，以第3列为例: j=3 I=I=y+y+y+y=65+74+62+67=268 j31278 II=II=y+y+y+y=71+73+70+73=287 j33456 K=k=8/2=4 j3 I/k=I/k=268/4=67.00 jj33 II/k=II/k=287/4=71.75 jj33 极差D=D=71.75-67.00=4.75 j3 2( 试验指标随各因素的变化趋势。 由表4-15中的第1(A)列:I/k=70.75， 1水平A=50? 111 II/k=68.00， 2水平A=70? 112 可见，反应温度(A)升高，收率下降。 同样可引出结论: 反应时间(B)加长，收率下降。 硫酸浓度(C)增大，收率增大。 操作方法()由搅拌改为不搅拌，收率增大。 D 3( 适宜的操作条件。 首先应搞清，所讨论问题的试验指标的数值是大好还是小好。很明显，本题的试验指标收 率是愈大愈好。 在确定适宜操作条件时，应优先考虑对试验指标影响大的试验因素和交互作用。也就是说 必须按对试验指标的影响从大到小的顺序，来确定适宜的操作条件。 ? 对于交互作用A×B: I/k=67.00时交互作用(A×B)列??A 、B的水平号码相同，即AB或AB。 3311122II/k=71.75时交互作用(A×B)列对应的因素A 、B的水平号码不同即AB或AB。 3321221可见，为提高收率，首先应让(A×B)=(A×B)，即应让A、B的水平号码不同。 2 ? 对于C因素，宜取2水平， ? 对于A因素，宜取1水平， ? 对于B因素，应让A、B的水平号码不同，故应取2水平， ? 对于D因素，宜取2水平。 所以，为提高收率，在本试验范围内，适宜的操作条件为: 反应温度A=A=50?;反应时间B=B=2h; 12 ==27%;操作方法==不搅拌。 硫酸浓度CCDD22 4( 对所得结论和进一步研究方向的讨论 这个问题涉及到正交试验设计方法的有关知识，也涉及到所研究的试验问题的理论知识和实践知识。 ? 由搅拌改为不搅拌，收率增大。这个结论很可能是错误的。因为加搅拌的唯一目的是使反应器内的浓度场和温度场更加均匀一致，这会使收率下降，是没有道理的。另外，因素的极D差=2.25，与“其他”因素的极差=1.75，相差不多，和的相对偏差 DDDD7676 可见，操作方法对指标的影响不算太大，可认为不搅拌有利于提高收率的结论是实验误差造成的。 也可这样认为，理论上搅拌有利于提高收率，但在本试验的范围内，实际上提高的不多。既然搅拌带来的好处不大，当然就不要搅拌了。 ? 在适宜操作条件中，取A=A=50?，B=B=2h, 与B=1h相比，时间长了1倍，这显然是121 不好的。在维持A×B=(A×B)的前提下，改取A=A=70?，B=B=1h，行吗, 221 由表4-15中的#3，#4号试验知，A=A，B=B时，收率y的平均值=(71+73)/2=72%。 12 由表4-15中的#5，#6号试验知，A=A，B=B时，收率y的平均值=(70+73)/2=71.5%。 21 可见，若改A=A=70?，B=B=1h, 收率只下降0.5%，而反应时间却缩短了一半。 21 所以，遇到这类情况，必要时，可重复进行试验，以便将适宜操作条件改为: A(70?)，B(1h)，C(27%), D(不搅拌)。 2122 ? 从二元图看交互作用A×B和A×C 交互作用A×B、A×C的二元表分别见表4-16和表4-17。 A×C的二元图分别见图4-4和图4-5。 交互作用A×B、 表4-16 交互作用A×B的二元表 收率 y/% 因素、水平 B=1h B=2h 12 A=50? (y+y)/2=(65+74)/2=69.5 (y+y)/2=(71+73)/2=72.0 11234 A=70? (y+y)/2=(70+73)/2=71.5 (y+y)/2=(62+67)/2=64.5 26678注: y~y的数据从表4-15中读取(下同)。 18 表4-17 交互作用A×C的二元表 收率 y/% 因素、水平 C=17% C=27% 12 =50? (y+y)/2=(65+71)/2=68.0 (y+y)/2=(74+73)/2=73.5 A11324 A=70? (y+y)/2=(70+62)/2=66.0 (y+y)/2=(73+67)/2=70.0 25768 由图4-4可见，与有效明显的交互作用，而因素与是否有交互作用(见图4-5)需ABAC要进一步作误差估算。 由表4-15得C=C=17%时直线的斜率 1 C=C=27%，A=70?时，y=70.0% 22实验 y对y的相对偏差 实验计算 -2e(y，y)= =-2.1×10， 实验 计算r 在常见的测试误差范围[?5%左右]之内，故可认为两直线不平行是由于数据误差所致, 交互作用A×C并不存在。