一种求解单配送中心物流运输调度问题的方法技术

本发明专利技术公开了一种求解单配送中心物流运输调度问题的方法。针对单配送中心物流运输调度问题，建立数学模型，提出并行混沌粒子群算法求解该模型。仿真分析表明，本发明专利技术提供的方法保持了群体的多样性，避免了过早陷入早熟，收敛速度快，能够有效地求解单配送中心物流运输调度问题。

【技术实现步骤摘要】

本专利技术属于人工智能物流路径规划控制领域，采用并行混沌粒子群算法来求解单配送中心物流运输调度问题的方法。
技术介绍
物流运输调度问题一直以来都是作为物流配送的热点研究问题，随着科学和社会的高速发展，它在国家经济发展和企业经营中扮演越来越重要的角色，被广泛应用于交通、工业管理、物流运输等领域。国内外学者分别从路径构造、局部搜索、数学规划等方向对其展开了研究。送货路径优化问题是指对产品的发货点和卸货点，通过选取适当的车辆行驶线路，在满足货物体积、需求量、送货时间、时间窗、车辆单次送货里程等约束条件下，以配送的车辆数最少、总费用尽量少、配送路径最短、配送所花时间最短等为目标。配送优化问题也称为车辆路径问题(VRP)被证明为NP难问题。粒子群优化算法是一种基于社会群体行为的演化算法，其源于鸟群和鱼群群体运动行为的研究。PSO算法与其他的算法一样，容易陷入局部最优解。粒子群算法出现早熟现象的主要原因四缺乏种群多样性，混沌是非线性系统中的普遍现象，它具有一些特异性质，如混沌运动能在一定的范围内按自身不重复地便利所有状态，而初始值条件极其微弱的变化会引起系统行为巨大的变化，利用上述两点，混沌游湖能表面粒子群算法陷入局部解，因此提出一种并行混沌粒子群算法，该算法保证了群体多样性从而粒子群过早陷入了早熟，仿真结果表明，该算法具有较高的收敛速度和收敛精度。
技术实现思路
本专利技术针对现有求解物流运输调度问题的不足，公开了一种求解单配送中心物流运输调度方法，建立一种更能表现物流运输实际问题的数学模型，并提出并行混沌粒子群算法对问题进行求解，对经典的车辆路径寻优问题。为实现上述目的，本专利技术提供的方案的具体步骤如下：1、一种求解单配送中心物流运输调度问题的方法，其特征在于，包括以下步骤：S1：建立单配送中心物流运输调度问题的数学模型：minz=Σk=1mΣi=1NΣj=1N(Cv+Cd)dijxijk]]>m=[Σi=1Ngi/αG]+1]]>满足：Σi=1Nqiyik≤Gk=1,2,...,m]]>Σk=1myik=1i=1,2,...,n]]>Σk=1my0k=m]]>Σi=1Nxijk=yjki,j=1,2,...,Nk=1,2,...,m]]>Σj=1Nxijk=yiki,j=1,2,...,Nk=1,2,...,m]]>其中，N表示要服务的客户数(N为正整数)，m表示配送中心发出的车辆数(m为正整数)，第i个客户的需求量为gi(i＝1，2，…，N)，式中[]代表不大于括号内的数字的最大正整数，α(0＜α＜1)是对装车和卸车的困难程度以及约束的估计，G表示每一辆车最大载重约束(吨)，Cv为单位距离车辆行驶费用(元/km)，Cd为司机单位费用(元/km)。dij表示车辆k行驶了客户i到j的路程(km)，目标函数公式(1)表示物流运输过程中的成本耗费，包括物流运输中车辆行驶费用和司机薪酬费用。式(3)表示各车辆承载的总重量不超过车辆的最大负载约束G；式(4)表示一个客户的货物只能由一辆车进行配送服务；式(5)表示物流公司车场的送货车辆的总数为m辆；式(6)表和式(7)说明了两个变量之间的关系。S2：采用并行混沌粒子群算法求解，进一步包括：S2-1:首先假设整个粒子群的最优解的个数为n，并且随机生成一列初始值：粒子的更迭在主进程和从进程同时进行，粒子群x，v是在主进程进行初始化，第i个变量xi所生成的初始值有ki个。S2-2:数据从主进程发送给从进程。S2-3:从进程中用混沌粒子群算子优化粒子群。S2-4:适应度评估。根据公式计算得到了每个粒子的适应值S2-5:通过统计多次迭代后每一代粒子个体的适应值得到局部最优适应值S2-6:将从进程得到的局部最优粒子值送往主进程，得出全局最优适应值gbest，如果当前进化代数<最大进化代数，或者，当前平均适应值-上一代平均适应值>允许值，主进程将gbest发送到从进程，完成指定的迭代次数后转向S2-3进行混沌优化；否则转向第S2-8。S2-7:对症粒子群进行优化，更新粒子群速度和位置；转向第S2-4。S2-8:退出循环。输入：①存储在计算机pi(i＝1,2,...,n)初始化数据Di②混沌生长模型。输出：全局最优解gbest。优化粒子群算法初始值公式如下：(1)将映射到区间[0,1]内，即：yj,i(0)=xj,i(0)-aibi-aii=1,2,...,n;j=1,2,...,k]]>(2)生成了本次迭代的混沌序列，以为初始值，根据logistic方程yn+1＝vyn(1-yn)求解得到序列其中(3)生成粒子成员作逆映射得到序列xj,i(0,l)=ai+(bi-ai)yj,i(0,l),i=1,2,...,nj=1,2,...,ki;l=1,2,...,Mj,i]]>粒子的速度和位置更新公式如下：vi+1=ω*vi+c1*rand()*(piBest-xi)+c2*rand()*(gbest-xi)]]>xi+1＝xi+vi+1其中，xi和xi+1分别代表第i代和i+1代粒子群的当前位置和更新位置，vi和vi+1分别代表第i代和i+1代粒子群的当前速度和更新速度。piBest是粒子本身所找到的最优解，gbest是这个整群目前找到的最优解，c1和c2是学习因子，rand()∈(0,1)，ω∈(0.1,0.9)是加权系数。附图说明图1本专利技术实施例的客户及配送中心位置；图2本专利技术实施例的最优配送轨迹图。具体实施方式为了验证本专利技术的优越性，下面将本专利技术通过1个配送中心对10个客户的配送进行实例分析，对本专利技术的的技术方案进行清楚、完整地描述。本专利技术的具体实施步骤如下：S1：获取但配送中心物流运输调度问题基本参数：表1表示的是客户的位置。车场的送货车辆最大运载能力为12吨，每次出行的最大里程为150km。不及空载费用，不考虑吃饭时间。车场坐标为(0，0)。车辆最早发车时间是8:00。S2：建立单配送中心物流运输调度问题的数学模型：minz=Σk=1mΣi=1NΣj=1N(Cv+Cd)dijxijk]]>m=[Σi=1Ngi/αG]+1]]>满足：Σi=1Nqiyik≤Gk=1,2,...,m]]>Σk=1myik=1i=1,2,...,n]]>Σk=1my0k=m]]>Σi=1Nxijk=yjki,j=1,2,...,Nk=1,2,...,m]]>Σj=1Nxijk=yiki,j=1,2,...,Nk=1,2,...,m]]>其中，N表示要服务的客户数(N为正整数)，m表示配送中心发出的车辆数(m为正整数)，第i个客户的需求量为gi(i＝1，2，…，N)，式中[]代表不大于括号内的数字的最大正整数，α(0＜α＜1)是对装车和卸车的困难程度以及约束的估计，G表示每一辆车最大载重约束(吨)，Cv为单位距离车辆行驶费用(元/km)，Cd为司机单位费用(元/km)。dij表示本文档来自技高网...

