正交实验设计四因素三水平

1、关于正交实验的设计四因素三水平第一张，PPT共一百二十三页，创作于2022年6月1.1 正交试验设计的基本概念 正交试验设计是利用正交表来安排与分析多因素试验的一种设计方法。它是由试验因素的全部水平组合中，挑选部分有代表性的水平组合进行试验的，通过对这部分试验结果的分析了解全面试验的情况，找出最优的水平组合。下一张 主 页 退 出 上一张 1 正交试验设计的概念及原理第二张，PPT共一百二十三页，创作于2022年6月 例如，要考察增稠剂用量、pH值和杀菌温度对豆奶稳定性的影响。每个因素设置3个水平进行试验 。 A因素是增稠剂用量，设A1、A2、A3 3个水平；B因素是pH值，设B1、B2、B3

2、 3个水平；C因素为杀菌温度，设C1、C2、C3 3个水平。这是一个3因素3水平的试验，各因素的水平之间全部可能组合有27种 。 全面试验：可以分析各因素的效应 ，交互作用，也可选出最优水平组合。但全面试验包含的水平组合数较多，工作量大 ，在有些情况下无法完成 。 若试验的主要目的是寻求最优水平组合，则 可利用正交表来设计安排试验。第三张，PPT共一百二十三页，创作于2022年6月 正交试验设计的基本特点是：用部分试验来代替全面试验，通过对部分试验结果的分析，了解全面试验的情况。 正因为正交试验是用部分试验来代替全面试验的，它不可能像全面试验那样对各因素效应、交互作用一一分析；当交互作用存在时

3、，有可能出现交互作用的混杂。虽然正交试验设计有上述不足，但它能通过部分试验找到最优水平组合 ，因 而 很 受实际工作者青睐。 下一张 主 页 退 出 上一张 第四张，PPT共一百二十三页，创作于2022年6月 如对于上述3因素3水平试验，若不考虑交互作用，可利用正交表L9(34)安排，试验方案仅包含9个水平组合，就能反映试验方案包含27个水平组合的全面试验的情况，找出最佳的生产条件。1.2 正交试验设计的基本原理 在试验安排中 ，每个因素在研究的范围内选几个水平，就好比在选优区内打上网格 ，如果网上的每个点都做试验，就是全面试验。如上例中，3个因素的选优区可以用一个立方体表示（图10-1），3

4、个因素各取 3个水平，把立方体划分成27个格点，反映在 图10-1上就是立方体内的27个“.”。若27个网格点都试验，就是全面试验，其试验方案如表10-1所示。 第五张，PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 表10-1第六张，PPT共一百二十三页，创作于2022年6月 3 因 素 3 水 平 的 全 面试验水平组合数为33=27，4 因素3水平的全面试验水平组合数为34=81 ，5因素3水平的全面试验水平组合数为35=243，这在科学试验中是有可能做不到的。 下一张 主 页 退 出 上一张 图10-1第七张，PPT共一百二十三页，创作于2022年6月 正交设

5、计就是从选优区全面试验点（水平组合）中挑选出有代表性的部分试验点（水平组合）来进行试验。图10-1中标有试验号的九个“()”，就是利用正交表L9(34)从27个试验点中挑选出来的9个试验点。即：(1)A1B1C1 (2)A2B1C2 (3)A3B1C3(4)A1B2C2 (5)A2B2C3 (6)A3B2C1(7)A1B3C3 (8)A2B3C1 (9)A3B3C2下一张 主 页 退 出 上一张 第八张，PPT共一百二十三页，创作于2022年6月 上述选择 ，保证了A因素的每个水平与B因素、C因素的各个水平在试验中各搭配一次 。对于A、B、C 3个因素来说 ， 是在27个全面试验点中选择9个试

6、验点 ，仅 是全面试验的 三分之一。 从图10-1中可以看到 ，9个试验点在选优区中分布是均衡的，在立方体的每个平面上 ，都恰是3个试验点；在立方体的每条线上也恰有一个试验点。 9个试验点均衡地分布于整个立方体内 ，有很强的代表性 ， 能 够比较全面地反映选优区内的基本情况。 下一张 主 页 退 出 上一张 第九张，PPT共一百二十三页，创作于2022年6月1.3 正交表及其基本性质1.3.1 正交表 由于正交设计安排试验和分析试验结果都要用正交表，因此，我们先对正交表作一介绍。 表10-2是一张正交表，记号为L8(27)，其中“L”代表正交表；L右下角的数字“8”表示有8行 ，用这张正交表安

7、排试验包含8个处理(水平组合) ；括号内的底数“2” 表示因素的水平数，括号内2的指数“7”表示有7列 ，用这张正交表最多可以安排7个2水平因素。 下一张 主 页 退 出 上一张 第十张，PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 表10-2第十一张，PPT共一百二十三页，创作于2022年6月 常用的正交表已由数学工作者制定出来，供进行正交设计时选用。2水平正交表除L8(27)外，还有L4(23)、L16(215)等；3水平正交表有L9(34)、L27(213)等（详见附表14及有关参考书）。1.3.2 正交表的基本性质 1.3.2.1 正交性 （1）任一列中，各

