单纯形优化法

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单纯形优化法,7-1 概述,单纯形法是一种优化设计方法,和正交试验相比的特点：,计算简便,不受因素数的限制,因素数的增加不会导致试验次数大量增加,它属于非线性动态调优过程,开展简史,1962年，Spendley提出根本单纯形法,1965年，Nelder等提出改进单纯形法,之后，Routh提出加权形心法与控制加权形心法,7-2 根本单纯形,一、双因素根本单纯形法,如果我们有一个试验设计，只选有两个影响因素，即因素数为2。分别取值a1和a2作为试验的初点。记为A(a1,a2)。对其余两个点分别设为B和C，再设三角形的边长为a(步长)。那么B、C点就可以计算出来,假设AB、,AC、BC,间距均为,，等边三角形可以算出B点为：,B=(a,1,+p, a,2,+q),根据对称性可知：,C=(a,1,+q, a,2,+p),可以根据等边三角形性质解得：,a,2,+p,a,2,+q,a,2,a,1,+p,a,1,+q,a,1,因素2,因素1,A,B,C,D,E,o,由A、B、C三点构成得单纯形称为初始单纯形,首先在A、B、C三点下分别试验，得出三个响应值，比较其大小，找出最坏响应值的点称为坏点,此处设A为坏点，去掉A点并取A的对称点D点作为新试验点，比较B、C、D三点响应值的好坏,此处设C为坏点，去点C点，取其反点E，此时C、D、E三点又构成新的单纯形,重复以上结果，最终到达优化试验的目的,二、新试验点的计算方法,以初始单纯形A、B、C为例，设A为坏点，A应该去掉，求其反射点D，此时,A(a,1,a,2,)、,B=(a,1,+p, a,2,+q)、C=(a,1,+q, a,2,+p),D=B+C-A=(a,1,+p+q,a,2,+p+q),E=B+D-C=(a,1,+2p,a,2,+2q),即：新试验点留下各点之和去掉点(9-8),三、多因素根本单纯形,设有n个因素n1个定点构成的n维空间单纯形，设有一点A=(a1, a2, a3, an)，步长为a,那么其余各点为：,B=(a1+p,a2+q,a3+q, an+q),C=(a1+q,a2+p,a3+q, an+q),(n)=(a1+q,a2+q, an-1+p, an+q),(n+1)=(a1+q,a2+q,a3+q, an+p),其中,新点计算,新坐标点2n留下点的坐标和/n,去掉点坐标,(9-11),四、n，p，q取值对应表,由(9-8) 我们可以算出n取不同值的p、q的取值,n,2,3,4,5,6,7,8,p,q,n,9,10,11,12,13,14,15,p,q,0.966,0.943,0.926,0.911,0.901,0.892,0.883,0.259,0.236,0.219,0.204,0.194,0.185,0.176,0.878,0.872,0.865,0.861,0.855,0.854,0.848,0.171,0.165,0.158,0.154,0.148,0.147,0.141,n、q、p取值对应表,五、小结,用前面的例子，对两因素问题A、B、C构成初始单纯形，在此三点上进展试验,规那么1：去掉最坏点，用其对称反射点作新试点,例A、B、C中，A为最坏点，去掉A点并取A的对称点D点作为新试验点。,D留下各点之和去掉点BCA,在B、C、D三角形中继续使用规那么1，如果C为坏点，去点C点，取其反点E，此时C、D、E三点又构成新的单纯形。,如果最坏点为D那么对称点就会返回到与A重合，此时改用规那么2,规那么2：去掉次坏点，用其对称反射点作新试点对称计算公式与前面一样,经过反复使用后，如果有一个点老是保存下来，必须使用规那么3,规那么3：重复、停顿和缩短步长,一般一个点劲3次单纯形后仍未被淘汰，它可能是一个很好点，也可能是偶然性或试验误差导致的假象。,此时需要重复试验：结果不好，淘汰；结果已很满意那么停顿试验,反之那么以它为起点缩短步长，继续试验,六、特殊方法,前面介绍的单纯形是正规的，任意两点间的距离一样，实际上，这个要求可以不要。尤其是由于各个因素所取的量纲不一样例如一个因素是温度，另一个因素是时间秒。即使量纲一样所取的单位也可以不一样。,一直角单纯形法,我们考虑双因素模型，开场不从正三角形出发，而是从一个直角三角形出发，其顶点取值如下：,=(a1,a2),=(a1+p1,a2),=(a1,a2+p2),用图表示如下,a,2,+p,2,a,2,a,1,+2p,1,a,1,+p,1,a,1,因素2,因素1,同样比较三个顶点响应值的结果，假设最坏，那么新点就用对称公式,=+-=(a1+p1,a2+p2),在得到点后，再用、三点试验，比较其结果，假设最坏，那么取其对称点做新试验点,=+-=(a1+2p1,a2),、构成一个新单纯形，比较其结果，假设最坏，那么用规那么2去掉次坏点，假设次坏点为，那么新点,=+-=(a1+2p1,a2-p2),如此等等，有时还会使用规那么3，直至结果满意为止。,一般在任意n个因素时,=(a,1, a,2,a,3,a,n,),=(a,1,+p,1,a,2,a,3, ,a,n,),=(a,1,a,2,+p,2,a,3, ,a,n,), ,(n)=(a,1,a,2,a,n-1,+p,n-1,a,n,),(n+1)=(a,1,a,2,a,3, ,a,n,+p,n,),二、双水平单纯形法,7-3 改进单纯形法,为了解决优化结果精度和优化速度的矛盾，可以采用可变步长推移单纯形，此即改进单纯形法，既能加快优化速度，又能获得较好的优化精度。,改进单纯形法是1965年JANelder等提出来的，它是在根本单纯形法的根底上引入了反射、扩大、收缩与整体收缩规那么，变固定步长为可变步长，较好地解决了优化速度与优化精度之间的矛盾，是各种单纯形优化法中应用最广泛的一种单纯形优化方法。,两因素单纯形的推移过程,因素1,因素2,B,A,C,D,E,O,N,A,改进单纯形,N,D,单纯形的整体收缩,因素1,因素2,A,B,C,C,A,在单纯形的推移过程中，新实验点在空间的位置坐标按以下方法计算：,讨论：,a1，此时(9-19)式变差根本单纯形中新点的计算公式，此时新试验点为去掉点的等距离反射点，这时改进单纯形又变成了根本单纯形,a1，按根本单纯形法a1计算出新点后，对新试验点做试验得出新试验点的响应值。如果新点的响应值好，说明我们搜索方向正确，可以进一步沿AD搜索。因此取a1，称为扩大。如果扩大点E不如反射点D好，那么“扩大失败，仍采用D，由反射点何留下点构成的单纯形BCD继续优化,-1a0，按a1计算出来的反射点D的响应值最坏，此时采用-1a0称为内收缩计算新试验点，此时形成新的单纯形BNAC,0a1，按根本单纯形法a1计算除反射点D响应值最坏。但比去掉点A响应值好。此时采用0a1，称为收缩，新试点仍按(9-19)式计算，此时形成新的单纯形BCND,如果去掉点与其反射点连线AD方向上所有点的响应值都比去掉点A坏，那么不能沿此方向搜索。这时应以单纯形中最好点为初点，到其它各点的一半为新点，构成新的单纯形BAC进展优化。此时步长减半，称为“整体收缩,7-4 加权形心法,根本单纯形和改进单纯形都是采用去掉点的反射方向为新试验点的搜索方向，这就意味着，去掉点的反射方向作为近似的优化方向，就是梯度变化最大的方向,实际上，这个方向是一个近似的梯度最大方向，这样的搜索结果可能导致搜索次数的增加和搜索结果精度的降低,为了解决这个问题，提出了加权形心法，加权形心法利用加权形心代替单纯的反射形心，使新点的搜索方向更接近实际的最优方向,因素1,因素2,B,C,O,E,E,O,形心点O和加权形心点O,如图，使W、B、C三个顶点组成的一个二因素的优化过程的一个单纯形，并知W点的响应最坏，B的响应最好。,如果搜索优化过程中函数不出现异常，那么搜索最优点的方向明显应当更靠近WB的方向，而不是靠近WC的方向。因此可以通过加权的方法来使搜索的方向由原来的WE反射方向变为WE方向加权方向，此时用加权形心点O代替反射形心点O,7-5 单纯形优化的参数选择,在试验中，我们只研究优化条件，可用根本单纯形法时，首先必须确定研究的因素,由于单纯形法不受因素的限制，考察的因素可以相对的多些,因素确定后，据分析仪器和试验要求，规定因素变化的上下限，据上下限的范围确定步长的大小。,步长较大，优化速度加快，精度较差；步长太小试验次数增多，优化速度变慢,一、试验指标,试验指标是用于衡量和考核试验响应的各种数值,在分析测试中可将仪器响应值作为试验指标，但有时须转换称其它的数量，试验指标是数量化的，以便直接比较结果的大小,二、初始单纯形的构成,本章第一节介绍的方法是根据初始点和步长来计算初始单纯形的各个顶点，各因素的步长是一样的,实际过程中，各因素步长和单位并不一样，利用这种方法会变得很麻烦，在实际应用中问题较多,我们介绍下述两个构成初始单纯形的方法,一long系数表法,提出一种用系数表构成初始单纯形各顶点的方法，可以解决试验设计中初始单纯形的构成问题,使用时把表中的对应值乘上该因素的步长后，再加到初始点坐标上,1,2,3,4,5,6,7,8,9,10,10,A B C D E F G H I J,0,1.00,0.50,0.50,0.50,0.50,0.50,0.50,0.50,0.50,0.50,0,0,0.866,0.289,0.289,0.289,0.289,0.289,0.289,0.289,0.289,0,0,0,0.817,0.204,0.158,0.204,0.204,0.204,0.204,0.204,0,0,0,0,0.791,0.158,0.158,0.158,0.158,0.158,0.158,0,0,0,0,0,0.775,0.129,0.129,0.129,0.129,0.129,0,0,0,0,0,0,0.764,0.109,0.109,0.109,0.109,0,0,0,0,0,0,0,0.756,0.094,0.094,0.094,0,0,0,0,0,0,0,0,0.750,0.083,0.083,0,0,0,0,0,0,0,0,0,0.745,0.075,0,0,0,0,0,0,0,0,0,0,0.742,因素,顶点,Long系数表,例：有一个二因素的设计过程，其初始点为(10.0,2.0)；步长为和，据Long系数表来计算其余两个顶点的坐标,顶点1： ,顶点2： ，,=,顶点3： ，,=，,二均匀设计表法,利用Long系数表法所构成的初始单纯形各顶点在空间的分布是不均匀的，因此进展的是不均匀优化,均匀设计表改变了这个缺点，使各顶点在空间均匀分布，这样进展的优化就是整体的均匀优化,据所选因素的因素数，确定一个比较适宜的均匀表，使用时把表中的对应数值乘以响应因素的步长，加到初始点坐标上即可,例：我们有一个四因素的优化过程，因此可以选用四因素的均匀设计表。设初点为(1.0,1.0,1,0,1.0);步长为。要求计算初始单纯形的各顶点,四因素均匀设表U,5,(5,4,),顶点1：,=,顶点2：,=,顶点3：,=,顶点4：,=,顶点5：,=,三单纯形的收敛,单纯形收敛的检验方法：,在n因素的单纯形中，如果有一个点经n1次单纯形仍为被淘汰，一般可以在此点收敛,这种检验方法未考虑到试验误差的存在，按助理统计或实际工作要求单纯形收敛准那么应为：,|R(B)-R()/R(B)|,式中R(B)和R()分别代表最好点B与最坏点的响应值， 为试验误差或预给定的允许误差,应用举例,三单纯形法与正交试验设计法比较,为了比较两种方法选出的最正确条件，用实验进展验证，实验结果利于表10-19,谢谢,