第一章试验设计与试验数据分析初步

1、第一章 试验设计与试验数据分析初步在科学实验与工农业生产中 ,经常要做实验 .如何安排实验， 使实验次数尽量少， 而又能达到好的实验效 果呢？这是经常会碰到的问题 ,解决这个问题有一门专门的学问 ,叫做“试验设计” 。 20 世纪 30 年代，由于 农业试验的需要， 在试验设计和统计分析方面做出了一系列先驱工作， 从此试验设计成为统计科 学的一个分支 .随后， F. Yates， R.C。 Bose， O。 Kempthome ， W. G. Cochran, D. R. Cox 和 G. E. P。 Box 对试验设计都作出了杰出的贡献 ,使该分支在理论上日趋完善 ,在应用上日趋广泛。 19

2、60 年代, 日本统计学家 田口玄一将试验设计中应用最广泛的正交设计表格化，在方法解说方面深入浅出为试验设计的更广泛普及 与应用做出了巨大的贡献。试验设计方法有很多种类 .国内方面 ,60 年代由华罗庚教授倡导与普及的“优选法” ，即国外的斐波那契 方法，与 70 年代我国的数理统计学者在工业部门中普及的“正交设计”法都是人们熟悉的试验设计法。 70 年代末期由方开泰教授和王元教授倡导和推广的均匀设计也是一种常用的试验设计方法。“优选法 " 是单变量的最优调试法， “正交设计”是基于拉丁方理论和群论的多因素试验设计方法， “均匀设计”则是基于数 论方法的另一种多因素试验设计方法。材料

3、科学与材料工业中经常会出现各种试验设计问题，例如配方试验或称混料试验 (Experiments with Mixtures ),就是在材料科学中经常遇到的问题之一。在各种介绍试验设计方法的书籍中，通常都会辟出专 门的章节来阐述配方试验设计的解决方案 .试验设计得好，会事半功倍 ,反之就会事倍功半了。好的试验设计方案可以大大减少试验次数，得到充 分的信息，简化数据处理过程，节省人力、物力和时间。正确合理的试验设计，可以使试验结果的可靠性 显著提高。试验设计还可以为迅速寻求参数的优化数值和选择最佳工艺方案指明方向 .1。 1. 试验设计的概念及发展历史1.1.1. 试验设计的基本概念? 指标：在试

4、验设计中，人们把判断试验效果好坏所采用的标准称为试验指标，或简称为指标 .? 因素 : 或称为因子， 有可能影响试验指标的条件， 称作因素。 通常情况下固定的因素在试验方案中并不 称为因素，只有变化的因素才称为因素。? 水平：或称为处理，能影响试验指标的因素 , 通常人为地给予控制、分组，在统计学上，统称其为因子 的水平。1.1.2. 试验设计的基本要求? 进行不同处理的实验单元间，要有相同的系统? 要有明确的试验目的、恰当的指标? 要挑选因子 ,适当确定水平。 使试验范围尽可能大一点， 试验范围太小的缺点是不容易获得比已有条件 有显著改善的结果。 每一个因素的水平个数最好适当多一些 ,水平的

5、间隔大小和生产控制精度是密切相 关的。同时，因素和水平的含意可以是广义的? 处理实验数据要能配上相应的数理统计方法，以达到预期的试验目的1.1.3. 因素的主效应和因素间的交互效应 各因素的水平所对应的目标值称为主效应，各个因素除了对目标值有独立的影响外,还可能共同对目标值产生作用， 即交互作用 .交互作用通常表现为因素的乘积对指标的影响,其系数为正称为正交互作用， 其系数为负称为负交互作用 .1.1.4. 试验设计方法的发展过程? 全面试验法： 将每一个因素的不同水平组合做同样数目的试验。 一般说 m 个因子 n 个水平的全面试验 需要做 nm 次试验。当因素的个数不多，每个因数的水平数也不

6、多时，用全面试验的方法，并且通过 数据分析可以获得较为丰富的结果， 结论也比较精确。 当因素较多， 水平数较大时 ,全面试验要求的试 验数目可能非常大，虽然最后能够早出最好的搭配方案，但费时费工，往往不可能实现，因此除了一 些比较简单的情况外，一般不进行全面试验 .简单对比法：又称孤立因素法，是将因子中只变化一个，其余的固定，然后逐步地得到好的搭配的方 法。'这种方法一般也能得到一定的效果，而且比全面试验的次数少，但也有缺点，就是对待各因子 和水平不是均等的。 并且先固定那些因子， 后变化那些因子 ,都会影响试验结果， 因此最后的结果是不 是最好的，还不能充分肯定。完全随机试验法：是一

