正交试验设计方法-讲义及举例.docx

《正交试验设计方法-讲义及举例.docx》由会员分享，可在线阅读，更多相关《正交试验设计方法-讲义及举例.docx（10页珍藏版）》请在三一文库上搜索。

1、精品资料可编辑修改正交试验设计方法讲义及举例第5章正交试验设计方法5.1 试验设计方法概述试验设计是数理统计学的一个重要的分支。多数数理统计方法主要用于分析已经得到 的数据，而试验设计却是用于决定数据收集的方法。试验设计方法主要讨论如何合理地安 排试验以及试验所得的数据如何分析等。例5-1 某化工厂想提高某化工产品的质量和产量，对工艺中三个主要因素各按三个水平进行试验（见表 5-1 ）。试验的目的是为提高合格产品的产量，寻求最适宜的操作条件。对此实例该如何进行试验方案水平因素温度C压力Pa加碱量kg符号TPm1(80 )P1(5.0)m 1(2.0)2T2(100)P2(6.0)m2(2.5)

2、3T3(120)P3(7.0)m3(3.0)表5- 1 因素水平的设计呢？很容易想到的是全面搭配法方案（如图5-1所示）：此方案数据点分布的均匀性极好，因素和水平的搭配十分全面，唯一的缺点是实验次数多达 3 3 = 27次（指数3代表3个因素，底数3代表每因素有3个水平）。因素、水平数愈多，则实验次数就愈多，例如，做一个 6因素3水平的试验，就需 36 = 729次实验, 显然难以做到。因此需要寻找一种合适的试验设计方法。试验设计方法常用的术语定义如下。试验指标：指作为试验研究过程的因变量，常为试验结果特征的量 （如得率、纯度等）。例1的试验指标为合格产品的产量。因素：指作试验研究过程的自变量

3、，常常是造成试验指标按某种规律发生变化的那些原因。如例1的温度、压力、碱的用量。水平：指试验中因素所处的具体状态或情况，又称为等级。如例1的温度有3个水平。温度用T表示，下标1、2、3表示因素的不同水平，分别记为 、T2、T3。常用的试验设计方法有：正交试验设计法、均匀试验设计法、单纯形优化法、双水平单纯形优化法、回归正交设计法、序贯试验设计法等。可供选择的试验方法很多，各种试 验设计方法都有其一定的特点。所面对的任务与要解决的问题不同，选择的试验设计方法 也应有所不同。由于篇幅的限制，我们只讨论正交试验设计方法。5. 2 正交试验设计方法的优点和特点用正交表安排多因素试验的方法，称为正交试验

4、设计法。其特点为：完成试验要求 所需的实验次数少。数据点的分布很均匀。可用相应的极差分析方法、方差分析方法、 回归分析方法等对试验结果进行分析，引出许多有价值的结论。从例1可看出，采用全面搭配法方案， 需做27次实验。那么采用简单比较法方案又如何呢？先固定Ti和pi,只改变m,观察因素m不同水平的影响，做了如图 2-2 (1)所示的 三次实验，发现 m = m2时的实验效果最好(好的用口 表示)，合格产品的产量最高，因此认为在后面的实验中因素 m应取m 2水平。固定Ti和m2,改变p的三次实验如图5-2 (2)所示，发现p=p 3时的实验效果最好， 因此认为因素 P应取P 3水平。固定P3和m

5、2,改变T的三次实验如图 5-2 (3)所示，发现因素 T宜取T2水平。因此可以引出结论： 为提高合格产品的产量， 最适宜的操作条件为 T2P 3m2。与全面搭 配法方案相比，简单比较法方案的优点是实验的次数少，只需做9次实验。但必须指出，简单比较法方案的试验结果是不可靠的。因为，在改变m值(或p值，或T值)的三次实验中，说 m2 (或p 3或T2 )水平最好是有条件的。在 Tw，pwpi时，m2水平不 是最好的可能性是有的。在改变m的三次实验中，固定 T =T2, p =p3应该说也是可以的，是随意的，故在此方案中数据点的分布的均匀性是毫无保障的。用这种方法比较 条件好坏时，只是对单个的试验

