正交实验设计及结果分析

1、-. z.正交试验设计对于单因素或两因素试验，因其因素少，试验的设计、实施与分析都比拟简单。但在实际工作中，常常需要同时考察3个或3个以上的试验因素，假设进展全面试验，则试验的规模将很大，往往因试验条件的限制而难于实施。正交试验设计就是安排多因素试验、寻求最优水平组合的一种高效率试验设计方法。1 正交试验设计的概念及原理1.1 正交试验设计的根本概念正交试验设计是利用正交表来安排与分析多因素试验的一种设计方法。它是由试验因素的全部水平组合中，挑选局部有代表性的水平组合进展试验的，通过对这局部试验结果的分析了解全面试验的情况，找出最优的水平组合。例如：设计一个三因素、3水平的试验 A因素，设A1

2、、A2、A3 3个水平；B因素，设B1、B2、B3 3个水平；C因素，设C1、C2、C3 3个水平，各因素的水平之间全部可能组合有27种。全面试验：可以分析各因素的效应，交互作用，也可选出最优水平组合。但全面试验包含的水平组合数较多图示的27个节点，工作量大，在有些情况下无法完成。假设试验的主要目的是寻求最优水平组合，则可利用正交表来设计安排试验。全面试验法示意图三因素、三水平全面试验方案正交试验设计的根本特点是：用局部试验来代替全面试验，通过对局部试验结果的分析，了解全面试验的情况。正因为正交试验是用局部试验来代替全面试验的，它不可能像全面试验那样对各因素效应、交互作用一一分析；当交互作用存

3、在时，有可能出现交互作用的混杂。虽然正交试验设计有上述缺乏，但它能通过局部试验找到最优水平组合，因而很受实际工作者青睐。如对于上述3因素3水平试验，假设不考虑交互作用，可利用正交表L9(34)安排，试验方案仅包含9个水平组合，就能反映试验方案包含27个水平组合的全面试验的情况，找出最正确的生产条件。1.2 正交试验设计的根本原理正交设计就是从选优区全面试验点水平组合中挑选出有代表性的局部试验点水平组合来进展试验。上图中标有试验号的九个()，就是利用正交表L9(34)从27个试验点中挑选出来的9个试验点。即：(1)A1B1C1 (2)A2B1C2 (3)A3B1C3(4)A1B2C2 (5)A2

4、B2C3 (6)A3B2C1(7)A1B3C3 (8)A2B3C1 (9)A3B3C2以上选择，保证了A因素的每个水平与B因素、C因素的各个水平在试验中各搭配一次。对于A、B、C 3个因素来说，是在27个全面试验点中选择9个试验点，仅是全面试验的三分之一。从上图中可以看到，9个试验点在选优区中分布是均衡的，在立方体的每个平面上，都恰是3个试验点；在立方体的每条线上也恰有一个试验点。 9个试验点均衡地分布于整个立方体，有很强的代表性，能够比拟全面地反映选优区的根本情况。1.3 正交表及其根本性质1.3.1 正交表由于正交设计安排试验和分析试验结果都要用正交表，因此，我们先对正交表作一介绍。下表是

5、一正交表，记号为L8(27)，其中L代表正交表；L右下角的数字8表示有8行，用这正交表安排试验包含8个处理(水平组合) ；括号的底数2表示因素的水平数，括号2的指数7表示有7列，用这正交表最多可以安排7个2水平因素。L8(27) 正交表常用的正交表已由数学工作者制定出来，供进展正交设计时选用。2水平正交表除L8(27)外，还有L4(23)、L16(215)等；3水平正交表有L9(34)、L27(213)等。1.3.2 正交表的根本性质1.3.2.1 正交性1任一列中，各水平都出现，且出现的次数相等例：L8(27)中不同数字只有1和2，它们各出现4次；L9(34)中不同数字有1、2和3，它们各出

6、现3次。2任两列之间各种不同水平的所有可能组合都出现，且对出现的次数相等例：L8(27)中(1, 1), (1, 2), (2, 1), (2, 2)各出现两次；L9(34) 中 (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3)各出现1次。即每个因素的一个水平与另一因素的各个水平所有可能组合次数相等，说明任意两列各个数字之间的搭配是均匀的。1.3.2.2 代表性一方面：1任一列的各水平都出现，使得局部试验中包括了所有因素的所有水平；2任两列的所有水平组合都出现，使任意两因素间的试验组合为全面试验。另一方

