《正交试验设计方法13192.pdf》由会员分享，可在线阅读，更多相关《正交试验设计方法13192.pdf（10页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-第 5 章 正交试验设计方法 51 试验设计方法概述 试验设计是数理统计学的一个重要的分支。多数数理统计方法主要用于分析已经得到的数据，而试验设计却是用于决定数据收集的方法。试验设计方法主要讨论如何合理地安排试验以及试验所得的数据如何分析等。例 5-1 *化工厂想提高*化工产品的质量和产量，对工艺中三个主要因素各按三个水平进展试验见表 5-1。试验的目的是为提高合格产品的产量，寻求最适宜的操作条件。对此实例该如何进展试验方案的设计呢？很容易想到的是全面搭配法方案如图 5-1 所示：此方案数据点分布的均匀性极好，因素和水平的搭配十分全面，唯一的缺点是实验次数多达3327次指数 3 代表 3 个
2、因素，底数 3 代表每因素有 3 个水平。因素、水平数愈多，则实验次数就愈多，例如，做一个 6 因素 3 水平的试验，就需 36729 次实验，显然难以做到。因此需要寻找一种适宜的试验设计方法。试验设计方法常用的术语定义如下。试验指标：指作为试验研究过程的因变量，常为试验结果特征的量 如得率、纯度等。例 1 的试验指标为合格产品的产量。因素：指作试验研究过程的自变量，常常是造成试验指标按*种规律发生变化的那些原因。如例 1 的温度、压力、碱的用量。水平：指试验中因素所处的具体状态或情况，又称为等级。如例 1 的温度有 3 个水平。温度用T表示，下标 1、2、3 表示因素的不同水平，分别记为T1
3、、T2、T3。常用的试验设计方法有：正交试验设计法、均匀试验设计法、单纯形优化法、双水平单纯形优化法、回归正交设计法、序贯试验设计法等。可供选择的试验方法很多，各种试表 51 因素水平 水平 因素 温度 压力 Pa 加碱量 kg 符号 T p m 1 2 3 T1(80)T2(100)T3(120)p1(5.0)p2(6.0)p3(7.0)m 1(2.0)m2(2.5)m3(3.0)图 51 全面搭配法方案-验设计方法都有其一定的特点。所面对的任务与要解决的问题不同，选择的试验设计方法也应有所不同。由于篇幅的限制，我们只讨论正交试验设计方法。52 正交试验设计方法的优点和特点 用正交表安排多因
4、素试验的方法，称为正交试验设计法。其特点为：完成试验要求所需的实验次数少。数据点的分布很均匀。可用相应的极差分析方法、方差分析方法、回归分析方法等对试验结果进展分析，引出许多有价值的结论。从例 1 可看出，采用全面搭配法方案，需做 27 次实验。则采用简单比拟法方案又如何呢？先固定T1和p1，只改变m，观察因素m不同水平的影响，做了如图 2-21所示的三次实验，发现 mm2时的实验效果最好好的用表示，合格产品的产量最高，因此认为在后面的实验中因素m应取m2水平。固定T1和m2，改变p的三次实验如图 5-22 所示，发现pp时的实验效果最好，因此认为因素p应取p水平。固定p和m2，改变T的三次实
5、验如图 5-23所示，发现因素T宜取T2水平。因此可以引出结论：为提高合格产品的产量，最适宜的操作条件为T2pm2。与全面搭配法方案相比，简单比拟法方案的优点是实验的次数少，只需做 9 次实验。但必须指出，简单比拟法方案的试验结果是不可靠的。因为，在改变m值或p值，或T值的三次实验中，说m2或p或T2 水平最好是有条件的。在TT，pp时，m2 水平不是最好的可能性是有的。在改变 m 的三次实验中，固定T T，p p 应该说也是可以的，是随意的，故在此方案中数据点的分布的均匀性是毫无保障的。用这种方法比拟条件好坏时，只是对单个的试验数据进展数值上的简单比拟，不能排除必然存在的试验数据误差的干扰。
6、运用正交试验设计方法，不仅兼有上述两个方案的优点，而且实验次数少，数据点分布均匀，结论的可靠性较好。正交试验设计方法是用正交表来安排试验的。对于例 1 适用的正交表是 L934，其试验安排见表 5-2。所有的正交表与 L934正交表一样，都具有以下两个特点：1 在每一列中，各个不同的数字出现的次数一样。在表 L934中，每一列有三个水平，水平 1、2、3 都是各出现 3 次。2 表中任意两列并列在一起形成假设干个数字对，不同数字对出现的次数也都一样。在表 L934 中，任意两列并列在一起形成的数字对共有 9 个：1,1，1,2，1,3，2,1，2,2，2,3，3,1，3,2，3,3，每一个数字
7、对各出现一次。图 52 简单比拟法方案-L 93 4 正交表的列数 每一列的水平数 实验的次数 正交表的代号 表 52 试验安排表 试验号 列号 1 2 3 4 因素 温度 压力 Pa 加碱量 kg 符号 T p m 1 2 3 4 5 6 7 8 9 1T1 1T1 1T1 2T2 2T2 2T2 3T3 3T3 3T3 1p1 2p2 3p3 1p1 2p2 3p3 1p1 2p2 3p3 1m1 2m2 3m3 2m2 3m3 1m1 3m3 1m1 2m2 1 2 3 3 1 2 2 3 1 这两个特点称为正交性。正是由于正交表具有上述特点，就保证了用正交表安排的试验方案中因素水平是均
8、衡搭配的，数据点的分布是均匀的。因素、水平数愈多，运用正交试验设计方法，愈发能显示出它的优越性，如上述提到的 6 因素 3 水平试验，用全面搭配方案需 729 次，假设用正交表 L27313来安排，则只需做 27 次试验。在化工生产中，因素之间常有交互作用。如果上述的因素T的数值和水平发生变化时，试验指标随因素p变化的规律也发生变化，或反过来，因素p的数值和水平发生变化时，试验指标随因素T变化的规律也发生变化。这种情况称为因素T、p间有交互作用，记为Tp。5 3 正交表 使用正交设计方法进展试验方案的设计，就必须用到正交表。正交表请查阅有关参考书。5.3.1 各列水平数均一样的正交表 各列水平
9、数均一样的正交表，也称单一水平正交表。这类正交表名称的写法举例如下：各列水平均为 2 的常用正交表有：L423，L827，L12211，L16215，L20219，L32231。各列水平数均为 3 的常用正交表有：L934，L27313。各列水平数均为 4 的常用正交表有：L1645 各列水平数均为 3 的常用正交表有：L2556-5.3.2 混合水平正交表 各列水平数不一样的正交表，叫混合水平正交表，下面就是一个混合水平正交表名称的写法：L 84124常简写为 L 8424。此混合水平正交表含有 1 个 4 水平列，4 个 2 水平列，共有 145 列。5.3.3 选择正交表的根本原则 一般
10、都是先确定试验的因素、水平和交互作用，后选择适用的 L 表。在确定因素的水平数时，主要因素宜多安排几个水平，次要因素可少安排几个水平。1先看水平数。假设各因素全是 2 水平，就选用 L(2)表；假设各因素全是 3 水平，就选 L(3)表。假设各因素的水平数不一样，就选择适用的混合水平表。2每一个交互作用在正交表中应占一列或二列。要看所选的正交表是否足够大，能否容纳得下所考虑的因素和交互作用。为了对试验结果进展方差分析或回归分析，还必须至少留一个空白列，作为误差列，在极差分析中要作为其他因素列处理。3要看试验精度的要求。假设要求高，则宜取实验次数多的 L 表。4假设试验费用很昂贵，或试验的经费很
11、有限，或人力和时间都比拟紧，则不宜选实验次数太多的 L 表。5按原来考虑的因素、水平和交互作用去选择正交表，假设无正好适用的正交表可选，简便且可行的方法是适当修改原定的水平数。6 对*因素或*交互作用的影响是否确实存在没有把握的情况下，选择 L 表时常为该选大表还是选小表而犹豫。假设条件许可，应尽量选用大表，让影响存在的可能性较大的因素和交互作用各占适当的列。*因素或*交互作用的影响是否真的存在，留到方差分析进展显著性检验时再做结论。这样既可以减少试验的工作量，又不致于漏掉重要的信息。5.3.4 正交表的表头设计 所谓表头设计，就是确定试验所考虑的因素和交互作用，在正交表中该放在哪一列的问题。
12、1有交互作用时，表头设计则必须严格地按规定办事。因篇幅限制，此处不讨论，请查阅有关书籍。2 假设试验不考虑交互作用，则表头设计可以是任意的。如在例 5-1 中，对 L 934表头设计，表 5-3 所列的各种方案都是可用的。但是正交表的构造是组合数学问题，必须满足 5.2 中所述的特点。对试验之初不考虑交互作用而选用较大的正交表，空列较多时，最好仍与有交互作用时一样，按规定进展表头设计。只不过将有交互作用的列先视为空列，L 84124 2 水平列的列数为 4 4 水平列的列数为 1 实验的次数 正交表的代号-待试验完毕后再加以判定。54 正交试验的操作方法 1分区组。对于一批试验，如果要使用几台
13、不同的机器，或要使用几种原料来进展，为了防止机器或原料的不同而带来误差，从而干扰试验的分析，可在开场做实验之前，用 L表中未排因素和交互作用的一个空白列来安排机器或原料。与此类似，假设试验指标的检验需要几个人或几台机器来做，为了消除不同人或仪器检验的水平不同给试验分析带来干扰，也可采用在 L 表中用一空白列来安排的方法。这样一种作法叫做分区组法。2 因素水平表排列顺序的随机化。如在例 5-1 中，每个因素的水平序号从小到大时，因素的数值总是按由小到大或由大到小的顺序排列。按正交表做试验时，所有的 1 水平要碰在一起，而这种极端的情况有时是不希望出现的，有时也没有实际意义。因此在排列因素水平表时
14、，最好不要简单地按因素数值由小到大或由大到小的顺序排列。从理论上讲，最好能使用一种叫做随机化的方法。所谓随机化就是采用抽签或查随机数值表的方法，来决定排列的别有顺序。3试验进展的次序没必要完全按照正交表上试验的顺序。为减少试验中由于先后实验操作熟练的程度不匀带来的误差干扰，理论上推荐用抽签的方法来决定试验的次序。4在确定每一个实验的实验条件时，只需考虑所确定的几个因素和分区组该如何取值，而不要其实也无法考虑交互作用列和误差列怎么办的问题。交互作用列和误差列的取值问题由实验本身的客观规律来确定，它们对指标影响的大小在方差分析时给出。5做实验时，要力求严格控制实验条件。这个问题在因素各水平下的数值
15、差异不大时更为重要。例如，例 5-1 中的因素 加碱量m的三个水平：m12.0，m2=2.5，m3=3.0，在以mm2=2.5 为条件的*一个实验中，就必须严格认真地让m2=2.5。假设因为粗心和不负责任，造成m2=2.2 或造成m2=3.0，那就将使整个试验失去正交试验设计方法的特点，使极差和方差分析方法的应用丧失了必要的前提条件，因而得不到正确的试验结果。55 正交试验结果分析方法 正交试验方法之所以能得到科技工作者的重视并在实践中得到广泛的应用，其原因不仅在于能使试验的次数减少，而且能够用相应的方法对试验结果进展分析并引出许多有价值的结论。因此，有正交试验法进展实验，如果不对试验结果进展
16、认真的分析，并引出应该引出的结论，那就失去用正交试验法的意义和价值。5.5.1 极差分析方法 表 53 L 93 4表头设计方案 列 号 1 2 3 4 方 案 1 2 3 4 T 空 m p p T 空 m m p T 空 空 m p T -下面以表 5-4 为例讨论 L423正交试验结果的极差分析方法。极差指的是各列中各水平对应的试验指标平均值的最大值与最小值之差。从表 5-4 的计算结果可知，用极差法分析正交试验结果可引出以下几个结论：1在试验围，各列对试验指标的影响从大到小的排队。*列的极差最大，表示该列的数值在试验围变化时，使试验指标数值的变化最大。所以各列对试验指标的影响从大到小的
17、排队，就是各列极差D的数值从大到小的排队。2试验指标随各因素的变化趋势。为了能更直观地看到变化趋势，常将计算结果绘制成图。3使试验指标最好的适宜的操作条件适宜的因素水平搭配。4可对所得结论和进一步的研究方向进展讨论。5.5.2 方差分析方法 5.5.2.1 计算公式和工程 试验指标的加和值niiy1，试验指标的平均值niiyny11，以第j列为例：j 1水平所对应的试验指标的数值之和 表 54 L423正交试验计算 列 号 1 2 3 试验指标 yi 试验号 1 2 3 n4 1 1 2 2 1 2 1 2 1 2 2 1 y1 y2 y3 y4 j j kj j/kj j/kj 极差Dj 1
18、y1y2 1y3y4 k12 1/k1 1/k1 ma*-min 2y1y3 2y2y4 k22 2/k2 2/k2 ma*-min 3y1y4 3y2y3 k32 3/k3 3/k3 ma*-min 注：j 第j列1水平所对应的试验指标的数值之和；j 第j列2水平所对应的试验指标的数值之和；kj 第j列同一水平出现的次数。等于试验的次数n除以第j列的水平数。j/kj 第j列1水平所对应的试验指标的平均值；j/kj 第j列1水平所对应的试验指标的平均值；Dj 第j列的极差。等于第j列各水平对应的试验指标平均值中的最大值减 最小值，即 Djma*j/kj，j/kj，-min j/kj，j/kj，
19、-j 2水平所对应的试验指标的数值之和 kj 同一水平出现的次数。等于试验的次数除以第j列的水平数 j/kj 1水平所对应的试验指标的平均值 j/kj 1水平所对应的试验指标的平均值 以上 7 项的计算方法同极差法见表 5-4。偏差平方和 fj自由度。fj 第j列的水平数1。Vj方差。Vj Sj fj。Ve误差列的方差。Ve Se fe。式中，e为正交表的误差列。Fj方差之比 Fj Vj Ve。查 F 分布数值表F 分布数值表请查阅有关参考书做显著性检验。总的偏差平方和 niiyyS12总 总的偏差平方和等于各列的偏差平方和之和。即 mjjSS1总 式中，m为正交表的列数。假设误差列由 5 个
20、单列组成，则误差列的偏差平方和Se等于 5 个单列的偏差平方和之和，即：Se Se1 Se2 Se3 Se4 Se5；也可用Se S总 S，来计算，其中S，为安排有因素或交互作用的各列的偏差平方和之和。5.5.2.2 可引出的结论 与极差法相比，方差分析方法可以多引出一个结论：各列对试验指标的影响是否显著，在什么水平上显著。在数理统计上，这是一个很重要的问题。显著性检验强调试验在分析每列对指标影响中所起的作用。如果*列对指标影响不显著，则，讨论试验指标随它的变化趋势是毫无意义的。因为在*列对指标的影响不显著时，即使从表中的数据可以看出该列水平变化时，对应的试验指标的数值与在以*种规律发生变化，
21、但那很可能是由于实验误差所致，将它作为客观规律是不可靠的。有了各列的显著性检验之后，最后应将影响不显著的交互作用列与原来的误差列合并起来。组成新的误差列，重新检验各列的显著性。56 正交试验方法在化工原理实验中的应用举例 例 5-2 为提高真空吸滤装置的生产能力，请用正交试验方法确定恒压过滤的最正确操作条件。其恒压过滤实验的方法、原始数据采集和过滤常数计算等见过滤实验局部。影响实验的主要因素和水平见表 5-5a。表中p为过滤压强差；T为浆液温度；w为浆液质量分数；M为过滤介质材质属多孔瓷。解：1试验指标确实定：恒压过滤常数Km2/s 2选正交表：根据表 5-5a的因素和水平，可选用 L 842
22、4表。3制定实验方案：按选定的正交表，应完成 8 次实验。实验方案见表 5-5b。-4实验结果：将所计算出的恒压过滤常数Km2/s列于表 5-5b。表 5-5a 过滤实验因素和水平 因素 压强差kPa 温度 质量分数 过滤介质 符号 p T w M 水平 1 2 3 4 2.94 3.92 4.90 5.88 室温18 室温1533 稀约 5%浓约 10%G2*G3*G2、G3为过滤漏斗的型号。过滤介质孔径：G2为3050 m、G3为1630 m。表 2-5b正交试验的试验方案和实验结果 列号 j1 2 3 4 5 6 因素 p T w M e Km2/s 试验号 水 平 1 2 3 4 5
23、6 7 8 1 1 2 2 3 3 4 4 1 2 1 2 1 2 1 2 1 2 1 2 2 1 2 1 1 2 2 1 1 2 2 1 1 2 2 1 2 1 1 2 4.01104 2.93104 5.21104 5.55104 4.83104 1.02103 5.11104 1.10103 5指标 K 的极差分析和方差分析：分析结果见表 5-5c。以第 2 列为例说明计算过程：24.011045.211044.831045.111041.92103 22.931045.551041.021031.101032.97103 k24 2/k21.9210-344.79104 2/k22.9
24、710347.42104 D27.42104-4.791042.63104 K4.88103K6.11104 S2k22/k2K2k22/k2K2 44.791046.11104 2 47.421046.11104 2 1.38107 f2第二列的水平数1211 V2S2f21.3810711.38107 SeS5k55/k5K2k55/k5K2 46.221046.11104 2 45.991046.11104 2 1.06109-fef51 VeSefe1.0610911.06109 F2=V2Ve1.381071.06109130.2 查F 分布数值表可知：F 0.01，f11，f214
25、052 F2 F 0.05，f11，f21161.4F2 F 0.10，f11，f2139.9 F2 F 0.25，f11，f215.83 F2 其中：f1为分子的自由度，f2 分母的自由度 所以第二列对试验指标的影响在0.10 水平上显著。其他列的计算结果见表 2-5c。表 55c K 的极差分析和方差分析 列号 j1 2 3 4 5 6 因素 p T w M e Km2/s 工程 j j j j kj j/kj j/kj j/kj j/kj Dj Sj fj Vj Fj F0.01 F0.05 F0.10 F0.25 显著性 6.94104 1.08103 1.50103 1.61103
26、2 3.47104 5.38104 7.52104 8.06103 4.59104 2.65107 3 8.84108 83.6 5403 215.7 53.6 8.20 2*(0.10)1.92103 2.97103 4 4.79104 7.42104 2.63104 1.38107 1 1.38107 130.2 4052 161.4 39.9 5.83 2*(0.10)3.04103 1.84103 4 7.61104 4.61104 3.00104 1.80107 1 1.80107 170.1 4052 161.4 39.9 5.83 3*(0.05)2.54103 2.35103
27、4 6.35104 5.86104 4.85105 4.70109 1 4.70109 4.44 4052 161.4 39.9 5.83 0*(0.25)2.49103 2.40103 4 6.22104 5.99104 2.30105 1.06109 1.06109 1.00 K 4.88103 m2/s K 6.11104 m2/s 6由极差分析结果引出的结论：请同学们自己分析。7由方差分析结果引出的结论。第 1、2 列上的因素 p、T 在0.10 水平上显著；第 3 列上的因素w在0.05水平上显著；第 4 列上的因素M 在0.25 水平上仍不显著。各因素、水平对K的影响变化趋势见图
28、5-3。图 5-3 是用表 5-5a的水平、因素和表-5-5c的j/kj 、j/kj、j/kj、j/k值来标绘的。从图中可看出：A过滤压强差增大，K值增大；B过滤温度增大，K值增大；C过滤浓度增大，K值减小；D过滤介质由 1 水平变为 2 水平，多孔瓷微孔直径减小，K值减小。因为第 4 列对K值的影响在0.25 水平上不显著，所以此变化趋势是不可信的。适宜操作条件确实定。由恒压过滤速率议程式可知，试验指标K值愈大愈好。为此，本例的适宜操作条件是各水平下K的平均值最大时的条件：过滤压强差为 4 水平，5.88kPa 过滤温度为 2 水平，33 过滤浆液浓度为 1 水平，稀滤液 过滤介质为 1 水平或 2 水平这是因为第 4 列对K值的影响在0.25 水平上不显著。为此可优先选择价格廉价或容易得到者。上述条件恰好是正交表中第 8 个试验号。图 53 指标随因素的变化趋势