6、数据进行数值上的简单比较，不能排除必然存在的试验数 据误差的干扰。运用正交试验设计方法，不仅兼有上述两个方案的优点，而且实验次数少，数据点分 布均匀，结论的可靠性较好。正交试验设计方法是用正交表来安排试验的。对于例1适用的正交表是 L9 (34),其试验安排见表5-2。所有的正交表与L9 (34)正交表一样，都具有以下两个特点：(1) 在每一列中，各个不同的数字出现的次数相同。在表 L9 (34)中，每一列有三个水平，水平1、2、3都是各出现3次。(2) 表中任意两列并列在一起形成若干个数字对，不同数字对出现的次数也都相同。在表L9 (34)中，任意两列并列在一起形成的数字对共有9个：(1,1

7、 ), (1,2), (1,3),(2,1 ), (2,2 ) , (2,3 ) , (3,1 ), (3,2 ), (3,3 ),每一个数字对各出现一次。表5-2 试验安排表试验号列号1234因素温度C压力Pa加碱量kg符号TPm11 (Ti)1 (p1)1 (m1)121 (Ti)2(p2)2(m2)231 (Ti)3(p3)3(m3)342 (T2)1 (p1)2 (ml2)352 (T2)2 (p2)3 (m3)162 (T2)3(p3)1(m1)273 (T3)1 (p1)3 (m3)283 (T3)2 (p2)1 (mh)393 (T3)3 (p3)2 (m2)1这两个特点称为正交

8、性。正是由于正交表具有上述特点，就保证了用正交表安排的试 验方案中因素水平是均衡搭配的，数据点的分布是均匀的。因素、水平数愈多，运用正交 试验设计方法，愈发能显示出它白优越性，如上述提到的6因素3水平试验，用全面搭配方案需729次，若用正交表L27 (313)来安排，则只需做 27次试验。在化工生产中， 因素之间常有交互作用。 如果上述的因素 T的数值和水平发生变化 时，试验指标随因素 p变化的规律也发生变化，或反过来，因素 p的数值和水平发生变化 时，试验指标随因素 T变化的规律也发生变化。 这种情况称为因素 T、p间有交互作用，记 为 Tx p。5. 3 正交表使用正交设计方法进行试验方案

9、的设计，就必须用到正交表。正交表请查阅有关参考书。5.1.1 各列水平数均相同的正交表各列水平数均相同的正交表，也称单一水平正交表。这类正交表名称的写法举例如下:L 32各列水平均为 (231)。2 的常用正交表有：L4 (2 3), L8 (27), L12 (211)L16 (215), L20 (219),各列水平数均为3的常用正交表有：L9 (34), L27 (313)。各列水平数均为4的常用正交表有：Li6 (4 5)各列水平数均为3的常用正交表有：L25 (56)5.1.2 混合水平正交表各列水平数不相同的正交表，叫混合水平正交表，下面就是一个混合水平正交表名称 的写法：L 8

10、(41X24)I I | 2水平列的列数为4 4水平列的列数为1|实验的次数|正交表的代号L 8 (41 X24)常简写为L 8 (4X24)。此混合水平正交表含有 1个4水平列，4个2 水平列，共有 1+4 = 5歹U。5.3.3 选择正交表的基本原则一般都是先确定试验的因素、水平和交互作用，后选择适用的L表。在确定因素的水平数时，主要因素宜多安排几个水平，次要因素可少安排几个水平。(1)先看水平数。若各因素全是 2水平，就选用L(2 *)表；若各因素全是 3水平，就 选L(3 *)表。若各因素的水平数不相同，就选择适用的混合水平表。(2)每一个交互作用在正交表中应占一列或二列。要看所选的正

11、交表是否足够大，能 否容纳得下所考虑的因素和交互作用。为了对试验结果进行方差分析或回归分析，还必须 至少留一个空白列，作为“误差”歹U,在极差分析中要作为“其他因素”列处理。(3)要看试验精度的要求。若要求高，则宜取实验次数多的L表。(4)若试验费用很昂贵，或试验的经费很有限，或人力和时间都比较紧张，则不宜选 实验次数太多的 L表。(5)按原来考虑的因素、水平和交互作用去选择正交表，若无正好适用的正交表可选，简便且可行的办法是适当修改原定的水平数。(6)对某因素或某交互作用的影响是否确实存在没有把握的情况下，选择L表时常为该选大表还是选小表而犹豫。若条件许可，应尽量选用大表，让影响存在的可能性

12、较大 的因素和交互作用各占适当的列。某因素或某交互作用的影响是否真的存在，留到方差分 析进行显著性检验时再做结论。这样既可以减少试验的工作量，又不致于漏掉重要的信息。5.3.4 正交表的表头设计所谓表头设计，就是确定试验所考虑的因素和交互作用，在正交表中该放在哪一列的问题。(1)有交互作用时，表头设计则必须严格地按规定办事。因篇幅限制，此处不讨论， 请查阅有关书籍。(2)若试验不考虑交互作用，则表头设计可以是任意的。如在例5-1中，对L 9 (3 4)表头设计，表 5-3所列的各种方案都是可用的。但是正交表的构造是组合数学问题，必须满足5.2中所述的特点。对试验之初不考虑交互作用而选用较大的正

