改进模拟退火算法及其在装填问题中应用.docVIP

改进模拟退火算法及其在装填问题中应用 　　摘要：NP难度问题一直是计算机科学研究的一个重要问题，具有很高的理论和实用价值。这篇文章主要研究利用模拟退火算法解决具有NP难度的装填问题，它的求解目标是寻求多个圆在一个矩形内的优良布局，使得这些圆两两互不嵌入地放置。通过将模拟退火算法与梯度法相结合，并且融入一些启发式的格局更新策略，提出了改进的模拟退火算法。计算结果显示，该算法对于解决圆形装填问题具有很高的效率。 　　关键词：模拟退火算法；装填问题；NP难度；梯度法 　　中图分类号：TP18 文献标识码：A 文章编号：1009-3044（2016）06-0181-03 　　1 概述 　　装填问题研究的是一组小物体在大区域内互相不重叠的布局方式，要求尽可能地提高空间或物体的利用率，并达到某些最优指标[1-3]，从而减少浪费。在物体堆放、运输以及原材料下料等领域有着广泛的应用，有效地解决这个问题可以增加资源的利用率，从而带来不可估量的经济效益。由于装填问题通常是NP难度的问题，所以目前学术界没有算法能保证在有效时间内求出其精确解，顺着应用行业的迫切要求，出现了很多有效的搜索策略来解决这一问题。 　　本文在模拟退火算法的基础上，提出了一个改进的模拟退火算法，它是在模拟退火算法能够有效地逃离局部最小值陷阱策略的基础上，又结合了梯度算法和空白点放置策略，有效地提高了算法的效率。 　　2 问题的描述 　　2.1 NP问题和NP完全问题介绍 　　首先我们来看下NP问题的定义，所谓N也就是Non-Deterministic，即不确定性，P则是多项式算法问题，所以NP就叫做非确定性多项式算法问题，有些多项式问题没有一套完整的求解公式，只能使用确定性的猜测的方法和验证的方法来求解，对于这一类求解的问题我们都可以叫它NP问题，如果来求解NP问题，我们一般使用一个合理的算法来对NP问题进行猜测和验证 　　如果以上的NP问题可以用一个多项式算法得到这个NP问题的最优解，就叫做NP完全问题的解，但是这个难度是很大的，也是目前世界数学界的难题。本文研究的是对于NP问题的相对最优解，采用改进模拟退火算法来对布局进行优化，从而得到优质的计算结果。 　　2.2 装填问题介绍 　　装填问题所探讨的是一组小物体放在一个大区域内并且让它们互相不重叠的布局方案，要求尽可能地提高空间和物体的利用率，并达到某些最优指标，从而减少浪费。下面来看下本文研究的装填问题的说明。 　　给定一个矩形区域和N个不同规格的圆，圆形装填问题就是问能否将这些圆互不重叠地放到矩形容器中，此问题更形式化的描述如下： 　　将二维笛卡尔坐标的原点取在矩形区域的左下角，如图1所示： 　　记矩形容器的长宽分别为L与W，圆i的圆心坐标为（xi， yi），半径为Ri，圆j的圆心坐标为（xj，yj） ， 半径为Rj，问是否存在2n个实数x1，y1，…，xn，yn满足以下两个约束条件： 　　[Ri≤xi≤L-Ri， Ri≤yi≤W-Ri（xi-xj）2+（yi-yj）2≥Ri+Rj] 　　如果存在，则给出一组合乎条件的解（x1，y1，…，xn，yn），这里i， j=1，2，…，N， and i≠j。 　　3 问题的转化 　　根据上面的问题描述，我们按照拟物方法[10， 11]，假设这N个圆为有弹性的光滑的圆形物体，将它们随机的放入这个矩形容器内，若物体间两两无挤压，则得到问题的解，否则根据弹性力学，受挤压的圆间将产生积压弹性力，并在此力的作用下产生一系列恢复原本形状的运动，直到物体间没有挤压则得到此问题的解，这样通过拟物的方法可以转化为一个势能函数的优化问题，根据弹性力学原理就可以得出势能函数。以下是势能函数的求得过程： 　　第一种情况，第i个圆和第j个圆相互嵌入时，这时嵌入深度的是它们半径之和减去圆心的距离；当两圆不嵌入时，则[dij]为0，见图2所示： 　　因此圆i与圆j的挤压深度为： 　　5 改进的模拟退火算法 　　5.1 挑圆策略 　　当模拟退火算法陷入局部最小值陷阱后，我们需要挑出一个圆重新放置来使算法跳出陷阱从而继续寻求最优解，那么选择挑出哪一个圆将是一个待解决的问题。我们借鉴黄文奇等人的拟人策略来解决这一问题，在热闹拥挤的公交车中，受挤压最甚者总能设法改变自己的位置，而处境宽松者往往也会让出一部分空间给予别人。所以当陷入局部最小值陷阱时，我们可以挑出受挤压最严重的圆饼，将它重新放置到矩形容器的某个地方去，也可以跳出那个挤压最宽松的圆放到矩形容器的某个地方去，挤压的严重程度我们可以用圆的势能与其半径比来表示，我们可以将挑出势能半径比最大的圆的策略称之为压力解除策略，将挑出势能半径比最小的圆的策略称之为资源让与策略。此两个策略统称为挑圆策略。 　　5.2 空