试验设计与数据处理第五章正交试验设计课件.ppt
试验设计与数据处理第试验设计与数据处理第五章正交试验设计五章正交试验设计第1页,此课件共34页哦第一节第一节 正交表介绍正交表介绍 列列号号试验号试验号1231111212232124221 正交表是已经制作好的规格化的表,是正交试验设计的基本工具.(见附表5)正交表L4(23)-最小的正交表第2页,此课件共34页哦第3页,此课件共34页哦常见的正交试验表常见的正交试验表:1.L1.L4 4(23 3)有有4 4个横行和个横行和3 3个纵列,由数码个纵列,由数码“1 1”和和“2”组成.它它 有两个特点:有两个特点:(1)(1)每纵列恰有两个“1 1”和两个和两个“2 2”;(2)任意两个纵列,其横方向形成的四个数字对中,恰任意两个纵列,其横方向形成的四个数字对中,恰 好(1(1,1)、(1,2)2)、(2,1)和(2(2,2)各出现一次.这就是这就是 说对于任意两个纵列,数码说对于任意两个纵列,数码“1”、“2”间的搭配是均衡间的搭配是均衡 的.2.L8(27)有有8 8个横行和个横行和7个纵列,由数码个纵列,由数码“1”和“2”组成组成.它 有两个特点:有两个特点:(1)(1)每纵列恰有四个“1”和四个和四个“2 2”;(2)任意两个纵列,其横方向形成的八个数字对中,恰任意两个纵列,其横方向形成的八个数字对中,恰 好(1,1)1)、(1,2)2)、(2,1)和(2(2,2)各出现两次各出现两次.这就是这就是 说对于任意两个纵列,数码说对于任意两个纵列,数码“1 1”、“2 2”间的搭配是均衡的.第4页,此课件共34页哦3.L3.L9 9(34)有9个横行和4个纵列,由数码“1 1”、“2”和”3”组成.它有两个特点:它有两个特点:(1)(1)每纵列“1”、“2”和和“3 3”出现的次数相同,都是三次;(2)任意两个纵列,其横方向形成的九个数字对中,任意两个纵列,其横方向形成的九个数字对中,(1(1,1)、(1,2)、(1(1,3)、(2,1)、(2,2)2)、(2(2,3)3)、(3(3,1)1)、(3,2)2)和和(3(3,3)3)出现的次数相同,都是一次,即任 意两列的数码“1 1”、“2”和和“3”间的搭配是均衡的.L n(t q)数码数(水平数)正交表的纵列数(最多安排的因素数)正交表的横行数(试验次数)正交表代号第5页,此课件共34页哦第二节第二节 二水平正交试验和分批试验二水平正交试验和分批试验 分批试验分批试验的思想:第一批撒大网第一批撒大网第二批撒小网第二批撒小网第6页,此课件共34页哦例例5.2.1 2,4二硝基苯肼的工艺改革试验目的试验目的:2,4二硝基苯肼是一种用于检验醛酮色层分析及肝功能试验的试剂产品.若产品外观出现紫色就为不合格品,桔黄色为合格品.北京化工厂决定采用水合肼与氯代苯在溶剂乙醇中直接合成新工艺.即在氯代苯与乙醇的混合溶液中,加入水合肼进行搅拌,反应中途再追加200ml乙醇.经小试验证,虽为合格品,但产率只有45%.通过正交试验,希望找到提高产率的工艺条件.考核指标考核指标:产率和颜色第7页,此课件共34页哦第一批撒大网第一批撒大网在多分批试验中,第一批试验可先分两个水平(二水平试二水平试验验).优点:节省试验次数,提高试验效果1.挑因素、选水平,制定因素水平表挑因素、选水平,制定因素水平表第8页,此课件共34页哦2.利用正交表,确定试验方案利用正交表,确定试验方案第9页,此课件共34页哦3.试验结果分析试验结果分析第10页,此课件共34页哦第二批撒小网第二批撒小网 对第一批试验确定的主次因素,若认为不再需要进一步考虑的因素,则主要因素固定在好水平上(本例如C2、F2),次要因素固定在适当水平上(本例如A2、E2).第二批试验中要考虑的因素包括第一批试验中确定的主要因素且认为有必要进一步详细考察的因素(本例如B1、D2)以及第一批试验中未及考虑到而认为有必要在第二批试验中补充考虑的因素(本例如加料速度).第11页,此课件共34页哦1.挑因素、选水平,制定因素水平表挑因素、选水平,制定因素水平表2.利用正交表,确定试验方案利用正交表,确定试验方案第12页,此课件共34页哦较优方案较优方案:用工业2,4二硝基氯代苯与粗品水合肼在溶剂乙醇中合成;水合肼用量为理论量的2.3倍;反应时间为2小时;反应温度掌握在6070之间;慢速加料;快速搅拌.效果效果:平均产率超过80%,均合格.优质、高产、低消耗 3.试验结果分析试验结果分析第13页,此课件共34页哦第三节第三节 多水平正交试验和水平趋势图多水平正交试验和水平趋势图例例例例5.3.1 晶体退火工艺改进试验目的试验目的试验目的试验目的:检查癌细胞用到一种碘化钠晶体f40,f40,要求应力越小越好,希望不超过2度度.晶体的退火工艺是影响晶体质量的一个重要环节.北京综合仪器厂经过北京综合仪器厂经过3030多炉试验多炉试验,采用的工艺如图所示:先将晶体放在炉内以5050/小时的速度升温至小时的速度升温至600,然后恒温6小时小时,再用用1.5A的电流使晶体降温至250250,最后断电,使其自使其自然降温.在此工艺下,其它指标都已合格,但应力未能低于7 7度度.通过正交试验,希望找到能降低应力的工艺条件.考核指标考核指标:应力应力(度度)第14页,此课件共34页哦2.利用正交表,确定试验方案利用正交表,确定试验方案1.挑因素、选水平,制定因素水平表挑因素、选水平,制定因素水平表(略略)第一批撒大网第一批撒大网第15页,此课件共34页哦第16页,此课件共34页哦(1)直观分析直观分析第5号试验最好,第7号试验次之.(2)极差分析极差分析A2B2C1D3称为全体水平组合关于应力的可能可能好的水平组合好的水平组合.本例因素的主次顺序为:升温速度A(极差21.5)恒温时间C(极差14)恒温速度B(极差12.5)降温速度D(极差11.5)3.试验结果分析试验结果分析第17页,此课件共34页哦(3)水平趋势图水平趋势图-定量因素的用量与试验结果之定量因素的用量与试验结果之 和的关系图和的关系图发现有用量选偏的因素发现有用量选偏的因素(本例中为恒温温度本例中为恒温温度B)是认识上是认识上的重大收获的重大收获的重大收获的重大收获.下批试验中把用量选准后,常能取得明显下批试验中把用量选准后,常能取得明显下批试验中把用量选准后,常能取得明显下批试验中把用量选准后,常能取得明显的进展的进展的进展的进展.第18页,此课件共34页哦第二批撒小网第二批撒小网 对第一批试验确定的主要因素A,认为不再需要进一步考虑,则固定在好水平A2上.第二批试验中要考虑的因素包括第一批试验中确定的主要因素且认为有必要进一步详细考察的因素B、C、D.1.挑因素、选水平,制定因素水平表挑因素、选水平,制定因素水平表 B:水平1-450,水平2-400;(从水平趋势图可看出呈单调递增性)C:水平1-3小时,水平2-5小时;(考查能否省电)D:水平1-15/小时,水平2-25/小时;(再考查一次)第19页,此课件共34页哦2.利用正交表,确定试验方案利用正交表,确定试验方案3.试验结果试验结果-这批四个试验基本都消除了应力这批四个试验基本都消除了应力.第20页,此课件共34页哦利用利用SAS进行方差分析进行方差分析n n正交试验采用极差分析或方差分析都可以,两者没有本质区别只是前者较方便而后者分析更精确,以下见例5.3.1用SAS进行方差方析 (E531).一、利用菜单系统对每个因素进行单因素方差分析,模型的平方和即该因素的平方和在正交试验的分析中与极差作用类似.亦可用因素的原始数据值(xa xb xc)求水平趋势图(plot选项中选取均值图)第21页,此课件共34页哦因素因素A的水平趋势图的水平趋势图:第22页,此课件共34页哦data E531;data E531;data E531;data E531;input A B C D y;input A B C D y;input A B C D y;input A B C D y;cards;cards;cards;cards;input A B C D xa xb xc y;input A B C D xa xb xc y;cards;cards;1 1 3 2 30 600 4 62 1 1 1 50 600 6 73 3 3 1 100 500 4 13;proc print;proc print;proc print;proc print;proc anova data=E531;class xa;model y=xa;means xa/snk;/*选项snk要求作各水平组均值间的两两比较,显著性水平默认为=0.05,若需修正则在snk后再加选项alpha=*/二、SAS编程第23页,此课件共34页哦proc anova data=E531;proc anova data=E531;class xb;class xb;model y=xb;model y=xb;means xb/snk;means xb/snk;proc anova data=E531;proc anova data=E531;class xc;class xc;model y=xc;model y=xc;means xc/snk;means xc/snk;proc anova data=E531;proc anova data=E531;class d;class d;model y=d;model y=d;means d/snk;means d/snk;run;run;proc anova data=E531;proc anova data=E531;class xa xb xc d;class xa xb xc d;model y=xa xb xc d;model y=xa xb xc d;/*/*当各因素当各因素 单因素方差分析单因素方差分析 同时进行时由于没有误差同时进行时由于没有误差 自由度所以无法进行显著性分析自由度所以无法进行显著性分析*/means xa xb xc d/snk;means xa xb xc d/snk;/*/*选项选项snksnk要求作要求作 各水平组均值间的两两比较各水平组均值间的两两比较*/run;run;第24页,此课件共34页哦说明:说明:由四个因素xa,xb,xc,D的单因素方差分析表(输出结果略)可得它们的p值分别为0.1073、0.6088、0.5425、0.6398,同样也可得因素的主次顺序为:升温速度A(p值0.1073)恒温时间C(p值0.5425)恒温速度B(p值0.6088)降温速度D(p值0.6398).当然比较各因素的SS模也能得此结果.第25页,此课件共34页哦Level of -y-Level of -y-xa N Mean Std Dev xa N Mean Std Dev 30 3 5.0000000 3.60555128 30 3 5.0000000 3.60555128 50 3 4.3333333 3.78593890 50 3 4.3333333 3.78593890 100 3 11.6666667 4.16333200 100 3 11.6666667 4.16333200Level of -y-Level of -y-xb N Mean Std Dev xb N Mean Std Dev 450 3 5.00000000 4.35889894 450 3 5.00000000 4.35889894 500 3 6.66666667 6.02771377 500 3 6.66666667 6.02771377 600 3 9.33333333 4.93288286 600 3 9.33333333 4.93288286 Level of -y-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四个因素各水平组均值输出结果:第26页,此课件共34页哦第四节第四节 正交试验的优良性正交试验的优良性单因素轮换法单因素轮换法:将其它因素固定在某一水平,每次只改变欲研究因素的水平,依次类推以研究A、B、C三因素三水平试验为例,在研究因素A的效应时,先将因素B、C固定在某一水平如B1、C1;如发现A3较好,固定B2,试验安排在A3B2C1、A3B2C2、A3B2C3;如发现C1较好,试验安排在A3C1B1、A3C1B2、A3C1B3;V如发现B2较好,则确定A3B2C1为最佳分析条件.(请看下图中的七个黑点七个黑点即为上述单因素轮换法安排试验时试验点的分布位置)缺点:实际考察的因素水平仅局限于欲考察的因素水平范围的局部区域,不能全面反映因素影响的情况,且不能考察因素之间的交互效应.全面试验法全面试验法:工作量大,如三因素三水平全面试验要进行33=27次.(请看下图中的27个交叉点个交叉点即为全面试验法安排试验时试验点的分布位置)第27页,此课件共34页哦第28页,此课件共34页哦正交试验法正交试验法正交试验法正交试验法:利用正交性的试验设计方法利用正交性的试验设计方法三因素三水平试验可选用L L9(34 4)正交表正交表,试验只需进行9 9次次.(请看上图中的请看上图中的9 9个红点个红点个红点个红点即为正交试验法安排试即为正交试验法安排试验时试验点的分布位置验时试验点的分布位置)第29页,此课件共34页哦正交性正交性包括:1.均衡分散性均衡分散性:使每个试验点有强烈的代表性正交试验的试验次数(红点总数):全面试验试验次数(交叉点数)=每个平面上红点数 :每个平面上交叉点数=每条棱上红点数 :每条棱上交叉点数即 9/27=3/9=1/3(三因素三水平场合)第30页,此课件共34页哦2.整齐可比性整齐可比性:使正交表的同一列可比较产生”好水平”.上图说明A1的三次试验(红点)中因素B的三个水平和因素C的三个水平各出现一次,A2和A3有同样情况,就是说,对因素A的三个不同水平,受因素B和因素C的水平影响处于对等状态.局部试验的最优点的平均名次局部试验的最优点的平均名次=(N+1)/(n+1)(N (N是全面试验次数是全面试验次数,n是局部试验次数是局部试验次数)如对于正交表如对于正交表L8(27)L8(27)则有(128+1)/(8+1)14,说明这8 8次次试验的好结果在全面试验中的排位大致在前试验的好结果在全面试验中的排位大致在前14名左右.第31页,此课件共34页哦第六节第六节 正交表的灵活安排正交表的灵活安排定量因素的试验用量或定性因素的状态叫做因素的定量因素的试验用量或定性因素的状态叫做因素的水平,在确定参加正交试验的因素后,各种类型的水平,在确定参加正交试验的因素后,各种类型的正交表的格式是固定的,有时候各因素的水平必须灵活处理以适应某张合适的正交表.1.活动水平2.拟水平和后备水平3.配比试验配比试验第32页,此课件共34页哦1.1.活动水平活动水平因素因素硼酸量硼酸量脱水量脱水量反应时间反应时间邻氨基邻氨基用量用量二甲苯二甲苯温度温度苹果酸苹果酸用量用量通二氧化碳时间通二氧化碳时间水平水平水平水平1.5g1.5g1.7g1.7g7ml7ml8.2ml8.2ml7 7小时小时小时小时45g45g50g50g7070606029g29g26g26g0.80.8反应时间反应时间 反应时间反应时间例北京印染厂生产白“的确良”,用到一种增白剂,提高该产品质量第33页,此课件共34页哦2.2.拟水平和后备水平拟水平和后备水平因素因素设备设备催化剂品种催化剂品种原料品种原料品种反应时间反应时间水平水平水平水平水平水平甲甲乙乙(甲)(甲)卯卯寅寅子子长春长春武汉武汉广州广州3030分分4040分分2020分分拟水平:作为水平的设备甲后备水平:催化剂品种丑催化剂品种丑第34页,此课件共34页哦