第7章正交试验设计的极差分析.docx

1、第7章正交试验设计的极差分析第7章正交试验设计的极差分析正交试验设计和分析方法大致分为二种：一种是极差分析法（又 称直观分析法）,另一种是方差分析法（又称统计分析法）。本章介绍极 差分析法，它简单易懂，实用性强，在工农业生产中广泛应用。7.1单指标正交试验设计及其极差分析极差分析法简称 R法。它包括计算和判断两个步骤，其内容如 图7-1所示。图7-1 R法示意图图中，Km为第j列因素m水平所对应的试验指标和，Kjm为Km的平均 值。由Km的大小可以判断j因素的优水平和各因素的水平组合，即 最优组合。R为第j列因素的极差，即第j列因素各水平下平均指标 值的最大值与最小值之差：R二max（心，耳,

2、，爲）-min（心,兀,爲）R反映了第j列因素的水平变动时，试验指标的变动幅度。 R越大，说明该因素对试验指标的影响越大，因此也就越重要。于是依据Rj的大小，就可以判断因素的主次。极差分析法的计算与判断，可直接在试验结果分析表上进行，现以例6 - 2来说明单指标正交试验结果的极差分析方法。一、确定因素的优水平和最优水平组合例6-2为提高山楂原料的利用率，某研究组研究了酶法液化工艺制 造山楂精汁。拟通过正交试验寻找酶法液化工艺的最佳工艺条件。在例6 -2中，不考虑因素间的交互作用（因例6 - 2是四因素三水平试验，故选用L9（34）正交表），表头设计如表6 - 5所示，试验方 案则示于表6 -

3、6中。试验结果的极差分析过程，如表7 - 1所示.表6-4 因素水平表加水量加酶量酶解温度酶解时间水平因素(ml/100g)(ml/100g)(9)(h)ABCD1101201.52504352.53907503.5表6-6 试验方案及结果试验号因 素试验结果液化率（%）ABCD11(10)1(1)1(20)1(1.5)0.002122(35)2(2.5)17.03133(50)3(3.5)24.042(50)12312.05223147.06231228.073(90)1321.008321318.09332142.0试验指标为液化率，用y表示，列于表6 - 6和表7 - 1的最后一一表7-

4、1试验方案及结果分析试验号因素试验结果液化率(%)ABCD11(10)1(1)1(20)1(1.5)0.002122(35)2(2.5)17.03133(50)3(3.5)24.042(50)12312.05223147.06231228.073(90)1321.008321318.09332142.0K41.013.046.089.0迟=189.0K287.082.071.046.0K361.094.072.054.0K；13.74.315.329.7K729.027.323.715.320.331.324.018.0优水平AC3DRj15.327.08.714.4主次顺序B A D C计算

5、示例：因素A的第1水平Ai所对应的试验指标之和及其平均值分别为:1Ki=yi+y2+y3=0+17+24=41, K ai = KAi=13.73同理，对因素A的第2水平A和第3水平A,有1 KA2=y4+y5+y6=12+47+28=87 Ka2 = - Ka2=2931KA3=y7+y8+y9=1 + 18+42=61 , KA3 Ka3=20.33由表7 - 1或表6 - 6可以看出，考察因素 A进行的三组试验中(A,A2,A3), B、C D各水平都只出现了一次，且由于 B C、D间无 交互作用，所以B、C D因素的各水平的不同组合对试验指标无影响， 因此，对A、A和A来说，三组试验的

6、试验条件是完全一样的。假如 因素A对试验指标无影响，那么 石,兀,心 应该相等，但由上面的 计标可知，实际上并不相等，显然，这是由于因素A的水 平变化引起的，因此，心,瓦2,瓦；的大小反映了 A、A和A对试验指 标影响的大小。由于液化率y越大越好，而 兀 k3 心，所以可判 断A为因素A的优水平。同理，可判断因素 B C、D的优水平分别为 B、C3、D。所以， 优水平组合为 ABGD，即最优工艺条件为加水量 A=50ml/100g、加 酶量B=7ml/100g、酶解温度G=50。C和酶解时间D=1.5小时。二、确定因素主次顺序极差R按定义计算，如Ra = K A2 一 KA厂 29.0 一 1