7、面：由于正交表的正交性，正交试验的试验点必然均衡地分布在全面试验点中，具有很强的代表性。因此，局部试验寻找的最优条件与全面试验所找的最优条件，应有一致的趋势。1.3.2.3 综合可比性1任一列的各水平出现的次数相等；2任两列间所有水平组合出现次数相等，使得任一因素各水平的试验条件一样。这就保证了在每列因素各水平的效果中，最大限度地排除了其他因素的干扰。从而可以综合比拟该因素不同水平对试验指标的影响情况。根据以上特性，我们用正交表安排的试验，具有均衡分散和整齐可比的特点。所谓均衡分散，是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是均匀的。整齐可比是指每一个因素的各水平间具有可比性。

8、因为正交表中每一因素的任一水平下衡地包含着另外因素的各个水平，当比拟*因素不同水平时，其它因素的效应都彼此抵消。如在A、B、C 3个因素中，A因素的3个水平 A1、A2、A3 条件下各有 B 、C 的 3个不同水平，即：在这9个水平组合中，A因素各水平下包括了B、C因素的3个水平，虽然搭配方式不同，但B、C皆处于同等地位，当比拟A因素不同水平时，B因素不同水平的效应相互抵消，C因素不同水平的效应也相互抵消。所以A因素3个水平间具有综合可比性。同样，B、C因素3个水平间亦具有综合可比性。正交表的三个根本性质中，正交性是核心，是根底，代表性和综合可比性是正交性的必然结果。1.4 正交表的类别 1、

9、等水平正交表各列水平数一样的正交表称为等水平正交表。如L4(23)、L8(27)、L12(211)等各列中的水平为2，称为2水平正交表；L9(34)、L27(313)等各列水平为3，称为3水平正交表。 2、混合水平正交表各列水平数不完全一样的正交表称为混合水平正交表。如L8(424)表中有一列的水平数为4，有4列水平数为2。也就是说该表可以安排一个4水平因素和4个2水平因素。再如L16(4423)，L16(4212)等都混合水平正交表。试验目的与要求试验指标选因素、定水平试验目的与要求试验指标选因素、定水平因素、水平确定选择适宜正交表列试验方案对于多因素试验，正交试验设计是简单常用的一种试验设

10、计方法，其设计根本程序如下图。正交试验设计的根本程序包括试验方案设计及试验结果分析两局部。 2.1 试验方案设计1明确试验目的，确定试验指标试验设计前必须明确试验目的，即本次试验要解决什么问题。试验目确实定后，对试验结果如何衡量，即需要确定出试验指标。试验指标可为定量指标，也可为定性指标。试验方案设计流程一般为了便于试验结果的分析，定性指标可按相关的标准打分或模糊数学处理进展数量化，将定性指标定量化。2选因素、定水平，列因素水平表根据专业知识、以往的研究结论和经历，从影响试验指标的诸多因素中，通过因果分析筛选出需要考察的试验因素。一般确定试验因素时，应以对试验指标影响大的因素、尚未考察过的因素

11、、尚未完全掌握其规律的因素为先。试验因素选定后，根据所掌握的信息资料和相关知识，确定每个因素的水平，一般以2-4个水平为宜。对主要考察的试验因素，可以多取水平，但不宜过多6，否则试验次数骤增。因素的水平间距，应根据专业知识和已有的资料，尽可能把水平值取在理想区域。四因素、三水平的试验因素水平表水平试验因素ABCD1233选择适宜的正交表正交表的选择是正交试验设计的首要问题。确定了因素及其水平后，根据因素、水平及需要考察的交互作用的多少来选择适宜的正交表。正交表的选择原则是在能够安排下试验因素和交互作用的前提下，尽可能选用较小的正交表，以减少试验次数。一般情况下，试验因素的水平数应等于正交表中的

12、水平数；因素个数包括交互作用应不大于正交表的列数；最低的试验次数(行数)(每列水平数一1)+l 等水平正交表 Labc试验结果分析试验结果分析例：选择一4个3水平因素试验的正交表可以选用L9(34)或L27(313) A不考察因素间的交互作用，宜选用L934。B考察交互作用，则应选用L27(313)。课堂练习：选择一5个3水平因子及一个2水平因子试验的正交表 L12(235)4表头设计表头设计，就是把试验因素和要考察的交互作用分别安排到正交表的各列中去的过程。在不考察交互作用时，各因素可随机安排在各列上；假设考察交互作用，就应按所选正交表的交互作用列表安排各因素与交互作用，以防止设计混杂。例：