7、种最基本的试验设计，即各因素的水平完全随机分配随机区试验法：在划分区间后，用随机数字表或抽签法来确定顺序的试验设计方案 拉丁方试验法： 将一定数的文字排成正方形， 每个文字在各行各列都出现一次而且只出现一次 方格称为拉丁方。第一行与第一列相同的顺序排列的拉丁方称为标准型拉丁方。对于因子数多于, 这样的3的实验， 实验数将随着因子数的增加而快速地增加， 采用拉丁方试验法， 可以大大减少实验数 .组成拉丁 方区的必要条件是 ,在拉丁方区内行数等于列数等于水平数。拉丁方是供不存在交互作用的因子的实验设计用的 。拉丁方设计的主要目的是研究单因素不同水平对实验结果的影响。正交拉丁方法：两个或两个阶数相同

8、的拉丁方之间呈正交关系而言的,正交拉丁方试验一般在5*5， 77， 8*8 中进行。正交拉丁方的个数不超过拉丁方字母的个数减 1，对于 6*6 的拉丁方，则不存在正 交拉丁方法，拉丁方或正交拉丁方试验，可用随机法安排试验。1.2. 正交试验设计与结果分析1。 2.1 。 正交试验设计正交试验设计方法，就是利用数理统计学与正交性原理，从大量的试验点挑选适量的具有代表性、典 型性的点，应用“正交表”合理安排试验的一种科学的试验设计方法 .统计学家将正交设计通过一系列表格来实现，这些表格叫做正交表，记为Lt (me), 其中 L 表示正交表，t表示总共做9次试验，m表示每个因素都有 m个水平，e表示

