正交试验设计原理

1、正交实验设计1 .概述任何生产部门，任何科学实验工作，为达到预期目的和效果都必须恰当地安排实验工作，力求通过次数不多的实验认识所研究课题的基本规律并取得满意的结果。例如为拟定一个正确而简便的分析方法，必然要研究影响这种分析方法效果的种种条件，诸如试剂浓度和用量、溶液酸度、反应时间以及共存组分的干扰等等。同时，对于影响分析效果的每一种条件，还应通过试验选择合理的范围。在这里，我们把受到条件影响的反系方法的准确度、精密度以及方法的效果等叫做指标；把试验中要研究的条件叫做因素；把每种条件在试验范围内的取值（或选取的试验点）叫做该条件的水平。这就是说我们常常遇到的问题可能包括多种因素，各种因素又有不同

2、的水平，每种因素可能对分析结果产生各自的影响，也可能彼此交织在一起而产生综合的效果。正交试验设计就是用于安排多因素实验并考察各因素影响大小的一种科学设计方法。它始于1942年,之后在各个领域里都得到很快的发展和广泛应用。这种科学设计方法是应用一套已规格化的表格一一正交表来安排实验工作，其优点是适合于多种因素的实验设计，便于同时考查多种因素各种水平对指标的影响通过较少的实验次数，选出最佳的实验条件，即选出各因素的某一水平组成比较合适的条件，这样的条件就所考查的因素和水平而言，可视为最佳条件。另一方面，还可以帮助我们在错综复杂的因素中抓住主要因素，并判断那些因素只起单独的作用，那些因素除自身的单独

3、作用外，它们之间还产生综合的效果。数理统计上的实验设计还能给出误差的估计。2 .试验设计的基本方法全面试验法正交设计的方法，首先应根据实验的目的，确定影响实验结果的各种因素，选择这些影响因素的试验点，进而拟出实验方案，之后按所拟方案进行实验并对实验结果作出评估。必要时再拟出进一步的实验方案，使实验工作更趋完善，所得结果也更为可靠。如在研究某一显色反应时，为选择合适的显色温度、酸度和显色完全的时间，可作如下的试验安排。首先确定上述三因素的实验范围：显色温度：25之一35C（温度以A表示）酸浓度：L（酸浓度以B表示）显色时间：1030min（时间以C表示）其次确定每种因素在上述实验范围内各取的水平

4、数（如各取三个水平）因素A的三个水平分别以Ai,A2,A表不；因素B的三个水平分别以Bi,B2,B3表小；因素C的三个水平分别以C,C2,C3表小；然后将显色试验的因素、水平列为下表。二、因素ABC水平f、温度（t/C）酸浓度(C/molL-1)时间(t/min)125102302033530这是一个三因素三水平的试验问题，对这样的试验工作可做如下的安排ABCABGABCABGABCABCABC3ABGA3B1C3AB2cAB2cA3B2CAB2C2AB2C2A362GAB2c3AB2c3A3B2C3ARGARCARCARGARGARCA1BC3ARGARC即三因素水平的试验共27种组合（33

5、=27）,按上组合方式做完27次试验后自然可得出在所确定的因素和水平下的最佳显色条件。这种全面试验的方法，对事物的内部规律剖析得十分清楚，但却费时费事。假如我们还需要对实验精密度，对试验误差的大小做出估计，则每一试验至少应重复一次。即应做54次实验。如果在讨论六因素而每种因素均取5个水平时，则全面试验的数目是56=15625次，这里还未包括为了给出误差估计所需的重复试验次数，显然这是难以付诸实施的。当考察的因素，水平数越多，在试验中所有可能的搭配也更多，要逐个地进行试验，显然是不可能的。这就提出了合理地设计和安排试验的问题。提出了通过较少量的试验次数以获得理想的实验条件取得最佳的试验效果，并对

