[置顶] 粒子群优化算法解决场地资源布局问题-算法原型探讨

版权声明：本文为博主原创文章，未经博主允许不得转载。

　　本篇文章的命题来源于博主参加工作后的第一项任务，用智能优化算法解决工程实际中关于某施工场地相关资源布局的问题。由于博主的硕士毕业论文（详见第3章 基于PSO的已知静态路径规划方法）涉及到了粒子群优化算法，有一定经验，固对于本命题欲采用粒子群优化算法进行求解。
　　对于本命题而言，与博主的硕士毕业论文中解决的问题有很大差别，可以说是应用在不同的工程领域。当时花了两周的时间讨论并实现算法的原型，通过算法原型求解出的结果验证了粒子群优化算法还是能有效解决该类问题的。本篇博文要分享的就是关于场地资源布局问题，算法原型的设计、结果的分析和验证。博主对于粒子群优化算法的研究还尚浅，也就刚刚好能解决一些问题，完全可以将本篇文章当做粒子群优化算法的一个小算例，如有纰漏恳请各路大仙不吝赐教。

1. 场地资源布局问题描述

　　由于现实中的场地资源布局问题涉及到相关的业务流程和专业术语，为了方便描述，将该问题抽象化。如图１所示，场地共分为两大区域：总装区域和装配区域。总装区域中各个矩形块即为该问题中的建造资源，而各个建造资源所组成的多边形即为最终的产品，每个矩形块在总装区域中的位置即为该矩形块在该问题中的目标位置（已知条件）。各个矩形块为了能运输到总装区域进行总装，需要预先在装配区域完成一系列的装配作业，每个矩形块因各自的装配作业不同，在装配区域中装配的周期（已知条件，矩形块中ｔ参数）也不尽相同；每个矩形块在装配区域中的布局（待求解）需要满足的条件包括：矩形块不与其他矩形块相干涉，不与装配区域中障碍物（现实中用来运输相关资源的绿色通道）相干涉，布局不超越装配区域。而矩形块从装配区域运输到总装区域，运输的方向仅能为水平的X方向和垂直的Y方向。
　　需要解决的优化问题可描述为：为每个在装配区域装配的矩形块分配有效位置，使各个矩形块从装配区域运输到总装区域的X方向和Y方向的距离总和最小。
　　实际问题中，需要考虑的优化目标不仅仅为距离，仍需考虑每个矩形块重量，运输矩形块的Ｘ、Ｙ方向权重。而在原型算法中为了能更快速的实现及更直观的对结果进行判断，仅考虑了Ｘ、Ｙ方向运输距离。

图１　场地资源布局问题示意图

2. 算法原型设计

　　对于该场地资源布局问题，每个矩形块均需要求解X、Y方向的布局位置，即每个矩形块包含2个维度，若某批装配包含N个矩形块，则该批的求解维度为2×N。该问题无法通过简单的数值求解进行计算，所以需采用智能优化算法求解该问题，以期得到全局最优解。

2.1 粒子群优化算法描述

　　关于粒子群优化算法描述详见第3章 基于PSO的已知静态路径规划方法。

2.2 初始化

　　对于标准粒子群优化算法，种群的初始化都是随机产生的，以保证初始化种群的分布性，以较大概率覆盖整个布局解空间。然而盲目的随机选择位置会降低求解效率，另一方面，装配区域空间是有限的，当某批装配中矩形块过多、区域利用率较高的情况下，随机分布可能会导致区域中矩形块摆放过于杂乱无章，无法利用的区域过多，后进入的矩形块没有位置摆放。
　　所以为减少粒子群算法初始化的随机性和盲目性，引入正态分布赋值，以使各矩形块在装配区域的初始化位置更接近其在总装区域中的位置；说白了就是让１号矩形块在初始化Ｘ方向位置时，使其落在总装区域的总装Ｘ方向位置概率最大，越向两边概率越小。同时，为了保证初始化种群的分布性，以较大概率覆盖整个布局解空间，仍需保留随机方法生成初始化种群。
　　算法原型采用随机初始化和正态分布赋值初始化相结合的策略对种群进行初始化。若某批装配中矩形块过多，即矩形块占装配区域可用面积较大的情况下，随机初始化会导致区域浪费空间较多，进而导致粒子易陷入初始化。正态分布赋值初始化种群占种群比例为：每天矩形块占用面积和／每天装配区域可用面积和。采用该粒子初始化策略，即保证了区域利用率较高情况下的求解效率，同时也保证了初始种群的分布性。

