第8章正交试验设计的方差分析例题.docx
《第8章正交试验设计的方差分析例题.docx》由会员分享，可在线阅读，更多相关《第8章正交试验设计的方差分析例题.docx（23页珍藏版）》请在冰豆网上搜索。
第8章正交试验设计的方差分析例题
832考虑交互作用的三水平正交试验的方差分析（因学时有
限和正交表太大L27（313），不讲解！
只讲解二水平情况，因为二水平会，三水平自然也会！
）
例8-4运动发酵单细胞菌是一种酒精生产菌。
为了确定其发酵培养
基的最佳配方，进行了四因素三水平正交试验，试验指标为酒精浓度
（g/ml）。
表8-12给出了因素水平表，要求考察交互作用AxB、AX
C和AXD。
查附表7可得，本试验应选用L27（313）正交表，表头设计应按照“L27（313）二列间的交互作用表”进行。
本例只考虑一级交互作用（p=1）,所以每个三水平交互作用应占（m-1）P=（3-1）"=2列，即A
13
XB、AXc,和AXD在L27（3）正交表中各占二列。
表8-12因素水平表
\\因素
水平'、、
葡萄糖浓度％
A
酵母膏浓度％
B
培养温度（C）
C
培养基pH
D
1
5
0
25
5.0
2
15
0.5
30
6.0
3
25
1.0
35
7.0
表头设计时应避免混杂，试验方案及试验结果见表8-13。
由交互作用表可知，将因素A、B安排在第1、2列之后，第3、4
列为AXB交互作用列；再将C安排在第5列后，AXC交互作用在第
6、7列；最后将D安排在第9列，则AXD交互作用类落在第&10
列（当然也可将D安排在第8列,则第9、10列为AXD交互作用列）。
表8-13试验方案及结果分析L27（313）
试
A
B
AXB
C
AXC
AXD
D
AXD
试验
验
号
1
2
3
4
5
6
7
8
9
10
11
12
13
结果
yi
1
1
1
1
1
1
1
1
1
1
1
1
1
1
0.20
2
1
1
1
1
2
2
2
2
2
2
2
2
2
0.50
3
1
1
1
1
3
3
3
3
3
3
3
3
3
0.50
4
1
2
2
2
1
1
1
2
2
2
3
3
3
1.50
5
1
2
2
2
2
2
2
3
3
3
1
1
1
1.10
6
1
2
2
2
3
3
3
1
1
1
2
2
2
1.20
7
1
3
3
3
1
1
1
3
3
3
2
2
2
1.60
8
1
3
3
3
2
2
2
1
1
1
3
3
3
1.60
9
1
3
3
3
3
3
3
2
2
2
1
1
1
1.20
10
2
1
2
3
1
2
3
1
2
3
1
2
3
0.40
11
2
1
2
3
2
3
1
2
3
1
2
3
1
0.50
12
2
1
2
3
3
1
2
3
1
2
3
1
2
0.20
13
2
2
3
1
1
2
3
2
3
1
3
1
2
6.30
14
2
2
3
1
2
3
1
3
1
2
1
2
3
2.70
15
2
2
3
1
3
1
2
1
2
3
2
3
1
4.20
16
2
3
1
2
1
2
3
3
1
2
2
3
1
5.90
17
2
3
1
2
2
3
1
1
2
3
3
1
2
7.70
18
2
3
1
2
3
1
2
2
3
1
1
2
3
6.15
19
3
1
3
2
1
3
2
1
3
2
1
3
2
0.40
20
3
1
3
2
2
1
3
2
1
3
2
1
3
0.30
21
3
1
3
2
3
2
1
3
2
1
3
2
1
0.30
22
3
2
1
3
1
3
2
2
1
3
3
2
1
1.75
23
3
2
1
3
2
1
3
3
2
1
1
3
2
4.75
24
3
2
1
3
3
2
1
1
3
2
2
1
3
5.30
25
3
3
2
1
1
3
2
3
2
1
2
1
3
2.90
26
3
3
2
1
2
1
3
1
3
2
3
2
1
7.30
27
3
3
2
1
3
2
1
2
1
3
1
3
2
2.80
K1j
9.40
3.30
32.75
26.40
19.95
26.20
22.60
28.30
16.65
22.90
19.70
24.20
22.45
68.25
Kj
33.05
27.80
17.90
24.55
26.45
23.20
18.80
20.00
23.45
25.00
22.40
21.90
24.45
Kj
25.80
37.15
17.60
17.30
21.85
18.85
26.85
19.95
28.15
20.35
26.15
22.15
21.35
灯
88.36
10.89
1072.6
696.96
398.00
685.44
510.76
800.89
277.22
524.41
388.09
585.64
504.00
K；j
1092.3
772.84
320.41
602.70
699.60
538.24
353.44
400.00
549.90
625.00
501.76
479.61
597.80
K3j
665.64
1380.1
309.76
299.29
477.42
355.32
720.92
398.00
792.42
414.12
683.82
490.62
455.82
S
32.62
67.90
21.81
2.48
6.64
6.34
7.43
6.34
3.23
一、计算（计算过程省略）
1.计算各列各水平的K值（KijgKj）和KWKjKij）
各列各水平对应的试验数据之和Kij,K2j,K3j,及其平方和K：
K2j,
K2j,列于表8-13中,例如
9
Ka八Xi=0.20+0.50*2+1.50+1.10+1.20*2+1.60*2=9.40二Kn，
i4
K21二88.36
9
Ka='Xi=0.40+0.50+,,+6.15=33.05=©,©=1092.30
i生
9
Ka=Xi=0.40+0.30+,,+2.80=25.80=K31，K21=665.64
i生
表示AXB的有两列，即第3,4列，计算后可知
K13=32.75,K23=17.90;K33=17.60
K14=26.40;K24=24.55,K34=17.30
2.计算各列的偏差平方和（Sj）及其自由度（fj）
由式（8-4）,可知:
-CT
1m
S=—'Kijrid
r=n/m=27/3=9;
CT=Tn=1/27X68.252=172.52
所以Sj=lJKij2—172.52二（Kij2+Kj2+Kj2）-172.539yj9
S^=^=l（Kii2+K2i2+K5i2）-172.52
9
i
=；（88.36+i092.30+665.64）-i72.52
=32.62
S=S=,,=67.90,S
3一,,
=i6.67,
=5.i4
所以SxB=S+S=2i.8i
所以Sxc=S+S=6.64
所以Saxd=S+So=6.34
Sd=S二,，=7.43
所以S=Si+S2+S3=3.23
A=fB=fC=fD=2
axb=f3+f4=2+2=4,faxc=f6+f7=2+2=4
fAXD=f8+f10=2+2=4,fe=f11+f12+f13=2X3=6验算：
①S的验算
n
QT='x：
=0.22+0.52+,+2.82=320.98
i=1
St=Qt—CT=320.98-172.52=148.46
k
另外S=\Sj=32.62+67.90+,+0.55=148.45
j=i
②fT的验算
ft=n-仁27-1=26
另夕卜ft=13fi=13x（m-1）=13x（3-1）=26
•••计算过程无误.
3.计算方差
Vj=Sfj
VA二S/fa=32.62/2=16.31
同理可得
VB=33.95,Vc=1.24,Vd=3.72
Vaxb=5.45,Vaxc=1.66,Vaxd=1.59
Ve=0.538
t—均大于2,且fe=6>1,•••无需校正Ve!
Ve
：
、显著性检验（计算过程省略）
1.计算Fj
Fa=VA/Ve=16.31/0.538=30.32,
tF二V/Ve,
二同理可得,Fb=63.10,Faxb=10.13,Fc=2.30,Faxc=3.09
2.查Fa
Fc（f因,fe）=fc（2,6）,Fc（f交,fe）二Fc（4,6）
当a=0.05时,查得F0.05（2,6）=5.14,F0.05（4,6）=4.53;
当a=0.01时,查得F0.01（2,6）=10.92,F0.01（4,6）=9.15.
3.显著性检验
vFa=30.32和Fb=63.10均大于Fo.oi（2,6）=5.14,
•••因素A和B均高度显著（用**表示）；
又vFaxb=10.13>F0.01（4,6）=9.15,
•••交互作用AXB也高度显著（用**表示）;
又vFd=6.91,v介于Fo.o5（2,6）和Fo.oi（2,6）之间，
•因素D显著（用*表示）;
又vFc=2.30vFo.o5（2,6）=5.14,以及Faxc=3.09和Faxd=2.96均小
于Fo.05（4,6）=4.53,
•因素C及交互作用AXC和AXD匀不显著.
根据F值大小，可知各因素及交互作用对试验指标影响的主次顺序为：
B,A,AXB,D,AXC,AXD,C.
4.列方差分析表
表8-14方差分析表
方差来源
偏差平方和
自由度
方差
F值
Fa
显著性
A
32.62
2
16.51
30.32
Fo.o5（2,6）=5.14
**
B
67.90
2
33.95
63.10
Fo.01（2,6）=10.92
**
AXB
21.81
4
5.45
10.13
**
C
2.48
2
1.24
2.30
Fo.05（4,6）=4.53
AXC
6.64
4
1.66
3.09
Fo.01（4,6）=9.15
D
7.43
2
3.72
6.91
*
AXD
6.34
4
1.59
2.96
误差e
3.23
6
0.538
总和
148.45
26
三、最优工艺条件确定
因素AB及交互作用AXB都高度显著，但因在主次顺序中，AXB排在A、B之后，因此应优先考虑AB的优水平。
A和B的优水平确定了，其搭配也就随之确定，不必再通过A、B的二元表确定A与B搭配。
通过比较试验指标和K值大小，可知A和B的优水平，分别为A和B
因素D作用显著，但D与A的交互作用AXD不显著,故可不考虑交互作用，通过比较K值可知D的优水平为D3。
因素C作用不显著，可以降低成本和操作方便等方面来考虑选取最适水平。
对本例通过比较K值确定C的优水平为G。
因此，最优水平组合为ARC2D，即最优工艺条件为葡萄糖浓度15%酵母膏浓度
1.0%、培养温度30C和培养基pH值7.0。
最后，最好能在最佳条件ABC2C3下，再实施一次试验，测定试验指标值（即酒精浓度），在L27（313）正交表中，没有ARCD3这一组试验。
在正常情况下，ARGD3组合条件下的试验指标值，应大于表8—13中的最大Xi值，即第17号试验的xi7=7.70。
8.4混合型正交表的方差分析
混合型正交表的方差分析与等水平正交表的方差分析无本质的
区别，只是用公式时，要注意各列水平数的差别。
例8-5试对例7-2试验数据进行方差分析。
课本中为简化计算，对表7-5（p144）的试验数据X作了线性变换，实际上没有必要。
在不对Xi作变换的情况，请同学们自己再做一次方差分析，作为课外作业去完成。
且求Sj时，对二水平因素用通式和简化式分别计算！
总的偏差平方和
n_n
St=二.（Xj~"X）2=二.xi-CT
i1i=1
ft=n-1
因素的偏差平方和S分两种类型进行计算：
1、对于四水平因素
1m2
Sj二八Kj-CT
rij
（对二水平因素，也可用这一通式计算，建议全部用通式计算，以免产
生混乱！
）
m=4,r二n/m=8/4=2
fj=m-1
2、对于二水平因素，简化计算公式为：
Sj=^（Kij-K2j）2，n=8
n
fj=m-1，m=2
（方差分析和显著性检验，见书上p181）
讨论：
（1）方差分析法与极差分析法得出的各因素主次顺序相同，都是A、C、B;
（2）由方差分析可知，因素A显著，因素C不显著，而因素B对试验结果无影响（即将Sb并入Se中,及Se，Se+S，fe：
二fe+fB）；
（3）主要因素A的优水平A;不显著因素C,可根据具体情况确定其水平，为缩短加工时间，可选G水平，但从指标值看，还的选C2为好；对试验结果无影响的因素B,选B或B均可，从试验的指标可知，ABC2为最佳工艺条件，（即5号试验）。
因此，此时指标值最大。
（251cnV100g）
（极差分析结果：
ASG或ABO
方差分析结果：
ABG或ABG）
8.5重复试验和重复取样的方差分析（因时间有限，不讲解！
）
在实际工作中，用正交表安排试验时，为了提高试验及其统计分析的精确性和可靠性，往往采取重复试验和重复取样，在安排试验时，将同号试验重复做若干次，从而得到在同一条件下若干次试验的数据，叫做重复试验，若在一个试验中，同时抽取若干个样品进行测试，则叫做重复取样。
8.5.1重复试验的方差分析
在用正交表安排试验时，若表上各列已被因素及交互作用占满，没有空列，也无经验误差。
这时，为了估计试验误差，一般选用更大的正交表以外，还可以重复试验，由于正交本身的需要，有时虽然正交表的所有列并未被因素及交互作用占满，但也要做重复试验。
重复试验的方差分析与无重复试验的方差分析比较，有以下几点
不同：
（1）假设每号试验重复数为S,在计算Kij，Kj,,Kj时，是以各号试验下“S个试验指标数据之和”进行计算；
（2）重复试验时，总偏差平方和Sr及其自由度fT按下式计算：
nsr2
smx/
itt壬ns
ft=ns-1
式中:
n---试验条件数，即正交表的总试验号;
s---各号试验重复数
Xit---第i号试验第t次重复试验数据（i=1,2,n；=1,2,,,
s）；
T----所有试验数据之和（包括重复试验）；
ns
T=二工乂止
i=1t4
（3）重复试验时，各列偏差平方和（S）计算公式中的“水平重复数”改为“水平重复数乘以试验重复数”，修正项CT也有变化，S的自由度fj仍为水平数减1。
1mt2n
S二丄-Kj2-CT，CT=—，r=-
rsi吕nsm
fj=m-1
（4）重复试验时，总误差平方和包括空列误差Si和重复试验误差Se2,即
S=Si+$2
其总的自由度fe等于Sei的自由度fei与Se2的自由度fe2之和，即：
fe=fel+fe2
Se2及fe2的计算公式如下
nsns
Se2二二二xit2-1-（、xio2
iWt：
!
si=1t=1
fe2=n（S-i）
（5）重复试验时，用Ve=Se/fe检验各因素及其交互作用的显著性。
当
正交表的各列都已排满因素及交互作用而无空列时（即Se=O和
fe1=0）用Ve2=S^/fe2来检验因素及交互作用的显著性。
例8-6（p183）四因素四水平正交试验，每号试验重复三次，由
附表7可知，对四因素四水平试验，选Li6（45）正交表最合适，本例不考虑因素间的交互作用，因素水平如表8-17所示，而表8-18为试验方案与试验结果计算表。
一、计算（简略）
1.计算各列水平Kj值（Kj，K2j，K3j，Kj）和&
如Kii=6+12.5+17.5+19.2=55.2
K2=55.22=3047.04
心=19.2+19.5+18.9+19.2=76.8
K45=76.82=5898.24
2计算各列偏差平方和（S）及其自由度（fj）
m2
Sj=—乞Kj2一ct，CT=—，r=—
rsynsm
=1'Kj-303二丄「Kj-1912.69
43i16312心
如
Sa二$=右（K:
+K21+K312+K412）-1912.69
1
二丄X（3047.04+6528.64+7656.25+6320.25）-1912.6=49.99
12
同理可得Sb=S=33.42Sc=S3=29.01
Sd=S=13.54Se1=S=9.65
nsns
Sq=X二xit2-丄v（'xit）2
idtdSi#td
163163
2I2
二工工Xit■（、Xit）
i4tSi4t4
22
=（2+2+
2222
+6.5+6.9）-1/3X（6+12.5+
+20.42）=2050.32-2048.31=2.01
所以Se二Se+Sa=9.65+2.01=11.66fj=m-1=4-1=3
fA=fB=fC=fD=3
fe1=f5=4-1=3
fe2=n（S-1）=16X（3-1）=32
fe=fe1+fe2=3+32=35
验算:
①Sr
nsr2
Sr=_Xi2丄
i=1t=1ns
3032
163
163
mXi2
i：
!
t=1
=2050.32-1912.69=137.62
5
另外ST=7SjSe2=SaSbScSdSe,Se2
j三
=49.99+33.42+29.01+13.54+9.65+2.01=137.62
②f
ft=ns-1=163-1=47
5
另外fTHfe^fAfBfcfDfe1fe2
j廿
3.计算方差VJ=
、显著性检验
Va
16.66
Ve
"0.33
Vb_
11.41
Ve
"0.33
Vc
9.67
Ve
0.33
Vd
4.51
Ve
0.33
j
Ve
F
=50.48
A
F
=33.76
B
-29.30
F
-13.67
D
Fc
三、确定最优条件
T四个因素的作用均高度显著，且由F值大小可知因素作用的主次顺序为A、B、C、D
•••通过比较K值，可知各因素的优水平为AB、G、D3,故最优水平组合为ABGD3，表8-18的试验方案中无该水平组合的试验，所以应在最优水平组合下，再安排实施一次试验，并且其试验指标值应大于表8-18中的最大指标值。
8.5.2重复取样的方差分析由于重复试验使试验次数成倍增加而增加试验费用，故在实际工作中，更常用的是采用重复取样方法来提高试验的可靠性，重复取样与重复试验在误差偏差平方和的计算上完全一样，但重复取样的误差，反映的是原材料和产品的不均匀性与试样的测量误差，即局部（试验）误差；而重复试验的误差，反映的是整个试验过程中的各种干扰引起的误差，即整体误差。
通常，局部误差比整体误差要小，原则上不能用来检验各因素水平间是否存在差异，否则，会得到几乎全部因素及交互作用都是显著的不正确结论。
但是，若符合下面两种情况，则可以把重复取样得到的局部误差Se2当作试验误差Se,进行统计检验。
（1）正交表中各列已排满，无空列提供一次误差（Sei），这时，为
了少做试验而用重复取样误差（也）作为试验误差（Se），检验各因
素交互作用的显著性，若检验结果有一半左右的因素及交互作用不显
著，就可以认为这种检验是合理的；
（2）若重复取样得到的局部（试验）误差（S2）与整体（试验）误
差（Si）相差不大，也就是说，要求两类误差的F值:
对于给定的信度a,有F不显著，这时，就可以将Se2和Se合并作为试验误差，即
Se=Sei+Sa;
Fe=fei+fe2
但是，若F>Fa（fei,fe2），则两类误差有显著差异，不能合并使用。
例8-7三因素三水平正交试验，不考虑交互作用，因此，选用L（34）正交表最合适。
因素水平表见表8-20,试验方案见表8-2i（seepi88）。
重复取样三次，即s=3.
解：
一、计算
i.计算各列水平的K值（Kj,K2j和Kj）和Kj2
如Ki=0.655+0.657+0.787=2.099,K：
=2.0992=4.406
Ki3=0.760+1.305+0.657=2.722,K23=2.7222=7.409
K和K的计算结果，列于表8-2i中.
2.计算各列偏差平方和（S）及自由度（fj）检验:
Sj」mkj二丄-込」kj-2.866jrsidjns33yj939yj
Sa=Si=1/9X（7.684+15.4
升级会员 