CN101739812A - 物流配送车辆路径优化方法 - Google Patents
物流配送车辆路径优化方法 Download PDFInfo
- Publication number
- CN101739812A CN101739812A CN200910155263A CN200910155263A CN101739812A CN 101739812 A CN101739812 A CN 101739812A CN 200910155263 A CN200910155263 A CN 200910155263A CN 200910155263 A CN200910155263 A CN 200910155263A CN 101739812 A CN101739812 A CN 101739812A
- Authority
- CN
- China
- Prior art keywords
- cost
- vehicle
- dispensing
- client
- route scheme
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 238000012546 transfer Methods 0.000 claims description 9
- 230000008521 reorganization Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000003780 insertion Methods 0.000 claims description 3
- 230000037431 insertion Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 abstract 2
- 230000000694 effects Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 239000012141 concentrate Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种物流配送车辆路径优化方法,包括下述步骤:步骤1,参数设置及初始化;步骤2:若Nc≤Nmax,判断迭代次数是否达到最大迭代次数,如果是,就结束循环,执行步骤7;否则,执行步骤3;步骤3,并购操作,再对物流配送路径进行调整操作;步骤4,重组操作,插入操作,调整操作;步骤5,成本评估;步骤6,选择最优物流配送路径方案;步骤7,获得最优物流配送车辆路径方案。本发明的物流配送路径优化方法是通过抽象企业并购行为,由并购操作,重组操作等组成的并购方法,提高物流配送路径优化的收敛性,可以高效地得到最优物流配送路径方案。
Description
(一)技术领域
本发明涉及一种物流配送车辆路径优化方法,具体涉及一种基于本专利
方法的有时间窗约束物流配送及车辆路径优化方法。
(二)背景技术
物流配送路径优化问题,即是车辆路径问题(Vehicle Routing Problem,VRP),是Dantzig和Ramser于1959年提出的。具体来说,就是根据物流配送的已知条件和目标约束条件,寻找一般最符合条件的配送路线,使物流配送成本最小化。随着市场经济的发展与繁荣,物流配送在货物运输成本中占比重也越来越高。因而,如何降低物流运输成本,提高物流运输效率等问题,日益凸显重要,成为影响物流业发展中的重要一环。特别是,商业竞争的残酷性,使得客户对时间的要求更为苛刻。因此,如何解决有时间窗物流配送路径优化问题,也相应地凸显迫切。目前,物流配送路径问题应用范围十分广泛,物流运输、交通建设、通信等。随着物流业的发展,根据现代物流业发展的需求,车辆路径问题又进一步发展了有时间窗车辆路径问题(Vehicle Routing Problem With Time Windows,VRPTW),或有时间窗物流配送路径优化问题,也即是在物流配送路径优化问题上,增加客户对配送时间的要求,目标仍然是求解问题的最优解,使得物流配送成本最小化。
目前,对物流配送路径优化问题的研究,部分趋于成熟,但多数还集中在利用传统算法来解决问题。其中,取得效果比较好的,有Thangiah等人采用遗传算法来解决物流配送路径优化问题,及有时间窗约束的物流配送路径优化问题。也有很多学者利用蚁群算法、禁忌算法等来求解同样问题,往往很多算法在求解该问题时,会遇到“陷入局部最优”等问题。在国内,研究比较有限,主要是在遗传算法基础上,提出各自的见解与主张。除此之外,也有很多学者在积极地研究新的算法来解决问题,效果也各有优劣,但还是没有传统算法那么成熟。因而,当前对物流配送路径优化问题的研究,处在一个比较缓慢的地步。
当前国内对物流配送路径优化问题研究,主要以传统算法为主。其主要存在以下不足之处:(1)、求解物流配送路径优化问题的算法研究的单一性,使得问题的解决也具体单一性;(2)、求解物流配送路径优化问题的算法中,有的效果好,但计算复杂度偏高;有的计算复杂度低,但效果不够理想;(3)、很多算法求解物流路径优化问题时,往往容易遇到“陷入局部最优”等问题。
(三)发明内容
本发明所要解决的技术问题是提供一种高效可行的物流配送路径优化方法。
一种物流配送车辆路径优化方法,包括下述步骤:
步骤1,参数设置及初始化:物流配送路径优化的方案数目NT;循环迭代次数Nc=1,最大迭代次数Nmax;当前参与并购操作的车辆标识CarMark=ones(NT,1),CarMark的选择顺序是根据步骤5中“单个车辆配送路径方案成本”从高到低排列;当前由车辆CarMark完成配送的客户标识CustomerMark=ones(NT,1),CustomerMark的选择顺序受步骤3中的并购操作中some函数影响;MACustSum表示由车辆CarMark完成车辆配送路径方案中临时选定参与并购操作客户的最大值,也受步骤3中的并购操作中some函数影响;Dist表示NT个物流配送方案;Cars表示参与物流配送方案的配送车辆的数目;MAType=zeros(NT,1)表示并购操作类型;先随机生成NT条物流配送方案,每一条物流配送方案有n个配送客户个体;其中,ones是赋值为1的函数,zeros是赋值为0的函数;
步骤2:若Nc≤Nmax,判断迭代次数是否达到最大迭代次数,如果是,就结束循环,执行步骤7;否则,执行步骤3;
步骤3,并购操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示车辆j的配送路径方案,则有:
若MAType(i)=0,执行以下步骤:
先对参与并购操作的配送对象的位置重定位:
若CustomerMark(i)>MACustSum(i),则CarMark(i)=CarMark(i)+1;
若CarMark(i)>Cars,则MAType(i)=1,跳出当前操作;
其中,MACustSum(i)=sum(some(CarDist(j))),表示客户个体集的数目;
再对物流配送路径进行调整操作:
物流配送方案i中当前参与并购操作的车辆CarMark(i),由它完成配送的客户个体CustomerMark(i),记作Temp,分别调至到除其以外的所有车辆中进行配送,分别计算出调整后的每一个车辆配送路径方案成本,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时消去车辆CarMark(i)对其的配送任务;
若MAType(i)=1,执行以下步骤:
提取由每一个车辆完成配送路径方案中的部分配送客户,组成一个待配送客户集合
选择其中一个客户个体Temp分别插入由每一个车辆配送路径方案中的适当位置,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时把客户个体Temp从待配送客户集合T删除,继续再选择一个客户个体,做同样操作,直到待配送客户集合T为空为止;
其中,some函数定义为:提取每一个车辆配送路径方案中的一定客户,组成一个临时客户个体集合,一定客户的选择原则是按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本mincost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集,而这个客户个体集的数目就是MACustSum;
步骤4,重组操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示j车辆的配送路径方案,则有:
插入操作:
对于每一个车辆j配送路径方案作如下操作:
由车辆j完成的配送路径方案中,按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本mincost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集T=some(CarDist(j)),依次把T中的每一个客户个体先后插入由车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
调整操作:
对于每一个车辆j配送路径方案作如下操作:
先后依次取出由车辆j完成配送路径方案中的每一个配送客户,再插入车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
步骤5,成本评估:
计算每一个物流配送路径方案的成本,成本包括三部分,分别是配送客户个体在单个车辆配送路径方案中的成本,单个车辆配送路径方案成本和物流配送方案成本;
配送客户个体在单个车辆配送路径方案中的成本包括距离成本cs、等待时间成本ce和迟到时间成本cl,即为(cs+ce+cl);
单个车辆配送路径方案成本,包括四部分,分别是又单个车辆完成配送任务所需的运输距离成本cs,运输超载成本cq,各个客户个体的等待时间成本之和ce和迟到时间成本之和cl,即为(cs+ce+cl)exp(cq);
物流配送路径方案成本,也包括四部分,分别是由完成所有客户配送所需的运输距离成本Cs、运输超载成本之和cq、各客户个体的等待时间成本之和ce和迟到时间成本之和cl,即为(Cs+Ce+Cl)exp(Cq);
步骤6,选择最优物流配送路径方案:
对于任意一个物流配送路径方案i,则有:
若MAType(i)=0,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,CarMark(i)=1,继续执行步骤2;否则,不更新,CustomerMark(i)=CustomerMark(i)+1,并在保存成本最小的物流配送路径方案,暂时作为当前最佳物流配送路径方案,继续执行步骤2;
若MAType(i)=1,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,MAType(i)=0,CarMark(i)=1,继续执行步骤2;否则,继续执行步骤2;
步骤7,获得最优物流配送车辆路径方案:
对于任意一个物流配送方案,按步骤5的“物流配送路径方案成本”方法求取最小成本minCost(1,2,…,NT),成本最低的物流配送路径方案即为最优物流配送车辆路径方案。
本发明中,配送客户个体在单个车辆配送路径方案中的距离成本,和其它运输距离成本的定义不同。运输距离成本就是完成配送任务所需的车程成本,而配送客户个体在单个车辆配送路径方案中的距离成本,则是当前配送客户个体和当前车辆配送路径方案中的其它客户个体之间距离,以及当前配送客户个体和当前车辆配送路径方案中的相邻客户个体之间距离的综合距离成本。
本发明的物流配送路径优化方法是通过抽象企业并购行为,由并购操作,重组操作等组成的并购方法,提高物流配送路径优化的收敛性,可以高效地得到最优物流配送路径方案。
(四)附图说明
图1为C101问题类型最优路线图;
图2为C201问题类型最优路线图;
图3为R101问题类型最优路线图;
图4为R201问题类型最优路线图;
图5为RC101问题类型最优路线图;
图6为RC201问题类型最优路线图。
(五)具体实施方式
下面通过附图和具体实施例对本发明的技术方案做进一步地详细地描
述,但本发明的保护范围并不限于此。
一种物流配送车辆路径优化方法,包括下述步骤:
步骤1,参数设置及初始化:物流配送路径优化的方案数目NT;循环迭代次数Nc=1,最大迭代次数Nmax;当前参与并购操作的车辆标识CarMark=ones(NT,1),CarMark的选择顺序是根据步骤5中“单个车辆配送路径方案成本”从高到低排列;当前由车辆CarMark完成配送的客户标识CustomerMark=ones(NT,1),CustomerMark的选择顺序受步骤3中的并购操作中some函数影响;MACustSum表示由车辆CarMark完成车辆配送路径方案中临时选定参与并购操作客户的最大值,也受步骤3中的并购操作中some函数影响;Dist表示NT个物流配送方案;Cars表示参与物流配送方案的配送车辆的数目;MAType=zeros(NT,1)表示并购操作类型;先随机生成NT条物流配送方案,每一条物流配送方案有n个配送客户个体;其中,ones是赋值为1的函数,zeros是赋值为0的函数;
步骤2:若Nc≤Nmax,判断迭代次数是否达到最大迭代次数,如果是,就结束循环,执行步骤7;否则,执行步骤3;
步骤3,并购操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示车辆j的配送路径方案,则有:
若MAType(i)=0,执行以下步骤:
先对参与并购操作的配送对象的位置重定位:
若CustomerMark(i)>MACustSum(i),则CarMark(i)=CarMark(i)+1;
若CarMark(i)>Cars,则MAType(i)=1,跳出当前操作;
其中,MACustSum(i)=sum(some(CarDist(j))),表示客户个体集的数目;
再对物流配送路径进行调整操作:
物流配送方案i中当前参与并购操作的车辆CarMark(i),由它完成配送的客户个体CustomerMark(i),记作Temp,分别调至到除其以外的所有车辆中进行配送,分别计算出调整后的每一个车辆配送路径方案成本,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时消去车辆CarMark(i)对其的配送任务;
若MAType(i)=1,执行以下步骤:
提取由每一个车辆完成配送路径方案中的部分配送客户,组成一个待配送客户集合
选择其中一个客户个体Temp分别插入由每一个车辆配送路径方案中的适当位置,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时把客户个体Temp从待配送客户集合T删除,继续再选择一个客户个体,做同样操作,直到待配送客户集合T为空为止;
其中,some函数定义为:提取每一个车辆配送路径方案中的一定客户,组成一个临时客户个体集合,一定客户的选择原则是按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本mincost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集,而这个客户个体集的数目就是MACustSum;
步骤4,重组操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示j车辆的配送路径方案,则有:
插入操作:
对于每一个车辆j配送路径方案作如下操作:
由车辆j完成的配送路径方案中,按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本mincost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集T=some(CarDist(j)),依次把T中的每一个客户个体先后插入由车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
调整操作:
对于每一个车辆j配送路径方案作如下操作:
先后依次取出由车辆j完成配送路径方案中的每一个配送客户,再插入车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
步骤5,成本评估:
计算每一个物流配送路径方案的成本,成本包括三部分,分别是配送客户个体在单个车辆配送路径方案中的成本,单个车辆配送路径方案成本和物流配送方案成本;
配送客户个体在单个车辆配送路径方案中的成本包括距离成本cs、等待时间成本ce和迟到时间成本cl,即为(cs+ce+cl);
单个车辆配送路径方案成本,包括四部分,分别是又单个车辆完成配送任务所需的运输距离成本cs,运输超载成本cq,各个客户个体的等待时间成本之和ce和迟到时间成本之和cl,即为(cs+ce+cl)exp(cq);
物流配送路径方案成本,也包括四部分,分别是由完成所有客户配送所需的运输距离成本Cs、运输超载成本之和Cq、各客户个体的等待时间成本之和Ce和迟到时间成本之和Cl,即为(Cs+Ce+Cl)exp(Cq);
步骤6,选择最优物流配送路径方案:
对于任意一个物流配送路径方案i,则有:
若MAType(i)=0,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,CarMark(i)=1,继续执行步骤2;否则,不更新,CustomerMark(i)=CustomerMark(i)+1,并在保存成本最小的物流配送路径方案,暂时作为当前最佳物流配送路径方案,继续执行步骤2;
若MAType(i)=1,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,MAType(i)=0,CarMark(i)=1,继续执行步骤2;否则,继续执行步骤2;
步骤7,获得最优物流配送车辆路径方案:
对于任意一个物流配送方案,按步骤5的“物流配送路径方案成本”方法求取最小成本minCost(1,2,…,NT),成本最低的物流配送路径方案即为最优物流配送车辆路径方案。
实验采用Solomon于1987年所设计的Solomon Benchmark problems的数据集。其中,每一个问题类型都有100个客户,一共有C1、C2、R1、R2、RC1、RC2等六大类的56种问题类型,每个类型的数据集合包含有7个到12个不等的带有100个配送点的问题数据。而且C问题类型的客户分布较聚群或集中,R问题类型的客户分布较均匀,RC问题类型是C问题类型和R问题类型的混合;C1、R1、RC1等I类问题的特点,车辆多,容量小,车程时间较短;C2、R2、RC2等II类问题的特点,则与前者恰好相反。
实验参数说明:(RT,EC,Wq,Wd,We,Wl)=(迭代次数,企业生态链数目,单位超载载量成本惩罚权重,单位距离成本惩罚权重,等待时间惩罚权重,迟到时间惩罚权重)。
实验结果与分析:
实验结果与目前已知的最优解(源自文献[Bent R,Van Hentenryck P. Atwo-stage hybrid local search for the vehicle routing problem with time windows[R].USA:Brown University,2001])进行比较,得到以下结果。其中,Data代表问题类型,NV代表需求的汽车数,TD代表汽车总的行驶路程,ET代表等待时间,LE代表迟到时间,N%表示第一目标与已经最优解的相对误差×100,T%表示第二目标与已经最优解的相对误差×100。
表1在Solomon数据集上的实验结果
实验结果分析:
C问题类型:特点是客户呈结构分布。可见,在本专利方法在解决C问题类型十分有效,第一目标和第二目标都获得了最优解,而且C104获得了更优解。
R问题类型和RC问题类型:R问题类型特点是客户分布比较均匀,RC问题类型是C问题类型和R问题类型的混合。R1和RC1问题类型都是有部分在第一目标上获得了最优解,但还是有几个问题类型不能达到最优解,不过在第二目标上效果也不错,接近最优解。主要原因是因为R和RC问题类型中客户分布比较均匀,所需车辆较多,而且时间窗宽窄不一,使得在本专利方法下所求得的解与目前最优解相比,还有进一步提高的空间。R2和RC2问题类型在第一目标上都获得到了最优解,而且第二目标也接近最优解,效果比较好。
上述实验结果充分证明,本专利方法在解决物流配送路径优化问题上具有可行性和有效性,而且可有效地解决其它方法在求解物流配送路径优化问题上遇到的“陷入局部最优”问题,同时也丰富了求解物流配送路径优化问题的方法。
Solomon几个常见问题类型详细实验结果
Solomon几个常见问题类型一般指C101、C201、R101、R201、RC101、RC201问题类型。
①C101问题类型:得到的最优解路程长度为TD=828.9396KM,等待时间ET=0min,迟到时间LT=0min,如图1所示;C201问题类型:得到的最优解路程长度为TD=591.5566KM,等待时间ET=0min,迟到时间LT=0min,如图2所示。
表2C101问题类型最优解的路线
1 | 0 20 24 25 27 29 30 28 26 23 22 21 0 |
2 | 0 67 65 63 62 74 72 61 64 68 66 69 0 |
3 | 0 57 55 54 53 56 58 60 59 0 |
4 | 0 43 42 41 40 44 46 45 48 51 50 52 49 47 0 |
5 | 0 81 78 76 71 70 73 77 79 80 0 |
6 | 0 5 3 7 8 10 11 9 6 4 2 1 75 0 |
7 | 0 13 17 18 19 15 16 14 12 0 |
8 | 0 98 96 95 94 92 93 97 100 99 0 |
9 | 0 32 33 31 35 37 38 39 36 34 0 |
10 | 0 90 87 86 83 82 84 85 88 89 91 0 |
表3C201问题类型最优解的路线
1 | 0 93 5 75 2 1 99 100 97 92 94 95 98 7 3 4 89 91 88 84 86 83 82 85 76 71 70 73 80 |
2 | 0 67 63 62 74 72 61 64 66 69 68 65 49 55 54 53 56 58 60 59 57 40 44 46 45 51 50 52 47 43 42 41 48 0 |
3 | 0 20 22 24 27 30 29 6 32 33 31 35 37 38 39 36 34 28 26 23 18 19 16 14 12 15 17 13 25 9 11 10 8 21 0 |
②R101问题类型:得到的最优解路程长度为TD=1739.8KM,等待时间ET=544.9547min,迟到时间LT=0min,如图3所示;R201问题类型:得到的最优解路程长度为TD=1277.8KM,等待时间ET=1077.0min,迟到时间LT=0min,如图4所示。
表4R101问题类型最优解的路线
编号 | 配送路径经过客户 | 编号 | 配送路径经过客户 |
1 | 0 52 99 94 97 0 | 11 | 0 31 30 81 34 35 0 |
2 | 0 14 44 38 43 13 0 | 12 | 0 92 95 98 16 85 37 93 0 |
3 | 0 5 82 18 84 17 0 | 13 | 0 33 29 79 50 0 |
4 | 0 36 47 7 8 46 60 89 0 | 14 | 0 45 83 61 86 91 100 0 |
5 | 0 11 19 49 48 0 | 15 | 0 2 75 22 55 4 25 0 |
6 | 0 65 71 20 32 70 0 | 16 | 0 72 21 73 41 26 0 |
7 | 0 63 64 90 10 0 | 17 | 0 62 88 6 96 0 |
8 | 0 27 69 51 9 66 1 0 | 18 | 0 59 42 15 87 57 58 0 |
9 | 0 28 12 76 78 3 68 77 0 | 19 | 0 39 23 67 56 74 0 |
10 | 0 40 53 54 24 80 0 |
表5R201问题类型最优解的路线
编号 | 配送路径经过客户 |
编号 | 配送路径经过客户 |
1 | 0 27 52 31 69 28 12 29 76 30 71 9 51 81 79 78 34 3 50 20 10 32 66 35 68 77 1 70 0 |
2 | 0 39 72 2 42 15 75 23 67 21 73 40 53 87 57 41 22 56 54 26 74 4 55 25 24 80 0 |
3 | 0 83 45 5 59 92 95 98 14 38 44 16 61 86 85 99 84 6 94 96 97 37 43 91 100 93 17 60 89 0 |
4 | 0 33 65 63 82 47 36 64 11 19 62 88 7 8 18 90 49 46 48 13 58 0 |
③RC101问题类型:得到的最优解路程长度为TD=1719.8KM,等待时间ET=131.8247min,迟到时间LT=0min,如图5所示;RC201问题类型:得到的最优解路程长度为TD=1503.3KM,等待时间ET=624.2545min,迟到时间LT=0min,如图6所示。
表6RC101问题类型最优解的路线
79 81 78 77 96 87 90 0 | |||
编号 | 配送路径经过客户 | 编号 | 配送路径经过客户 |
1 | 0 63 76 51 85 89 91 0 | 9 | 0 12 11 16 9 10 13 17 0 |
2 | 0 69 98 88 53 78 60 0 | 10 | 0 23 21 18 49 22 20 48 25 0 |
3 | 0 65 59 75 87 97 58 77 0 | 11 | 0 92 33 28 30 34 50 0 |
4 | 0 14 47 15 73 79 46 4 100 0 | 12 | 0 95 62 67 71 94 56 66 0 |
5 | 0 5 45 2 7 6 8 3 1 70 0 | 13 | 0 82 52 99 86 57 74 24 0 |
6 | 0 61 81 90 55 68 0 | 14 | 0 83 64 19 84 0 |
7 | 0 72 38 44 41 54 96 93 0 | 15 | 0 31 29 27 26 32 80 0 |
8 | 0 42 39 36 40 43 37 35 0 |
表7RC201问题类型最优解的路线。
编号 | 配送路径经过客户 |
编号 | 配送路径经过客户 |
1 | 0 72 92 95 63 33 28 27 29 31 30 62 67 71 90 6 7 8 46 3 4 1 55 68 0 |
2 | 0 64 83 52 75 23 21 18 19 76 85 84 51 49 22 57 20 66 56 50 32 89 48 24 25 77 58 0 |
3 | 0 39 36 42 45 5 2 69 82 98 88 61 44 40 38 41 81 94 34 26 96 54 43 35 37 93 91 80 0 |
4 | 0 65 59 14 47 16 15 11 12 73 79 78 53 99 86 87 9 10 97 74 13 17 60 100 70 0 |
Claims (1)
1.一种物流配送车辆路径优化方法,其特征在于包括下述步骤:
步骤1,参数设置及初始化:物流配送路径优化的方案数目NT;循环迭代次数Nc=1,最大迭代次数Nmax;当前参与并购操作的车辆标识CarMark=ones(NT,1),CarMark的选择顺序是根据步骤5中“单个车辆配送路径方案成本”从高到低排列;当前由车辆CarMark完成配送的客户标识CustomerMark=ones(NT,1),CustomerMark的选择顺序受步骤3中的并购操作中some函数影响;MACustSum表示由车辆CarMark完成车辆配送路径方案中临时选定参与并购操作客户的最大值,也受步骤3中的并购操作中some函数影响;Dist表示NT个物流配送方案;Cars表示参与物流配送方案的配送车辆的数目;MAType=zeros(NT,1)表示并购操作类型;先随机生成NT条物流配送方案,每一条物流配送方案有n个配送客户个体;其中,ones是赋值为1的函数,zeros是赋值为0的函数;
步骤2:若Nc≤Nmax,判断迭代次数是否达到最大迭代次数,如果是,就结束循环,执行步骤7;否则,执行步骤3;
步骤3,并购操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示车辆j的配送路径方案,则有:
若MAType(i)=0,执行以下步骤:
先对参与并购操作的配送对象的位置重定位:
若CustomerMark(i)>MACustSum(i),则CarMark(i)=CarMark(i)+1;
若CarMark(i)>Cars,则MAType(i)=1,跳出当前操作;
其中,MACustSum(i)=sum(some(CarDist(j))),表示客户个体集的数目;
再对物流配送路径进行调整操作:
物流配送方案i中当前参与并购操作的车辆CarMark(i),由它完成配送的客户个体CustomerMark(i),记作Temp,分别调至到除其以外的所有车辆中进行配送,分别计算出调整后的每一个车辆配送路径方案成本,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时消去车辆CarMark(i)对其的配送任务;
若MAType(i)=1,执行以下步骤:
提取由每一个车辆完成配送路径方案中的部分配送客户,组成一个待配送客户集合
选择其中一个客户个体Temp分别插入由每一个车辆配送路径方案中的适当位置,按步骤5的“单个车辆配送路径方案成本”方法求取最小成本minCost(1,2,…,Cars),选择配送成本最小的车辆k,最终确认把配送客户个体Temp调至到由车辆k完成配送,同时把客户个体Temp从待配送客户集合T删除,继续再选择一个客户个体,做同样操作,直到待配送客户集合T为空为止;
其中,some函数定义为:提取每一个车辆配送路径方案中的一定客户,组成一个临时客户个体集合,一定客户的选择原则是按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本min cost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集,而这个客户个体集的数目就是MACustSum;
步骤4,重组操作:
对于任意一个物流配送方案i,令CarDist=Dist(i),CarDist(j)表示j车辆的配送路径方案,则有:
插入操作:
对于每一个车辆j配送路径方案作如下操作:
由车辆j完成的配送路径方案中,按步骤5的“配送客户个体在单个车辆配送路径方案中的成本”方法求取最小成本min cost(1,2,…,x),选择最小的成本隶属客户个体p,以p为临界,选择包括p内的左侧(或右侧)所以客户中的不大于二分之一个客户个体集T=some(CarDist(j)),依次把T中的每一个客户个体先后插入由车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
调整操作:
对于每一个车辆j配送路径方案作如下操作:
先后依次取出由车辆j完成配送路径方案中的每一个配送客户,再插入车辆j完成配送路径方案中,同时保证当前车辆j的“单个车辆配送路径方案成本”为最小成本min(Cost(CarDist(j))),使得当前由车辆j完成配送路径方案的成本也是最小的;
步骤5,成本评估:
计算每一个物流配送路径方案的成本,成本包括三部分,分别是配送客户个体在单个车辆配送路径方案中的成本,单个车辆配送路径方案成本和物流配送方案成本;
配送客户个体在单个车辆配送路径方案中的成本包括距离成本cs、等待时间成本ce和迟到时间成本cl,即为(cs+ce+cl);
单个车辆配送路径方案成本,包括四部分,分别是又单个车辆完成配送任务所需的运输距离成本cs,运输超载成本cq,各个客户个体的等待时间成本之和ce和迟到时间成本之和cl,即为(cs+ce+cl)exp(cq);
物流配送路径方案成本,也包括四部分,分别是由完成所有客户配送所需的运输距离成本Cs、运输超载成本之和Cq、各客户个体的等待时间成本之和Ce和迟到时间成本之和Cl,即为(Cs+Ce+Cl)exp(Cq);
步骤6,选择最优物流配送路径方案:
对于任意一个物流配送路径方案i,则有:
若MAType(i)=0,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,CarMark(i)=1,继续执行步骤2;否则,不更新,CustomerMark(i)=CustomerMark(i)+1,并在保存成本最小的物流配送路径方案,暂时作为当前最佳物流配送路径方案,继续执行步骤2;
若MAType(i)=1,执行以下步骤:
根据得到最新物流配送路径方案成本,如果低于旧的物流配送路径方案成本,则更物流配送路径方案,MAType(i)=0,CarMark(i)=1,继续执行步骤2;否则,继续执行步骤2;
步骤7,获得最优物流配送车辆路径方案:
对于任意一个物流配送方案,按步骤5的“物流配送路径方案成本”方法求取最小成本minCost(1,2,…,NT),成本最低的物流配送路径方案即为最优物流配送车辆路径方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910155263A CN101739812A (zh) | 2009-12-10 | 2009-12-10 | 物流配送车辆路径优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910155263A CN101739812A (zh) | 2009-12-10 | 2009-12-10 | 物流配送车辆路径优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101739812A true CN101739812A (zh) | 2010-06-16 |
Family
ID=42463241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910155263A Pending CN101739812A (zh) | 2009-12-10 | 2009-12-10 | 物流配送车辆路径优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101739812A (zh) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445208A (zh) * | 2011-09-21 | 2012-05-09 | 清华大学深圳研究生院 | 一种从地图数据中获取多车辆导航路径的方法 |
CN103364005A (zh) * | 2012-03-30 | 2013-10-23 | 本田技研工业株式会社 | 集散支援系统 |
CN103426074A (zh) * | 2013-08-22 | 2013-12-04 | 深圳市华傲数据技术有限公司 | 一种计算物流收发路径的方法及系统 |
CN103440524A (zh) * | 2013-04-03 | 2013-12-11 | 天津大学 | 一种基于电子地图的物流路径优化方法 |
WO2015027734A1 (zh) * | 2013-09-02 | 2015-03-05 | 深圳市华傲数据技术有限公司 | 一种物流派件调度方法、平台和系统 |
CN104504453A (zh) * | 2014-11-20 | 2015-04-08 | 昆明理工大学 | 一种多目标优化的军事运输过程的优化调度方法 |
CN105631530A (zh) * | 2014-11-24 | 2016-06-01 | Sap欧洲公司 | 时间可划分资源的多顺序规划和分配 |
CN106600180A (zh) * | 2016-11-04 | 2017-04-26 | 中联物流(中国)有限公司 | 物流管理系统的运输路径优化方法和系统 |
CN107194575A (zh) * | 2017-05-19 | 2017-09-22 | 重庆大学 | 一种处理快递新增取件需求的车辆自治调度方法 |
CN107749010A (zh) * | 2016-08-29 | 2018-03-02 | 北京三快在线科技有限公司 | 订单处理方法、装置及服务器 |
CN108253988A (zh) * | 2017-12-29 | 2018-07-06 | 重庆大学 | 基于电动汽车的循环取货路径规划方法 |
CN108805488A (zh) * | 2018-05-21 | 2018-11-13 | 湖北迪天信息科技有限公司 | 一种物流承运社交管理平台 |
CN109978471A (zh) * | 2019-04-11 | 2019-07-05 | 聊城大学 | 一种带时间窗的冷链物流路径优化方法 |
CN109990795A (zh) * | 2019-04-01 | 2019-07-09 | 杉数科技(北京)有限公司 | 确定车辆开行距离的方法及装置 |
CN111288991A (zh) * | 2018-12-06 | 2020-06-16 | 北京京东尚科信息技术有限公司 | 路径规划方法、装置、机器人和计算机可读存储介质 |
CN111325371A (zh) * | 2018-12-13 | 2020-06-23 | 顺丰科技有限公司 | 一种运输线路规划方法及系统 |
CN112288347A (zh) * | 2020-02-21 | 2021-01-29 | 北京京东振世信息技术有限公司 | 冷链配送的路线确定方法、装置、服务器及存储介质 |
CN112862212A (zh) * | 2021-03-09 | 2021-05-28 | 无锡弘宜智能科技有限公司 | 基于改进麻雀搜索算法的多agv调度方法、装置及设备 |
CN113361739A (zh) * | 2020-03-06 | 2021-09-07 | 北京京东振世信息技术有限公司 | 一种生成拣货路径的方法和装置 |
CN113723804A (zh) * | 2021-08-30 | 2021-11-30 | 合肥工业大学 | 考虑多无人机站点的车机协同配送方法及系统 |
-
2009
- 2009-12-10 CN CN200910155263A patent/CN101739812A/zh active Pending
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445208A (zh) * | 2011-09-21 | 2012-05-09 | 清华大学深圳研究生院 | 一种从地图数据中获取多车辆导航路径的方法 |
CN103364005A (zh) * | 2012-03-30 | 2013-10-23 | 本田技研工业株式会社 | 集散支援系统 |
CN103440524A (zh) * | 2013-04-03 | 2013-12-11 | 天津大学 | 一种基于电子地图的物流路径优化方法 |
CN103426074A (zh) * | 2013-08-22 | 2013-12-04 | 深圳市华傲数据技术有限公司 | 一种计算物流收发路径的方法及系统 |
WO2015027734A1 (zh) * | 2013-09-02 | 2015-03-05 | 深圳市华傲数据技术有限公司 | 一种物流派件调度方法、平台和系统 |
CN104504453A (zh) * | 2014-11-20 | 2015-04-08 | 昆明理工大学 | 一种多目标优化的军事运输过程的优化调度方法 |
CN105631530A (zh) * | 2014-11-24 | 2016-06-01 | Sap欧洲公司 | 时间可划分资源的多顺序规划和分配 |
CN107749010A (zh) * | 2016-08-29 | 2018-03-02 | 北京三快在线科技有限公司 | 订单处理方法、装置及服务器 |
CN106600180A (zh) * | 2016-11-04 | 2017-04-26 | 中联物流(中国)有限公司 | 物流管理系统的运输路径优化方法和系统 |
CN106600180B (zh) * | 2016-11-04 | 2018-01-02 | 中联物流(中国)有限公司 | 物流管理系统的运输路径优化方法和系统 |
CN107194575A (zh) * | 2017-05-19 | 2017-09-22 | 重庆大学 | 一种处理快递新增取件需求的车辆自治调度方法 |
CN107194575B (zh) * | 2017-05-19 | 2020-07-31 | 重庆大学 | 一种处理快递新增取件需求的车辆自治调度方法 |
CN108253988A (zh) * | 2017-12-29 | 2018-07-06 | 重庆大学 | 基于电动汽车的循环取货路径规划方法 |
CN108805488A (zh) * | 2018-05-21 | 2018-11-13 | 湖北迪天信息科技有限公司 | 一种物流承运社交管理平台 |
CN111288991A (zh) * | 2018-12-06 | 2020-06-16 | 北京京东尚科信息技术有限公司 | 路径规划方法、装置、机器人和计算机可读存储介质 |
CN111325371B (zh) * | 2018-12-13 | 2023-04-18 | 顺丰科技有限公司 | 一种运输线路规划方法及系统 |
CN111325371A (zh) * | 2018-12-13 | 2020-06-23 | 顺丰科技有限公司 | 一种运输线路规划方法及系统 |
CN109990795A (zh) * | 2019-04-01 | 2019-07-09 | 杉数科技(北京)有限公司 | 确定车辆开行距离的方法及装置 |
CN109990795B (zh) * | 2019-04-01 | 2020-01-21 | 杉数科技(北京)有限公司 | 确定车辆开行距离的方法及装置 |
CN109978471A (zh) * | 2019-04-11 | 2019-07-05 | 聊城大学 | 一种带时间窗的冷链物流路径优化方法 |
CN112288347A (zh) * | 2020-02-21 | 2021-01-29 | 北京京东振世信息技术有限公司 | 冷链配送的路线确定方法、装置、服务器及存储介质 |
WO2021164390A1 (zh) * | 2020-02-21 | 2021-08-26 | 北京京东振世信息技术有限公司 | 冷链配送的路线确定方法、装置、服务器及存储介质 |
KR20220059531A (ko) * | 2020-02-21 | 2022-05-10 | 베이징 징동 젼스 인포메이션 테크놀로지 씨오., 엘티디. | 콜드 체인 배송 노선 확정 방법, 장치, 서버 및 저장매체 |
KR102655411B1 (ko) * | 2020-02-21 | 2024-04-04 | 베이징 징동 젼스 인포메이션 테크놀로지 씨오., 엘티디. | 콜드 체인 배송 노선 확정 방법, 장치, 서버 및 저장매체 |
CN113361739A (zh) * | 2020-03-06 | 2021-09-07 | 北京京东振世信息技术有限公司 | 一种生成拣货路径的方法和装置 |
CN113361739B (zh) * | 2020-03-06 | 2024-04-12 | 北京京东振世信息技术有限公司 | 一种生成拣货路径的方法和装置 |
CN112862212A (zh) * | 2021-03-09 | 2021-05-28 | 无锡弘宜智能科技有限公司 | 基于改进麻雀搜索算法的多agv调度方法、装置及设备 |
CN113723804A (zh) * | 2021-08-30 | 2021-11-30 | 合肥工业大学 | 考虑多无人机站点的车机协同配送方法及系统 |
CN113723804B (zh) * | 2021-08-30 | 2023-08-22 | 合肥工业大学 | 考虑多无人机站点的车机协同配送方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101739812A (zh) | 物流配送车辆路径优化方法 | |
Sun et al. | Behaviorally stable vehicle platooning for energy savings | |
Defryn et al. | Multi-objective optimisation models for the travelling salesman problem with horizontal cooperation | |
CN109919376A (zh) | 多车场与多车型车辆路径调度控制方法 | |
CN108764634A (zh) | 一种考虑充电需求增长的电动汽车充电站动态规划方法 | |
Rimmel et al. | Optimization of the nested monte-carlo algorithm on the traveling salesman problem with time windows | |
CN107301504A (zh) | 基于混合蛙跳—路径重连的生产运输协同调度方法和系统 | |
CN105489002A (zh) | 一种基于智能匹配和路径优化的拼车方法及系统 | |
CN108564211B (zh) | 物流运输路径规划方法及系统 | |
CN108253988A (zh) | 基于电动汽车的循环取货路径规划方法 | |
CN104111642A (zh) | 设备预防性维护与柔性作业车间控制集成优化方法 | |
Li et al. | Deploying autonomous mobile lockers in a two-echelon parcel operation | |
CN106779570A (zh) | 一种智能冷链物流路径多目标优化算法 | |
CN104504453A (zh) | 一种多目标优化的军事运输过程的优化调度方法 | |
CN101650805B (zh) | 用于路径规划的搜索方法 | |
Shi et al. | A home health care routing problem with stochastic travel and service time | |
CN114444843A (zh) | 一种基于大规模变邻域搜索策略的农产品绿色物流配送车辆调度方法及系统 | |
CN108416482B (zh) | 一种基于区域性多策略的物流配送路径规划方法 | |
CN109635998A (zh) | 一种求解带时间窗车辆路径问题的自适应多目标优化方法 | |
CN104537425A (zh) | 一种汽车空调出风口的生产装配过程的优化调度方法 | |
Yang et al. | User-based battery swapping strategy in an electric bike-sharing system | |
CN113344267A (zh) | 一种基于合作的物流网络资源配置优化方法 | |
Nguyen et al. | GRASP with learning process for a two-echelon location routing problem | |
Yang et al. | Comparison study on algorithms for vehicle routing problem with time windows | |
CN107423810A (zh) | 基于骆驼群算法的作业车间调度排产方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20100616 |