2.3 参数设置

1）惯性权重ω的取值

　　算法原型采用线性递减公式取权值，即：
　　 ω=ωmax−(ωmax−ωmin)∙runrunmax
　　式中， ωmax 为最大惯性权重，取1.2； ωmin 为最小惯性权重，取0.4；run是当前迭代次数； runmax 为算法迭代的总次数。
　　这种迭代方法在初期时局部搜索能力较弱，后期全局搜索能力变弱。

2）学习因子 c1 和 c2 的取值

　　 c1 和 c2 通常为相等的常数，取值在0～4之间时，可搜索到较优的解。经过多次试验，本文 c1 和 c2 值均取为2.0。

2.4 位置优化算法

　　由装配区域与总装区域布局形式可知，较优的布局方案会紧贴距总装区域较近的一边，如图1中装配区域的下边；且朝向矩形块在总装区域中位置紧密靠拢。因此，为将布局结果与区域布局形式相结合，得到更符合实际的最优结果，加入位置优化算法策略（以图1中装配区域为例）：
　　Y轴坐标临位移动优化算法：将矩形块沿Y轴方向紧贴装配区域下边移动。当然算法原型由于目的在于验证，仅考虑了Ｙ轴方向的优化，后续发现与Ｘ方向进行优化（X轴坐标临位移动优化算法、X轴坐标正位搜索移动优化算法，实际算法中用到，算法原型未加入）配合，计算效果更好。

2.5 适应度函数

　　装配区域布局优化要求布局方案中所有矩形块移动距离之和最小，适应度函数如下：
　　 F=∑ni=1(LXi+LYi)

2.6 算法实现过程

　　本文提出的针对场地资源布局问题算法原型流程图见图2所示：

图2　算法原型流程图

3. 计算结果及探讨

　　采用上述描述的算法原型，对模拟数据进行计算。由于算法原型程序未进行优化，后续实践中发现仅对算法进行代码级别的优化，即可实现计算耗时的大幅缩短，因此计算结果中的计算耗时仅代表相对计算耗时值。

3.1 模拟数据计算-5个矩形块

　　模拟数据计算-5个矩形块数据及甘特图详见表1。表中时间单位为天，代表在起始和终止时间内矩形块都在装配区域；装配位置即为优化计算结果；计算的原点（X=0，Y=0）为总装区域左下角；通过甘特图可知只需第2天（t=2）的布局图即可反应整个优化计算结果。经过多次计算，可得出最优计算结果唯一，最佳适应度为124，见图3所示。从计算结果可知通过算法原型优化计算出的结果，符合约束条件，且得到了最优解，可以尝试更为复杂的模型进行测试。本次计算的粒子群优化算法参数见图3中所示。

表1　模拟数据计算-5个矩形块数据表及甘特图

图3　模拟数据计算-5个矩形块计算结果（t=2）

　　图４为该次计算出的适应度值随迭代次数变化情况，可知初始化的第一代适应度值为141，初始化没有过于盲目的随机赋值，即引入的正态分布赋值和位置优化算法显现成效。在第25代适应度值达到了130，与最优值的差距在5%以内，可以说收敛速度也较快。

图4　模拟数据计算-5个矩形块，计算适应度值随迭代次数变化情况

