试验设计与数据处理第五章-正交试验设计的课件

1、第一节 正交表介绍 列号试验号1231111212232124221 正交表是已经制作好的规格化的表, 是正交试验设计的基本工具. (见附表5) 正交表L4 (23 )-最小的正交表常见的正交试验表:1. L4(23)有4个横行和3个纵列，由数码“1”和“2”组成. 它 有两个特点： (1)每纵列恰有两个“1”和两个“2”； (2)任意两个纵列，其横方向形成的四个数字对中，恰 好(1，1)、(1，2)、(2，1)和(2，2)各出现一次. 这就是 说对于任意两个纵列，数码“1”、“2”间的搭配是均衡 的. 2. L8(27)有8个横行和7个纵列，由数码“1”和“2”组成. 它 有两个特点： (1

2、)每纵列恰有四个“1”和四个“2”； (2)任意两个纵列，其横方向形成的八个数字对中，恰 好(1，1)、(1，2)、(2，1)和(2，2)各出现两次. 这就是 说对于任意两个纵列，数码“1”、“2”间的搭配是均衡的. 3. L9(34)有9个横行和4个纵列，由数码“1”、“2”和”3”组成. 它有两个特点： (1)每纵列“1”、“2”和“3”出现的次数相同，都是三次; (2)任意两个纵列，其横方向形成的九个数字对中， (1，1)、(1，2)、(1，3)、(2，1)、(2，2)、(2，3)、(3， 1)、(3，2)和(3，3)出现的次数相同，都是一次，即任 意两列的数码“1”、“2”和“3”间的

3、搭配是均衡的. L n( t q)数码数（水平数） 正交表的纵列数(最多安排的因素数)正交表的横行数（试验次数）正交表代号第二节 二水平正交试验和分批试验 分批试验的思想: 第一批撒大网第二批撒小网例5.2.1 2,4二硝基苯肼的工艺改革试验目的: 2,4二硝基苯肼是一种用于检验醛酮色层分析及肝功能试验的试剂产品. 若产品外观出现紫色就为不合格品, 桔黄色为合格品. 北京化工厂决定采用水合肼与氯代苯在溶剂乙醇中直接合成新工艺. 即在氯代苯与乙醇的混合溶液中, 加入水合肼进行搅拌, 反应中途再追加200ml乙醇. 经小试验证, 虽为合格品, 但产率只有45%.通过正交试验, 希望找到提高产率的工

4、艺条件.考核指标: 产率和颜色第一批撒大网在多分批试验中, 第一批试验可先分两个水平 ( 二水平试验 ).优点: 节省试验次数, 提高试验效果1. 挑因素、选水平，制定因素水平表2. 利用正交表，确定试验方案3. 试验结果分析第二批撒小网 对第一批试验确定的主次因素，若认为不再需要进一步考虑的因素，则主要因素固定在好水平上(本例如C2、F2 )，次要因素固定在适当水平上(本例如A2、E2 ). 第二批试验中要考虑的因素包括第一批试验中确定的主要因素且认为有必要进一步详细考察的因素(本例如B1、D2 )以及第一批试验中未及考虑到而认为有必要在第二批试验中补充考虑的因素(本例如加料速度). 1.

5、挑因素、选水平，制定因素水平表2. 利用正交表，确定试验方案较优方案:用工业2,4二硝基氯代苯与粗品水合肼在溶剂乙醇中合成;水合肼用量为理论量的2.3倍;反应时间为2小时;反应温度掌握在6070之间;慢速加料;快速搅拌.效果: 平均产率超过80%, 均合格. 优质、高产、低消耗 3. 试验结果分析第三节 多水平正交试验和水平趋势图例5.3.1 晶体退火工艺改进试验目的: 检查癌细胞用到一种碘化钠晶体f40,要求应力越小越好,希望不超过2度.晶体的退火工艺是影响晶体质量的一个重要环节.北京综合仪器厂经过30多炉试验,采用的工艺如图所示:先将晶体放在炉内以50/小时的速度升温至600, 然后恒温6

6、小时,再用1.5A的电流使晶体降温至250, 最后断电, 使其自然降温. 在此工艺下, 其它指标都已合格, 但应力未能低于7度. 通过正交试验, 希望找到能降低应力的工艺条件.考核指标: 应力(度)2. 利用正交表，确定试验方案1. 挑因素、选水平，制定因素水平表 ( 略 )第一批撒大网(1)直观分析第5号试验最好, 第7号试验次之.(2)极差分析A2B2C1D3称为全体水平组合关于应力的可能好的水平组合.本例因素的主次顺序为:升温速度A(极差21.5) 恒温时间C(极差14)恒温速度B(极差12.5) 降温速度D(极差11.5) 3. 试验结果分析(3)水平趋势图-定量因素的用量与试验结果之

7、 和的关系图发现有用量选偏的因素(本例中为恒温温度B)是认识上的重大收获. 下批试验中把用量选准后，常能取得明显的进展.第二批撒小网 对第一批试验确定的主要因素A，认为不再需要进一步考虑，则固定在好水平A2上. 第二批试验中要考虑的因素包括第一批试验中确定的主要因素且认为有必要进一步详细考察的因素B、C、D.挑因素、选水平，制定因素水平表 B: 水平1- 450, 水平2- 400; (从水平趋势图可看出呈单调递增性)C: 水平1- 3小时, 水平2- 5小时;(考查能否省电)D: 水平1- 15/小时, 水平2- 25/小时; (再考查一次)2. 利用正交表，确定试验方案3. 试验结果- 这