8、水平都出现，且出现的次数相等 例如L8(27)中不同数字只有1和2，它们各出现4次；L9(34)中不同数字有1、2和3，它们各出现3次 。下一张 主 页 退 出 上一张 第十二张，PPT共一百二十三页，创作于2022年6月（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次。即每个因素的一个水平与另一因素的各个水平所有可

9、能组合次数相等，表明任意两列各个数字之间的搭配是均匀的。 下一张 主 页 退 出 上一张 第十三张，PPT共一百二十三页，创作于2022年6月1.3.2.2 代表性 一方面： （1）任一列的各水平都出现，使得部分试验中包括了所有因素的所有水平； （2）任两列的所有水平组合都出现，使任意两因素间的试验组合为全面试验。另一方面：由于正交表的正交性，正交试验的试验点必然均衡地分布在全面试验点中，具有很强的代表性。因此，部分试验寻找的最优条件与全面试验所找的最优条件，应有一致的趋势。第十四张，PPT共一百二十三页，创作于2022年6月1.3.2.3 综合可比性 （1）任一列的各水平出现的次数相等；（2

10、）任两列间所有水平组合出现次数相等，使得任一因素各水平的试验条件相同。这就保证了在每列因素各水平的效果中，最大限度地排除了其他因素的干扰。从而可以综合比较该因素不同水平对试验指标的影响情况。第十五张，PPT共一百二十三页，创作于2022年6月 根据以上特性，我们用正交表安排的试验，具有均衡分散和整齐可比的特点。 所谓均衡分散，是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是均匀的 。 由 图10-1可以看出，在立方体中 ，任一平面内都包含 3 个“()”， 任一直线上都包含1个“()” ，因此 ，这些点代表性强 ，能够较好地反映全面试验的情况。 下一张 主 页 退 出 上一张 第十

11、六张，PPT共一百二十三页，创作于2022年6月 整齐可比是指每一个因素的各水平间具有可比性。因为正交表中每一因素的任一水平下都均衡地包含着另外因素的各个水平 ，当比较某因素不同水平时，其它 因素的效应都彼此抵消。如在A、B、C 3个因素中，A因素的3个水平 A1、A2、A3 条件下各有 B 、C 的 3个不同水平，即： 第十七张，PPT共一百二十三页，创作于2022年6月 在这9个水平组合中，A因素各水平下包括了B、C因素的3个水平，虽然搭配方式不同，但B、C皆处于同等地位，当比较A因素不同水平时，B因素不同水平的效应相互抵消，C因素不同水平的效应也相互抵消。所以A因素3个水平间具有综合可比

12、性。同样，B、C因素3个水平间亦具有综合可比性。下一张 主 页 退 出 上一张 第十八张，PPT共一百二十三页，创作于2022年6月正交表的三个基本性质中，正交性是核心，是基础，代表性和综合可比性是正交性的必然结果第十九张，PPT共一百二十三页，创作于2022年6月1.4 正交表的类别 1、等水平正交表 各列水平数相同的正交表称为等水平正交表。如L4(23)、L8(27)、L12(211)等各列中的水平为2，称为2水平正交表；L9(34)、L27(313)等各列水平为3，称为3水平正交表。 2、混合水平正交表 各列水平数不完全相同的正交表称为混合水平正交表。如L8(424)表中有一列的水平数为

13、4，有4列水平数为2。也就是说该表可以安排一个4水平因素和4个2水平因素。再如L16(4423)，L16(4212)等都混合水平正交表。下一张 主 页 退 出 上一张 第二十张，PPT共一百二十三页，创作于2022年6月2 正交试验设计的基本程序 对于多因素试验，正交试验设计是简单常用的一种试验设计方法，其设计基本程序如图所示。正交试验设计的基本程序包括试验方案设计及试验结果分析两部分。第二十一张，PPT共一百二十三页，创作于2022年6月试验目的与要求试验指标选因素、定水平因素、水平确定选择合适正交表表头设计列试验方案试验方案设计：试验结果分析第二十二张，PPT共一百二十三页，创作于2022

14、年6月进行试验，记录试验结果试验结果极差分析计算K值计算k值计算极差R绘制因素指标趋势图优水平因素主次顺序优组合结 论试验结果分析：试验结果方差分析列方差分析表，进行F 检验计算各列偏差平方和、自由度分析检验结果，写出结论第二十三张，PPT共一百二十三页，创作于2022年6月2.1 试验方案设计实例：为提高山楂原料的利用率，研究酶法液化工艺制造山楂原汁，拟通过正交试验来寻找酶法液化的最佳工艺条件。 试验设计前必须明确试验目的，即本次试验要解决什么问题。试验目的确定后，对试验结果如何衡量，即需要确定出试验指标。试验指标可为定量指标，如强度、硬度、产量、出品率、成本等；也可为定性指标如颜色、口感、

15、光泽等。一般为了便于试验结果的分析，定性指标可按相关的标准打分或模糊数学处理进行数量化，将定性指标定量化。（1） 明确试验目的，确定试验指标 对本试验而言，试验目的是为了提高山楂原料的利用率。所以可以以液化率液化率=(果肉重量-液化后残渣重量)/果肉重量100%为试验指标，来评价液化工艺条件的好坏。液化率越高，山楂原料利用率就越高。第二十四张，PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 根据专业知识、以往的研究结论和经验，从影响试验指标的诸多因素中，通过因果分析筛选出需要考察的试验因素。一般确定试验因素时，应以对试验指标影响大的因素、尚未考察过的因素、尚未完全掌握