　　表2为使用不同的迭代次数、粒子数参数计算结果的汇总，从平均计算耗时和平均最佳适应度值来看，越大的迭代次数和粒子数，计算结果越优，而计算耗时越长。从相同条件下的计算耗时可知，该算法原型计算并不是很稳定，迭代次数100、粒子数50的计算数据中，甚至计算耗时能相差10倍以上，究其原因是因为初始化的随机性造成的。

表2 模拟数据计算-5个矩形块，计算结果汇总

3.2 模拟数据计算-9个矩形块

　　模拟数据计算-9个矩形块数据及甘特图详见表3。即在上节中5个矩形块的基础上再新增4个矩形块。通过甘特图可知需展示第2、4、5天的布局图反应整个优化计算结果。经过多次计算，亦得出最优计算结果唯一，最佳适应度为230，见图5所示。从计算结果可知通过算法原型优化计算出的结果，符合约束条件，且得到了最优解。本次计算的粒子群优化算法参数见图5中所示。

表3　模拟数据计算-9个矩形块数据表及甘特图

　　从图5中可以看到与上节计算结果不同的地方在4和5的位置，再通过图6、7可知，4和5的布局方式是为了在第5天，9号矩形块进入装配区域时能落位到最好的位置，即4和5为9进行了让位。

图5　模拟数据计算-9个矩形块计算结果（t=2）

图6　模拟数据计算-9个矩形块计算结果（t=4）

图7　模拟数据计算-9个矩形块计算结果（t=5）

　　然而你可能会质疑4和5的摆放位置，如图8中A即为求出的最优解，B为另外一种方案，看似结果相似；但经分析可知：若采用B方案9号（X方向增加4步），4号（X方向减少8步，Y方向增加6步；共减少2步），固B方案比A方案多耗费2步。

图8　9个矩形块方案对比

　　表4为使用不同的迭代次数、粒子数参数计算结果的汇总，总结同上节。

表4 模拟数据计算-9个矩形块，计算结果汇总

4. 算法优势

　　算法原型考虑的问题较简单，涉及因素较少，虽然能求解出简单计算模型的最优解，但未能充分提现出该算法的优势。结合后续的工作提前提出该算法解决该类问题的优势：
　　1）矩形块优先占用距离该矩形块在总装区位置最近的区域，靠近总装区域的矩形块排列紧密，场地利用率最高。
　　2）考虑全局最优，周转率、优先级较低的矩形块会为周转率、优先级较高的矩形块让位置。
　　3）充分有效利用装配区域中距离总装区域位置最近区域，最近区域变化最快，各矩形块尽量向总装位置靠拢，占地面积较大及占用场地周期长的矩形块尽量放置在距总装位置较远区域。

5. 仍需探讨的问题

　　通过算法原型的计算结果，可以初步判定该算法还是可以解决场地资源布局问题的，但算法原型设计过程中抽象并省略了很多约束条件及相关业务限制。在真正设计处理工程实际的算法中还有很多需要探讨的问题：
　　1) 原型算法中装配区域仅为一个，而实际中可能会存在两个及以上个装配区域，如何协调多个装配区域？
　　2) 原型算法初始化是随机的，可能每次计算结果都不一样，如何让每次计算结果相同，以便对结果进行重现？
　　3) 若因矩形块位置分布过于随机，且矩形块占装配区域面积较大，而导致后续入场的矩形块没有位置可放无法计算下去时，该如何处理？
　　4) 原型算法基于资源为矩形，若资源为不规则多边形如何处理？
　　5) 原型算法是基于矩形块周期不会改变的情况下的，但现实中该周期可能会经常改变，如何协调？
　　6) 若矩形块周期可前后浮动，如何与位置分布结合，计算出更优的方案？

==============文档信息===============
版权声明：本文为博主原创文章，未经博主允许不得转载
署名(BY) ：dkjkls（dkj卡洛斯）
文章出处：http://blog.csdn.net/dkjkls