混沌自适应水波算法在包装配送问题中的应用
彭 维
(重庆城市管理职业学院工商管理学院 重庆 401331)
摘 要 追求高效的启发式算法是包装配送问题的重要研究方向。对此，设计一种求解包装配送问题的混沌自适应水波算法(CAWWA)。根据包装配送问题特征设计水波算法(WWA)个体表达式；利用混沌系统生成WWA算法初始种群，并提出一种自适应调整的碎波系数，根据进化代数调整算法搜索范围；利用包装配送实例和6个国际算例进行仿真实验。结果表明，该算法能够适用于包装配送问题求解，且与GA算法、TS算法和ACO算法相比，其全局搜索能力更强、收敛速度更快。
关键词 包装配送 混沌 自适应 水波算法
0 引 言
包装配送问题是典型的组合优化问题[1]，通常可以描述为：从一个配送中心出发，安排若干车辆为客户配送包装，要求在满足交货时间、车辆载重、客户需求量等条件下，合理安排配送线路，使配送的路程、时间或者费用达到最优。近年来，随着电子商务深入发展，包装需求量和配送量井喷式增长，包装配送问题愈发受到重视，其求解算法也逐渐成为了学者们的研究热点。
分析大量文献可知，包装配送问题的求解算法可以分为精确算法和启发式算法，其中启发式算法占比高达85%左右。这主要是由于包装配送问题属于NP-hard难题，可行解数量会随着问题规模的增大而发生“组合爆炸”，计算开销也随之呈指数式增长[2-3]。精确算法在求解该类问题时，求解效率低、运行速度慢，往往无法取得令人满意的结果。相比之下，启发式算法具有快速寻优能力，可以在较短时间内求得大规模问题的较好解，因此成为了包装配送问题的主要求解算法。目前，应用于包装配送问题的主要启发式算法包括：遗传算法[4]、萤火虫算法[5]、禁忌搜索算法和其他算法[6]。
虽然已有大量优秀算法，但追求更高效的算法一直都是包装配送及其他组合优化问题的重要研究方向[7]。对此，本文引入一种新型启发式算法—水波算法(WWA)，并将其应用于包装配送问题中。WWA算法由郑宇军教授于2014年首次提出，具有参数较少、实现简单、计算开销小等优势[8]。在标准WWA算法基础上，本文对算法个体编码方式进行重新设计，并引入混沌机制和碎波系数自适应调整机制，提出了基于包装配送问题的CAWWA算法。实验表明，CAWWA算法能够很好地适用于包装配送问题的求解。
该词常用作为调侃自嘲自己要倒霉了、很衰的意思。《凉凉》原本是电视剧《三生三世十里桃花》片尾曲。“一首凉凉送给某某”源自于有些网络主播在直播玩游戏失败时，观众就会在弹幕上刷“一首凉凉送给主播”，从而走红网络。
1 包装配送问题的数学模型
假设配送中心安排K 辆车为L 个客户配送包装，车辆最大载重量为Q ，客户需求量为q i (i =1,2,…,L )，用i =0表示配送中心，c ij 表示客户i 到j 的配送成本(如距离、费用等)，定义决策变量为：
则包装配送问题的数学模型可表示为[9]：
(1)
(2)
(3)
(4)
(5)
(6)
x ijk =0或1 ∀i ,j ,k
(7)
y ik =0或1 ∀i ,k
(8)
其中:式(2)表示车辆最大载重约束；式(3)-式(5)表示每个客户都有且仅有一辆车配送；式(6)表示消除子回路；式(7)-式(8)表示变量取值范围。
2 标准水波算法
2.1 基本原理
WWA算法是模拟浅水波运动模型求解问题的过程而设计的一种新型启发式算法[10]。算法将问题求解空间看作为海床，每个可能解看作一个波高为h 和波长为λ 的水波。水波的适应度值与水波到海床的垂直距离成反比，即水波距海床越远，适应度值越小，内含能量也越小，且波长λ 越长，波高h 越低，这使得优质水波可以在最优解附近进行局部搜索，而劣质水波则能够跳出局部最优进行大范围搜索，进而求得问题最优解。WWA算法包括传播、折射和碎波3个算子。不同深度水域的波形如图1所示。
图1 不同深度水域的波形
2.2 传播算子
假设问题维度为D ，水波每一维位置为x (d )(1≤d ≤D )，每次迭代中水波通过传播公式对自身位置进行更新。
1.2.2.1 建立患者管理档案，与社区护士结对 出院前由专职小组人员与患者沟通，建立造口患者管理档案，内容包括患者姓名、地址、电话、疾病名称和造口情况等。在征得患者同意的情况下转诊社区医院，与社区护士结对，并提供病区咨询服务电话以及可联系的护理人员姓名等内容，方便护患联系和沟通。观察组和对照组出院后定期门诊随访。
x ′(d )=x (d )+rand (-1,1)·λ ·L (d )
(9)
式中:rand (-1,1)表示[-1,1]之间的均匀随机数；L (d )表示第d 维搜索空间的长度。位置更新后，如果适应度函数值f (x ′)>f (x )，则用x ′代替x ，水波波高h 重置为最大值h max；反之，保留x ，水波波高h 伴随能量的损耗而衰减。水波位置更新后，水波波长也发生改变，公式如下：
λ =λ ·α -(f(x )-f min+ε )/(f max-f min+ε )
(10)
式中：α 表示最大波长减少率；ε 表示为保证分母不为0而设置的极小正数；f max、f min分别表示本次迭代中最优水波和最差水波对应的适应度值。分析式(10)可知，波长与适应度值成反比。适应度值越大，水波波长越短，传播范围越窄；反之，水波波长越长，传播范围越大。
2.3 折射算子
步骤2 生成配送路径。根据M 2生成客户序列，并按照车辆载重约束生成配送路径。图2中客户序列为x =3-1-4-2。假设客户1到客户4的包装需求量分别为50、20、40和20，车辆最大载重量为70。则可得两条配送子路径，分别为x 1=0-3-1-0(0表示配送中心)，x 2=0-4-2-0。每条子路径的配送成本(距离、费用等)之和即为该方案的配送成本。
(11)
步骤2 根据第3.3节生成Num 个初始随机矩阵。
(12)
2.4 碎波算子
随着能量的不断增加，水波波峰变得越来越陡峭，最终破碎成大量独立的小波浪。WWA算法只对最优水波进行碎波操作，公式如下：
x ′(d )=x *(d )+N (0,1)·β ·L (d )
(13)
在WWA算法中，碎波算子主要用于对最优水波附近的空间进行搜索，而碎波系数β 是控制搜索范围的一个重要参数。β 越大，搜索范围越广，反之则搜索范围越小。传统WWA算法采用固定的碎波系数，明显无法满足算法前期和后期的不同搜索需求，无法达到全局搜索和局部搜索的平衡，不利于算法求解质量的提高。对此，本文提出了基于对数递减的碎波系数自适应调整机制。在算法前期，碎波系数较大，算法可以大范围地进行全局搜索，从而得到更优的水波。而在算法后期，碎波系数较小，算法能够更加精准地进行局部搜索，进而提高算法求解精度。公式如下：
3 混沌自适应水波算法
3.1 个体表达式构造
步骤1：生成随机矩阵。随机生成一个L ×L (L 表示客户数量)的客户矩阵M ，每个元素m ij (i 、j 为横坐标和纵坐标)的值介于0～1之间。
步骤1 设置算法参数，包括初始种群规模Num 、初始水波波高h max、波长λ 、波长减少率α 、最大碎波系数β max、最小碎波系数β min、对数调整因子γ 和算法最大迭代次数iterMax ，令当前迭代次数iter =0。
3.2 个体表达式解析
启发式算法对初始种群具有严重依赖性，如果初始种群个体在最优解邻域内，那么算法能够快速收敛到较好解甚至最优解。如果初始种群个体离最优解较远，则算法收敛速度较慢，且求解精度会大幅降低。混沌是一种广泛存在于社会和自然界中的非周期性运动，具有随机性、遍历性和规律性等特征，能够在一定范围内不重复地遍历所有状态[11-12]。对此，本文采用混沌机制来克服启发式算法初始种群分布不理想的缺点，以提高初始种群质量。本文采用混沌Logistic映射生成初始随机矩阵，公式如下：
水波传播过程中，有的水波会在浅水域聚集，最终在深水域发散反射。具体来说，当水波波高h 变为0时，水波将进行反射，对应位置根据下式进行更新：
随着科学技术发展和产品检验检测需要，检验方法和各项标准的迅速更新以及各类产品的日益多样化，要求检测机构各类人员不断学习，各检测机构应尽可能给各类人员提供适当的培训条件和机会[13]，有计划、有步骤、分层次地进行继续教育和学历教育，加快知识更新，提升学历层次，适应科学技术发展和产品检验检测需要[14]，以提高管理人员的管理水平和检测人员的专业技术能力和科学素养[15]，逐步建立并完善人员培训制度，从整体上提高员工队伍素质。
图2 个体表达式构成及解析
3.3 混沌初始化
步骤1 生成0-1矩阵。将概率矩阵M 1元素与波长λ 进行比较，若小于λ ，则相应位置为1，反之为0，得到0-1矩阵M 2，如图2所示(假设λ =0.2)。进行矩阵调整，确保M 2每行每列都只有一个1。调整原则为：如果第i 列只有一个1，则保留该1，如图2的第2列；如果第i 列不止一个1，假设所有1后面的列都还有1出现，则优先保留第一个1，假设有的1后面的列全是0，则优先保留该1，其他1变为0，如图2中的第1列和第3列。
住院患者服务中心由8人组成，其中3人为心理学专业。部门开展患者心理教育，服务对象覆盖从婴儿到80岁老人的全人群。刚刚降临人世的孩子需要人文关怀吗？很多人摇头，但每个孩子分明都是哲学家，他们的心理与成人不同，但不代表没有心理。
M n+1 =μM n (1-M n )
(14)
n =0,1,2,… 0≤μ ≤4
(15)
3.4 碎波系数自适应调整
式中:β 为碎波系数。碎波操作后，如果f (x ′)>f (x *)，则用x ′代替x *；否则，不更新。
β =β max-γ (β max-β min)×logiterMax iter
(16)
式中:β max、β min分别表示最大最小碎波系数；γ 表示对数调整因子；iter 表示算法当前迭代次数，iterMax 表示最大迭代次数。
3.5 CAWWA 算法求解包装配送问题
本文用
表示包装配送问题的适应度函数，算法求解过程如下：
步骤2：生成概率矩阵。对矩阵M 每一行进行归一化处理，得到概率矩阵M 1。比如，假设M 第一行为[0.20.30.80.9]，对应四个元素相加之和为0.2+0.3+0.8+0.9=2.2，各元素分别除以2.2，可得归一化处理结果为[0.090.140.360.41]，将其作为M 1的第一行。同理，可生成M 1其他数据。可知，M 1中每一行元素之和均为1，且每个元素表示该客户在客户序列中某个位置的概率。比如，假设M 1的第一行为[0.090.140.360.41]，则表示第一个客户排在第一位的概率为0.09，排在第二位的概率为0.14，排在第三位的概率为0.36，排在第四位的概率为0.41。
英国的培根说：“美犹如盛夏的水果，容易腐烂而难能保持。世上有许多美人，他们有过放荡的青春，却迎受着愧悔的晚年。因此，把美的形貌与美的德行结合起来吧。只有这样，美才会放射出真正的光辉。”这段话在学生时代就给我以深刻的印象。是的，在这个“看脸的社会”，一些人把爹妈给的容貌太当回事，傲娇自大，放弃努力，从而不学无术，活在虚荣和浮华里，而不懂得奋斗的意义和内涵的价值。“金玉其外，败絮其中”难免让人遗憾，是浮躁的价值观误导蒙蔽了他们。其实榜样不远，周恩来、宋庆龄就是高颜值，他们何曾因颜值而成就事业，他们最闪耀的是人格和才能。
式中：x *(d )表示当前种群的最优水波，N (μ ,δ 2)表示均值为μ 、标准差为δ 的高斯随机数。折射后的波长为：
步骤3 判断iter 是否大于或等于iterMax 。若是，算法停止，输出最优水波x *及f (x *)；反之，算法进入步骤5。
步骤4 按照第3.1节对个体矩阵进行归一化处理，并根据第3.2节解析个体并计算适应度值，记录最优水波x *及适应度值f (x *)。
步骤5 利用式(9)对每个水波x 进行传播，得到新的水波x ′，利用式(10)更新波长λ 。如果f (x ′)>f (x )，则用x ′代替x ，重置波高为h max；反之，保留x ，令波高h =h -1。
步骤6 判断水波波高h 是否等于0。若是，则利用式(11)-式(12)进行折射操作。
（5）区块写入账本。将对所有节点成功解出数学难题的广播答案进行验证，如果正确，它会将该区块纳入自己的账本中，每个节点同步进行；否则，将丢弃该区块。
步骤7 利用式(13)对最优水波x *进行碎波操作，得到新的水波。如果f (x ′)>f (x *)，则用x ′代替x *；反之，保留x *。
步骤8 令iter =iter +1，返回步骤3。
从上文分析可以看到，现在适用权利用尽理论的主要经济体，除了美国已经否定了该规则的区别适用理论之外，德国、日本及我国都还多多少少保留着区别适用的情况。
4 仿真实验
为验证算法有效性，本文在8×quad core 2.3 MHz CPU、64 GB内存、Window10系统环境下，利用CAWWA算法对包装配送实例进行仿真测试。其中，配送中心位置为(30 km,30 km)，30个客户需求量及位置信息如表1所示，车辆最大载重为8 t。
表1 客户信息
如表2所示，完成此次配送任务所需车辆数为8，最优配送距离为809.59 km。由此可知，CAWWA算法能够有效应用于包装配送问题的求解。为进一步验证算法的求解性能，分别利用CAWWA算法、GA算法[13]、ACO算法[14]和TS算法[15]对配送问题的国际通用算例进行50次求解，统计求得最好解(BS)、平均运行时间(AT)等指标如表3所示。各算例最优配送路径如表4所示。
表2 最优配送方案
表3 GA算法、ACO算法、TS算法和 CAWWA算法仿真结果对比
表4 最优配送方案
续表4
由表3可知， TS算法只能求得B-n31-k5算例的已知最优解，GA算法和ACO算法分别还能求得E-n23-k3、E-n33-k4算例的已知最优解，CAWWA算法不但能求得所有算例已知最优解，还更新了B-n31-k5、B-n51-k7和P-n55-k8的已知最优解。由此可见，CAWWA算法全局寻优能力最强，GA算法和ACO算法次之，而TS算法最弱。同时，CAWWA算法求解各算例的平均运行时间最短，说明该算法运行速度最快、搜索效率最高。这主要得益于混沌机制提高了算法初始种群质量，自适应碎波系数增强了算法前期全局搜索能力，并提高了算法后期的搜索精度，进而优化了算法综合求解性能。
在传统互联网上，数据主要用来表示信息，核心要义是传播、复制和分享。随着互联网与实体经济融合的不断深化、大数据的兴起，数据正在资产化，资产正在数据化。现在的数据，很多已经是用来表示价值而不是信息了。数据表示价值，核心要义是所有权、控制和交易。
综上分析，CAWWA算法能够很好地适应包装配送问题的求解，且求解性能优于GA算法、ACO算法和TS算法。
5 结 语
本文主要对求解包装配送问题的CAWWA算法进行了研究。首先，设计了基于包装配送问题的个体编码方式；然后，引入混沌机制生成初始种群，并提出了自适应调整的碎波系数，扩大算法前期搜索范围，提高后期搜索精度。仿真结果表明，CAWWA算法求解性能优于GA算法、TS算法和ACO算法，能够适应于包装配送问题的求解。
梁爽(1987-)，女，汉族，新疆乌鲁木齐人，经济系硕士，研究方向：国际收支，国际投资，国际货币流动，国际储备。
参考文献
[1] 樊贵香.混合模拟植物生长算法在包装件配送中的应用[J].包装工程,2016(13):43-49.
[2] Uchoa E, Pecin D, Pessoa A, et al. New benchmark instances for the capacitated vehicle routing problem[J]. European Journal of Operational Research, 2017, 257(3):845-858.
[3] 王帅,蒋华伟.遗传-蚁群算法在灾后应急物资路径规划问题中的应用研究[J].计算机应用与软件, 2018,35(9): 99-103,131.
[4] Wang K , Lan S , Zhao Y . A genetic-algorithm-based approach to the two-echelon capacitated vehicle routing problem with stochastic demands in logistics service[J]. Journal of the Operational Research Society,2017,68(11):1409-1421.
[5] Marinaki M, Marinakis Y. A glowworm swarm optimization algorithm for the vehicle routing problem with stochastic demands [J]. Expert Systems with Applications, 2016, 46(C):145-163.
[6] 裴小兵, 贾定芳. 基于模拟退火算法的城市物流多目标配送车辆路径优化研究[J]. 数学的实践与认识, 2016, 46(2):105-113.
[7] Teymourian E, Kayvanfar V, Komaki G M, et al. Enhanced intelligent water drops and cuckoo search algorithms for solving the capacitated vehicle routing problem [J]. Information Sciences, 2016, 334-335(C):354-378.
[8] 吴秀丽, 周永权. 一种基于混沌和单纯形法的水波优化算法[J]. 计算机科学, 2017, 44(5):218-225.
[9] Ai T J, Kachitvichyanukul V. Particle swarm optimization and two solution representations for solving the capacited vehicle routing problem [J]. Computers & Industrial Engineering, 2009, 56(1):380-387.
[10] 王万良, 陈超, 李笠,等. 基于模拟退火的自适应水波优化算法[J]. 计算机科学, 2017, 44(10):216-221.
[11] Costa D D A, Savi M A. Chaos control of an SMA-pendulum system using thermal actuation with extended time-delayed feedback approach [J]. Nonlinear Dynamics, 2018, 93(2):571-583.
[12] Mariotti C. A new leapfrog scheme for rotational motion in 3D [J]. International Journal for Numerical Methods in Engineering, 2016, 107(4):273-289.
[13] 周生伟,蒋同海,张荣辉.改进遗传算法求解VRP问题[J].计算机仿真,2013,30(12):140-143,157.
[14] Gianluca. Ant colony optimization for capacitated vehicle routing problem[J]. Journal of Computer Science, 2012, 8(6):846-852.
[15] Ren C Y. Tabu search algorithm for capacitated vehicle routing problem [J]. Advanced Materials Research, 2013, 753-755:3060-3063.
APPLICATION OF CHAOTIC ADAPTIVE WATER WAVE ALGORITHM IN PACKAGING AND DISTRIBUTION PROBLEM
Peng Wei
(College of Business Administration ,Chongqing City Management College ,Chongqing 401331 ,China )
Abstract The efficient heuristic algorithm is an important research direction of packaging and distribution problem. So this paper designed the chaotic adaptive water wave algorithm (CAWWA) to solve packaging and distribution problem. The individual expression of WWA was designed according to the characteristics of packaging and distribution problem. Then, the initial population of WWA was generated by the chaotic system, and an adaptive wave breaking coefficient was proposed, which was used to adjust the search range of algorithm on evolutionary algebra. Simulation experiments were carried out by using packaging distribution examples and 6 international examples. The results show that this algorithm can be applied to solve packaging and distribution problems. Compared with GA, TS and ACO, CAWWA has stronger global search ability and faster convergence speed.
Keywords Packaging and distribution Chaotic Adaptive Water wave algorithm
中图分类号 TP301
文献标识码 A
DOI: 10.3969/j.issn.1000-386x.2019.10.010
收稿日期： 2019-01-20。
重庆市教育委员会2017年度科学技术研究项目(1609155488)。彭维 ，讲师，主研领域：物流与供应链管理。
标签：包装配送论文; 混沌论文; 自适应论文; 水波算法论文; 重庆城市管理职业学院工商管理学院论文;