16、其规律的因素为先。试验因素选定后，根据所掌握的信息资料和相关知识，确定每个因素的水平，一般以2-4个水平为宜。对主要考察的试验因素，可以多取水平，但不宜过多（6），否则试验次数骤增。因素的水平间距，应根据专业知识和已有的资料，尽可能把水平值取在理想区域。（2） 选因素、定水平，列因素水平表 第二十五张，PPT共一百二十三页，创作于2022年6月 对本试验分析，影响山楂液化率的因素很多，如山楂品种、山楂果肉的破碎度、果肉加水量、原料pH 值、果胶酶种类、加酶量、酶解温度、酶解时间等等。经全面考虑，最后确定果肉加水量、加酶量、酶解温度和酶解时间为本试验的试验因素，分别记作A、B、C和D，进行四因素

17、正交试验，各因素均取三个水平，因素水平表见表10-3所示。 水平试验因素加水量（mL/100g）A加酶量（mL/100g）B酶解温度（）C酶解时间（h）D1101201.52504352.53907503.510-3 因素水平表第二十六张，PPT共一百二十三页，创作于2022年6月 正交表的选择是正交试验设计的首要问题。确定了因素及其水平后，根据因素、水平及需要考察的交互作用的多少来选择合适的正交表。正交表的选择原则是在能够安排下试验因素和交互作用的前提下，尽可能选用较小的正交表，以减少试验次数。 一般情况下，试验因素的水平数应等于正交表中的水平数；因素个数（包括交互作用）应不大于正交表的列数

18、；各因素及交互作用的自由度之和要小于所选正交表的总自由度，以便估计试验误差。若各因素及交互作用的自由度之和等于所选正交表总自由度，则可采用有重复正交试验来估计试验误差。（3） 选择合适的正交表第二十七张，PPT共一百二十三页，创作于2022年6月La（bc）正交设计试验总次数，行数因素水平数因素个数，列数等水平正交表 La（bc）第二十八张，PPT共一百二十三页，创作于2022年6月列：正交表的列数c因素所占列数+交互作用所占列数+空列。自由度：正交表的总自由度（a-1）因素自由度+交互作用自由度+误差自由度。正交表选择依据：第二十九张，PPT共一百二十三页，创作于2022年6月 此例有4个3

19、水平因素，可以选用L9(34)或L27(313) ；因本试验仅考察四个因素对液化率的影响效果，不考察因素间的交互作用，故宜选用L9（34）正交表。若要考察交互作用，则应选用L27(313)。 下一张 主 页 退 出 上一张 第三十张，PPT共一百二十三页，创作于2022年6月 所谓表头设计，就是把试验因素和要考察的交互作用分别安排到正交表的各列中去的过程。 在不考察交互作用时，各因素可随机安排在各列上；若考察交互作用，就应按所选正交表的交互作用列表安排各因素与交互作用，以防止设计“混杂” 。 此例不考察交互作用，可将加水量(A)、加酶量(B)和酶解温度 (C)、酶解时间（D）依次安排在L9(3

20、4)的第1、2、3、4列上，见表10-4所示。（4） 表头设计列号1234因素ABCD表10-4 表头设计第三十一张，PPT共一百二十三页，创作于2022年6月 把正交表中安排各因素的列（不包含欲考察的交互作用列）中的每个水平数字换成该因素的实际水平值，便形成了正交试验方案（表10-5）。下一张 主 页 退 出 上一张 （5）编制试验方案，按方案进行试验，记录试验结果。第三十二张，PPT共一百二十三页，创作于2022年6月试验号因 素ABCD111112122231333421235223162312731328321393321表10-5 试验方案及试验结果说明：试验号并非试验顺序，为了排除

21、误差干扰，试验中可随机进行； 安排试验方案时，部分因素的水平可采用随机安排。1（10）2（50）3（90）2（4）3（7）1（1）2（35）1（20）3（50）3(3.5)2(2.5)1(1.5)试验结果（液化率 %）0172412472811842第三十三张，PPT共一百二十三页，创作于2022年6月例10-2 鸭肉保鲜天然复合剂的筛选。试验以茶多酚作为天然复合保鲜剂的主要成分，分别添加不同增效剂、被膜剂和不同的浸泡时间，进行4因素4水平正交试验。试设计试验方案。（西南农业大学）有机酸和盐处理对鸭肉保鲜有明显效果，但大部分属于合成的化学试剂，在卫生安全上得不到保证，并且不符合满足消费者纯天然

22、、无污染的要求。 明确目的，确定指标。本例的目的是通过试验，寻找一个最佳的鸭肉天然复合保鲜剂。 选因素、定水平。根据专业知识和以前研究结果，选择4个因素，每个因素定4个水平，因素水平表见表10-6。第三十四张，PPT共一百二十三页，创作于2022年6月 选择正交表。此试验为4因素4水平试验，不考虑交互作用，4因素共占4列，选L16（45）最合适，并有1空列，可以作为试验误差以衡量试验的可靠性。 表头设计。4因素任意放置。 编制试验方案。试验方案见表10-7。水平因素A茶多酚浓度/B增效剂种类C被膜剂种类D浸泡时间/min10.10.5维生素C0.5海藻酸钠120.20.1柠檬酸0.8海藻酸钠2