6、试验结果做出科学评估的问题。对于上述试验，一种习惯的试验方法是简单比较法。简单比较法这种方法首先固定因素A、B为某一水平（如Ai、Bi）,改变C以获得在Ai、Bi时C的最佳水平（设为在其下以“-。C；AiBiCC3然后固定A为Ai,C为G,改变B以获得在Ai、G时B的最佳水平（设为B3）AiC23再固定B为B3,C为G,改变A以获得在B、G时的最佳水平（设为A）这样可以认为A2B3c为较佳的显色条件，即简单比较法经过9次试验也能获得较佳的试验条件，但却存在以下缺点：当各因素之间交互影响较大时，ABG不认为是最佳试验条件。它未能保证三因素中任何两因素的不同水平之间相碰一次因而上不均衡的，它提供的

7、信息也是不丰富的。在不做重复试验的情况下，不能给出误差的估计。如何保持这种方法试验次数少的优点而又能避免上述缺点呢，可采用正交设计的方法来解决。在这9次试验中实际上有两次试验是在相同条件下的重复试验（ARG和ARG）,所以只有7次属不同条件下的实验，另一方面还可看出各因素、各水平出现的机会是不均衡的，其中AnC各出现了7次；R、Ci各出现了4次；而入、A、C、C、B2却只出现了一次，显然，它们的出现的机会是很不均衡的。简单比较法认为最佳的分析条件是ABsG,但在试验过程中G是在AB条件下与C和C3相比，是最佳的一个条件水平，至于因素A、B取其他水平时是否也得出同样的结论，却未做过实验，也不能得

8、出同样的结论，故上述的条件不能视为最佳的显色条件，而只能是最佳条件的一种估计。导致上述几种问题的原因是简单比较法中各因素各水平的搭配不是均衡分散的，只能在同一批试验中做单因素比较，而在不同批数的试验之间却无法进行比较。2. 3正交设计法试验设计是数理统计中的一个重要内容，正交设计是利用预先编制好的正交表来合理的安排多因素试验，以便通过少量的试验次数来获得满意的结果，同时对试验数据进行统计分析。现在对三因素三水平的试验做如下的安排，首先只考虑A、B两因素，起全面实验应作9次，如下表所示。BiB2B3AAiBiAB2A1B3AABAB2A2B3AA3B1AB2A3B3这时两因素的三水平相互各碰一次

9、，它反映的情况全面，现在将因素C考虑进去，也同样希望在任何两个因素的不同水平之间各相碰一次而有不增加试验的次数，可做如下按排.OCB、AB1B2B3AAABCAB1C2AB2C2AB2c3ABC3A2RCAAB2cA3B3C2按上表安排的9次试验与简单比较法相比，试验次数相同但却克服了简单比较法的不均衡性，A的每个水平和B、C的三个水平分别各碰一次，B的每个水平和A、C的三个水平分别各碰一次，对C也是类似的情况。即三因素中任何两因素的不同水平均相碰一次因而试验是均衡的，上述9次试验可视为三因素三水平的全面试验的代表。为了书写方便，上述试验设计可简化为下表：CB123A112322313312表

10、中右下角部分的每一行和每一列中，1,2,3正好各出现一次，我们把具有这样的性质方块叫拉丁方,在排这种方块时常用拉丁字母，故有拉丁方之称。3正交设计法的基本特征3. 1均衡分散性在正交设计的试验安排中，各因素之间的搭配是均匀的，这种因素间搭配的均匀性一一试验点分布的均衡性成为正交设计的均衡分散性。或者说，正交试验设计把各试验条件均衡地分散在排列完全的水平组合之中，是之更具有代表性，更易于通过最少的试验次数来寻求最佳的试验条件，正交设计的这种性质，可以从试验结果的平均值中消除由于非均衡所引起的误差，有利于提高测定结果的可靠信。整齐可比性正交试验设计中，各因素各水平之间不仅搭配均匀，而且变化很有规律