13、不考察交互作用，可将因素(A)、(B)和(C)、D依次安排在L9(34)的第1、2、3、4列上，见下表所示。表头设计列号1234因素ABCD5编制试验方案，按方案进展试验，记录试验结果。把正交表中安排各因素的列不包含欲考察的交互作用列中的每个水平数字换成该因素的实际水平值，便形成了下表中的正交试验方案。下表说明：试验号并非试验顺序，为了排除误差干扰，试验中可随机进展；安排试验方案时，局部因素的水平可采用随机安排。试验方案及试验结果表试验号因素试验结果ABCD111112122231333421235223162312731328321393321作业1、正交表有哪些类型？它们的核心性质是什么？

14、2、写出正交表的表达式，并简述正交试验设计的根本程序。3、不考虑交互作用，设计一个4水平的3因素正交试验方案极差分析方差分析2.2 试验结果分析分清各因素及其交互作用的主次顺序，分清哪个是主要因素，哪个是次要因素；判断因素对试验指标影响的显著程度；找出试验因素的优水平和试验围的最优组合，即试验因素各取什么水平时，试验指标最好；分析因素与试验指标之间的关系，即当因素变化时，试验指标是如何变化的。找出指标随因素变化的规律和趋势，为进一步试验指明方向；了解各因素之间的交互作用情况；估计试验误差的大小。3 正交试验的结果分析3.1 直观分析法极差分析法计算简便，直观，简单易懂，是正交试验结果分析最常用

15、方法。以下说明极差分析过程。Rj为第j列因素的极差，反映了第j列因素水平波动时，试验指标的变动幅度。Rj越大，说明该因素对试验指标的影响越大。根据Rj大小，可以判断因素的主次顺序。Kjm为第j列因素m水平所对应的试验指标和，kjm为Kjm平均值。由kjm大小可以判断第j列因素优水平和优组合。3.1.1 不考察交互作用的试验结果分析1确定试验因素的优水平和最优水平组合分析A因素各水平对试验指标的影响。根据正交设计的特性，对A1、A2、A3来说，三组试验的试验条件是完全一样的综合可比性，可进展直接比拟。如果因素A对试验指标无影响时，则kA1、kA2、kA3应该相等，不相等时说明，A因素的水平变动对

16、试验结果有影响。根据kA1、kA2、kA3的大小可以判断A1、A2、A3对试验指标的影响大小。kA值愈接近要求值的水平是A因素的优水平。同理，可以计算并确定B、C、D因素的优水平。四个因素的优水平组合为试验的最优水平组合。例1：分析下表中温度、时间、加碱量对转化率影响试验中各条件的最优值和最正确的工艺条件。转化率试验数据表2确定因素的主次顺序根据极差Rj的大小，可以判断各因素对试验指标的影响主次。比拟各R值大小，R值愈大的表示因素对指标的影响大，因素越重要，R值愈小因素的影响较小。3绘制因素与指标趋势图以各因素水平为横坐标，试验指标的平均值kjm为纵坐标，绘制因素与指标趋势图。由因素与指标趋势

17、图可以更直观地看出试验指标随着因素水平的变化而变化的趋势，可为进一步试验指明方向。以上即为正交试验极差分析的根本程序与方法极差R：表示该因素在其取值围试验指标变化的幅度。 R = ma*Ki- minKi例2：根据转化率试验结果计算极差R，并分析影响转化率因素的主次顺序。解例：计算的k值和R值如下表：温度时间加碱量以上计算后分析得到下面的试验结论各条件的最优值：温度390，时间2120分钟，加碱量26%。最正确工艺条件是以上三个最优水平的组合。对转化率影响最大的因素是温度，其次是加碱量，时间的影响最小。附1 多指标正交试验极差分析对于多指标试验，方案设计和实施与单指标试验一样，不同在于每做一次

18、试验，都需要对考察指标一一测试，分别记录。试验结果分析时，也要对考察指标一一分析，然后综合评衡，确定出优条件。1试验方案设计确定试验指标*、Y、Z挑因素，选水平，列因素水平表见下表选正交表、设计表头、编制试验方案。下表为四因素三水平试验，不考虑交互作用，选L934安排试验。试验号因素试验结果ABCD*YZ1111121222313334212352231623127313283213933212试验结果分析计算各因素各水平下每种试验指标的数据和以及平均值，并计算极差R。根据极差大小列出各指标下的因素主次顺序。试验指标：主次顺序优化水平组合 * ACDB A3B3C1D2 Y CDAB A1B2

