页数:60页文件大小:1.82 MB来源: 食品试验设计优化ppt课件
本课件通过阿魏酸制备等实例演示均匀设计的完整步骤:列出因素水平表、选择均匀设计表、安排试验方案、用回归模型匹配数据并优化因素水平组合,并介绍混合型水平均匀设计、含定性因素的伪变量处理与混料配方均匀设计。
Uniform Design andits Application
沐浴在改革开发的阳光下,神州大地生机盎然,新生事物层出不穷。在科教兴国建设社会主义的过程中,人们所熟悉的那些传统的试验设计方法(如对比试验设计、全面试验设计、正交试验设计等),已不能充分满足快节奏高效率的要求。新时期呼唤新思维,新方法。中国科学家巧妙的将“数论方法”和“统计试验设计”相结合,发明了一种全新的试验设计方法,这就是均匀设计法。
均匀设计法诞生于1978年。由中国著名数学家方开泰教授和王元院士合作共同发明。前言
正交设计可使试验点“均匀分散、整齐可比”,为保证“整齐可比性”,使试验设计的均匀性受到了一定限制,使试验点的代表性还不够强,试验次数不能充分地少。均匀设计是另一种部分实施的试验设计方法。它可以用较少的试验次数,安排多因素、多水平的析因试验,是在均匀性的度量下最好的析因试验设计方法。
它可以使试验点在试验范围内充分地均匀分散,不仅可大大减少试验点,而且仍能得到反映试验体系主要特征的试验结果。
均匀设计方法
下面通过制药工业中的一个实例来说明均匀试验设计方法。例1.1:阿魏酸的制备根据试验目的,确定以阿魏酸产量作为试验指标Y。阿魏酸是某些药品的主要成分,在制备过程中,我们想提高阿魏酸产量。
全面交叉试验要N=73=343次,太多了。建议使用均匀设计。查阅均匀设计表。经过资料查阅,分析研究,选出影响阿魏酸产量的试验因素,确定试验因素水平为:原料配比:1.0---3.4吡啶总量:10----28反应时间:0.5---3.5确定每个因素相应的水平数为7。
如何安排试验?“方开泰,均匀设计与均匀设计表,科学出版社(1994)”之附表1网络地址:http://www.math.hkbu.edu.hk/UniformDesing也可以浏览如下网页:
第1步:列出试验因素水平表表1.1.1试验因素水平表
第2步:选择相应的均匀设计表均匀设计表格式见下,其含义为:Un(qs)均匀设计试验总次数因素水平数因素数
本页源文件仅提取到少量文字(多为页码或符号残片),建议对照原课件图片补充本页要点说明。
每个均匀设计表都有一个使用表,它将建议我们如何选择适当的列安排试验因素,进行试验设计,这样可以减少“试验偏差”。其中‘偏差’为均匀性的度量值,数值小的设计表示均匀性好。例如U7 (74)的使用表为:表1.1.4:表1.1.2:
第3步:应用选择的UD-表安排试验,设计试验方案1. 将x1, x2和x3放入均匀设计表的1,2和3列;x1 x2 x32.用x1的7个水平(值)替代第一列的1到7;
| 试验号 | x1 | x2 | x3 | 试验结果Y |
|---|---|---|---|---|
| 1 | 1.0 | 13 | 1.5 | 源文未给出 |
| 2 | 1.4 | 19 | 3.0 | 源文未给出 |
| 3 | 1.8 | 25 | 1.0 | 源文未给出 |
| 4 | 2.2 | 10 | 2.5 | 源文未给出 |
| 5 | 2.6 | 16 | 0.5 | 源文未给出 |
| 6 | 3.0 | 22 | 2.0 | 源文未给出 |
| 7 | 3.4 | 28 | 3.5 | 源文未给出 |
源文仅说明“得到7个结果,将其放入最后一列”,未给出具体数值,建议对照原课件图片补全。
对第二列,第三列做同样的替代;13 1.519 3.025 1.010 2.516 0.522 2.028 3.54. 按设计的方案进行试验,得到7个结果,将其放入最后一列。表1.1.5:试验方案
第4步:用回归模型匹配数据首先,考虑线性回归模型:这个结果与人们的经验不符。使用回归分析中变量筛选的方法,比如‘向后法’,得到推荐的模型为:
然后,我们尝试用二次回归模型来匹配这些数据:使用‘向前’的变量选择法,我们发现适宜的模型:14表1.1.6: 方差分析(ANOVA)表
状态是正常的,所以模型(1.1.4)是可接受的。图1.1.1:模型中的三项,在5%的水平下都是显著的。
第5步:优化--寻找最佳的因素水平组合表1.1.5的设计是73=343个全面试验的部分实施,其中最好的试验点是值为Y= 48.2%的#7。它不一定是全局最好的。人们想找到满足下式的x1*和x3*:这里求取max的区域为:
x1x3的回归系数是正的,x3的回归系数也是正的, x1* = 3.4。在x3*= 2.7575达到最大值。图1.1.3等值线图(x1*,x3*)在x1*= 3.4和x3*= 2.7575处估计响应的最大值是51.85% 。它比7个试验点的最好值48.2%还大。18
讨论:因素x2没有给响应Y予显著的贡献,我们可以选x2为其中点x2= 19 ml.求出的x1*= 3.4 在边界上, 我们需要扩大x1的试验上限。在x1= 3.4和x3= 2.7575的邻域,追加一些试验是必要的。在第5步,一些优化算法是很有用的。
混合型水平的均匀设计试验中各因素若有不同水平数,比如,其水平数分别为q1,…,qk。这时应使用相应的混合均匀设计表。见“方开泰,均匀设计与均匀设计表,科学出版(1994).”之附表2
每个混合水平表有一个记号,含义为:Un(q1 ×… ×qk)均匀设计试验次数各定量因素之水平数定量因素的最大数
下表是一个混合水平均匀设计表:此均匀设计表试验总数n为12,用它可以安排水平数为6、4、3的因素各一个。
U12(624)此表也是混合水平均匀设计表。23它的试验数n为12。可以安排二个6水平因素和一个4水平因素的设计。
混合因素试验设计
考虑4个因素:平均施肥量X,分为12个水平(70,74,78,82,86,90,94,98,102,106,110,114);种子播种前浸种时间T,分为6个水平(1,2,3,4,5,6);土壤类型B,分4种B1,B2,B3,B4;种子品种A,分3个A1,A2,A3;对某农作物产量的影响。可以看出前两个为定量因素,后两个为定性因素。
例2.1:在农业试验中如何进行试验安排?
混合型因素混合型水平的均匀设计一般情况下试验中既有定量型连续变化因素,又有定性型状态变化因素。假设有k个定量因素X1,…,Xk;这k个因素可化为k个连续变量,其水平数分别为q1,…,qk。又有t个定性因素G1,…,Gt,这t个定性因素分别有d1,…,dt个状态。可以使用“拟水平法”,或用优化方法计算,求出相应的均匀设计表。
混合因素混合水平表有如下的记号和含义:Un(q1×… ×qk×d1×…×dt)均匀设计试验次数各定性因素之水平数定性因素的最大数各定量因素之水平数定量因素的最大数
U12(12×6×4×32×22)1 2 3 4 5 6 7 1 1 1 1 2 3 1 2 2 2 2 2 3 2 2 1 3 3 3 3 2 1 1 2 4 4 4 4 3 1 2 1 5 5 5 1 1 2 2 2 6 6 6 2 3 2 1 1 7 7 1 3 1 1 1 1 8 8 2 4 3 3 2 1 9 9 3 1 1 3 2 2 10 10 4 2 2 2 1 2 11 11 5 3 1 1 1 1 12 12 6 4 2 3 2 2例:12次试验。
可以安排2个水平数为12、6的定量因素,以及总数为5的一个水平为4、两个水平为3和两个水平为2的定性因素的设计。
U12(12×6×4×3)表2.1.1选混合均匀设计表2.1.1安排此试验第1列安排平均施肥量X,分为12个水平;第2列安排种子播种前浸种时间T,分为6个水平;第3列安排土壤类型B,分4种B1,B2,B3,B4;第4列安排种子品种A,分3个A1,A2,A3。
试验安排及结果如表2.1.2
为了进行分析,我们引进5个‘伪变量’。它们的记号和取值如下:它们和X、T一起进行回归分析。B因素的A因素的回归方程如下:
本页源文件仅提取到少量文字(多为页码或符号残片),建议对照原课件图片补充本页要点说明。
不显著。需进一步考虑高阶回归项。若我们考虑除主效应外,再多考虑一个2次效应和一个交互效应。这时回归方程化为解得回归系数的最小二乘估计及其R和F值为:
本页源文件仅提取到少量文字(多为页码或符号残片),建议对照原课件图片补充本页要点说明。
解得非常显著回归系数的最小二乘估计及其R和F值为:
方程为:其中1.含变量x 的两项与其它是分离的(即可加的),最大值点在x=100.127 。2.含变量z41z42的两项与其它是分离的,最大值点在z41=0 z42=0,即品种3为好。
3.含变量Tz31z32 z33 的四项与其它是分离的,最大值点可能在z31=1 z32=0 z33=0 类型为1,T=6或z31=0 z32=1 z33=0 类型为2,T=6比较后知道为后者。
所以得到最佳状态组合为施肥量X=100.127,浸种时间T=6,土壤类型B取2,种子品种A取3,此时最大值估计为
一、表的选择,因素及水平的安排若试验中有k个定量因素和t个定性因素时,我们从混合型均匀设计表中选出带有s=k+t列的Un(q1×…×qk×d1×…×dt)表。这里要求n≥k+d+1,其中d=(d1+…+dt-t). 为了给误差留下自由度,其中的n最好不取等号。表中前k列对应k个连续变量,表中后t列可安排定性因素。
安排n个试验,得到n个结果y1,y2,…,yn。下面综述应注意的事项:
为了分析,首先要将定性因素之状态,依照伪变量法,将第i个因素分别化成(di-1)个相对独立的n维伪变量Zi1,Zi2,…,Zi(di-1)。将这总共d=(d1+…dt-t)个伪变量与相应的k个连续变量X1,…,Xk一起进行建模分析。为了保证主效应不蜕化,要对混合型均匀设计表进行挑选。
二、试验结果的回归建模分析如果不理想,则首先考察它们的一阶回归模型:再考虑一些交互效应,和一些连续变量的高次效应。显然最多可考虑的附加效应数为m个,这里m≤n-(k+d-2)
值得指出的是,由于Zij*Zij=Zij,因此无需考虑伪变量的高阶效应,只考虑连续变量的高次效应即可. 又因为Zij1*Zij2=0,j1≠j2时,因此也无需考虑同一状态因素内的伪变量间的交互效应。只有i1≠i2时,才有可能使Zi1j1*Zi2j2≠0,即不同状态因素间的交互效应可能要考虑.。此外,不要忘记考虑连续变量与伪变量的交互效应。
至于三个以上的状态因素间的交互效应项Zi1j1*Zi2j2*Zi3j3≠0的可能性就更少了。
混料配方试验
许多产品都是混合多种成分在一起形成的。面粉水糖蔬菜汁椰子汁盐发酵粉乳酸钙咖啡粉香料色素咖啡面包怎样确定各种成分的比例呢?经验试验混料试验混料配方均匀设计43
有s个因素: X1, , Xs满足Xi0, i= 1, , s和X1 + + Xs= 1.试验区域为单纯形Ts= {(x1, , xs): xi0, i= 1, , s, x1 + + xs= 1.}单纯形格子点设计(Scheffe, 1958).单纯形重心设计(Scheffe, 1963).轴设计(Cornell, 1975)人们提出了许多混料设计方法,如
例如, 成分数s = 3单纯形格子点设计单纯形重心设计d轴设计这些设计的全面评价请参考:Cornell, J. A. (1990). Experiments with Mixtures: Designs, Models and the Analysis of Mixture Data.Wiley, New York.45
混料均匀设计上述设计的弱点:许多点在Ts的边界上;给用户设计的选择不多。混料均匀设计是要寻找在Ts上均匀散布的试验点。问题:怎样设计这些试验点呢?变换方法46
给定s-1维单位立方体Cs-1上的均匀设计,且用{Ck= (ck1, ,ck,s-1), k= 1, ,n} 表示。则进行下列必要的变换:(3.1.1){xk= (xk1, ,xks), k= 1, ,n} 是Ts.上的均匀设计。47变换方法
例3.1构造T3上带有11个(配方)试验点的均匀设计。假设我们选用U11(112) 和相关的Ck, k= 1, ,11:48变换公式(4.1) 现在成为:(3.1.2)
用这个变换公式, 正方形[0,1]2上的均匀设计Ck= (ck1, ck2), k= 1, ,11 导出T3上的均匀设计Xk= (xk1, xk2, xk3), k= 1, ,11 如下:49
区域T3是一个边长为的等边三角形,用V2表示。111T3x1x2x3T3可以证明:V2上的任何点(z1, z2) 到V2的三条边之距离d1, d2和d3,满足d1+d2+d3= 1.d1d2d3因此,V2上任何点(z1, z2) 都对应一个T3上的点(x1, x2, x3),如果我们像这样在V2上建立一个新坐标系统的话。x1x2x350
给定点(x1, x2, x3),计算点(z1, z2)的公式是:图3.1.2ac1c2图3.1 .2b51
-4-均匀设计软件52介绍
均匀设计软件均匀设计软件有中、英文两个版本。该软件中列举了许多均匀的设计表,并给出了数据分析方法。53
本页源文件仅提取到少量文字(多为页码或符号残片),建议对照原课件图片补充本页要点说明。
程序设计者杜明亮和方法指导者方开泰教授在一起
均匀设计软件可用于*与试验设计相关的大学本科或研究生课程*自然科学研究的试验设计*工业试验和国防科研*系统工程、仿真试验
(一)设计•均匀设计表-用好格子点法生成-用拉丁方生成-用优化方法生成•带拟水平的均匀设计表•带约束的配方设计•无约束的配方设计
(二)数据分析•建模-简单线性模型-二次模型-自选模型•选择变量-前进法-逐步回归法-最优子集法-自选变量•统计诊断-残差点图、正态Q-Q点图-偏回归点图-拟合比较图-等高线图-三维图•综合分析-多响应模型分析-优化-预报-相关系数-各类诊断点图
小结我们介绍了有关均匀设计的一些知识:均匀设计表的构造和用法;介绍了有关均匀设计软件的一些内容。我们强调的是正确使用均匀设计表。即:能确定试验目标,能找出影响因素及其变化范围,合理确定水平数及其值,正确安排试验,对试验结果进行适当的分析,得出恰当的认识。幸好,均匀设计分会研制的“均匀设计软件V3.0”可以帮助你完成这些工作。想了解原理者请接看下面的附录。
主要参考文献Uniform Design and its Application,PPT,kaitai,Fang含有定性因素的均匀设计,PPT, 王柱均匀设计八讲,拷贝,方开泰均匀设计与均匀设计表,书,方开泰均匀设计与正交设计的关联和比较,文章,方开泰,马长兴含有定性因素均匀设计均匀性的度量,文章,王柱,方开泰均匀设计理论及其应用研讨会,论文集均匀设计论文选,第一集均匀设计V3.0,软件,方开泰,杜明亮均匀设计论文选,第二集
下载地址:{{downurl}}
{{{a}}}
{{#b}}{{{b}}}
{{/b}} {{/title}} {{#nolog}}登陆注册
{{/nolog}} {{#djsx}}点我刷新
{{/djsx}} {{#gocz}} {{/gocz}} {{#ffqr}} {{/ffqr}}