ch1正交试验设计2幻灯片

2021/8/21第一章第一章 正交试验设计正交试验设计 目的与要求：目的与要求： 学会使用正交设计法优选影响实验的各个学会使用正交设计法优选影响实验的各个因素，减少实验次数，节约人力物力。因素，减少实验次数，节约人力物力。2021/8/22 正交试验设计（Orthogonal Design）是于二十世纪50年代初期，由日本质量管理专家田口玄一（Tachugi）博士提出的在多因素试验设计方法的基础上，进一步研究开发出来的一种试验设计技术。 正交试验设计法使用一种规范化的表格（正交表）进行试验设计，可以用较少的试验次数，取得较为准确、可靠的优选结论。正交试验设计主要可以完成： 确定出各因素对试验指标的影响规律，得知哪些因素的确定出各因素对试验指标的影响规律，得知哪些因素的影响是主要的、哪些因素的影响是次要的、哪些因素之间影响是主要的、哪些因素的影响是次要的、哪些因素之间存在相互影响；存在相互影响； 选出各因素的一个水平组合来确定最佳生产条件。选出各因素的一个水平组合来确定最佳生产条件。 正交试验设计的基础是正交表。正交试验设计的基础是正交表。 2021/8/23试验设计例n 为提高某化工产品的转化率，选择了三个有关因素进行条件试验，反应温度(A)，反应时间(B)，用碱量(C)，并确定了它们的试验范围：A：80-90B：90-150分钟C：5-7试验目的是搞清楚因子A、B、C对转化率有什么影响，哪些是主要的，哪些是次要的，从而确定最适生产条件，即温度、时间及用碱量各为多 少才能使转化率高。试制定试验方案。2021/8/24n这里，对因子A，在试验范围内选了三个水平；因子B和C也都取三个水平：A：Al80，A285，A3=90B：Bl90分，B2120分，B3=150分C：Cl5，C26%，C37%当然，在正交试验设计中，因子可以是定量的，也可以是定性的。而定量因子各水平间的距离可以相等，也可以不相等。这个三因子三水平的条件试验，通常有两种试验进行方法： 2021/8/25n()取三因子所有水平之间的组合，即AlBlC1，A1BlC2，A1B2C1， ，A3B3C3，共有33=27次试验。用图表示就是图1 立方体的27个节点。这种试验法叫做全面试验法全面试验法。 2021/8/26全面试验法的优缺点：n优点：对各因素于试验指标之间的关系剖析得比较清楚n缺点：n试验次数太多，费时、费事，当因素水平比较多时，试验无法完成。n不做重复试验无法估计误差。n无法区分因素的主次。n例如选六个因素，每个因素选五个水平时，全面试验的数目是56 15625次。2021/8/27n()简单对比法简单对比法n 变化一个因素而固定其他因素，如首先固定B、C于Bl、Cl，使A变化之： A1B1C1 A2 A3 (好结果) 如得出结果A3最好，则固定A于A3，C还是Cl，使B变化之： B1A3C1 B2 (好结果) B3得出结果以B2为最好，则固定B于B2，A于A3，使C变化之： C1A3B2C2 (好结果) C3试验结果以C2最好。于是就认为最好的工艺条件是A3B2C2。 2021/8/28简单比较法的优缺点：n优点：试验次数少n缺点：n （1）试验点不具代表性。考察的因素水平仅局限于局部区域，不能全面地反映因素的全面情况。n （2）无法分清因素的主次。n （3）如果不进行重复试验，试验误差就估计不出来，因此无法确定最佳分析条件的精度。n （4）无法利用数理统计方法对试验结果进行分析，提出展望好条件。2021/8/29n首先这种方法的选点代表性很差，如按上述方法进行试验，试验点完全分布在一个角上，而在一个很大的范围内没有选点。因此这种试验方法不全面，所选的工艺条件A3B2C2不一定是27个组合中最好的。其次，用这种方法比较条件好坏时，是把单个的试验数据拿来，进行数值上的简单比较，而试验数据中必然要包含着误差成分，所以单个数据的简单比较不能剔除误差的干扰，必然造成结论的不稳定。 2021/8/210n 正交试验的提出：n考虑兼顾全面试验法和简单比较法的优点，利用根据数学原理制作好的规格化表正交表来设计试验不失为一种上策。n如上例， 对应于A有Al、A2、A3三个平面，对应于B、C也各有三个平面，共九个平面。则这九个平面上的试验点都应当一样多，即对每个因子的每个水平都要同等看待。具体来说，每个平面上都有三行、三列，要求在每行、每列上的点一样多。这样，作出如图2所示的设计，试验点用 表示。我们看到，在9个平面中每个平面上都恰好有三个点而每个平面的每行每列都有一个点，而且只有一个点，总共九个点。这样的试验方案，试验点的分布很均匀，试验次数也不多。 2021/8/2111.1 正交试验设计的基本概念 正交试验设计是利用正交表来安排与正交试验设计是利用正交表来安排与分析多因素试验的一种设计方法。它是由试验分析多因素试验的一种设计方法。它是由试验因素的全部水平组合中，挑选部分有代表性的因素的全部水平组合中，挑选部分有代表性的水平组合进行试验的，通过对这部分试验结果水平组合进行试验的，通过对这部分试验结果的分析了解全面试验的情况，找出最优的水平的分析了解全面试验的情况，找出最优的水平组合。组合。 2021/8/212 例如，要考察增稠剂用量、例如，要考察增稠剂用量、pHpH值和杀菌温值和杀菌温度对豆奶稳定性的影响。每个因素设置度对豆奶稳定性的影响。每个因素设置3 3个水个水平进行试验平进行试验 。A A因素是增稠剂用量，设因素是增稠剂用量，设A A1 1、A A2 2、A A3 3 3 3个水平；个水平；B B因素是因素是pHpH值，设值，设B B1 1、B B2 2、B B3 3 3 3个个水平；水平；C C因素为杀菌温度，设因素为杀菌温度，设C C1 1、C C2 2、C C3 3 3 3个水个水平。这是一个平。这是一个3 3因素因素3 3水平的试验，各因素的水水平的试验，各因素的水平之间全部可能组合有平之间全部可能组合有2727种。种。2021/8/213 全面试验：可以分析各因素的效应全面试验：可以分析各因素的效应 ，交互，交互作用，也可选出最优水平组合。但全面试验包作用，也可选出最优水平组合。但全面试验包含的水平组合数较多，工作量大含的水平组合数较多，工作量大 ，在有些情，在有些情况下无法完成。况下无法完成。 若试验主要目的是寻求最优水平组合，则可若试验主要目的是寻求最优水平组合，则可利用正交表来设计安排试验。利用正交表来设计安排试验。 正交试验设计的基本特点是：用部分试验来正交试验设计的基本特点是：用部分试验来代替全面试验，通过对部分试验结果的分析，代替全面试验，通过对部分试验结果的分析，了解全面试验的情况。了解全面试验的情况。2021/8/214 如对于上述如对于上述3 3因素因素3 3水平试验，若不考水平试验，若不考虑交互作用，可利用正交表虑交互作用，可利用正交表L L9 9(3(34 4) )安排，试验安排，试验方案仅包含方案仅包含9 9个水平组合，就能反映试验方案个水平组合，就能反映试验方案包含包含2727个水平组合的全面试验的情况，找出最个水平组合的全面试验的情况，找出最佳的生产条件。佳的生产条件。2021/8/2151.21.2正交试验设计的基本原理正交试验设计的基本原理 在试验安排中在试验安排中 ，每个因素在研究的范围，每个因素在研究的范围内选几个水平，就好比在选优区内打上网格内选几个水平，就好比在选优区内打上网格 ，如果网上的每个点都做试验，就是全面试验。如果网上的每个点都做试验，就是全面试验。如上例中，如上例中，3 3个因素的选优区可以用一个立方个因素的选优区可以用一个立方体表示（图体表示（图1-11-1），），3 3个因素各取个因素各取 3 3个水平，把个水平，把立方体划分成立方体划分成2727个格点，反映在个格点，反映在 图图1-11-1上就是上就是立方体内的立方体内的2727个个“.”.”。若。若2727个网格点都试验，个网格点都试验，就是全面试验，其试验方案如表就是全面试验，其试验方案如表1-11-1所示所示。2021/8/2162021/8/217 3 因因 素素 3 水水 平平 的的 全全 面试验水平组合面试验水平组合数为数为33=27，4 因素因素3水平的全面试验水平组水平的全面试验水平组合数为合数为34=81 ，5因素因素3水平的全面试验水平水平的全面试验水平组合数为组合数为35=243，这在科学试验中是有可能，这在科学试验中是有可能做不到的。做不到的。2021/8/218n正交设计就是从选优区全面试验点（水平组合）中挑正交设计就是从选优区全面试验点（水平组合）中挑选出有代表性的部分试验点（水平组合）来进行试验。选出有代表性的部分试验点（水平组合）来进行试验。图图1-11-1中标有试验号的九个中标有试验号的九个“()”()”，就是利用正交，就是利用正交表表L L9 9(3(34 4) )从从2727个试验点中挑选出来的个试验点中挑选出来的9 9个试验点。即：个试验点。即：n(1)A(1)A1 1B B1 1C C1 1 (2)A (2)A2 2B B1 1C C2 2 (3)A (3)A3 3B B1 1C C3 3n(4)A(4)A1 1B B2 2C C2 2 (5)A (5)A2 2B B2 2C C3 3 (6)A (6)A3 3B B2 2C C1 1n(7)A(7)A1 1B B3 3C C3 3 (8)A (8)A2 2B B3 3C C1 1 (9)A (9)A3 3B B3 3C C2 22021/8/219 上述选择上述选择 ，保证了，保证了A A因素的每个水平与因素的每个水平与B B因素、因素、C C因素的各个水平在试验中各搭配一次。对于因素的各个水平在试验中各搭配一次。对于A A、B B、C C 3 3个因素来说，是在个因素来说，是在2727个全面试验点中选择个全面试验点中选择9 9个试验点，个试验点，仅是全面试验的三分之一。仅是全面试验的三分之一。 从图中可以看到，从图中可以看到，9 9个试验点在选优区中分布是个试验点在选优区中分布是均衡的，在立方体的每个平面上，都恰是均衡的，在立方体的每个平面上，都恰是3 3个试验点；个试验点；在立方体的每条线上也恰有一个试验点。在立方体的每条线上也恰有一个试验点。 9 9个试验点均衡地分布于立方体内，有很强的代个试验点均衡地分布于立方体内，有很强的代表性，能够比较全面地反映选优区内基本情况。表性，能够比较全面地反映选优区内基本情况。2021/8/2201.3 1.3 正交表及其基本性质正交表及其基本性质1.3.1 1.3.1 正交表正交表 由于正交设计安排试验和分析试验结果都要用正由于正交设计安排试验和分析试验结果都要用正交表，因此先对正交表作介绍。交表，因此先对正交表作介绍。正交表正交表 在设计安排正交试验时制作好的标准化的表格。 表表1-21-2是一正交表，记号为是一正交表，记号为L L8 8(2(27 7) )，其中，其中“L L”代代表正交表；数字表正交表；数字“8”8”表示有表示有8 8行，用此正交表安排试行，用此正交表安排试验包含验包含8 8个处理个处理( (水平组合水平组合) ) ；底数；底数“2” 2” 表示因素的表示因素的水平数，指数水平数，指数“7”7”表示有表示有7 7列，用这张正交表最多可列，用这张正交表最多可以安排以安排7 7个个2 2水平因素。水平因素。2021/8/2212021/8/222行（试验）号列号1231111221231224221L4（23）正交表）正交表L4 （23）正交表代号正交表行数正交表字码总数（可安排的水平数）正交表列数（可安排的因素数）常用的正交表已由数学工作者制定出来，供进行正交设常用的正交表已由数学工作者制定出来，供进行正交设计时选用。计时选用。2水平正交表除水平正交表除L8(27)外，还有外，还有L4(23)、L16(215)等；等；3水平正交表有水平正交表有L9(34)、L27(213)等。等。2021/8/223 1.3.2 正交表的特性正交表的特性 任何一张正交表都有如下两个特性：任何一张正交表都有如下两个特性： 正交性正交性 （ 1）任一列中，不同数字出现的次数相等任一列中，不同数字出现的次数相等 例如例如L8(27)中不同数字只有中不同数字只有1和和2，它们各出现，它们各出现4次；次；L9(34)中不同数字有中不同数字有1、2和和3，它们各出现，它们各出现3次次 。2021/8/224 （ 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次。即每个因素次。即每个因素的一个水平与另一因素的各个水平互碰次数相等，表明任意两的一个水平与另一因素的各个水平互碰次数相等，表明任意两列各个数字之间的搭配是均匀的。列各个数字之间的搭配是均匀的。 2021/8/2252021/8/2262021/8/227 根据以上两个特性，我们用正交表安排根据以上两个特性，我们用正交表安排的试验，具有的试验，具有均衡分散均衡分散和和整齐可比整齐可比的特点。的特点。 所谓均衡分散，是指用正交表挑选出来所谓均衡分散，是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是的各因素水平组合在全部水平组合中的分布是均匀的均匀的 。 由由 图图1-1可以看出，在立方体中可以看出，在立方体中 ，任一平面内都包含，任一平面内都包含 3 个个“()”， 任一直线任一直线上都包含上都包含1个个“()” ，因此，因此 ，这些点代表性，这些点代表性强强 ，能够较好地反映全面试验的情况。，能够较好地反映全面试验的情况。 2021/8/228整齐可比是指每一个因素的各水平间具有可比性。因为正交表中每一因素的任一水平下都均衡地包含着另外因素的各个水平，当比较某因素不同水平时，其它因素的效应都彼此抵消。如在A、B、C 3个因素中，A因素的3个水平 A1、A2、A3 条件下各有 B 、C 的 3 个不同水平，即： 2021/8/229 在这在这9个水平组合中，个水平组合中，A因素各水平下包括了因素各水平下包括了B、C因素的因素的3个水平，虽然搭配方式不同，但个水平，虽然搭配方式不同，但B、C皆处于同等地位，当比较皆处于同等地位，当比较A因素不同水平时因素不同水平时，B因素不同水平的效应相互抵消，因素不同水平的效应相互抵消，C因素不同因素不同水平的效应也相互抵消。所以水平的效应也相互抵消。所以A因素因素3个水平间个水平间具有可比性。同样，具有可比性。同样，B、C因素因素3个水平间亦具个水平间亦具有可比性。有可比性。2021/8/2302021/8/2311.4 1.4 正交表的类别正交表的类别n1 1、等水平正交表、等水平正交表 各列水平数相各列水平数相同的正交表称为等水平正交表。如同的正交表称为等水平正交表。如L L4 4(2(23 3) )、L L8 8(2(27 7) )、L L1212(2(21111) )等各列中等各列中的水平为的水平为2 2，称为，称为2 2水平正交表；水平正交表；L L9 9(3(34 4) )、L L2727(3(31313) )等各列水平为等各列水平为3 3，称为称为3 3水平正交表。水平正交表。2021/8/2321.4 1.4 正交表的类别正交表的类别 2 2、混合水平正交表、混合水平正交表 各列水平数不完全各列水平数不完全相同的正交表称为混合水平正交表。如相同的正交表称为混合水平正交表。如L L8 8(4(42 24 4) )表中有一列的水平数为表中有一列的水平数为4 4，有，有4 4列水平数为列水平数为2 2。也就是说该表可以安排一。也就是说该表可以安排一个个4 4水平因素和水平因素和4 4个个2 2水平因素。再如水平因素。再如L L1616(4(44 42 23 3) )，L L1616(4(42 21212) )等都混合水平正等都混合水平正交表。交表。2021/8/2331.51.52021/8/2342021/8/235进行试验，记录试验结果进行试验，记录试验结果列方差分析表，列方差分析表，进行进行F F 检验检验计算各列偏差平方和、计算各列偏差平方和、自由度自由度2021/8/236 （1 1） 明确试验目的，确定试验指标明确试验目的，确定试验指标 试验设计前必须明确试验目的，即本次试试验设计前必须明确试验目的，即本次试验要解决什么问题。试验目的确定后，对试验验要解决什么问题。试验目的确定后，对试验结果如何衡量，即需要确定出结果如何衡量，即需要确定出试验指标试验指标。试验。试验指标可为指标可为定量指标定量指标，如强度、硬度、产量、出，如强度、硬度、产量、出品率、成本等；也可为品率、成本等；也可为定性指标定性指标如颜色、口感、如颜色、口感、光泽等。光泽等。2021/8/237 根据专业知识、以往的研究结论和经验，根据专业知识、以往的研究结论和经验，从影响试验指标的诸多因素中，通过因果分析从影响试验指标的诸多因素中，通过因果分析筛选出需要考察的试验因素。一般确定试验因筛选出需要考察的试验因素。一般确定试验因素时，应以对试验指标影响大的因素、尚未考素时，应以对试验指标影响大的因素、尚未考察过的因素、尚未完全掌握其规律的因素为先。察过的因素、尚未完全掌握其规律的因素为先。试验因素选定后，确定每个因素的水平，一般试验因素选定后，确定每个因素的水平，一般以以2 24 4个水平为宜。对主要考察的试验因素可个水平为宜。对主要考察的试验因素可多取水平，但不宜过多（多取水平，但不宜过多（66），否则试验次），否则试验次数骤增。数骤增。2021/8/238 实例：为提高山楂原料的利用率，研究酶法实例：为提高山楂原料的利用率，研究酶法液化工艺制造山楂原汁，拟通过正交试验来寻液化工艺制造山楂原汁，拟通过正交试验来寻找酶法液化的最佳工艺条件。找酶法液化的最佳工艺条件。 对本试验而言，试验目的是为了提高山楂对本试验而言，试验目的是为了提高山楂原料的利用率。所以可以以液化率原料的利用率。所以可以以液化率 液化率液化率=(=(果肉重量果肉重量- -液化后残渣重量液化后残渣重量)/)/果肉重果肉重量量 100%100%为试验指标，来评价液化工艺条件为试验指标，来评价液化工艺条件的好坏。液化率越高，山楂原料利用率就越高。的好坏。液化率越高，山楂原料利用率就越高。2021/8/239对本试验分析，影响山楂液化率的因素很对本试验分析，影响山楂液化率的因素很多，如山楂品种、山楂果肉的破碎度、果肉加多，如山楂品种、山楂果肉的破碎度、果肉加水量、原料水量、原料pH pH 值、果胶酶种类、加酶量、酶值、果胶酶种类、加酶量、酶解温度、酶解时间等等。经全面考虑，最后确解温度、酶解时间等等。经全面考虑，最后确定定果肉加水量、加酶量、酶解温度和酶解时间果肉加水量、加酶量、酶解温度和酶解时间为本试验的试验因素，分别记作为本试验的试验因素，分别记作A A、B B、C C和和D D，进行进行四因素正交试验，各因素均取三个水平四因素正交试验，各因素均取三个水平，因素水平表见表因素水平表见表1-31-3所示。所示。2021/8/2401-3 1-3 因素水平表因素水平表水平水平试验因素试验因素加水量加水量（mL/100gmL/100g）A A加酶量加酶量（mL/100gmL/100g）B B酶解温度酶解温度（）C C酶解时间酶解时间（h h）D D1 110101 120201.51.52 250504 435352.52.53 390907 750503.53.52021/8/241 正交表的选择是正交试验设计的首要问题。正交表的选择是正交试验设计的首要问题。确定因素及其水平后，根据因素、水平及需要确定因素及其水平后，根据因素、水平及需要考察的交互作用的多少来选择合适的正交表。考察的交互作用的多少来选择合适的正交表。正交表的选择原则是在能够安排下试验因素和正交表的选择原则是在能够安排下试验因素和交互作用的前提下，尽可能选用较小的正交表，交互作用的前提下，尽可能选用较小的正交表，以减少试验次数。以减少试验次数。 （3 3） 选择合适的正交表选择合适的正交表2021/8/242 一般情况下，试验因素的水平数应等于正一般情况下，试验因素的水平数应等于正 表中的水平数；因素个数（包括交互作用）应表中的水平数；因素个数（包括交互作用）应不大于正交表的列数；各因素及交互作用的自不大于正交表的列数；各因素及交互作用的自由度之和要小于所选正交表的总自由度，以便由度之和要小于所选正交表的总自由度，以便估计试验误差。若各因素及交互作用的自由度估计试验误差。若各因素及交互作用的自由度之和等于所选正交表总自由度，则可采用有重之和等于所选正交表总自由度，则可采用有重复正交试验来估计试验误差。复正交试验来估计试验误差。 2021/8/2432021/8/244列：列：正交表的列数正交表的列数cc因素所占列数因素所占列数+ +交互作用所占列数交互作用所占列数+ +空列。空列。自由度：自由度：正交表的总自由度（正交表的总自由度（a-1a-1）因素自由度因素自由度+ +交互作用自由度交互作用自由度+ +误差自误差自由度。由度。正交表选择依据：正交表选择依据：2021/8/245 此例有此例有4 4个个3 3水平因素，可以选用水平因素，可以选用L L9 9(3(34 4) )或或L L2727(3(31313) ) ；因本试验仅考察四个因素对液化；因本试验仅考察四个因素对液化率的影响效果，不考察因素间的交互作用，故率的影响效果，不考察因素间的交互作用，故宜选用宜选用L L9 9（3 34 4）正交表。若要考察交互作用，）正交表。若要考察交互作用，则应选用则应选用L L2727(3(31313) )。 （4 4）表头设计）表头设计所谓表头设计，就是把试验因素和要考察所谓表头设计，就是把试验因素和要考察的交互作用分别安排到正交表的各列中去的过的交互作用分别安排到正交表的各列中去的过程。程。2021/8/246在不考察交互作用时，各因素可随机安排在各在不考察交互作用时，各因素可随机安排在各列上；若考察交互作用，就应按所选正交表的交互列上；若考察交互作用，就应按所选正交表的交互作用列表安排各因素与交互作用，以防止设计作用列表安排各因素与交互作用，以防止设计“混混杂杂” ” 。 此例不考察交互作用，可将加水量此例不考察交互作用，可将加水量(A)(A)、加酶量、加酶量(B)(B)和酶解温度和酶解温度 (C)(C)、酶解时间（、酶解时间（D D）依次安排在）依次安排在L L9 9(3(34 4) )的第的第1 1、2 2、3 3、4 4列上，见表列上，见表1-41-4所示。所示。列号列号1234因素因素ABCD表表1-4 1-4 表头设计表头设计2021/8/247 把正交表中安排各因素的列（不包把正交表中安排各因素的列（不包含欲考察的交互作用列）中的每个水平含欲考察的交互作用列）中的每个水平数字换成该因素的实际水平值，便形成数字换成该因素的实际水平值，便形成了正交试验方案（表了正交试验方案（表1-51-5）。）。（5 5）编制试验方案，按方案进行试验，）编制试验方案，按方案进行试验，记录试验结果。记录试验结果。2021/8/248试验号试验号因因 素素ABCD111112122231333421235223162312731328321393321表表1-5 1-5 试验方案及试验结果试验方案及试验结果试验结果试验结果（液化率（液化率 %）01724124728118422021/8/249例1-2 鸭肉保鲜天然复合剂的筛选。试验以茶多酚作为天然复合保鲜剂的主要成分，分别添加不同增效剂、被膜剂和不同的浸泡时间，进行4因素4水平正交试验。试设计试验方案。（西南农业大学）n 有机酸和盐处理对鸭肉保鲜有明显的效果，但大部分属于合成的化学试剂，在卫生安全上得不到保证，并且不符合满足消费者纯天然、无污染的要求。2021/8/250n 明确目的，确定指标。明确目的，确定指标。本例目的是通过本例目的是通过试验，寻找最佳鸭肉天然复合保鲜剂。试验，寻找最佳鸭肉天然复合保鲜剂。n 选因素、定水平。选因素、定水平。根据专业知识和以前根据专业知识和以前研究结果，选择研究结果，选择4 4个因素，每个因素定个因素，每个因素定4 4个水个水平，因素水平表见表平，因素水平表见表1-61-6。n 选择正交表。选择正交表。此试验为此试验为4 4因素因素4 4水平试验，水平试验，不考虑交互作用，不考虑交互作用，4 4因素共占因素共占4 4列，选列，选L L1616（4 45 5）最合适，并有最合适，并有1 1空列，可以作为试验误差以空列，可以作为试验误差以衡量试验的可靠性。衡量试验的可靠性。n表头设计。表头设计。4 4因素任意放置。因素任意放置。n编制试验方案。编制试验方案。试验方案见表试验方案见表1-71-7。2021/8/251表表1-6 1-6 天然复合保鲜剂筛选试验因素水平表天然复合保鲜剂筛选试验因素水平表水平水平因因 素素A A茶多酚浓度茶多酚浓度/ /B B增效剂种类增效剂种类C C被膜剂种类被膜剂种类D D浸泡时间浸泡时间/min/min1 10.10.10.50.5维生素维生素C C0.50.5海藻酸海藻酸钠钠1 12 20.20.20.10.1柠檬酸柠檬酸0.80.8海藻酸海藻酸钠钠2 23 30.30.30.20.2-CD-CD1.01.0海藻酸海藻酸钠钠3 34 40.40.4生姜汁生姜汁1.01.0葡萄糖葡萄糖4 42021/8/252试验号试验号A A茶多酚浓度茶多酚浓度/ /B B增效剂种类增效剂种类C C被膜剂种类被膜剂种类D D浸泡时间浸泡时间/min/minE E 空列空列结果结果1 11 12 23 33 32 236.2036.202 22 24 41 12 22 231.5431.543 33 34 43 34 43 330.0930.094 44 42 21 11 13 329.3229.325 51 13 31 14 44 431.7731.776 62 21 13 31 14 435.0235.027 73 31 11 13 31 132.3732.378 84 43 33 32 21 132.6432.649 91 11 14 42 23 338.7938.7910102 23 32 23 33 330.9030.9011113 33 34 41 12 232.8732.8712124 41 12 24 42 234.5434.5413131 14 42 21 11 138.0238.0214142 22 24 44 41 135.6235.6215153 32 22 22 24 434.0234.0216164 44 44 43 34 432.8032.80表表1-7 天然复合保鲜剂筛选试验方案天然复合保鲜剂筛选试验方案2021/8/2531.5.2 1.5.2 试验结果分析试验结果分析n分清各因素及其交互作用的主次顺序，分清哪个是主要分清各因素及其交互作用的主次顺序，分清哪个是主要因素，哪个是次要因素；因素，哪个是次要因素；n判断因素对试验指标影响的显著程度；判断因素对试验指标影响的显著程度；n找出试验因素的优水平和试验范围内的最优组合，即试找出试验因素的优水平和试验范围内的最优组合，即试验因素各取什么水平时，试验指标最好；验因素各取什么水平时，试验指标最好；n分析因素与试验指标之间的关系，即当因素变化时，试分析因素与试验指标之间的关系，即当因素变化时，试验指标是如何变化的。找出指标随因素变化的规律和趋验指标是如何变化的。找出指标随因素变化的规律和趋势，为进一步试验指明方向；势，为进一步试验指明方向；n了解各因素之间的交互作用情况；了解各因素之间的交互作用情况；n估计试验误差的大小。估计试验误差的大小。2021/8/254直观分析法极差分析法直观分析法极差分析法Kjm为第为第j列因素列因素m水平所水平所对应的试验指标和，对应的试验指标和，kjm为为Kjm平均值。由平均值。由kjm大大小可以判断第小可以判断第j列因素优列因素优水平和优组合。水平和优组合。Rj为第为第j列因素的极差，反映了第列因素的极差，反映了第j列列因素水平波动时，试验指标的变动幅因素水平波动时，试验指标的变动幅度。度。Rj越大，说明该因素对试验指越大，说明该因素对试验指标的影响越大。根据标的影响越大。根据Rj大小，可以大小，可以判断因素的主次顺序。判断因素的主次顺序。2021/8/255（1 1）确定试验因素的最优水平组合）确定试验因素的最优水平组合（2 2）确定因素的主次顺序）确定因素的主次顺序（3 3）绘制因素与指标趋势图）绘制因素与指标趋势图 以各因素水平为横坐标，试验指标的以各因素水平为横坐标，试验指标的平均值为纵坐标，绘制因素与指标趋势图。平均值为纵坐标，绘制因素与指标趋势图。由因素与指标趋势图可以更直观地看出试由因素与指标趋势图可以更直观地看出试验指标随着因素水平的变化而变化的趋势，验指标随着因素水平的变化而变化的趋势，可为进一步试验指明方向。可为进一步试验指明方向。不考察交互作用的试验结果分析不考察交互作用的试验结果分析2021/8/256表表1-8 试验结果分析试验结果分析2021/8/257（2 2）计算各因素同一水平的平均值）计算各因素同一水平的平均值K Ki i。 K K1 1=36.20=36.20，K K2 2=33.27=33.27，K K3 3=32.34=32.34，K K4 4=31.83=31.831 1）计算）计算K Ki i值。值。K Ki i为同一水平之和。以第一列为同一水平之和。以第一列A A因素因素为例：为例：K K1 1=36.20+31.77+38.79+38.02=144.78=36.20+31.77+38.79+38.02=144.78 K K2 2=31.54+35.02+30.90+35.62=133.08=31.54+35.02+30.90+35.62=133.08 K K3 3=30.09+32.37+32.87+34.02=129.35=30.09+32.37+32.87+34.02=129.35 K K4 4=29.32+32.64+34.54+32.80=129.30=29.32+32.64+34.54+32.80=129.302021/8/258（3 3）计算各因素的极差）计算各因素的极差R R，R R表示该因素在其取值表示该因素在其取值范围内试验指标变化的幅度。范围内试验指标变化的幅度。 R R=max=max（K Ki i）-min-min（K Ki i）2021/8/2596 6）选优组合，即根据各因素各水平的平均）选优组合，即根据各因素各水平的平均值确定优水平，进而选出优组合。值确定优水平，进而选出优组合。 本例本例A A、B B、C C为主要因素，按照平均值大小为主要因素，按照平均值大小选取优水平为选取优水平为A A1 1B B1 1C C4 4，即茶多酚用量取，即茶多酚用量取0.1%0.1%水水平；以平；以0.5%0.5%维生素维生素C C作为增效剂；作为增效剂；1.0%1.0%葡萄糖葡萄糖液为被膜剂为形成的鸭肉保鲜复合剂为优组液为被膜剂为形成的鸭肉保鲜复合剂为优组合，而浸泡时间为次要因素，选取操作时间合，而浸泡时间为次要因素，选取操作时间1-3min1-3min即可。即可。2021/8/260表表1-9 鸭肉保鲜天然复合剂筛选试验结果鸭肉保鲜天然复合剂筛选试验结果2021/8/261思考题n 为提高某化工产品的转化率，选择了三个有关因素进行条件试验，反应温度(A)，反应时间(B)，用碱量(C)，并确定了它们的试验范围：A：80-90B：90-150分钟C：5-7试验目的是搞清楚因子A、B、C对转化率有什么影响，哪些是主要的，哪些是次要的，从而确定最适生产条件，即温度、时间及用碱量各为多 少才能使转化率高。试制定试验方案。2021/8/262n这里，对因子A，在试验范围内选了三个水平；因子B和C也都取三个水平：A：Al80，A285，A3=90B：Bl90分，B2120分，B3=150分C：Cl5，C26%，C37%这个三因子三水平的条件试验，2021/8/2632021/8/2642021/8/265n由此分别得出结论：温度越高转化率越好，以90为最好，但可以进一步探索温度更好的情况。反应时间以120分转化率最高。用碱量以6转化率最高。所以最适水平是A3B2C2。2021/8/2661.6 正交试验设计的方差分析正交试验设计的方差分析极差分析法简单明了，通俗易懂，计算工作量少便于推极差分析法简单明了，通俗易懂，计算工作量少便于推广普及。但这种方法广普及。但这种方法不能将试验中由于试验条件改变引不能将试验中由于试验条件改变引起的数据波动同试验误差引起的数据波动区分开来，起的数据波动同试验误差引起的数据波动区分开来，也也就是说，不能区分因素各水平间对应的试验结果的差异就是说，不能区分因素各水平间对应的试验结果的差异究竟是由于因素水平不同引起的，还是由于试验误差引究竟是由于因素水平不同引起的，还是由于试验误差引起的，起的，无法估计试验误差的大小无法估计试验误差的大小。此外，各因素对试验。此外，各因素对试验结果的影响大小无法给以精确的数量估计，不能提出一结果的影响大小无法给以精确的数量估计，不能提出一个标准来判断所考察因素作用是否显著。为了弥补极差个标准来判断所考察因素作用是否显著。为了弥补极差分析的缺陷，可采用方差分析。分析的缺陷，可采用方差分析。2021/8/267n方差分析基本思想是将数据的总变异分解成因方差分析基本思想是将数据的总变异分解成因素引起的变异和误差引起的变异两部分，构造素引起的变异和误差引起的变异两部分，构造F统计量，作统计量，作F检验，即可判断因素作用是否检验，即可判断因素作用是否显著。显著。2021/8/268表表1-10 L9(34)正交表正交表处理号处理号 第第1 1列（列（A A） 第第2 2列列 第第3 3列列 第第4 4列列 试验结果试验结果yiyi1 11 11 11 11 1y1y12 21 12 22 22 2y2y23 31 13 33 33 3y3y34 42 21 12 23 3y4y45 52 22 23 31 1y5y56 62 23 31 12 2y6y67 73 31 13 32 2y7y78 83 32 21 13 3y8y89 93 33 32 21 1y9y9因素因素重复重复1 1重复重复2 2重复重复3 3A A1 1y1y1y2y2y3y3A A2 2y4y4y5y5y6y6A A3 3y7y7y8y8y9y9和和y1+y2+y3y1+y2+y3K K1 1y4+y5+y6y4+y5+y6K K2 2y7+y8+y9y7+y8+y9K K3 39T-KKK31 9y9.y2y1y9y8y7y6y5y4y3y2y131SS22322212222A）（修正项）（）（）（）（2021/8/269表头设计表头设计A AB B试验数据试验数据列号列号1 12 2k kx xi ix xi i2 2试验号试验号1 11 1x x1 1x x1 12 22 21 1x x2 2x x2 22 2n nm mx xn nx xn n2 2K K1j1jK K1111K K1212K K1k1kK K2j2jK K2121K K2222K K2k2kK KmjmjK Km1m1K Km2m2K KmkmkK K1j1j2 2K K11112 2K K12122 2K K1k1k2 2K K2j2j2 2K K21212 2K K22222 2K K2k2k2 2K Kmjmj2 2K Km1m12 2K Km2m22 2K Kmkmk2 2SSSSj jSSSS1 1SSSS2 2SSSSk kCTQSSr1QCTQSSxQnT CTxTjjijTTniiT2niim12ij121K 表表1-11 Ln（mk）正交表及计算表格）正交表及计算表格2021/8/270nxxSSn1i2in1i2iT）（），（）（k.21j nxKr1SSn1i2im1i2ijj），（）（k.21j K-K n1SS22j1jj2021/8/271n-1dfT为因素水平个数，m j1mdf2021/8/272进行方差分析时，可制订如下的方差分析表：进行方差分析时，可制订如下的方差分析表： 总和总和 T误差误差 e（组间）（组间）因素因素 A（组内）（组内）显著性显著性结论结论F值值平均偏差平方平均偏差平方和和（方差）（方差）VS / f自由度自由度f偏差平方偏差平方和和S方差种类方差种类11 mf)1(2 kmf1 mkfTCTQS 1112SSST CTQSTT 111/ fSV 222/ fSV TTTfSV/ 221121/fSfSVVF 方差分析表方差分析表 2021/8/273不考虑交互作用等水平正交试验方差分析不考虑交互作用等水平正交试验方差分析 例：自溶酵母提取物是一种多用途食品配料。为探例：自溶酵母提取物是一种多用途食品配料。为探讨啤酒酵母的最适自溶条件，安排三因素三水平正讨啤酒酵母的最适自溶条件，安排三因素三水平正交试验。试验指标为自溶液中蛋白质含量（）。交试验。试验指标为自溶液中蛋白质含量（）。试验因素水平表见表试验因素水平表见表1-12，试验方案及结果分析，试验方案及结果分析见表见表1-13。试对试验结果进行方差分析。试对试验结果进行方差分析。水 平试验因素温度（）ApH值B加酶量（）C1506.52.02557.02.43587.52.8表表1-12 因素水平表因素水平表2021/8/274处理号处理号 A AB BC C空列空列试验结果试验结果yiyi1 11 1（5050）1 1（6.56.5）1 1（2.02.0）1 16.256.252 21 12 2（7.07.0）2 2（2.42.4）2 24.974.973 31 13 3（7.57.5）3 3（2.82.83 34.544.544 42 2（5555）1 12 23 37.537.535 52 22 23 31 15.545.546 62 23 31 12 25.55.57 73 3（5858）1 13 32 211.411.48 83 32 21 13 310.910.99 93 33 32 21 18.958.95K K1j1j15.76 15.76 25.18 25.18 22.65 22.65 20.74 20.74 K K2j2j18.57 18.57 21.41 21.41 21.45 21.45 21.87 21.87 K K3j3j31.25 31.25 18.99 18.99 21.48 21.48 22.97 22.97 K K1j1j2 2248.38 248.38 634.03 634.03 513.02 513.02 430.15 430.15 K K2j2j2 2344.84 344.84 458.39 458.39 460.10 460.10 478.30 478.30 K K3j3j2 2976.56 976.56 360.62 360.62 461.39 461.39 527.62 527.62 58.65T2021/8/275（1）计算）计算计算各列各水平的计算各列各水平的K值值 计算各列各水平对应数据之和计算各列各水平对应数据之和K1j、K2j、K3j及其平及其平方方K1j2、K2j2、K3j2。计算各列偏差平方和及自由度计算各列偏差平方和及自由度miijjCTKrss12186.477958.6522nTCT2221121311()31(248.38344.84976.56)477.8645.43AssKKKCT同理，同理，SSB=6.49，SSC=0.31 SSe=0.83（空列）（空列）2021/8/276自由度：自由度：dfAdfBdfCdfe3-1=2计算方差计算方差7 .2224 .45AAAdfSSV155. 0231. 0CCCdfSSV23. 3249. 6BBBdfSSV415. 0283. 0eeedfSSV（2）显著性检验）显著性检验根据以上计算，进行显著性检验，列出方差分析表，结果见表根据以上计算，进行显著性检验，列出方差分析表，结果见表1-14变异来源变异来源 平方和平方和 自由度自由度 均方均方 F F值值 Fa显著水平显著水平 A45.4045.402 222.7022.7079.6F0.05(2,4) =6.94*B6.496.492 23.243.2411.4F0.01(2,4)=18.0*C0.310.312 20.160.16误差误差e0.830.832 20.410.41误差误差e 1.141.144 40.2850.285总和总和 53.0353.03表表1-14 方差分析表方差分析表2021/8/277因素A高度显著，因素B显著，因素C不显著。因素主次顺序A-B-C。3）优化工艺条件的确定）优化工艺条件的确定本试验指标越大越好。对因素A、B分析，确定优水平为A3、B1；因素C的水平改变对试验结果几乎无影响，从经济角度考虑，选C1。优水平组合为A3B1C1。即温度为58，pH值为6.5，加酶量为2.0%。2021/8/278 自由度自由度 总的自由度总的自由度 1 nfT误差的自由度误差的自由度 )( BATeffff因素的自由度因素的自由度 11 mAfA的水平数的水平数11 mBfB的水平数的水平数 若存在交互作用（例如若存在交互作用（例如 ），那么交互作用的自），那么交互作用的自由度由度 ，则有，则有 BA BABAfff )(BABATefffff 2021/8/279同水平正交试验设计的方差分析同水平正交试验设计的方差分析 分析思路分析思路 ), 2 , 1(nixi 设正交试验共有设正交试验共有A、B、C、D、等实验因素（其中可等实验因素（其中可能包括交互作用），试验共进行能包括交互作用），试验共进行n次，因此有次，因此有n个试验结果，个试验结果，记为记为 。 求出总的偏差平方和求出总的偏差平方和 、各因素偏差平方和、各因素偏差平方和 以及误差偏差平方和以及误差偏差平方和 之后，根据它们各自的自由度之后，根据它们各自的自由度 、 以及以及 ，求出方差，求出方差 、 以及以及 ，然后计，然后计算相对于算相对于 的的F值值 ，再分别，再分别与各自的与各自的F临界值临界值 等进行比较，对各因素影响等进行比较，对各因素影响的显著性进行分析、判断。的显著性进行分析、判断。 TS,BASSeSTf,BAffefTV,BAVVeVeV,/,/BBBeAAfVFVVF ,BAFF 2021/8/280 计算计算 设试验因素为设试验因素为A、B、，总的试验次数为，总的试验次数为n，每个因，每个因素的水平数相同均为素的水平数相同均为m，水平的重复次数为，水平的重复次数为r，试验结果为，试验结果为 或表示成或表示成 （其中（其中 为水平序号，为水平序号， 为重复序号，为重复序号， 为总的试验次数）。为总的试验次数）。 ijx), 2 , 1(nixi mi, 2 , 1 rj, 2 , 1 mrn 总的偏差平方和总的偏差平方和ST CTQSTT 其中其中 mirjijTxQ112 所有所有n个数据的平方和个数据的平方和 nTCT2 校正系数校正系数 mirjijxT11总和总和 2021/8/281 各因素的偏差平方和（组内，各列）各因素的偏差平方和（组内，各列） SA、SB CTQSCTQSBBAA2111, mirjijBAxrQQ其中其中注意：注意： 计算计算 等时虽用同一公式，但括号中的等时虽用同一公式，但括号中的 依据依据不同的因素而不同（参见后面的示例不同的因素而不同（参见后面的示例 ）。）。,BAQQijx2021/8/282 试验误差的偏差平方和试验误差的偏差平方和 Se（与正交表中的空列相对应（与正交表中的空列相对应 ） eBATSSSS 一般根据关系一般根据关系 反求得到：反求得到： )( BATeSSSSeBABATSSSSS 注意：注意： 若试验存在因素的交互作用（例如若试验存在因素的交互作用（例如 ），则将这），则将这种交互作用视为一个因素，在数据计算表中单独占据一列种交互作用视为一个因素，在数据计算表中单独占据一列并计算并计算 （计算方法见后），则有（计算方法见后），则有BA BAS BAS 反求反求Se时也必须将时也必须将 考虑进去，估计显著性时也要对考虑进去，估计显著性时也要对其做相应的考虑。其做相应的考虑。 2021/8/283进行方差分析时，可制订如下的方差分析表：进行方差分析时，可制订如下的方差分析表： 总和总和 T误差误差 e（组间）（组间）因素因素 A（组内）（组内）显著性显著性结论结论F值值平均偏差平方平均偏差平方和和（方差）（方差）VS / f自由度自由度f偏差平方偏差平方和和S方差种类方差种类11 mf)1(2 kmf1 mkfTCTQS 1112SSST CTQSTT 111/ fSV 222/ fSV TTTfSV/ 221121/fSfSVVF 方差分析表方差分析表 2021/8/284 方差（平均偏差平方和）计算方差（平均偏差平方和）计算 组内偏差平方和的自由度组内偏差平方和的自由度 )1(2 kmf组间偏差平方和的自由度组间偏差平方和的自由度 11 mf总的偏差平方和的自由度总的偏差平方和的自由度 11 nmkfT三者关系：三者关系： 21fffT 因此因此 组内方差（组内平均偏差平方和）组内方差（组内平均偏差平方和） 222/ fSV 组间方差（组间平均偏差平方和）组间方差（组间平均偏差平方和） 111/ fSV 总总 方方 差（总的平均偏差平方和）差（总的平均偏差平方和） TTTfSV/ 2021/8/285 显著性检验显著性检验 可以证明，可以证明， 、 分别是自由度为分别是自由度为 、 的的 变量，变量， 是第一自由度为是第一自由度为 、第二自由度为、第二自由度为 的的F变量。根据上述基础可确定显著性检验的步骤如下：变量。根据上述基础可确定显著性检验的步骤如下： 1S2S1f2f2 221121/fSfSVVF 1f2f 计算试验的实际计算试验的实际F值：值： 221121/fSfSVVF 根据给定的显著性水平根据给定的显著性水平 及自由度及自由度 、 查查F分布表，确分布表，确定定F的临界值的临界值 。 1f2f),(21ffF （ 常取常取0.01、0.05、0.10、0.25） 将试验得到的实际将试验得到的实际F值与各种值与各种 下的临界值下的临界值 进行进行比较，得到因素水平变化对试验结果影响的显著性结论。比较，得到因素水平变化对试验结果影响的显著性结论。 ),(21ffF 2021/8/286F分布表（临界值分布表（临界值 表）表） F)25. 0( 2021/8/287F分布表（临界值分布表（临界值 表）表） F)10. 0( 2021/8/288F分布表（临界值分布表（临界值 表）表） F)05. 0( 2021/8/289F分布表（临界值分布表（临界值 表）表） F)01. 0( 部分资料从网络收集整理而来，供大家参考，感谢您的关注！