正交实验设计

1、正交实验设计当析因设计要求的实验次数太多时，一个非常自然的想法就是从析因设计的水平组合中，选择一部分有代表性水平组合进行试验。因此就出现了分式析因设计（fractionalfactorialdesigns）,但是对于试验设计知识较少的实际工作者来说，选择适当的分式析因设计还是比较困难的。正交试验设计（Orthogonalexperimentaldesign）是研究多因素多水平的又一种设计方法，它是根据正交性从全面试验中挑选出部分有代表性的点进行试验，这些有代表性的点具备了均匀分散，齐整可比”的特点，正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家

2、田口玄一将正交试验选择的水平组合列成表格，称为正交表。例如作一个三因素三水平的实验，按全面实验要求，须进行33=27种组合的实验，且尚未考虑每一组合的重复数。若按L9（3）3正交表按排实验，只需作9次，按Li8（3）7正交表进行18次实验，显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。1 .正交表正交表是一整套规则的设计表格，用工式/）。L为正交表的代号，n为试验的次数，t为水平数，c为列数，也就是可能安排最多的因素个数。例如L9（34）,（表11）,它表示需作9次实验，最多可观察4个因素，每个因素均为3水平。一个正交表中也可以各列的水平数不相等，我们称它为混合型正

3、交表，如L8（4>24）（表12）,此表的5列中，有1列为4水平，4列为2水平。根据正交表的数据结构看出，正交表是一个n行c列的表，其中第j列由数码1,2,S组成，这些数码均各出现N/S次，例如表11中，第二列的数码个数为3,S=3,即由1、2、3组成，各数码均出现明3=邛=3次。表nl式厂）表引号13411i11212213iJ33.4212352A3I52?上73133R333P3331表广X公）表民喳号-班片123451111112122223211224222IL531A1W6331Q17A122LS42112表134（2）正交表遍彳列EI2345d711111112111222

4、2312211224122221152111316212212172211221£221211口正交表具有以下两项性质：（1）每一列中，不同的数字出现的次数相等。例如在两水平正交表中，任何一列都有数码“1与“2；且任何一列中它们出现的次数是相等的；如在三水平正交表中，任何一列都有“17”27”3；且在任一列的出现数均相等。（2）任意两列中数字的排列方式齐全而且均衡。例如在两水平正交表中，任何两列（同一横行内）有序对子共有4种：（1,1）、（1,2）、（2,1）、（2,2）。每种对数出现次数相等。在三水平情况下，任何两列（同一横行内）有序对共有9种，1.1、1.2、1.3、2.1、2.

5、2、2.3、3.1、3.2、3.3,且每对出现数也均相等。以上两点充分的体现了正交表的两大优越性，即均匀分散性，整齐可比通俗的说，每个因素的每个水平与另一个因素各水平各碰一次，这就是正交性。2 .交互作用表每一张正交表后都附有相应的交互作用表，它是专门用来安排交互作用试验。表14就是L8（27）表的交互作用表。表145（2，）表的交互作用表号列4、12345672325416167457654123321安排交互作用的试验时，是将两个因素的交互作用当作一个新的因素，占用一列，为交互作用列，从表14中可查出L8（27）正交表中的任何两列的交互作用列。表中带（）的为主因素的列号，它与另一主因素的交

6、互列为第一个列号从左向右，第二个列号顺次由下向上，二者相交的号为二者的交互作用列。例如将A因素排为第（1）列，B因素排为第（2）列，两数字相交为3,则第3列为AXB交互作用列。又如可以看到第4列与第6列的交互列是第2歹U,等等。3 .正交实验的表头设计表头设计是正交设计的关键，它承担着将各因素及交互作用合理安排到正交表的各列中的重要任务，因此一个表头设计就是一个设计方案。表头设计的主要步骤如下：(1)确定列数根据试验目的，选择处理因素与不可忽略的交互作用，明确其共有多少个数，如果对研究中的某些问题尚不太了解，列可多一些，但一般不宜过多。当每个试验号无重复，只有1个试验数据时，可设2个或多个空白

7、列，作为计算误差项之用。(2)确定各因素的水平数根据研究目的，一般二水平(有、无)可作因素筛选用；也可适用于试验次数少、分批进行的研究。三水平可观察变化趋势，选择最佳搭配；多水平能以一次满足试验要求。(3)选定正交表根据确定的列数(c)与水平数(t)选择相应的正交表。例如观察5个因素8个一级交互作用，留两个空白列，且每个因素取2水平，则适宜选Li6(215)表。由于同水平的正交表有多个，如L8(27)、Li2(211)、Li6(215),一般只要表中列数比考虑需要观察的个数稍多一点即可，这样省工省时。(4)表头安排应优先考虑交互作用不可忽略的处理因素，按照不可混杂的原则，将它们及交互作用首先在