23、30.30.2-CD1.0海藻酸钠340.4生姜汁1.0葡萄糖4表10-6 天然复合保鲜剂筛选试验因素水平表第三十五张，PPT共一百二十三页，创作于2022年6月试验号A茶多酚浓度/B增效剂种类C被膜剂种类D浸泡时间/minE 空列结果11233236.2022412231.5433434330.0944211329.3251314431.7762131435.0273113132.3784332132.6491142338.79102323330.90113341232.87124124234.54131421138.02142244135.62153222434.02164443432.8

24、0表10-7 天然复合保鲜剂筛选试验方案第三十六张，PPT共一百二十三页，创作于2022年6月2.2 试验结果分析分清各因素及其交互作用的主次顺序，分清哪个是主要因素，哪个是次要因素；判断因素对试验指标影响的显著程度；找出试验因素的优水平和试验范围内的最优组合，即试验因素各取什么水平时，试验指标最好；分析因素与试验指标之间的关系，即当因素变化时，试验指标是如何变化的。找出指标随因素变化的规律和趋势，为进一步试验指明方向；了解各因素之间的交互作用情况；估计试验误差的大小。极差分析方差分析第三十七张，PPT共一百二十三页，创作于2022年6月Kjm，kjm 计算简便，直观，简单易懂，是正交试验结果

25、分析最常用方法。以上例为实例来说明极差分析过程。 3 正交试验的结果分析3.1 直观分析法极差分析法极差分析法R法1. 计算2. 判断Rj因素主次优水平优组合Kjm为第j列因素m水平所对应的试验指标和，kjm为Kjm平均值。由kjm大小可以判断第j列因素优水平和优组合。Rj为第j列因素的极差，反映了第j列因素水平波动时，试验指标的变动幅度。Rj越大，说明该因素对试验指标的影响越大。根据Rj大小，可以判断因素的主次顺序。第三十八张，PPT共一百二十三页，创作于2022年6月（1） 确定试验因素的优水平和最优水平组合 分析A因素各水平对试验指标的影响。由表3可以看出，A1的影响反映在第1、2、3号

26、试验中，A2的影响反映在第4、5、6号试验中，A3的影响反映在第7、8、9号试验中。A因素的1水平所对应的试验指标之和为KA1=y1+y2+y3=0+17+24=41，kA1= KA1/3=13.7；A因素的2水平所对应的试验指标之和为KA2=y4+y5+y6=12+47+28=87，kA2=KA2/3=29；A因素的3水平所对应的试验指标之和为KA3=y7+y8+y9=1+18+42=61，kA3=KA3/3=20.3。3.1.1 不考察交互作用的试验结果分析第三十九张，PPT共一百二十三页，创作于2022年6月 根据正交设计的特性，对A1、A2、A3来说，三组试验的试验条件是完全一样的（综

27、合可比性），可进行直接比较。如果因素A对试验指标无影响时，那么kA1、kA2、kA3应该相等，但由上面的计算可见，kA1、kA2、kA3实际上不相等。说明，A因素的水平变动对试验结果有影响。因此，根据kA1、kA2、kA3的大小可以判断A1、A2、A3对试验指标的影响大小。由于试验指标为液化率，而kA2kA3kA1，所以可断定A2为A因素的优水平。 第四十张，PPT共一百二十三页，创作于2022年6月同理，可以计算并确定B3、C3、D1分别为B、C、D因素的优水平。四个因素的优水平组合A2B3C3D1为本试验的最优水平组合，即酶法液化生产山楂清汁的最优工艺条件为加水量50mL/100g，加酶量

28、7mL/100g，酶解温度为50,酶解时间为1.5h。第四十一张，PPT共一百二十三页，创作于2022年6月 根据极差Rj的大小，可以判断各因素对试验指标的影响主次。本例极差Rj计算结果见表10-8，比较各R值大小，可见RBRARDRC,所以因素对试验指标影响的主次顺序是BADC。即加酶量影响最大，其次是加水量和酶解时间，而酶解温度的影响较小。（2） 确定因素的主次顺序 以各因素水平为横坐标，试验指标的平均值（kjm）为纵坐标，绘制因素与指标趋势图。由因素与指标趋势图可以更直观地看出试验指标随着因素水平的变化而变化的趋势，可为进一步试验指明方向。（3） 绘制因素与指标趋势图以上即为正交试验极差

29、分析的基本程序与方法第四十二张，PPT共一百二十三页，创作于2022年6月表10-8 试验结果分析试验号因素液化率ABCD1111102122217313332442123125223147623122873132183213189332142K141134689K287827146K361947254k113.74.315.329.7k229.027.323.715.3k320.331.324.018.0极差R15.327.08.714.3主次顺序BADC优水平A2B3C3D1优组合A2B3C3D1第四十三张，PPT共一百二十三页，创作于2022年6月表10-8 试验结果分析第四十四张，PP

