正交实验法原理与应用

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,正交实验法原理与应用,杨明军,2014-1-24,目录,一、正交实验法意义和原理,二、正交表,三、,正交实验设计的基本步骤,四、结果分析,一、正交实验法意义和原理,凡是要做实验，就存在着如何安排实验和如何分析实验结果的问题。这就是做实验的方法问题。一项科学的安排实验的方法应能做到以下两点：,1,）在实验安排上尽可能地减少实验次数；,2,）在进行较少次数实验的基础上，能够利用所得到的实验数据，分析出指导实验的正确结论，并得到较好的结果。因此，科学的实验方法是使我们的工作达到快捷高效的一个工具。,1.1,、正交实验法意义,在实验研究中，我们一般要考察原料种类、原料数量、反应温度、反应时间、反应压力等多个影响因素，若进行全面实验 ，则实验的规模将很大 ，往往因实验条件的限制而难于实施 。例如，一实验有,6,个因素： 每因素取,5,个水平，全面实验就需要,5,6,=15625,个组合。因素是指参与实验并对其结果有影响的要素或对象，水平是指因素的变化状态或用量。,1.2,、,正交实验法原理,19,世纪,20,年代，英国统计学家,R. A. Fisher,首先在马铃薯肥料实验当中，运用排列均衡的拉丁方，解决了实验时的不均匀实验条件，获得成功，并创立了,“,实,验设计,”,这一新兴学科。,“,均衡分布,”,思想在,20,世纪,50,年代应用于工业领域，,60,年代应用于农业领域，使正交实验在科研生产实际中得到推广,。,1.2.1,正交拉丁方,用从,1,开始的,n,个连续正整数排成,n,行,n,列的方阵，如果每一行和每一列都 没有重复的数，就称为一个,n,阶拉丁方。两个,n,阶拉丁方在同一位置上的数依次配置成对时，如果这两个有序数对恰好各不相对（一般处理方法把当中某些行或列对调）（这种相同即经过有限次旋转和镜像对称后不重合）。,在实验安排中 ，每个因素在研究的范围内选几个水平，就好比在选优区内打上网格 ，如果网上的每个点都做实验，就是全面实验。,3,个因素的选优区可以用一个立方体表示,(,图,1),，,3,个因素各取,3,个水平，把立方体划分成,27,个格点，反映在 图,1,上就是立方体内的,27,个“,.”,。若,27,个网格点都实验，就是全面实验，其实验方案如表,1,所示。,表,1,3,因素,3,水平全面实验方案,C1,C2,C3,A1,B1,A1B1C1,A1B1C2,A1B1C3,B2,A1B2C1,A1B2C2,A1B2C3,B3,A1B3C1,A1B3C2,A1B3C3,A2,B1,A2B1C1,A2B1C2,A2B1C3,B2,A2B1C1,A2B1C2,A2B1C3,B3,A2B2C1,A2B2C2,A2B2C3,A3,B1,A3B1C1,A3B1C2,A3B1C3,B2,A3B2C1,A3B2C2,A3B2C3,B3,A3B3C1,A3B3C2,A3B3C3,正交实验法就是从选优区全面实验点（水平组合）中挑选出有代表性的部分实验点（水平组合）来进行实验。图,1,中标有实验号的九个“,()”,，就是利用正交表,L,9,(3,4,),从,27,个实验点中挑选出来的,9,个实验点。即：,关于正交的直观印象,数据点分布是均匀的,每一个面都有,3,个点,每一条线都有,1,个点,图,1,1.2.2,正交实验法,正交实验法是用来科学地设计多因素实验的一种方法。它利用一套规格化的正交表安排实验，得到的实验结果再用数理统计方法进行处理，使之得出科学结论。正交表是实验设计的基本工具，它是根据均衡分布的思想，运用组合数学理论构造的一种数学表格，均衡分布性是正交表的核心。,2., 正交,表 正交拉丁方的自然推广,由于正交设计安排实验和分析实验结果都要 用 正交表，因此，我们先对正交表作一介绍。,安排,4,因素,3,水平的实验，编上实验号，列成另外一种形式，见正交表,L,9,(,3,4,),（表,2,）,。,可以由此得到系列正交表。,二、正交表,表,2,因素,实验号,A,B,C,D,1,1,1,1,1,2,1,2,2,2,3,1,3,3,3,4,2,1,2,3,5,2,2,3,1,6,2,3,1,2,7,3,1,3,2,8,3,2,1,3,9,3,3,2,1,常用的正交表已由数学工作者制定出来，供进行正交设计时选用。,2,水平正交表除,L,8,(2,7),外，还有,L,4,(2,3,),、,L,16,(2,15,),等；,3,水平正交表有,L,9,(3,4,),、,L,27,(,3,13,),等。,表,2,是一张正交表，记号为,L,9,(3,4,),，其中“,L”,代表正交表；,L,右下角的数字“,9”,表示有,9,行 ，用这张正交表安排实验包含,3,个处理,(,水平组合,),；括号内的底数“,3”,表示因素的水平数，括号内,3,的指数“,4”,表示有,4,列 ，也指安排的因素数，用这张正交表最多可以安排,4,个因素,3,个水平。,2.2,正交表的表示符号,正交表记号所表示的含义归纳如下：,L,n (t,q,）,式中：,L,为正交表符号，是,Latin,的第一个字母；,n,为实验次数，即正交表行数；,t,为因素的水平数，即,1,列中出现不同数字的个数,;q,为最多能安排的因素数，即正交表的列数。,正交表中,1,列可以安排,1,个因素，因此它可安排的因素数可以小于或等于,q,，但不能大于,q,。,括号内的,t,q,表示,q,个因素、每个因素,t,个水平全面实验的水平组合数（即处理数）。因为安排因素个数不能大于,q,，所以,n /t,q,为最小部分实施。,显然,，,L,4,(2,3,),是最简单的正交表，有,4,列,3,行用它最多能安排,3,个因素,2,水平的实验。部分实验为,4,次，全面实验为,8,次，最小部分实施为,1/2,，即用它安排实验可比全面实验少做,1/2,。所以，当实验因素数,q,及每个因素的水平数,t,增加时,n/t,q,则下降，节省实验,次数的效果更明显。,2.3,常用正交表的分类及特点,1,、标准表（,相同水平正交表）,2,水平,:L,4,(2,3,）,L,8,(2,7,）,L,16,(2,15,），,3,水平,:L,9,(3,4,），,L,27,（,3,13,），,L,81,(3,40,），,4,水平：,L,16,(4,5,）,L,64,(4,21,）,L,256,(4,85,），,5,水平：,L,25,（,5,6,）,L,125,(5,31,）,L,625,(5,156,），,各列中出现的最大数字相同的正交表称为相同水平正交表。如,L,4,(2,3,),、,L,8,(2,7,),、,L,12,(2,11,),等各列中最大数字为,2,，称为两水平正交表；,L,9,(3,4,),、,L,27,(3,13,),等各列中最大数字为,3,，称为,3,水平正交表。,2,、非标准表（,混合水平正交表）,各列中出现的最大数字不完全相同的正交表称为混合水平正交表。如,L,8,(42,4,),表中有一列最大数字为,4,，有,4,列最大数字为,2,。也就是说该表可以安排一个,4,水平因素和,4,个,2,水平因素。再如,L,16,(4,4,2,3,),，,L,16,(42,12,),等都混合水平正交表。,2.4,正交表的基本性质,任何一张正交表都有如下三个特性：,2.4.1,正交性,1,、,任一列中，不同数字出现的次数相等,例如,L,8,(2,7,),中不同数字只有,1,和,2,，它们各出现,4,次；,L,9,(3,4,),中不同数字有,1,、,2,和,3,，它们各出现,3,次,。,2,、,任两列中，同一横行所组成的数字对出现的次数相等,例如,L,8,(2,7,),中,(1, 1), (1, 2), (2, 1), (2, 2),各出现两次；,L,9,(3,4,),中,(1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3),各出现,1,次。即每个因素的一个水平与另一因素的各个水平互碰次数相等，表明任意两列各个数字之间的搭配是均匀的。,由正交表的正交性可以看出：,正交表各列的地位是平等的，表中各列之间可以互相置换，称为列间置换；,正交表各行之间也可相互置换，称行间置换；,正交表中同一列的水平数字也可以相互置换，称水平置换。,上述,3,种置换即正交表的,3,种初等置换。经过初等置换所能得到的一切正交表，称为原正交表的同构表或等价表，显然，实际应用时，可以根据不同需要进行变换。,2.4.2,代表性。,代表性的含义之一，在于正交表的正交性中：,任一列的各水平都出现，使得部分实验中包含所有因素的所有水平。,任意,2,列间的所有组合全部出现，使任意两因素间都是全面实验。因此，在部分实验中，所有因素的所有水平信息及两因素间的所有组合信息都无一遗漏。这样，虽然安排的是部分实验，却能够了解全面实验的情况，从这个意义上,讲可以代表全面实验。,因为正交性，使部分实验点必然均衡地分布在全面实验的实验点中。,所谓均衡分散，是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是均匀的 。 由 图,1,可以看出，在立方体中 ，任一平面内都包含,3,个“,()”,， 任一直线上都包含,1,个“,()”,，因此 ，这些点代表性强 ，能够较好地反映全面实验的情况。,2.4.3,综合可比性。,反映在正交性当中：,任一列各水平出现的次数都相等。,任,2,列间所有可能的组合出现的次数都相等。因此使任一因素各水平的,实验条件相同。这就保证了在每列因素各个水平的效果中，最大限度地排除其他因素的干扰，突出本列因素的作用，从而可以综合比较该因素不同水平对实验指标的影响。这种性质称为综合可比性或整齐可比性。,如在,A,、,B,、,C 3,个因素中，,A,因素的,3,个水平,A1,、,A2,、,A3,条件下各有,B,、,C,的,3,个不同水平，即：,在这,9,个水平组合中，,A,因素各水平下包括了,B,、,C,因素的,3,个水平，虽然搭配方式不同，但,B,、,C,皆处于同等地位，当比较,A,因素不同水平时，,B,因素不同水平的,在这,9,个水平组合中，,A,因素各水平下包括了,B,、,C,因素的,3,个水平，虽然搭配方式不同，但,B,、,C,皆处于同等地位，当比较,A,因素不同水平时，,B,因素不同水平的效应相互抵消，,C,因素不同水平的效应也相互抵消。所以,A,因素,3,个水平间具有可比性。同样，,B,、,C,因素,3,个水平间亦具有可比性。,根据以上两个特性，我们用正交表安排的实验，具有均衡分散和整齐可比的特点。,正交表的,3,个基本性质中，正交性即均衡性是核心，是基础，代表性和综合可比性是正交性的必然结果，从而使正交表得以具体应用。,正交表集其,3,个性质于一体，成为正交实验设计的有效工具，用它来安排实验，也必然具有,“,均衡分散，整齐可比,”,的特性，代表性强，效率也高。因而，实际应用越来越广。,三、,正交实验法的基本步骤,正交实验法的基本步骤主要有又下几步：,第一步，明确实验目的，确定考核指标。,第二步，挑因素，选水平。,第三步，选择合适的正交表。,第四步，进行表头设计。,第五步，确定实验方案。,为了更好地说明问题，我们结合一个例子来说明。,例,1,用乙酸和环己醇酯化来制备乙酸环己酯的反应,简单实验步骤是：将原料加入四口瓶中加热回流，反应结束冷至室温，反应液转移至分流漏斗中碱洗分层，有机层经柱色谱法纯化得到乙酸环己酯。,3.1,明确实验目的，确定考核指标,实验目的，就是通过正交实验要想解决什么问题。,我们做实验的目的一般就是提高产品质量或收率，降低成本 或者查找生产中出现问题的原因,。,3.2,挑因素，选水平,实验指标确定之后，就可着手分析影响指标的各种因素。影响实验结果的因素往往是很多的，究竟挑选哪些因素作为考察因素，需要根据以往的经验和实验条件来确定。排除那些对指标影响不大，或已经掌握得比较好的因素，挑选出那些对指标可能影响较大，但又没有把握的因素来进行考察。特别注意不能把重要因素固定，否则，由于重要因素固定在不适当的水平上，使实验得不到应有效果。,结合例,1,根据文献报道，在实验前分析影响反应的可能因子为：催化剂的种类（,A,）、反应温度（,B/ ),、反应时间（,C/h),、带水剂的用量（,D/ml),及反应原料比（,E,）,5,个因素,.,每个因素在其实验范围内取的实验点叫做该因素的水平。对挑出的因素，各水平选在何处是值得很好研究的。水平选得好，通过较少的实验次数，就能迅速地找到最优的实验条件。该实验中每个因素取,4,个水平,(,数据如表,3),。,表,3,3.3,、选出合适的正交表,选用哪种正交表，一般根据因素和水平多少以及实验工作量大小而定。,总原则：能容纳所有考察因素，又使实验号最小。例如，若因素都是二水平的，当有,3,个因素时，可用,L4,（,2,3,）（当然也可用,L8,（,2,7,），但实验工作量大），当有,4,7,个因素时，一般用,L8,（,2,7,）（也可用,L16,（,2,15,）等等。,但必须注意，因素水平表中的水平数与所选正交表中的水平数要完全一致，因素水平表中的因素个数要小于或等于正交表中的列数。,在讨论实验方案时，挑选因素和水平及选用正交表有时是结合进行的，例如，原考虑四因素四水平实验，应选,L16,（,4,5,），需要做,16,次实验。为减少实验次数，改为四因素三水平，选,L9,（,3,4,），做,9,次实验就够了。例,1,为,5,因素,4,水平实验，可选用,L16,（,4,5,）（如表,4,）（若全面实验则要,4,5,=1024,次实验）。,表,4,表,5,3.4,、作表头设计,把五个因素,A,、,B,、,C,、,D,、,E,分别列在正交表,L16,（,4,5,）表头的任意五列上。例如，依次放在第,1,，,2,，,3,，,4,，,5,列上。这一步通常叫做“表头设计”，列成表,5,。,列号,1,2,3,4,5,因素,A,B,C,D,E,3.5,列出实验方案,按因素水平表把,A,、,B,、,C,、,D,、,E,占有的各列中对应的“,1”,、“,2”,、“,3”,、“,4”,、“,5”,换成因素的具体水平。总共,16,个实验（如表,6,）。,因素,催化剂的种类（,A,）,反应温度（,B/ ),反应时间（,C/h),带水剂的用量（,D/ml),反应原料比（,E,）,实验号,1,浓硫酸（,F,）,65,1,4,1,2,浓硫酸（,F,）,70,1.5,5,1.5,3,浓硫酸（,F,）,75,2,6,2,4,浓硫酸（,F,）,80,2.5,7,2.5,5,维生素,C,片（,G,）,65,1.5,6,2.5,6,维生素,C,片（,G,）,70,1,7,2,7,维生素,C,片（,G,）,75,2.5,4,1.5,8,维生素,C,片（,G,）,80,2,5,1,9,十二水合硫酸铁铵（,H,）,65,2,7,1.5,10,十二水合硫酸铁铵（,H,）,70,2.5,6,1,11,十二水合硫酸铁铵（,H,）,75,1,5,2.5,12,十二水合硫酸铁铵（,H,）,80,1.5,4,2,13,催化剂,I,65,2.5,5,2,14,催化剂,I,70,2,4,2.5,15,催化剂,I,75,1.5,7,1,16,催化剂,I,80,1,6,1.5,表,6,四、,正交实验的结果分析,现在我们还结合例,1,来进行结果分析。实验结果如表,7,：,因素,催化剂的种类（,A,）,反应温度（,B/ ),反应时间（,C/h),带水剂的用量（,D/ml),反应原料比（,E,）,收率（,%,）,实验号,1,浓硫酸（,F,）,65,1,4,1,47.4,2,浓硫酸（,F,）,70,1.5,5,1.5,55.6,3,浓硫酸（,F,）,75,2,6,2,57.6,4,浓硫酸（,F,）,80,2.5,7,2.5,57.2,5,维生素,C,片（,G,）,65,1.5,6,2.5,72.6,6,维生素,C,片（,G,）,70,1,7,2,71.4,7,维生素,C,片（,G,）,75,2.5,4,1.5,72.6,8,维生素,C,片（,G,）,80,2,5,1,73.6,9,十二水合硫酸铁铵（,H,）,65,2,7,1.5,80.6,10,十二水合硫酸铁铵（,H,）,70,2.5,6,1,80.5,11,十二水合硫酸铁铵（,H,）,75,1,5,2.5,82.5,12,十二水合硫酸铁铵（,H,）,80,1.5,4,2,83.6,13,催化剂,I,65,2.5,5,2,91.4,14,催化剂,I,70,2,4,2.5,92.5,15,催化剂,I,75,1.5,7,1,95.8,16,催化剂,I,80,1,6,1.5,95.3,4.1,、计算综合平均值,按例,1,计算的结果如表,8,A,B,C,D,E,以表,8,为例。用,a,表示因素,A,取第一水平时相应的实验结果之和，用,a,表示因素,A,取第,2,水平时相应的实验结果之和，用,a,、,a,分别表示因素,A,取第,3,、,4,水平时相应的实验结果之和。即,a=47.4+55.6+57.6+57.2=217.8,a=72.6+71.4+72.6+73.6=290.2,a=80.6+80.5+82.5+83.6=327.2,a=91.4+92.5+95.8+95.3=375,为比较因素,A,不同水平的好坏，求平均值,K,A1,=a/4=217.8/4=54.45,K,A2,=a/4=290.2/4=72.55,K,A3,=a/4=327.2/4=81.8,K,A4,=a/4=375/4=93.75,K,A1,、,K,A2,、,K,A3,、,K,A4,分别称为因素,A,相应水平的综合平均值。表,8,中第一列中的数值就是这样算的。根据此法分别算出其它因素相应水平的综合平均值。,4.2,、,画出指标与因素关系图,为了更直观地反映各因素对乙酸环己酯合成的实验指标的影响规律和趋势，分别以因素的水平作为横坐标，综合平均值为纵坐标，得到了各因素对收率的影响趋势图，详见图,2,6,，由表,8,或图,2,6,可看出催化剂,I,效果好，反应温度越高越好，反应时间,1.5,小时最好，带水剂的用量,6ml,最好，反应原料比环己醇：乙酸,=1,：,2.5,为最好，确定最优水平组合为,A4B4C2D3E4,。,图,2,图,3,图,4,图,5,图,6,4.3,、分析因素的主次,从图,2,6,又可看出，当因素取不同水平时，点子散布的范围大，即点子上升（或下降）的幅度大，该因素就是影响指标的主要因素，反之，即为次要因素。我们用一个数（极差）来描述分散程度的大小。因素,A,的各个综合平均值中最大的减去最小的称为因素,A,的极差，用,R,表示。,计算如下：,R,A,= K,A1,、,K,A2,、,K,A3,、,K,A4,中最大值减去最小值,=93.750-54.450=39.300,R,B,=77.425-73.000=4.425,R,C,=76.900-74.150=2.75,R,D,=76.500-74.025=2.475,R,E,=76.200-74.325=1.875,表,8,的最末一行就是各因素的极差。极差大者，为重要因素,;,极差小者，为次要因素。依极差大小，区分因素主次如下：,ABCDE,4.4,、选取最优水平组合,如果要寻找使指标越大越好的条件，就选取使各因素综合平均值最大的水平组合，作为最优水平组合。如果要寻找使指标越小越好的条件，就选取使各因素综合平均值为最小的水平组合，作为最优水平组合。在例,1,中已提到，最优水平组合是,A4B4C2D3E4,。值得注意的是，这个条件并不在做过的,16,个实验之内。这是正交实验的优点，它可以指出寻找最优实验条件的途径。,4.5,、验证实验,对最优水平组合应该做验证实验。在例,1,中将最优水平组合,A4B4C2D3E4,和,16,个实验中最好的,A4B3C2D4E1,（第,15,号实验）做了对比实验。,对比实验的结果是，,A4B4C2D3E4,的收率为,96.5%,，,A4B3C2D4E1,的收率为,95.8%,，证实选出的实验条件确实是最好的。,对最优水平组合应做验证实验。,报告结束，谢谢,！,