单纯形优化法.ppt

《单纯形优化法.ppt》由会员分享，可在线阅读，更多相关《单纯形优化法.ppt（55页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、关于单纯形优化法现在学习的是第1页，共55页7-1 概述 单纯形法是一种优化设计方法单纯形法是一种优化设计方法 和正交试验相比的特点：和正交试验相比的特点：计算简便计算简便不受因素数的限制不受因素数的限制因素数的增加不会导致试验次数大量增加因素数的增加不会导致试验次数大量增加 它属于非线性动态调优过程它属于非线性动态调优过程现在学习的是第2页，共55页 发展简史发展简史19621962年，年，SpendleySpendley提出基本单纯形法提出基本单纯形法19651965年，年，NelderNelder等提出改进单纯形法等提出改进单纯形法之后，之后，RouthRouth提出加权形心法与控制加权

2、形心法提出加权形心法与控制加权形心法现在学习的是第3页，共55页7-2 7-2 基本单纯形基本单纯形 一、双因素基本单纯形法一、双因素基本单纯形法如果我们有一个试验设计，只选有两个影响因素，即如果我们有一个试验设计，只选有两个影响因素，即因素数为因素数为2 2。分别取值。分别取值a a1 1和和a a2 2作为试验的初点。记为作为试验的初点。记为A(aA(a1 1,a,a2 2) )。对其余两个点分别设为。对其余两个点分别设为B B和和C C，再设三角形，再设三角形的边长为的边长为a(a(步长步长) )。那么。那么B B、C C点就可以计算出来点就可以计算出来现在学习的是第4页，共55页假设假

3、设ABAB、 ACAC、BCBC间距均为间距均为 ，等边三角形可以算出，等边三角形可以算出B B点为：点为： B=(aB=(a1 1+p, a+p, a2 2+q)+q)根据对称性可知：根据对称性可知： C=(aC=(a1 1+q, a+q, a2 2+p)+p)可以根据等边三角形性质解得：可以根据等边三角形性质解得：312 2(91)312 2qapa现在学习的是第5页，共55页a2+pa2+qa2a1+pa1+qa1因素2因素1ABCDEo现在学习的是第6页，共55页 由由A A、B B、C C三点构成得单纯形称为初始单纯形三点构成得单纯形称为初始单纯形首先在首先在A A、B B、C C三

4、点下分别试验，得出三个响应值，比较其大小，找三点下分别试验，得出三个响应值，比较其大小，找出最坏响应值的点称为坏点出最坏响应值的点称为坏点 此处设此处设A A为坏点，去掉为坏点，去掉A A点并取点并取A A的对称点的对称点D D点作为新试验点，比点作为新试验点，比较较B B、C C、D D三点响应值的好坏三点响应值的好坏此处设此处设C C为坏点，去点为坏点，去点C C点，取其反点点，取其反点E E，此时，此时C C、D D、E E三点又构成三点又构成新的单纯形新的单纯形 重复以上结果，最终达到优化试验的目的重复以上结果，最终达到优化试验的目的现在学习的是第7页，共55页 二、新试验点的计算方法

5、二、新试验点的计算方法以初始单纯形以初始单纯形A A、B B、C C为例，设为例，设A A为坏点，为坏点，A A应该应该去掉，求其反射点去掉，求其反射点D D，此时，此时A(aA(a1 1,a,a2 2) )、B=(aB=(a1 1+p, a+p, a2 2+q)+q)、C=(aC=(a1 1+q, a+q, a2 2+p)+p)D=B+C-A=(aD=B+C-A=(a1 1+p+q,a+p+q,a2 2+p+q)+p+q)E=B+D-C=(aE=B+D-C=(a1 1+2p,a+2p,a2 2+2q)+2q)即：即： 新试验点新试验点 留下各点之和留下各点之和 去掉点去掉点(9-8)(9-8

6、)现在学习的是第8页，共55页 三、多因素基本单纯形三、多因素基本单纯形 设有设有n n个因素个因素n n1 1个定点构成的个定点构成的n n维空间单纯形，设有一维空间单纯形，设有一点点A=(aA=(a1 1, a, a2 2, , a a3 3, , a an n) )，步长为，步长为a a 则其余各点为：则其余各点为： B=(aB=(a1 1+p,a+p,a2 2+q,a+q,a3 3+q, a+q, an n+q)+q) C=(a C=(a1 1+q,a+q,a2 2+p,a+p,a3 3+q, a+q, an n+q)+q) (n)=(a (n)=(a1 1+q,a+q,a2 2+q,

7、 a+q, an-1n-1+p, a+p, an n+q)+q) (n+1)=(a (n+1)=(a1 1+q,a+q,a2 2+q,a+q,a3 3+q, a+q, an n+p)+p)现在学习的是第9页，共55页 其中其中112(98)112nnqannpan 新点计算新点计算新坐标点新坐标点2n留下点的坐标和留下点的坐标和/n 去掉点坐标去掉点坐标 (9-11)现在学习的是第10页，共55页 四、四、n n，p p，q q取值对应表取值对应表 由由(9-8) (9-8) 我们可以算出我们可以算出n n取不同值的取不同值的p p、q q的取值的取值112(98)112nnqannpan 现

