固体退火过程组合优化问题物理系统中的一个状态组合优化Read课件.ppt

1、l局部搜索方法l模拟退火算法l遗传算法l很多问题属于优化问题，或者可以转化为优化问题l如TSP问题，皇后问题l设x是决策变量，D是x的定义域，f(x)是指标函数，g(x)是约束条件集合。则优化问题可以表示为，求解满足g(x)的f(x)最小值问题。l如果在定义域D上，满足条件g(x)的解是有限的，则优化问题称为组合优化问题。)(|)(minxgxfDxl对于组合优化问题，由于其可能的解是有限的，当问题的规模比较小时，总可以通过枚举的方法获得问题的最优解，但当问题的规模比较大时，就难于求解了。l常用的算法复杂度函数)(),!(),(),2(),(),log(),(),(loglog2nnnnOnO

2、nOOnOnnOnOnO 输入量输入量n复杂性函数复杂性函数10203040100n10ns20ns30ns40ns100nsnlogn10ns26.0ns44.3ns64.1ns200nsn2100ns400ns900ns1.6us10us2n1.0us1.0ms1.1s18.3min4.0世纪世纪n!3.6ms77.1年年8.41013世纪世纪2.61029世纪世纪3.010139世纪世纪时间复杂性函数比较（10亿次/秒）l旅行商问题l背包问题l装箱问题l.l寻求在可以接受的时间内得到满意解的方法l邻域，简单的说就是一个点附近的其他点的集合。l在距离空间，邻域就是以某一点为中心的圆。l组合

3、优化问题的定义：l设D是问题的定义域，若存在一个映射N，使得：则称N(S)为S的邻域。DSNDSN2)(:lS=Si表示一个可能解，其中Si表示在第i行，第Si列有一个皇后。l如四皇后问题的一个解：S=(2,4,1,3)Q QQ Q l定义映射N为棋盘上任意两个皇后的所在行或列进行交换，即S中任意两个元素交换位置。l例：当S=(2,4,1,3)时，其邻域为：lN(S)=(4,2,1,3),(1,4,2,3),(3,4,1,2),(2,1,4,3),(2,3,1,4),(2,4,3,1)l用一个城市的序列表示一个可能的解。l通过交换两个城市的位置获取S的邻居 l例：简单交换方法 设S=(x1,x

4、2,xi-1,xi,xi+1,xj-1,xj,xj+1,xn)则通过交换xi和xj两个城市的位置可以得到S的一个邻居：S=(x1,x2,xi-1,xj,xi+1,xj-1,xi,xj+1,xn)x1x2xnxj+1xjxj-1xi-1xixi+1x1x2xnxj+1xjxj-1xi-1xixi+1l例：逆序交换方法 设xi、xj是选取的两个城市，所谓的逆序交换方式是指，通过逆转xi、xj两个城市之间的城市次序来得到S的邻居。设：S=(x1,x2,xi-1,xi,xi+1,xj-1,xj,xj+1,xn)则：S=(x1,x2,xi-1,xi,xj-1,x j-2,xi+1,xj,xj+1,xn)

5、x1x2xnxj+1xjxj-1xi-1xixi+1x1x2xnxj+1xjxj-1xi-1xixi+1l基本思想：在搜索过程中，始终向着离目标最接近的方向搜索。l目标可以是最大值，也可以是最小值。l在后面的介绍中，如果没有特殊说明，均假定是最小值。局部搜索算法（Local Search）1，随机的选择一个初始的可能解x0D，xb=x0，P=N(xb)；2，如果不满足结束条件，则3，Begin4，选择P的一个子集P，xn为P中的最优解5，如果f(xn)f(xb),P=P xn =(a,d,c,b,e),(a,e,c,d,b),(a,b,d,c,e),(a,b,e,d,c),(a,b,c,e,d

6、)从P中选择一个元素，假设xn=(a,d,c,b,e),f(xn)=45,f(xn)f(xb),P=P xn =(a,e,c,d,b),(a,b,d,c,e),(a,b,e,d,c),(a,b,c,e,d)从P中选择一个元素，假设xn=(a,e,c,d,b),f(xn)=44,f(xn)f(xb),P=P xn =(a,b,d,c,e),(a,b,e,d,c),(a,b,c,e,d)从P中选择一个元素，假设xn=(a,b,d,c,e),f(xn)=44,f(xn)f(xb),P=P xn=(a,b,e,d,c),(a,b,c,e,d)从P中选择一个元素，假设xn=(a,b,e,d,c),f(x

7、n)=34,f(xn)f(xb),P=P xn =(a,d,e,b,c),(a,c,e,d,b),(a,b,d,e,c),(a,b,c,d,e),(a,b,e,c,d)从P中选择一个元素，假设xn=(a,d,e,b,c),f(xn)=39,f(xn)f(xb),P=P xn =(a,c,e,d,b),(a,b,d,e,c),(a,b,c,d,e),(a,b,e,c,d)从P中选择一个元素，假设xn=(a,c,e,d,b),f(xn)=38,f(xn)f(xb),P=P xn =(a,b,d,e,c),(a,b,c,d,e),(a,b,e,c,d)从P中选择一个元素，假设xn=(a,b,d,e,

8、c),f(xn)=38,f(xn)f(xb),P=P xn=(a,b,c,d,e),(a,b,e,c,d)从P中选择一个元素，假设xn=(a,b,c,d,e),f(xn)=38,f(xn)f(xb),P=P xn=(a,b,e,c,d)从P中选择一个元素，假设xn=(a,b,e,c,d),f(xn)=41,f(xn)f(xb),P=P xn=P等于空，算法结束，得到结果为xb=(a,b,e,d,c),f(xb)=34。l局部最优问题l l每次并不一定选择邻域内最优的点，而是依据一定的概率，从邻域内选择一个点，指标函数优的点，被选中的概率比较大，而指标函数差的点，被选中的概率比较小。l设求最大值

9、：)(max)()()(xNxjiijxfxfxPl当求最小值时：)(11)(1)(1()(1)(max)(maxmaxminixNxjiixPxNxPxPxPj局部搜索算法1（Local Search 1）1，随机的选择一个初始的可能解x0D，xb=x0，P=N(xb)2，如果不满足结束条件，则3，Begin4，对于所有的xP计算指标函数f(x)，并按照式（3）或者式（4）计算每一个点 x的概率5，依计算的概率值，从P中随机选择一个点 xn，xb xn，P=N(xb)，转26，End7，输出计算结果8，结束l步长问题 初始值搜索到的最优解l变步长初始值搜索到的最优解局部搜索算法2（Local

10、 Search 2）1，随机的选择一个初始的可能解x0D，xb=x0，确定一个初始步长计算P=N(xb)2，如果不满足结束条件，则3，Begin4，选择P的一个子集P，xn为P中的最优解5，如果f(xn)f(xb)，则xb xn6，按照某种策略改变步长，计算P=N(xb)，转27，否则P=P P，转2。8，End9，输出计算结果10，结束l起始点问题AB全局最大值局部最大值l随机的生成一些初始点，从每个初始点出发进行搜索，找到各自的最优解。再从这些最优解中选择一个最好的结果作为最终的结果。局部搜索算法3（Local Search 3）1，k=02，随机的选择一个初始的可能解x0D，xb=x0，

11、P=N(xb)3，如果不满足结束条件，则4，Begin5，选择P的一个子集P，xn为P中的最优解6，如果f(xn)f(xb)，则xb xn，P=N(xb)，转37，否则P=P P，转3。8，End9，k=k+110，如果k达到了指定的次数，则从k个结果中选 择一个最好的结果输出，否则转（2）11，输出结果12，结束 l以上几种解决方法可以结合在一起使用，比如第一、第二种方法的结合，就产生了我们将在后面介绍的模拟退火方法。皇后搜索算法（Queen Search）1，随机地将n个皇后分布在棋盘上，使得棋盘 的每行、每列只有一个皇后。2，计算皇后间的冲突数conflicts。3，如果冲突数confl

12、icts等于0，则转（6）4，对于棋盘上的任意两个皇后，交换他们的行 或者列，如果交换后的冲突数conflicts减少，则接受这种交换，更新冲突数conflicts，转3。5，如果陷入了局部极小，既交换了所有的皇后 后，冲突数仍然不能下降，则转1。6，输出结果7，结束。皇皇 后后 数数1005001000200050001000030000平均时间平均时间（秒）（秒）55122817090010000l是局部搜索算法的一种扩展l最早由Metropolis在1953年提出，Kirkpatrick等人在1983年成功地将模拟退火算法用于求解组合优化问题。l基本思想是借用金属的退化过程改进局部搜索算

13、法l溶解过程：随着温度的不断上升，粒子逐渐脱离开其平衡位置，变得越来越自由，直到达到固体的溶解温度，粒子排列从原来的有序状态变为完全的无序状态。l退火过程：随着温度的下降，粒子的热运动逐渐减弱，粒子逐渐停留在不同的状态，其排列也从无序向有序方向发展，直至到温度很低时，粒子重新以一定的结构排列。l l粒子不同的排列结构，对应着不同的能量水平。如果退火过程是缓慢进行的，也就是说，温度的下降如果非常缓慢的话，使得在每个温度下，粒子的排列都达到一种平衡态，则当温度趋于0（绝对温度）时，系统的能量将趋于最小值。l如果以粒子的排列或者相应的能量来表达固体所处的状态，在温度T下，固体所处的状态具有一定的随机

14、性。一方面，物理系统倾向于能量较低的状态，另一方面，热运动又妨碍了系统准确落入低能状态。l从状态i转换为状态j的准则：l如果E(j)E(i)，则状态转换被接受；l如果E(j)E(i)，则状态转移被接受的概率为：l其中E(i)、E(j)分别表示在状态i、j下的能量，T是温度，K0是波尔兹曼常数。KTjEiEe)()(l在给定的温度T下，当进行足够多次的状态转换后，系统将达到热平衡。此时系统处于某个状态i的概率由波尔兹曼（Boltzmann）分布给出：l （6）l其中 为归一化因子，S是所有可能状态的集合。TKTiEiZeTP)(SjKTjETeZ)(l考察一下式（6）随温度T的变化情况：同一温度

15、下，两个能量不同的状态高温下的情况低温下的情况 当温度下降时的情况l在给定的温度T下，设有i、j两个状态，E(i)E(j)：l即在任何温度T下，系统处于能量低的状态的概率大于处于能量高的状态的概率。KTiEjEKTiETKTiEKTjEKTiETTKTjETKTiEjieeZeeeZZeZeTPTP)()()()()()()()(1111)()(由于E(i)E(j)，所以该项小于1 l当温度趋于无穷时：其中|S|表示系统所有可能的状态数。l当温度很高时，系统处于各个状态的概率基本相等，接近于平均值，与所处状态的能量几乎无关。SeeTPSjKTjEKTiETiT1lim)(lim)()(l当温度

16、趋于0时：l设Sm表示系统最小能量状态的集合，Em是系统的最小能量。上式分子、分母同乘以 SjKTjEKTiETiTeeTP)()(00lim)(limKTmEemmmSjKTEjEKTEiETSjKTEjESjKTEjEKTEiETSjKTEjEKTEiETiTSiSiSeeeeeeeTPmmmmmmmmmm如果如果01limlimlim)(lim)()(0)()()(0)()(00l当温度趋近于0时，系统以等概率趋近于几个能量最小的状态之一，而系统处于其他状态的概率为0。以概率1达到能量最小的状态。l当温度上升或下降时：TiSjjiSjTKTjEiSjKTjETiiTTiiTKTiETTK

17、TiETKTiEiEiEKTTPTPjEiEKTTPZejEiEKTTPejEKTZTPTPKTiETZZTPTPKTiETZeZZeKTiEZeTTTP)()()()()()()()()()(1)()()()()()(1)()(22)(2)(222)(2)(2)(l系统落入低能量状态的概率随着温度的下降单调上升，而系统落入高能量状态的概率随着温度的下降单调下降。TTiEiEEiETTP)(0)(0)(如果如果l在高温下，系统基本处于无序的状态，基本以等概率落入各个状态。在给定的温度下，系统落入低能量状态的概率大于系统落入高能量状态的概率，这样在同一温度下，如果系统交换的足够充分，则系统会趋向

18、于落入较低能量的状态。随着温度的缓慢下降，系统落入低能量状态的概率逐步增加，而落入高能量状态的概率逐步减少，使得系统各状态能量的期望值随温度的下降单调下降，而只有那些能量小于期望值的状态，其概率才随温度下降增加，其他状态均随温度下降而下降。因此，随着能量期望值的逐步下降，能量低于期望值的状态逐步减少，当温度趋于0时，只剩下那些具有最小能量的状态，系统处于其他状态的概率趋近于0。因此最终系统将以概率1处于具有最小能量的一个状态。l（1）初始温度必须足够高；l（2）在每个温度下，状态的交换必须足够充分；l（3）温度T的下降必须足够缓慢。固体退火过程固体退火过程组合优化问题组合优化问题物理系统中的一

19、个状态物理系统中的一个状态组合优化问题的解组合优化问题的解状态的能量状态的能量解的指标函数解的指标函数能量最低状态能量最低状态最优解最优解温度温度控制参数控制参数1，随机选择一个解i，k=0，t0=Tmax（初始温度），计算指标函数f(i)。2，如果满足结束条件，则转（15）。3，Begin4，如果在该温度内达到了平衡条件，则转（13）。5，Begin6，从i的邻域N(i)中随机选择一个解j。7，计算指标函数f(j)。8，如果f(j)j)Random(0,1)，则i=j，f(i)=f(j)。11，转（4）12，End13，tk+1=Drop(tk)，k=k+1。14，End15，输出结果。16

20、，结束。tifjftejiP)()()(l初始温度的选取l内循环的结束条件，即每个温度状态交换何时结束l外循环的结束条件，即温度下降到什么时候结束l温度的下降方法l在模拟退火过程中，给定温度下状态（解）的转移可以看作是一个马尔可夫链。对于任意两个状态i和j，我们用Pt(i,j)表示温度t下，从状态i转移到状态j的一步转移概率，则有：l其中：Gt(i,j)是产生概率，表示从状态i产生状态j的概率。At(i,j)是接受概率，表示在状态i产生状态j后，接受状态j的概率。jiliAliGjijiAjiGjiPiliNlttttt如果如果),(),(),(1),(),(),(定理1满足条件的Gt(i,j

21、)、At(i,j)举例：其他如果0)(1),(iNjNjiGt其他如果tifjfteifjfjiA)()()()(1),(说明：条件2的后半部分除外，该条件与具体的问题有关。定理2：在定理1的条件下，如果对于任意两个状态 有：则有：Sji,0),(lim)()(0jiAjfiftt其他如果01)(lim0mmttSiSiP定理3（放宽了定理1的条件）Gt(i,j)、At(i,j)满足定理1中除条件（2）以外的所有其他条件，并且：1，对于任意两个状态i、j，它们相互为邻居或者相互都不为邻居；2，对于任意i，Gt(i,j)满足：3，状态空间S对于邻域是连通的；则与模拟退火算法相伴的时齐马尔可夫链存

22、在平稳分布，其分布概率为：其他如果0)()(1),(iNjiNjiGtSijiAjNiiAiNiPSjmtmtt,),()(),()()(l（1）初始温度t0；l（2）温度t的衰减函数，即温度的下降 方法；l（3）算法的终止准则，用终止温度tf或 者终止条件给出；l（4）每个温度t下的马尔可夫链长度Lk。l一个合适的初始温度，应保证平稳分布中每一个状态的概率基本相等，也就是接受概率P0近似等于1。在Metropolis准则下，即要求：10),(tjifel如果我们给定一个比较大的接受概率P0，则：100ln),(Pjiftl其中，可以有以下估计方式：),(jif)(min)(max),(ifi

23、fjifSiSi2,)()(),(SjfifjifSji)1()(),(10SiSfiSfjifSil假设在t0下随机的生成一个状态序列，分别用m1和m2表示指标函数下降的状态数和指标函数上升的状态数，表示状态增加的平均值。则m2个状态中，被接受的个数为：),(jif0),(2tjifeml所以平均接受率为：l求解有：21),(2100mmemmPtjif0102201ln),(PmPmmjiftl模拟固体的升温过程：（1）给定一个希望的初始接受概率P0，给定一个较低的初始温度t0，比如t01；（2）随机的产生一个状态序列，并计算该序列的接收率：如果接收率大于给定的初始接受概率P0，则转（4）

24、；（3）提高温度，更新t0，转（2）；（4）结束。产生的状态总数接受的状态数l等比例下降 kktt110l等值下降 tttkk1l由定理1我们知道，在一定的条件下，与模拟退火算法相伴的时齐马尔可夫链存在平稳分布。如果温度每次下降的幅度比较小的话，则相邻温度下的平稳分布应该变化不大，也就是说，对于任意一个状态i，相邻温度下的平稳分布应满足：1)()(111iPiPkkttl一个充分条件是：11)()(kmkmtfiftfifee1)()()1(kmkmtfiftfifeel两边取对数，并整理得：l用 代替 可得温度的衰减函数：mkkkfifttt)()1ln(11kt3mfif)(ktkkktt

25、t3)1ln(11l固定长度方法 l在每一个温度下，都使用相同的Lk。Lk的选取与具体的问题相关，一般与邻域的大小直接关联，通常选择为问题规模n的一个多项式函数。l基于接受率的停止准则：规定一个接受次数R，在某一温度下，只有被接受的状态数达到R时，在该温度下的迭代才停止，转入下一个温度。规定一个状态接受率R，R等于该温度下接受的状态数除以总生成的总状态数。如果接受率达到了R，则停止该温度下的迭代，转入下一个温度。在迭代的过程中，若干相邻的状态称为“一代”，如果相邻两代的解的指标函数差值小于规定的值的话，则停止该温度下的迭代。l零度法 设定一个正常数e，当时tke时，算法结束。l循环总控制法 给

26、定一个指定的温度下降次数K，当温度的迭代次数达到K次时，则算法停止。l无变化控制法 如果在相邻的n个温度中，得到的解的指标函数值无任何变化，则说明算法已经收敛。l接受概率控制法 给定一个小的概率值p，如果在当前温度下除了局部最优状态外，其他状态的接受概率小于p值，则算法结束。l领域平均概率控制法 设大小为N的一个领域，在邻域内一个状态被接受的平均概率为1/N。设f0、f1为该领域中的局部最优值和局部次最优值。则次最优解是除了局部最优解以外接受概率最大的，其接受概率为：tffe01l如果该概率值小于平均值1/N时，即：可以认为从局部最优解跳出的可能性已经很小了，因此可以终止算法。此时的终止温度t

27、f为：Netff101Nfftfln21l相对误差估计法 设温度t时指标函数的期望值为：则当终止温度1时，由泰勒展开近似有：SitiPiftf)()()(fttfmfdttfdtftf)()(l由于：l所以可用下式估计当前解与最优解之间的误差：mtftf)(lim0fttfmfdttfdtftf)()(l或者使用相对于 的相对误差：fttfdttfdft)()()(f)()()(2ftdttfdftftttfffl实际计算时：l其中：)()()(022tfttftffffniiXfntf1)(1)(niiXfntf122)(1)(l解的表示：n个城市的任何一种排列均是问题的一个可能解，表示为:

28、l指标函数(能量函数)其中),.,(1nniniidf111),.,(11nl新解的产生 采用第一节介绍的两个城市间的逆序交换方式得到问题的一个新解。设当前解是 ，被选中要逆序交换的城市是第u和第v个到访的城市，uv。则逆序排列u和v之间的城市，得到问题的新解为：则两个路径的距离差为：),.,(1n),.,.,.,(1111nvvuvu)()(1111vvuuvuvuddddfl新解的接受准则 其他如果tftefA01l初始参数的确定l康立山等人的方法：初始温度t0=280；在 每 个 温 度 下 采 用 固 定 的 迭 代 次 数，Lk=100n，n为城市数；温度的衰减系数0.92，即tk+

29、1=0.92tk；算法的停止准则为：当相邻两个温度得到的解无任何变化时算法停止。lNirwan Ansari和Edwin Hou的方法：初始温度t0是这样确定的：从t0=1出发，并以t0=1.05t0对t0进行更新，直到接受概率大于等于0.9时为止，此时得到的温度为初始温度；在每个温度下采用固定的迭代次数，Lk=10n，n为城市数；温度的衰减系数0.95，即tk+1=0.95tk；路径长度出现次数平均转移次数路径最优2.6919063952BCADEFGHIJ次优2.752464056BCADEGFHIJ第三2.769104053DEFGHIJCBA最差2.89854497ABCDEFHIJG 路径长度出现次数平均转移次数路径最优24.387928740ACLBIQFTMEPRGSOJHDKN次优24.621678638ADCLBIQFTMEPRGSOJHKN第三25.17399902ANKDHIOJSGRPEMTFQBLC最差25.5015794AQFTMEPRGSJOIBLCDHKN