30、T共一百二十三页，创作于2022年6月（2）计算各因素同一水平的平均值Ki。K1=36.20，K2=33.27，K3=32.34，K4=31.83例10-2试验结果极差分析（1）计算Ki值。Ki为同一水平之和。以第一列A因素为例：K1=36.20+31.77+38.79+38.02=144.78K2=31.54+35.02+30.90+35.62=133.08K3=30.09+32.37+32.87+34.02=129.35K4=29.32+32.64+34.54+32.80=129.30第四十五张，PPT共一百二十三页，创作于2022年6月（3）计算各因素的极差R，R表示该因素在其取值范围内

31、试验指标变化的幅度。 R=max（Ki）-min（Ki）（4）根据极差大小，判断因素的主次影响顺序。R越大，表示该因素的水平变化对试验指标的影响越大，因素越重要。由以上分析可见，因素影响主次顺序为A-C-B-D，A因素影响最大，为主要因素，D因素为不重要因素。（5）做因素与指标趋势图，直观分析出指标与各因素水平波动的关系。第四十六张，PPT共一百二十三页，创作于2022年6月（6）选优组合，即根据各因素各水平的平均值确定优水平，进而选出优组合。 本例A、B、C为主要因素，按照平均值大小选取优水平为A1B1C4，即茶多酚用量取0.1%水平；以0.5%维生素C作为增效剂；1.0%葡萄糖液为被膜剂为

32、形成的鸭肉保鲜复合剂为优组合，而浸泡时间为次要因素，选取操作时间1-3min即可。第四十七张，PPT共一百二十三页，创作于2022年6月表10-9 鸭肉保鲜天然复合剂筛选试验结果第四十八张，PPT共一百二十三页，创作于2022年6月附1 ： 多指标正交试验极差分析 对于多指标试验，方案设计和实施与单指标试验相同，不同在于每做一次试验，都需要对考察指标一一测试，分别记录。试验结果分析时，也要对考察指标一一分析，然后综合评衡，确定出优条件。 油炸方便面生产中，主要原料质量和主要工艺参数对产品质量有影响。通过试验确定最佳生产条件。第四十九张，PPT共一百二十三页，创作于2022年6月（1）试验方案设

33、计确定试验指标。本试验目的是探讨方便面生产的最佳工艺条件，以提高方便面的质量。试验以脂肪含量、水分含量和复水时间指标。脂肪含量越低越好，水分含量越高越好，复水时间越短越好。挑因素，选水平，列因素水平表。根据专业知识和实践经验，确定试验因素和水平见表10-10。表10-10 因素水平表第五十张，PPT共一百二十三页，创作于2022年6月选正交表、设计表头、编制试验方案。本试验为四因素三水平试验，不考虑交互作用，选L9（34）安排试验。表头设计和试验方案以及试验结果记录见表。（2）试验结果分析计算各因素各水平下每种试验指标的数据和以及平均值，并计算极差R。根据极差大小列出各指标下的因素主次顺序。试

34、验指标： 主次顺序脂肪含量（）：ACDB水分含量（）：CDAB复水时间（s）：ADBC第五十一张，PPT共一百二十三页，创作于2022年6月表10-11 试验结果极差分析表第五十二张，PPT共一百二十三页，创作于2022年6月初选优化工艺条件。根据各指标不同水平平均值确定各因素的优化水平组合。脂肪含量（）：A3B3C1D2水分含量（）：A1B2C1D1复水时间（ s ）：A2B2C2D3综合平衡确定最优工艺条件。以上三指标单独分析出的优化条件不一致，必须根据因素的影响主次，综合考虑，确定最佳工艺条件。 对于因素A，其对粗脂肪影响大小排第一位，此时取A3；其对复水时间影响也排第一位，取A2；而其

35、对水分影响排次要第三位，为次要因素，因此A可取A2或A3，但取A2时，复水时间比取A3缩短了14%，而粗脂肪增加了11.3%，且由水分指标看，取A2比A3水分高，故A因素取A2。同理可分析B取B2，C取C1，D取D3。优组合为A2B2C1D3.第五十三张，PPT共一百二十三页，创作于2022年6月附2 ： 混合型正交表试验设计与极差分析试验设计与结果分析同前。 某油炸膨化食品的体积与油温、物料含水量及油炸时间有关，为确保产品质量，现通过正交试验来寻求理想的工艺参数。表10-12 因素水平表第五十四张，PPT共一百二十三页，创作于2022年6月表10-13 试验方案及结果分析结论：油炸温度对油炸

36、食品的体积影响最大，其次是油炸时间，而物料含水量影响最小。优化组合为A3B2C2或A3B1C2，即理想工艺参数为油炸温度230，油炸时间40s，物料含水量可取2%或4%。r 为因素每个水平试验重复数d 折算系数，与因素水平有关。表10-14 折算系数表第五十五张，PPT共一百二十三页，创作于2022年6月（1） 交互作用 在多因素试验中，不仅因素对指标有影响，而且因素之间的联合搭配也对指标产生影响。因素间的联合搭配对试验指标产生的影响作用称为交互作用。因素之间的交互作用总是存在的，这是客观存在的普遍现象，只不过交互作用的程度不同而异。一般地，当交互作用很小时，就认为因素间不存在交互作用。对于交