8、在学习的是第11页，共55页n2345678pqn9101112131415pq0.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取值对应表现在学习的是第12页，共55页五、小结五、小结 用前面的例子，对两因素问题用前面的例子，对两因素问题A A、B B、C C构成初始单纯形，在此三点上构成初始单纯形，

9、在此三点上进行试验进行试验 规则规则1 1：去掉最坏点，用其对称反射点作新试点：去掉最坏点，用其对称反射点作新试点 例例A A、B B、C C中，中，A A为最坏点，去掉为最坏点，去掉A A点并取点并取A A的对称点的对称点D D点作为新试验点。点作为新试验点。 D D 留下各点之和留下各点之和 去掉点去掉点 B BC CA A 在在B B、C C、D D三角形中继续使用规则三角形中继续使用规则1 1，如果，如果C C为坏点，去点为坏点，去点C C点，取其反点点，取其反点E E，此时此时C C、D D、E E三点又构成新的单纯形。三点又构成新的单纯形。 如果最坏点为如果最坏点为D D那么对称点

10、就会返回到与那么对称点就会返回到与A A重合，此时改用规则重合，此时改用规则2 2现在学习的是第13页，共55页规则规则2 2：去掉次坏点，用其对称反射点作新试点对称计算公式与前面：去掉次坏点，用其对称反射点作新试点对称计算公式与前面相同相同 经过反复使用后，如果有一个点老是保留下来，必须使用规则经过反复使用后，如果有一个点老是保留下来，必须使用规则3 3 规则规则3 3：重复、停止和缩短步长：重复、停止和缩短步长 一般一个点劲一般一个点劲3 3次单纯形后仍未被淘汰，它可能是一个很好点，也可能次单纯形后仍未被淘汰，它可能是一个很好点，也可能是偶然性或试验误差导致的假象。是偶然性或试验误差导致的

11、假象。 此时需要重复试验：结果不好，淘汰；结果已很满意则停止试验此时需要重复试验：结果不好，淘汰；结果已很满意则停止试验 反之则以它为起点缩短步长，继续试验反之则以它为起点缩短步长，继续试验现在学习的是第14页，共55页 六、特殊方法六、特殊方法前面介绍的单纯形是正规的，任意两点间的距离前面介绍的单纯形是正规的，任意两点间的距离一样，实际上，这个要求可以不要。尤其是由于各个一样，实际上，这个要求可以不要。尤其是由于各个因素所取的量纲不一样（例如一个因素是温度（因素所取的量纲不一样（例如一个因素是温度（），），另一个因素是时间（秒）。即使量纲一样所取的单位另一个因素是时间（秒）。即使量纲一样所取

12、的单位也可以不一样。也可以不一样。现在学习的是第15页，共55页 （一）直角单纯形法（一）直角单纯形法我们考虑双因素模型，开始不从正三角形出发，而是从一个直我们考虑双因素模型，开始不从正三角形出发，而是从一个直角三角形出发，其顶点取值如下：角三角形出发，其顶点取值如下： =(a=(a1 1,a,a2 2) ) =(a=(a1 1+p+p1 1,a,a2 2) ) =(a=(a1 1,a,a2 2+p+p2 2) )用图表示如下用图表示如下现在学习的是第16页，共55页a2+p2a2a1+2p1a1+p1a1因素2因素1现在学习的是第17页，共55页同样比较三个顶点响应值的结果，若同样比较三个顶