13、交表，空列较多时,表53 L 9(3 4)表头设计方案列号12341TPm空方案2空Tpm3m空Tp4pm空T最好仍与有交互作用时一样，按规定进行表头设计。只不过将有交互作用的列先视为空列， 待试验结束后再加以判定。5. 4 正交试验的操作方法(1)分区组。对于一批试验，如果要使用几台不同的机器，或要使用几种原料来进行，为了防止机器或原料的不同而带来误差，从而干扰试验的分析， 可在开始做实验之前， 用L表中未排因素和交互作用的一个空白列来安排机器或原料。与此类似，若试验指标的检验需要几个人(或几台机器)来做，为了消除不同人(或 仪器)检验的水平不同给试验分析带来干扰，也可采用在L表中用一空白列

14、来安排的办法。这样一种作法叫做分区组法。(2)因素水平表排列顺序的随机化。如在例5-1中，每个因素的水平序号从小到大时，因素的数值总是按由小到大或由大到小的顺序排列。按正交表做试验时，所有的1水平要碰在一起，而这种极端的情况有时是不希望出现的，有时也没有实际意义。因此在排列因 素水平表时，最好不要简单地按因素数值由小到大或由大到小的顺序排列。从理论上讲， 最好能使用一种叫做随机化的方法。所谓随机化就是采用抽签或查随机数值表的办法，来 决定排列的别有顺序。(3)试验进行的次序没必要完全按照正交表上试验号码的顺序。为减少试验中由于先 后实验操作熟练的程度不匀带来的误差干扰，理论上推荐用抽签的办法来

15、决定试验的次序。(4)在确定每一个实验的实验条件时，只需考虑所确定的几个因素和分区组该如何取 值，而不要(其实也无法)考虑交互作用列和误差列怎么办的问题。交互作用列和误差列 的取值问题由实验本身的客观规律来确定，它们对指标影响的大小在方差分析时给出。(5)做实验时，要力求严格控制实验条件。这个问题在因素各水平下的数值差别不大时更为重要。例如，例5-1中的因素(加碱量)m的三个水平：mi = 2.0 , m2=2.5 , m3=3.0 , 在以m = m2=2.5为条件的某一个实验中， 就必须严格认真地让 m2=2.5。若因为粗心和不 负责任，造成 m2=2.2或造成m2=3.0 ,那就将使整个

16、试验失去正交试验设计方法的特点， 使极差和方差分析方法的应用丧失了必要的前提条件，因而得不到正确的试验结果。5. 5正交试验结果分析方法正交试验方法之所以能得到科技工作者的重视并在实践中得到广泛的应用，其原因不 仅在于能使试验的次数减少，而且能够用相应的方法对试验结果进行分析并引出许多有价 值的结论。因此，有正交试验法进行实验，如果不对试验结果进行认真的分析，并引出应该引出的结论，那就失去用正交试验法的意义和价值。5.1.1 极差分析方法下面以表5-4为例讨论L4 (23)正交试验结果的极差分析方法。极差指的是各列中各 水平对应的试验指标平均值的最大值与最小值之差。从表5-4的计算结果可知，用

17、极差法分析正交试验结果可引出以下几个结论：(1)在试验范围内，各列对试验指标的影响从大到小的排队。某列的极差最大，表示 该列的数值在试验范围内变化时，使试验指标数值的变化最大。所以各列对试验指标的影 响从大到小的排队，就是各列极差 D的数值从大到小的排队。(2)试验指标随各因素的变化趋势。为了能更直观地看到变化趋势，常将计算结果绘 制成图。(3)使试验指标最好的适宜的操作条件(适宜的因素水平搭配)。(4)可对所得结论和进一步的研究方向进行讨论。表5-4 L4 (23)正交试验计算歹1J号123试验指标yi1111y1试验号2122y23212y3n= 4221y4I ji 1=y1 + y21

18、2=y1 + y313=y1 + y4口 jn 1= y3 + y4n 2 = y2 + y4n 3= y2+ y3kjk1=2k2 = 2k3 = 2I j/ kjI 1/ k1I 2/ k2I 3/ k3 j/ kjn 1/ k1n 2/ k2n 3/ k3极差(Dj)max -min max -min max -min 注:Ij 第j列“1”水平所对应的试验指标的数值之和；L 第j列“2”水平所对应的试验指标的数值之和；K 第j列同一水平出现的次数。等于试验的次数(n)除以第j列的水平数I j/ kj 第j列“ 1”水平所对应的试验指标的平均值；口 j/ kj 第j列“ 1 ”水平所对应

