模拟退火算法 第一节优秀课件.ppt

《模拟退火算法 第一节优秀课件.ppt》由会员分享，可在线阅读，更多相关《模拟退火算法 第一节优秀课件.ppt（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、模拟退火算法 第一节第1 页，本讲稿共22 页一固体退火过程退火是一种物理过程，固体退火是先将固体加热至熔化，再徐徐冷却使之凝固成规整晶体的热力学过程退火过程中，系统在每一温度下达到平衡态，系统状态的分布满足一定的概率分布，即在温度 T，系统达到平衡态后，分子停留在状态 r 满足波兹曼(Boltzmann)概率分布2.1模拟退火算法及模型第2 页，本讲稿共22 页其中，E(r)为状态 r 的能量，kB 0为波兹曼常数，为分子能量的一个随机变量，称为波兹曼因子Z(T)为概率分布的标准化因子，先研究由(2.1)确定的函数随 T 变化的趋势选定两个能量 E1 E2，在同一个温度 T，有D 为状态空间

2、第3 页，本讲稿共22 页在同一个温度，(2.2)表示分子停留在能量小状态的概率比停留在能量大状态的概率要大当温度相当高时，(2.1)的概率分布使得每个状态的概率基本相同，接近平均值1D，D 为状态空间 D中状态的个数此时，具有最低能量状态的波兹曼概率接近并超出平均值1D 第4 页，本讲稿共22 页当 rmin 是 D 中具有最低能量的状态时，得由第5 页，本讲稿共22 页所以，关于温度 T 是单调下降的又有其中，D0是具有最低能量的状态集合，第6 页，本讲稿共22 页因此得到，当 T 趋向于 0 时，当温度趋向于 0时，(2.1)决定的概率渐近由此可以得到，在温度趋向于 0时，分子停留在最低

3、能量状态的概率趋向1综合上面的讨论，分子在最低能量状态的概率变化趋势由图(a)表示第7 页，本讲稿共22 页对于非能量最小的状态，由(2.2)和分子在能量最小状态的概率是单调减小的事实，在温度较高时，分子在这些状态的概率在 附近，依赖于状态的不同，使(2.1)决定的概率在(0,t)是单调升的；再由(2.4)可知，当温度趋于 0时，(2.1)定义的概率趋于 0概率变化曲线见图(b)可能超过由(2.3)和(2.4)可知存在一个温度t，第8 页，本讲稿共22 页从上面的讨论得到，在温度很低时，能量越低的状态的概率值越高，在极限状况，只有能量最低的点概率不为即有1.系统在 T 平衡时，系统状态的概率分

4、布趋于(2.1)式，第9 页，本讲稿共22 页0.002 0.016 0.117 0.865 t=0.50.181 0.221 0.269 0.325 t=50.232 0.243 0.256 0.269 t=20例2.1简化概率分布(2.1)为其中q(t)为标准化因子设共有四个能量点x=1,2,3,4,在此观察 t=20,5,0.5,三个温度点概率分布变化第10 页，本讲稿共22 页二 Metropolis 准则固体在恒定温度下达到热平衡的过程可以进行模拟1953年，Metropolis 等提出重要性采样法他们用下述方法产生固体的状态序列：先给定以粒子相对位置表征的初始状态 i，作为固体的当

5、前状态，该状态的能量是 Ei 然后用摄动装置使随机选取的某个粒子的位移随机地产生一微小变化，得到一个新状态 j，新状态的能量是Ej 如Ej Ei，则考虑热运动的影响，该新状态是否重要状态，要依据固体处于该状态的几率来第1 1 页，本讲稿共22 页判断由(2.1)知，固体处于 i 和 j 的概率的比值等于相应Boltzmann 因子的比值，即r 是一个小于1的数用随机数发生器产生一个0,1)区间的随机数，若r，则新状态 j 作为重要状态，否则舍去若新状态 j 是重要状态，就以 j 取代 i 成为当前状态，否则仍以 i 为当前状态，再重复以上新状态产生过程在大量固体状态的变换后，系统趋于能量较低的