13、点响应值的结果，若最坏，则新点最坏，则新点就就用对称公式用对称公式 = =+ +- -=(a=(a1 1+p+p1 1,a,a2 2+p+p2 2) )在得到在得到点后，再用点后，再用、三点试验，比较其结果，若三点试验，比较其结果，若最最坏，则取其对称点坏，则取其对称点做新试验点做新试验点 = =+ +- -=(a=(a1 1+2p+2p1 1,a,a2 2) )、构成一个新单纯形，比较其结果，若构成一个新单纯形，比较其结果，若最坏，则用规则最坏，则用规则2 2去掉次坏点，若次坏点为去掉次坏点，若次坏点为，则新点，则新点 = =+ +- -=(a=(a1 1+2p+2p1 1,a,a2 2-p

14、-p2 2) )如此等等，有时还会使用规则如此等等，有时还会使用规则3 3，直至结果满意为止。，直至结果满意为止。现在学习的是第18页，共55页一般在任意一般在任意n n个因素时个因素时 =(a=(a1 1, a, a2 2, , a a3 3, , a an n) ) =(a=(a1 1+p+p1 1,a,a2 2,a,a3 3, a, an n) ) =(a=(a1 1,a,a2 2+p+p2 2,a,a3 3, a, an n) ) (n)=(a (n)=(a1 1,a,a2 2, a, an-1n-1+p+pn-1n-1, a, an n) ) (n+1)=(a (n+1)=(a1 1

15、,a,a2 2,a,a3 3, a, an n+p+pn n) )现在学习的是第19页，共55页 （二）、双水平单纯形法（二）、双水平单纯形法现在学习的是第20页，共55页7-3 7-3 改进单纯形法改进单纯形法 为了解决优化结果精度和优化速度的矛盾，可以采用可变步为了解决优化结果精度和优化速度的矛盾，可以采用可变步长推移单纯形，此即改进单纯形法，既能加快优化速度，又长推移单纯形，此即改进单纯形法，既能加快优化速度，又能获得较好的优化精度。能获得较好的优化精度。 改进单纯形法是改进单纯形法是19651965年年J JA ANelderNelder等提出来的，它等提出来的，它是在基本单纯形法的基

16、础上引入了是在基本单纯形法的基础上引入了反射反射、扩大扩大、收缩收缩与与整体收缩整体收缩规则，变固定步长为规则，变固定步长为可变步长可变步长，较好地解决了，较好地解决了优化速度与优化精度之间的矛盾，是各种单纯形优化优化速度与优化精度之间的矛盾，是各种单纯形优化法中应用最广泛的一种单纯形优化方法。法中应用最广泛的一种单纯形优化方法。现在学习的是第21页，共55页 两因素单纯形的推移过程两因素单纯形的推移过程因素1因素2BACDEONA改进单纯形ND现在学习的是第22页，共55页 单纯形的整体收缩单纯形的整体收缩因素1因素2ABCCA现在学习的是第23页，共55页 在单纯形的推移过程中，新实验点在

17、空间的位置坐标在单纯形的推移过程中，新实验点在空间的位置坐标按以下方法计算：按以下方法计算：(919)0naa留下各点的坐标和新试点的坐标(1+a)去掉点的坐标式中 是大于 的系数现在学习的是第24页，共55页 讨论：讨论： a a1 1，此时，此时(9-19)(9-19)式变差基本单纯形中新点的计算公式，此时新试式变差基本单纯形中新点的计算公式，此时新试验点为去掉点的等距离反射点，这时改进单纯形又变成了基本单纯验点为去掉点的等距离反射点，这时改进单纯形又变成了基本单纯形形 a1a1，按基本单纯形法（，按基本单纯形法（a a1 1）计算出新点后，对新试验点做试验）计算出新点后，对新试验点做试验

