第五章 优选法.pptVIP

第五章 优选法 5.1单因素优选法 Y＝f（x），x在（a，b）内求y的极值的问题。但f（x）形式未知。 5.1.1来回调试法 假定f（x）在（a，b）内仅有一个峰， 取ax1x2b，看试验结果，有三种情况： y1y2、y1＝y2和y1y2 y1y2 时 说明峰不在（x2，b）范围内，再在（a，x2）内取x3，对应试验值y3。 y1y2 时 说明峰不在（a，x1）范围内再在（x1，b）内取x3，对应试验值y3。 y1＝y2时 说明峰在（x1，x2）范围内，再在（x1，x2）内取x3、x4，对应试验值y3、y4。 每次都是比较中间两点对应的y值，舍弃y值较低的x及其相邻的边界之间的范围，n次后y达到最优，或三次y值相差不大为止。 但是如何选点x1、x2的位置，才能用最少的试验次数最快达到最优结果？ 5.1.2黄金分割法（0.618法） x1、x2的位置取（a，b）的0.618和0.382处。 例题：钢中加入某种元素能够提高钢的硬度，1吨钢中该元素的加入量在1000～2000g之间，如何通过试验尽快找到最佳值？ 5.1.3二分法 当我们做一次试验就能得到下一步试验的方向时，可以用二分法更快地接近目标。 例：某厂生产溴樟脑，原来按1：1氯仿：樟脑使樟脑溶解后反应，产品合格，但需要蒸馏回收溶剂氯苯，才能使产物结晶。故减少氯仿用量0.5：1,合格；0.25：1合格；0.13：1合格；0.06：1溶解不完全，所以用0.13：1投产，大大降低成本，而且去掉了回收工序，产品自然结晶。 5.1.4抛物线法 先做三次试验x1、x2、x3对应得到y1、y2、y3，用y＝a0＋a1x＋a2x2，回归得到a0、a1、a2，求出极值对应的x4，预测y4’，由x4点试验得y4，再用y1～y4中最大的三个y重新拟合，求极值，直到y达到要求或y与预测值基本吻合。 需要注意的问题 某些因素A不能大规模调整时，可以先在原来的生产条件下在A减小的方向上跨一小步，看效果，好就继续；不好就反向跨一小步；遵循先小步再大步再小步的原则，逐渐接近最优。 多峰情况：做分布较均匀的疏松试验，看是否多峰，把每个峰都找出来。 5.2多因素优选法 5.2.1单因素轮换法 固定x2、x3……xn，只改变x1，按单因素寻优法找到x1的最佳值，x1固定在最佳条件下，改变x2，找x2的最佳……，找到xn的最佳后，返回来再找x1的最佳…… 直到x1……xn的条件都不再改变。 5.2.2对开法 axb，在(a+b)/2的位置用单因素找最大值P， cyd，在(c+d)/2的位置用单因素找最大值Q，有三种情况：PQ,PQ,P=Q PQ舍弃下半部分 。 PQ,舍弃左半边 P=Q，只留PQ共同所在象限 如此再找剩下部分的中线，找最佳值，再比较，舍弃一半……直到P、Q非常接近，可以认为是同一点。 5.2.3旋升法 5.2.4平行线法 两个因素中，有一个不易调整，将不易调整的Y分别取0.618和0.382处，找到各自条件下的x的最佳点P(0.618)、Q(0.382)，比较P、Q两点对应的Y值，舍弃小值以外的部分，再取剩余部分的0.618处…… 5.3单纯形法 单纯形：在n维空间里，由n+1个顶点构成的凸多面体。若凸多面体的边长相等，则为正规单纯形。如：二维空间的单纯形为三角形，正规单纯形为等边三角形；三维空间的单纯形为四面体，正规单纯形为正四面体；n维空间内的单纯形为超四面体。 单纯形原理 如二因素实验，首先选三个点，并比较其响应值，响应最差的以W表示；次差的为N；最好的记为B。P为NB的中心，称为重心。称R为W关于P点的反射点。这种做法叫反射，将点W去掉，得到一个新的单纯形BNR，仍按好、次坏、最坏记为BNW重复操作。可使单纯形移向响应最大的区域。 P=(N+B)/k k一般取2 R=P+(P-W)×α α一般取1 5.3.1初始单纯形的构造 a.给定某一顶点和步长，构造正规初始单纯形 首先根据化学知识和经验确定初始顶点x1(x11、x12…x1n)，并对每一个因素选一步长aj，为保证后继实验中各单纯形边长相等，这个步长与各因素中的变化范围的比值应相同。 b.黄金分割法构造初始单纯形 前提：知道各因素的范围。首先将各因素按黄金分割分为两个水平：0.382水平和0.618水平，计算出A1/2=(A1+A2)/2。按表构造初始单纯形。 c.根据均匀设计表构造初始单纯形。 例题 设有一光度分析体系，现需研究A：盐酸浓度0.1-1.0mol/L,B：温度：25-85℃，C:反应时间：10-70min，D:显色剂浓度1.0-4.0对该反应的影响，用上述三种方法分别构造初始单纯形。 构造正规单纯形 设步长取因素水平范围的2/3，则aj＝（0.6，40，40，2.0） 选初始顶点为（0.2