在线正交试验设计网站-正交实验

正交实验怎么做?

正交试验设计法,就是使用已经造好了的表格--正交表--来安排试验并进行资料分析的一种方法。它简单易行,计算表格化,使用者能够迅速掌握。下边通过一个例子来说明正交试验设计法的基本思想。

[例1]为提高某化工产品的转化率,选择了三个有关因素进行条件试验,反应温度(A),反应时间(B),用碱量(C),并确定了它们的试验范围:

A:80-90℃

B:90-150分钟

C:5-7%

试验目的是搞清楚因子A、B、C对转化率有什么影响,哪些是主要的,哪些是次要的,从而确定最适生产条件,即温度、时间及用碱量各为多少才能使转化率高。试制定试验方案。

这里,对因子A,在试验范围内选了三个水平;因子B和C也都取三个水平:

A:A1=80℃,A2=85℃,A3=90℃

B:B1=90分,B2=120分,B3=150分

C:C1=5%,C2=6%,C3=7%

当然,在正交试验设计中,因子可以是定量的,也可以是定性的。而定量因子各水平间的距离可以相等,也可以不相等。

正交实验是怎么做的

我的实验是三因素三水平,做正交的话就是9次实验(我做一次实验得出一组资料就是一个月啊)。 我看不少论文上都是分两大步骤: 一、正交设计 1、进行正交设计,设计出后续的实验方案。 2、根据正交设计的实验方案进行试验(不写出过程,不做图表,因为各单因素的搭配是随机的,没法画曲线图分析),只给出结果。 3、根据实验结果,进行正交分析,得出各单因素的影响大小。 (在正交设计这部分,我就得做九个实验啊) 二、单因素实验 对正交设计的分析结果进行验证所补充的实验,这次所做的实验具有可比性,可以对单个因素影响效果做曲线图分析。 如此一来,那我的实验岂不是就做了两遍(虽然这两遍实验的控制因素水平可能不一样)。 做这么多实验就啥时候毕业呀。 我的问题是: 有人说先单因素方便确定正交设计的因素水平。 2 ,我不想做正交设计。想直接用简单比较法得出最优实验方案算了,不过这样做出的论文不够炫,感觉没有技术含量,和人家做正交设计的相比就差了一大截。 很是头疼,定不下来方案。希望各位走过路过的朋友帮帮忙,

正交实验怎么设计?

正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。例如作一个三因素三水平的实验,按全面实验要求,须进行3的3次方=27 种组合的实验,且尚未考虑每一组合的重复数。若按L9(3)3 正交表按排实验,只需作9 次,按L18(3)7 正交表进行18 次实验,显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。(汗,这里不能打出来正确的表达,反正学这个的都知道具体的写法)

正交表是一整套规则的设计表格,L 为正交表的代号,n 为试验的次数,t为水平数,c 为列数,也就是可能安排最多的因素个数。例如L9(34),它表示需作9次实验,最多可观察4 个因素,每个因素均为3 水平。一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表,如L8(4×24) ,此表的5 列中有1 列为4 水平,4 列为2水平。根据正交表的资料结构看出,正交表是一个n 行c 列的表,其中第j 列由数码1,2,… Sj 组成,这些数码均各出现N/S 次,例如表11 中,第二列的数码个数为3,S=3 ,即由1、2、3 组成,各数码均出现N/3=9/3=3次。

正交实验

正交实验设计