19、C1D1Z ADBC A2B2C2D3初选优化工艺条件：根据各指标不同水平平均值确定各因素的优化水平组合。综合平衡确定最优工艺条件。以上三指标单独分析出的优化条件不一致，必须根据因素的影响主次，综合考虑，确定最正确工艺条件。例：如下表数据，对于因素A，其对*影响大小排第一位，此时取A3；其对Z影响也排第一位，取A2；而其对Y影响排次要第三位，为次要因素，因此A可取A2或A3，但取A2时，Z比取A3减小了14.8%，而*增加了12%，且由Y指标看，取A2比A3的Y值高，故A因素取A2。同理可分析B取B2，C取C1，D取D3。优组合为A2B2C1D3.K1K2K3*23.621.819.4Y2.6

20、2.42.3Z3.42.73.1附2 混合型正交表试验设计与极差分析试验设计与结果分析同前，不同的是将极差R进展调整，用调整后的R进展比拟。 r r 为因素每个水平试验重复数 d 折算系数，与因素水平有关。3.1.2 考察交互作用的试验设计1交互作用在多因素试验中，不仅因素对指标有影响，而且因素之间的联合搭配也对指标产生影响。因素间的联合搭配对试验指标产生的影响作用称为交互作用。因素之间的交互作用总是存在的，这是客观存在的普遍现象，只不过交互作用的程度不同而异。一般地，当交互作用很小时，就认为因素间不存在交互作用。对于交互作用，设计时应引起高度重视。在试验设计中，表示A、B间的交互作用记作AB

21、，称为1级交互作用；表示因素A、B、C之间的交互作用记作ABC，称为2级交互作用；依此类推，还有3级、4级交互作用等。2交互作用的处理原则试验设计中，交互作用一律当作因素对待，这是处理交互作用问题的总原则。作为因素，各级交互作用都可以安排在能考察交互作用的正交表的相应列上，它们对试验指标的影响情况都可以分析清楚，而且计算非常简单。但交互作用又与因素不同，表现在：用于考察交互作用的列不影响试验方案及其实施；一个交互作用并不一定只占正交表的一列，而是占有m-1p列。表头设计时，交互作用所占列数与因素的水平m有关，与交互作用级数p有关。2水平因素的各级交互作用均占1列；对于3水平因素，一级交互作用占

22、两列，二级交互作用占四列，可见，m和p越大，交互作用所占列数越多。例：对一个25因素试验，表头设计时，如果考虑所有各级交互作用，则连同因素本身，总计应占列数为： C51 + C52 +C53 +C54 +C555+10+10+5+131，则此试验必选L3225正交表进展设计。一般对于多因素试验，在满足试验要求的条件下，要有选择地、合理地考察*些交互作用。综合考虑试验目的、专业知识、以往的经历及现有试验条件等多方面情况进展交互作用选择。一般原则是：忽略高级交互作用有选择地考察一级交互作用。通常只考察那些作用效果较明显的，或试验要求必须考察的。试验允许的条件下，试验因素尽量取2水平。3有交互作用的

23、试验表头设计表头设计时，各因素及其交互作用不能任意安排，必须严格按交互作用列表进展安排。这是有交互作用正交试验设计的一个重要特点，也是关键的一步。在表头设计中，为了防止混杂，那些主要因素，重点要考察的因素，涉及交互作用较多的因素，应该优先安排，次要因素，不涉及交互作用的因素后安排。所谓混杂，就是指在正交表的同列中，安排了两个或两个以上的因素或交互作用，这样，就无法区分同一列中这些不同因素或交互作用对试验指标的影响效果。4有交互作用的正交设计在实际研究中，有时试验因素之间存在交互作用。对于既考察因素主效应又考察因素间交互作用的正交设计，除表头设计和结果分析与前面介绍略有不同外，其它根本一样。1、

24、右表是一正交试验方案和结果表，请用极差分析法对试验结果进展分析，确定各因素对指标Y的影响大小和优水平，明确最优的工艺组合。3.2 正交试验结果的方差分析极差分析法简单明了，通俗易懂，计算工作量少。但这种方法不能将试验中由于试验条件改变引起的数据波动同试验误差引起的数据波动区分开来，也就是说，不能区分因素各水平间对应的试验结果的差异终究是由于因素水平不同引起的，还是由于试验误差引起的，无法估计试验误差的大小。此外，各因素对试验结果的影响大小无法给以准确的数量估计，不能提出一个标准来判断所考察因素作用是否显著。为了弥补极差分析的缺陷，可采用方差分析。3.2.1 正交试验结果的方差分析方差分析根本思