11、。在考虑某因素的每一水平的试验中，其他各因素各水平出现的次数都相同，所作的贡献也认为是一致的。这样在比较各因素的每一水平对指标生产的影响时，就能最大限度地排除其他因素的干扰，突出本因素的作用，也就将各因素的效应清楚地加以区别并估计其大小，这就是正交试验设计的整齐可比性。在数学上把均衡分散性和整齐可比性称为正交性，凡具有这特性的试验设计方法都称为正交设计法。正是由于正交试验设计最大限度地排除了其他因素的干扰并消除了非均匀分散性可能造成的误差，因而只要比较因素各水平的试验指标的平均植，就能估计各因素对试验指标的影响大小，这在后面将作具体的介绍。两拉丁方的叠合在上述三因素三水平的基础上，如果还需同时

12、考虑第四个因素D,且因素D也取三个水平（D1,D2,D3）,那么能否在不增加试验次数而又能保持前述的要求呢？这首先应将D的三个水平拼成拉丁方，其次D的拉丁方和C的拉丁方不一样。对于前着，是使D也能与AB均衡搭配；对于后者，是使D与C之间也能均衡，既无重复，又无遗漏。若用（1）,（2）,（3）表示D的三个水平，而D的拉丁方与C的拉丁方相同时，其9次试验安排为：C(D)A-BX12311(1)2(2)3(3)22(2)3(3)1(1)33(3)1(1)2(2)这时A、B和D间是均衡的搭配，但C和D的搭配却不均衡，C的（1）水平和D的（1）水平相碰三次而不与D的（2）、（3）水平相碰，C的其他水平也

13、有类似的情况。所以上述的试验安排是不妥的，当试验的结果表明C的（1）水平最好，而在C取（1）水平时总是伴随着D的（1）水平的出现，自然也可以认为是D的（1）水平也最好，导致C和D的作用混杂。改进上述试验设计时，只需使D的拉丁方和C的拉丁方不同，两拉丁方具有均匀的搭配。按此原则可作如下的设计:C、DBA12311(1)21322(3)3(1)1(2)33(2)1(3)2(1)这时D的三个水平组成的是拉丁方，它和A、B及C之间的搭配都是均衡的，D的每一水平和C的1、2、3水平各碰一次，C的每一水平也和D的（1）、（2）、（3）水平各碰一次，既无重复，也无遗漏。现将两个拉丁方叫正交拉丁方。两个拉丁方

14、叠合在一起，就获得上述的试验设计，习惯上把具有这种性质的1(1)2(2)3(3)2(3)3(1)1(2)3(2)1(3)2(1)正交拉方设计因其搭配均衡，在分析试验数据时可以把每个因素的作用剖析得十分清楚而不致混杂，同时还可简便地寻求到最优的测量条件，达到预期的效果。第一部分正交试验结果的直观分析1,正交表及其使用正交表它是一种预先编制好的表格，根据这种表可合理安排试验并对试验数据作出判断对于前述的三因素三水平试的设计安排，可米用L9(3)正交表来完成。L9(3)表见表1.表1L9(34)正交表水平试验号112341111121222313334212352231623127313283213

15、93321表L9(34)读作L934,符号L表示正交表，L右下角的数字“9”表示此表有9行，即需安排9个实验，括号内数字的指数“4”表示有4歹U,即最多能安排四个因素；括号内数字的底数“3”表示每个因素取三个水平。表头的列号是置放试验中的因素(因素常记为A、B、CD)，表中列号1、2、3、4是在不考虑交互作用时最多可置放四个因素(因素少于四时，可只用其中几列)，表的左侧为试验号，表内的1、2、3是因素在试验中应分别取的水平，故称作水平号。L9(34)正交表可解决四因素(或少于四因素)的三水平试验设计问题，是一种较为简单的正交表。当试验因素及所取水平数更多时，则应选择其它种类的正交表，如Lie(

16、45)、L27(313)、L25(56)、Lie(42X29)等，其中Lie(42X29)表示作16个试验，可安两个四水平的因素和9个二水平的因素。.正交表的选择选择正交表时可考虑以下几点：()根据试验目的确定要考查的因素，如对试验的变化规律有大致的了解，有把握判断出影响试验效果的主要因素，可少取些因素，也可多取些因素，总之不能将主要影响因素漏掉。()确定各因素的变化范围和水平数，每个因素的水平数可以相等，也可以不等，一般地说，重要因素或者特别希望详细考查的因素，其变化范围可宽些水平数可多些，其余的因素所取水平数则可少些。()根据试验者进行试验时一次能平行完成的试验次数而选择正交表。()选用正

