第五讲，最优化方法.pptVIP

最优化方法 最优化方法 改变单因子法 在一实际测试体系中，有许多影响因素(称之为因子)，习惯的做法是每次改变一种因子，而其他因子固定不变，以寻找这种因子对于体系响应的最适范围。 改变单因子法 如用分光光度法测定铜，影响测定的主要因素，如Cu的浓度、发色试剂的浓度、酸度及温度等。若观察溶液酸度的影响，则可固定Cu的浓度和溶液温度，并用某一固定浓度的发色试剂测试不同酸度下的吸光度。以类似的方法可以测试Cu浓度和温度的最适范围等，从而得到Cu的最佳测定条件。 改变单因子法的前提 各因子间彼此是不相关的。 假若因子间相关，运用这种方法所得最佳测试条件有可能不是实际上最佳的。 单纯形法 基本单纯形法 单纯形为一几何图形，在n维空间中，它由n+1个点所构成。 当n＝2时，即在二维平面上，单纯形为一三角形。 当n＝3时，即在三维空间，单纯形为一四面体。 当n＞4时，即在高维空间，单纯形为一超四面体。 基本单纯形的操作方法 以二因子单纯形为例 首先选定3个点，并且比较在这3个点处系统的响应(函数值)，我们的目的是寻求系统的最佳响应，故系统的响应值越大，此点越好；反之，响应值越小，此点越坏。 响应最差的点以W表示，次差点以N表示，最好点以B表示。如下图所示，即为一单纯形BNW。 根据这3点的好坏情况，可以取R点，使 WR＝2WP (P为NB中点，称为重心) 基本单纯形的操作方法 我们称点R为W关于P点的反射点，这种做法叫反射。 将点W去掉，则得一新的单纯形为BNR。比较3点的“好”，“坏”，仍可记得新的单纯形为BAW。 重复如上操作，则可使单纯形逐步移向系统响应的最大区域。 坐标计算 P点坐标 P=（N+B）/K (K为因子数) R点坐标 R=P+(P-W) N点 4.5 6.33 B点 ＋ 7 2 P点 （11.5 8.33）÷2=（5.75 4.165） 5.75 4.165 W点 －2 2 3.75 2.165 ＋5.75 4.165 R点 9.5 6.33 基本单纯形的移动规则 除初始的单纯形外，每次单纯形的移动均是在观察各点处的响应之后进行。 单纯形的移动均是在舍弃系统的响应最坏点之后进行。 在新单纯形中，若反射点为最坏点，此时将反射点保留，而将次最坏点舍弃，即取次最好点N的反射点R，以构成新的单纯形。 基本单纯形的移动应遵循如下规则 若某一顶点在运用规则2前得以保留，则在保留的点处重复进行实验 若一新的顶点超出允许的范围，则将此点赋一最坏的响应值，而强使此点回到正常范围。 基本单纯形法有如下不足： 单纯形在移动过程中不能够加速，当单纯形的起始位置远离最优处时，此缺点则更明显。另外当沿着响应曲面的“脊”移动时，仍需进行大量实验，因为单纯形走过的为一锯齿形。 运用规则4时是将一最坏的值赋予给超出允许范围的顶点而强行使之回到正常范围，并且舍弃次最坏点。若次最坏点的反射点处系统的响应更差，则有可能使单优形移不到最优区域，此时，单纯形被迫停止。 改良单纯形(变步长)法 在改良单纯形中引入两个新概念，即扩展和压缩 它们可使单纯形在响应曲面的陡峭处改变自己移动的步长。 在基本单纯形中，有的反射点系统的响应改进并不大(尚不知是否由于噪声所致)，则也要构成一新的单纯形。在改良单纯形中，则沿着反射点R的方向进行扩展以试探系统响应的好坏，由此，扩展可以增强单纯形的抗噪声能力。 改良单纯形的移动规则 规则1 若在反射点R处系统的响应优于在B点处的响应，则将反射点扩展到E。顶点E的坐标为 E＝P+2×(P-W) 若在E处系统的响应优于在B点处的响应，则保留点E，并构成一新的单纯形BNE。若在E点系统的响应差于B点，此时我们说扩展失败，则BNR为新的单纯形。 改良单纯形的移动规则 若在顶点R处的响应优于在N点的响应，但差于在B处的响应，则此时既不扩展也不压缩，以BNR为单纯形重新计算。 若反射点R处系统的响应次于N点，此时反射点R不能构成新单纯形，应该压缩。如右图，压缩分为两种情况。 改良单纯形的移动规则 若在反射点R处系统的响应差于在N点处的响应，但优于在W点处的响应，则压缩应该靠近R点处，此时压缩点称为CR。CR点为正压缩，坐标为 CR＝P+0.5（P - W) 改良单纯形的移动规则 若在反射点R处系统的响应差于在W处，则压缩应靠近W点处。 此时压缩点称为CW。Cw的坐标为 CW＝P - 0.5(P - W) 改良单纯形的移动规则 若在反射点CW