正交实验设计及结果分析

正交试验设计对于单因素或两因素试验，因其因素少，试验的设计、实施与分析都比较简单。但在实际工作中，常常需要同时考察3个或3个以上的试验因素，若进行全面试验，则试验的规模将很大，往往因试验条件的限制而难于实施。正交试验设计就是安排多因素试验、寻求最优水平组合的一种高效率试验设计方法。1正交试验设计的概念及原理11正交试验设计的基本概念正交试验设计是利用正交表来安排与分析多因素试验的一种设计方法。它是由试验因素的全部水平组合中，挑选部分有代表性的水平组合进行试验的，通过对这部分试验结果的分析了解全面试验的情况，找出最优的水平组合。例如设计一个三因素、3水平的试验A因素，设A1、A2、A33个水平；B因素，设B1、B2、B33个水平；C因素，设C1、C2、C33个水平，各因素的水平之间全部可能组合有27种。全面试验可以分析各因素的效应，交互作用，也可选出最优水平组合。但全面试验包含的水平组合数较多（图示的27个节点），工作量大，在有些情况下无法完成。若试验的主要目的是寻求最优水平组合，则可利用正交表来设计安排试验。全面试验法示意图三因素、三水平全面试验方案正交试验设计的基本特点是用部分试验来代替全面试验，通过对部分试验结果的分析，了解全面试验的情况。正因为正交试验是用部分试验来代替全面试验的，它不可能像全面试验那样对各因素效应、交互作用一一分析；当交互作用存在时，有可能出现交互作用的混杂。虽然正交试验设计有上述不足，但它能通过部分试验找到最优水平组合，因而很受实际工作者青睐。如对于上述3因素3水平试验，若不考虑交互作用，可利用正交表L934安排，试验方案仅包含9个水平组合，就能反映试验方案包含27个水平组合的全面试验的情况，找出最佳的生产条件。12正交试验设计的基本原理正交设计就是从选优区全面试验点（水平组合）中挑选出有代表性的部分试验点（水平组合）来进行试验。上图中标有试验号的九个“”，就是利用正交表L934从27个试验点中挑选出来的9个试验点。即1A1B1C12A2B1C23A3B1C34A1B2C25A2B2C36A3B2C17A1B3C38A2B3C19A3B3C2以上选择，保证了A因素的每个水平与B因素、C因素的各个水平在试验中各搭配一次。对于A、B、C3个因素来说，是在27个全面试验点中选择9个试验点，仅是全面试验的三分之一。从上图中可以看到，9个试验点在选优区中分布是均衡的，在立方体的每个平面上，都恰是3个试验点；在立方体的每条线上也恰有一个试验点。9个试验点均衡地分布于整个立方体内，有很强的代表性，能够比较全面地反映选优区内的基本情况。13正交表及其基本性质131正交表由于正交设计安排试验和分析试验结果都要用正交表，因此，我们先对正交表作一介绍。下表是一张正交表，记号为L827，其中“L”代表正交表；L右下角的数字“8”表示有8行，用这张正交表安排试验包含8个处理水平组合；括号内的底数“2”表示因素的水平数，括号内2的指数“7”表示有7列，用这张正交表最多可以安排7个2水平因素。L827正交表常用的正交表已由数学工作者制定出来，供进行正交设计时选用。2水平正交表除L827外，还有L423、L16215等；3水平正交表有L934、L27213等。132正交表的基本性质1321正交性（1）任一列中，各水平都出现，且出现的次数相等例L827中不同数字只有1和2，它们各出现4次；L934中不同数字有1、2和3，它们各出现3次。（2）任两列之间各种不同水平的所有可能组合都出现，且对出现的次数相等例L827中1,1,1,2,2,1,2,2各出现两次；L934中1,1,1,2,1,3,2,1,2,2,2,3,3,1,3,2,3,3各出现1次。即每个因素的一个水平与另一因素的各个水平所有可能组合次数相等，表明任意两列各个数字之间的搭配是均匀的。1322代表性一方面（1）任一列的各水平都出现，使得部分试验中包括了所有因素的所有水平；（2）任两列的所有水平组合都出现，使任意两因素间的试验组合为全面试验。另一方面由于正交表的正交性，正交试验的试验点必然均衡地分布在全面试验点中，具有很强的代表性。因此，部分试验寻找的最优条件与全面试验所找的最优条件，应有一致的趋势。1323综合可比性（1）任一列的各水平出现的次数相等；（2）任两列间所有水平组合出现次数相等，使得任一因素各水平的试验条件相同。这就保证了在每列因素各水平的效果中，最大限度地排除了其他因素的干扰。从而可以综合比较该因素不同水平对试验指标的影响情况。根据以上特性，我们用正交表安排的试验，具有均衡分散和整齐可比的特点。所谓均衡分散，是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是均匀的。整齐可比是指每一个因素的各水平间具有可比性。因为正交表中每一因素的任一水平下都均衡地包含着另外因素的各个水平，当比较某因素不同水平时，其它因素的效应都彼此抵消。如在A、B、C3个因素中，A因素的3个水平A1、A2、A3条件下各有B、C的3个不同水平，即在这9个水平组合中，A因素各水平下包括了B、C因素的3个水平，虽然搭配方式不同，但B、C皆处于同等地位，当比较A因素不同水平时，B因素不同水平的效应相互抵消，C因素不同水平的效应也相互抵消。所以A因素3个水平间具有综合可比性。同样，B、C因素3个水平间亦具有综合可比性。试验目的与要求试验指标选因素、定水平因素、水平确定选择合适正交表列试验方案正交表的三个基本性质中，正交性是核心，是基础，代表性和综合可比性是正交性的必然结果。14正交表的类别1、等水平正交表各列水平数相同的正交表称为等水平正交表。如L423、L827、L12211等各列中的水平为2，称为2水平正交表；L934、L27313等各列水平为3，称为3水平正交表。2、混合水平正交表各列水平数不完全相同的正交表称为混合水平正交表。如L8424表中有一列的水平数为4，有4列水平数为2。也就是说该表可以安排一个4水平因素和4个2水平因素。再如L164423，L164212等都混合水平正交表。2正交试验设计的基本程序对于多因素试验，正交试验设计是简单常用的一种试验设计方法，其设计基本程序如图所示。正交试验设计的基本程序包括试验方案设计及试验结果分析两部分。21试验方案设计（1）明确试验目的，确定试验指标试验设计前必须明确试验目的，即本次试验要解决什么问题。试验目的确定后，对试验结果如何衡量，即需要确定出试验指标。试验指标可为定量指标，也可为定性指标。试验方案设计流程一般为了便于试验结果的分析，定性指标可按相关的标准打分或模糊数学处理进行数量化，将定性指标定量化。（2）选因素、定水平，列因素水平表根据专业知识、以往的研究结论和经验，从影响试验指标的诸多因素中，通过因果分析筛选出需要考察的试验因素。一般确定试验因素时，应以对试验指标影响大的因素、尚未考察过的因素、尚未完全掌握其规律的因素为先。试验因素选定后，根据所掌握的信息资料和相关知识，确定每个因素的水平，一般以24个水平为宜。对主要考察的试验因素，可以多取水平，但不宜过多（6），否则试验次数骤增。因素的水平间距，应根据专业知识和已有的资料，尽可能把水平值取在理想区域。四因素、三水平的试验因素水平表试验因素水平ABCD123（3）选择合适的正交表正交表的选择是正交试验设计的首要问题。确定了因素及其水平后，根据因素、水平及需要考察的交互作用的多少来选择合适的正交表。正交表的选择原则是在能够安排下试验因素和交互作用的前提下，尽可能选用较小的正交表，以减少试验次数。一般情况下，试验因素的水平数应等于正交表中的水平数；因素个数（包括交互作用）应不大于正交表的列数；最低的试验次数行数每列水平数一1L等水平正交表LA（BC）例选择一4个3水平因素试验的正交表可以选用L934或L27313（A）不考察因素间的交互作用，宜选用L9（34）。（B）考察交互作用，则应选用L27313。课堂练习试验结果分析选择一5个3水平因子及一个2水平因子试验的正交表L12235（4）表头设计表头设计，就是把试验因素和要考察的交互作用分别安排到正交表的各列中去的过程。在不考察交互作用时，各因素可随机安排在各列上；若考察交互作用，就应按所选正交表的交互作用列表安排各因素与交互作用，以防止设计“混杂”。例不考察交互作用，可将因素A、B和C、（D）依次安排在L934的第1、2、3、4列上，见下表所示。表头设计列号1234因素ABCD（5）编制试验方案，按方案进行试验，记录试验结果。把正交表中安排各因素的列（不包含欲考察的交互作用列）中的每个水平数字换成该因素的实际水平值，便形成了下表中的正交试验方案。下表说明试验号并非试验顺序，为了排除误差干扰，试验中可随机进行；安排试验方案时，部分因素的水平可采用随机安排。试验方案及试验结果表试验号因素试验结果ABCD111112122231333421235223162312731328321393321作业1、正交表有哪些类型它们的核心性质是什么2、写出正交表的表达式，并简述正交试验设计的基本程序。3、不考虑交互作用，设计一个4水平的3因素正交试验方案极差分析方差分析22试验结果分析分清各因素及其交互作用的主次顺序，分清哪个是主要因素，哪个是次要因素；判断因素对试验指标影响的显著程度；找出试验因素的优水平和试验范围内的最优组合，即试验因素各取什么水平时，试验指标最好；分析因素与试验指标之间的关系，即当因素变化时，试验指标是如何变化的。找出指标随因素变化的规律和趋势，为进一步试验指明方向；了解各因素之间的交互作用情况；估计试验误差的大小。3正交试验的结果分析31直观分析法极差分析法计算简便，直观，简单易懂，是正交试验结果分析最常用方法。以下说明极差分析过程。RJ为第J列因素的极差，反映了第J列因素水平波动时，试验指标的变动幅度。RJ越大，说明该因素对试验指标的影响越大。根据RJ大小，可以判断因素的主次顺序。KJM为第J列因素M水平所对应的试验指标和，KJM为KJM平均值。由KJM大小可以判断第J列因素优水平和优组合。311不考察交互作用的试验结果分析（1）确定试验因素的优水平和最优水平组合分析A因素各水平对试验指标的影响。根据正交设计的特性，对A1、A2、A3来说，三组试验的试验条件是完全一样的（综合可比性），可进行直接比较。如果因素A对试验指标无影响时，那么KA1、KA2、KA3应该相等，不相等时说明，A因素的水平变动对试验结果有影响。根据KA1、KA2、KA3的大小可以判断A1、A2、A3对试验指标的影响大小。KA值愈接近要求值的水平是A因素的优水平。同理，可以计算并确定B、C、D因素的优水平。四个因素的优水平组合为试验的最优水平组合。例1分析下表中温度、时间、加碱量对转化率影响试验中各条件的最优值和最佳的工艺条件。转化率试验数据表（2）确定因素的主次顺序根据极差RJ的大小，可以判断各因素对试验指标的影响主次。比较各R值大小，R值愈大的表示因素对指标的影响大，因素越重要，R值愈小因素的影响较小。（3）绘制因素与指标趋势图以各因素水平为横坐标，试验指标的平均值（KJM）为纵坐标，绘制因素与指标趋势图。由因素与指标趋势图可以更直观地看出试验指标随着因素水平的变化而变化的趋势，可为进一步试验指明方向。以上即为正交试验极差分析的基本程序与方法极差R表示该因素在其取值范围内试验指标变化的幅度。RMAX（KI）MIN（KI）例2根据转化率试验结果计算极差R，并分析影响转化率因素的主次顺序。解例计算的K值和R值如下表温度时间加碱量以上计算后分析得到下面的试验结论各条件的最优值温度3（90），时间2（120分钟），加碱量2（6）。最佳工艺条件是以上三个最优水平的组合。对转化率影响最大的因素是温度，其次是加碱量，时间的影响最小。附1多指标正交试验极差分析对于多指标试验，方案设计和实施与单指标试验相同，不同在于每做一次试验，都需要对考察指标一一测试，分别记录。试验结果分析时，也要对考察指标一一分析，然后综合评衡，确定出优条件。（1）试验方案设计确定试验指标X、Y、Z挑因素，选水平，列因素水平表（见下表）选正交表、设计表头、编制试验方案。下表为四因素三水平试验，不考虑交互作用，选L9（34）安排试验。因素试验结果试验号ABCDXYZ111112122231333421235223162312731328321393321ABCDXYZ1281005380215524821351200751701150225383713010275316023620302321212382830223322417222312193272833611318425303222190202733312072336K1709670602670K2655639664631K3581636679644K1236223201223K2218213221210K3194212226215极差R43112613K179749089K272756868K368706162K126253030K224252323K323232021极差R04021009K11029595103K280868790K393949382K134323234K227292930K331313127极差R07030307试验因素XY试验号Z123456789（2）试验结果分析计算各因素各水平下每种试验指标的数据和以及平均值，并计算极差R。根据极差大小列出各指标下的因素主次顺序。试验指标主次顺序优化水平组合XACDBA3B3C1D2YCDABA1B2C1D1ZADBCA2B2C2D3初选优化工艺条件根据各指标不同水平平均值确定各因素的优化水平组合。综合平衡确定最优工艺条件。以上三指标单独分析出的优化条件不一致，必须根据因素的影响主次，综合考虑，确定最佳工艺条件。例如下表数据，对于因素A，其对X影响大小排第一位，此时取A3；其对Z影响也排第一位，取A2；而其对Y影响排次要第三位，为次要因素，因此A可取A2或A3，但取A2时，Z比取A3减小了148，而X增加了12，且由Y指标看，取A2比A3的Y值高，故A因素取A2。同理可分析B取B2，C取C1，D取D3。优组合为A2B2C1D3K1K2K3X236218194Y262423Z342731附2混合型正交表试验设计与极差分析试验设计与结果分析同前，不同的是将极差R进行调整，用调整后的R进行比较。R为因素每个水平试验重复数DRD折算系数，与因素水平有关。水平数M2345678910折算系数D07105204504037035034032031312考察交互作用的试验设计（1）交互作用在多因素试验中，不仅因素对指标有影响，而且因素之间的联合搭配也对指标产生影响。因素间的联合搭配对试验指标产生的影响作用称为交互作用。因素之间的交互作用总是存在的，这是客观存在的普遍现象，只不过交互作用的程度不同而异。一般地，当交互作用很小时，就认为因素间不存在交互作用。对于交互作用，设计时应引起高度重视。在试验设计中，表示A、B间的交互作用记作AB，称为1级交互作用；表示因素A、B、C之间的交互作用记作ABC，称为2级交互作用；依此类推，还有3级、4级交互作用等。（2）交互作用的处理原则试验设计中，交互作用一律当作因素看待，这是处理交互作用问题的总原则。作为因素，各级交互作用都可以安排在能考察交互作用的正交表的相应列上，它们对试验指标的影响情况都可以分析清楚，而且计算非常简单。但交互作用又与因素不同，表现在用于考察交互作用的列不影响试验方案及其实施；一个交互作用并不一定只占正交表的一列，而是占有（M1）P列。表头设计时，交互作用所占列数与因素的水平M有关，与交互作用级数P有关。2水平因素的各级交互作用均占1列；对于3水平因素，一级交互作用占两列，二级交互作用占四列，可见，M和P越大，交互作用所占列数越多。例对一个25因素试验，表头设计时，如果考虑所有各级交互作用，那么连同因素本身，总计应占列数为C51C52C53C54C55510105131，那么此试验必选L32（25）正交表进行设计。一般对于多因素试验，在满足试验要求的条件下，要有选择地、合理地考察某些交互作用。综合考虑试验目的、专业知识、以往的经验及现有试验条件等多方面情况进行交互作用选择。一般原则是忽略高级交互作用有选择地考察一级交互作用。通常只考察那些作用效果较明显的，或试验要求必须考察的。试验允许的条件下，试验因素尽量取2水平。（3）有交互作用的试验表头设计表头设计时，各因素及其交互作用不能任意安排，必须严格按交互作用列表进行安排。这是有交互作用正交试验设计的一个重要特点，也是关键的一步。在表头设计中，为了避免混杂，那些主要因素，重点要考察的因素，涉及交互作用较多的因素，应该优先安排，次要因素，不涉及交互作用的因素后安排。所谓混杂，就是指在正交表的同列中，安排了两个或两个以上的因素或交互作用，这样，就无法区分同一列中这些不同因素或交互作用对试验指标的影响效果。（4）有交互作用的正交设计在实际研究中，有时试验因素之间存在交互作用。对于既考察因素主效应又考察因素间交互作用的正交设计，除表头设计和结果分析与前面介绍略有不同外，其它基本相同。1、右表是一正交试验方案和结果表，请用极差分析法对试验结果进行分析，确定各因素对指标Y的影响大小和优水平，明确最优的工艺组合。32正交试验结果的方差分析极差分析法简单明了，通俗易懂，计算工作量少。但这种方法不能将试验中由于试验条件改变引起的数据波动同试验误差引起的数据波动区分开来，也就是说，不能区分因素各水平间对应的试验结果的差异究竟是由于因素水平不同引起的，还是由于试验误差引起试验号ABC试验结果Y111111210212222208321122215422211230531212251632121247741221238842112230K1418914902921925K2445915927908904K3498K4468K120902285225523032313K222252288231822702260K32490K42340极差R40025625325525调整R2546036888462746主次顺序优组合ACBA3B2C2或A3B1C2的，无法估计试验误差的大小。此外，各因素对试验结果的影响大小无法给以精确的数量估计，不能提出一个标准来判断所考察因素作用是否显著。为了弥补极差分析的缺陷，可采用方差分析。321正交试验结果的方差分析方差分析基本思想是将数据的总变异分解成因素引起的变异和误差引起的变异两部分，构造F统计量，作F检验，即可判断因素作用是否显著。（1）偏差平方和分解总偏差平方和各列因素偏差平方和误差偏差平方和（2）自由度分解若计算出的F值F0FA，则拒绝原假设，认为该因素或交互作用对试验结果有显著影响；若F0FA，则认为该因素或交互作用对试验空列（误差）因素SSST空列（误差）因素TDFDF结果无显著影响。（6）正交试验方差分析说明由于进行F检验时，要用误差偏差平方和SSE及其自由度DFE，因此，为进行方差分析，所选正交表应留出一定空列。当无空列时，应进行重复试验，以估计试验误差。误差自由度一般不应小于2，DFE很小，F检验灵敏度很低，有时即使因素对试验指标有影响，用F检验也判断不出来。为了增大DFE，提高F检验的灵敏度，在进行显著性检验之前，先将各因素和交互作用的方差与误差方差比较，若MS因（MS交）2MSE，可将这些因素或交互作用的偏差平方和、自由度并入误差的偏差平方和、自由度，这样使误差的偏差平方和和自由度增大，提高了F检验的灵敏度。L934正交表处理号第1列（A）第2列第3列第4列试验结果YI11111Y121222Y231333Y342123Y452231Y562312Y673132Y783213Y893321Y9分析第1列因素时，其它列暂不考虑，将其看做条件因素。单因素试验数据资料格式LN（MK）正交表及计算表格表头设计AB试验数据列号12KXIXI2试验号9T）K（31（修正项）9）Y2Y1（）9Y87（）Y654（）Y2（S232222A11X1X1221X2X22NMXNXN2K1JK11K12K1KK2JK21K22K2KKMJKM1KM2KMKK1J2K112K122K1K2K2J2K212K222K2K2KMJ2KM12KM22KMK2SSJSS1SS2SSKCTQSR1CTQSXNTCTJIJTNIT2NIM2J121K总偏差平方和列偏差平方和），（）（K21JNXKR1SN1I2IMI2IJJ试验总次数为N，每个因素水平数为M个，每个水平作R次重复RN/M。试验方案及结果分析表处理号ABC空列试验结果11111625212224973133345442123753522315546231255731321148321310993321895K1J1576251822652074K2J1857214121452187K3J3125189921482297K1J224838634035130243015K2J234484458394601047830K3J297656360624613952762（1）计算计算各列各水平的K值计算各列各水平对应数据之和K1J、K2J、K3J及其平方K1J2、K2J2、K3J2。自由度DFADFBDFCDFE312计算方差（2）显著性检验根据以上计算，进行显著性检验，列出方差分析表，结果见下表方差分析表7245ADFSV232496BBDFSV150230CCDFSV41502830EEDFSV变异来源平方和自由度均方F值FA显著水平A454022270796F0052,4694B6492324114F0012,4180C0312016误差E0832041误差E11440285总和5303结论因素A高度显著，因素B显著，因素C不显著。因素主次顺序ABC。（3）优化工艺条件的确定本试验指标越大越好。对因素A、B分析，确定优水平为A3、B1；因素C的水平改变对试验结果几乎无影响，可从经济、操作等角度考虑取值，选C1。优水平组合为A3B1C1。323考虑交互作用正交试验方差分析例用石墨炉原子吸收分光光度法测定锌锭中的铅，为了提高测定灵敏度，希望吸光度越大越好，现研究影响吸光度的因素，确定最佳测定条件（1）计算计算各列各水平对应数据之和K1J、K2J及（K1JK2J）；计算各列偏差平方和及自由度（见吸光度测定试验方案及结果分析表）。（2）显著性检验根据方差分析表，因素B高度显著，因素A、C及交互作用AB、AC、BC均不显著。各因素对试验结果影响的主次顺序为B、A、AC、C、AB、BC。吸光度测定试验方案及结果分析表试验号ABABCACBC空列吸光度111111112422111222222431221122266412222112585212121223662122121247