17、交表除考虑因素水平及试验条件外，还应考虑对试验结果精度的要求。当对试验结果的精度要求高时，宜取试验次数多的正交表，试验费用贵或试验周期长的，可取试验次数少的正交表。当存在交互作用时，应选用具交互作用的正交表。，般情况下，若因素全为二水平时，可选用L4(23)、L8(27)、Lie(215)等正交表；因素全是三水平时，可选用L9(34)、Li8(2X37)、L27(318)等正交表；若因素全为四水平的，可选用Lie(45)正交表；因素全为五水平的则选用Lie(45)正交表。当因素取不同水平时，一方面可采用下面即将介绍的拟水平法，一方面可直接套用L8(4X*28)、Li2(3X28)、Lie(4X

18、2i2)、Li8(42X29)等混合水平正交表。在三水平实验种选Li8(2X厂)，其中2水平所在的列，不做安排。三水平因素可在其它7列选用。正交试验的工作程序及几点说明在选择所需要的正交表后，将已确定的因素放置在表的任意列上，并把每一列的i、2、3填入具体水平，即得出试验方案。今仍以前述三因素三水平的显色反应为例，其试验方案如下表所示。表：三因素三水平正交试验表水平123试验结果A(t/oCB(mol/L)C(t/min)试验号11(25)101(10)21(25)202(20)31(25)303(30)42(30)102(20)52(30)203(30)62(30)301(10)73(35)

19、103(30)83(35)201(10)93(35)302(20)表中每一横行表示一次试验及进行该试验时所取的条件，按上安排作完实验后并将所测结果填入最后一列内，至于试验结果的分析，将在以后再作讨论。上面的试验设计表未考虑因素之间的交互作用，故选用L9(34)正交表，三因素在表上所处的列可任意选择而且可将因素的次序进行交换。如在1、2、3列可依次排列A、B、C三因素，也可安排为AC、B三因素，在把因素及水平排入正交表后而获得一张试验设计表，这过程叫表头设计。L9(34)表所安排的9次试验，不一定按表上的试验号码排列，也可按抽签的方法来决定，这样处理是为了减少试验中由于先后掌握不匀所带来的影响，

20、但对有些试验，其次序却不宜随意变更。对于每个因素的水平并不一定总是由小到大(或由大到小)按顺序排列，一般采用随机化方法来处理，即对部分因素的水平作随机的排列。常用的正交表三因素二水平正交表试验号正交表为L4(23),表头设计为：1111212232124221七因素二水平正交表正交表为L8(27),表头设计为:列号试验号12345671111111121112222312211224122221152121212621221217221122182212112更多因素二水平的正交法612221221211721221122121821212221112921122212211102221111

21、2212112212121112222112121221四因素三水平正交表正交表为L9(34),表头设计在前已述及，当为三因素时，此三因素可在表头上占取任意三列，如三因素三水平在选用L9(34)时，表头设计可为：试验号列号123111121223133421352216232731283239331七因素三水平正交表试'111111111212222221313333331421122331522233111623311221731213231832321311933132122101133221211121133221213221132132123132214223121321523

22、123212163132312217321312321833212312*：若把二水平的列1排进L18(37)表中，便得到混合型L18(21X37)表。更多因素的三水平正交表可选用L27(313)、L36(313)正交表。五因素四水平正交表正交表为L16(45),表头设计为:列号12345试验号'-、一、.1111112122223133334144445212346221437234128243219313421032431113312412342131341423144231415432411644132更多因素的四水平，可选用L32(49)正交表。六因素五水平正交表正交表为L25

23、(56),表头设计为:列号123456实验号、111111121222223133333414444451555556212345722345182345129245123102512341131352412324135133352411434135215352413164142531742531418431425194425312045314221515432225215432353215424543215255543212,二列间交互作用正交表二列间指两因素之间（因为因素占歹U）交互作用正交表除能对因素的主效应进行考查外，有时还能简便地考查各因素之间的交互作用并给出交互效应的大小。所谓交互

24、作用，是指在某些试验中，不仅因素自身对实验结果产生影响，而且因素之间产生协同的影响这种协同作用叫交互作用。如考查氮肥（N）和磷肥（P）对豆类增产效果，可在四块土质情况基本相同的土地上做四个试验，试验中施肥情况及产量如表所示.表：氮肥,磷肥对豆类产量的影响试验号N量(n/kg)P量(n/kg)产量（n/kg）100200230215302225432275由表知，单施氮肥3kg增产豆类15kg；单施磷肥2kg增产豆类25kg；同时施加了3kg氮肥和2kg磷，豆类增产量不是把两种肥料单独使用时增产豆类量的加和，而是增产了75kg,说明两种肥料对豆类增产起了协同的效果，这种作用叫氮肥和磷肥的交互作用

25、，以NXP表示。对于其它的因素，则记作因素1X因素2,或AXRAXC等。二列间交互作用正交表试验设计时，要考虑各因素间有无交互作用，这既可从专业本身加以判断，也可对一定的试验方案下的实验数据经统计分析来加以确定。在常用正交表中，有的只能考查因素本身的效应，不能用以考查因素间的交互作用；有的则可以分析因素间的交互作用，很多正交表都附有相应的二列间的交互作用表。在作表头设计时，若不考虑因素间的交互作用，则因素置那一列上可任意选取，若因素间存在交互作用，则因素的置放要根据一定的规则，应利用有交互作用的表来设计表头。今以L8（27）正交表来安排具有二列间交互作用的试验工作时，可由表2对因素及交互列在表

26、头中所处的列号作出安排。表2:心（27）二列间交互作用表列号列号、1(A)2(B)3(A*B)4(C)5(A*C)6(B*C)71(A)(1)3254762(B)(2)167453(A*B)(3)76544(C)(4)1235(A*C)(5)326(B*C)(6)17(7)表2中最上一行和最左侧一列数字以及括号(呈对角线)内的数字是列号，其余数字均为交互作用的列号。对于三因素而言，先将因素置放在表的第1、2歹U,则A和B相交的位置上的数字为3。即A*B应置放在第3列上，再将因素c置放于第4歹U,则A和C相交位置上的数字是5,B和C相交位置上的数字是6,这样A和C及B和C的交互作用列应分别为第5

27、列和第6歹U。如果考查时还有第四个因素D,并将它置放于第6歹U,根据上表可得如下的表头设计。列号1234567因素ABC*DA*BCB*DA*CDB*CA*D这样的设计中，虽有B和CXDC与BXDD与BXC的混杂，但如果已知BC、D之间的交互作用很小。故不致影响试验结果的分析，仍可进引因素A、B、C及交互作用AXBAXC及AXD的考查。如果要对四个因素及其两两之间的交互作用都作全面的考查，不允许上述存在的几种混杂，故此时不能选用L8(27)表，而选用Li6(215)二列向的交互作用表，见表3。.15表3:L16(2)二列向的交互作用表列123456789101112131415列'号1

28、(1)325476981110131215142(2)16745101189141512133(3)7654111098151413124(4)123121314158910115(5)32131215149811106(6)1141512131011897(7)151413121110988(8)12345679(9)32547610(1016745)11(117654)12(12123)13(1332)14(141)15(15)这样，对于四因素的表头设计为:列号123456789101112131415因素ABAXBCAxcBxcDAXDBXDCXD表3中，D未置入第7歹U。原因是D置于7列后，AXD应置第6歹U,导致与BXC的混杂。对于五因素。二水平的试验，在同时考虑各因素之间的交互作用时，因五因素自身及它们之间的两两交互作用共有15项，仍可用L16(215)二列间交互作用表，其表头设计为:列号12345678因素ABA*BCA*CB*CD*ED列号9101112131415因素A*DB*DC*EC*DB*EA*EE如果考查一个四因素三水平的问题，在只考虑因素主效应时，选用L8(27)正交表，让因素顺序上列，L8(27)