响应面分析实用举例.pptx

会计学1响应面分析实用举例响应面分析实用举例一、为什么要做响应面分析一、为什么要做响应面分析一、为什么要做响应面分析一、为什么要做响应面分析n n许多试验设计与优化方法，特别是在做回归分析过程中，都未能给出直观的许多试验设计与优化方法，特别是在做回归分析过程中，都未能给出直观的图形，因而也不能凭直觉观察其最优化点，虽然能找出最优值，但难以直观图形，因而也不能凭直觉观察其最优化点，虽然能找出最优值，但难以直观地判别优化区域。地判别优化区域。n n响应面分析是将体系的响应（如萃取化学中的萃取率）作为一个或多个因素响应面分析是将体系的响应（如萃取化学中的萃取率）作为一个或多个因素（如萃取剂浓度、酸度等）的函数，运用图形技术将这种函数关系显示出来，（如萃取剂浓度、酸度等）的函数，运用图形技术将这种函数关系显示出来，以供我们凭借直觉的观察来选择试验设计中的最优化条件。以供我们凭借直觉的观察来选择试验设计中的最优化条件。第1页/共21页二、如何做响应面分析二、如何做响应面分析二、如何做响应面分析二、如何做响应面分析n n要构造响应面并进行分析以确定最优条件或寻找最优区域，首先必须通过试要构造响应面并进行分析以确定最优条件或寻找最优区域，首先必须通过试验获取大量的测量数据，并建立一个合适的数学模型（建模），然后再用此验获取大量的测量数据，并建立一个合适的数学模型（建模），然后再用此数学模型作图。数学模型作图。n n建模最常用和最有效的方法之一就是多元线性回归方法。对于非线性体系可建模最常用和最有效的方法之一就是多元线性回归方法。对于非线性体系可作适当处理化为线性形式。作适当处理化为线性形式。第2页/共21页二、如何做响应面分析二、如何做响应面分析二、如何做响应面分析二、如何做响应面分析n n设有个因素影响指标取值，通过试验测量，得到组设有个因素影响指标取值，通过试验测量，得到组试验数据。假设指标与因素之间的关系可用线性模型表试验数据。假设指标与因素之间的关系可用线性模型表示，则可将各系数写成矩阵式。示，则可将各系数写成矩阵式。n n应用最小二乘法即可求出模型参数矩阵，将矩阵代入原应用最小二乘法即可求出模型参数矩阵，将矩阵代入原假设的回归方程，就可得到响应关于各因素水平的数学假设的回归方程，就可得到响应关于各因素水平的数学模型，进而可以图形方式绘出响应与因素的关系图。模型，进而可以图形方式绘出响应与因素的关系图。n n模型中如果模型中如果只有一个因素（或自变量），响应（曲）面只有一个因素（或自变量），响应（曲）面是二维空间中的一条曲线是二维空间中的一条曲线；当有二个因素时，响应面是当有二个因素时，响应面是三维空间中的曲面三维空间中的曲面。第3页/共21页三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析 n n在化学量测实践中，一般不考虑三因素及三因素以上间在化学量测实践中，一般不考虑三因素及三因素以上间的交互作用。因此假设二因素响应（曲）面的数学模型的交互作用。因此假设二因素响应（曲）面的数学模型为二次多项式模型。为二次多项式模型。nn n n通过次测量试验（试验次数应大于参数个数，一般认通过次测量试验（试验次数应大于参数个数，一般认为至少应是它的倍），以最小二乘法估计模型各参数，为至少应是它的倍），以最小二乘法估计模型各参数，从而建立模型；从而建立模型；nn n n求出模型后，以两因素水平为坐标和坐标，以相应求出模型后，以两因素水平为坐标和坐标，以相应的响应为坐标作出三维空间的曲面（这就是因素响的响应为坐标作出三维空间的曲面（这就是因素响应曲面）。应曲面）。第4页/共21页三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析 n n应当指出，上述求出的模型只是最小二乘解，不一定与应当指出，上述求出的模型只是最小二乘解，不一定与实际体系相符，也即，计算值与试验值之间的差异不一实际体系相符，也即，计算值与试验值之间的差异不一定符合要求。因此，求出系数的最小二乘估计后，应进定符合要求。因此，求出系数的最小二乘估计后，应进行检验。行检验。n n一个简单实用的方法就是以响应的计算值与试验值之间一个简单实用的方法就是以响应的计算值与试验值之间的相关系数是否接近于或观察其相关图是否所有的点的相关系数是否接近于或观察其相关图是否所有的点都基本接近直线进行判别。都基本接近直线进行判别。第5页/共21页三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析三、二因素响应面分析 n n应当指出，上述求出的模型只是最小二乘解，不一定与应当指出，上述求出的模型只是最小二乘解，不一定与实际体系相符，也即，计算值与试验值之间的差异不一实际体系相符，也即，计算值与试验值之间的差异不一定符合要求。因此，求出系数的最小二乘估计后，应进定符合要求。因此，求出系数的最小二乘估计后，应进行检验。行检验。n n一个简单实用的方法就是以响应的计算值与试验值之间一个简单实用的方法就是以响应的计算值与试验值之间的相关系数是否接近于或观察其相关图是否所有的点的相关系数是否接近于或观察其相关图是否所有的点都基本接近直线进行判别。都基本接近直线进行判别。第6页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n在多因素数量处理试验的分析中，可以分析试验指标在多因素数量处理试验的分析中，可以分析试验指标(依依变量变量)与多个试验因素与多个试验因素(自变量自变量)间的回归关系，这种回归间的回归关系，这种回归可能是曲线或曲面的关系，因而称为响应面分析。可能是曲线或曲面的关系，因而称为响应面分析。n n例如农作物产量与例如农作物产量与N N、P P、K K的施肥量有关，可以通过回归的施肥量有关，可以通过回归分析建立产量与施肥要素间的回归关系，从而求得最佳分析建立产量与施肥要素间的回归关系，从而求得最佳施肥配方。施肥配方。第7页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n在回归分析中，观察值在回归分析中，观察值在回归分析中，观察值在回归分析中，观察值y y可以表述为：可以表述为：可以表述为：可以表述为：其中其中其中其中 是是是是 自变量的函数，自变量的函数，自变量的函数，自变量的函数，是是是是误差项。误差项。误差项。误差项。n n在响应面分析中，首先要得到回归方程在响应面分析中，首先要得到回归方程在响应面分析中，首先要得到回归方程在响应面分析中，首先要得到回归方程 然后通过对自变量然后通过对自变量然后通过对自变量然后通过对自变量 的合理取值，求得最优的合理取值，求得最优的合理取值，求得最优的合理取值，求得最优的值的值的值的值 ，这就是响应面分析的目的。，这就是响应面分析的目的。，这就是响应面分析的目的。，这就是响应面分析的目的。第8页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n 例例例例 有一个大麦氮磷肥配比试验，施氮肥量为每亩尿有一个大麦氮磷肥配比试验，施氮肥量为每亩尿有一个大麦氮磷肥配比试验，施氮肥量为每亩尿有一个大麦氮磷肥配比试验，施氮肥量为每亩尿素素素素0 0，3 3，6 6，9 9，1212，1515，18kg 718kg 7个水平，施磷肥量为每个水平，施磷肥量为每个水平，施磷肥量为每个水平，施磷肥量为每亩过磷酸钙亩过磷酸钙亩过磷酸钙亩过磷酸钙0 0，7 7，1414，2121，2828，3535，42kg 742kg 7个水平，共个水平，共个水平，共个水平，共4949个处理组合，试验结果列于下表个处理组合，试验结果列于下表个处理组合，试验结果列于下表个处理组合，试验结果列于下表1 1，试作产量对于氮、，试作产量对于氮、，试作产量对于氮、，试作产量对于氮、磷施肥量的响应面分析。磷施肥量的响应面分析。磷施肥量的响应面分析。磷施肥量的响应面分析。第9页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 第10页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n对于表对于表1 1的数据可以采用二元二次多项式拟合，那么产的数据可以采用二元二次多项式拟合，那么产量可表示为：量可表示为：n n其中其中N Ni i、P Pj j、分别表示分别表示N N、P P施用量和误差，按此模型施用量和误差，按此模型的的方差分析见表方差分析见表2 2。n n从表从表2 2结果看，结果看，b b2 2和和b b3 3这两个偏回归系数不显著，应该这两个偏回归系数不显著，应该将模型缩减，逐步去掉不显著的回归系数，将模型缩减，逐步去掉不显著的回归系数，结果见表结果见表3 3。得到的模型为：得到的模型为：第11页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n使用该模型分析的结果为表使用该模型分析的结果为表3 3，从，从表表3 3中可以看出，中可以看出，b b1 1，b b4 4，b b5 5达到极显著水平，达到极显著水平，b b2 2接近达到显著性，只有接近达到显著性，只有b3b3达达不到显著水平。不到显著水平。n n该模型的回归变异占总变异的该模型的回归变异占总变异的98%98%，因此可以较好地说，因此可以较好地说明施用明施用N N、P P对产量的影响。对产量的影响。n n二元二次多项式回归的回归系数二元二次多项式回归的回归系数及其显著性测验见表及其显著性测验见表4 4。第12页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n 由表由表4 4，可以列出产量对，可以列出产量对N N、P P施用量的回归方程为：施用量的回归方程为：n n由回归方程，可以作出产量对由回归方程，可以作出产量对N N、P P施用量的施用量的响应曲面响应曲面图图1 1。第13页/共21页四、响应面分析实例四、响应面分析实例四、响应面分析实例四、响应面分析实例 n n分别对回归方程求对分别对回归方程求对N N和和P P的偏导数，并令偏导数等于的偏导数，并令偏导数等于0 0，可以求得极值：，可以求得极值：N=13.87(kg)P=21.61(kg)l因而由回归方程估计得尿素施用量为13.87kg，过磷酸钙施用量为21.61kg时产量最高。l响应面分析中通过回归方程进行预测时一般不能超过自变量的取值范围，例如氮肥的取值范围为0至18kg/亩，而磷肥的取值范围为0至42kg/亩。第14页/共21页第15页/共21页第16页/共21页第17页/共21页第18页/共21页第19页/共21页图1 大麦产量对于氮、磷肥的响应面图 第20页/共21页