响应面法及软件中文教程ppt课件

响应曲面法（响应曲面法（RSM）响响应应曲面法（曲面法（RSM）学习目标学习目标描述为何使用描述为何使用RSM及什么是及什么是RSM解释响应曲面法设计的常用类型解释响应曲面法设计的常用类型用用minitab实施实施RSM方法方法掌握掌握RSM设计数据分析设计数据分析了解最快上升路线法了解最快上升路线法学学习习目目标标RSM之起源与背景英国学者英国学者Box&Wilson（1951年）正式提出响年）正式提出响应曲面方法论应曲面方法论目的：探究多个输入变量与化学制程产出值之目的：探究多个输入变量与化学制程产出值之间关系。间关系。在实验设计规划范围内，如何寻找实验因子最在实验设计规划范围内，如何寻找实验因子最佳组合，以达到最佳反应值。佳组合，以达到最佳反应值。系列化实验的最佳规划。系列化实验的最佳规划。Minitab使分析变成更容易。使分析变成更容易。RSM之起源与背景英国学者之起源与背景英国学者Box&Wilson（1951年）年）What Is RSM?什么是响应面方法（RSM）When doing DOE to maximize yield,which plot do you prefer to see?Why?当实施当实施DOEDOE把良率提到最高，你希望看到那个图？为什么？把良率提到最高，你希望看到那个图？为什么？Plot APlot BWhat Is RSM?什么是响什么是响应应面方法（面方法（RSM）WhenWhat Is RSM?什么是响应面方法（RSM）This plot indicates there is opportunity for higher yield.此图显示良率还有再提高的机会此图显示良率还有再提高的机会良率温度时间What Is RSM?什么是响什么是响应应面方法（面方法（RSM）ThisWhat is RSM?什么是响应面方法（RSM）?Yield良率良率Temp温度温度Time时间时间Optimal Area（Highest Yield）最佳区域（最高良率）最佳区域（最高良率）What is RSM?什么是响什么是响应应面方法（面方法（RSM）?YieWhat is RSM?什么是响应面方法（RSM）?RSMRSM有如正在爬山而看不见山顶。有如正在爬山而看不见山顶。What is RSM?什么是响什么是响应应面方法（面方法（RSM）?RSMWhat is RSM?什么是响应面方法（RSM）?当到达山顶时，用当到达山顶时，用RSMRSM方法对周围区域进行勘查。方法对周围区域进行勘查。What is RSM?什么是响什么是响应应面方法（面方法（RSM）?当到达当到达What is RSM?什么是响应面方法（RSM）?然后对过程制订规格界限然后对过程制订规格界限What is RSM?什么是响什么是响应应面方法（面方法（RSM）?然后然后对对Path of Steepest Ascent 最陡的上升路线How can I move to the top the fastest?我怎样能更快到达山顶？我怎样能更快到达山顶？良率温度时间Path of Steepest Ascent 最陡的上最陡的上Path of Steepest Ascent 最陡的上升路线90Path of steepest ascent最陡上升路线最陡上升路线Path of Steepest Ascent 最陡的上最陡的上Path of Steepest Ascent 最陡的上升路线Optimum 最佳条件最佳条件Path of Steepest Ascent 最陡的上最陡的上RSM的使用时机寻找因子参数设定使反应值得到最佳结果寻找因子参数设定使反应值得到最佳结果确认新的操作条件能使产品质量获得提升确认新的操作条件能使产品质量获得提升建构因子与反应值之间的关系式建构因子与反应值之间的关系式当不确定曲线关系是否存在时当不确定曲线关系是否存在时 当当DOE中发现有曲率（中发现有曲率（Factorial+Ct Point）系列化实验系列化实验-中央复合设计（中央复合设计（Central Composite Design，CCD）当事先已知有曲线当事先已知有曲线 3k全因子全因子 CCD Box-Benhnken设计设计RSM的使用的使用时时机机寻寻找因子参数找因子参数设设定使反定使反应值应值得到最佳得到最佳结结果果RSM二级模型的设计类型二级模型的设计类型1.3k全因子全因子2.中心组合（复合）设计中心组合（复合）设计（CCD）3.Box-Behnken设计（设计（BBD）RSM二二级级模型的模型的设计类设计类型型1.3k全因子全因子1.31.3k k全因子实验全因子实验全因子实验全因子实验K个因子，每个因子取三个水平个因子，每个因子取三个水平优点：能够估计所有主效果（线性的和二次的）和交互作用优点：能够估计所有主效果（线性的和二次的）和交互作用缺点：实验次数过多缺点：实验次数过多 1.3k全因子全因子实验实验K个因子，每个因子取三个水平个因子，每个因子取三个水平1.31.33 3全因子设计全因子设计全因子设计全因子设计1.33全因子全因子设计设计2.2.中心组合设计（中心组合设计（中心组合设计（中心组合设计（CCDCCD）中心复合设计是在中心复合设计是在2水平全因子和分部试验水平全因子和分部试验设计的基础上发展出来的一种试验设计方法，它是设计的基础上发展出来的一种试验设计方法，它是2水平全因子和分部试验设计的拓展。通过对水平全因子和分部试验设计的拓展。通过对2水平水平试验增加一个设计点（相当于增加了一个水平），试验增加一个设计点（相当于增加了一个水平），从而可以对评价指标（输出变量）和因素间的非线从而可以对评价指标（输出变量）和因素间的非线性关系进行评估。它常用于在需要对因素的非线性性关系进行评估。它常用于在需要对因素的非线性影响进行测试的试验。影响进行测试的试验。2.中心中心组组合合设计设计（CCD）中心复合中心复合设设中心复合设计的特点1、可可以以进进行行因因素素数数在在26个个范范围围内内的的试试验验。2、试试验验次次数数一一般般为为1490次次：2因因素素12次次，3因因素素20次次，4 因因素素30次次，5因因素素54次次，6因因素素90次次。3、可可以以评评估估因因素素的的非非线线性性影影响响。4、适适用用于于所所有有试试验验因因素素均均为为计计量量值值数数末末尾尾的的试试验验。5、在在使使用用时时，一一般般按按三三个个步步骤骤进进行行试试验验。（1）先先进进行行2水水平平全全因因子子或或分分部部试试验验设设计计。（2）再再加加上上中中心心点点进进行行非非线线性性测测试试。（3）如如果果发发现现非非线线性性影影响响为为显显著著影影响响，则则加加上上轴轴向向点点进进行行补补充充试试验验以以得得到到非非线线性性预预测测方方程程。6、中中心心复复合合试试验验也也可可一一次次进进行行完完毕毕，（在在确确信信有有非非线线性性影影响响的的情情况况下下）。中心复合中心复合设计设计的特点的特点1、可以、可以进进行因素数在行因素数在26个范个范围围内的内的试验试验。中心复合设计（CCD）优点：优点：1）能够预估所有主效果，双向交互作用和四分条件）能够预估所有主效果，双向交互作用和四分条件2）可以通过增加轴向点，从一级筛选设计转化而来（即中）可以通过增加轴向点，从一级筛选设计转化而来（即中心复合法）心复合法）缺点：缺点：1）轴向点的选择也许会造成在非理想条件下进行实验）轴向点的选择也许会造成在非理想条件下进行实验中心复合中心复合设计设计（CCD）优优点：点：立方点立方点 轴向点轴向点 中心点中心点 区组区组 序贯试验序贯试验 旋转性旋转性基本概念中心复合试验设计中心复合试验设计中心复合试验设计中心复合试验设计 立方点基本概念中心复合立方点基本概念中心复合试验设计试验设计中心复合试验中的立方点、轴向点和中心点中心复合试验设计由立方点、轴向点和中心点试验三部分组中心复合试验设计由立方点、轴向点和中心点试验三部分组成，下面以成，下面以2因子中心复合试验设计为例分别对三种点加以因子中心复合试验设计为例分别对三种点加以说明。说明。立方点立方点 立方点即全因子设计或分部试验设计中的立方点即全因子设计或分部试验设计中的2水平对应的水平对应的“-1”和和“+1”点，表示如下图：点，表示如下图：中心复合中心复合试验试验中的立方点、中的立方点、轴轴向点和中心点中心复合向点和中心点中心复合试验设计试验设计由立方由立方轴向点轴向点又称始点、星号点，分布在轴向上。除一个坐标为又称始点、星号点，分布在轴向上。除一个坐标为+或或-外，其余外，其余坐标皆为坐标皆为0。在。在k个因素的情况下，共有个因素的情况下，共有2k个轴向点。个轴向点。记为（记为（+a，0）、）、（-a，0）、（）、（0，+a）、（）、（0，-a），如下图表示。），如下图表示。轴向点轴向点a=n1/4,如：如：81/4=1.68，41/4=1.414轴轴向点向点轴轴向点向点a=n1/4,如：如：81/4=1.68，41/4中心点中心点中心点亦即设计中心，在坐标轴上表示为（中心点亦即设计中心，在坐标轴上表示为（0，0），表示在），表示在图上，坐标皆为图上，坐标皆为0。即（。即（0，0）点。将三种点集成在一个图上表）点。将三种点集成在一个图上表示如下：示如下：中心点中心点三因素下的立方点、轴向点和中心点三因素下的立方点、轴向点和中心点三因素下的立方点、三因素下的立方点、轴轴向点和中心点向点和中心点 序贯试验（顺序试验）先后分几段完成试验，前次试验设计的先后分几段完成试验，前次试验设计的点上做过的试验结果，在后续的试验设计中点上做过的试验结果，在后续的试验设计中继续有用。继续有用。序序贯试验贯试验（顺顺序序试验试验）先后分几段完成）先后分几段完成试验试验，前次，前次试验设计试验设计的点的点 旋转性(rotatable)设计旋转设计具有在设计中心旋转设计具有在设计中心等距点上预测等距点上预测方差恒定方差恒定的性质，这改善了预测精度。的性质，这改善了预测精度。旋旋转转性性(rotatable)设计设计旋旋转设计转设计具有在具有在设计设计中心等中心等的选取在在的选取上可以有多种出发点，旋转性是的选取上可以有多种出发点，旋转性是个很有意义的考虑。在个很有意义的考虑。在k k个因素的情况下，应个因素的情况下，应取取=2 =2 k/4k/4当当k=2k=2，=1.414=1.414；当；当k=3k=3，=1.682=1.682；当当k=4k=4，=2.000=2.000；当；当k=5k=5，=2.378=2.378的的选选取在取在的的选选取上可以有多种出取上可以有多种出发发点，旋点，旋转转性是个很有意性是个很有意义义的考的考 按上述公式选定的按上述公式选定的值来安排中心复合试验设计值来安排中心复合试验设计(CCD)(CCD)是最典型的情形，它可以实现试验的序贯性，这种是最典型的情形，它可以实现试验的序贯性，这种CCDCCD设计特称设计特称中心复合序贯设计中心复合序贯设计(central composite(central composite circumscribed design,CCC)circumscribed design,CCC)，它是，它是CCDCCD中最常用的一种。中最常用的一种。按上述公式按上述公式选选定的定的值值来安排中心复合来安排中心复合试验设计试验设计(CCD)对于对于值选取的另一个出发点也是有意义的，就是值选取的另一个出发点也是有意义的，就是取取=1=1，这意味着将轴向点设在立方体的表面上，这意味着将轴向点设在立方体的表面上，同时不改变原来立方体点的设置，这样的设计称为同时不改变原来立方体点的设置，这样的设计称为中心复合表面设计中心复合表面设计 (central composite face-(central composite face-centered design,CCF)centered design,CCF)。这样做，每个因素的取值水平只有这样做，每个因素的取值水平只有3 3个个(-1,0,1)(-1,0,1)，而，而一般的一般的CCDCCD设计，因素的水平是设计，因素的水平是5 5个个(-,-(-,-1,0,1,),1,0,1,),这在更换水平较困难的情况下是有意义这在更换水平较困难的情况下是有意义的。的。这种设计失去了旋转性。但保留了序贯性，即前一次在立方一次在立方点上已经做过的试验结果，在点上已经做过的试验结果，在后续的后续的CCFCCF设计中可以继续使用设计中可以继续使用，可以在二阶回归中采用。对对于于值选值选取的另一个出取的另一个出发发点也是有意点也是有意义义的，就是取的，就是取=1，这这意味意味中心点的个数选择在满足旋转性的前提下，如果适当选择在满足旋转性的前提下，如果适当选择NcNc，则可以使整个试验区域内的预测值都有一致均匀则可以使整个试验区域内的预测值都有一致均匀精度精度(uniform precision)(uniform precision)。见下表：。见下表：中心点的个数中心点的个数选择选择在在满满足旋足旋转转性的前提下，如果适当性的前提下，如果适当选择选择Nc，则则可可但有时认为，这样做的试验次数多，代价但有时认为，这样做的试验次数多，代价太大，太大，NcNc其实取其实取2 2以上也可以；如果中心点的以上也可以；如果中心点的选取主要是为了估计试验误差，选取主要是为了估计试验误差，NcNc取取4 4以上也以上也够了。够了。总之，当时间和资源条件都允许时，应尽总之，当时间和资源条件都允许时，应尽可能按推荐的可能按推荐的NcNc个数去安排试验，设计结果和个数去安排试验，设计结果和推测出的最佳点都比较可信。实在需要减少试推测出的最佳点都比较可信。实在需要减少试验次数时，中心点至少也要验次数时，中心点至少也要2-52-5次。次。但有但有时认为时认为，这样这样做的做的试验试验次数多，代价太大，次数多，代价太大，Nc其其实实取取2以上以上首先建立一个首先建立一个23因子设计因子设计统计统计DOE修改设计修改设计即一个单位的面上，即一个单位的面上，当轴向点太远时，实当轴向点太远时，实验条件达不到情况验条件达不到情况当轴向点太远时，实当轴向点太远时，实验条件达不到情况，验条件达不到情况，可以自己定义可以自己定义首先建立一个首先建立一个23因子因子设计设计即一个即一个单单位的面上，当位的面上，当轴轴向点太向点太远时远时，实实3.Box-Behnken试验设计(BBD)Box-Behnken试验设计是可以评价指标和试验设计是可以评价指标和因素间的非线性关系的一种试验设计方法。和中心因素间的非线性关系的一种试验设计方法。和中心复合设计不同的是它不需连续进行多次试验，并且复合设计不同的是它不需连续进行多次试验，并且在因素数相同的情况下，在因素数相同的情况下，Box-Behnken试验的试试验的试验组合数比中心复合设计少因而更经济。验组合数比中心复合设计少因而更经济。Box-Behnken试验设计常用于在需要对因素的非线性试验设计常用于在需要对因素的非线性影响进行研究时的试验。影响进行研究时的试验。3.Box-Behnken试验设计试验设计(BBD)Box-Behnken试验设计的特点试验设计的特点1、可以进行因素数在、可以进行因素数在37个范围内的试验。个范围内的试验。2、试验次数一般为、试验次数一般为15-62次。在因素数相同时比中心复合次。在因素数相同时比中心复合设计所需的试验次数少，比较如下设计所需的试验次数少，比较如下Box-Behnken试验设计试验设计的特点的特点1、可以、可以进进行因素数在行因素数在33、可以评估因素的非线性影响。、可以评估因素的非线性影响。4、适用于所有因素均为计量值的试验。、适用于所有因素均为计量值的试验。5、使用时无需多次连续试验。、使用时无需多次连续试验。6、Box-Behnken试验方案中没有将所有试验因素同时安试验方案中没有将所有试验因素同时安排为高水平的试验组合，对某些有特别需要或安全要求的排为高水平的试验组合，对某些有特别需要或安全要求的试验尤为适用。试验尤为适用。和中心复合试验相比，和中心复合试验相比，Box-Behnken试验设计试验设计不存在轴向点，因而在实际操作时其水平设置不会超出安不存在轴向点，因而在实际操作时其水平设置不会超出安全操作范围。而存在轴向点的中心复合试验却存在生成的全操作范围。而存在轴向点的中心复合试验却存在生成的轴向点可能超出安全操作区域或不在研究范围之列考虑的轴向点可能超出安全操作区域或不在研究范围之列考虑的问题。问题。3、可以、可以评评估因素的非估因素的非线线性影响。性影响。一个一个一个一个k=3 Box-Behnkenk=3 Box-Behnken的图像分析的图像分析的图像分析的图像分析注意：加入了一引进中心点，注意：加入了一引进中心点，并未增加轴向点，因而更完并未增加轴向点，因而更完全。设计并不包括全。设计并不包括任何极限任何极限值值，当因子在极限的组合因，当因子在极限的组合因为太昂贵，或根本无法进行为太昂贵，或根本无法进行实验时，这是一个有利的特实验时，这是一个有利的特性。性。当一个实验设计需要推倒从当一个实验设计需要推倒从来时，可以选择来时，可以选择BBD设计设计一个一个k=3 Box-Behnken的的图图像分析注意：加入了一引像分析注意：加入了一引统计统计DOE响应曲面响应曲面创建响应曲面设计创建响应曲面设计统计统计DOE响响应应曲面曲面创创建响建响应应曲面曲面设计设计中心复合法CCD中心复合法中心复合法CCD实验设计指南RSM1.1.问题的认知及陈述问题的认知及陈述2.2.反应变量的选择反应变量的选择3.3.因子选择与水平个数及范围的选择因子选择与水平个数及范围的选择4.4.选择合适的实验设计选择合适的实验设计5.5.进行试验收集数据进行试验收集数据实验设计实验设计指南指南RSM问题问题的的认认知及知及陈陈述述实验设计指南RSM6.资料分析资料分析为整个模型建立为整个模型建立Anova表表模式精简：去除不显著项（模式精简：去除不显著项（P-value高）或平方和影响低的高）或平方和影响低的项次（在项次（在Pareto图或常态图）后，进行模型的简化。切记：图或常态图）后，进行模型的简化。切记：一次删一项，重新分析再评估。一次删一项，重新分析再评估。注意注意Lack of fit问题是否显著问题是否显著解释能力是否足够：解释能力是否足够：R2值要大于值要大于80%。残差分析，确认模型的前提假设是否成立：四合一残差图残差分析，确认模型的前提假设是否成立：四合一残差图研究显著的交互作用研究显著的交互作用/主效应主效应(P-value小于小于0.05）-从高阶从高阶着手着手7.结论与建议结论与建议列出数学模型列出数学模型评估各方差源实际的重要性评估各方差源实际的重要性将模型转换为实际的流程设置（优化器）将模型转换为实际的流程设置（优化器）实验设计实验设计指南指南RSM6.资资料分析料分析例题一位化学工程师想了解使制程产能为最大的操作条件，有两个可控一位化学工程师想了解使制程产能为最大的操作条件，有两个可控因子会影响制程能力：反应时间和反应温度；因子会影响制程能力：反应时间和反应温度；工程师决定讨论制程在反应时间为（工程师决定讨论制程在反应时间为（80,90）分钟与反应温度在）分钟与反应温度在（170,180）F之范围的变化；之范围的变化；因为事前没有任何实验上的证据，而且因为时间上的急迫，所以工因为事前没有任何实验上的证据，而且因为时间上的急迫，所以工程师决定直接用一阶的实验来找到最佳化的条件，所以设计了一个两因程师决定直接用一阶的实验来找到最佳化的条件，所以设计了一个两因子两水平与一个区组化的响应曲面法；子两水平与一个区组化的响应曲面法；反应变量为产能（最低反应变量为产能（最低75，目标，目标80，望大），产品粘度，望大），产品粘度（60,65,70），分子量），分子量Molecular Weight（3000,3200,3400）RSM-CCD1.mtx例例题题一位化学工程一位化学工程师师想了解使制程想了解使制程产产能能为为最大的操作条件，有两个可最大的操作条件，有两个可试验数据5 1-1177.9289 175.00075.61 2 0185.0000 175.00080.06 3-1192.0711 175.00078.42 41190.0000 170.00078.08 5-1185.0000 182.07178.51 61180.0000 170.00076.512 70185.0000 175.00079.713 80185.0000 175.00079.87 9-1185.0000 167.92977.010 100185.0000 175.00080.34 111190.0000 180.00079.53 121180.0000 180.00077.09 130185.0000 175.00079.2StdOrder RunOrder PtType Blocks Time Temp Productivity试验试验数据数据StdOrder RunOrder PtType 完整模型之ANOVAProductivity 的估计回归系数的估计回归系数 项项 系数系数 系数标准误系数标准误 T P常量常量 79.8000 0.1642 486.087 0.000Time 0.9950 0.1298 7.666 0.000Temp 0.5152 0.1298 3.969 0.005Time*Time -1.3062 0.1392 -9.385 0.000Temp*Temp -0.9312 0.1392 -6.691 0.000Time*Temp 0.2500 0.1835 1.362 0.215S=0.367091 PRESS=3.04577R-Sq=96.53%R-Sq（预测）（预测）=88.80%R-Sq（调整）（调整）=94.06%可以简化哪项？解释能力是否足够？完整模型之完整模型之ANOVAProductivity 的估的估计计回回归归系数系数Productivity 的方差分析的方差分析来源来源 自由度自由度 Seq SS Adj SS Adj MS F P回归回归 5 26.2598 26.2598 5.2520 38.97 0.000 线性线性 2 10.0430 10.0430 5.0215 37.26 0.000 Time 1 7.9198 7.9198 7.9198 58.77 0.000 Temp 1 2.1232 2.1232 2.1232 15.76 0.005 平方平方 2 15.9668 15.9668 7.9834 59.24 0.000 Time*Time 1 9.9339 11.8698 11.8698 88.08 0.000 Temp*Temp 1 6.0329 6.0329 6.0329 44.77 0.000 交互作用交互作用 1 0.2500 0.2500 0.2500 1.86 0.215 Time*Temp 1 0.2500 0.2500 0.2500 1.86 0.215残差误差残差误差 7 0.9433 0.9433 0.1348 失拟失拟 3 0.2833 0.2833 0.0944 0.57 0.663 纯误差纯误差 4 0.6600 0.6600 0.1650合计合计 12 27.2031是不是缺失度的问题是不是缺失度的是不是缺失度的问题问题产能最优化-结论与建议产产能最能最优优化化-结论结论与建与建议议重要度与权重重要度（重要度（0.1-10），又称相对重要度，优化器优先满足重要度高的），又称相对重要度，优化器优先满足重要度高的响应输出。响应输出。权重（权重（0.1-10），又称满足目标的程度。权重越高，必须要在），又称满足目标的程度。权重越高，必须要在target越接近才能满意。越接近才能满意。结论：结论：两个因子对产能的影响都是重要的，而且其效应是有二次式存在两个因子对产能的影响都是重要的，而且其效应是有二次式存在最佳设计条件：目标为望大最佳设计条件：目标为望大时间：时间：86.92min温度：温度：176.35度度产能：预计可以达到产能：预计可以达到80.06重要度与重要度与权权重重要度（重重要度（0.1-10），又称相），又称相对对重要度，重要度，优优化器化器优优练习接续前一个情况：接续前一个情况：针对黏度（针对黏度（Viscosity）（望目：）（望目：60,65,70）与分子量（）与分子量（Molecular）（望目：）（望目：3000,3200，3400）分布讨论其数学式与结论）分布讨论其数学式与结论练习练习接接续续前一个情况：前一个情况：StdOrder RunOrder PtType Blocks Time Temp Productivity Viscosity Molecular试验数据StdOrder RunOrder PtType BlockRSM-CCD10.mtx例题-重叠等值线图RSM-CCD10.mtx例例题题-重叠等重叠等值线图值线图RSM-CCD1.mtx响应优化器RSM-CCD1.mtx响响应优应优化器化器例题-step1一位制程工程师想要了解一个化学制程的最佳化设定，其目的是要产能一位制程工程师想要了解一个化学制程的最佳化设定，其目的是要产能为最大化。所以该工程师要先进行一个两因子（时间和压力）加中心点的实验为最大化。所以该工程师要先进行一个两因子（时间和压力）加中心点的实验设计，以确定是否有曲率现象存在。设计，以确定是否有曲率现象存在。因子因子:时间：时间：80min，100min 温度：温度：140度，度，150度度反应变数：产率（反应变数：产率（80，95，望大），望大）RSM-CCD2-step1例例题题-step1一位制程工程一位制程工程师师想要了解一个化学制程的最佳化想要了解一个化学制程的最佳化设设例题-step2该工程师发现线性效应不成立，决定进行第二阶段实验，以取得最该工程师发现线性效应不成立，决定进行第二阶段实验，以取得最佳化的模型与参数设定。佳化的模型与参数设定。第二阶段以轴点第二阶段以轴点+中心点，实验数据已经在中心点，实验数据已经在RSM-CCD-step2中中将第二阶段的数据复制到第一阶段的后半部，形成一个具有两因子将第二阶段的数据复制到第一阶段的后半部，形成一个具有两因子两个区组化的两个区组化的CCD设计，设计，请进行分析，获得最佳输出结果。请进行分析，获得最佳输出结果。RSM-CCD2-step2.mtwRSM-CCD2.mtw例例题题-step2该该工程工程师发现线师发现线性效性效应应不成立，决定不成立，决定进进行第二行第二阶阶段段加入新增轴向点+中心点加入新增加入新增轴轴向点向点+中心点中心点输出结果结果结果:RSM_CCD2.mtw 响应曲面回归响应曲面回归:Yield 与区组与区组 ,Time,Temp 分析是使用已编码单位进行的。分析是使用已编码单位进行的。Yield 的估计回归系数的估计回归系数项项 系数系数 系数标准误系数标准误 T P常量常量 87.7667 0.7179 122.251 0.000区组区组 0.6214 0.4700 1.322 0.228Time -1.3837 0.6217 -2.226 0.061Temp 0.3620 0.6217 0.582 0.579Time*Time -2.3396 0.6471 -3.615 0.009Temp*Temp -3.2896 0.6471 -5.083 0.001Time*Temp -4.8750 0.8793 -5.544 0.001S=1.75854 PRESS=138.459R-Sq=91.37%R-Sq（预测）（预测）=44.78%R-Sq（调整）（调整）=83.97%区组非显区组非显著，可以著，可以进一步简进一步简化化 输输出出结结果果结结果果:RSM_CCD2.mtw区区组组非非显显著，可以著，可以进进一一等值线图与曲面图等等值线图值线图与曲面与曲面图图等值线图等等值线图值线图结论与建议Yield 的估计回归系数，使用未编码单位的数据的估计回归系数，使用未编码单位的数据项项 系数系数常量常量 -4138.70Time 18.2104Temp 47.0066Time*Time -0.0233958Temp*Temp -0.131583Time*Temp -0.0975000结论结论与建与建议议Yield 的估的估计计回回归归系数，使用未系数，使用未编码单编码单位的数据位的数据最优化工具最最优优化工具化工具结论试验的目的是要确定可以让得率为最大的因子设定。由此序列化试验的目的是要确定可以让得率为最大的因子设定。由此序列化实验发现：实验发现：在时间在时间=76min和温度和温度=150度，产能可以达到度，产能可以达到89以上以上由等值线图我们可以发现将时间下降与温度提高时，似乎有较高由等值线图我们可以发现将时间下降与温度提高时，似乎有较高的得率，然而，如果还有其它质量特性（例如成本）必须同时考虑时，的得率，然而，如果还有其它质量特性（例如成本）必须同时考虑时，就不得如此贸然行事；因为虽然可以蔣得率提升，但相对是否会让其就不得如此贸然行事；因为虽然可以蔣得率提升，但相对是否会让其它质量特性恶化。它质量特性恶化。结论试验结论试验的目的是要确定可以的目的是要确定可以让让得率得率为为最大的因子最大的因子设设定。由此序列化定。由此序列化例题一位黑带应邀帮助把一个试产制程推广到更大范围，该一位黑带应邀帮助把一个试产制程推广到更大范围，该制程的目的是把一间金属加工厂之废水中金属污染消除掉。制程的目的是把一间金属加工厂之废水中金属污染消除掉。处理废水的下一步要用一种昂贵的有机树脂，所以他们要在处理废水的下一步要用一种昂贵的有机树脂，所以他们要在离子交换程序前尽可能将废水中金属污染去除。离子交换程序前尽可能将废水中金属污染去除。分析两个因子：时间和温度分析两个因子：时间和温度输出：除去金属的重量；输出：除去金属的重量；RSM-Exercise1.mtw例例题题一位黑一位黑带应带应邀帮助把一个邀帮助把一个试产试产制程推广到更大范制程推广到更大范围围，该该制程的目制程的目输出结果-注意需要选中区组选项Weight 的估计回归系数的估计回归系数 系数系数 标标项项 系数系数 准误准误 T P常量常量 760.721 6.207 122.564 0.000区组区组 10.716 4.516 2.373 0.049Time -8.405 7.259 -1.158 0.285Temp 20.804 5.375 3.870 0.006Time*Time -38.488 7.620 -5.051 0.001Temp*Temp -14.440 5.854 -2.467 0.043S=15.2033 PRESS=5296.43R-Sq=88.44%R-Sq（预测）（预测）=62.15%R-Sq（调整）（调整）=80.18%输输出出结结果果-注意需要注意需要选选中区中区组选项组选项例题-结果例例题题-结结果果练习提高烧碱纯度问题。在烧碱生产中，经过因子的筛选，提高烧碱纯度问题。在烧碱生产中，经过因子的筛选，最后得知反应炉内压力及温度是两个关键因子。再改进最后得知反应炉内压力及温度是两个关键因子。再改进阶段先进行全因子实验：阶段先进行全因子实验：A 压力：压力：-50Bar，60BarB 温度：温度：-260度，度，320度度中心点也做了三次实验，实验结果见：中心点也做了三次实验，实验结果见：DOE-烧碱纯度烧碱纯度（反应（反应1）.mtw烧碱纯度（烧碱纯度（Larger is better):LSL=80,Target=100.练习练习提高提高烧烧碱碱纯纯度度问题问题。在。在烧烧碱生碱生产产中，中，经过经过因子的因子的筛选筛选，最后得知，最后得知练习（续）经第一步的分析得知，的确存在有曲率的现象。因此规划经第一步的分析得知，的确存在有曲率的现象。因此规划中加入四个轴点位置的实验，构成一个响应曲面设计。中加入四个轴点位置的实验，构成一个响应曲面设计。全部资料见：全部资料见：DOE-烧碱纯度（反应烧碱纯度（反应2）.mtw.请确认请确认:最佳华模式最佳华模式最佳参数设定。最佳参数设定。练习练习（续续）经经第一步的分析得知，的确存在有曲率的第一步的分析得知，的确存在有曲率的现现象。因此象。因此规规划划最陡上升路线法-寻找最佳区间寻找最佳区间最陡上升路最陡上升路线线法法-寻寻找最佳区找最佳区间间Opportunity f如正在爬山而看不见山顶。如正在爬山而看不见山顶。如正在爬山而看不如正在爬山而看不见见山山顶顶。最陡上升路线法最陡上升路线法 -寻找最佳区间寻找最佳区间寻找最佳区间寻找最佳区间1）第一步：确定有可能含有最佳条件的区）第一步：确定有可能含有最佳条件的区域域2）第二步：建立一个描述响应与重要变量）第二步：建立一个描述响应与重要变量之间关系的模型之间关系的模型3）第三步：利用该模型进行过程优化，并）第三步：利用该模型进行过程优化，并且确认结果且确认结果最陡上升路最陡上升路线线法法-寻寻找最佳区找最佳区间间1）第一步：确定有可能含有）第一步：确定有可能含有最陡上升路线最陡上升路线假设已经做了一个实验，并且得到以下模型：假设已经做了一个实验，并且得到以下模型：Y=0 0+1 1x x1 1+2 2x x2 2+（同主效果有关）（同主效果有关）+1212x x1212+（同双向交互作用有关）（同双向交互作用有关）我们可以利用这个模型确定一条最陡上升路线，帮助我们我们可以利用这个模型确定一条最陡上升路线，帮助我们更接近最佳条件。更接近最佳条件。最陡上升路最陡上升路线线假假设设已已经经做了一个做了一个实验实验，并且得到以下模型：，并且得到以下模型：坡度坡度根据微积分学，以下函数的坡度根据微积分学，以下函数的坡度Y=f（x1,x2,xn)微分：微分：坡度根据微坡度根据微积积分学，以下函数的坡度分学，以下函数的坡度寻找最佳上升路线举例寻找最佳上升路线举例6sigma小组打算优化某一制程的良率。分析阶段过后，识别出两个小组打算优化某一制程的良率。分析阶段过后，识别出两个因子（因子（KPIV）因子因子A:温度（温度（F）因子因子B：压力（：压力（psi）文件在文件在RSM01.mtw寻寻找最佳上升路找最佳上升路线举线举例例6sigma小小组组打算打算优优化某一制程的良率。化某一制程的良率。因因子子图图及及交交互互作作用用因子因子图图及交互作用及交互作用举例举例举例举例-分析结果分析结果分析结果分析结果Yield 的效应和系数的估计（已编码单位）的效应和系数的估计（已编码单位）项项 效应效应 系数系数 系数标准误系数标准误 T P 常量常量 74.450 0.2121 350.96 0.002Temp 9.100 4.550 0.2121 21.45 0.030Pres -2.900 -1.450 0.2121 -6.84 0.092Temp*Pres -3.300 -1.650 0.2121 -7.78 0.081Ct Pt 1.450 0.3674 3.95 0.158来源来源 自由度自由度 Seq SS Adj SS Adj MS F P主效应主效应 2 91.220 91.2200 45.6100 253.39 0.0442因子交互作用因子交互作用 1 10.890 10.8900 10.8900 60.50 0.081 弯曲弯曲 1 2.803 2.8033 2.8033 15.57 0.158残差误差残差误差 1 0.180 0.1800 0.1800 纯误差纯误差 1 0.180 0.1800 0.1800合计合计 5 105.093曲率效果不明显曲率效果不明显关系式：关系式：Yield=74.45+4.55*Temp-1.45*Press-1.65*Temp*PressH0和和H1分别是？分别是？举举例例-分析分析结结果果Yield 的效的效应应和系数的估和系数的估计计（已（已编码单编码单位）曲位）曲3D曲面图曲面图图形图形3D曲面图曲面图3D曲面曲面图图图图形形3D曲面曲面图图图形图形等值线图等值线图图图形等形等值线图值线图选选出出显显著著效效果果通过方差分析（通过方差分析（ANOVA）表中，选出显著因子（）表中，选出显著因子（a风险风险是是0.1）A：温度：温度B：压力：压力AB：AB交互作用交互作用模型：模型：Yield=74.45+4.55*Temp-1.45*Press-1.65*Temp*Press选选出出显显著效果通著效果通过过方差分析（方差分析（ANOVA）表中，）表中，选选出出显显著因子（著因子（a最陡的上升路线最陡的上升路线最陡的上升路线的方向由以下得出：最陡的上升路线的方向由以下得出：x*i的起点被确定的起点被确定将设计中心（将设计中心（0，0）设为起点，得到：）设为起点，得到：最陡的上升路最陡的上升路线线最陡的上升路最陡的上升路线线的方向由以下得出：的方向由以下得出：沿着路线前进沿着路线前进沿着路线前进沿着路线前进 X*T X*P 0 03.14 -16.28 -29.42 -3沿着路沿着路线线前前进进 X*T X*P将已解码变为未解码将已解码变为未解码 XT50X*T+150,XP25X*P+150将已解将已解码变为码变为未解未解码码新的试验结果新的试验结果温度温度压力压力 150 307 464 621 150 125 100 75新的新的试验结试验结果温度果温度压压力力 150 307 RSM分析步骤分析步骤步骤：步骤：1）求）求p，确定在山坡位置，确定在山坡位置2）写出关系式）写出关系式3）偏导）偏导4）出发点，确定方向（一般以中心点出发）出发点，确定方向（一般以中心点出发）5）确定步幅（确定做实验的规格点），写出实验方案（代）确定步幅（确定做实验的规格点），写出实验方案（代码化的）码化的）6）代码化转化成实际方案）代码化转化成实际方案RSM分析步分析步骤骤步步骤骤：练习（练习（练习（练习（2 2水平水平水平水平+5+5中心点）中心点）中心点）中心点）40.340.340.540.540.740.740.240.240.640.616015515030 35 40练习练习（2水平水平+5中心点）中心点）40.316030 练练习习-结结果果Yield=40.425+0.775*Temp-0.325*Press-0.025*Temp*Press练习练习-结结果果Yield=40.425+0.775*Temp-0练习练习-结果结果将已解码变为未解码将已解码变为未解码 XT5X*T+35,XP5X*P+155练习练习-结结果将已解果将已解码变为码变为未解未解码码练练练练习习习习（2 2水水水水平平平平+4 4中中中中心心心心点点点点）40.340.540.740.240.6605040100 150 200100.00102.0098.0099.00练习练习（2水平水平+4中心点）中心点）40.360100 选出新的试验水平（见选出新的试验水平（见P17）见数据：见数据：RSM02.mtw选选出新的出新的试验试验水平（水平（见见P17）见见数据：数据：RSM02.mtw当发现一级模型无效时，该怎么做？当发现一级模型无效时，该怎么做？当发现一级模型无效时，该怎么做？当发现一级模型无效时，该怎么做？拟合因子拟合因子:Yield 与与 Temp,Press Yield 的效应和系数的估计（已编码单位）的效应和系数的估计（已编码单位）项项 效应效应 系数系数 系数标准误系数标准误 T P常量常量 83.175 0.3182 261.39 0.002Temp 3.750 1.875 0.3182 5.89 0.107Press -3.650 -1.825 0.3182 -5.74 0.110Temp*Press 7.550 3.775 0.3182 11.86 0.054Ct Pt 7.375 0.5511 13.38 0.047S=0.636396 PRESS=*R-Sq=99.74%R-Sq（预测）（预测）=*%R-Sq（调整）（调整）=98.71%对于对于 Yield 方差分析（已编码单位）方差分析（已编码单位）来源来源 自由度自由度 Seq SS Adj SS Adj MS F P主效应主效应 2 27.385 27.3850 13.6925 33.81 0.1212因子交互作用因子交互作用 1 57.003 57.0025 57.0025 140.75 0.054 弯曲弯曲 1 72.521 72.5208 72.5208 179.06 0.047残差误差残差误差 1 0.405 0.4050 0.4050 纯误差纯误差 1 0.405 0.4050 0.4050合计合计 5 157.313中心点显著，曲面中心点显著，曲面效益显著。效益显著。当当发现发现一一级级模型无效模型无效时时，该该怎么做？怎么做？拟拟合因子合因子:Yield 与与 主主效效果果分分析析图图主效果分析主效果分析图图需需要要一一个个更更高高级级模模型型曲面曲面P-数值的统计显著性表明一级模型已经不够用了数值的统计显著性表明一级模型已经不够用了 Y=0 0+1 1x x1 1+2 2x x2 2+（同主效果有关）（同主效果有关）+1212x x1212+（同双向交互作用有关）（同双向交互作用有关）下面我们该怎么做呢？快找出最佳条件的时候，往往需要下面我们该怎么做呢？快找出最佳条件的时候，往往需要二次模型二次模型Y=0 0+1 1x x1 1+2 2x x2 2+（同主效果有关）（同主效果有关）+1212x x1212+（同双向交互作用有关）（同双向交互作用有关）+1111x x1 12 2+(同纯二次项有关）同纯二次项有关）需要一个更高需要一个更高级级模型曲面模型曲面P-数数值值的的统计显统计显著性表明一著性表明一级级模型已模型已经经不不RSMRSM设计设计设计设计最佳区域，最佳区域，最高良率最高良率当曲面（或中心点）变为显著时，我们已经达到顶部，需要当曲面（或中心点）变为显著时，我们已经达到顶部，需要RSM方法进一步方法进一步勘查最佳区域。勘查最佳区域。RSM设计设计最佳区域，最高良率当曲面（或中心点）最佳区域，最高良率当曲面（或中心点）变为显变为显著著时时，我，我二阶模型之二阶模型之CCD法法回顾之间的例题回顾之间的例题见数据：见数据：RSM03.mtw二二阶阶模型之模型之CCD法法回回顾顾之之间间的例的例题见题见数据：数据：RSM03.mtw统计统计DOE响应曲面响应曲面创建响应曲面设计创建响应曲面设计统计统计DOE响响应应曲面曲面创创建响建响应应曲面曲面设计设计MinitabMinitab输出输出输出输出分析是使用已编码单位进行的。分析是使用已编码单位进行的。Yield 的估计回归系数的估计回归系数 系数标系数标项项 系数系数 准误准误 T P常量常量 90.5500 2.134 42.431 0.000Temp 2.3164 1.067 2.171 0.096Press -0.3645 1.067 -0.342 0.750Temp*Temp -3.7875 1.412 -2.683 0.055Press*Press -5.0625 1.412 -3.587 0.023Temp*Press 3.7750 1.509 2.502 0.067S=3.01799 PRESS=257.819R-Sq=86.47%R-Sq（预测）（预测）=4.24%R-Sq（调整）（调整）=69.55%对于对于 Yield 的方差分析的方差分析来源来源 自由度自由度 Seq SS Adj SS Adj MS F P回归回归 5 232.808 232.808 46.5616 5.11 0.070 线性线性 2 43.987 43.987 21.9935 2.41 0.205 平方平方 2 131.818 131.818 65.9092 7.24 0.047 交互作用交互作用 1 57.003 57.003 57.0025 6.26 0.067残差误差残差误差 4 36.433 36.433 9.1083 失拟失拟 3 36.028 36.028 12.0093 29.65 0.134 纯误差纯误差 1 0.405 0.405 0.4050合计合计 9 269.241Minitab输输出分析是使用已出分析是使用已编码单编码单位位进进行的。行的。等高等高/表面图表面图统计统计DOE响应曲面响应曲面等值线等值线/表面图表面图情况一等高等高/表面表面图统计图统计DOE响响应应曲面等曲面等值线值线/表面表面图图情况一情况一选择初步最佳条件选择初步最佳条件选择初步最佳条件选择初步最佳条件找到具体最佳条件需要通过解开一系列公式，当因子数量找到具体最佳条件需要通过解开一系列公式，当因子数量超过两个的时候，这将是一个非常复杂的工作超过两个的时候，这将是一个非常复杂的工作使用使用“多反应优化法多反应优化法”（Multiple Responses Optimism）找到大约最佳条件通常已足够）找到大约最佳条件通常已足够选择选择初步最佳条件找到具体最佳条件需要通初步最佳条件找到具体最佳条件需要通过过解开一系列公式，当因解开一系列公式，当因选择初步最佳条件选择初步最佳条件选择初步最佳条件选择初步最佳条件选择选择初步最佳条件初步最佳条件CCD练习练习Data：CDD-1.mtw课题：粘合剂生产条件优化课题：粘合剂生产条件优化粘合剂生产条件优化问题。在粘合剂生产中，经过因子的筛选，最后得知，反应罐粘合剂生产条件优化问题。在粘合剂生产中，经过因子的筛选，最后得知，反应罐内温度及反应时间是两个关键因子。在本阶段的最初全内温度及反应时间是两个关键因子。在本阶段的最初全 因子实验时因子实验时，因子，因子 A（Temp）的）的 低水平及高水平取为低水平及高水平取为200度及度及300度，因子度，因子B（Time）的低水平及高）的低水平及高水平分别取为水平分别取为40s和和70s，在中心处也做了三次实验，实验结果如表所示：，在中心处也做了三次实验，实验结果如表所示：CCD练习练习Data：CDD-1.mtw课题课题：粘合：粘合剂剂生生产产条件条件优优增加试验数据增加试验数据数据：CDD-2.mtw增加增加试验试验数据数据：数据数据：CDD-2.mtw响应优化器响应优化器目标目标目标目标100100，最小为，最小为，最小为，最小为1010分析：分析：当温度（temp）取259.5281，时间（time)取67.6777）时，所获得的黏度最大，最佳值可以达到50.6927.计算机提供了自动求最优解的功能，利用“响应变量优化器”可以直接获得最佳点的设置及最佳值，同时可以用人工进行调整，对最优点取整等等。响响应优应优化器化器目目标标100，最小，最小为为10分析：分析：试验策划试验策划建立试验目标建立试验目标本试验的目标是确定化工原料化学反应工序的温度、压力和本试验的目标是确定化工原料化学反应工序的温度、压力和反应时间与原料关键参数反应时间与原料关键参数Y之间的关系，并对三种因素的水之间的关系，并对三种因素的水平设置进行优化。平设置进行优化。确定测量指标确定测量指标小组确定以该原料的关键参数小组确定以该原料的关键参数Y为测量指标，其目标值为为测量指标，其目标值为873mg确定影响因素确定影响因素XS小组确定的影响因素为小组确定的影响因素为1、反应压力、反应压力2、反应温度、反应温度3、反应时间、反应时间Box-Behnken试验设计（试验设计（BBD）例）例试验试验策划策划Box-Behnken试验设计