19、的试验指标的平均值；Dj 第j列的极差。等于第j列各水平对应的试验指标平均值中的最大值减最小值，即Dj=max I j/ kj , 口 j/ kj ,-min j kj , U j/ kj ,5.1.2 方差分析方法5.1.2.1 计算公式和项目试验指标的加和值=yi ,以第j列为例:i，一 1yi ,试验指标的平均值 y -n .iI j “1”水平所对应的试验指标的数值之和Ilj “2”水平所对应的试验指标的数值之和 kj同一水平出现的次数。等于试验的次数除以第j列的水平数I j/ kj “1 ”水平所对应的试验指标的平均值 n j/ kj “1”水平所对应的试验指标的平均值 偏差平方和S

20、j以上7项的计算方法同极差法（见表 5-4 ）。22Ij 一II j III jkj y k j y k j kjkjkjfj 自由度。fj =第j列的水平数一1。(10)Vj方差。Vj = Sj / fj （11） Ve 误差列的方差。 Ve =Se /fe。式中，e为正交表的误差列。Fj 方差之比Fj =Vj /Ve。 查F分布数值表（F分布数值表请查阅有关参考书）做显著性检验。 n（14）总的偏差平方和s总1 y2i 1 m（15）总的偏差平方和等于各列的偏差平方和之和。即 5总= Sj j 1式中，m为正交表的列数。若误差列由5个单列组成，则误差列的偏差平方和 Se等于5个单列的偏差平

21、方和之和， 即：Se =Se1 +Se2 +Se3 +Se4 + Se5 ；也可用Se =S总+ S，来计算，其中 S，为安排有因素或交互作用的各列的偏差平方和之和。5.5.2.2 可引出的结论与极差法相比，方差分析方法可以多引出一个结论：各列对试验指标的影响是否显著， 在什么水平上显著。在数理统计上，这是一个很重要的问题。显著性检验强调试验在分析 每列对指标影响中所起的作用。如果某列对指标影响不显著，那么，讨论试验指标随它的 变化趋势是毫无意义的。因为在某列对指标的影响不显著时，即使从表中的数据可以看出 该列水平变化时，对应的试验指标的数值与在以某种“规律”发生变化，但那很可能是由 于实验误

22、差所致，将它作为客观规律是不可靠的。有了各列的显著性检验之后，最后应将影响不显著的交互作用列与原来的“误差列”合并起来。组成新的“误差列”，重新检验各列的显著性。5. 6 正交试验方法在化工原理实验中的应用举例例5-2 为提高真空吸滤装置的生产能力，请用正交试验方法确定恒压过滤的最佳操作 条件。其恒压过滤实验的方法、原始数据采集和过滤常数计算等见过滤实验部分。影 响实验的主要因素和水平见表5-5 ( a)。表中A p为过滤压强差；T为浆液温度；w为浆液质量分数；M为过滤介质(材质属多孔陶瓷)。解：(1)试验指标的确定：恒压过滤常数K (m2/s )(2)选正交表：根据表 5-5 (a)的因素和

23、水平，可选用 L 8 (4X24)表。(3)制定实验方案：按选定的正交表，应完成8次实验。实验方案见表 5-5 (b)。(4)实验结果：将所计算出的恒压过滤常数K (m2/s )列于表5-5 (b)。分析结果见表5-5 (c)。以第2列为例说明计算过程:12 =4.01 X 10112 = 2.93 X 104 + 5.2144 + 5.55X 10 4+4.83X 10 4+1.02X 10x 104 + 5.113+1.10X 10 4X 10 3= 1.92 X10 3= 2.97 X 10 3表5-5 (a)过滤实验因素和水平因素压强差/ kPaP温度/cT质量分数w过滤介质M符号12

24、.94(室温)18稀(约5% )G2*23.92(室温+15) 33浓(约10% )G3*水平34.9045.88* G2、G3为过滤漏斗的型号。过滤介质孔径：G2为3050 “m、G3为1630 gm表2-5 (b)正交试验的试验方案和实验结果列号j=123456因素 PTwMeK (m2/s )试验号水平111111一一 一44.01 X 10 421222242.93 X 10 432112245.21 X 10 442221145.55 X 10 453121244.83 X 10 463212131.02 X 10 374122145.11 X 10 4842112一一 31.10