18、得出新试验点的响应值。如果新点的响应值好，说明我们搜索方向得出新试验点的响应值。如果新点的响应值好，说明我们搜索方向正确，可以进一步沿正确，可以进一步沿ADAD搜索。因此取搜索。因此取a1a1，称为，称为扩大扩大。如果扩大点。如果扩大点E E不如反射点不如反射点D D好，则好，则“扩大扩大”失败，仍采用失败，仍采用D D，由反射点何留下点，由反射点何留下点构成的单纯形构成的单纯形BCDBCD继续优化继续优化现在学习的是第25页，共55页 -1a0-1a0，按（，按（a a1 1）计算出来的反射点）计算出来的反射点D D的响应值最坏，此时采用的响应值最坏，此时采用- -1a01a0（称为（称为内

19、收缩内收缩）计算新试验点，此时形成新的单纯形）计算新试验点，此时形成新的单纯形BNBNA AC C 0a10a1，按基本单纯形法（，按基本单纯形法（a a1 1）计算除反射点）计算除反射点D D响应值最坏。但比去响应值最坏。但比去掉点掉点A A响应值好。此时采用响应值好。此时采用0a10a1，称为，称为收缩收缩，新试点仍按，新试点仍按(9-19)(9-19)式计式计算，此时形成新的单纯形算，此时形成新的单纯形BCNBCND D现在学习的是第26页，共55页 如果去掉点与其反射点连线如果去掉点与其反射点连线ADAD方向上所有点的响应值都比去掉点方向上所有点的响应值都比去掉点A A坏，则不能沿此方

20、向搜索。这时应以单纯形中最好点为初点，到其它各点的坏，则不能沿此方向搜索。这时应以单纯形中最好点为初点，到其它各点的一半为新点，构成新的单纯形一半为新点，构成新的单纯形BACBAC进行优化。此时步长减半，称为进行优化。此时步长减半，称为“整体收缩整体收缩”现在学习的是第27页，共55页7-4 7-4 加权形心法加权形心法 基本单纯形和改进单纯形都是采用去掉点的反射方向基本单纯形和改进单纯形都是采用去掉点的反射方向为新试验点的搜索方向，这就意味着，去掉点的反射为新试验点的搜索方向，这就意味着，去掉点的反射方向作为近似的优化方向，就是梯度变化最大的方向方向作为近似的优化方向，就是梯度变化最大的方向

21、 实际上，这个方向是一个近似的梯度最大方向，这样实际上，这个方向是一个近似的梯度最大方向，这样的搜索结果可能导致搜索次数的增加和搜索结果精度的搜索结果可能导致搜索次数的增加和搜索结果精度的降低的降低 为了解决这个问题，提出了加权形心法，加权形心法利为了解决这个问题，提出了加权形心法，加权形心法利用加权形心代替单纯的反射形心，使新点的搜索方向更用加权形心代替单纯的反射形心，使新点的搜索方向更接近实际的最优方向接近实际的最优方向现在学习的是第28页，共55页因素1因素2BCOEEO形心点O和加权形心点O现在学习的是第29页，共55页 如图，使如图，使WW、B B、C C三个顶点组成的一个二因素的优

22、化过三个顶点组成的一个二因素的优化过程的一个单纯形，并知程的一个单纯形，并知WW点的响应最坏点的响应最坏，B B的响应最好的响应最好。 如果搜索优化过程中函数不出现异常，那么搜索最优如果搜索优化过程中函数不出现异常，那么搜索最优点的方向明显应当更靠近点的方向明显应当更靠近WBWB的方向，而不是靠近的方向，而不是靠近WCWC的方向。因此可以通过加权的办法来使搜索的方向由原的方向。因此可以通过加权的办法来使搜索的方向由原来的来的WEWE（反射方向）（反射方向）变为变为WEWE方向（加权方向）方向（加权方向），此，此时用加权形心点时用加权形心点O O 代替反射形心点代替反射形心点O O现在学习的是第

23、30页，共55页3131 ()( )(920) ()( )i()i( )iiiiiiiR pRpOR pRpR pR加权形心点其中为第 点的坐标为的 点的响应值。为最坏点的响应值。现在学习的是第31页，共55页11()11()n ()( )(921) ()( ) niiiiniiiR pRpOR pROO同样对于 因素的加权形心点计算如下：然后将代替改进单纯形法中的形心点，即成为加权形心法。现在学习的是第32页，共55页7-5 7-5 单纯形优化的参数选择单纯形优化的参数选择 在试验中，我们只研究优化条件，可用基本单纯形法时，首在试验中，我们只研究优化条件，可用基本单纯形法时，首先必须确定研究