37、互作用，设计时应引起高度重视。 在试验设计中，表示A、B间的交互作用记作AB，称为1级交互作用；表示因素A、B、C之间的交互作用记作ABC，称为2级交互作用；依此类推，还有3级、4级交互作用等。3.1.2 考察交互作用的试验设计与结果分析第五十六张，PPT共一百二十三页，创作于2022年6月（2）交互作用的处理原则 试验设计中，交互作用一律当作因素看待，这是处理交互作用问题的总原则。作为因素，各级交互作用都可以安排在能考察交互作用的正交表的相应列上，它们对试验指标的影响情况都可以分析清楚，而且计算非常简单。但交互作用又与因素不同，表现在： 用于考察交互作用的列不影响试验方案及其实施； 一个交互

38、作用并不一定只占正交表的一列，而是占有（m-1）p列。表头设计时，交互作用所占列数与因素的水平m有关，与交互作用级数p有关。第五十七张，PPT共一百二十三页，创作于2022年6月 2水平因素的各级交互作用均占1列；对于3水平因素，一级交互作用占两列，二级交互作用占四列，可见，m和p越大，交互作用所占列数越多。 例如，对一个25因素试验，表头设计时，如果考虑所有各级交互作用，那么连同因素本身，总计应占列数为： C51 + C52 +C53 +C54 +C55 5+10+10+5+131，那么此试验必选L32（24）正交表进行设计。一般对于多因素试验，在满足试验要求的条件下，有选择地、合理地考察某

39、些交互作用。第五十八张，PPT共一百二十三页，创作于2022年6月 综合考虑试验目的、专业知识、以往的经验及现有试验条件等多方面情况进行交互作用选择。一般原则是： 忽略高级交互作用 有选择地考察一级交互作用。通常只考察那些作用效果较明显的，或试验要求必须考察的。 试验允许的条件下，试验因素尽量取2水平。第五十九张，PPT共一百二十三页，创作于2022年6月（3）有交互作用的试验表头设计 表头设计时，各因素及其交互作用不能任意安排，必须严格按交互作用列表进行安排。这是有交互作用正交试验设计的一个重要特点，也是关键的一步。 在表头设计中，为了避免混杂，那些主要因素，重点要考察的因素，涉及交互作用较

40、多的因素，应该优先安排，次要因素，不涉及交互作用的因素后安排。 所谓混杂，就是指在正交表的同列中，安排了两个或两个以上的因素或交互作用，这样，就无法区分同一列中这些不同因素或交互作用对试验指标的影响效果。第六十张，PPT共一百二十三页，创作于2022年6月 在实际研究中，有时试验因素之间存在交互作用。对于既考察因素主效应又考察因素间交互作用的正交设计，除表头设计和结果分析与前面介绍略有不同外，其它基本相同。 【例】 某一种抗菌素的发酵培养基由A、B、C 三种成分组成，各有两个水平，除考察A、B、C三个因素的主效外，还考察A与B、B与C的交互作用。试安排一个正交试验方案并进行结果分析。 下一张

41、主 页 退 出 上一张 （4）有交互作用的正交设计与分析实例第六十一张，PPT共一百二十三页，创作于2022年6月 选用正交表，作表头设计 由于本试验有3个两水平的因素和两个交互作用需要考察，各项自由度之和为：3(2-1)+2(2-1)(2-1)=5，因此可选用L8(27)来安排试验方案。 正交表L8(27)中有基本列和交互列之分，基本列就是各因素所占的列，交互列则为两因素交互作用所占的列。可利用L8(27)二列间交互作用列表来安排各因素和交互作用。下一张 主 页 退 出 上一张 第六十二张，PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 第六十三张，PPT共一百二十

42、三页，创作于2022年6月 如果将A因素放在第1列 ，B 因素 放在第 2列，查表可知，第1列与第2列的交互作用列是第3列 ，于是将 A与B 的交互作用 AB放在第3列。这样第3列不能再安排其它因素 ，以免出现“混杂”。然后将C放在第4列， 查表 12-30 可知，BC应放在第6列，余下列为空列 ，如此可得表头设计，见表10-15。 下一张 主 页 退 出 上一张 第六十四张，PPT共一百二十三页，创作于2022年6月 列出试验方案 根据表头设计，将A、B、C各列对应的数字“1”、“2”换成各因素的具体水平，得出试验方案列于表10-16。下一张 主 页 退 出 上一张 表10-15第六十五张，

43、PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 表10-16第六十六张，PPT共一百二十三页，创作于2022年6月 结果分析 按表所列的试验方案进行试验，其结果分析与前面并无本质区别，只是：应把互作当成因素处理进行分析； 应根据互作效应，选择优化组合。 下一张 主 页 退 出 上一张 第六十七张，PPT共一百二十三页，创作于2022年6月下一张 主 页 退 出 上一张 *试验结果以对照为100计。试验号ABABC空列BC空列试验结果111111115521112222383122112297412222118952121212122621221211247221122

44、1798221211261K1279339233353337327347K2386326432312328338318k169.75 84.75 58.25 88.25 84.25 81.75 86.75 k296.50 81.50 108.00 78.00 82.00 84.50 79.50 极差R 26.75 3.25 49.75 10.25 2.25 2.75 7.25 主次顺序ABACBBC优水平A2B1C1优组合A2B1C1表10-17 极差分析结果因素主次顺序为ABACBBC，表明AB交互作用、 A因素影响最大，因素C影响次之，因素B影响最小。优组合为A2B1C1。二元表B1B2A

