1、第章正交试验设计的极差分析第7章 正交试验设计的极差分析正交试验设计和分析方法大致分为二种:一种是极差分析法(又称直观分析法),另一种是方差分析法(又称统计分析法)。本章介绍极差分析法,它简单易懂,实用性强,在工农业生产中广泛应用。7.1 单指标正交试验设计及其极差分析极差分析法简称R法。它包括计算和判断两个步骤,其内容如图7-1所示。图7-1 R法示意图图中,Kjm为第j列因素m水平所对应的试验指标和, jm为Kjm的平均值。由Kjm的大小可以判断j因素的优水平和各因素的水平组合,即最优组合。Rj 为第j列因素的极差,即第j列因素各水平下平均指标值的最大值与最小值之差:Rj=max()-mi
2、n()Rj反映了第j列因素的水平变动时,试验指标的变动幅度。Rj越大,说明该因素对试验指标的影响越大,因此也就越重要。于是依据Rj的大小,就可以判断因素的主次。极差分析法的计算与判断,可直接在试验结果分析表上进行,现以例-来说明单指标正交试验结果的极差分析方法。一、 确定因素的优水平和最优水平组合例6-2 为提高山楂原料的利用率,某研究组研究了酶法液化工艺制造山楂精汁。拟通过正交试验寻找酶法液化工艺的最佳工艺条件。在例-中,不考虑因素间的交互作用(因例-是四因素三水平试验,故选用L9(34)正交表),表头设计如表-所示,试验方案则示于表-中。试验结果的极差分析过程,如表-所示.表6-4 因素水
3、平表水平 因素加水量(ml/100g)A加酶量(ml/100g)B酶解温度(C)C酶解时间(h)D1231050901472035501.52.53.5表6-6 试验方案及结果试验号因 素试验结果液化率(%)ABCD1234567891(10)112(50)223(90)331(1)2(4)3(7)1231231(20)2(35)3(50)2313121(1.5)2(2.5)3(3.5)3122310.0017.024.012.047.028.01.0018.042.0试验指标为液化率,用yi表示,列于表-和表-的最后一列。表7-1 试验方案及结果分析试验号因 素试验结果液化率(%)ABCD1
4、234567891(10)112(50)223(90)331(1)2(4)3(7)1231231(20)2(35)3(50)2313121(1.5)2(2.5)3(3.5)3122310.0017.024.012.047.028.01.0018.042.0K1K2K341.087.061.013.082.094.046.071.072.089.046.054.0=189.013.729.020.34.327.331.315.323.724.029.715.318.0优水平A2B3C3D1Rj15.327.08.714.4主次顺序B A D C计算示例:因素A的第水平A1所对应的试验指标之和及其
5、平均值分别为:KA1=y1+y2+y3=0+17+24=41, KA1=13.7同理,对因素A的第水平A2和第水平A3,有KA2=y4+y5+y6=12+47+28=87, KA2=29KA3=y7+y8+y9=1+18+42=61, KA3=20.3由表-或表-可以看出,考察因素A进行的三组试验中(A1,A2,A3),B、C、D各水平都只出现了一次,且由于B、C、D间无交互作用,所以B、C、D因素的各水平的不同组合对试验指标无影响,因此,对A1、A2和A3来说,三组试验的试验条件是完全一样的。假如因素A对试验指标无影响,那么应该相等,但由上面的计标可知,实际上并不相等,显然,这是由于因素A的
6、水平变化引起的,因此,的大小反映了A1、A2和A3对试验指标影响的大小。由于液化率y越大越好,而,所以可判断A2为因素A的优水平。同理,可判断因素B、C、D的优水平分别为B3、C3、D1。所以,优水平组合为A2B3C3D1,即最优工艺条件为加水量A2=50ml/100g、加酶量B3=7ml/100g、酶解温度C3=50。C和酶解时间D1=1.5小时。二、确定因素主次顺序极差Rj按定义计算,如,同理可求出RC和RD. 计算结果列于表7-1中。比较Rj值可知RBRARDRC,所以试验因素对试验指标的影响的主次顺序为BADC。即加酶量影响最大,其次是加水量和酶解时间,而酶解温度的影响最小。三、绘制因
7、素与指标趋势图为了更直观地反映因素对试验指标的影响规律和趋势,用因素的水平作横坐标,试验指标的平均值()作纵坐标,画出因素与指标的关系图(即趋势图),如图7-2所示. (p137)趋势图可为进一步试验时选择因素水平指明方向.如对因素A,由图7-2可见,A2水平时,指标最高,但若能在A2附近再取一些水平(如40、60)作进一步试验,则有可能取得更高的指标;对D因素,若能取一些比D1更小的水平(如1.0和0.5)作进一步试验,也有可能得到更好的结果.以上三个步骤即为极差分析的基本程序与方法.四、说明与讨论1、计算结果的检验: 每一列的Kj之和应等于全部试验结果(即指标值)之和, 即,m为水平数,n
8、为试验总实施次数.2.因素的最优水平组合,在实际处理中是灵活的,即对于主要因素,一定要选最优水平;而对次要因素,则应权衡利弊,综合考虑其它条件进行水平选取,从而得到最符合实际生产的最优或较优生产工艺条件.3.例6-2的最优工艺条件A2B3C3D1并不在实施的9个试验之中.这表明优化结果不仅反映了已做的试验信息,而且反映了全面试验信息.因此,正交试验设计的部分实施方案反映了全面试验信息.4.例6-2得出的最优工艺条件,只有在试验所考察的范围内才有意义,超出这个范围,情况就可能发生变化。另外,只能说是“较优工艺条件”,而不能说是“最优工艺条件”.最好能根据趋势图做进一步试验,找出最靠近最优的工艺条
9、件.5.对已确定的最优工艺条件(如例6-2的A2B3C3D1)进行重复试验,验证其试验指标是否最优.7.2 多指标正交试验设计及其极差分析 在实际生产和科研试验中,所要考察的指标往往不止一个,这一类的试验设计叫做多指标试验设计.在多指标试验设计中,各指标之间可能存在一定的矛盾,如何兼顾各个指标,找出使每个试验都尽可能好的试验条件呢?换言之,应如何分析多指标试验设计的结果呢?常用的有两种方法:综合平衡法和综合评分法.下面举例说明综合平衡法的分析方法. 这种方法在试验方案安排和各指标计算分析方法上,与单指标试验完全一样.其步骤是先分别找出各个指标最优或较优的生产条件,然后将这些生产条件综合平衡,找
10、出兼顾每个指标都尽可能好的生产条件.例7-1 在油炸方便面的生产中,主要原料质量和主要工艺参数对产品的质量有影响。今欲通过正交试验确定最佳生产条件。一. 试验方案设计 1.确定试验指标评价方便面质量好坏的主要指标是: 脂肪含量(越低越好),水分含量(越高越好)和复水时间(越短越好)。 2挑因素,选水平,列出因素水平表根据专业知识和实际经验,确定试验因素和水平,如表7-2所示。表7-2 因素水平表水平 因素湿面筋值(%)A改良剂用量(%)B油炸时间(s)C油炸温度(C)D1232832360.050.0750.10707580150155160 3选正交表,设计表头,编制试验方案本试验是四因素三
11、水平试验,不考虑因素间的交互作用,因此,可应选L9(34)安排试验,表头设计和试验方案见表7-3(p140)。按上述方案实施后,将每一项试验指标都记录下来,见表7-3。注:对极差分析可以这样选正交表,但对方差分析应留有空列,以便估计试验误差.表7-3 试验方案及结果分析试验号因 素试验结果ABCD脂肪(%)水分(%)复水时间(s)1234567891(28)112(32)223(36)331(0.05)2(0.075)3(0.10)1231233(80)1(70)2(75)2311232(155)1(150)3(160)13232124.822.523.623.822.419.318.419.
12、020.72.13.82.02.81.72.72.52.02.33.53.73.03.02.22.83.02.73.6脂肪含量K1K2K370.965.558.167.063.063.660.266.467.967.063.164.4=194.523.621.819.422.321.321.220.122.122.622.321.021.5R4.21.12.51.3水分含量K1K2K37.97.26.87.47.56.99.06.86.18.96.86.2=21.92.632.402.272.472.502.303.002.272.032.972.272.07R0.360.200.970.90
13、复水时间K1K2K310.28.09.39.58.69.49.58.79.310.39.08.2=27.53.402.673.103.172.873.133.172.903.103.433.002.73R0.730.300.270.70二试验结果分析1计算每列各水平下每种试验指标的数据和(K1,K2,K3),及其平均值(),并计算极差R,填入表7-3中。2画出因素与各种指标的趋势图,如图7-3所示(p140)。3按极差大小列出各指标下各因素主次顺序:各因素主次顺序表试 验 指 标主-次脂肪含量(%)A C D B水分含量(%)C D A B复水时间(s)A D B C4初选最优工艺条件根据各指
14、标下的平均数据和,初步确定各因素的最优水平组合为: 对脂肪含量(%):A3B3C1D2 (脂肪含量越低越好) 对水分含量(%):A1B2C1D1 (水分含量越高越好)对复水时间(s):A2B2C2D3 (复水时间越短越好)5综合平衡确定最优工艺条件(难点)! 由于三个指标单独分析出来的最优条件并不一致,所以必须根据因素对三个指标影响的主次顺序,综合考虑,确定出最优条件。首先,把水平选取上没有矛盾的因素的水平定下来,即如果对三个指标影响都重要的某一因素,都是取某一水平时最好,则该因素就是选这一水平。在本试验中无这样的因素,因此我们只能逐个考察每一因素。对因素A:从主次顺序来看,对脂肪含量和复水时
15、间的影响都排在第一位为主要因素,而对水分含量的影响则排在第三位,属次要因素,因此,应以主要因素为主选因素的水平。从初选的最优水平组合中可以看出,对脂肪含量选A3为好,而对复水时间,则选A2为好。因为二者不一致,所以还须根据试验结果分析确定选A2还是A3。从表7-3可知,当取A2时,复水时间比取A3时缩短16.1%(有利),即(2.67-3.10)2.67100%=-16.1%,而脂肪含量只比取A3时增加11.0%(不利),即 (21.8-19.4)21.8100%=11.0%,且从水分含量指标来看,取A2也比取A3时更好,因此,应选取A2水平。注: 当取A3时,脂肪含量比取A2时降低12.4%
16、(有利),即(19.4-21.8)/19.4100%=12.4%,复水时间比取A2时增加13.9%(不利),即(3.10-2.67)/3.10100%=13.9%。综合平衡A不 利有 利A211.0%16.1%A313.9%12.4%对 “有利”部分,A2A3;对 “不利”部分,A2“不利”;选D3时,“不利”“有利”.并且D1 (有利)D3(有利之和绝对值),D1 (不利之和) D3(不利绝对值).因此,从定量分析来看,D应取D1,而不是取D3.那么,究竟如何决定D的水平呢?最后,应该再进行A2B2C1D1和A2B2C1D3两次试验,由试验结果决定D1好还是D3好!实践是检验真理的唯一标准!
17、7.3 混合型正交表的试验设计极差分析前面讨论的都是水平数相同的正交试验设计.但在实际工作中,有些试验受到设备、原材料和生产条件等限制.某些因素的水平选择受到制约,或者在有些试验中,要重点考察某个(或某些)因素需要多取几个水平,这时就会遇到水平数不同的正交试验设计.在这种情况下,通常有三种解决方法:一是直接选用合适的混合型正交表;二是采用拟水平法;三是采用拟因素法.我们现在只讨论第一种方法,即使用混合型正交表进行正交试验设计.例7-2 某油炸膨化食品的体积与油温、物料含水量及油炸时间有关,为确保产品质量,提出工艺要求。现通过正交试验设计寻求理想的工艺条件。一. 试验方案设计1.确定试验指标 本
18、试验的指标为油炸膨化食品的体积,体积越大越好.2.挑因素、选水平、制定因素水平表根据专业知识,制定因素水平表如7-4所示,因素A取4个水平,因素B和C各取2个水平,所以属于水平数不相等的正交试验设计.表7-4 因素水平表水平 因素油炸温度(C)A物料含水量(%)B油炸时间(s)C12342102202302402.04.030403.选正交表、设计表头、编制试验方案本试验宜选用L8(4124)正交表安排试验,表头设计时,把A因素放在第一列,其余两个因素可随意安排在四个二水平列中,比如依次排在第二、三列中,把所安排因素的各列的水平数字后标上相应因素的具体水平值,即得出试验方案,如表7-5所示.按
19、表7-5试验方案实施后,所得试验结果列于表7-5中的最后一列.表7-5 试验方案及结果分析试验号油温A含水量B时间C体积xi(cm3/100g)12345123456781(210)12(220)23(230)34(240)41(2.0)2(4.0)1212121(30)2(40)1221211221122112212112210.0208.0215.0230.0251.0247.0238.0230.0K1K2K3K4418.0445.0498.0468.0914.0915.0902.0927.0=1829.0209.0222.5249.0234.0228.5228.75225.5231.75
20、R40.00.256.2525.460.3558.875二. 试验结果分析1. 计算各列各水平下的K、及R由于各列的水平数不完全相同,所以K和的计算略有差异.第1列: 由于有四个水平数,所以要计算四个K与,每个K由二个数据相加得到,因此=K/2.例如:第2、3列:由于只有两个水平,所以只要计算两个K与,每个K由四个数据相加得到,因此=K/4.例如:按上述方法计算出各列各水平下的K、以及R值,列于表7-5中.2. 计算R的折算值R(极差R的折算)当因素的水平数相同时,因素的主次顺序完全由R决定.但当因素的水平数不同时,直接比较R是不行的.这是因为,若两个因素对试验指标有影响,一般来说,水平数多的
21、因素极差可能大一些.因此,要用一个系数把极差R折算后才能作比较.极差的折算公式如下:式中 -折算后的极差; R-因素的极差; r-该因素每个水平试验的重复数,r=; d-折算系数,与因素的水平数有关,其值见表7-6。表7-6 折算系数表 水平数m 2折算系数d0.71 0.52 0.45 0.40 0.37 0.35 0.34 0.32 0.31本例中,的折算如下:计算结果列于表-中.3.根据R大小确定因素的主次顺序主 -次A C B即油炸温度对实验指标的影响最大,其次是油炸时间,而物料含水量的影响最小。4.画出因素指标趋势图,如图7-4所示(p146)5.选各因素的最优水平及最优水平组合 比
22、较各因素各水平下的值(本例中越大越好),并参考因素指标趋势图,得出最优水平组合为A3B2C2或A3B1C2,即油炸温度230摄氏度,油炸时间40秒,物料含水量对试验指标影响很小,故取2%或4%都可以,视具体情况而定。由表7-5可见,若最优水平组合A3B1C2,则该试验即表中的第5号试验,实验指标值即膨化体积为251.03/100g,为表中所列最大值;若最优水平组合为A3B2C2,则需再实施一次该水平组合下的试验,作为验证。7.4 考察交互作用的正交试验设计及极差分析一、交互作用的概念 前面介绍的正交试验设计与试验结果的分析方法,都是指因素间没有(或不考虑)交互作用的情况,实际上,在许多试验中,
23、不仅因素对指标有影响,而且因素之间还会联合搭配起来对指标产生影响。所以,因素对试验产生的总效果,是由每一个因素对试验的单独作用再加上各个因素之间的搭配作用决定的。这种因素间的联合搭配对试验指标产生的影响作用,称为交互作用。 例如,我们要考虑化学反应的温度(A)与时间(B)对产品收率的影响,温度和时间都取二个水平,即和。在各AiBj组合条件的平均产品收率,可能有如下三种情况:(1)不论B因素取哪个水平,A2水平下收率总比A1水平高10;同样,不论A因素取哪个水平,B2水平下的收率总比B1水平下高5。在这种情况下,一个水平的好坏或好坏程度不受另一个因素水平的影响,这种情况称为因素A与B之间无交互作
24、用。(2)在B1水平下A2比A1的收率高,但在B2水平下,A1比A2的收率高。这种一个因素水平的好坏或好坏程度受到另一因素水平制约的情况,称为因素A由于因素B存在交互作用,一般用AB表示。(3)不论B因素取哪个水平,A2水平的收率总比A1水平下高,但高的程度不等,这也说明因素A与B存在交互作用。A1A2B17585B28090(1) A与B间无交互作用(平行线)A1A2B17585B28065(2) A与B间有交互作用(AB)A1A2B17585B28095(3)A与B间存在交互作用(AB)图7-4 A与B 间的交互作用情况 事实上,因素之间总是存在着交互作用的,这是客观存在的普遍现象,只不过
25、交互作用的程度不同而已。一般的,当交互作用很小时,就认为不存在交互作用。因素间的交互作用对试验指标的影响,可能是正的,也可能是负的。有人说:“中国人一个人像一条龙,三个人像一条虫;日本人一个人像一条虫,三个人像一条龙。”这说明中国人之间的交互作用常常产生负面效应。(一个和尚挑水喝,二个和尚抬水喝,三个和尚没水喝。团结就是力量,集体主义精神)在试验设计中,表示因素A、B间的交互作用记作AB,称作一级交互作用;表示因素A、B、C之间的交互作用记作ABC,称作二级交互作用;依次类推,还有三级、四级交互作用。二级和二级以上的交互作用称为高级交互作用。在试验设计中,通常忽略高级交互作用。2交互作用的处理原则 处理交互作用的总原则是,将交互作用当作因素看待,并将交互作用安排在能考察交互
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1