当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。因此就出现了分式析因设计(fractional factorial designs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。

正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。例如作一个三因素三水平的实验,按全面实验要求,须进行33=27种组合的实验,且尚未考虑每一组合的重复数。若按L9(3)3正交表按排实验,只需作9次,按L18(3)7正交表进行18次实验,显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。

1.正交表

正交表是一整套规则的设计表格,用 。L为正交表的代号,n为试验的次数,t为水平数,c为列数,也就是可能安排最多的因素个数。例如L9(34), (表11),它表示需作9次实验,最多可观察4个因素,每个因素均为3水平。一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表,如L8(4×24) (表12),此表的5列中,有1列为4水平,4列为2水平。根据正交表的资料结构看出,正交表是一个n行c列的表,其中第j列由数码1,2,… Sj 组成,这些数码均各出现N/S 次,例如表11中,第二列的数码个数为3,S=3 ,即由1、2、3组成,各数码均出现 次。

正交表具有以下两项性质:

(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,且每对出现数也均相等。

以上两点充分的体现了正交表的两大优越性,即“均匀分散性,整齐可比”。通俗的说,每个因素的每个水平与另一个因素各水平各碰一次,这就是正交性。

2. 互动作用表 每一张正交表后都附有相应的互动作用表,它是专门用来安排互动作用试验。表14就是L8(27)表的互动作用表。

安排互动作用的试验时,是将两个因素的互动作用当作一个新的因素,占用一列,为互动作用列,从表14中可查出L8(27)正交表中的任何两列的互动作用列。表中带( )的为主因素的列号,它与另一主因素的互动列为第一个列号从左向右,第二个列号顺次由下向上,二者相交的号为二者的互动作用列。例如将A因素排为第(1)列,B因素排为第(2)列,两数字相交为3,则第3列为A×B互动作用列。又如可以看到第4列与第6列的互动列是第2列,等等。

3.正交实验的表头设计 表头设计是正交设计的关键,它承担著将各因素及互动作用合理安排到正交表的各列中的重要任务,因此一个表头设计就是一个设计方案。

表头设计的主要步骤如下:

(1)确定列数 根据试验目的,选择处理因素与不可忽略的互动作用,明确其共有多少个数,如果对研究中的某些问题尚不太了解,列可多一些,但一般不宜过多。当每个试验号无重复,只有1个试验资料时,可设2个或多个空白列,作为计算误差项之用。

(2)确定各因素的水平数 根据研究目的,一般二水平(有、无)可作因素筛选用;也可适用于试验次数少、分批进行的研究。三水平可观察变化趋势,选择最佳搭配;多水平能以一次满足试验要求。

(3)选定正交表 根据确定的列数(c)与水平数(t)选择相应的正交表。例如观察5个因素8个一级互动作用,留两个空白列,且每个因素取2水平,则适宜选L16(215)表。由于同水平的正交表有多个,如L8(27)、L12(211)、L16(215),一般只要表中列数比考虑需要观察的个数稍多一点即可,这样省工省时。

(4)表头安排 应优先考虑互动作用不可忽略的处理因素,按照不可混杂的原则,将它们及互动作用首先在表头排妥,而后再将剩余各因素任意安排在各列上。例如某专案考察4个因素A、B、C、D及A×B互动作用,各因素均为2水平,现选取L8(27)表,由于AB两因素需要观察其互动作用,故将二者优先安排在第1、2列,根据互动作用表查得A×B应排在第3列,于是C排在第4列,由于A×C互动在第5列,B×C互动作用在第6列,虽然未考查A×C与B×C,为避免混杂之嫌,D就排在第7列。

(5)组织实施方案 根据选定正交表中各因素占有列的水平数列,构成实施方案表,按实验号依次进行,共作n次实验,每次实验按表中横行的各水平组合进行。例如L9(34)表,若安排四个因素,第一次实验A、B、C、D四因素均取1水平,第二次实验A因素1水平,B、C、D取2水平,……第九次实验A、B因素取3水平,C因素取2水平,D因素取1水平。实验结果资料记录在该行的末尾。因此整个设计过程我们可用一句话归纳为:“因素顺序上列、水平对号入座,实验横著作”。

4.二水平有互动作用的正交实验设计与方差分析

例8 某研究室研究影响某试剂回收率的三个因素,包括温度、反应时间、原料配比,每个因素都为二水平,各因素及其水平见表16。选用L8(27)正交表进行实验,实验结果见表17。

首先计算Ij 与IIj ,Ij为第j列第1水平各试验结果取值之和,IIj为第j列第2水平各试验结果取值之和。然后进行方差分析。过程为:

求:总离差平方和

各列离差平方和 SSj=

本例各列离均差平方和见表10最底部一行。即各空列SSj之和。即误差平方和

自由度v为各列水平数减1,互动作用项的自由度为相交因素自由度的乘积。

分析结果见表18。

从表18看出,在α=0.05水准上,只有C因素与A×B互动作用有统计学意义,其余各因素均无统计学意义,A因素影响最小,考虑到互动作用A×B的影响较大,且它们的二水平为优。在C2的情况下, 有B1A2和B1,A1两种组合状况下的回收率最高。考虑到B因素影响较A因素影响大些,而B中选B1为好,故选A2B1。这样最后决定最佳配方为A2B1C2,即80℃,反应时间2.5h,原料配比为1.2:1。

如果使用计算机进行统计分析,在资料是只需要输入试验因素和实验结果的内容,互动作用界的内容不用输入,然后按照表头定义要分析的模型进行方差分析。

帮忙 做个正交实验表

L44 正交实验表

所在列 1 2 3 4

因素 A B C D

实验01 1 1 1 1

实验02 1 2 2 2

实验03 1 3 3 3

实验04 1 4 4 4

实验05 2 1 2 3

实验06 2 2 1 4

实验07 2 3 4 1

实验08 2 4 3 2

实验09 3 1 3 4

实验10 3 2 4 3

实验11 3 3 1 2

实验12 3 4 2 1

实验13 4 1 4 2

实验14 4 2 3 1

实验15 4 3 2 4

实验16 4 4 1 3

L46 实验计划表

所在列 1 2 3 4 5 6

因素 A B C D E F

实验01 1 1 1 1 1 1

实验02 1 2 2 2 2 2

实验03 1 3 3 3 3 3

实验04 1 4 4 4 4 4

实验05 2 1 1 2 2 3

实验06 2 2 2 1 1 4

实验07 2 3 3 4 4 1

实验08 2 4 4 3 3 2

实验09 3 1 2 3 4 1

实验10 3 2 1 4 3 2

实验11 3 3 4 1 2 3

实验12 3 4 3 2 1 4

实验13 4 1 2 4 3 3

实验14 4 2 1 3 4 4

实验15 4 3 4 2 1 1

实验16 4 4 3 1 2 2

实验17 1 1 4 1 4 2

实验18 1 2 3 2 3 1

实验19 1 3 2 3 2 4

实验20 1 4 1 4 1 3

实验21 2 1 4 2 3 4

实验22 2 2 3 1 4 3

实验23 2 3 2 4 1 2

实验24 2 4 1 3 2 1

实验25 3 1 3 3 1 2

实验26 3 2 4 4 2 1

实验27 3 3 1 1 3 4

实验28 3 4 2 2 4 3

实验29 4 1 3 4 2 4

实验30 4 2 4 3 1 3

实验31 4 3 1 2 4 2

实验32 4 4 2 1 3 1

求正交实验表

L44

正交实验表

所在列

1 2 3 4

因素

A B C D

实验01

1 1 1 1

实验02

1 2 2 2

实验03

1 3 3 3

实验04

1 4 4 4

实验05

2 1 2 3

实验06

2 2 1 4

实验07

2 3 4 1

实验08

2 4 3 2

实验09

3 1 3 4

实验10

3 2 4 3

实验11

3 3 1 2

实验12

3 4 2 1

实验13

4 1 4 2

实验14

4 2 3 1

实验15

4 3 2 4

实验16

4 4 1 3

L46 实验计划表

所在列

1 2 3 4 5 6

因素

A B C D E F

实验01

1 1 1 1 1 1

实验02

1 2 2 2 2 2

实验03

1 3 3 3 3 3

实验04

1 4 4 4 4 4

实验05

2 1 1 2 2 3

实验06

2 2 2 1 1 4

实验07

2 3 3 4 4 1

实验08

2 4 4 3 3 2

实验09

3 1 2 3 4 1

实验10

3 2 1 4 3 2

实验11

3 3 4 1 2 3

实验12

3 4 3 2 1 4

实验13

4 1 2 4 3 3

实验14

4 2 1 3 4 4

实验15

4 3 4 2 1 1

实验16

4 4 3 1 2 2

实验17

1 1 4 1 4 2

实验18

1 2 3 2 3 1

实验19

1 3 2 3 2 4

实验20

1 4 1 4 1 3

实验21

2 1 4 2 3 4

实验22

2 2 3 1 4 3

实验23

2 3 2 4 1 2

实验24

2 4 1 3 2 1

实验25

3 1 3 3 1 2

实验26

3 2 4 4 2 1

实验27

3 3 1 1 3 4

实验28

3 4 2 2 4 3

实验29

4 1 3 4 2 4

实验30

4 2 4 3 1 3

实验31

4 3 1 2 4 2

实验32

4 4 2 1 3 1

以上仅供参考

希望能解决您的问题。

如何设计正交实验?通过正交实验如何优化实验结果?

正交实验设计

当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。因此就出现了分式析因设计(fractional factorial designs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。

正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。例如作一个三因素三水平的实验,按全面实验要求,须进行33=27种组合的实验,且尚未考虑每一组合的重复数。若按L9(3)3正交表按排实验,只需作9次,按L18(3)7正交表进行18次实验,显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。

1.正交表

正交表是一整套规则的设计表格,用 。L为正交表的代号,n为试验的次数,t为水平数,c为列数,也就是可能安排最多的因素个数。例如L9(34), (表11),它表示需作9次实验,最多可观察4个因素,每个因素均为3水平。一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表,如L8(4×24) (表12),此表的5列中,有1列为4水平,4列为2水平。根据正交表的资料结构看出,正交表是一个n行c列的表,其中第j列由数码1,2,… Sj 组成,这些数码均各出现N/S 次,例如表11中,第二列的数码个数为3,S=3 ,即由1、2、3组成,各数码均出现 次。

正交表具有以下两项性质:

(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,且每对出现数也均相等。

以上两点充分的体现了正交表的两大优越性,即“均匀分散性,整齐可比”。通俗的说,每个因素的每个水平与另一个因素各水平各碰一次,这就是正交性。

2. 互动作用表 每一张正交表后都附有相应的互动作用表,它是专门用来安排互动作用试验。表14就是L8(27)表的互动作用表。

安排互动作用的试验时,是将两个因素的互动作用当作一个新的因素,占用一列,为互动作用列,从表14中可查出L8(27)正交表中的任何两列的互动作用列。表中带( )的为主因素的列号,它与另一主因素的互动列为第一个列号从左向右,第二个列号顺次由下向上,二者相交的号为二者的互动作用列。例如将A因素排为第(1)列,B因素排为第(2)列,两数字相交为3,则第3列为A×B互动作用列。又如可以看到第4列与第6列的互动列是第2列,等等。

3.正交实验的表头设计 表头设计是正交设计的关键,它承担著将各因素及互动作用合理安排到正交表的各列中的重要任务,因此一个表头设计就是一个设计方案。

表头设计的主要步骤如下:

(1)确定列数 根据试验目的,选择处理因素与不可忽略的互动作用,明确其共有多少个数,如果对研究中的某些问题尚不太了解,列可多一些,但一般不宜过多。当每个试验号无重复,只有1个试验资料时,可设2个或多个空白列,作为计算误差项之用。

(2)确定各因素的水平数 根据研究目的,一般二水平(有、无)可作因素筛选用;也可适用于试验次数少、分批进行的研究。三水平可观察变化趋势,选择最佳搭配;多水平能以一次满足试验要求。

(3)选定正交表 根据确定的列数©与水平数(t)选择相应的正交表。例如观察5个因素8个一级互动作用,留两个空白列,且每个因素取2水平,则适宜选L16(215)表。由于同水平的正交表有多个,如L8(27)、L12(211)、L16(215),一般只要表中列数比考虑需要观察的个数稍多一点即可,这样省工省时。

(4)表头安排 应优先考虑互动作用不可忽略的处理因素,按照不可混杂的原则,将它们及互动作用首先在表头排妥,而后再将剩余各因素任意安排在各列上。例如某专案考察4个因素A、B、C、D及A×B互动作用,各因素均为2水平,现选取L8(27)表,由于AB两因素需要观察其互动作用,故将二者优先安排在第1、2列,根据互动作用表查得A×B应排在第3列,于是C排在第4列,由于A×C互动在第5列,B×C互动作用在第6列,虽然未考查A×C与B×C,为避免混杂之嫌,D就排在第7列。

(5)组织实施方案 根据选定正交表中各因素占有列的水平数列,构成实施方案表,按实验号依次进行,共作n次实验,每次实验按表中横行的各水平组合进行。例如L9(34)表,若安排四个因素,第一次实验A、B、C、D四因素均取1水平,第二次实验A因素1水平,B、C、D取2水平,……第九次实验A、B因素取3水平,C因素取2水平,D因素取1水平。实验结果资料记录在该行的末尾。因此整个设计过程我们可用一句话归纳为:“因素顺序上列、水平对号入座,实验横著作”。

4.二水平有互动作用的正交实验设计与方差分析

例8 某研究室研究影响某试剂回收率的三个因素,包括温度、反应时间、原料配比,每个因素都为二水平,各因素及其水平见表16。选用L8(27)正交表进行实验,实验结果见表17。

首先计算Ij 与IIj ,Ij为第j列第1水平各试验结果取值之和,IIj为第j列第2水平各试验结果取值之和。然后进行方差分析。过程为:

求:总离差平方和

各列离差平方和 SSj=

本例各列离均差平方和见表10最底部一行。即各空列SSj之和。即误差平方和

自由度v为各列水平数减1,互动作用项的自由度为相交因素自由度的乘积。

分析结果见表18。

从表18看出,在α=0.05水准上,只有C因素与A×B互动作用有统计学意义,其余各因素均无统计学意义,A因素影响最小,考虑到互动作用A×B的影响较大,且它们的二水平为优。在C2的情况下, 有B1A2和B1,A1两种组合状况下的回收率最高。考虑到B因素影响较A因素影响大些,而B中选B1为好,故选A2B1。这样最后决定最佳配方为A2B1C2,即80℃,反应时间2.5h,原料配比为1.2:1。

如果使用计算机进行统计分析,在资料是只需要输入试验因素和实验结果的内容,互动作用界的内容不用输入,然后按照表头定义要分析的模型进行方差分析。

正交实验方差分析怎么用SPSS做?

需要资料分析+qq

正交实验法中的F比怎么算

F值是均方与自由度的比值,F比值是F值与其在相应显著性水平下的F临界值的比值,二者都可以作为反应显著性差异的引数

分类: 教育/科学 >> 科学技术

问题描述:

正交实验的目的是什么?如果最佳参数不在你的设计之内怎么确定实验的最佳参数?

解析:

正交实验设计

当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。因此就出现了分式析因设计(fractional factorial designs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。

正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。例如作一个三因素三水平的实验,按全面实验要求,须进行33=27种组合的实验,且尚未考虑每一组合的重复数。若按L9(3)3正交表按排实验,只需作9次,按L18(3)7正交表进行18次实验,显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。

1.正交表

正交表是一整套规则的设计表格,用 。L为正交表的代号,n为试验的次数,t为水平数,c为列数,也就是可能安排最多的因素个数。例如L9(34), (表11),它表示需作9次实验,最多可观察4个因素,每个因素均为3水平。一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表,如L8(4×24) (表12),此表的5列中,有1列为4水平,4列为2水平。根据正交表的数据结构看出,正交表是一个n行c列的表,其中第j列由数码1,2,… Sj 组成,这些数码均各出现N/S 次,例如表11中,第二列的数码个数为3,S=3 ,即由1、2、3组成,各数码均出现 次。

正交表具有以下两项性质:

(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,且每对出现数也均相等。

以上两点充分的体现了正交表的两大优越性,即“均匀分散性,整齐可比”。通俗的说,每个因素的每个水平与另一个因素各水平各碰一次,这就是正交性。

2. 交互作用表 每一张正交表后都附有相应的交互作用表,它是专门用来安排交互作用试验。表14就是L8(27)表的交互作用表。

安排交互作用的试验时,是将两个因素的交互作用当作一个新的因素,占用一列,为交互作用列,从表14中可查出L8(27)正交表中的任何两列的交互作用列。表中带( )的为主因素的列号,它与另一主因素的交互列为第一个列号从左向右,第二个列号顺次由下向上,二者相交的号为二者的交互作用列。例如将A因素排为第(1)列,B因素排为第(2)列,两数字相交为3,则第3列为A×B交互作用列。又如可以看到第4列与第6列的交互列是第2列,等等。

3.正交实验的表头设计 表头设计是正交设计的关键,它承担着将各因素及交互作用合理安排到正交表的各列中的重要任务,因此一个表头设计就是一个设计方案。

表头设计的主要步骤如下:

(1)确定列数 根据试验目的,选择处理因素与不可忽略的交互作用,明确其共有多少个数,如果对研究中的某些问题尚不太了解,列可多一些,但一般不宜过多。当每个试验号无重复,只有1个试验数据时,可设2个或多个空白列,作为计算误差项之用。

(2)确定各因素的水平数 根据研究目的,一般二水平(有、无)可作因素筛选用;也可适用于试验次数少、分批进行的研究。三水平可观察变化趋势,选择最佳搭配;多水平能以一次满足试验要求。

(3)选定正交表 根据确定的列数(c)与水平数(t)选择相应的正交表。例如观察5个因素8个一级交互作用,留两个空白列,且每个因素取2水平,则适宜选L16(215)表。由于同水平的正交表有多个,如L8(27)、L12(211)、L16(215),一般只要表中列数比考虑需要观察的个数稍多一点即可,这样省工省时。

(4)表头安排 应优先考虑交互作用不可忽略的处理因素,按照不可混杂的原则,将它们及交互作用首先在表头排妥,而后再将剩余各因素任意安排在各列上。例如某项目考察4个因素A、B、C、D及A×B交互作用,各因素均为2水平,现选取L8(27)表,由于AB两因素需要观察其交互作用,故将二者优先安排在第1、2列,根据交互作用表查得A×B应排在第3列,于是C排在第4列,由于A×C交互在第5列,B×C交互作用在第6列,虽然未考查A×C与B×C,为避免混杂之嫌,D就排在第7列。

(5)组织实施方案 根据选定正交表中各因素占有列的水平数列,构成实施方案表,按实验号依次进行,共作n次实验,每次实验按表中横行的各水平组合进行。例如L9(34)表,若安排四个因素,第一次实验A、B、C、D四因素均取1水平,第二次实验A因素1水平,B、C、D取2水平,……第九次实验A、B因素取3水平,C因素取2水平,D因素取1水平。实验结果数据记录在该行的末尾。因此整个设计过程我们可用一句话归纳为:“因素顺序上列、水平对号入座,实验横着作”。

4.二水平有交互作用的正交实验设计与方差分析

例8 某研究室研究影响某试剂回收率的三个因素,包括温度、反应时间、原料配比,每个因素都为二水平,各因素及其水平见表16。选用L8(27)正交表进行实验,实验结果见表17。

首先计算Ij 与IIj ,Ij为第j列第1水平各试验结果取值之和,IIj为第j列第2水平各试验结果取值之和。然后进行方差分析。过程为:

求:总离差平方和

各列离差平方和 SSj=

本例各列离均差平方和见表10最底部一行。即各空列SSj之和。即误差平方和

自由度v为各列水平数减1,交互作用项的自由度为相交因素自由度的乘积。

分析结果见表18。

从表18看出,在α=0.05水准上,只有C因素与A×B交互作用有统计学意义,其余各因素均无统计学意义,A因素影响最小,考虑到交互作用A×B的影响较大,且它们的二水平为优。在C2的情况下, 有B1A2和B1,A1两种组合状况下的回收率最高。考虑到B因素影响较A因素影响大些,而B中选B1为好,故选A2B1。这样最后决定最佳配方为A2B1C2,即80℃,反应时间2.5h,原料配比为1.2:1。

如果使用计算机进行统计分析,在数据是只需要输入试验因素和实验结果的内容,交互作用界的内容不用输入,然后按照表头定义要分析的模型进行方差分析。

正交试验设计(Orthogonal experimental design)是研究多因素多水平的又一种设计方法,它是根据正交性从全面试验中挑选出部分有代表性的点进行试验,这些有代表性的点具备了“均匀分散,齐整可比”的特点,正交试验设计是分式析因设计的主要方法。是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交试验选择的水平组合列成表格,称为正交表。

当析因设计要求的实验次数太多时,一个非常自然的想法就是从析因设计的水平组合中,选择一部分有代表性水平组合进行试验。因此就出现了分式析因设计(fractional factorial designs),但是对于试验设计知识较少的实际工作者来说,选择适当的分式析因设计还是比较困难的。 例如作一个三因素三水平的实验,按全面实验要求,须进行3^3=27种组合的实验,且尚未考虑每一组合的重复数。若按L9(3^4)正交表安排实验,只需作9次,按L15(3^7)正交表进行15次实验,显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛套用。

基本介绍 中文名 :正交试验 外文名 :Orthogonal experimental design 类别 :一种实验设计方法 特点 :“均匀分散,齐整可比” 基本思想,正交表,方案设计,数据分析, 基本思想 正交试验设计法,就是使用已经造好了的表格--正交表--来安排试验并进行数据分析的一种方法。它简单易行,计算表格化,使用者能够迅速掌握。下边通过一个例子来说明正交试验设计法的基本思想。 [例1]为提高某化工产品的转化率,选择了三个有关因素进行条件试验,反应温度(A),反应时间(B),用碱量(C),并确定了它们的试验范围:

A:80-90℃

B:90-150分钟

C:5-7%

试验目的是搞清楚因子A、B、C对转化率有什么影响,哪些是主要的,哪些是次要的,从而确定最适生产条件,即温度、时间及用碱量各为多少才能使转化率高。试制定试验方案。

这里,对因子A,在试验范围内选了三个水平;因子B和C也都取三个水平:

A:A1=80℃,A2=85℃,A3=90℃

B:B1=90分,B2=120分,B3=150分

C:C1=5%,C2=6%,C3=7%

当然,在正交试验设计中,因子可以是定量的,也可以是定性的。而定量因子各水平间的距离可以相等,也可以不相等。

这个三因子三水平的条件试验,通常有两种试验进行方法: (Ⅰ)取三因子所有水平之间的组合,即A1B1C1,A1B1C2,A1B2C1, ……,A3B3C3,共有

3^3=27次

试验。用图表示就是图1 立方体的27个节点。这种试验法叫做全面试验法。

全面试验对各因子与指标间的关系剖析得比较清楚。但试验次数太多。特别是当因子数目多,每个因子的水平数目也多时。试验量大得惊人。如选六个因子,每个因子取五个水平时,如欲做全面试验,则需5^6=15625次试验,这实际上是不可能实现的。如果套用正交实验法,只做25次试验就行了。而且在某种意义上讲,这25次试验代表了15625次试验。图1 全面试验法取点。 (Ⅱ)简单对比法,即变化一个因素而固定其他因素,如首先固定B、C于B1、C1,使A变化之:

图1

JA1

B1C1 →A2

KA3 (好结果)

如得出结果A3最好,则固定A于A3,C还是C1,使B变化之:

JB1

A3C1 →B2 (好结果)

KB3

得出结果以B2为最好,则固定B于B2,A于A3,使C变化之:

JC1

A3B2→C2 (好结果)

KC3

图2

试验结果以C2最好。于是就认为最好的工艺条件是A3B2C2。

这种方法一般也有一定的效果,但缺点很多。首先这种方法的选点代表性很差,如按上述方法进行试验,试验点完全分布在一个角上,而在一个很大的范围内没有选点。因此这种试验方法不全面,所选的工艺条件A3B2C2不一定是27个组合中最好的。其次,用这种方法比较条件好坏时,是把单个的试验数据拿来,进行数值上的简单比较,而试验数据中必然要包含着误差成分,所以单个数据的简单比较不能剔除误差的干扰,必然造成结论的不稳定。

简单对比法的最大优点就是试验次数少,例如六因子五水平试验,在不重复时,只用5+(6-1)×(5-1)=5+5×4=25次试验就可以了。 考虑兼顾这两种试验方法的优点,从全面试验的点中选择具有典型性、代表性的点,使试验点在试验范围内分布得很均匀,能反映全面情况。但我们又希望试验点尽量地少,为此还要具体考虑一些问题。 如上例,对应于A有A1、A2、A3三个平面,对应于B、C也各有三个平面,共九个平面。则这九个平面上的试验点都应当一样多,即对每个因子的每个水平都要同等看待。具体来说,每个平面上都有三行、三列,要求在每行、每列上的点一样多。这样,作出如图2所示的设计,试验点用⊙表示。我们看到,在9个平面中每个平面上都恰好有三个点而每个平面的每行每列都有一个点,而且只有一个点,总共九个点。这样的试验方案,试验点的分布很均匀,试验次数也不多。

当因子数和水平数都不太大时,尚可通过作图的办法来选择分布很均匀的试验点。但是因子数和水平数多了,作图的方法就不行了。

试验工作者在长期的工作中总结出一套办法,创造出所谓的正交表。按照正交表来安排试验,既能使试验点分布得很均匀,又能减少试验次数,图2正交试验设计图例而且计算分析简单,能够清晰地阐明试验条件与指标之间的关系。用正交表来安排试验及分析试验结果,这种方法叫

正交试验设计法

。

正交表的性质

(1)每列中不同数字出现的次数是相等的,如L9( ),每列中不同的数字是1,2,3,它们各出现3次; (2)在任意两列中,将同一行的两个数字看成有序数对时,每种数对出现的次数是相等的,如L9( ),有序数对共有9个:(1,1),(1,2),(1,3),(2,1),(2,2),(2,3),(3,1),(3,2),(3,3),它们各出现一次。 由于正交表有这两条性质,用它来安排试验时,各因素的各种水平的搭配是均衡的。 正交表 为了叙述方便,用L代表正交表,常用的有L8( ),L9( ),L16( ),L8(4× ),L12( ),等等。此符号各数字的意义如下:

L8( )

7为此表列的数目(最多可安排的因子数)

2为因子的水平数

8为此表行的数目(试验次数)

L16(2× )

有7列是3水平的

有1列是2水平的

L16(2× )的数字告诉我们,用它来安排试验,做16个试验最多可以考察一个2水平因子和7个3水平因子。

在行数为mn型的正交表中(m,n是正整数),试验次数(行数)=Σ(每列水平数-1)+ 1

如L8( ),

8=7×(2-1)+l

利用上述关系式可以从所要考察的因子水平数来决定最低的试验次数,进而选择合适的正交表。比如要考察五个3水平因子及一个2水平因子,则起码的试验次数为5×(3-1)+1×(2-1)+1=12(次)

这就是说,要在行数不小于13,既有2水平列又有3水平列的正交表中选择,L18(2× )适合。

正交表具有两条性质:(1)每一列中各数字出现的次数都一样多。(2)任何两列所构成的各有序数对出现的次数都一样多。所以称之谓正交表。

例如在L9( )中(见表1),各列中的1、2、3都各自出现3次;任何两列,例如第3、4列,所构成的有序数对从上向下共有九种,既没有重复也没有遗漏。其他任何两列所构成的有序数对也是这九种各出现一次。这反映了试验点分布的均匀性。

表1

方案设计 安排试验时,只要把所考察的每一个因子任意地对应于正交表的一列(一个因子对应一列,不能让两个因子对应同一列),然后把每列的数字"翻译"成所对应因子的水平。这样,每一行的各水平组合就构成了一个试验条件(不考虑没安排因子的列)。 对于[例1],因子A、B、C都是三水平的,试验次数要不少于

3×(3-1)+1=7(次)

可考虑选用L9( )。因子A、B、C可任意地对应于L9( )的某三列,例如A、B、C分别放在l、2、3列,然后试验按行进行,顺序不限,每一行中各因素的水平组合就是每一次的试验条件,从上到下就是这个正交试验的方案,见表2。这个试验方案的几何解释正好是图2。

表2

三个3水平的因子,做全面试验需要3*3*3=27次试验,现用L9( )来设计试验方案,只要做9次,工作量减少了2/3,而在一定意义上代表了27次试验.。

再看一个用L9( )安排四个3水平因子的例子。

[例2]某矿物气体还原试验中,要考虑还原时间(A)、还原温度(B)、气体流速(C)、还原气体比例(D)这四个因子对全铁含量X〔越高越好)、金属化率Y(越高越好)、二氧化钛含量Z(越低越好)这三项指标的影响。希望通过试验找出主要影响因素,确定最适工艺条件。

首先根据专业知识以确定各因子的水平:

时间:A1=3(小时),A2=4(小时),A3=5(小时)

温度:B1=1000(℃),B2=1100(℃),B3=1200(℃)

流速:Cl=600(毫升/分),C2=400(毫升/分),

C3=800(毫升/分)

CO:H2:D1=1:2,D2=2:1,D3=1:1

这是四因子3水平的多指标(X、Y、Z)问题,如果做全面试验需3^4=81次试验,而用L9( )来做只要9次。具体安排如表3。同全面试验比较,工作量少了8/9。由于缩短了试验周期,可以提高试验精度,时间越长误差干扰越大。并且对于多指标问题,采用简单对比法,往往顾此失彼,最适工艺条件很难找;而套用正交表来设计试验时可对各指标通盘考虑,结论明确可靠。

表3

数据分析 正交表的另一个好处是简化了试验数据的计算分析。还是以[例1]为例来说明。按照表2的试验方案进行试验,测得9个转化率数据,见表4。 通过9次试验,我们可以得两类收获。 第一类收获是拿到手的结果。第9号试验的转化率为64,在所做过的试验中最好,可取用之。因为通过L9( )已经把试验条件均衡地打散到不同的部位,代表性是好的。假如没有漏掉另外的重要因素,选用的水平变化范围也合适的话,那么,这9次试验中最好的结果在全体可能的结果中也应该是相当好的了,所以不要轻易放过。

第二类收获是认识和展望。9次试验在全体可能的条件中(远不止3^3=27个组合,在试验范围内还可以取更多的水平组合)只是一小部分,所以还可能扩大。精益求精。寻求更好的条件。利用正交表的计算分折,分辨出主次因素,预测更好的水平组合,为进一步的试验提供有份量的依据。

其中I、Ⅱ、Ⅲ分别为各对应列(因子)上1、2、3水平效应的估计值,其计算式是:

Ⅰi(Ⅱi,Ⅲi)=第i列上对应水平1(2,3)的数据和

K1 为1水平数据的综合平均=Ⅰ/水平1的重复次数

Si为变动平方和=

[例1]的转化率试验数据与计算分析见表4。

先考虑温度对转比率的影响。但单个拿出不同温度的数据是不能比较的,因为造成数据差异的原因除温度外还有其他因素。但从整体上看,80℃时三种反应时间和三种用碱量全遇到了,85℃时、90℃时也是如此。这样,对于每种温度下的三个数据的综合数来说,反应时间与加碱量处于完全平等状态,这时温度就具有可比性。所以算得三个温度下三次试验的转化率之和:

80℃: ⅠA=x1+x2+x3=31+54+38=123;

85℃: ⅡA=x4+x5+x6=53+49+42=144;

90℃: ⅢA=x7+x8+x9=57+62+64=183。

分别填在A列下的Ⅰ、Ⅱ、Ⅲ三行。再分别除以3,表示80℃、85℃、90℃时综合平均意义下的转化率,填入下三行K1、K2、K3。R行称为极差,表明因子对结果的影响幅度。

同样地,为了比较反应时间;用碱量对转化率的影响,也先算出同一水平下的数据和IB、ⅡB、ⅢB,IC、ⅡC、ⅢC,再计算其平均值和极差。都填入表4中;

由此分别得出结论:温度越高转化率越好,以90℃为最好,但可以进一步探索温度更好的情况。反应时间以120分转化率最高。用碱量以6%转化率最高。所以最适水平是A3B2C2。

正交试验的方差分析

(一)假设检验 在数理统计中假设检验的思想方法是:提出一个假设,把它与数据进行对照,判断是否舍弃它。其判断步骤如下:

(1)设假设H 0 正确,得到一个理论结论,设此结论为R 0 ;

(2)再根据试验得出一个试验结论,与理论结论相对应,设为R 1 ;

(3)比较R 0 与R 1 :若R 0 与R 1 没有大的差异,则没有理由怀疑H 0 ,从而判定为:"不舍弃H。"(采用H。);若R 0 与R 1 有较大差异,则可以怀疑H 0 ,此时判定为:"舍弃H 0 "。

但是,R 1 /R 0 比值为多少才能舍弃H 0 呢?为确定这个量的界限,需要利用数理统计中F分布的理论。

若yl服从自由度为φ1的χ2分布,y2服从自由度为φ2的χ2分布,并且yl、y2相互独立,则(y1/φ1)/(y2/φ2)服从自由度为(φ1,φ2)的F分布。F分布是连续分布,分布模数是两个自由度(φ1,φ2)。称φ1为分子自由度,称φ2为分母自由度。在自由度为(φ1,φ2)的F分布中,某点右侧面积为p,也就是F比此值大的机率为p,把这个值写为 (p)。若检验的显著性水平(或危险率)给定为α时,则可以把 (α)作为临界值来检验假设。

这里,Se/σ2服从自由度为φe,的χ2分布;当H。成立,σ2=0时,SA/σ2也服从自由度为φA的χ2分布;又SA与Se相互成立,所以(SA/(φAσ2)/ Se/(φeσ2))=VA/Ve服从自由度为(φA,φe)的F分布。这就是假定H。正确时的理论结论R。。而试验结论Rl要与理论结论R。相比较。由给定的显著性水平,通常是α=0.05;分子自由度φ1=φA=a-1,分母自由度φ2=φe=a(n-1);查F分布表得出 (α)。所以H。:α1=α2=……=αa=0(σA2=0)的检验是:(显著性水平α)

FA=VA/Ve> (α) → 舍弃H 0;

FA=VA/Ve≤ (α) → 不舍弃H 0;

通常, (α)一般性地表示成Fα(φA,φB)。

假设因子A对试验结果的影响不显著,那么A的两个水平的效应该表现为相等或相近,即假设H 0 :α1=α2=0。如果因子A显著,则舍弃假设。

为了判断因子A是否显著,首先要计算比值显然,这个比值越大,因子A对指标的影响越显著;反之,因子A就不显著。在给定置信度α后,如α=0.05,查F分布表,自由度φA是因子A的,自由度φe是误差的,其临界值Fα(φA,φe),如果FA>Fα(φA,φe)就舍弃假设,可以认为因子A是显著的;如果FA≤Fα(φA,φe)就没有理由否定假设,而只能认为因子A是不显著的。因为按照F分布表的物理念义,F值小于Fα(φA,φe)的机率是95%,即有95%的机会出现小于Fα(φA,φe)的F值,既然出现了这种情况,就有了95%的把握,所以就没有理由否定假设,只能接受假设,认为因子A不显著。另一方面,F值大于Fα(φA,φe)的机率是5%,也就是只有5%的机会出现大于Fα(φA,φe)的F值,这是小机率事件,如果小机率事件居然发生了,则可认为情况异常,假设不可信,必须否定假设,因子A是显著的。对其他因子的显著性检验完全类似。

(二)方差分析表

由总平方和与各因素平方和即可求得误差平方和,亦称剩余平方和。是总平方和减各因素平方和所得。如正交表有一空列,则该列的平方和就是误差平方和。但在正交表饱和试验的情况下,即所有各列全部排满时,误差平方和一般用各因素平方和中几个最小的平方和之和来代替,同时,这几个因素不再作进一步的分析。

自由度:φT=试验次数一1

φA,B…=水平数一1

φA×B=φA×φB

φe=φT-φA-φB-……-φD

表5