7、3.7 二 15.3 ,Rb 尺-爲31.27.0同理可求出RC和FD.计算结果列于表7-1中。比较R值可知 FBFAFDFC,所以试验因素对试验指标的影响的主次顺序为 BADC即 加酶量影响最大，其次是加水量和酶解时间，而酶解温度的影响最小。三、 绘制因素与指标趋势图为了更直观地反映因素对试验指标的影响规律和趋势 ，用因素的 水平作横坐标，试验指标的平均值（Kj）作纵坐标，画出因素与指标的 关系图（即趋势图）,如图7-2所示.（p137）趋势图可为进一步试验时选择因素水平指明方向.如对因素A,由 图7-2可见,A2水平时,指标最咼,但若能在A附近再取一些水平（如 40、60）作进一步试验，则

8、有可能取得更高的指标；对D因素,若能取一 些比D更小的水平（如1.0和0.5）作进一步试验，也有可能得到更好 的结果.以上三个步骤即为极差分析的基本程序与方法.四、 说明与讨论1. 计算结果的检验：每一列的K之和应等于全部试验结果（即指标值）m n之和，即Kj yj , m为水平数,n为试验总实施次数.t jm2. 因素的最优水平组合，在实际处理中是灵活的，即对于主要因素，一定要选最优水平；而对次要因素，则应权衡利弊，综合考虑其它条件进 行水平选取，从而得到最符合实际生产的最优或较优生产工艺条件 .3. 例6-2的最优工艺条件ABGD并不在实施的9个试验之中.这表明 优化结果不仅反映了已做的试

9、验信息，而且反映了全面试验信息.因 此,正交试验设计的部分实施方案反映了全面试验信息 .4. 例6-2得出的最优工艺条件，只有在试验所考察的范围内才有意义， 超出这个范围 ,情况就可能发生变化。另外 , 只能说是“较优工艺条 件”,而不能说是“最优工艺条件 ”. 最好能根据趋势图做进一步试验 找出最靠近最优的工艺条件 .5. 对已确定的最优工艺条件（如例6-2的ABGD）进行重复试验，验证 其试验指标是否最优 .7.2 多指标正交试验设计及其极差分析在实际生产和科研试验中 , 所要考察的指标往往不止一个 , 这一类 的试验设计叫做多指标试验设计 . 在多指标试验设计中 , 各指标之间 可能存在

10、一定的矛盾 , 如何兼顾各个指标 , 找出使每个试验都尽可能 好的试验条件呢 ?换言之，应如何分析多指标试验设计的结果呢 ?常用 的有两种方法 : 综合平衡法和综合评分法 . 下面举例说明综合平衡法 的分析方法 .这种方法在试验方案安排和各指标计算分析方法上 , 与单指标试 验完全一样 . 其步骤是先分别找出各个指标最优或较优的生产条件 , 然后将这些生产条件综合平衡 , 找出兼顾每个指标都尽可能好的生产 条件.例 7-1 在油炸方便面的生产中，主要原料质量和主要工艺参数对产 品的质量有影响。今欲通过正交试验确定最佳生产条件。一 . 试验方案设计1. 确定试验指标评价方便面质量好坏的主要指标是

11、 : 脂肪含量 （越低越好 ），水分含量(越高越好)和复水时间(越短越好)2 .挑因素，选水平，列出因素水平表根据专业知识和实际经验，确定试验因素和水平，如表7-2所示表7-2 因素水平表水平因素湿面筋值(%)A改良剂用量(%)B油炸时间(s)C油炸温度(9)D1280.05701502320.075751553360.10801603 .选正交表，设计表头，编制试验方案本试验是四因素三水平试验，不考虑因素间的交互作用，因此， 可应选L9(34)安排试验，表头设计和试验方案见表 7-3 (p140)。按上述方案实施后，将每一项试验指标都记录下来，见表 7-3注：对极差分析可以这样选正交表，但对

12、方差分析应留有空列，以便 估计试验误差.表7-3试验方案及结果分析因素试验结果试验号ABCD脂肪(%)水分(%)复水时间(s)11(28)1(0.05)3(80)2(155)24.82.13.5212(0.075)1(70)1(150)22.53.83.7313(0.10)2(75)3(160)23.62.03.042(32)12123.82.83.05223322.41.72.26231219.32.72.873(36)11318.42.53.08322219.02.02.79333120.72.33.6KK270.965.567.063.060.266.467.063.1送=194.5旨

