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

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

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

3、如表7-1所示.表6-4因素水平表X加水量(ml/10加酶量(ml/ 1酶解温度酶解时间水車、素0 g)00g)(C)(h)XABcD11 01201 .52504352.539 07so3. 5表6-6 试验方案及结果试验号因素试验结果液化率(%)ABcD11( 1 0)1(1)1 (20)1(1. 5)0.00212(4)2(35)2(2. 5)17. 0313 (7)3 (50)3 (3. 5)24 . 042(50)12312.0522314 7.06231228.073 (90)1321.008321318.09332142. 0试验指标为液化率，用y,表示，列于表6-6和表7-1的

4、最后一 列。表7-1 试验方案及结果分析试验号因素试验结果 液化率（）ABcD11(10)1 (1)1 ( 20)1(1.5)0.00212(4)2(35)2 (2.5)17. 0313(7)3(50)3(3. 5)24.042(5 0)1231 2.05223147. 06231228.073 (90)1321. 008321318.09332142.0K4 1.01 3.04 6. 089. 0为=189. 0K:87. 08 2. 071. 04 6.0K361.094. 072 . 05 4.0K】374.315. 32 9. 7K；29.02 7.323.715.320.331. 3

5、24. 018. 0优水平A2B,C3D：Rj15. 327. 08.71 4.4主次顺序B AD C计算示例:因素A的第1水平A】所对应的试验指标之和及其平均值分别为: 1KAi=y i+y+y 3=0+ 1 7 + 2 4二 4 1, =矗讦 1 3. 7同理，对因素A的第2水平A 2和第3水平A3,有Ka2= y 4+ y s+y6= 1 2 +47+28=87, 7=1ka2=2 9K.=y7+ys+y9= 1 +18+42 = 61, F = KA3=20. 3由表71或表6-6可以看出，考察因素A进行的三组试验中(AbA2, A3),B. C、D各水平都只出现了一次,且由于B、C、

6、D间无 交互作用，所以B、C、D因素的各水平的不同组合对试验指标无影响, 因此，对入、A：和乩来说，三组试验的试验条件是完全一样的。假如 因素A对试验指标无影响，那么石，石，石 应该相等，但由上而的 计标可知，石，石，石实际上并不相等，显然，这是由于因素A的水平 变化引起的，因此，石，兀，石的大小反映了 A、A?和码对试验指 标影响的大小。由于液化率y越大越好，而兀兀兀，所以可判 断A：为因素A的优水平。同理，可判断因素B、C、D的优水平分别为Bs、C3、D1O所以， 优水平组合为A泯CD,即最优工艺条件为加水量A2= 5 0ml/10 0 g、 加酶量B3=7ml/1 OOg、酶解温度C,=

7、5 0 - C和酶解时间21. 5小时。二、确定因素主次顺序极差出按定义计算，如Ra =兀一石= 29.0 13.7 = 15.3,念=石石= 31.3-4.3 = 27.0同理可求出Rc和Rd.计算结果列于表7- 1中。比较Rj值可知 RbRaRdRc,所以试验因素对试验指标的影响的主次顺序为BAD Co 即加酶量影响最大，其次是加水量和酶解时间，而酶解温度的影响最三、 绘制因素与指标趋势图为了更直观地反映因素对试验指标的影响规律和趋势，用因素的 水平作横坐标，试验指标的平均值（心）作纵坐标，画出因素与指标的 关系图（即趋势图），如图72所示.（pl 3 7）趋势图可为进一步试验时选择因素水

8、平指明方向.如对因素扎由 图7- 2可见，A：水平时，指标最高，但若能在A：附近再取一些水平（如 4 0、6 0 ）作进一步试验，则有可能取得更高的指标;对D因素,若能取 一些比D更小的水平（如1.0和0. 5）作进一步试验，也有可能得到 更好的结果.以上三个步骤即为极差分析的基木程序与方法.四、 说明与讨论1.计算结果的检验：每一列的K,之和应等于全部试验结果（即指标值） 之和，即士K产士兀，m为水平数,n为试验总实施次数.J-1 /】2.因素的最优水平组合，在实际处理中是灵活的，即对于主要因素， 一定要选最优水平;而对次要因素，则应权衡利弊，综合考虑其它条件 进行水平选取，从而得到最符合实

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

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

