《均匀设计与分析》课件

均匀试验设计正交设计是利用正交表的均衡分散性和整齐可比性，以较少的实验次数获得基本上能反映全面情况的试验结果的一种优化试验设计方法．为了保证整齐可比和搭配均衡的特点，简化数据处理，实验点应在试验范围内充分地均衡分散，因此试验点不能过少．当欲考察的因素较多，特别是因素水平数较多时，需要的试验次数仍然很多，例如要考察9个水平试验，用正交表安排试验，至少要进行92次试验．均匀试验设计正交设计是利用正交表的均衡分散性和整齐可比性，以1为此，寻找一种适用于多因素、多水平而试验次数更少的试验设计方法是有意义的，我国数学家方开泰和王元等利用数论方法构造了均匀试验设计（uniformdesign）表．如果不考察试验数据的整齐可比性，而让试验点在试验范围内充分地均衡分散，则可以从全面试验中挑选比正交试验设计更少的实验点作为代表进行试验，这种着眼于实验点充分地均衡分散的试验方法，称为均匀试验设计方法．为此，寻找一种适用于多因素、多水平而试验次数更少的试验设计方2《均匀设计与分析》课件32008年度国家科学技术奖励大会日前在北京举行，香港浸会大学荣休教授方开泰和中国科学院数学与系统科学研究院王元院士合作研究逾30年的“均匀试验设计的理论、方法及其应用”，获颁国家自然科学奖二等奖。2008年度国家科学技术奖励大会日前在北京举行，香港浸会大学4均匀设计最先运用在军事工业上（我国导弹设计），后来在石油、化工、生物以及科学计算等高新产业上也获得成功应用。著名汽车品牌福特汽车在开发6汽缸汽车引擎时，便应用了均匀设计，其后该理论更成为福特汽车计算机仿真实验的常规方法。东北制药总厂为了使数理统计方法在工业参数优化中发挥更多的作用，成立了优化技术应用研究室，将各种实用的数学方法在计算机上实现，供科研和生产应用。为此研制出“均匀设计与统计调优软件包”，用这一技术完成十几项科研和生产课题，创百万元以上的经济效益。1993年“均匀设计与统计调优技术应用”通过国家医药管理局的技术鉴定，1994年列入全国医药行业“八五”科技推广项目。均匀设计最先运用在军事工业上（我国导弹设计），后来在石油、化5《均匀设计与分析》课件6和正交试验设计需要正交表一样，均匀试验设计也需要用规格化的表格来安排实验，这种表格称均匀设计表，简称U表（uniform）。通常只列出试验次数为奇数的表，对于偶数次数试验可以用试验次数多一次的奇数表划去最后一行来安排．和正交试验设计需要正交表一样，均匀试验设计也需要用规格化的表7列数（最多能安排的因素个数，包括交互作用、误差等）行数（需要做的试验次数）各因素的水平数（各因素的水平数相等）均匀设计表的代号均匀设计表的记号及含义列数（最多能安排的因素个数，包括交互作用、误差等）行数（需要8均匀试验设计的突出优点是试验工作量很少，特别适用于水平数较多时的试验安排．但它与正交表是不同的，不仅表中各列的地位不平等，而且各因素安排在表中的位置也是不能随便变换的，需根据试验中欲考察的实际因素数，依照附在每一张均匀设计表后的使用表来确定因素所对应的列号．例如用U11(1110）安排2因素11水平的试验，因素安排在第1列与第7列；5因素11水平试验则安排在第1，2，3，5，7列.均匀试验设计的突出优点是试验工作量很少，特别适用于水平数较9均匀设计的特点1）每个因素的每个水平做一次且仅做一次试验。2）任两个因素的试验点点在平面的格子点上，每行每列有且仅有一个试验点。性质1）和2）反映了试验安排的“均衡性”，即对各因素，每个因素的每个水平一视同仁。3）均匀设计表任两列组成的试验方案一般并不等价。均匀设计表的这一性质和正交表有很大的不同，因此，每个均匀设计表必须有一个附加的使用表。4）当因素的水平数增加时，试验数按水平数的增加量在增加。而正交设计当水平增加时，试验数按水平数的平方的比例在增加。由于这个特点，使均匀设计更便于使用。均匀设计的特点1）每个因素的每个水平做一次且仅做一次试验。10试验安排的特点使试验数据失去了整齐可比性，数据一般应采用回归分析法进行分析．由于实验次数较少，试验精度较差，为了提高其精度，可采用试验次数较多的均匀设计表来重新安排因素各水平的试验．试验安排的特点使试验数据失去了整齐可比性，数据一般应采用回归11均匀设计表的构造如均匀表U7(76)，第1次（表中第2行）6个因素分别以1，2，3，4，5和6的水平进行组合试验；下一次试验各因素水平分别在前一次水平基础上加1，2，3，4，5和6，并以7进制进位取余数（当余数为0时，水平号取7）例如在表中第3列的7次试验的水平号分别为3，3＋3=6，6＋3=9→2（取余数），2＋3=5，5＋3=8→1，1＋3=4和4＋3=7，即为3，6，2，5，1，4和7。同理可构造出其他均匀表的因素水平安排．均匀设计表的构造如均匀表U7(76)，第1次（表中第2行）12《均匀设计与分析》课件13例8-10

:阿魏酸的制备阿魏酸是某些药品的主要成分，为了在制备过程中能增加其产量。经过分析研究，挑选出因素和试验区域，为：原料配比:1.0---3.4吡啶总量:10----28反应时间:0.5---3.5确定了每个因素相应的水平数为7。例8-10:阿魏酸的制备14表1.第1步:将试验因素的水平列成下表：表1.第1步:将试验因素的水平列成下表：15例如:表1.2:表1.3:16例如:表1.2:表1.3:1616每个表还有一个使用表，将建议我们如何选择适当的列。其中‘偏差’为均匀性的度量值，数值小的设计表示均匀性好。例如

U7(74)的使用表为：表1.1.4:表1.1.2:每个表还有一个使用表，将建议我们如何选择适当的列。其中‘偏差17第3步:

应用选择的

UD-表,做出试验安排。1.将

x1,x2和

x3放入列1,２和3.x1x2x3

2．用x1的７个水平替代第一列的1到7.1.01.41.82.22.63.03.43.对第二列，第三列做同样的替代.

131.5193.0251.0102.5160.5222.0283.54.完成该设计对应的试验，得到７个结果，将其放入最后一列.表1.1.5:第3步:应用选择的UD-表,做出试验安排。1.将18第4步:

用回归模型匹配数据首先，考虑线性回归模型:dataex;inputx1x2x3y;cards;1.0131.50.3301.4193.00.3661.8251.00.2942.2102.50.4762.6160.50.2093.0222.00.4513.4283.50.482;proc

glm;modely=x1x2x3;run;第4步:用回归模型匹配数据首先，考虑线性回归模型:da19

TheSASSystem23:13Tuesday,October26,20041TheGLMProcedureNumberofobservations7TheSASSystem23:13Tuesday,October26,20042TheGLMProcedureDependentVariable:ySumofSourceDFSquaresMeanSquareFValuePr>FModel30.048770100.016256703.290.1773Error30.014837610.00494587CorrectedTotal60.06360771R-SquareCoeffVarRootMSEyMean0.76673218.876080.0703270.372571SourceDFTypeISSMeanSquareFValuePr>Fx110.010452890.010452892.110.2420x210.000324110.000324110.070.8145x310.037993100.037993107.680.0695SourceDFTypeIIISSMeanSquareFValuePr>Fx110.004542330.004542330.920.4086x210.002195720.002195720.440.5529x310.037993100.037993107.680.0695StandardParameterEstimateErrortValuePr>|t|Intercept0.20236417750.099329342.040.1344x10.03718344160.038800000.960.4086x2-.00344696970.00517333-0.670.5529x30.07694805190.027763022.770.0695TheSASSystem23:13Tu20dataex;inputx1x2x3y;X4=x3*x3;x5=x1*x3;cards;1.0131.50.3301.4193.00.3661.8251.00.2942.2102.50.4762.6160.50.2093.0222.00.4513.4283.50.482;proc

glm;modely=x3x4x5;run;dataex;21

TheSASSystem23:13Tuesday,October26,20049TheREGProcedureModel:MODEL1DependentVariable:yAnalysisofVarianceSumofMeanSourceDFSquaresSquareFValuePr>FModel30.062190.0207343.880.0056Error30.001420.00047244CorrectedTotal60.06361RootMSE0.02174R-Square0.9777DependentMean0.37257AdjR-Sq0.9554CoeffVar5.83398ParameterEstimatesParameterStandardVariableDFEstimateErrortValuePr>|t|Intercept10.062320.034961.780.1727x310.251130.039186.410.0077x41-0.060000.01063-5.640.0110x510.023470.004814.880.0164

22rcoplot(r,rint)作残差分析图（如残差置信区间不包含中心线，则可视为异常点）rcoplot(r,rint)作残差分析图（如残差置信区23第5步:

优化--寻找最佳的因素水平组合表1.1.5的设计是73=343个全面试验的部分实施,

其中最好的试验点是值为Y=48.2%的#7。它不一定是全局最好的。如果想找到满足下式的x1*和

x3*:这里求取max的区域为：第5步:优化--寻找最佳的因素水平组合表1.1.5的设24bint=-0.04890.17360.12640.3758-0.0938-0.02620.00820.0388stats= 0.977743.87860.0056《均匀设计与分析》课件25可以看到，结果分析与正交设计有所不同，一般应采用多元线性回归分析法．指标最佳试验点所对应的试验条件，即使不是全面试验中最好的条件，往往也是接近于全面试验中最佳条件的试验条件．由于均匀试验中试验点在整个试验区域内分布均匀，因此采用多元线性回归对实验数据进行解析以确定最佳实验条件是可行的．可以看到，结果分析与正交设计有所不同，一般应采用多元线性回归26例8-11石墨炉原子吸收分析钯试验因素：灰化温度、灰化时间、原子化温度，原子化时间例8-11石墨炉原子吸收分析钯试验因素：灰化温度、灰化时间27《均匀设计与分析》课件28进一步去掉x4

进一步去掉x429Y=-0.05345-3.045×10-3x2-3.138×10-8x12

+3.530×10-5x22

+3.419×10-8

x32Y=-0.05345-3.045×10-3x2-30最优试验条件将计算得的b阵代入可算得最优条件为 TC=139.6℃、tc=41.2s、Ta=1791.3℃，ta=6.6s 由于该条件并不在设计表内，应再安排试验得出其指标，判断其是否最优。最优试验条件31例8-12.均匀设计法在全光亮镀镍研究中的应用1.均匀设计表的选取 本实验的目的是提高镀层光亮性。经初步研究，取其固定组成为硫酸镍25g/L，次磷酸钠25g/L，醋酸钠25g/L。考察因素为稳定剂，主光亮剂，辅助光亮剂，润湿剂4个因素，每个因素取值范围为t个水平（t为实验次数），4个因素的一次项及二次项各有4项，4项因素间的两两交互作用设有6项，共14项，实验数不能小于14，本实验选用U17(178)表。例8-12.均匀设计法在全光亮镀镍研究中的应用1.均匀32均匀表U17(178)均匀表U17(178)33U17(178)表的使用表本实验为4因素，这4个因素安排在均匀表的1，5，7，8列，去掉U17(178)的最后一行，将实验方案及结果见下表。U17(178)表的使用表本实验为4因素，这4个因素安排在均34《均匀设计与分析》课件352.指标的选择和优化指标是回归方程中的响应函数，在本实验中即是镀件质量。根据我们对镀件的要求，定义一个综合指标z，z的分值由外观评分R，沉积速度评分V，耐腐蚀性评分Q乘以不同的权重构成，z=0.5R+0.2V+0.3Q。R，V，Q的分值分别为0－100。2.指标的选择和优化指标是回归方程中的响应函数，在本实验中即363.实验方法试样为10cm×5cm×0.2cm的低碳钢板，在88－90℃的恒温水浴槽内施镀，镀液pH值控制在4.5-5.0。镀前处理按常规进行，按均匀设计表中确定的组成分别配成16种化学镀液，挂镀法施镀1h，清洗，晾干，对试样进行外观的评定。沉积速度测定：沉积速度，样片增加的重量/样片的面积(g/cm2)耐腐蚀性测定：10％硫酸浸泡24h，根据失重及腐蚀后外观评分3.实验方法试样为10cm×5cm×0.2cm的低碳钢板，在374.结果处理及分析实验结果用计算机处理，主要运用软件为SPSS和Matlab。

4.1建立数学模型及筛选变量 考虑到可能有的数学关系，将各因素的一次项，二次项，两因子间的交互作用项均作为考察对象，回归方程模型为： R=b0+∑bixi+∑bijxixj+∑biixi2(i=1,2,3,4;i≠j)b为各项系数。将给因素的值及综合指标输入计算机，用自后淘汰变量法(backwardselection)进行回归分析和变量筛选，sigF＞0.10的变量被淘汰，最后得到指标与相关组成的回归方程。4.结果处理及分析实验结果用计算机处理，主要运用软件为SPS38Z=86.726+6.555×d－4.554×p2＋1.384×c2＋0.01641×ω2－3.177×p×c＋0.1932×p×ω－0.1209×c×ω－0.3779×d×ω

c为主光亮剂；d为辅助光亮剂；ω为润湿剂；p为稳定剂。4.2对回归方程的优化处理用求条件极值的强约束优化法对回归方程进行优化，用Matlab语言编程，用BFGS拟牛顿(Quasi-Newton)算法及最小二乘法寻优，本实验找到的最优解为：主光亮剂C3.7mg/L，辅助光亮剂HD1.1ml/g，稳定剂0.2mg/L，润湿剂19.7mg/L，乳酸6mol/L。Z=86.726+6.555×d－4.554×p2＋1.3839§8-6

单纯形试验法（simplex）单纯形优化法最初由Erns引入化学研究，它是一种按黑箱方式工作的试验设计方法。正交设计试验是通过在多维空间中均衡地布置试验点，并比较它们的优劣来寻优单纯形试验法是基于在多维空间中的几何图形的变换来寻找试验点，并逐一比较试验结果后逐步搜索寻优的。§8-6

单纯形试验法（simplex）单纯形优化法40这里所说的图形是指的在n维空间中，具有n+1个顶点的图形，例如，两维空间使用三角形，三维空间使用有四个顶点的多面体。这些图形的顶点即是试验所安排的试验点，通过比较各试验点的结果并丢弃最差点而代之以新点，从而再构成新的单纯形，这样就可逐步逼近最优试验点。这里所说的图形是指的在n维空间中，具有n+1个顶点的图形，例41谢谢！谢谢！42均匀试验设计正交设计是利用正交表的均衡分散性和整齐可比性，以较少的实验次数获得基本上能反映全面情况的试验结果的一种优化试验设计方法．为了保证整齐可比和搭配均衡的特点，简化数据处理，实验点应在试验范围内充分地均衡分散，因此试验点不能过少．当欲考察的因素较多，特别是因素水平数较多时，需要的试验次数仍然很多，例如要考察9个水平试验，用正交表安排试验，至少要进行92次试验．均匀试验设计正交设计是利用正交表的均衡分散性和整齐可比性，以43为此，寻找一种适用于多因素、多水平而试验次数更少的试验设计方法是有意义的，我国数学家方开泰和王元等利用数论方法构造了均匀试验设计（uniformdesign）表．如果不考察试验数据的整齐可比性，而让试验点在试验范围内充分地均衡分散，则可以从全面试验中挑选比正交试验设计更少的实验点作为代表进行试验，这种着眼于实验点充分地均衡分散的试验方法，称为均匀试验设计方法．为此，寻找一种适用于多因素、多水平而试验次数更少的试验设计方44《均匀设计与分析》课件452008年度国家科学技术奖励大会日前在北京举行，香港浸会大学荣休教授方开泰和中国科学院数学与系统科学研究院王元院士合作研究逾30年的“均匀试验设计的理论、方法及其应用”，获颁国家自然科学奖二等奖。2008年度国家科学技术奖励大会日前在北京举行，香港浸会大学46均匀设计最先运用在军事工业上（我国导弹设计），后来在石油、化工、生物以及科学计算等高新产业上也获得成功应用。著名汽车品牌福特汽车在开发6汽缸汽车引擎时，便应用了均匀设计，其后该理论更成为福特汽车计算机仿真实验的常规方法。东北制药总厂为了使数理统计方法在工业参数优化中发挥更多的作用，成立了优化技术应用研究室，将各种实用的数学方法在计算机上实现，供科研和生产应用。为此研制出“均匀设计与统计调优软件包”，用这一技术完成十几项科研和生产课题，创百万元以上的经济效益。1993年“均匀设计与统计调优技术应用”通过国家医药管理局的技术鉴定，1994年列入全国医药行业“八五”科技推广项目。均匀设计最先运用在军事工业上（我国导弹设计），后来在石油、化47《均匀设计与分析》课件48和正交试验设计需要正交表一样，均匀试验设计也需要用规格化的表格来安排实验，这种表格称均匀设计表，简称U表（uniform）。通常只列出试验次数为奇数的表，对于偶数次数试验可以用试验次数多一次的奇数表划去最后一行来安排．和正交试验设计需要正交表一样，均匀试验设计也需要用规格化的表49列数（最多能安排的因素个数，包括交互作用、误差等）行数（需要做的试验次数）各因素的水平数（各因素的水平数相等）均匀设计表的代号均匀设计表的记号及含义列数（最多能安排的因素个数，包括交互作用、误差等）行数（需要50均匀试验设计的突出优点是试验工作量很少，特别适用于水平数较多时的试验安排．但它与正交表是不同的，不仅表中各列的地位不平等，而且各因素安排在表中的位置也是不能随便变换的，需根据试验中欲考察的实际因素数，依照附在每一张均匀设计表后的使用表来确定因素所对应的列号．例如用U11(1110）安排2因素11水平的试验，因素安排在第1列与第7列；5因素11水平试验则安排在第1，2，3，5，7列.均匀试验设计的突出优点是试验工作量很少，特别适用于水平数较51均匀设计的特点1）每个因素的每个水平做一次且仅做一次试验。2）任两个因素的试验点点在平面的格子点上，每行每列有且仅有一个试验点。性质1）和2）反映了试验安排的“均衡性”，即对各因素，每个因素的每个水平一视同仁。3）均匀设计表任两列组成的试验方案一般并不等价。均匀设计表的这一性质和正交表有很大的不同，因此，每个均匀设计表必须有一个附加的使用表。4）当因素的水平数增加时，试验数按水平数的增加量在增加。而正交设计当水平增加时，试验数按水平数的平方的比例在增加。由于这个特点，使均匀设计更便于使用。均匀设计的特点1）每个因素的每个水平做一次且仅做一次试验。52试验安排的特点使试验数据失去了整齐可比性，数据一般应采用回归分析法进行分析．由于实验次数较少，试验精度较差，为了提高其精度，可采用试验次数较多的均匀设计表来重新安排因素各水平的试验．试验安排的特点使试验数据失去了整齐可比性，数据一般应采用回归53均匀设计表的构造如均匀表U7(76)，第1次（表中第2行）6个因素分别以1，2，3，4，5和6的水平进行组合试验；下一次试验各因素水平分别在前一次水平基础上加1，2，3，4，5和6，并以7进制进位取余数（当余数为0时，水平号取7）例如在表中第3列的7次试验的水平号分别为3，3＋3=6，6＋3=9→2（取余数），2＋3=5，5＋3=8→1，1＋3=4和4＋3=7，即为3，6，2，5，1，4和7。同理可构造出其他均匀表的因素水平安排．均匀设计表的构造如均匀表U7(76)，第1次（表中第2行）54《均匀设计与分析》课件55例8-10

:阿魏酸的制备阿魏酸是某些药品的主要成分，为了在制备过程中能增加其产量。经过分析研究，挑选出因素和试验区域，为：原料配比:1.0---3.4吡啶总量:10----28反应时间:0.5---3.5确定了每个因素相应的水平数为7。例8-10:阿魏酸的制备56表1.第1步:将试验因素的水平列成下表：表1.第1步:将试验因素的水平列成下表：57例如:表1.2:表1.3:58例如:表1.2:表1.3:1658每个表还有一个使用表，将建议我们如何选择适当的列。其中‘偏差’为均匀性的度量值，数值小的设计表示均匀性好。例如

U7(74)的使用表为：表1.1.4:表1.1.2:每个表还有一个使用表，将建议我们如何选择适当的列。其中‘偏差59第3步:

应用选择的

UD-表,做出试验安排。1.将

x1,x2和

x3放入列1,２和3.x1x2x3

2．用x1的７个水平替代第一列的1到7.1.01.41.82.22.63.03.43.对第二列，第三列做同样的替代.

131.5193.0251.0102.5160.5222.0283.54.完成该设计对应的试验，得到７个结果，将其放入最后一列.表1.1.5:第3步:应用选择的UD-表,做出试验安排。1.将60第4步:

用回归模型匹配数据首先，考虑线性回归模型:dataex;inputx1x2x3y;cards;1.0131.50.3301.4193.00.3661.8251.00.2942.2102.50.4762.6160.50.2093.0222.00.4513.4283.50.482;proc

glm;modely=x1x2x3;run;第4步:用回归模型匹配数据首先，考虑线性回归模型:da61

TheSASSystem23:13Tuesday,October26,20041TheGLMProcedureNumberofobservations7TheSASSystem23:13Tuesday,October26,20042TheGLMProcedureDependentVariable:ySumofSourceDFSquaresMeanSquareFValuePr>FModel30.048770100.016256703.290.1773Error30.014837610.00494587CorrectedTotal60.06360771R-SquareCoeffVarRootMSEyMean0.76673218.876080.0703270.372571SourceDFTypeISSMeanSquareFValuePr>Fx110.010452890.010452892.110.2420x210.000324110.000324110.070.8145x310.037993100.037993107.680.0695SourceDFTypeIIISSMeanSquareFValuePr>Fx110.004542330.004542330.920.4086x210.002195720.002195720.440.5529x310.037993100.037993107.680.0695StandardParameterEstimateErrortValuePr>|t|Intercept0.20236417750.099329342.040.1344x10.03718344160.038800000.960.4086x2-.00344696970.00517333-0.670.5529x30.07694805190.027763022.770.0695TheSASSystem23:13Tu62dataex;inputx1x2x3y;X4=x3*x3;x5=x1*x3;cards;1.0131.50.3301.4193.00.3661.8251.00.2942.2102.50.4762.6160.50.2093.0222.00.4513.4283.50.482;proc

glm;modely=x3x4x5;run;dataex;63

TheSASSystem23:13Tuesday,October26,20049TheREGProcedureModel:MODEL1DependentVariable:yAnalysisofVarianceSumofMeanSourceDFSquaresSquareFValuePr>FModel30.062190.0207343.880.0056Error30.001420.00047244CorrectedTotal60.06361RootMSE0.02174R-Square0.9777DependentMean0.37257AdjR-Sq0.9554CoeffVar5.83398ParameterEstimatesParameterStandardVariableDFEstimateErrortValuePr>|t|Intercept10.062320.034961.780.1727x310.251130.039186.410.0077x41-0.060000.01063-5.640.0110x510.023470.004814.880.0164

64rcoplot(r,rint)作残差分析图（如残差置信区间不包含中心线，则可视为异常点）rcoplot(r,rint)作残差分析图（如残差置信区65第5步:

优化--寻找最佳的因素水平组合表1.1.5的设计是73=343个全面试验的部分实施,

其中最好的试验点是值为Y=48.2%的#7。它不一定是全局最好的。如果想找到满足下式的x1*和

x3*:这里求取max的区域为：第5步:优化--寻找最佳的因素水平组合表1.1.5的设66bint=-0.04890.17360.12640.3758-0.0938-0.02620.00820.0388stats= 0.977743.87860.0056《均匀设计与分析》课件67可以看到，结果分析与正交设计有所不同，一般应采用多元线性回归分析法．指标最佳试验点所对应的试验条件，即使不是全面试验中最好的条件，往往也是接近于全面试验中最佳条件的试验条件．由于均匀试验中试验点在整个试验区域内分布均匀，因此采用多元线性回归对实验数据进行解析以确定最佳实验条件是可行的．可以看到，结果分析与正交设计有所不同，一般应采用多元线性回归68例8-11石墨炉原子吸收分析钯试验因素：灰化温度、灰化时间、原子化温度，原子化时间例8-11石墨炉原子吸收分析钯试验因素：灰化温度、灰化时间69《均匀设计与分析》课件70进一步去掉x4

进一步去掉x471Y=-0.05345-3.045×10-3x2-3.138×10-8x12

+3.530×10-5x22

+3.419×10-8

x32Y=-0.05345-3.045×10-3x2-72最优试验条件