13、方 J 二八 -_ 二K358.163.667.964.4K23.622.320.122.3含量21.821.322.121.019.421.222.621.5R4.21.12.51.3K7.97.49.08.9迟=21.9K27.27.56.86.8水 分K36.86.96.16.2K2.632.473.002.97含量2.402.502.272.272.272.302.032.07R0.360.200.970.90K10.29.59.510.3Z=27.5K28.08.68.79.0复K39.39.49.38.2水K3.403.173.173.43时 间2.672.872.903.00K3

14、3.103.133.102.73R0.730.300.270.70二.试验结果分析1. 计算每列各水平下每种试验指标的数据和（ K, K2, K3）,及 其平均值（心,心,心），并计算极差R,填入表7-3中。2. 画出因素与各种指标的趋势图，如图 7-3所示（p140）3. 按极差大小列出各指标下各因素主次顺序:各因素主次顺序表主-次脂肪含量（%） 水分含量（咧 复水时间（s）4. 初选最优工艺条件根据各指标下的平均数据和K1,k；,K；，初步确定各因素的最优水平组合为:5. 综合平衡确定最优工艺条件（难点）！ 由于三个指标单独分析出来的最优条件并不一致， 所以必须根据 因素对三个指标影响的主

15、次顺序，综合考虑，确定出最优条件。首先，把水平选取上没有矛盾的因素的水平定下来， 即如果对三 个指标影响都重要的某一因素，都是取某一水平时最好，则该因素就 是选这一水平。在本试验中无这样的因素，因此我们只能逐个考察每 一因素。对因素A:从主次顺序来看，对脂肪含量和复水时间的影响都排在第一位为主要因素，而对水分含量的影响则排在第三位，属次要因 素，因此，应以主要因素为主选因素的水平。从初选的最优水平组合 中可以看出，对脂肪含量选 A3为好，而对复水时间，则选 A为好。 因为二者不一致，所以还须根据试验结果分析确定选 A还是A从表 7-3可知，当取 A时，复水时间比取 A时缩短16.1%(有利)，

16、即 (2.67-3.10) -2.67 X 100%=-16.1% 而脂肪含量只比取 A时增加11.0% (不利)，即(21.8-19.4) - 21.8 X 100%=11.0% 且从水分含量指标来看，取A也比取A3时更好，因此，应选取 A水平。注： 当取A时，脂肪含量比取A时降低12.4%(有利)，即 (19.4-21.8)/19.4 X 100%=- 12.4%,复水时间比取 A时增加 13.9(不 利)，即(3.10-2.67 ) /3.10 X 100%=13.9%综合平衡a不利有利A11.0%16.1%A13.9%12.4%对“有利”部分,A2Aj；对“不利”部分,A2 “不利”；

17、选D3时，“不利” “有 利”.并且D （有利）Da（有利之和绝对值）,D i （不利之和） D 3（不利 绝对值）.因此,从定量分析来看，D应取D,而不是取D3.那么,究竟如 何决定D的水平呢？最后,应该再进行 ABCDi和ABCC3两次试验，由 试验结果决定D好还是D好!实践是检验真理的唯一标准！7.3混合型正交表的试验设计极差分析前面讨论的都是水平数 Ln（mk）相同的正交试验设计.但在实际工 作中,有些试验受到设备、原材料和生产条件等限制某些因素的水平 选择受到制约，或者在有些试验中，要重点考察某个（或某些）因素需 要多取几个水平，这时就会遇到水平数不同的正交试验设计 .在这种 情况下

18、,通常有三种解决方法：一是直接选用合适的混合型正交表；二 是采用拟水平法；三是采用拟因素法 .我们现在只讨论第一种方法， 即使用混合型正交表Ln（mik1 m22）进行正交试验设计.例7-2 某油炸膨化食品的体积与油温、物料含水量及油炸时间有关， 为确保产品质量，提出工艺要求。现通过正交试验设计寻求理想的工 艺条件。一. 试验方案设计1. 确定试验指标本试验的指标为油炸膨化食品的体积，体积越大越好.2. 挑因素、选水平、制定因素水平表根据专业知识，制定因素水平表如7-4所示,因素A取4个水平, 因素B和C各取2个水平,所以属于水平数不相等的正交试验设计.表7-4 因素水平表水平因素油炸温度（9

