响应曲面设计培训资料课件.ppt

1、响应面试验设计Response SurfaceMethodology,RSMYouthRSM LIULI关于RSMRSMRSMRSM是利用合理的实验设计方法合理的实验设计方法并通过实验的到一定的数据，采用多元二次回归多元二次回归方程方程来拟合因素与响应值之间的函数关系，通过对回归方程的分析来寻求最优最优的工艺参数的工艺参数，解决多变量问题的一种统计方法。RSM可用三维效应面三维效应面或二二维等高线图维等高线图表示。由此可以直观的看出自变量去不同值时的效应值。反过来也可通过效应面上选取一定的效应值也可以找出对应的自变量取值，即在选定的最佳值范围内可以找到最佳的试验条件最佳的试验条件。什么是什么是

2、RSMRSM？目录目录响应曲面设计概论响应曲面设计概论BoxBox-Behnken-Behnken试验设计试验设计中心复合设计中心复合设计CCDCCD实例和总结实例和总结响应曲面设计概论响应曲面设计概论一、响应曲面设计概论1 1、什么是响应曲面设计？、什么是响应曲面设计？通过对响应的曲面图形进行分析，寻找最佳响应的设计方法。2 2、包含二次项的回归方程、包含二次项的回归方程一般的形式如下：Y=bY=bo o+b+b1 1x x1 1+b+b2 2x x2 2+b+b1111x x1 12 2+b+b2222x x2 22 2+b+b1212x x1 1x x2 2+由于增加了两个因子各自的平方

3、项，需要增加试验点。先后分为几个阶段完成全部试验的策略，称为序贯试验策略一、响应曲面设计概论3 3、怎样获得响应的曲面图形？、怎样获得响应的曲面图形？大概的步骤如下：先用包含中心点的2水平因子试验的数据，拟合一个线性回归方程（可以包含交叉乘积项）如果发现曲面有弯曲的趋势，则应拟合一个含有二次项的回归方程；如果无，且Y没有达到目标，则用最速下降法最速下降法寻找最优的区域，直到达成目标或发现弯曲再拟合含有二次项的回归方程。一、响应曲面设计概论响应曲面设计适用范围响应曲面设计适用范围:确信或怀疑因素对指标存在非线性非线性影响因素个数2-72-7个，一般在4个以内所有因素均为计量值数据试验区域已接近最

4、优最优区域基于2水平的全因子正交试验一、响应曲面设计概论响应曲面设计的方法分为两类：中心复合序贯设计（CCC）中心复合试验设计（CCD）中心复合有界设计（CCI）中心复合表面设计（CCF）Box-Behnken试验设计（BBD）中心复合试验设计中心复合试验设计central composite desingncentral composite desingn中心复合试验设计CCDCCDCCD由以下由以下3 3类点构成的实验设计称为类点构成的实验设计称为中心复合设计中心复合设计。立方点 CCD由3部分组成：轴向点 中心点 CCD的决定因素：序贯性设计 旋转性设计 中心复合试验设计中心复合试验设计

5、 CCDCCD立方点立方点(cube point)(cube point)又称为立方体点、角点，即时2水平对应的“-1-1”和“+1+1”点。各点坐标皆为+1或-1.在k个因素的情况下，共有2 2k k个立方点。轴向点轴向点（axial pointaxial point）又称为始点、星号点，分布在轴向上。除一个坐标为“+”或-“”外，其余坐标皆为0,。在K个因素的情况下，共有2k2k个轴向点。中心复合试验设计中心复合试验设计 CCDCCD中心点中心点（center pointcenter point）中心点，即时设计中心，表示在图上，坐标皆为“0 0”.序贯试验设计（顺序试验）序贯试验设计（顺

6、序试验）线后分段完成试验，前次试验设计的点上做过的试验结果，在后续的试验设计中继续有用。中心复合试验设计中心复合试验设计 CCDCCD旋转性旋转性(rotatable)(rotatable)设计设计 将来在某点处预报值的方差仅与该店到试验中心的距离有关，即响应变量的预测精度在意设计中心为球心的球面上是相同的，可保证均匀一致的精度。旋转性旋转性具有在设计中心等距点上预测设计中心等距点上预测方差恒定方差恒定的性质。中心复合试验设计中心复合试验设计 CCDCCD中心复合设计一般步骤：中心复合设计一般步骤：确定因素及水平，注意水平数为水平数为2 2，因素数一般不超过4个，因素均为计量数据.创建“中心复

7、合中心复合”或“Box-Behnken Box-Behnken 试验设计试验设计”确定试验运行顺序（Display Design）进行试验并收集数据分析数据优化因素的设置水平中心复合设计实验方案的确定中心复合设计实验方案的确定第一步：第一步：如果只选择立方点和中心点立方点和中心点，则构成一般3 3因子因子2 2水平水平的全因子设计，可以拟合各因子的主效应和二阶、三阶主效应和二阶、三阶交互效应交互效应。第二步第二步:如果拟合模型三维图形三维图形出现弯弯曲曲的情况，可在上述全因子设计的基础上增加增加6 6个轴点个轴点，完成第二阶段的试验，可拟合各因子的二阶项，即为序贯试序贯试验验。中心复合设计实验

