6
表2.实验方案及实验结果的直观分析
列号 实验号
1 A wH2SO4 (%) 1 B mCuSO4· 5H2O(g) 1 C mZn (g) 1 空白列 2 10min内H2的 产率 32.62
2
3 4
2
3 1
1
1 2
2
3 3
1
3 1
40.40
41.07 34.97
5
6 7 8 9
2
3 1 2 3
此，因素水平变化所引起的波动，即因素A的偏差平方和SA应为：
SA=∑(yi-y总)2= (34.74-39.08)2+(38.71-39.08)2+(43.78-39.08)2
=123.37 上述计算结果我们可以通过S总=SA+ Se式来检验SA和 Se 计算正确与否。
(4)自由度和平均偏差平方和的计算 为了消除个数不同对实验指标所产生的影响，
S总 ( yi y总 )2
i 1
n
由表3知:
y总=1/9(32.62+34.97+36.62+40.40+…+44.53)=39.08
则：S总=(32.62-39.08)2+(34.97-39.08)2+…+(44.5339.08)2=151.08 S总反映了实验数据总的波动情况，如果硫酸质量分

i 1
i
即除以(n-1)，就得到平均偏差平方和。
平均偏差平方和 nS 1
为什么不除以n而要除以(n-1)呢？这是因为n个 数(y1, y2, y3, ……yn)之间并非彼此毫无关系，它们满 足的关系是： 1 n y yi n i 1 即n个数之和的均值为一定值，因此，n个数中 只有(n-1)个可“自由”变动，所以，求平均偏差平 方和时除以(n-1)，数学上将这个(n-1)称为S的自由 度。
4.因素的显著性判断 设因素A的F比为FA：
当FA >F0. 01 (n1, n2 )时，说明该因素水平的改变
对实验结果有很显著的影响，记作**。
当FA >F0. 05 (n1, n2 )时，说明该因素水平的改变
对实验结果有显著的影响，记作*。 当FA >F0. 10 (n1, n2 )时，说明该因素水平的改变 对实验结果有一定的影响，记作O。
114.09 117.25 120.34 38.03 39.08 40.11 2.08
C mZn (g)
122.77 115.23 113.68 40.92 38.41 37.89 3.03
空白列
10min内H2的 产率
119.9 117.56 114.22 39.96 39.18 38.07 1.89
为了弥补直观分析方法的不足，可采用方差分析 方法对实验结果进行计算分析。所谓方差分析就是将 因素水平(或交互作用)的变化引起的实验结果间的差 异与误差的波动所引起的实验结果间的差异区分开来 的一种数学方法。 方差分析的中心要点是：把实验数据总的波动分 解成两部分，一部分反映因素水平变化引起的波动， 另一部分反映实验误差引起的波动。即把数据总的偏 差平方和(S总)分解为因素的偏差平方和(SA、SB、SC ……)与误差的偏差平方和(Se)，并计算它们的平均偏 差平方和(也称均方和，或均方)，然后进行检验，最 后得出方差分析表。
二.方差分析中的一些基本概念
1.偏差平方和
方差分析的关键是对偏差平方和的分解，因此，
充分理解这一概念是至关重要的。
所谓偏差平方和是指一组数据中，各个数(y1, y2, y3……yn)与它们的算术平均数y之差的平方和。用符号 S来表示。即：
1 1 y ( y1 y2 ...... yn ) yi n n i 1
若以S1表示A1水平下实验误差所引起的波动，其值应 为：S1=(32.62-34.74)2+(34.97-34.74)2+(36.62-34.74)2 =8.0870。同理可以求出A2 、A3水平下实验误差所引 起的波动，其值分别为S2=7.8389，S3=11.7875 则，A因素的各个水平下总的偏差平方和应为： Se= S1+ S2+ S3=8.0870+7.8389+11.7875=27.71 (2) S总的计算 总的偏差平方和S总是指全部实验数据中，每个数据(yi) 与总平均值(y总)之差的平方和，即：
应采用平均偏差平方和，其计算公式为：
因素A的平均偏差平方和=SA/fA 误差的平均偏差平方和=Se/fe 式中SA、Se分别代表因素A和误差的偏差平方和 fA=A因素的水平数-1，它代表SA的自由度
fe=f总-fA，它代表Se的自由度
f总=总的实验次数-1，它代表S总的自由度
在本例中f总=9-1=8，fA=3-1=2， fe=8-2=6
F比 = S
f因素
误差
f 误差
为了判断F比值的大小所表明的物理意义(即F比值多大 时，可以认为实验结果的差异主要是由因素水平的 改变所引起的；其值多小时，可以认为实验结果的 差异主要是由实验误差所引起的)，这就需要有一个 标准来衡量F比值，此标准就是根据统计数学原理编 制的F分布表，F分布表列出了各种自由度情况下F比 的临界值。
数水平的改变对实验指标不发生影响，而且实验中
也没有误差产生的话，那么全部实验数据理应都一
样，即S总应等于零，但情况并非如此。
(3) S wH2SO4 (SA)的计算 对于因素A来讲，当它取一水平时，3次实验(即1、4、7实验)结 果的均值(y)应为： y1=1/3(y1+y4+y7)=1/3(32.62+34.93+36.62)=34.74 y1代表了3次一水平实验对H2产率的影响。同理： y2=1/3(y2+y5+y8)=38.71 y3=1/3(y3+y6+y9)=43.78 y2、y3分别代表了3次二水平和三水平实验对H2产率的影响。因
则
1 n 1 n X xi yi C n i 1 n i 1 X y C
于是
S ( xi x)2 [( yi C ) ( y C )]2 ( yi y)2
i 1 i 1 i 1
n
n
n
3. F比与F分布表 (1) F比
F比是指因素水平的改变引起的平均偏差平方和与误 S因素 差的平均偏差平方和的比值。即： (2) F分布表及其查阅方法
当实验所测得的n个数(y1, y2, y3, ……yn)数值较 大时，为了简化计算，可将每一个原始数据yi(i=1, 2, 3……n)都减去同一个常数C，这并不影响偏差平方 和的计算结果，但计算的工作量却简化了许多。
上述推论可通过以下简单换算予以证明。 若令Xi=yi-C (i=1, 2, ……n)
但是，极差值仅仅反映了各因素影响实验指标的主次 关系，它不能告诉我们各个因素对实验指标影响的程 度。也就是说，它既不能指明这些因素中哪个是影响 实验指标的关键因素，也不能提供一个标准，用来考 察、判断各个因素的作用是否显著。
第二：就因素A而言(因素B、C也类同)，其中k1、k2、 k3值之间的差异是如何产生的？是由于A因素水平不 同引起的呢？还是由于实验误差所造成的呢？还是 两者综合作用的结果？从直观分析角度是无法说清 楚的。 正是由于直观分析存在着上述的缺点，所以需 要采用方差分析的方法来弥补上述的不足。 1.单因素实验的方差分析 为了便于讨论，我们仍以实验室制取H2的因素 之一------A因素(硫酸的质量分数)为例，来说明单个 因素的实验数据的方差分析方法。
2 ( y y ) i i 1 n
n
则
S
为了计算方便，上式可简化为一种更常见的形式：
S yi 2 yi y y yi 2 ny 2
2 2 i 1 n i 1 i 1 i 1 n n n n
若令：
G yi
i 1
CT
ቤተ መጻሕፍቲ ባይዱ
G2 n
n
则
S yi 2 CT
方差分析是把实验数据总的波动(即数据的总的偏差平方 和S总)分解成两部分：一部分反映因素水平变化引起的波动 (即因素的偏差平方和)，对本例而言仅为S wH2SO4；另一部分 反映实验误差引起的波动(即误差的偏差平方和Se)。即： (1) Se的计算
表3.实验结果分析 参与wH2SO4某一水平的实验编号 A1(20%) 1 4 7 A2 (25%) 2 5 8 平均值y A3 (30%) 3 6 9 10minH2产率 A1(20%) 32.62 34.97 36.62 34.74 A2 (25%) 40.40 36.53 39.19 38.71 A3 (30%) 41.07 45.75 44.53 43.78
三.正交试验设计的方差分析 现以实验室制取H2为例，来说明正交设计的方 差分析的基本方法。若该实验所考察的因素、水平 如表1和表2所示。
表1. 因素水平
因素 水平 一 二 A wH2SO4 (%) 20 25 B mCuSO4· 5H2O(g) 0.4 0.5 C mZn (g) 4 5
三
30
0.6
2
2 3 3 3
1
2 2 3 1
3
2 3 2 1
36.53
45.75 36.62 39.19 44.53
列号 实验号
K1 K2 K3 k1 k2 k3 R
A wH2SO4 (%)
104.21 116.12 131.35 34.78 38.70 43.78 9.05
B mCuSO4· 5H2O(g)
正交试验设计的方差分析
一.方差分析的意义
前面我们介绍了正交设计方案及其结果的直 观分析，该方法简单明了，通俗易懂，计算工作 量少，便于普及和推广。但直观分析方法不能把 实验中由于实验条件的改变而引起的数据波动同 实验误差引起的数据波动区分开来，也就是说， 不能区分因素各水平所对应的实验结果间的差异， 究竟是由于因素水平不同引起的，还是由于实验 误差引起的。