24、的因素先必须确定研究的因素 由于单纯形法不受因素的限制，考察的因素可以相对的多些由于单纯形法不受因素的限制，考察的因素可以相对的多些 因素确定后，据分析仪器和试验要求，规定因素变化因素确定后，据分析仪器和试验要求，规定因素变化的上下限，据上下限的范围确定步长的大小。的上下限，据上下限的范围确定步长的大小。 步长较大，优化速度加快，精度较差；步长太小试验次数增步长较大，优化速度加快，精度较差；步长太小试验次数增多，优化速度变慢多，优化速度变慢现在学习的是第33页，共55页 一、试验指标一、试验指标 试验指标是用于衡量和考核试验响应的各种数值试验指标是用于衡量和考核试验响应的各种数值 在分析测试中

25、可将仪器响应值作为试验指标，但有时须转换称其在分析测试中可将仪器响应值作为试验指标，但有时须转换称其它的数量，试验指标是数量化的，以便直接比较结果的大小它的数量，试验指标是数量化的，以便直接比较结果的大小现在学习的是第34页，共55页 二、初始单纯形的构成二、初始单纯形的构成 本章第一节介绍的方法是根据本章第一节介绍的方法是根据初始点初始点和和步长步长来计算初始来计算初始单纯形的各个顶点，各因素的步长是相同的单纯形的各个顶点，各因素的步长是相同的 实际过程中，各因素步长和单位并不相同，利用这种方实际过程中，各因素步长和单位并不相同，利用这种方法会变得很麻烦，在实际应用中问题较多法会变得很麻烦，

26、在实际应用中问题较多 我们介绍下述两个构成初始单纯形的方法我们介绍下述两个构成初始单纯形的方法现在学习的是第35页，共55页（一）（一）longlong系数表法系数表法D.E.LongD.E.Long提出一种用系数表构成初始单纯形各顶点的方法，可提出一种用系数表构成初始单纯形各顶点的方法，可以解决试验设计中初始单纯形的构成问题以解决试验设计中初始单纯形的构成问题使用时把表中的对应值乘上该因素的步长后，再加到初始点坐标使用时把表中的对应值乘上该因素的步长后，再加到初始点坐标上上现在学习的是第36页，共55页1234567891010A B C D E F G H I J01.000.500.50

27、0.500.500.500.500.500.500.50000.8660.2890.2890.2890.2890.2890.2890.2890.2890000.8170.2040.1580.2040.2040.2040.2040.20400000.7910.1580.1580.1580.1580.1580.158000000.7750.1290.1290.1290.1290.1290000000.7640.1090.1090.1090.10900000000.7560.0940.0940.094000000000.7500.0830.0830000000000.7450.07500000000

28、000.742因素顶点Long系数表现在学习的是第37页，共55页 例：有一个二因素的设计过程，其初始点为例：有一个二因素的设计过程，其初始点为(10.0,2.0)(10.0,2.0)；步长为步长为1.01.0和和0.50.5，据，据LongLong系数表来计算其余两个顶点的系数表来计算其余两个顶点的坐标坐标 顶点顶点1 1： （10.0,2.010.0,2.0） 顶点顶点2 2： （10.0+1.0010.0+1.001.01.0，1.0+01.0+00.50.5） = =（11.0,2.011.0,2.0） 顶点顶点3 3： （10.0+0.510.0+0.51.01.0，2.0+0.86

29、62.0+0.8660.50.5） = =（10.510.5，2.4332.433）现在学习的是第38页，共55页 （二）均匀设计表法（二）均匀设计表法 利用利用LongLong系数表法所构成的初始单纯形各顶点在空间的系数表法所构成的初始单纯形各顶点在空间的分布是不均匀的，因此进行的是不均匀优化分布是不均匀的，因此进行的是不均匀优化 均匀设计表改变了这个缺点，使各顶点在空间均匀分布，这样均匀设计表改变了这个缺点，使各顶点在空间均匀分布，这样进行的优化就是进行的优化就是整体的均匀优化整体的均匀优化 据所选因素的因素数，确定一个比较合适的均匀表，使用时把表据所选因素的因素数，确定一个比较合适的均匀