25、X 10 3(5)指标K的极差分析和方差分析:k2 = 4I 2/ k2= 1.92X 10-3/4 = 4.79 X 10n 2/ k2=2.97D2 = 7.42 X 10X 10 3/4= 7.42 X 104 - 4.79 X 10 4=2.63X10 410-72 K = 4.88 X 10 3 K 6.11 X 10 4S2 = k 2 (I 2/ k 2 K)+k2 ( n 2 / k 2 K)=4 (4.79 X 10 4-6.11 X 10 4 ) 2 +4 (7.42 X 10 4-6.11 X 10 4 ) 2= 1.38 XX 10f2=第二列的水平数一1 =2 1 =

26、 1V2=S2/f2= 1.38 x 10 7/1=1.38 x 10 7Se=S5=k5 (I 5/ k 5 K ) 2 + k5 (II5/ k 5 K ) 2=4 (6.22 X 10 4-6.11 X 10 4 ) 2 +4 (5.99 X 10 4-6.11 X 10 4 ) 2 = 1.06fe= f5 = 1Ve=Se/ fe= 1.06 X 10 9/1 = 1.06 X 10 9F2 =V2/Ve= 1.38 X 10 7/1.06 X 10 9= 130.2查F分布数值表可知：F( a= 0.01,f1 = 1 ,f2= 1 )= 4052 F2F( a= 0.05,f =

27、 1,f2=1)= 161.4F2F(a= 0.10,f=1,f2=1)= 39.9 F2F( a= 0.25,f1=1,fz=1)= 5.83 F2(其中：f1为分子的自由度，f2分母的自由度)所以第二列对试验指标的影响在=0.10水平上显著。其他列白计算结果见表2-5 (c)。表55 (c) K的极差分析和方差分析列号j=123456因素 pTwMeK (m2/s )项目I j_433一336.94 X 10 41.92 X 10 33.04 X 10 3 2.54 X 10 3 2.49 X 10 32 K =口 j3333_31.08 X 10 32.97 X 10 31.84 X 1

28、0 3 2.35 X 10 3 2.40 X 10 34.88 X 10Wj31.50 X 10 3-3IVj一31.61 X 10 3(m2/s )kj24444i j/ kj_4444_43.47 X 10 44.79 X 10 47.61 X 10 46.35 X 10 46.22 X 10 4nj/ kj_444445.38 X 10 4 7.42 X 10 4 4.61 X 10 45.86 X 10 45.99 X 10 4叫/ kj _ 47.52 X 10 4IVj/ kj _ 38.06 X 10 3KDj444554.59 X 10 42.63 X 10 43.00 X 1

29、0 4 4.85 X 10 5 2.30 X 10 56.11 X 10Sj _ _7_ _ 7_ _ 7 _ _ 9_ _92.65 X 101.38 X 101.80 X 104.70 X 101.06 X 10-4fj3111(m2/s )Vj877998.84 X 10 81.38 X 101.80 X 104.70 X 101.06 X 10Fj83.6130.2170.14.441.00F0.015403405240524052Fg.05215.7161.4161.4161.4F0.1053.639.939.939.9Fg.258.205.835.835.83显著性2*(0.10)

30、2*(0.10)3*(0.05)0*(0.25)(6)由极差分析结果引出的结论：请同学们自己分析。(7)由方差分析结果引出的结论。 第1、2列上的因素 A p、T在 =0.10水平上显著；第3列上的因素 w在 = 0.05水平上显著；第 4列上的因素 M在 =0.25水平上仍不显著。 各因素、水平对 K的影响变化趋势见图 5-3。图5-3是用表5-5 (a)的水平、因 素和表5-5 (c)的I j/ kj 、Hj/ kj、nij/ kj、iVj/ k值来标绘的。从图中可看出：A.过滤压强差增大，K值增大；B.过滤温度增大，K值增大；C.过滤浓度增大，K值减小；D.过滤介质由1水平变为2水平，多孔陶瓷微孔直径减小，K值减小。因为第4列对K值的影响在=0.25水平上不显著，所以此变化趋势是不可信的。不rQIxx图5 3指标随因素的变化趋势K值愈大愈好。为适宜操作条件的确定。由恒压过滤速率议程式可知，试验指标 此，本例的适宜操作条件是各水平下K的平均值最大时的条件：过滤压强差为4水平，5.88kPa过滤温度为2水平，33 C过滤浆液浓度为1水平，稀滤液= 0.25水平上不显过滤介质为1水平或2水平(这是因为第 4列对K值的影响在 著。为此可优先选择价格便宜或容易得到者)。上述条件恰好是正交表中第8个试验号。