45、146.593A212370第六十八张，PPT共一百二十三页，创作于2022年6月例：p348 要生产每种食品添加剂，根据试验发现影响添加剂得率的因素有4个，每个因素设置2水平。因素水平表见表10-18。试验中可考虑交互作用AB、AC、BC。水平试验因素温度A /时间B /h配比C（两种原料）真空度C /kPa17522:0153.3229033:0166.65表10-18 某种食品添加剂得率试验因素水平表正交表的选择：自由度：dfT 因素+交互作用+空列4*（2-1）+3*1+17+18那么正交表的行数a dfT +19 无空列时a 8，选L8（27）即可。第六十九张，PPT共一百二十三页，

46、创作于2022年6月列：c因素所占列+交互作用所占列+误差列（空列）因素列：各因素各占一列，共计4列（4个因素）交互作用列：因试验因素为2水平因素，其1级交互作用分占1列，共计3列（3组交互作用）。误差列：0或1列c4+3+07，因素水平为2，列为7的最小正交表即L8(27)。可以看出尚无空列估计试验误差，应做重复试验或忽略某些交互作用。第七十张，PPT共一百二十三页，创作于2022年6月试验号ABABCACBCD试验结果11111111862111222295312211229141222211945212121291621221219672211221838221211288K136636

47、8352351361359359K2358356372373363365365k191.5 92.0 88.0 87.8 90.3 89.8 89.8 k289.5 89.0 93.0 93.3 90.8 91.3 91.3 极差R 2.0 3.0 5.0 5.5 0.5 1.5 1.5 主次顺序CABBABC、DAC优水平A2B1C2D1或D2优组合A2B1C2D1或D2表10-19 食品添加剂得率试验结果极差分析因素主次顺序为CABBABC、D AC ，表明C影响最大，AB交互作用影响其次，为重要考察因素；AC、BC、D等影响小，为次要因素， AC、BC交互作用是由误差引起的，可以忽略。表

48、10-16 二元表A1A2B190.593.5B292.585.5结论：优组合为A2B1C2D1或A2B1C2D2第七十一张，PPT共一百二十三页，创作于2022年6月 极差分析法简单明了，通俗易懂，计算工作量少便于推广普及。但这种方法不能将试验中由于试验条件改变引起的数据波动同试验误差引起的数据波动区分开来，也就是说，不能区分因素各水平间对应的试验结果的差异究竟是由于因素水平不同引起的，还是由于试验误差引起的，无法估计试验误差的大小。此外，各因素对试验结果的影响大小无法给以精确的数量估计，不能提出一个标准来判断所考察因素作用是否显著。为了弥补极差分析的缺陷，可采用方差分析。下一张 主 页 退

49、 出 上一张 3.2 正交试验结果的方差分析第七十二张，PPT共一百二十三页，创作于2022年6月 3.2.1 正交试验结果的方差分析 方差分析基本思想是将数据的总变异分解成因素引起的变异和误差引起的变异两部分，构造F统计量，作F检验，即可判断因素作用是否显著。正交试验结果的方差分析思想、步骤同前！第七十三张，PPT共一百二十三页，创作于2022年6月总偏差平方和各列因素偏差平方和+误差偏差平方和（1）偏差平方和分解：（2）自由度分解：（3）方差：第七十四张，PPT共一百二十三页，创作于2022年6月（4）构造F统计量：（5）列方差分析表，作F检验若计算出的F值F0Fa，则拒绝原假设，认为该因

50、素或交互作用对试验结果有显著影响；若F0Fa，则认为该因素或交互作用对试验结果无显著影响。第七十五张，PPT共一百二十三页，创作于2022年6月（6）正交试验方差分析说明由于进行F检验时，要用误差偏差平方和SSe及其自由度dfe，因此，为进行方差分析，所选正交表应留出一定空列。当无空列时，应进行重复试验，以估计试验误差。误差自由度一般不应小于2，dfe很小，F检验灵敏度很低，有时即使因素对试验指标有影响，用F检验也判断不出来。为了增大dfe，提高F检验的灵敏度，在进行显著性检验之前，先将各因素和交互作用的方差与误差方差比较，若MS因（MS交） 2MSe，可将这些因素或交互作用的偏差平方和、自由

51、度并入误差的偏差平方和、自由度，这样使误差的偏差平方和和自由度增大，提高了F检验的灵敏度。第七十六张，PPT共一百二十三页，创作于2022年6月表10-20 L9(34)正交表处理号 第1列（A） 第2列 第3列 第4列 试验结果yi11111y121222y231333y342123y452231y562312y673132y783213y893321y9分析第1列因素时，其它列暂不考虑，将其看做条件因素。因素A第1水平3次重复测定值因素A第2水平3次重复测定值因素A第3水平3次重复测定值因素重复1重复2重复3A1y1y2y3A2y4y5y6A3y7y8y9单因素试验数据资料格式和y1+y2

52、+y3K1y4+y5+y6K2y7+y8+y9K3第七十七张，PPT共一百二十三页，创作于2022年6月表头设计AB试验数据列号12kxixi2试验号11x1x1221x2x22nmxnxn2K1jK11K12K1kK2jK21K22K2kKmjKm1Km2KmkK1j2K112K122K1k2K2j2K212K222K2k2Kmj2Km12Km22Kmk2SSjSS1SS2SSk表10-21 Ln（mk）正交表及计算表格第七十八张，PPT共一百二十三页，创作于2022年6月总偏差平方和：列偏差平方和： 试验总次数为n，每个因素水平数为m个，每个水平作r次重复rn/m。当m2时，第七十九张，P