30、表，使用时把表中的对应数值乘以响应因素的步长，加到初始点坐标上即可中的对应数值乘以响应因素的步长，加到初始点坐标上即可现在学习的是第39页，共55页 例：我们有一个四因素的优化过程，因此可以选用四因素的例：我们有一个四因素的优化过程，因此可以选用四因素的均匀设计表。设初点为均匀设计表。设初点为(1.0,1.0,1,0,1.0);(1.0,1.0,1,0,1.0);步长为步长为0.5,1.0,1.5,2.00.5,1.0,1.5,2.0。要求计算初始单纯形的各顶点。要求计算初始单纯形的各顶点现在学习的是第40页，共55页四因素均匀设表四因素均匀设表U5(54)现在学习的是第41页，共55页 顶点

31、顶点1 1：（1.0+1 1.0+1 0.5,1.0+2 0.5,1.0+2 1.0, 1.0+3 1.0, 1.0+3 0.5,1.0+4 0.5,1.0+4 2.02.0） = =（1.5,3.0,5.5,9.01.5,3.0,5.5,9.0） 顶点顶点2 2：（1.0+2 1.0+2 0.5,1.0+4 0.5,1.0+4 1.0, 1.0+1 1.0, 1.0+1 0.5,1.0+30.5,1.0+32.02.0） = =（2.0,5.0,2.5,7.02.0,5.0,2.5,7.0） 顶点顶点3 3：（1.0+3 1.0+3 0.5,1.0+10.5,1.0+11.0, 1.0+4

32、1.0, 1.0+4 0.5,1.0+2 0.5,1.0+2 2.02.0） = =（2.5,2.0,7.0,5.02.5,2.0,7.0,5.0）现在学习的是第42页，共55页 顶点顶点4 4：（1.0+4 1.0+4 0.5,1.0+3 0.5,1.0+3 1.0, 1.0+2 1.0, 1.0+2 0.5,1.0+1 0.5,1.0+1 2.02.0） = =（3.0,4.0,4.0,3.03.0,4.0,4.0,3.0） 顶点顶点5 5：（1.0+5 1.0+5 0.5,1.0+5 0.5,1.0+5 1.0, 1.0+5 1.0, 1.0+5 0.5,1.0+50.5,1.0+52.

33、02.0） = =（3.5.6.0,8.5.11.03.5.6.0,8.5.11.0）现在学习的是第43页，共55页 （三）单纯形的收敛（三）单纯形的收敛 单纯形收敛的检验办法：单纯形收敛的检验办法：在在n n因素的单纯形中，如果有一个点经因素的单纯形中，如果有一个点经n n1 1次单纯形仍为被淘汰，次单纯形仍为被淘汰，一般可以在此点收敛一般可以在此点收敛 这种检验方法未考虑到试验误差的存在，按助理统计或实际工作要这种检验方法未考虑到试验误差的存在，按助理统计或实际工作要求单纯形收敛准则应为：求单纯形收敛准则应为：|R(B)-R(|R(B)-R( )/R(B)|)/R(B)| 式中式中R(B)

34、R(B)和和R(R( ) )分别代表最好点分别代表最好点B B与最坏点与最坏点 的响应值，的响应值， 为试验为试验误差或预给定的允许误差误差或预给定的允许误差现在学习的是第44页，共55页应用举例应用举例现在学习的是第45页，共55页现在学习的是第46页，共55页现在学习的是第47页，共55页现在学习的是第48页，共55页现在学习的是第49页，共55页现在学习的是第50页，共55页现在学习的是第51页，共55页现在学习的是第52页，共55页（三）单纯形法与正交试验设计法比较（三）单纯形法与正交试验设计法比较 为了比较两种方法选出的最佳条件，用实验进行验证，为了比较两种方法选出的最佳条件，用实验进行验证，实验结果利于表实验结果利于表10-1910-19现在学习的是第53页，共55页现在学习的是第54页，共55页感谢大家观看感谢大家观看现在学习的是第55页，共55页