25、想是将数据的总变异分解成因素引起的变异和误差引起的变异两局部，构造F统计量，作F检验，即可判断因素作用是否显著。1偏差平方和分解：总偏差平方和各列因素偏差平方和+误差偏差平方和2自由度分解：假设计算出的F值F0Fa，则拒绝原假设，认为该因素或交互作用对试验结果有显著影响；假设F0Fa，则认为该因素或交互作用对试验结果无显著影响。6正交试验方差分析说明由于进展F检验时，要用误差偏差平方和SSe及其自由度dfe，因此，为进展方差分析，所选正交表应留出一定空列。当无空列时，应进展重复试验，以估计试验误差。误差自由度一般不应小于2，dfe很小，F检验灵敏度很低，有时即使因素对试验指标有影响，用F检验也

26、判断不出来。为了增大dfe，提高F检验的灵敏度，在进展显著性检验之前，先将各因素和交互作用的方差与误差方差比拟，假设MS因MS交2MSe，可将这些因素或交互作用的偏差平方和、自由度并入误差的偏差平方和、自由度，这样使误差的偏差平方和和自由度增大，提高了F检验的灵敏度。L9(34) 正交表处理号第1列A第2列第3列第4列试验结果yi11111y121222y231333y342123y452231y562312y673132y783213y893321y9分析第1列因素时，其它列暂不考虑，将其看做条件因素。单因素试验数据资料格式Lnmk正交表及计算表格表头设计AB试验数据列号12k*i*i2试验

27、号11*1*1221*2*22nm*n*n2K1jK11K12K1kK2jK21K22K2kKmjKm1Km2KmkK1j2K112K122K1k2K2j2K212K222K2k2Kmj2Km12Km22Kmk2总偏差平方和：列偏差平方和： 总偏差平方和：列偏差平方和： 试验总次数为n，每个因素水平数为m个，每个水平作r次重复 rn/m。SS1SS2SSk试验方案及结果分析表处理号ABC空列试验结果1 1 1 116.252 1 2 224.973 1 3 334.544 2 1 237.535 2 2 315.546 2 3 125.57 3 1 3211.48 3 2 1310.99 3

28、3 218.95K1j15.76 25.18 22.65 20.74 K2j18.57 21.41 21.45 21.87 K3j31.25 18.99 21.48 22.97 K1j2248.38 634.03 513.02 430.15 K2j2344.84 458.39 460.10 478.30 K3j2976.56 360.62 461.39 527.62 1计算计算各列各水平的K值计算各列各水平对应数据之和K1j、K2j、K3j及其平方K1j2、K2j2、K3j2。自由度：dfAdfBdfCdfe3-1=2计算方差：2显著性检验根据以上计算，进展显著性检验，列出方差分析表，结果见下

29、表方差分析表变异来源平方和自由度均方F值Fa显著水平A45.40222.7079.6F0.05(2,4) =6.94*B6.4923.2411.4F0.01(2,4)=18.0*C0.3120.16误差e0.8320.41误差e1.1440.285总和53.03结论因素A高度显著，因素B显著，因素C不显著。因素主次顺序A-B-C。3优化工艺条件确实定本试验指标越大越好。对因素A、B分析，确定优水平为A3、B1；因素C的水平改变对试验结果几乎无影响，可从经济、操作等角度考虑取值，选C1。优水平组合为A3B1C1。3.2.3 考虑交互作用正交试验方差分析例：用石墨炉原子吸收分光光度法测定锌锭中的铅

30、，为了提高测定灵敏度，希望吸光度越大越好，现研究影响吸光度的因素，确定最正确测定条件计算计算各列各水平对应数据之和K1j、K2j及K1j-K2j；计算各列偏差平方和及自由度见吸光度测定试验方案及结果分析表。2显著性检验根据方差分析表，因素B高度显著，因素A、C及交互作用AB、AC、BC均不显著。各因素对试验结果影响的主次顺序为：B、A、AC、C、AB、BC。吸光度测定试验方案及结果分析表试验号ABABCACBC空列吸光度111111112.42211122222.24312211222.66412222112.58521212122.36621221212.4722112212.79822121122.76K1j9.99.4210.2110.2310.2410.1210.19K2j10.3110.79109.989.9710.0910.02K1j-K2j-0.41-1.370.210.250.270.030.17SSj0.0210.2350.00550.00780.00