8、批四个试验基本都消除了应力.利用SAS进行方差分析正交试验采用极差分析或方差分析都可以, 两者没有本质区别只是前者较方便而后者分析更精确, 以下见例5.3.1用SAS进行方差方析 (E531).一、利用菜单系统对每个因素进行单因素方差分析, 模型的平方和即该因素的平方和在正交试验的分析中与极差作用类似. 亦可用因素的原始数据值(xa xb xc)求水平趋势图(plot选项中选取均值图)因素A的水平趋势图:data E531;input A B C D y;cards;input A B C D xa xb xc y;cards;1 1 3 2 30 600 4 62 1 1 1 50 600

9、6 73 3 3 1 100 500 4 13;proc print;proc anova data=E531; class xa;model y=xa;means xa/snk ; /*选项snk要求作各水平组均值间的两两比较,显著性水平默认为=0.05,若需修正则在snk后再加选项alpha= */二、SAS编程proc anova data=E531; class xb;model y=xb;means xb/snk ;proc anova data=E531; class xc;model y=xc;means xc/snk ; proc anova data=E531; class

10、d;model y=d;means d/snk ;run;proc anova data=E531; class xa xb xc d;model y= xa xb xc d ; /*当各因素 单因素方差分析 同时进行时由于没有误差 自由度所以无法进行显著性分析*/means xa xb xc d/snk ; /*选项snk要求作 各水平组均值间的两两比较*/run;说明： 由四个因素xa, xb, xc, D的单因素方差分析表（输出结果略）可得它们的p值分别为0.1073 、0.6088、 0.5425、0.6398，同样也可得因素的主次顺序为:升温速度A (p值0.1073 ) 恒温时间C

11、 (p值0.5425)恒温速度B (p值0.6088) 降温速度D (p值0.6398). 当然比较各因素的SS模也能得此结果.Level of -y- xa N Mean Std Dev 30 3 5.0000000 3.60555128 50 3 4.3333333 3.78593890 100 3 11.6666667 4.16333200Level of -y- xb N Mean Std Dev 450 3 5.00000000 4.35889894 500 3 6.66666667 6.02771377 600 3 9.33333333 4.93288286 Level of -y

12、- xc N Mean Std Dev 2 3 9.66666667 4.72581563 4 3 6.33333333 6.50640710 6 3 5.00000000 3.46410162Level of -y- D N Mean Std Dev 1 3 9.33333333 3.21455025 2 3 6.33333333 0.57735027 3 3 5.33333333 8.38649708四个因素各水平组均值输出结果：第四节 正交试验的优良性单因素轮换法: 将其它因素固定在某一水平, 每次只改变欲研究因素的水平, 依次类推以研究A、B、C三因素三水平试验为例, 在研究因素A的效

13、应时, 先将因素B、C固定在某一水平如B1、C1; 如发现A3较好, 固定B2, 试验安排在A3B2C1、A3B2C2、A3B2C3; 如发现C1较好, 试验安排在A3C1B1、A3C1B2、A3C1B3; V如发现B2较好, 则确定A3B2C1为最佳分析条件. (请看下图中的七个黑点即为上述单因素轮换法安排试验时试验点的分布位置)缺点：实际考察的因素水平仅局限于欲考察的因素水平范围的局部区域, 不能全面反映因素影响的情况, 且不能考察因素之间的交互效应.全面试验法: 工作量大, 如三因素三水平全面试验要进行33=27次. (请看下图中的27个交叉点即为全面试验法安排试验时试验点的分布位置)正

14、交试验法: 利用正交性的试验设计方法三因素三水平试验可选用L9(34)正交表, 试验只需进行9次. (请看上图中的9个红点即为正交试验法安排试验时试验点的分布位置)正交性包括：1. 均衡分散性: 使每个试验点有强烈的代表性正交试验的试验次数(红点总数) ：全面试验试验次数(交叉点数) =每个平面上红点数 ：每个平面上交叉点数=每条棱上红点数 ：每条棱上交叉点数即 9/27=3/9=1/3 (三因素三水平场合)2. 整齐可比性: 使正交表的同一列可比较产生”好水平”. 上图说明A1的三次试验(红点)中因素B的三个水平和因素C的三个水平各出现一次，A2和A3有同样情况，就是说，对因素A的三个不同水平，受因素B和因素C的水平影响处于对等状态. 局部试验的最优点的平均名次=(N+1)/(n+1) ( N是全面试验次数, n是局部试验次数)如对于正交表L8(27)则有(128+1)/(8+1)14, 说明这8次试验的好结果在全面试验中的排位大致在前14名左右.第六节 正交表的灵活安排定量因素的试验用量或定性因素的状态叫做因素的水平，在确定参加正交试验的因素后，各种类型的正交表的格式是固定的，有时候各因素的水平必须灵活处理以适应