为什么要用正交试验
我们知道如果有很多的因素变化制约着一个事件的变化,那么为了弄明白哪些因素重要,哪些不重要,什么样的因素搭配会产生极值,必须通过做实验验证(仿真也可以说是试验,只不过试验设备是计算机),如果因素很多,而且每种因素又有多种变化(专业称法是:水平),那么试验量会非常的大,显然是不可能每一个试验都做的。
就影响主轴温升的试验来讲,影响主轴温升的因素很多,比如转速、预紧力、油气压力、喷油间隙时间、油品等等;每种因素的水平也很多,比如转速从8Krpm到20Krpm,等等,计算一下,所有因素都做,大概一共要900次试验,按一天3次试验计,要不停歇的做10个月,显然是不可能的。
能够大幅度减少试验次数而且并不会降低试验可行度的方法就是使用正交试验法。
首先需要选择一张和你的试验因素水平相对应的正交表,已经有数学家制好了很多相应的表,你只需找到对应你需要的就可以了。
所谓正交表,也就是一套经过周密计算得出的现成的试验方案,他告诉你每次试验时,用那几个水平互相匹配进行试验,这套方案的总试验次数是远小于每种情况都考虑后的试验次数的。
比如3水平4因素表就只有9行,远小于遍历试验的81次;我们同理可推算出如果因素水平越多,试验的精简程度会越高。
正交试验设计介绍
当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。
因此就出现了分式析因设计(fractional factorial de signs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。
正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。
是一种高效率、快速、经济的实验设计方法。
日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。
例如作一个三因素三水平的实验,按全面实验要求,须进行3³=27种组合的实验,且尚未考虑每一组合的重复数。
若按L9(3)³正交表按排实验,只需作9次,7
按L18(3) 正交表进行18次实验,显然大大减少了工作量。
因而正交实验设计在很多领域的研究中已经得到广泛应用。
1.正交表
正交表是一整套规则的设计表格。
L 为正交表的代号,n 为试验的次数,t 为水平数,c 为列数,也就是可能安排最多的因素个数。
例如L 9(3 ) , (表11),它表示需作9次实验,
最多可观察4个因素,每个因素均为3水平。
一个正交表中也可以各列的水平数不相等,我
们称它为混合型正交表,如L 8(4×2 ) (表12),此表的5列中,有1列为4水平,4列为2
水平。
根据正交表的数据结构看出,正交表是一个n 行c 列的表,其中第j 列由数码1,2,… Sj 组成,这些数码均各出现N/S 次,例如表11中,第二列的数码个数为3,S=3 ,即由1、2、3组成,各数码均出现1次。
正交表具有以下两项性质:
(1) 每一列中,不同的数字出现的次数相等。
例如在两水平正交表中,任何一列都有数码“1”与“2”,且任何一列中它们出现的次数是相等的;如在三水平正交表中,任何一列都有“1”、“2”、“3”,且在任一列的出现数均相等。
(2) 任意两列中数字的排列方式齐全而且均衡。
例如在两水平正交表中,任何两列(同一横行内)有序对子共有4种:(1,1)、(1,2)、(2,1)、(2,2)。
每种对数出现次数相等。
在三水平情况下,任何两列(同一横行内)有序对共有9种,1.1、1.2、1.3、2.1、2.2、2.3、3.1、
3.2、3.3,且每对出现数也均相等。
以上两点充分的体现了正交表的两大优越性,即“均匀分散性,整齐可比”。
通俗的说,每个因素的每个水平与另一个因素各水平各碰一次,这就是正交性。
4 4
表1 L 4 (2 ) 正交表
2. 交互作用表 每一张正交表后都附有相应的交互作用表,它是专门用来安排交互作用试验。
表14就是L 8(2 )表的交互作用表。
安排交互作用的试验时,是将两个因素的交互作用当作一个新的因素,占用一列,为交互作用列,从表14中可查出L 8(2 )正交表中的任何两列的交互作用列。
表中带( )的为主因素的列号,它与另一主因素的交互列为第一个列号从左向右,第二个列号顺次由下向上,二者相交的号为二者的交互作用列。
例如将A 因素排为第(1)列,B 因素排为第(2)列,两数字相交为3,则第3列为A×B 交互作用列。
又如可以看到第4列与第6列的交互列是第2列,等等。
3.正交实验的表头设计
表头设计是正交设计的关键,它承担着将各因素及交互作用合理安排到正交表的各列中的重要任务,因此一个表头设计就是一个设计方案。
表头设计的主要步骤如下:
(1) 确定列数
根据试验目的,选择处理因素与不可忽略的交互作用,明确其共有多少个数,如果对研究中的某些问题尚不太了解,列可多一些,但一般不宜过多。
当每个试验号无重复,只有1个试验数据时,可设2个或多个空白列,作为计算误差项之用。
(2) 确定各因素的水平数
根据研究目的,一般二水平(有、无)可作因素筛选用。
也可适用于试验次数少、分批进行的研究。
三水平可观察变化趋势,选择最佳搭配,多水平能以一次满足试验要求。
(3) 选定正交表
根据确定的列数c ,与水平数(t)选择相应的正交表。
例如观察5个因素8个一级交互作用,留两个空白列,且每个因素取2水平,则适宜选L 16(2 )表。
由于同水平的正交表有多
个,如L 8(2 )、L 12(2 )、L 16(2 ),一般只要表中列数比考虑需要观察的个数稍多一点即可,这样省工省时。
(4) 表头安排
应优先考虑交互作用不可忽略的处理因素,按照不可混杂的原则,将它们及交互作用首先在表头排妥,而后再将剩余各因素任意安排在各列上。
例如某项目考察4个因素A 、B 、C 、7 7 7 11 15 15
3
D 及A×B 交互作用,各因素均为2水平,现选取L 8(2 )表,由于A 、B 两因素需要观察其交互作用,故将二者优先安排在第1、2列,根据交互作用表查得A×B 应排在第3列,于是C 排在第4列,由于A×C 交互在第5列,B×C 交互作用在第6列,虽然未考查A×C 与B×C ,为避免混杂之嫌,D 就排在第7列。
(5) 组织实施方案
根据选定正交表中各因素占有列的水平数列,构成实施方案表,按实验号依次进行,共作n 次实验,每次实验按表中横行的各水平组合进行。
例如L 9(3 )表,若安排四个因素,第一次实验A 、B 、C 、D 四因素均取1水平,第二次实验A 因素1水平,B 、C 、D 取2水平,……第九次实验A 、B 因素取3水平,C 因素取2水平,D 因素取1水平。
实验结果数据记录在该行的末尾。
因此整个设计过程我们可用一句话归纳为:“因素顺序上列、水平对号入座,实验横着作”。
4.二水平有交互作用的正交实验设计与方差分析
例8 某研究室研究影响某试剂回收率的三个因素,包括温度、反应时间、原料配比,每个因素都为二水平,各因素及其水平见表16。
选用L 8(2 )正交表进行实验,实验结果见表 17。
首先计算Ij 与IIj ,Ij 为第j 列第1水平各试验结果取值之和,IIj 为第j 列第2水平各试验结果取值之和。
然后进行方差分析。
过程为:
求:总离差平方和
各列离差平方和 SSj=
本例各列离均差平方和见表10最底部一行。
即各空列SSj 之和。
即误差平方和 自由度v 为各列水平数减1,交互作用项的自由度为相交因素自由度的乘积。
分析结果见表18。
从表18看出,在α=0.05水准上,只有C 因素与A×B 交互作用有统计学意义,其余各因素均无统计学意义,A 因素影响最小,考虑到交互作用A×B 的影响较大,且它们的二水平为优。
在C 2的情况下, 有B 1、A 2和B 1、A 1两种组合状况下的回收率最高。
考虑到B 因素影响较A 因素影响大一些,而B 中选B 1为好,故选A 2、B 1。
这样最后决定最佳配方为A 2、B 1、C 2,即80℃,反应时间2.5h ,原料配比为1.2:1。
如果使用计算机进行统计分析,在数据是只需要输入试验因素和实验结果的内容,交互作用界的内容不用输入,然后按照表头定义要分析的模型进行方差分析。
4 7。