【技术保护点】
一种求解单配送中心物流运输调度问题的方法，其特征在于，包括以下步骤：S1：建立单配送中心物流运输调度问题的数学模型：满足：其中，N表示要服务的客户数(N为正整数)，m表示配送中心发出的车辆数(m为正整数)，第i个客户的需求量为gi(i＝1，2，…，N)，式中[]代表不大于括号内的数字的最大正整数，α(0＜α＜1)是对装车和卸车的困难程度以及约束的估计，G表示每一辆车最大载重约束(吨)，Cv为单位距离车辆行驶费用(元/km)，Cd为司机单位费用(元/km)。dij表示车辆k行驶了客户i到j的路程(km)，目标函数公式(1)表示物流运输过程中的成本耗费，包括物流运输中车辆行驶费用和司机薪酬费用；式(3)表示各车辆承载的总重量不超过车辆的最大负载约束G；式(4)表示一个客户的货物只能由一辆车进行配送服务；式(5)表示物流公司车场的送货车辆的总数为m辆；式(6)表和式(7)说明了两个变量之间的关系；S2：采用并行混沌粒子群算法求解，进一步包括：S2‑1:首先假设整个粒子群的最优解的个数为n，并且随机生成一列初始值：i＝1,2,...,n。粒子的更迭在主进程和从进程同时进行，粒子群x，v是在主进程进行初始化，第i个变量xi所生成的初始值有ki个；S2‑2:数据从主进程发送给从进程；S2‑3:从进程中用混沌粒子群算子优化粒子群；S2‑4:适应度评估。根据公式计算得到了每个粒子的适应值S2‑5:通过统计多次迭代后每一代粒子个体的适应值得到局部最优适应值S2‑6:将从进程得到的局部最优粒子值送往主进程，得出全局最优适应值gbest，如果当前进化代数<最大进化代数，或者，当前平均适应值‑上一代平均适应值>允许值，主进程将gbest发送到从进程，完成指定的迭代次数后转向S2‑3进行混沌优化；否则转向第S2‑8；S2‑7:对症粒子群进行优化，更新粒子群速度和位置；转向第S2‑4；S2‑8:退出循环，输入：①存储在计算机pi(i＝1,2,...,n)初始化数据Di②混沌生长模型，输出：全局最优解gbest。...