6、平衡状态，固体状态的概率分布趋于(2.1)式的Boltzmann 概率分布第12 页，本讲稿共22 页由()式可知，高温下可接受与当前状态能差较大的新状态为重要状态而在低温下只能接受与当前状态能差较小的新状态为重要状态.这与不同温度下热运动的影响完全一致在温度趋与零时，就不能接受任一 Ej Ei 的新状态 j 了上述接受新状态的准则称为Metropolis 准则，相应的算法称为Metropolis 算法这种算法的计算量显著减少第13 页，本讲稿共22 页三模拟退火算法对固体退火过程的研究给人们以新的启示1982年，Kirkpatrick 等首先意识到固体退火过程与组合优化问题之间存在的类似性，

7、Metropolis等对固体在恒定温度下达到热平衡的模拟也给他们以启迪：应该把Metropolis 准则引入到过程中来最终他们得到一种对Metropolis 算法进行迭代的组合优化算法，这种算法模拟固体退火过程，称之为模拟退火算法第14 页，本讲稿共22 页我们可以将组合优化问题同金属物体退火进行类比：组合优化问题金属物体假设算法用以解决如下组合优化问题：解费用(目标)函数最优解状态能量能量最低的状态第15 页，本讲稿共22 页模拟退火算法Step1 任选一个初始解 x0；xi:=x0；k:=0；Step2 若在该温度达到内循环条件，则到step3;Step3 tk+1:=d(tk)；k:=k

8、+1；若满足停止条件，终t0:=tmax；(初始温度)；否则，从邻域 N(xi)中随机选一xj，计算fij=f(xj)f(xi)；若fij 0，则 xi:=xj；否则若exp(fij/tk)random(0,1)时，则 xi:=xj；重复step2；止计算；否则，回到step2 产生一个0与1之间的一个随机数第16 页，本讲稿共22 页(1)随算法进程递减其值的控制参数 t 担当固体退火过程中的温度 T 的角色(2)对 t 的每一取值，算法持续进行“产生新解判断接受舍弃”的迭代过程就对应着固体在某一恒定温度下趋于热平衡的过程也就是执行了一次Metropolis 算法模拟退火算法从某个初始解出发

9、，经过大量解的变换后，可以求得给定控制参数值时组合优化问题的相对最优解然后减小t 的值，重复执行Metropolis 算法，就可以在 t 趋于零时，最终求得整体最优解(3)由于退火必须“徐徐”降温，t 也必须缓慢衰减，才能确保模拟退火算法最终趋于组合优化问题的整体最优解集第17 页，本讲稿共22 页(4)模拟退火算法依据Metropolis 准则接受新解，因此除接受优化解外，还在一个限定范围内接受恶化解，这正是模拟退火算法与局部搜索算法的本质区别所在开始时 t 值大，可能接受较差的恶化解；随着 t 的减小，只能接受较好的恶化解；最后在 t 值趋于零时，就不再接受任何恶化解了这就使算法既可以从局

10、部最优的陷阱中跳出，更有可能求得组合优化问题的整体最优解；又不失简单性和通用性因此对大多数组合优化问题而言，模拟退火算法要优于局部搜索算法第18 页，本讲稿共22 页模拟退火算法的数学模型可以描述为：在给定邻域结构后，模拟退火过程是从一个状态到另一个状态不断地随机游动我们可用马尔可夫链描述这一过程当温度 t 为一确定值时，两个状态的转移概率定义为：Gij(t)称为从 i 到 j 的产生概率，Gij(t)表示在状态i 时，j 状态被选取的概率比较容易理解的是 j 是 i 的邻居，如果在邻域中等概率选取，则 j 被第19 页，本讲稿共22 页选中的概率为Aij(t)称为接受概率，Aij(t)表示在

11、状态 i 产生 j 后，接受 j 的概率，如在模拟退火算法中接受的概率是第20 页，本讲稿共22 页(2.5),(2.6),(2.7)为模拟退火算法的主要模型模拟退火算法主要可以分为两类第一类为齐次算法，在(2.5)中对每一个固定的 t，计算对应的马尔可夫链变化，直至达到一个稳定状态，然后再使温度下降第二类是非齐次算法，由一个马尔可夫链组成，要求在两个相邻的转移中，温度 t 是下降的第21 页，本讲稿共22 页描述模拟退火算法过程的马尔可夫链应满足下列条件：(1)可达性无论起点如何，任何一个状态都可以到达(2)渐进不依赖起点(3)分布稳定性包含两个内容：一是当温度不变时，其马尔可夫链的极限分布存在；二是当温度渐近 0时，其马尔可夫链也有极限分布(4)收敛到最优解当温度渐近 0时，最优状态的极限分布和为1第22 页，本讲稿共22 页