试验设计与数据分析幻灯片.ppt
《试验设计与数据分析幻灯片.ppt》由会员分享,可在线阅读,更多相关《试验设计与数据分析幻灯片.ppt(79页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、试验设计与数据分析第1页,共79页,编辑于2022年,星期二目 录 第一章第一章 绪论绪论 第二章第二章 常用常用统计统计分布分布 第三章第三章 参数参数估估计计 第四章第四章 假设检验假设检验 第五章第五章 方差分析方差分析 第六章第六章 试验设计试验设计 第七章第七章 回归分析回归分析 第八章第八章 常用统计软件常用统计软件第2页,共79页,编辑于2022年,星期二第六章 试验设计第3页,共79页,编辑于2022年,星期二主要内容6.0 简单的试验设计技术6.1 正交试验设计6.2 响应曲面试验设计6.3 均匀试验设计第4页,共79页,编辑于2022年,星期二6.3 均匀试验设计第5页,共
2、79页,编辑于2022年,星期二引子:黎元生留学加拿大随笔我于1996年受国家教委公派去加拿大滑铁卢大学化工系做访问学者,从事重油乳液的现场制氢破乳加氢改质课题研究工作。加拿大有很多重油,开采过程中严重乳化;加拿大又有全世界最大的沥青砂开采工业,沥青砂在开采和水蒸气抽提过程中也产生大量的沥青乳液。这些乳状液不仅破乳困难,而且破乳脱水后还需再加氢处理才能作为合成原油出售。指导教授想在乳液中通入一氧化碳,在催化剂的作用下使一氧化碳和乳液中的水反应生成氢,氢再与重油或沥青中的含杂原子化合物反应,起到一步过程既破乳脱水,又对重油或沥青加氢处理的作用。这无疑是个好想法。在我去之前的10年中,陆续已有好几
3、个研究生、博士后和访问学者在这个实验室做过这个课题了。第6页,共79页,编辑于2022年,星期二 他们的研究方法是用含硫模型化合物和溶剂与水混合,然后在高他们的研究方法是用含硫模型化合物和溶剂与水混合,然后在高压反应釜中通一氧化碳反应,考察脱硫效果。由于重点放在考察脱硫压反应釜中通一氧化碳反应,考察脱硫效果。由于重点放在考察脱硫上,他们并没有用真正意义上的乳状液做过试验。听说我研究过沥青上,他们并没有用真正意义上的乳状液做过试验。听说我研究过沥青乳液,教授给我的任务就是制备出稳定的含模型化合物的甲苯乳液,教授给我的任务就是制备出稳定的含模型化合物的甲苯水乳水乳液。我以前并没有研究过轻油的乳化,
4、当时心里一点底也没有,但我液。我以前并没有研究过轻油的乳化,当时心里一点底也没有,但我仍信心十足地答应了下来。仍信心十足地答应了下来。我一面订购乳化剂,一面把从国内带去的我一面订购乳化剂,一面把从国内带去的“均匀设计与统计均匀设计与统计调优软件包调优软件包”安装在实验室的微机里,当然在安装之前要经过系安装在实验室的微机里,当然在安装之前要经过系里的电脑管理员对我的软件进行查毒,并把试验方案设计好。由里的电脑管理员对我的软件进行查毒,并把试验方案设计好。由于心里没底,我计划从于心里没底,我计划从5 5种乳化剂中筛选乳化配方,加上油和水的种乳化剂中筛选乳化配方,加上油和水的比例,一共是比例,一共是
5、6 6个自变量。同时我又设计出一种快速评价乳液稳定个自变量。同时我又设计出一种快速评价乳液稳定性的方法,将稳定性试验结果作为考察变量。订的乳化剂来了以性的方法,将稳定性试验结果作为考察变量。订的乳化剂来了以后,我就开始了紧张的试验和评价工作。后,我就开始了紧张的试验和评价工作。第7页,共79页,编辑于2022年,星期二 第一轮试验安排做12个样,评价以后进行回归处理,从中剔除两个对乳化影响不大的乳化剂,再安排第二轮7次试验。在第二轮试验中就出现了稳定性较好的样品。第三轮试验下来,整个稳定区间就出来了。将乳化剂加入量少而又能得到稳定乳状液的配方算出,验证之,又存放两天观察,得到了看起来像雪花膏一
6、样的雪白的含苯并噻吩的甲苯水乳化液。又按对水含量变化的要求,制备出从10 25不同含水量的稳定乳液。做完这些以后,我又观察三天,确信乳液稳定后,将乳化条件、配方变化和稳定性变化关系图整理出来,然后向指导教授汇报。第8页,共79页,编辑于2022年,星期二 指导教授看到我在两个星期内就拿出了雪白又细腻均匀的乳液样品,而且还有配方变化后的稳定区间图,简直不敢相信这是事实。当得知我的“秘密武器”后,又让我给详细介绍和解释软件的使用方法和功能以及均匀设计的数学依据。由于我不能用英语将均匀设计的数学原理讲明白,他又派题目组内一位曾在数学院修过三门研究生课程的数学功底很深的博士生专门去数学院的统计和优化系
7、请教。尽管他们没能在数学院得到满意的解释,但由于亲眼见到均匀设计和统计调优能快速解决问题,还是对它产生了极大的兴趣。接下来又让我用需要加氢脱硫的直馏柴油做成乳化液,由于不需要新订乳化剂,又有了甲苯的经验,一个星期我就拿出了稳定的乳液样品和配方数据。这样,到滑铁卢的第一个月,我就得到了同事风趣相送的“Double E”(Emulsion Expert)的外号。(本节资料来自互联网:黎元生 留学加拿大随笔 1997年6月)第9页,共79页,编辑于2022年,星期二前言均匀设计是一种试验设计方法。它可以用较少的试验次数,安排多因素、多水平的试验,是在均匀性的度量下最好的试验设计方法。第10页,共79
8、页,编辑于2022年,星期二q所有的试验设计方法本质上就是在试验的范围内给出挑选代表点的方法。q正交设计是根据正交性准则来挑选代表点,使得这些点能反映试验范围内各因素和试验指标的关系。q正交设计在挑选代表点时有两个特点:均匀分散,整齐可比。第11页,共79页,编辑于2022年,星期二q“均匀分散”使试验点均衡地布在试验范围内,让每个试验点有充分的代表性,。q“整齐可比”使试验结果的分析十分方便,易于估计各因素的主效应和部分交互效应,从而可分析各因素对指标的影响大小和变化规律。q为了照顾“整齐可比”,它的试验点并没有能做到充分“均匀分散”;为了达到“整齐可比”,试验点的数目就必须比较多。第12页
9、,共79页,编辑于2022年,星期二q若在一项试验中有s 个因素,每个因素各有q 水平,用正交试验安排试验,则至少要作个q2试验,当q较大时,将更大,使实验工作者望而生畏。例如,当 q=12 时,q2=144,对大多数实际问题,要求做144 次试验是太多了!q每一个方法都有其局限性,正交试验也不例外,它只宜于用于水平数不多的试验中。第13页,共79页,编辑于2022年,星期二q为了保证“整齐可比”的特点,正交设计必须至少要求做q2次试验。若要减少试验的数目,只有去掉整齐可比的要求。q均匀设计就是只考虑试验点在试验范围内均匀散布的一种试验设计方法。第14页,共79页,编辑于2022年,星期二中国
10、科学家巧妙的将“数论方法”和“统计试验设计”相结合,发明了均匀设计法。均匀设计法诞生於年。由中国著名数学家方开泰方开泰教授和王元王元院士合作共同发明。第15页,共79页,编辑于2022年,星期二华罗庚王元第16页,共79页,编辑于2022年,星期二第17页,共79页,编辑于2022年,星期二第18页,共79页,编辑于2022年,星期二q均匀设计只考虑试验点在试验范围内充分“均匀散布”而不考虑“整齐可比”,因此它的试验布点的均匀性会比正交设计试验点的均匀性更好,使试验点具有更好的代表性。q由于这种方法不再考虑正交设计中为“整齐可比”而设置的实验点,因而大大减少了试验次数,这是它与正交试验设计法的
11、最大不同之处。第19页,共79页,编辑于2022年,星期二均匀设计的特点第20页,共79页,编辑于2022年,星期二1)均均匀匀设设计计的的最最大大特特点点是是试试验验次次数数等等于于因因素素的最大水平数,而不是平方的关系的最大水平数,而不是平方的关系。如如当当水水平平数数从从9水水平平增增加加到到10水水平平时时,试试验验数数n 也也从从9增增加加到到10。而而正正交交设设计计当当水水平平增增加加时时,试试验验数数按按水水平平数数的的平平方方的的比比例例在在增增加加;当当水水平数从平数从9到到10时,试验数将从时,试验数将从81增加到增加到100。由于这个特点,使均匀设计更便于使用。由于这个
12、特点,使均匀设计更便于使用。第21页,共79页,编辑于2022年,星期二2)每个因素的每个水平做一次且仅做一次试验。)每个因素的每个水平做一次且仅做一次试验。3)任任两两个个因因素素的的试试验验,画画在在平平面面的的格格子子点点上上,每行每列有且仅有一个试验点。每行每列有且仅有一个试验点。特特点点2)和和3)反反映映了了试试验验安安排排的的“均均衡衡性性”,即即对对各各因因素素,每每个个因因素素的每个水平一视同仁。的每个水平一视同仁。第22页,共79页,编辑于2022年,星期二5)均均匀匀设设计计表表任任两两列列组组成成的的试试验验方方案案一一般般并并不等价。不等价。例例如如用用U6(66)的
13、的1,3 和和1,6列列分分别别画画图图,得得图图8(a)和和图图8(b)。我我们们看看到到,(a)的的点点散散布布比比较均匀,而较均匀,而(b)的点散布并不均匀。的点散布并不均匀。第23页,共79页,编辑于2022年,星期二均匀设计表的这一性质和正交表有很大的均匀设计表的这一性质和正交表有很大的不同,因此,每个均匀设计表必须有一个不同,因此,每个均匀设计表必须有一个附加的附加的使用表使用表。第24页,共79页,编辑于2022年,星期二4)把把奇奇数数表表划划去去最最后后一一行行就就得得到到比比原原奇奇数数表表少少一一个个水水平平的的偶偶数数表表,相相应应地地,试试验验次次数数也也少少一,而使
14、用表不变。一,而使用表不变。例例如如,把把U7(76)划划去去最最后后一一行行即即得得U6(66)。)。第25页,共79页,编辑于2022年,星期二6)由由于于均均匀匀设设计计不不再再考考虑虑正正交交试试验验的的整整齐齐可可比比性性,因因此此其其试试验验安安排排既既不不能能考考虑虑交交互互作作用用,也不能估计试验误差;也不能估计试验误差;试试验验结结果果的的分分析析只只能能采采用用直直观观分分析析法法和和回回归归分分析析方方法法,根根据据回回归归系系数数的的绝绝对对值值大大小小,得得出出试试验验因因素素对对指指标标影影响响的的主主次次顺顺序序;根根据据方方程程极值点得出最佳工艺条件。极值点得出
15、最佳工艺条件。第26页,共79页,编辑于2022年,星期二均匀设计及其应用均匀设计及其应用第27页,共79页,编辑于2022年,星期二我们通过制药工业中的一个实例,来看均匀设计表的使用方法。这就是说以阿魏酸的产量作为目标 Y。阿魏酸是某些药品的主要成分,在制备过程中,我们想增加其产量。例例例例6.3-16.3-16.3-16.3-1 :阿魏酸的制备阿魏酸的制备第28页,共79页,编辑于2022年,星期二全面交叉试验要N=73=343次,太多了。建议使用均匀设计。有现成的均匀设计表,提供使用。参见:经过分析研究,挑选出因素和试验区域,为原料配比:1.0-3.4吡啶总量:10-28反应时间:0.5
16、-3.5确定了每个因素相应的水平数为7。如何安排试验呢?“方开泰,均匀设计与均匀设计表,科学出版社(1994).”之附表1网络地址:http:/www.math.hkbu.edu.hk/UniformDesing也可以浏览如下网页第29页,共79页,编辑于2022年,星期二第第1步步:将试验因素的水平列成下表:表表 6.3.1:第30页,共79页,编辑于2022年,星期二第第2步步:选择相应的均匀设计表.每个均匀设计表有一个记号,它有如下的含义:Un(qs)均匀设计试验次数水平数因素的最大数第31页,共79页,编辑于2022年,星期二例如:表表 6.3.2:表表 6.3.3:第32页,共79页
17、,编辑于2022年,星期二每个表还有一个使用表,将建议我们如何选择适当的列。其中偏差为均匀性的度量值,数值小的设计表示均匀性好。例如U7(74)的使用表为,表表 6.3.4:表表6.3.2:第33页,共79页,编辑于2022年,星期二第第3步步:应用选择的UD-表,做出试验安排。1.将x1,x2和x3放入列1,和3.x1 x2 x3 2用x1的个水平替代第一列的1到7.1.01.41.82.22.63.03.43.对第二列,第三列做同样的替代.131.5193.0251.0102.5160.5222.0283.54.完成该设计对应的试验,得到个结果,将其放入最后一列.表表 6.3.5:第34页
18、,共79页,编辑于2022年,星期二第第 4步步:直观分析从表 6.3.5中试验数据可见,第7号试验的指标值最大,第7号试验对应的条件即为较优的工艺条件,即:原料配比3.4,吡啶总量28,反应时间3.5。试验区域,为试验区域,为 原原料配比料配比:1.0-3.4:1.0-3.4吡啶吡啶总量总量:10-28:10-28反应时反应时间间:0.5-3.5:0.5-3.5这些条件都是试验条件的上限,可见,还需要进行进一步的试验以寻求更佳的工艺参数。第35页,共79页,编辑于2022年,星期二第第 5步步:用回归模型匹配数据首先,考虑线性回归模型:Regression Analysis:y versus
19、 x1,x2,x3The regression equation isy=0.202+0.0372 x1-0.00345 x2+0.0769 x3Predictor Coef SE Coef T PConstant 0.20236 0.09933 2.04 0.134x1 0.03718 0.03880 0.96 0.409x2 -0.003447 0.005173 -0.67 0.553x3 0.07695 0.02776 2.77 0.069S=0.07033 R-Sq=76.7%R-Sq(adj)=53.3%Analysis of VarianceSource DF SS MS F PR
20、egression 3 0.048770 0.016257 3.29 0.177Residual Error 3 0.014838 0.004946Total 6 0.063608第36页,共79页,编辑于2022年,星期二线性回归效果不佳,可能存在非线性影响,用逐步回归法拟合非线性方程:Stepwise Regression:y versus x1,x2,x3,x11,x22,x33,x12,x13,x23Stepwise Regression:y versus x1,x2,x3,x11,x22,x33,x12,x13,x23 Alpha-to-Enter:0.15 Alpha-to-Rem
21、ove:0.15 Response is y on 9 predictors,with N=7 Step 1Constant 0.2141x3 0.079T-Value 3.34P-Value 0.021S 0.0627R-Sq 69.06R-Sq(adj)62.87拟合效果不好,包括的自拟合效果不好,包括的自变量太少。增大变量太少。增大Alpha-to-到0.3第37页,共79页,编辑于2022年,星期二 Step 1 2 3 4 5 Step 1 2 3 4 5Constant 0.21414 0.10457 0.06232 0.08483 0.06689Constant 0.21414
22、0.10457 0.06232 0.08483 0.06689x3 0.0792 0.2253 0.2511 0.2318 0.2400 x3 0.0792 0.2253 0.2511 0.2318 0.2400T-Value 3.34 2.24 6.41 11.47 25.73T-Value 3.34 2.24 6.41 11.47 25.73P-Value 0.021 0.089 0.008 0.008 0.025P-Value 0.021 0.089 0.008 0.008 0.025x33 -0.0365 -0.0600 -0.0503 -0.0464x33 -0.0365 -0.06
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 试验 设计 数据 分析 幻灯片
限制150内