【技术特征摘要】
1.一种求解单配送中心物流运输调度问题的方法，其特征在于，包括以下步骤：S1：建立单配送中心物流运输调度问题的数学模型：满足：其中，N表示要服务的客户数(N为正整数)，m表示配送中心发出的车辆数(m为正整数)，第i个客户的需求量为gi(i＝1，2，…，N)，式中[]代表不大于括号内的数字的最大正整数，α(0＜α＜1)是对装车和卸车的困难程度以及约束的估计，G表示每一辆车最大载重约束(吨)，Cv为单位距离车辆行驶费用(元/km)，Cd为司机单位费用(元/km)。dij表示车辆k行驶了客户i到j的路程(km)，目标函数公式(1)表示物流运输过程中的成本耗费，包括物流运输中车辆行驶费用和司机薪酬费用；式(3)表示各车辆承载的总重量不超过车辆的最大负载约束G；式(4)表示一个客户的货物只能由一辆车进行配送服务；式(5)表示物流公司车场的送货车辆的总数为m辆；式(6)表和式(7)说明了两个变量之间的关系；S2：采用并行混沌粒子群算法求解，进一步包括：S2-1:首先假设整个粒子群的最优解的个数为n，并且随机生成一列初始值：i＝1,2,...,n。粒子的更迭在主进程和从进程同时进行，粒子群x，v是在主进程进行初始化，第i个变量xi所生成的初始值有ki个；S2-2:数据从主进程发送给从进程；S2-3:从进程中用混沌粒子群算子优化粒子群；S2-4:适应度评估。根据公式计算得到了每个粒子的适应值S2-5:通过统计多次迭代后每一代粒子个体的适应值得到局部最优适应值S2-6:将从进程得到的局部最优粒子值送往主进程，得...

【专利技术属性】
技术研发人员：蔡延光，戚远航，蔡颢，马骏，
申请(专利权)人：广东工业大学，
类型：发明
国别省市：广东;44