9、表中有e列，最多可以安排e个因素。常用的二水平正交表有L4(23)，L8(27),Li6 (215)，L32( 231);三水平正交表有L9( 3勺，L27(313)；四 水平正交表有 L16( 45)；五水平正交表有L25(56)等。还有一批混合水平的表在实际中也十分有用，如L8(4 24),L12(23 31)， L16(44 23),L16(43 26)， L16(42 29)， L16(41 212)， L18(81 28)， L18(21 37)等。正交设计是利用数学上的正交性确定的设计方法，本质上具有“均匀分散、整齐可比" 的特点。具体表现为：1，水平均匀性：每个因子和因

10、子的每个水平都是均匀分配的；2,搭配的均匀性每个因子的各个水平出现的次数都是相同的，任何两个因子的搭配也都以相同的次数出现。正交试验过程 , 就是在确定指标、因子和水平后，用正交表安排试验方案 . 它主要要求解决三个方面的 问题：1，分析因子与指标的关系，即当因子变化时，指标是怎样变化。找出这种变化的规律，可以利用它 能动地指导生产 ;2 ，分析因子影响指标的主次， 即分析哪个因子是影响指标的主要因素，哪些是次要因素。 找出主要影响因素常常是生产中关键问题之一； 3，寻找好的生产工艺，即找到每个因子各取什么水平，会 得到最好的指标。也就是选出最优方案，这是生产中最需要解决的问题.用正交表安排试

11、验的步骤： 1，根据因子数和水平数以及试验条件的限制，选择合适的正交表； 2, 将各 因子放到表头的各列中 ;3 ，将各因子的水平安排到相应的表格中； 4, 形成试验方案。正交试验设计缺点在于它只适合于水平数不多的试验。通常情况每个因子有q 个水平，用正交表安排试验，至少要作q2次试验，当q较多时,q2非常大，很多情况无法进行试验.正交试验设计只考虑全局平衡， 只得到可能的优化方向。1.2 。 2。 直观分析利用正交表进行试验设计并分析结果一般有两种方法，即直观分析与方差分析。直观分析就是通过计算将各个因子、水平对试验结果质量指标影响的大小，用图形表示出来，通过直 观分析，综合比较 , 以确定

12、最优化试验方案的方法。直观分析的目的是 :? 因子与指标的变化规律，从 k1， k2,k3 与因子的关系图来考察?因子影响指标的主次顺序，从极差R来考察，R越大影响越大?选顶最优方案，看各个因子中哪些水平的平均指标最高，来获得最优方案?以实验验证最优方案举例如下：例：研究温度、压力、配比及时间四个因子对某产品质量指标的影响：温度取430、450、470°C，压力取10、20、30kg，配比取3 %、5% 7%,时间取1、2、3小时，质量指标越大越好。1,根据因子和水平数选用可以安排四因子、三水平的正交表L9（34）2，将温度、压力、时间和配比随机安排在正交表的各个列上3，安排各因素的

13、水平4，得到试验方案经过实验，得到各个试验方案的质量指标如下表：试验号列 号A（温度）1B （压力）2C （配比）3D （时间）4试验方案质量指标11111A1B1C1D12221222A1B2C2D25231333A1B3C3D34342123A2B1C2D35852231A2B2C3D16162312A2B3C1D26173132A3B1C3D25583213A3B2C1D37093321A3B3C2D164K1117135153147总和K2180183174168486K3189168159171k139455149总平均k26061585654k363565357R2416785，对

14、各个因子的每个水平所对应的质量指标求和得到Ki， K2， K36，在每一列中计算的 Ki, K2, K3平均值也，k2, k37，在每一列中计算极差，即用每一列的ki， k2， k3中最大者减去最小者所得到的值，记为R8，用每一个因子的k值对该因子的水平作图，观察各个因子与指标的变化规律9, 由极差R决定因子影响指标的主次顺序：温度压力时间配比10, 选顶最优方案：由各个因子中平均指标最高的水平组合而成试验方案11, 实验验证优选方案上述分析可以看岀，用正交表直观分析岀来的最优方案并不包含在正交表的九次试验之中，者说明利 用正交表还没能从全面试验中挑选岀好的方案来，但是，进行这九次有代表性的试

15、验也是完全必要的.在实际问题中，用来衡量试验效果的指标往往不止一个，而是多个，这类试验叫做多指标试验，进行多指标试验设计的分析方法有综合评分法和综合平衡法。综合评分法，是将多指标化为一个评分指标来进行直 观分析的方法；综合平衡法是先分别将各个指标按单指标进行计算和分析，再将各指标的分析结果进行综 合平衡，以得到“最优”试验方案。综合评分法和综合平衡法所得到的结论，一般是一致的，在实际应用 中综合评分法用得较多，因为计算上比较方便。有些试验中有的因子可以取较多的水平，但是也有的因子受条件的限制，不能多选水平，这就遇到水 平数不同的情况。解决不同水平数问题的方法：可以直接用不同水平混合型的正交表；

16、也可以采用拟水平方法，即将水平一律补齐为相同水平再加以解决.对于因子间有相互作用，即存在交互作用的正交试验设计，采用专门的有交互作用的正交表 .交互作用放 在哪一列上，由二列间交互作用表决定.1。23方差分析直观分析方法，其优点是简单、直观、计算量较小.但是，直观分析不能给岀误差大小的估计，也就不能知道结果的精度。方差分析可以弥补直观分析的不足之处在一批实验数据中，数据的算术平均值代表了数据的平均水平 ，反映了数据的集中性；而数据的方差，反 映了数据的波动性，即数据的分散性，方差大小表明数据变化的显著程度，而数据变化的显著程度，又反映了 因素对指标影响的大小。例 T8钢淬火试验(四因素二水平)

17、因素A温度B时间A * B交互作用C冷却液D操作方法水平1800°C15 minA1*B1，A2*B2油D12820°C11 minA1*B2，A2*B 1水D2其中A和B有交互作用，测试淬火硬度，硬度越大越好 选用L8(27)进行表头设计误差误差ABA * BCeeD1234567正交表及试验结果如下:列号A(0C)1B(min)2A * B3C4e误差5e误差6D7硬度硬度-55试验号1111111150-521112222594312211225614122221158352121212550621221215837221122147882212112523l32-1

18、2-124-57ll8-777-102m13/41/2-331-5/4-7/4m2-27/47/47/41/401/2T111/89/8-19/819/83/8-5/8-9/8T211/89/819/819/83/85/89/8l294144144162549ll264494949104s121/881/8361/8361/89/825/881/8表中 Ti=m (T/8)(i=1 , 2)以A因素为例，令硬度值分别为Y1 , Y2 , Y3 , Y4 , Y5 , Y6 , Y7,Y8,则|a=Y1+Y2+Y3+Y4 ，IIa=Y5+Y6+Y7+Y8T= Y1+Y2+Y3+Y4+ Y5+Y6

19、+Y7+Y8Sa= 4(Ia/4) - (T/8) 2+( (I I a /4)(T/8) 2=4: (Ia2+IIa 2)/16 T ( Ia +1I a)/16+T2/32=4: (Ia2+|Ia 2)/16 T2/32 =(I a 2+| I a 2)/4T2/8即： S = (l2+|2) /4 -T2/8由此计算出：Sa=121/8 ,Sb=81/8, Sa * b=361/8, Sc=81/8 ,Sd=81/8,Se= 9/8 + 25/8 =34/8用F检验法检验各个因素对指标影响的显著性：(各因素方差/因素方差的自由度)：(误差方差/误差方差自由度)=F即：F =(S/f) /

20、 (S0fo)由此计算结果如下表：因素方差S自由度fS/fF值显著性A（温度）121/81121/8(121/8): (17/8) =7.1B （时间）81/8181/8(81/8): (17/8)=4.8C （冷却液）361/81361/8(361/8): (17/8)=21。2*D （操作方法）81/8181/8(81/8):(17/8)=4.8A*B (交互作用)361/81361/8(361/8) :(17/8)=21。2*e (误差)34/8217/8查表可以得到F检验值：Fo。25(1,2)=2.57, Fo。05(1,2)=18.51, F0.01 ( 1, 2)=98。5F值大

21、的因素对指标的影响大,F值小的因素对指标的影响小。故此，可以得到，C和A*B对硬度指标的影响较显著，A和B、D的影响依次减小.筛选最优方案,C是最显著因素，先选 C，C列的m2m1，故此选C2方案;A * B也是显著因素，A * B 的m2m1,选二水平得到 A、B的搭配A1B2或A2B1，A比B重要，先选 A，A的m1>m2,故选定A1B2; D 不显著，可以选 D2,因此最优方案是 A1B2 C2D2。1。3。均匀试验设计所有的试验设计方法本质上就是在试验的范围内给岀挑选代表点的方法。正交设计是根据正交性准则 来挑选代表点，使得这些点能反映试验范围内各因素和试验指标间的关系。正交设计

22、在挑选代表点时有两 个特点：均匀分散，整齐可比。均匀分散时试验点有代表性；整齐可比便于试验数据的分析。为了保证整 齐可比的特点，正交设计至少要求做q2次试验.若要减少试验的数目，只有去掉整齐可比的要求。均匀设计就是只考虑试验点在试验范围内均匀散布的一种试验设计方法。每一个均匀设计表都有一个代号Un ( qs)或Un* (qs)，其中“ U"表示均匀设计，“n”表示要做n次试验，q表示每个因素有q个水平，s表示该表有s列；U的右上角加“ *"和不加“ *"代表两种不同类型的均匀 设计表，通常加“ *”的均匀设计表有更好的均匀性，应优先选用。但是不加“*”的均匀设计表

23、表能安排更多的因素，所以当因素数s较大，且超过加“ *”的均匀设计表的使用范围时，可使用不加“*"的均匀设计表每一个均匀设计表都附有一个使用表，它指示我们如何从均匀设计表中选用适当的列，以及由这些列 所组成的试验方案的均匀度。使用表最后一列D表示刻划均匀度的偏差(discrepancy)，偏差值越小，表示均匀度越好均匀设计有其独特的布点方式，其特点表现在:1，每个因素的每个水平做一次且只做一次试验2，任何两个因素的试验点点在平面的格子上，每行每列有且仅有一个试验点 性质(1 )和反映了试验安排的均衡性，即对各因素，每个因素的每个水平一视同仁3，均匀设计表的任何两列组成的试验方案一般并

24、不等价，因此每个均匀设计表必须有一个使用表 4，当因素的水平数增加时，试验数按水平数的增加量在增加均匀设计表U6* ( 64)1234112362246533624441535531266541均匀设计表U7 ( 74)123411236224653362444153553126654177777均匀设计表U7* (74)123411357226263317544444557136626277531U6* ( 64)的使用表s列号D2130。187531230.2656412340。2990U7(74)s列号D2130.239831230.3721412340.4760U7* ( 74)s列

25、号D2130.158232340.2132例题：在阿魏酸的合成工艺考察中，为了提高产量，选取了原料配比(A )、吡啶量(B)和反应时间(C)三个因素，它们各取了 7个水平如下：原料配比(A) :1。0, 1。4，1.8, 2.2, 2。6， 3.0, 3。4吡啶量(B)(ml) :10，13, 16,19, 22, 25，28反应时间（C） （h） :0。5,1.0,1。5, 2.0, 2.5, 3。0, 3.5根据因素和水平，我们选取均匀设计表U7*（74）或U7 （74），由它们的使用表可以查到，当s=3时，两个表的偏差分别为 0.2132和0。3721，故应当选用 5”（74）来安排试验

26、，其试验方案列于下表。该方案是 将A、B、C分别放在U7*（74）表的后3列而获得的。No。配比（A）吡啶量（B）反应时间（C）收率（Y）11230.33022460.33633620。29444150.47655310。20966540。45177770.482No.配比（A）吡啶量（B）反应时间（C）13: 1.85: 227: 3。526: 3。02: 136: 3.031: 1。07: 285: 2。544:2。24: 194:2.057:3。41: 103: 1。562: 1。46:252 : 1。075: 2。63: 161: 0.5试验数相同时的偏差的比较OD& UDS=

27、2S=3S=4S=5L8 (27)0.43750。 57810.6836Us*(88)0。14450。 20000。 2709L9(34)0。30560。 42130.5177U9(95)0.19440.31020。 4066L 12 (211)0.43750。 57810.68360。 7627U12* ( 1210)0.11630.18380。 22330.2272L16 (25)0。43750.57810.68360。 7627U16*(1612)0。09080。 12620.17050。2070L16 (45)0.23440.33010。 41380.4871U16* ( 1612)0。09080。 12620.17050。2070L25 ( 56)