正交试验方差分析通俗易懂.docx
《正交试验方差分析通俗易懂.docx》由会员分享,可在线阅读,更多相关《正交试验方差分析通俗易懂.docx(20页珍藏版)》请在冰豆网上搜索。
正交试验方差分析通俗易懂
第十一章正交设计试验资料的方差分析
在实际工作中,常常需要同时考察3个或3个以上的试验因素,若进行全面试验,则试验的规模将很大,往往因试验条件的限制而难于实施。
正交设计是安排多因素试验、寻求最优水平组合的一种高效率试验设计方法。
第一节、正交设计原理和方法
(一)正交设计的基本概念
正交设计是利用正交表来安排多因素试验、分析试验结果的一种设计方法。
它从
多因素试验的全部水平组合中挑选部分有代表性的水平组合进行试验,通过对这部分试验结
果的分析了解全面试验的情况,找出最优水平组合。
例如,研究氮、磷、钾肥施用量对某小麦品种产量的影响:
A因素是氮肥施用量,设Ai、A2、A33个水平;
B因素是磷肥施用量,设Bi、B2、B33个水平;
C因素是钾肥施用量,设Ci、C2、C33个水平。
这是一个3因素每个因素3水平的试验,各因素的水平之间全部可能的组合有27种。
如果进行全面试验,可以分析各因素的效应,交互作用,也可选出最优水平组合。
但全面试验包含的水平组合数较多,工作量大,由于受试验场地、经费等限制而难于
实施。
如果试验的主要目的是寻求最优水平组合,则可利用正交设计来安排试验。
正交设计的基本特点是:
用部分试验来代替全面试验,通过对部分试验结果的分析,了解全
面试验的情况。
正交试验是用部分试验来代替全面试验,它不可能像全面试验那样对各因素效应、交互
作用一一分析;当交互作用存在时,有可能出现交互作用的混杂。
如对于上述3因素每个因素3水平试验,若不考虑交互作用,可利用正交表L9(34)安排,试
验方案仅包含9个水平组合,就能反映试验方案包含27个水平组合的全面试验的情况,找
出最佳的生产条件。
一、正交设计的基本原理
表11-133试验的全面试验方案
5GG
AjB2C2
AjB2C3
4
b2
AjB7Cj
AjB.C.
AjB3Cf
/屈G
ajb^j3
A尹G
A尹2s
Jii圧4
/岂q
島
/尹心
/尹£
A^B2C2
A尹&3
A3B3Cj
a3b3c2
A3B3C3
图14」3因素侮个因素3水平试验点的均衡分布图
正交设计就是从全面试验点(水平组合)中挑选出有代表性的部分试验点(水平组合)来进行试验。
图1中标有9'个试验点,就是利用正交表L%34)从27个试验点中挑选出来的9个试验点。
即:
("A1B1C1
(2)A1B2C2
(3)A1B3C3
(4)A2B1C2
(5)A2B2C3
(6)A2B3C1
(7)A3B1C3
(8)A3B2C1
(9)A3B3C2
上述选择
,保证了A
因素的每个水平与B因素
、C因素的各个水平在试验中各搭
配一次。
从图1
中可以看到,
9个试验点分布是均衡的
,在立方体的每个平面上有且仅有3
个试验点;每两个平面的交线上有且仅有1个试验点。
9个试验点均衡地分布于整个立方体内,有很强的代表性,能够比较全面地反映全面
试验的基本情况。
二、正交表及其特性
(一)正交表
表11-2是L8(27)正交表,其中“L代表正交表;L右下角的数字“8表示有8行,
用这张正交表安排试验包含8个处理(水平组合);括号内的底数“2”表示因素的水平数,
括号内2的指数“7表示有7列,用这张正交表最多可以安排7个2水平因素。
表11-2L8(27)正交表
试验号
列
号
1
2
3
4
5
6
7
1
1
i
1
1
1
1
1
2
1
1
1
2
2
2
2
3
1
2
2
1
1
2
2
4
1
2
2
2
2
1
1
5
2
1
2
1
2
1
2
6
2
j
2
2
1
2
i
3
2
2
1
1
2
2
1
8
2
2
1
2
1
1
2
2水平正交表还有L4(23)、L16(215)等;
3水平正交表有L9(34)、L27(313)、…、等。
(2)正交表的特性
1、任一列中,不同数字出现的次数相同
例如L8(27)中不同数字只有1和2,它们各出现4次;L9(34)中不同数字有1、2和3,它们各出现3次。
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次。
即每个因素的一个水平与另一因素的各个水平互碰次数相等,表明任意两列各个数字之间的搭配是均匀的。
用正交表安排的试验,具有均衡分散和整齐可比的特点。
均衡分散,是指用正交表挑选出来的各因素水平组合在全部水平组合中的分布是均
衡的。
由图11-1可以看出,在立方体中,任一平面内都包含3个试验点,任两平面
的交线上都包含1个试验点。
整齐可比是指每一个因素的各水平间具有可比性。
因为正交表中每一因素的任一水平下都均衡地包含着另外因素的各个水平,当比较某因
素不同水平时,其它因素的效应都彼此抵消。
如在A、B、C3个因素中,A因素的3个水
平A1、A2、A3条件下各有B、C的3个不同水平,即:
B1C]
B1C2
B1C3
Ai
B2C2
勉
B2C3
A3
B2CI
B3C3
B3C1
B3C2
在这9个水平组合中,A因素各水平下包括了B、C因素的3个水平,虽然搭配方式不同,但B、C皆处于同等地位,当比较A因素不同水平时,B因素不同水平的效应相互抵消,C因素不同水平的效应也相互抵消。
所以A因素3个水平间具有可比性。
同样,B、C因素
3个水平间亦具有可比性。
(3)正交表的类别
1、相同水平正交表各列中出现的最大数字相同的正交表称为相同水平正交表。
L4(23)、L8(27)、L12(211)等各列中最大数字为2,称为两水平正交表;L9(34)、L27(313)等各列中最大数字为3,称为3水平正交表。
2、混合水平正交表各列中出现的最大数字不完全相同的正交表称为混合水平正交
表。
L8(41X24)表中有一列最大数字为4,有4列最大数字为2。
也就是说该表可以安
排1个4水平因素和4个2水平因素。
L16(44X23),L16(4X212)等都混合水平正交表。
三、正交设计方法
【例111】某水稻栽培试验选择了3个水稻优良品种(A):
二九矮、高二矮、窄叶
青,3种密度(B):
15、20、25(万苗/666.7m2);3种施氮量(C):
3、5、8(kg/666.7m2),试采用正交设计安排一个试验方案。
(一)确定试验因素及其水平,列出因素水平表
表11-3因素水平表
水平
因
素
品种
(A)
密度但)
施氮量
(C)
1
二九嫌AJ
3(E)
2
高二煤
20(B2>
5(C2)
3
窄叶青(AJ
25(引
8(CJ
(二)选用合适的正交表
根据因素、水平及需要考察的交互作用的多少来选择合适的正交表。
选用正交表的原则是:
既要能安排下试验的全部因素(包括需要考查的交互作用),又
要使部分水平组合数(处理数)尽可能地少。
一般情况下,试验因素的水平数应恰好等于正交表记号中括号内的底数;因素的个数(包
括需要考查交互作用)应不大于正交表记号中括号内的指数;各因素及交互作用的自由度之
和要小于所选正交表的总自由度,以便估计试验误差。
若各因素及交互作用的自由度之和等于所选正交表总自由度,则可采用有重复正交试验
来估计试验误差。
此例有3个3水平因素,若不考察交互作用,则各因素自由度之和为因素个数X(水平数-1)
=3X(3-1)=6,小于L9(34)总自由度9-1=8,故可以选用L9(34);
若要考察交互作用,则应选用L27(313),此时所安排的试验方案实际上是全面试验方案。
(三)表头设计
表头设计就是把挑选出的因素和要考察的交互作用分别排入正交表的表头适当的列上。
在不考察交互作用时,各因素可随机安排在各列上;若考察交互作用,就应按该正交表的交互作用列表安排各因素与交互作用。
此例不考察交互作用,可将品种(A)、密度(B)和施氮量(C)依次安排在L9(34)的第1、2、
3列上,第4列为空列,见表2-4。
表11-4表头设计
列号1234因素ABC空
L9(34)表头设计
因素数
列
号
1
2]
|3|
4
2
A
B
AXBj
3
A
B
AXG
C
axb2AxCjBxC3
4
A
BXCt
BXD,
B
AXQAXDiCxDa
C
AXB1axd2BXD?
Daxb2AXCaBXQa
L8(27)表头设计
列号
1
2
■
4
5
6
7
3
A
B
AXB
AXC
BXC
4
A
AXBCXD
L
AXC
BXD
BXCAXD
1>
4
A'
BCXD
rAXB
pedTJXD
AXC|
D
BXC
AXD
5
A
DXE
B
CXD
AXSCXE
CBXD
AXCHXE
D
AXE
BXC
EAXB
(四)列出试验方案
把正交表中安排因素的各列(不包含欲考察的交互作用列)中的每个数字依次换成该因
素的实际水平,就得到一个正交试验方案。
表11-5正交试验方案
因
素
试輪号
A
B
C
1
2
3
1
1(二九煙)
1(询
2
1(二九瘫)
2(20)
2(5)
3
1〔二九矮〕
3(25)
3⑻
4
2〔高—燼
1(15)
2(5>
5
2(高二矮)
2(20)
3⑻
6
2〔高二矮)
3(25)
1(3>
1
3(窄叶青)
1(15)
3(8>
8
肌窄叶青)
2(20^
1⑶
9
3(宕叶青)
3(?
5)
2⑸
第二节正交试验资料的方差分析
若各号试验处理都只有一个观测值,则称之为单个观测值正交试验;
若各号试验处理都有两个或两个以上观测值,则称之为有重复观测值正交试验。
一、单个观测值正交试验资料的方差分析
对【例11-1】用L9(34)安排试验方案后,各号试验只进行一次,试验结果列于表试对其进行方差分析。
表11-6正交试验结果计算表
皿号
因素
~产量
A
(1)
Bm
C
P)
1
1
i
1
34O.O(xt)
2
1
7
422*5(X2)
3
1
3
3
4
2
1
2
dr
3百0咻)
5
2
2
3
492怒)
6
2
3
1
4矽.0(阳
7
3
1
3
392.0(x7)
8
3
2
1
9
3
3