一种改进模拟退火组合优化法.pdfVIP

信息与控制 1990年 第 3期 一 种改进的模拟退火组合优化法 徐 雷 HE京大学数学系) 摘要 本文将模拟退、火组台优化珐的搜索控错过程与最优■的更新过程加 以分离．提出了一种有教的改进方珐 ． 一 … 它得到的解总是优于(至少不劣于)原有模拟退火珐的解．尤其是在固有关参数选辑不台适使原模撤退火法获得的解较 差时 改进{击仍能得到较椭意的鼻．另外．为碱少计算量，改进珐迁撮供 了一种判断M~~op01i$抽样停止和退火过程 结束的有技方案．最后．培出了改进珐在两十模式讽别典型问题上的应用．以及计算机实验结果． 美键谓：模拯退火．属性匿匹配．商傲优化．囊类分析 1 引言 近年来 ，S．Kirkpatrick等人 提出的用于求解组台优化问题 的模拟退火方法，引起了不 同领域中的许多人的注意和兴趣。该法将组合优化问题与统计力学中的热平衡问题类比，另辟 了 一 条求解组合优化问题的新途径，它通过模拟退火过程 ，可找到全局{或近似的)最优解。目前 ， 该法已被成功地应用于 TSP问题、VLSI设计和计算机设计中的体设计、结构宏配置及布线等问 题 。例如，用该方法可得到规模达 10000个城市的TSP问题的近似最优解 ．而用现有的其它方 法能谚求解的问题规模却只约有 300个城市。事实上，模拟退火方法是一种非常简单而有效的通 用随机优化方法，不仅在离散组合优化方面有广泛的应用前景，而且可同样有效地应用于求解连 续多变量擐优化 问题 ” 。 本文通过对模拟退火法的算法过程加 以剖析，提 出将算法的搜索控制过程与最优解的更新过 程加 以分离，给出了一种改进的模拟退火法。它得到的解总是优于(至少不劣于)原模拟退火法的 解，可克服原模拟退火法有时最终结果还不如中间结果的缺点 ．尤其是在因参数选择不合适 使原方法获得的解较差时，改进法仍能得到较满意的解。另外，改进法还提供了一种判断Me． tJopolis抽样停止和退火过程结束的有效方案，可使计算量大大减少，本文还将改进法应用于求 解两个模式识别 的典型问题，并以计算机实验结果表明了改进法的有效性． 本文只讨论用于离散问题的改进法。但本文的思路和主要做法，同样适用于连续多变量优化 问题 ”。 2 模拟退火法简介 设 = ． ，…，s 为所有可能的组台(珥下简称状态)所构成的集，C：S~R为非负代价 函数，即C(ss)1反映取状态 为解的代价，则组合优化问题可形式地表述为寻找 3．∈S，使 C(s．)=minC(s sf∈ (】) 模拟退火法求解组合优化 问题的基本思想是：把每种组台状态 ，看成某一物质体系的微观状 态，而 c()看成该物质体系在状态 下的内能，并用控制参数 r类比温度。让 r从一个足够高 的值慢慢下降，对于每个 用 Metropolis抽样法在计算机上模拟该体系在此 r下 的热平衡 态．即对 当前状态 做随机扰动 产生一个新状态 ，，计算增量 AC’=c(’一c()，并 以概率 exp卜AC／T3接受 作 为新 的当前状态。当重复地如此随机扰动足够 多次数后，各状态 出现 为当前状态 的概率将服从 Boltzmann分布，即 厂=zf P—co Z(T1= 1／ P ” (2) 若 7’下降的足够慢，且 一0，由式 (2)可见，当前状态将是具有最小 Cfs)的状态。 将上述思想写成算法形式，将包括Metropolis抽样算法和退火过程实现算法两部分． 收到本文的时间是1988年 l1月I5口．收到修改稿的时间是l990年2月25日 信息与控翩 1990年 第 3期 (1)退火过程实现算法(AnnealingProoedure，简称 AP算法) 步一： 任选一初始状态 so作为初始的当前解 s(O)=os，设初始温度为 ，并令 f=O 步二： 令 = ，