8、表头排妥，而后再将剩余各因素任意安排在各列上。例如某项目考察4个因素A、B、C、D及AXB交互作用，各因素均为2水平，现选取L8(27)表，由于AB两因素需要观察其交互作用，故将二者优先安排在第1、2歹U,根据交互作用表查得AXB应排在第3歹U,于是C排在第4歹U,由于AXC交互在第5歹U,BXC交互作用在第6歹U,虽然未考查AXC与BXC,为避免混杂之嫌，D就排在第7歹U。表151。丁)表头设计列号123453因素与交互作用ABCD(5)组织实施方案根据选定正交表中各因素占有列的水平数列，构成实施方案表，按实验号依次进行，共作n次实验，每次实验按表中横行的各水平组合进行。例如L9(34)表，

9、若安排四个因素，第一次实验A、B、C、D四因素均取1水平，第二次实验A因素1水平，B、C、D取2水平，第九次实验A、B因素取3水平，C因素取2水平,D因素取1水平。实验结果数据记录在该行的末尾。因此整个设计过程我们可用一句话归纳为：因素顺序上列、水平对号入座，实验横着作4 .二水平有交互作用的正交实验设计与方差分析例8某研究室研究影响某试剂回收率的三个因素，包括温度、反应时间、原料配比，每个因素都为二水平，各因素及其水平见表16。选用L8(27)正交表进行实验，实验结果见表17。5表1。因素与水平因素水平12A温度七60SOB反应时间h253：.5C原料配比11:112:1表17票试剂回收率的

10、正交实嘘口（2，）去结果试虻号kk2B3AXB4-C563试脂给果回收率了乂i11Ii111U2L1r：2.22.餐953122品.122914L222211945)12ia1291621221219672211.221£3822鸟1；11'128S3663633523513593SP3581(53723733(533fi53fi5?1220尔聋-a-6,占J喀VS.&50,。缸50.54.54.5首先计算Ij与山，Ij为第j列第1水平各试验结果取值之和，山为第j列第2水平各试验结果取值之和。然后进行方差分析。过程为：/=58-=M6.0求：总离差平方和-3TI/各列

11、离差平方和SS=，本例各列离均差平方和见表10最底部一行。即各空列SS之和。即误差平方和=£&&至负=05+45+45=95自由度v为各列水平数减1,交互作用项的自由度为相交因素自由度的乘积。分析结果见表18表三种因素对某试剂回收率哥晌的正较实验方差分析表变异来源离差平方和S写自由度U均方MSF值P值A801S.02.2530.2102B1S.0112.05右80.0973AXB5k0150.015.390.02S5C60U160.519.10.0222误差Q.53316总变异146.0?从表18看出，在a=0.05水准上，只有C因素与AXB交互作用有统计学意义，其余

12、各因素均无统计学意义，A因素影响最小，考虑到交互作用AXB的影响较大，且它们的二水平为优。在C2的情况下，有B1A2和Bi,Ai两种组合状况下的回收率最高。考虑到B因素影响较A因素影响大些，而B中选Bi为好，故选A2B10这样最后决定最佳配方为A2B1C2,即80C,反应时间2.5h,原料配比为121。如果使用计算机进行统计分析，在数据是只需要输入试验因素和实验结果的内容，交互作用界的内容不用输入，然后按照表头定义要分析的模型进行方差分析。附录1:常用正交表(1)L4(23)列号试验号123111121223二2124221L8(27)列号1234567试验号11111111211122122

13、31221122412222115212121262口221217221122182212112(3)L 12(211)列号1234567891011试验号111111111111211111222222J311222111222412122122112512212|21212161I2I212212117212211221218212122211129211222122111022211112212112212121112212221121212214(4)L9(3)列号1234试验号1111121222313334212352231623127313283213933216(5)L16(4

14、5)列号试验号1234511111121222231333341444452123462214372341218243219313421032431113312412342131341423144231415432411644132(6)L25(56)列号试验号123456111111121222223133333414444451555556212345722一3451823451292451231025123411313524123241351333524114341352153524131641425317425314184314251944253120453142215154322252一1543235321542454321525554321，、，4、L8(4X2)列号实验号123451|1|1|1112|12222321122422211531212632121741221842112(8)L12(3X24)列号试验号123451111112111