11、方案设计1确定试验指标评价方便而质量好坏的主要指标是：脂肪含量(越低越好)，水分含量(越高越好)和复水时间(越短越好)。2.挑因素，选水平，列岀因素水平表根据专业知识和实际经验,确定试验因素和水平，如表7- 2所示。表7-2 因素水平表湿面筋值(%)A改良剂用量 (%)B油炸时间 (s)C油炸温度(。0D1280. 05701 502320. 0757515 53360. 1 08 01603.选正交表，设计表头,编制试验方案本试验是四因素三水平试验，不考虑因素间的交互作用，因此, 可应选L9(3)安排试验，表头设计和试验方案见表7-3 ( p 1 4 0 ) o按上述方案实施后，将每一项试验

12、指标都记录下来，见表7-3o 注：对极差分析可以这样选正交表,但对方差分析应留有空列，以便 估计试验误差.表7-3试验方案及结果分析试验号因素试验结果ABCD脂肪(%)水分(%)复水时间(s)11(28)1 (0.03(80)2 (155)24.82 . 13.5215)1(7 0)1(150)22.53.83. 7312(0. 075)2(75)3(1 6 0)2 3.62 . 03. 042(32)3 (0. 10)2123.82.83 . 0521332 2.41. 72. 26221219.2. 72.873 (36)31332. 53.08312218.42. 02.79323119

13、.02.33.6320. 7K：70.967.060.267.0=194.5K：6 5.563.066.463. 1&58. 163.667.964.4/J FI 肪K23.622.32 0. 12 2. 3含量瓦2 1.821.322.12 1. 01 9.421. 222.621. 5R4 . 21.12. 51. 3Ki7.97.49.08. 9工= 21.9K27.27.56.86.8水k36.86.96 . 16.2分K2. 6 32 . 473. 002.97含疋2. 402. 502. 2 72. 27量2.2 72. 302. 0 32.07R0. 3 60. 200. 9 7

14、0. 90K：10.29.59. 510.3工二27. 5k28.08.6& 79.0k39 . 39.49. 38. 2夂水K3. 4 03. 173. 1 73.4 3时疋2.672.872. 903.00间K3. 1 03. 133. 1 02. 7 3R0.7 30.300 . 270. 70二.试验结果分析1.计算每列各水平下每种试验指标的数据和（Ki，K:, KJ,及其 平均值（石，心,氐），并计算极差R,填入表7-3中。2.画出因素与各种指标的趋势图，如图7-3所示（pl40）。3.按极差大小列出各指标下各因素主次顺序：各因素主次顺序表试验指标主一次脂肪含量（%）A C D B水

15、分含量（）C D A B复水时间（S）A D B C4初选最优工艺条件根据各指标下的平均数据和KK 尺,初步确定各因素的最优水 平组合为：对脂肪含量（%）： A3B3C】 （脂肪含量越低越好）对水分含量（%） : A iB： C ! Di （水分含量越高越好）对复水时间（s）: A:B2C：D：5 （复水时间越短越好）5.综合平衡确定最优工艺条件（难点）！由于三个指标单独分析出来的最优条件并不一致，所以必须 根据因素对三个指标影响的主次顺序，综合考虑，确定出最优条件。首先，把水平选取上没有矛盾的因素的水平定下来，即如果对三 个指标影响都重要的某一因素，都是取某一水平时最好，则该因素就 是选这一

16、水平。在本试验中无这样的因素，因此我们只能逐个考察每 一因素。对因素A:从主次顺序来看，对脂肪含量和复水时间的影响都排在 第一位为主要因素，而对水分含量的影响则排在第三位,属次要因素, 因此,应以主要因素为主选因素的水平。从初选的最优水平组合中可 以看出，对脂肪含量选As为好，而对复水时间，则选A?为好。因为 二者不一致，所以还须根据试验结果分析确定选短还是乩。从表7-3 可知，当取航时,复水时间比取仏时缩短16. 1%（有利），即（2.6 7 -3. 10）4-2. 6 7 X1 0 0 %=- 1 6. 1%,而脂肪含量只比取仏时增 加 11.0%（不利），即 （21. 8 -19 4）+

17、2 1 . 8 X 1 00%=11.0%, 且从水分含量指标来看，取A 2也比取乩时更好，因此，应选取比水平。 注：当取时，脂肪含量比取乩时降低1 2. 4%（有利），即（19. 4 一2 18 ） /19.4X 1 0 0%二一 1 2. 4%,复水时间比取A?时增加1 3 . 9%（不利）,即（3. 10-2.67）/3. 1 0 X 100% 二 13. 9%。综合平衡A不利有利a211. 0%1 6 . 1%A313. 9%12.4%对“有利”部分,A2As；对“不利”部分,A2“不利”；选6时，“不利” “有 利”.并且D】（有利）D3（有利之和绝对值），Dx （不利之和） D,不

18、利绝对值）.因此从定量分析来看,D应取D而不是取D3.那么, 究竟如何决定D的水平呢?最后，应该再进行A：B：CD】和两 次试验，由试验结果决定D好还是6好!实践是检验真理的唯一标准!7. 3混合型正交表的试验设计极差分析前而讨论的都是水平数M）相同的正交试验设计.但在实际工 作中，有些试验受到设备、原材料和生产条件等限制.某些因素的水平 选择受到制约，或者在有些试验中，要重点考察某个（或某些）因素需 要多取几个水平，这时就会遇到水平数不同的正交试验设计.在这种 情况下，通常有三种解决方法:一是直接选用合适的混合型正交表; 二是采用拟水平法；三是采用拟因素法.我们现在只讨论第一种方法, 即使用

19、混合型正交表厶0伸X今）进行正交试验设计.例7-2 某油炸膨化食品的体积与油温、物料含水量及油炸时间有 关，为确保产品质量，提出工艺要求。现通过正交试验设计寻求理想 的工艺条件。一. 试验方案设计1.确定试验指标木试验的指标为油炸膨化食品的体积，体积越大越好.2.挑因素、选水平、制定因素水平表根据专业知识，制定因素水平表如7-4所示，因素A取4个水平, 因素B和C各取2个水平，所以属于水平数不相等的正交试验设计.表7-4 因素水平表因油炸温度（。C ）A物料含水量（）B油炸时间（s）C12102. 03022204.04032 3042403.选正交表、设计表头、编制试验方案本试验宜选用Ls（

20、4】X2J正交表安排试验，表头设计时，把A因素 放在第一列，其余两个因素可随意安排在四个二水平列中，比如依次 排在第二、三列中，把所安排因素的各列的水平数字后标上相应因素的具体水平值，即得岀试验方案，如表7-5所示.按表7-5试验方案实施后,所得试验结果列于表7- 5中的最后一列.表7-5试验方案及结果分析试验号油温A含水量B时间C体积&(c m3/100 g )1234511(210)1 (2. 0)1(30)11210.0212 (4. 0)2 (40)222 0 8. 032 (2 2 0)1122215.0422211230.053(230)1212251.0632122 47. 07

21、4(240)122123 8.0812112230.0K,418.0914.09 02.0工=1k24 45.09 15. 0927. 0829.0K,4 9 8-0K.4 6 8 0K2 09. 022 8. 5225.5K；222. 5228. 7523 1 . 75&24 9. 0&2 34. 0R40. 00.256. 2 5R25. 460.3558. 87 5二.试验结果分析1.计算各列各水平下的K、斤及R由于各列的水平数不完全相同，所以K和尺的计算略有差异.第1列：由于有四个水平数，所以要计算四个K与斤，每个K由 二个数据相加得到，因此斤=K/2.例如：K、=210.0 + 20

22、8.0 = 418.0, Ka, = K*/2 = 418.0一 2 = 209.0/? = 249.0 - 209.0 = 40.0第2、3列：由于只有两个水平，所以只要计算两个K与瓦每个K 由四个数据相加得到，因此K=K/4.例如：0 =210.0 + 215.0 + 251.0 + 23&0 = 914.0Kb. =Kb/4 = 914.0/4 = 228.5R = 228.75 - 22&5 = 0.25按上述方法计算出各列各水平下的K、&以及R值，列于表7-5 中.2.计算R的折算值R（极差R的折算）当因素的水平数相同时，因素的主次顺序完全由R决定但当因 素的水平数不同时，直接比较R

23、是不行的.这是因为，若两个因素对试 验指标有影响，一般来说，水平数多的因素极差可能大一些.因此，要 用一个系数把极差R折算后才能作比较.极差的折算公式如下：R = dRr式中 疋一折算后的极差；R-因素的极差；r该因素每个水平试验的重复数,r二丄;Hid-折算系数，与因素的水平数有关，其值见表7-60水平数m2 3 4 5 6 7 8 910折算系数d0. 7 1 0. 52 0. 45 0.40 0. 37 0.3 50. 34 0. 32 0. 31木例中，R的折算如下：Ra = 0.45 x 40 x = 25.46Rb =0.71x0.25x74 = 0.355 =0.71x6.25x = 8.875计算结果列于表7-5中.3.根据R大小确定因素的主次顺序主一次A C B即油炸温度对实验指标的影响最大，其次是油炸时间，而物料含 水量的影响最小。4.画岀因素指标趋势图，如图7-4所示(pl 4 6 )5.选各因素的最优水平及最优水平组合比较各因素各水平下的斤值(木例中斤越大越好)，并参考因素指 标趋势图，得出最优水平组合为ADG或AECz,即油炸温度230 摄氏度，油炸时间40秒，物料含水量对试验指标影