美文网首页
模拟退火算法 Simulated Annealing

模拟退火算法 Simulated Annealing

作者: Cloud_J | 来源:发表于2019-05-25 16:31 被阅读0次

算法思想

自然中的退火(不理解问题不大)

模拟退火算法的思想受启发于自然界中固体由高温到低温的过程中其内部分子状态及内部能量的变化规律。

退火指物体逐渐降温冷却的物理现象。温度越低,物体的能量越低,在结晶状态是系统的能量状态到达最低。在自然中,缓慢降温(退火)可以导致结晶,而与之相对的快速降温(淬火)会导致不是最低能态的非晶体形态。

退火的过程可以表示为下图,左边为最初的非晶态状态;经过升温,系统能量增大后到达中间的状态;再缓慢降温到达晶体态,此时能量最小

自然界退火

模拟退火

我们用一个搜索函数最优解来直观表示:C为函数的全局最优解,在只采用贪心策略的情况下,如果从A点开始搜索,最终得到的解为B点,然而这只是一个局部的较好解。

搜索最优解

为了避免陷入局部的最优解,模拟退火算法在搜索过程中加入了一个随机因素,会以一定的概率接收一个比当前解较差的解,因此就有可能越过B与C之间的高峰,到达全局最优解

在这里,可以将解(横坐标的值)理解为固体的状态,函数值理解为系统的内能。算法以固体所处的温度T为控制参数,随着T的下降使固体内能(目标函数值)也逐渐下降,直至趋于全局最小。

根据Metropolis准则,在温度为T时,接受能量从E(old)到E(new)的概率为P:
P=\begin{cases} e^{\frac{E(x_{new}) - E(x_{old})}{kT}},&E(x_{new}) > E(x_{old})\\ 1,&E(x_{new}) < E(x_{old}) \end{cases}
在特定温度下,经过充分转换,材料达到热平衡。这时材料处于状态 i 的概率为
P(x=i) = \dfrac{e^{-\frac{E(i)}{KT}}}{\sum_{j\in S}\;{e^{-\frac{E(j)}{KT}}}}
其中x表示材料当前的状态,S表示材料的状态集合,k 为玻尔兹曼常数。根据上式可以得到以下结论

  • T 趋向无穷大时,任意状态的 P 都相同
  • T 下降到极低时,材料进入最低能量的概率为1

因此,如果我们运用退火思想放在优化问题上,在降温过程中问题的解进行充分地“热交换”,即进行充分地重新排列,同样可以帮助我们寻找最优解,理论上也会具有达到全局最优解的性能!

算法步骤

  1. 设置初始温度T_0,令当前温度T_{current}=T_0,并设置一个终止温度T_f和一个随机的初始解x_0,令当前解x_i=x_0
  2. 使当前温度发生变化,T_{current} = \alpha T_{current}\alpha的值取0.5~0.99之间
  3. 计算当前状态x_i(解)下的内能E(x_i)(函数值), 根据当前状态进行扰动,产生一个新的解x_j,计算\Delta E
  4. \Delta E<0,则接受该解,否则以概率P接受该解
  5. 在当前温度下,循环k次执行3和4
  6. 比较T_{current}T_f的大小,若T_{current}<T_f,则终止,否则执行2

可以看出,算法实际上是两层循环嵌套,外层循环控制温度,内层循环来进行扰动产生新解。

随着温度逐渐降低,算法最终由可能收敛到全局最优,这里说有可能的原因是因为,在温度很低时,虽然从地内能状态跳到高内能状态的可能性不大,但是也有可能发生。

应用场景

在已知的某个定义域内求函数最优值的问题通常有以下三种情况,以求最小值问题为例(求最大值可以转换为求最小值)

  • 如果是离散取值,则可以用穷举法获得问题最优解
  • 如果取值连续,且函数是凸函数,则可以用梯度下降等方法获得最优解
  • 如果连续非凸,虽说根据已有的近似求解法能够找到问题解,可解是否是最优的还有待考量,很多时候若初始值选择的不好,非常容易陷入局部最优值
三种场景

模拟退火算法就是为了应对第三种情况而提出的

参考资料:
https://www.cnblogs.com/ranjiewen/p/6084052.html
https://blog.csdn.net/weixin_40562999/article/details/80853418

相关文章

网友评论

      本文标题:模拟退火算法 Simulated Annealing

      本文链接:https://www.haomeiwen.com/subject/edlzzqtx.html