19、）A物料含水量（%）B油炸时间（s）C12102.03022204.040323042403. 选正交表、设计表头、编制试验方案本试验宜选用L8（41 X 24）正交表安排试验，表头设计时，把A因素 放在第一列，其余两个因素可随意安排在四个二水平列中 ，比如依次排在第二、三列中，把所安排因素的各列的水平数字后标上相应因素 的具体水平值，即得出试验方案，如表7-5所示.按表7-5试验方案实施后，所得试验结果列于表7-5中的最后一 列.表7-5试验方案及结果分析试验号油温A含水量B时间C体积Xi(cm3/100g)1234511(210)1(2.0)1(30)11210.0212(4.0)2(40

20、)22208.032(220)1122215.0422211230.053(230)1212251.0632121247.074(240)1221238.0842112230.0K1418.0914.0902.01=1829.0&445.0915.0927.0K3498.0K4468.0209.0228.5225.5心222.5228.75231.75249.0石234.0R40.00.256.251R25.460.3558.875二.试验结果分析1. 计算各列各水平下的K、K及R由于各列的水平数不完全相同，所以K和K的计算略有差异.第1列：由于有四个水平数，所以要计算四个K与Q每个K由二 个

21、数据相加得到，因此K=K/2.例如：KA =210.0 208.0 =418.0, K凡=KA/2 = 418.0 2 = 209.0R =249.0 -209.0 =40.0第2、3列:由于只有两个水平，所以只要计算两个K与k_,每个K 由四个数据相加得到，因此K=K/4.例如：KBi 210.0 215.0 251.0 238.914.0Kbi = Kb/4 = 914.0/4 = 228.5R =228.75-228.5 =0.25按上述方法计算出各列各水平下的 K、K以及R值，列于表7-5中.2. 计算R的折算值R （极差R的折算）当因素的水平数相同时，因素的主次顺序完全由 R决定.但

22、当因 素的水平数不同时，直接比较R是不行的.这是因为，若两个因素对试 验指标有影响，一般来说,水平数多的因素极差可能大一些.因此，要 用一个系数把极差R折算后才能作比较.极差的折算公式如下：R = dR.r式中 R -折算后的极差；R-因素的极差；r-该因素每个水平试验的重复数，r=-；md-折算系数，与因素的水平数有关，其值见表 7-6。表7-6 折算系数表 m r dR -水平数m2345678910折算系数d0.71 0.52 0.45 0.40 0.37 0.35 0.34 0.32 0.31本例中，R的折算如下：Ra =0.45 40 2 = 25.46RB =0.71 0.25 4

23、 =0.355Rc =0.71 6.25 一4=8.875计算结果列于表7 - 5中.3. 根据R大小确定因素的主次顺序主 - 次A C B即油炸温度对实验指标的影响最大，其次是油炸时间，而物料含 水量的影响最小。4. 画出因素指标趋势图，如图 7-4 所示（ p146）5. 选各因素的最优水平及最优水平组合比较各因素各水平下的 K 值（本例中 K 越大越好），并参考因素指 标趋势图，得出最优水平组合为A3B2C2或ABG,即油炸温度230摄氏 度，油炸时间 40秒，物料含水量对试验指标影响很小，故取 2%或 4% 都可以，视具体情况而定。由表7-5可见，若最优水平组合 ABC2,则该试验即表

24、中的第 5 号试验，实验指标值即膨化体积为251.0 cm 3/100g ,为表中所列最大 值；若最优水平组合为ABG,则需再实施一次该水平组合下的试验， 作为验证。7.4 考察交互作用的正交试验设计及极差分析 一、交互作用的概念前面介绍的正交试验设计与试验结果的分析方法，都是指因素间 没有（或不考虑）交互作用的情况，实际上，在许多试验中，不仅因 素对指标有影响，而且因素之间还会联合搭配起来对指标产生影响。 所以，因素对试验产生的总效果， 是由每一个因素对试验的单独作用 再加上各个因素之间的搭配作用决定的。 这种因素间的联合搭配对试验指标产生的影响作用，称为交互作用。例如，我们要考虑化学反应的温度（A）与时间（B）对产品收率 的影响，温度和时间都取二个水平，即 A A2和B B2。在各AB组合 条件的平均产品收率，可能有如下三种情况：（1） 不论B因素