8、方案的确定中心复合设计实验方案的确定第三步第三步:如果确定试验区域已经接近最优区域接近最优区域，则可选择三类点直接进行中心复合设计。需要考虑的问题如下：1 1、如何选择全因子设计部分、如何选择全因子设计部分 2 2、如何确定星号点的位置（即确定、如何确定星号点的位置（即确定值）值）3 3、如何确定中心点的个数、如何确定中心点的个数中心复合设计实验方案的确定中心复合设计实验方案的确定1 1、如何选择全因子设计部分、如何选择全因子设计部分一般选择全因子设计全因子设计（因子数在2-42-4之间），因子数 5 5时考虑采用部分因子设计部分因子设计。2 2、如何确定星号点的位置（多考虑旋转性）、如何确定

9、星号点的位置（多考虑旋转性）F=KF=K2 2 或 F=1/2KF=1/2K2 2(5个因素）=(F)=(F)1/41/4F为因子试验点的总数，K为因子的个数即可满足旋转性旋转性，又可满足序贯性序贯性，称为中心复合贯序设计中心复合贯序设计（CCC）。=2K/4中心复合设计实验方案的确定中心复合设计实验方案的确定3 3、中心点个数的选择、中心点个数的选择在满足旋转性的前提下，适当的选择中心点数，可以使整个试验区域内的预测值具有一致均匀精度。一般至少选2-52-5次。因子数立方点星号点中心点合计244513386620416863053210105251610733Box-BehnkenBox-B

10、ehnken试验设计试验设计 Box-Behnken desingn Box-Behnken desingnBox-BehnkenBox-Behnken试验设计（试验设计（BBDBBD）将各试验点取在立方棱的中点上各试验点取在立方棱的中点上，所需要的点数比CCD少，试验区域是球形的，有近似旋转性，无序贯性，3 3个因子个因子需要12+312+3次试验，4 4个因子个因子需要24+324+3次试验。Box-BehnkenBox-Behnken试验设计特点试验设计特点1、在因素相同时，比中心复合设计的试验次数少试验次数少2、没有将所有试验因素同时安排为高水平的试验组合，对某些有安全要求或特别需求的

11、试验尤为适用3、具有近似旋转性，无序贯性。有近似旋转性，无序贯性。3 3因子4 4种响应曲面试验点计划表响应面设计实例和总结响应面设计实例和总结响应面设计实例（响应面设计实例（CCDCCD）因素因素水平水平-1.732-10+1+1.732X1乙醇浓度50.0058.4570.0081.5590.00X2提取时间4073.81120.00166.19200.00X3溶剂量6.007.6910.0012.3114.00因素水平表因素水平表模型拟合 以远志皂苷元为因变量，对各因素进行多元线性回归和二项式拟合，模型如下：多元线性回归：Y=b0+b1x1+b2x2+b3x3二项式：Y=b0+b1x1+

12、b2x2+b3x3+b4x12+b5x22+b6x32+b7x1x2+b8x1x3+b9x2x3星点试验设计与结果试验号试验号X1X1X2X2X3X3远志皂苷含量远志皂苷含量远志皂苷元提取率远志皂苷元提取率1-1-1-123.3863.502+1-1-119.3945.443-1+1-126.3778.094+1+1-122.5257.445-1-1+126.5172.116+1-1+119.0145.817-1+1+130.9983.668+1-1+126.7673.249-1.7320028.7982.4310+1.7320016.9332.33110-1.732020.9048.3712

13、0+1.732028.9481.541300-1.73225.6869.021400+1.73229.9081.7615-2000028.5478.04 多元线性回归：r=0.8526 F(3,16)=14.194 P0.01 虽然通过试验，但拟合度不佳，预测性较差，因此线性模型不合适。数学模型：Y=160.4345160.4345+6.28596.2859x1+6.28596.2859x20.05210.0521x12-0.0021x22+0.0136x2x3复相关系数r=0.9790，相对线性拟合有大幅度的提高，方程删项简化后，r值降幅很小，表明该方程有较大的可信度。b0b1b2b4b5b9估计值-160.43456.28596.2859-0.0521-0.00210.0136标准误33.98940.94570.10830.00670.00040.0033t-4.72016.64665.0910-7.73272-4.92254.0886P0.00030.00000.00020.00000.00020.0011工艺优化与预测提取浓度：提取浓度：150-200min150-200min提取时间：提取时间：55%-65%55%-65%溶剂量：溶剂量：12-1412-14倍倍响应面设计总结响应面设计总结谢谢大家谢谢大家