53、PT共一百二十三页，创作于2022年6月总自由度：因素自由度：第八十张，PPT共一百二十三页，创作于2022年6月 3.2.2 不考虑交互作用等水平正交试验方差分析 例：自溶酵母提取物是一种多用途食品配料。为探讨啤酒酵母的最适自溶条件，安排三因素三水平正交试验。试验指标为自溶液中蛋白质含量（）。试验因素水平表见表10-22，试验方案及结果分析见表10-23。试对试验结果进行方差分析。水 平试验因素温度（）ApH值B加酶量（）C1506.52.02557.02.43587.52.8表10-22 因素水平表第八十一张，PPT共一百二十三页，创作于2022年6月处理号 ABC空列试验结果yi11（5

54、0）1（6.5）1（2.0）16.25212（7.0）2（2.4）24.97313（7.5）3（2.834.5442（55）1237.53522315.54623125.573（58）13211.48321310.9933218.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 表10-23

55、 试验方案及结果分析表第八十二张，PPT共一百二十三页，创作于2022年6月（1）计算计算各列各水平的K值 计算各列各水平对应数据之和K1j、K2j、K3j及其平方K1j2、K2j2、K3j2。计算各列偏差平方和及自由度同理，SSB=6.49，SSC=0.31 SSe=0.83（空列）第八十三张，PPT共一百二十三页，创作于2022年6月自由度：dfAdfBdfCdfe3-1=2计算方差（2）显著性检验根据以上计算，进行显著性检验，列出方差分析表，结果见表10-24变异来源 平方和 自由度 均方 F值 Fa显著水平 A45.40222.7079.6F0.05(2,4) =6.94*B6.492

56、3.2411.4F0.01(2,4)=18.0*C0.3120.16误差e0.8320.41误差e 1.1440.285总和 53.03表10-24 方差分析表第八十四张，PPT共一百二十三页，创作于2022年6月因素A高度显著，因素B显著，因素C不显著。因素主次顺序A-B-C。（3）优化工艺条件的确定 本试验指标越大越好。对因素A、B分析，确定优水平为A3、B1；因素C的水平改变对试验结果几乎无影响，从经济角度考虑，选C1。优水平组合为A3B1C1。即温度为58，pH值为6.5，加酶量为2.0%。第八十五张，PPT共一百二十三页，创作于2022年6月 3.2.3 考虑交互作用正交试验方差分析

57、 例： 用石墨炉原子吸收分光光度法测定食品中的铅，为了提高测定灵敏度，希望吸光度越大越好，今欲研究影响吸光度的因素，确定最佳测定条件。（1）计算 计算各列各水平对应数据之和K1j、K2j及（K1j-K2j）；计算各列偏差平方和及自由度。第八十六张，PPT共一百二十三页，创作于2022年6月表10-25 试验方案及结果分析表试验号ABABCACBC空列吸光度111111112.42211122222.24312211222.66412222112.58521212122.36621221212.4722112212.79822121122.76K1j9.99.4210.2110.2310.241

58、0.1210.19K2j10.3110.79109.989.9710.0910.02K1j-K2j-0.41-1.370.210.250.270.030.17SSj0.0210.2350.00550.00780.00910.00010.0036第八十七张，PPT共一百二十三页，创作于2022年6月变异来源 平方和 自由度 均方 F值 临界值Fa显著水平 A0.0210 10.021 6.82F0.05(1,3)=10.13B0.2346 10.235 76.19F0.01(1,3)=34.12*AB0.0055 10.006 C0.0078 10.008 2.53AC0.0091 10.009

59、 2.96BC 0.0001 10.000 误差e0.0036 10.004 误差e 0.0923 30.00308 总 和 0.2818 表10-26 方差分析表（2）显著性检验因素B高度显著，因素A、C及交互作用AB、AC、BC均不显著。各因素对试验结果影响的主次顺序为：B、A、AC、C、AB、BC。第八十八张，PPT共一百二十三页，创作于2022年6月（3）优化条件确定 交互作用均不显著，确定因素的优水平时可以不考虑交互作用的影响。对显著因素B，通过比较K1B和K2B的大小确定优水平为B2；同理A取A2，C取C1或C2。优组合为A2B2C1或A2B2C2。方差分析可以分析出试验误差的大小

60、，从而知道试验精度；不仅可给出各因素及交互作用对试验指标影响的主次顺序，而且可分析出哪些因素影响显著，哪些影响不显著。对于显著因素，选取优水平并在试验中加以严格控制；对不显著因素，可视具体情况确定优水平。但极差分析不能对各因素的主要程度给予精确的数量估计。第八十九张，PPT共一百二十三页，创作于2022年6月 3.2.4 混合型正交试验方差分析 混合型正交试验方差分析与等水平正交试验方差分析没有本质区别。（1）计算二水平列：第九十张，PPT共一百二十三页，创作于2022年6月试验号油温A含水量B油炸时间s C空列空列试验指标11111112122220.83211221.54222113531