CN105260785B - Logistics distribution vehicle path optimization method based on improved cuckoo algorithm - Google Patents
Logistics distribution vehicle path optimization method based on improved cuckoo algorithm Download PDFInfo
- Publication number
- CN105260785B CN105260785B CN201510521483.4A CN201510521483A CN105260785B CN 105260785 B CN105260785 B CN 105260785B CN 201510521483 A CN201510521483 A CN 201510521483A CN 105260785 B CN105260785 B CN 105260785B
- Authority
- CN
- China
- Prior art keywords
- path
- parasitic
- algorithm
- vehicle
- distribution
- 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.)
- Active
Links
Images
Abstract
The invention relates to a logistics distribution vehicle path optimization method based on an improved cuckoo algorithm, which is characterized by comprising the following steps of: step 1: setting parameters for improving the cuckoo algorithm; step 2: initializing a population and calculating a fitness value; wherein the population randomly generates size number of parasitic nests with the nd dimension search space range of [1, number ] by adopting a uniform distribution function; calculating the fitness value, firstly, locally optimizing the path in each line by adopting a 2-opt algorithm, and solving the sum of the path values of each optimized line by adopting a piecewise penalty function method to serve as the fitness value of the scheme; and step 3: executing levy walking operation of the parasitic nest; and 4, step 4: performing a parasitic nest discovery operation; and 5: performing a parasitic nest mutation operation; step 6: the discovery probability is dynamically adjusted. The method is applied to the solution of the logistics distribution vehicle path problem by combining the improved cuckoo algorithm and the 2-opt algorithm, provides a new feasible and effective solution for solving the logistics distribution vehicle optimization problem, and enriches the method for solving the logistics distribution path optimization problem.
Description
Technical Field
The invention relates to a logistics distribution vehicle path optimization method, in particular to a logistics distribution vehicle path optimization method based on an improved cuckoo algorithm.
Background
With the rapid development of modern economy and network technology, the logistics industry has become the "third source of profit" for modern enterprises. The total logistics cost of China in 2011 is up to 8.5 trillion yuan, which accounts for 17.8% of GDP. In each link of logistics, the transportation and distribution cost accounts for about 60% of the total logistics cost, and the overhigh logistics cost restricts the development of national economy and simultaneously weakens the market competitiveness of enterprises. The logistics distribution vehicle path optimization problem (VRP) is one of the most critical links in the logistics distribution optimization. The logistics distribution vehicle path optimization problem is a typical NP-hard problem, which was proposed by DantZig and Ramser in 1959, and mainly researches how to realize optimization of vehicle distribution cost under the condition of meeting customer requirements and other constraint conditions (such as vehicle maximum load, vehicle maximum forming distance and the like) in the logistics vehicle distribution process, such as: shortest path, least cost and the like. The problem has become a research hotspot in the field of operational research and combinatorial optimization.
In recent years, research on the logistics vehicle distribution path problem has mainly focused on solving the problem by adopting various heuristic algorithms. The King-iron monarch and the like provide a chaotic particle swarm optimization algorithm for solving the optimization problem of the logistics distribution path by combining the advantages of chaotic ergodicity and rapidity of particle swarm; the method for solving the logistics distribution optimization problem by improving the genetic algorithm proposed by Audort and the like provides an effective way for solving the optimization problem related to the logistics distribution optimization problem; wuyuanchun provides a method for applying a self-adaptive variation particle swarm algorithm to the problem optimization of a logistics distribution path; wang Huadong, Li Wei and the like propose a logistics distribution path optimization method of a particle swarm algorithm. The research on the logistics distribution optimization problem is limited to the application of traditional intelligent colony algorithms such as a genetic algorithm, a particle swarm algorithm, an ant colony algorithm and the like in the logistics distribution optimization problem. However, in solving the problem of the logistics distribution vehicle path, the single algorithm is prone to be trapped in a local optimal value, and optimization accuracy is low. The cuckoo algorithm adopts a Levy searching mode, so that the implementation is simple, the required set parameters are few, and the optimization precision and the convergence speed are superior to those of a particle swarm algorithm and a genetic algorithm.
Disclosure of Invention
The invention aims to solve the problem of realizing the route optimization of the logistics distribution vehicle by adopting a method based on an improved cuckoo algorithm. The technical scheme is as follows:
a method for solving the problem of improving the route of a logistics distribution vehicle based on a cuckoo algorithm comprises the following steps:
step 1: setting parameters of an algorithm, setting a population size, the number nd of service customers, the number of vehicles required for solving the problem, the load CarrayCarCan of the vehicles, and the discovery probability of the bird eggs in the parasitic nestsAlgorithm search spatial Range [1, number ]]The number of times of population iteration Max _ iter, and the iteration counter N _ iter = 1;
step 2: initializing population and calculating fitness value, and randomly generating size nd dimension search space range of [1, number ] by adopting uniform distribution function]The position of the ith parasitic nest is recorded asAnd rounding each nest with a ceil function (i.e.) To ensure that each customer site is served by a vehicle and calculate its corresponding fitness value;
And step 3: performing levy walk operations of parasitic nests, usingAnd the updating mode generates a new parasitic nest position, compares the position with the position before levy migration is executed, and selects a parasitic nest with a better position to be reserved to the next generation.
And 4, step 4: performing a parasitic nest discovery operation to generate random numbersIf, ifAnd disturbing the parasitic nest to generate a new parasitic nest, comparing the position with the corresponding position before disturbance, and reserving the parasitic nest with a better position.
And 5: performing parasitic nest mutation operations usingAnd (5) carrying out mutation operation on the parasitic nests.
Step 6: dynamically adjusting discovery probability, usingDynamically adjusting parasitic nest discovery probabilityIn the middle of the termRepresenting the probability of finding that the t-th iteration is carried out;、maximum discovery probability and minimum discovery probability respectively;is the maximum iteration number;is the current iteration number. This method ensures that at the early stage of the algorithm, a greater rate of change of position is required due to the individual being further from the optimum; at the later stages of the algorithm, a smaller rate of change of position is required, since most of the solution is concentrated around the optimal position.
And 7: keeping the optimal position and the fitness of each search, judging whether the search result meets the requirements, and if so, judging whether the search result meets the requirementsGo to Step 8, otherwise go to Step 3.
And 8: and outputting the optimal parasitic nest position and the corresponding fitness value thereof to obtain the optimal logistics distribution vehicle path scheme.
Further, the calculation steps of the fitness value are as follows:
step 1: determine the vehicles required for any logistics distribution scheme x by roadindex = unique (x), by [ row, car _ index [ ]]= size (roadindex) determine vehicle number car _ index, set fitness valueThe optimal distribution path total _ road of the current scheme;
step 2: determining a customer point served by each vehicle by subRout = find (Zx = = roadindex (i)), and forming a distribution path subRout formed by the current vehicle according to the size of a distribution number;
and step 3: optimizing each distribution path subRout by adopting a 2-opt algorithm so as to obtain an optimal path road for single vehicle distribution, and recording the sequence of the vehicle distribution paths: total _ road (subrout) = road;
and 4, step 4: calculating the length of the path which is passed by the current vehicle optimal path road;
and 5: calculating the sum _ q of the customer demands of the current vehicle optimal path road;
step 6: punishing the load exceeding the maximum load of the current vehicle by adopting a sectional penalty function method, and setting the punishment amount as follows: fp;
and step 9: outputting fitness value of scheme xAnd the current scheme optimal delivery path total _ road.
The improved cuckoo algorithm is adopted for planning vehicle dispatching for grouping, and then 2-opt algorithm is adopted for optimizing each line so as to obtain an optimal distribution line in each group.
The method provides a feasible and effective solution for solving the logistics vehicle distribution optimization problem, and enriches the method for solving the logistics distribution path optimization problem.
Drawings
FIG. 1 is a flow chart of the algorithm for solving the logistics distribution vehicle routing problem;
fig. 2 is a flow chart for solving the fitness value for each logistics distribution scheme.
Detailed Description
The technical solutions of the present invention are further explained in detail by the drawings and the specific examples, but the scope of the present invention is not limited thereto.
The implementation method for solving the logistics distribution vehicle path problem based on the improved cuckoo algorithm is disclosed by the embodiment. Referring to fig. 1, fig. 1 is a schematic flow chart of the algorithm, which includes the following steps:
① set the parameters of the algorithm, set the size of the group, the number of service clients nd, the number of vehicles required to solve the problem at present, set the load of the vehicles carraycaran. the probability of finding a bird egg in a nest in the basic cuckoo algorithm is usually set to a constant value, so that when the algorithm executes the operation of finding a nest in the basic cuckoo algorithm, the nest in either the superior or inferior position is replaced with the same probability, if the cluster is in the superior or inferior positionIf the value is set too small, the poor solution convergence in the optimization process is slow, if it is too smallThe value is set too large, and the solution of a better position is difficult to converge to the optimal solution, so the invention adopts a dynamic discovery mechanism:in the middle of the termAlgorithm search spatial Range [1, number ]]The number of times of population iteration Max _ iter, and the iteration counter t = 1;
② initialize the population and calculate the fitness value the cost of transportation between the distribution center and each customer demand point is calculated, in this example, by targeting the distance between the points as pseudo code:
for i=1:nd
for j=1:nd
cost(i,j)=sqrt((xy(i,1)-xy(j,1))^2+(xy(i,2)-xy(j,2))^2);
end
end
where xy is the coordinates of each point.
Randomly generating size nd dimension search space range of [1, number ] by using uniform distribution function]The position of the ith parasitic nest is recorded asAnd rounding each nest with a ceil function (i.e.) (ii) a To ensure that each customer site is served by a vehicle and calculate its corresponding fitness value;
③ carry out levy walk operations of parasitic nestsAnd a new parasitic nest position is generated in an updating mode, and because Levy distribution integral is difficult, equivalent calculation is realized by adopting a Mantegana algorithm, and the equivalent calculation is compared with the position before Levy migration is executed, so that the parasitic nest with a better position is selected and reserved to the next generation.
④ perform a nest discovery operationIf, ifAnd disturbing the parasitic nest to generate a new parasitic nest, comparing the position with the corresponding position before disturbance, and reserving the parasitic nest with a better position.
⑤ executing the parasitic nest mutation operationTo hostThe living nests are subjected to mutation operation.
⑥ dynamically adjusting the probability of discoveryDynamically adjusting parasitic nest discovery probabilityIn the middle of the termRepresenting the discovery probability of performing the Nth iter iteration;、maximum discovery probability and minimum discovery probability respectively;is the maximum iteration number;is the current iteration number. This method ensures that at the early stage of the algorithm, a greater rate of change of position is required due to the individual being further from the optimum; at the later stages of the algorithm, a smaller rate of change of position is required, since most of the solution is concentrated around the optimal position.
⑦ keeping the optimal position and its fitness for each search, determining whether the search result meets the requirement, and if so, determining whether the result meets the requirementGo to ⑧, otherwise go to ③.
⑧ outputs the optimal parasitic nest position and the corresponding fitness value to obtain the optimal logistics distribution vehicle path scheme.
The fitness value involved in the above solving stepReferring to fig. 2, fig. 2 is a schematic flowchart of the fitness value solving, including the following steps:
step 1: determining vehicles required by any logistics distribution scheme x by roadindex = unique (x), and setting a fitness valueThe optimal distribution path total _ road of the current scheme;
step 2: determining a customer point served by the ith vehicle by subRout = find (Zx = = roadindex (i)), and forming a distribution path subRout formed by the current vehicle according to the size of a distribution number;
and step 3: optimizing each distribution path subRout by adopting a 2-opt algorithm so as to obtain an optimal path road for the ith vehicle distribution, and recording the sequence of the vehicle distribution paths: total _ road (subrout) = road;
and 4, step 4: calculating the length of the path which is passed by the current vehicle optimal path road;
and 5: calculating the sum _ q of the customer demands of the current vehicle optimal path road;
step 6: punishing the load exceeding the maximum load of the current vehicle by adopting a sectional penalty function method, and setting the punishment amount as follows: fp, the specific calculation method is as follows:
the violation degree of the solution on the constraint condition is represented and calculated by the following formula:
determining different punishment degrees for the multi-segment mapping function according to different violation degrees, wherein the punishment mode is calculated by adopting the following formula:
and 8: judging whether the vehicle is the last vehicle, if so, turning to the step 2, otherwise, turning to the step 9;
Claims (2)
1. A logistics distribution vehicle path optimization method based on an improved cuckoo algorithm is characterized by comprising the following steps of:
solving:
step 1: setting parameters of an algorithm, setting a population size, the number nd of service customers, the number of vehicles required for solving the problem, the load CarrayCarCan of the vehicles, and the discovery probability of the bird eggs in the parasitic nestsAlgorithm search spatial Range [1, number ]]The number of times of population iteration Max _ iter, and the iteration counter N _ iter = 1;
step 2: initializing population and calculating fitness value, and randomly generating size nd dimension search space range of [1, number ] by adopting uniform distribution function]The position of the ith parasitic nest is recorded asAnd rounding each nest with a ceil function, i.e.To ensure that each customer site is served by a vehicle and calculate its corresponding fitness value
And step 3: performing levy walk operations of parasitic nests, usingThe updating mode generates a new parasitic nest position, compares the position with the position before levy walk execution, and selects the parasitic nest with better position to be reserved to the next generation, whereinWhich represents a step-size control factor,as a function of column dimension distributionThe parameters of (1);
and 4, step 4: performing a parasitic nest discovery operation to generate random numbersDisturbing the parasitic nest to generate a new parasitic nest, comparing the new parasitic nest with the corresponding position before disturbance, and reserving the parasitic nest with a better position;
and 5: performing parasitic nest mutation operations usingCarrying out mutation operation on the parasitic nests;
step 6: dynamically adjusting discovery probability, usingDynamically adjusting parasitic nest discovery probabilityIn the middle of the termRepresenting the probability of finding that the t-th iteration is carried out;maximum discovery probability and minimum discovery probability respectively;is the maximum iteration number;the current iteration number is; this method ensures that at the early stage of the algorithm, a greater rate of change of position is required due to the individual being further from the optimum; at the later stages of the algorithm, a smaller rate of change of position is required, since most of the solution is concentrated around the optimal position;
and 7: keeping the optimal position and the fitness of each search, judging whether the search result meets the requirements, and if so, judging whether the search result meets the requirementsThen, the process goes to Step 8,otherwise, turning to Step 3;
and 8: and outputting the optimal parasitic nest position and the corresponding fitness value thereof to obtain the optimal logistics distribution vehicle path scheme.
2. The logistics distribution vehicle path optimization method based on the improved cuckoo algorithm as claimed in claim 1, wherein the calculation steps of the involved fitness values are as follows:
step 1: determine the vehicles required for any logistics distribution scheme x by roadindex = unique (x), by [ row, car _ index [ ]]= size (roadindex) determine vehicle number car _ index, set fitness valueThe optimal distribution path total _ road of the current scheme;
step 2: determining a customer point served by each vehicle by subRout = find (Zx = = roadindex (i)), and forming a distribution path subRout formed by the current vehicle according to the size of a distribution number;
and step 3: optimizing each distribution path subRout by adopting a 2-opt algorithm so as to obtain an optimal path road for single vehicle distribution, and recording the sequence of the vehicle distribution paths: total _ road (subrout) = road;
and 4, step 4: calculating the length of the path which is passed by the current vehicle optimal path road;
and 5: calculating the sum _ q of the customer demands of the current vehicle optimal path road;
step 6: punishing the load exceeding the maximum load of the current vehicle by adopting a sectional penalty function method, and setting the punishment amount as follows: fp;
and 7: calculating the fitness value of the current optimal vehicle path road as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510521483.4A CN105260785B (en) | 2015-08-24 | 2015-08-24 | Logistics distribution vehicle path optimization method based on improved cuckoo algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510521483.4A CN105260785B (en) | 2015-08-24 | 2015-08-24 | Logistics distribution vehicle path optimization method based on improved cuckoo algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260785A CN105260785A (en) | 2016-01-20 |
CN105260785B true CN105260785B (en) | 2020-06-23 |
Family
ID=55100464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510521483.4A Active CN105260785B (en) | 2015-08-24 | 2015-08-24 | Logistics distribution vehicle path optimization method based on improved cuckoo algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260785B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105975342B (en) * | 2016-04-29 | 2019-02-15 | 广东工业大学 | Based on the cloud computing method for scheduling task and system for improving cuckoo searching algorithm |
CN106127295A (en) * | 2016-06-21 | 2016-11-16 | 湘潭大学 | A kind of Optimal Design of Pressure Vessel method based on self adaptation cuckoo Yu fireworks hybrid algorithm |
CN106651086B (en) * | 2016-09-07 | 2020-06-26 | 河南科技学院 | Automatic stereoscopic warehouse scheduling method considering assembly process |
CN107909228B (en) * | 2017-12-23 | 2021-10-29 | 深圳大学 | Dynamic vehicle goods receiving and dispatching path planning method and device based on modular factor calculation |
CN108388250B (en) * | 2018-03-30 | 2021-03-05 | 哈尔滨工程大学 | Water surface unmanned ship path planning method based on self-adaptive cuckoo search algorithm |
CN108921338B (en) * | 2018-06-20 | 2022-04-19 | 广东工业大学 | Multi-vehicle workshop logistics transportation scheduling method |
CN110334853A (en) * | 2019-06-10 | 2019-10-15 | 福建工程学院 | A kind of imitative nature body optimization method of logistics distribution center Warehouse Location |
CN111178596A (en) * | 2019-12-12 | 2020-05-19 | 浙江浙大网新国际软件技术服务有限公司 | Logistics distribution route planning method and device and storage medium |
CN111178730A (en) * | 2019-12-24 | 2020-05-19 | 中国航空工业集团公司西安飞机设计研究所 | Method and device for planning supply of oiling machine |
CN112232602A (en) * | 2020-11-19 | 2021-01-15 | 湘潭大学 | Logistics distribution path optimization method and system for large-scale nodes |
CN116502473B (en) * | 2023-06-27 | 2024-01-12 | 中科航迈数控软件(深圳)有限公司 | Wire harness electromagnetic compatibility optimization method and device, electronic equipment and storage medium |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10606962B2 (en) * | 2011-09-27 | 2020-03-31 | Autodesk, Inc. | Horizontal optimization of transport alignments |
-
2015
- 2015-08-24 CN CN201510521483.4A patent/CN105260785B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105260785A (en) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105260785B (en) | Logistics distribution vehicle path optimization method based on improved cuckoo algorithm | |
Xu et al. | Service offloading with deep Q-network for digital twinning-empowered internet of vehicles in edge computing | |
CN110705742A (en) | Logistics distribution method based on improved ant colony algorithm | |
Wang et al. | Ant colony optimization with an improved pheromone model for solving MTSP with capacity and time window constraint | |
CN103971530A (en) | Real-time route planning method | |
CN105302858B (en) | A kind of the cross-node enquiring and optimizing method and system of distributed data base system | |
CN110516871B (en) | Dynamic vehicle path optimization method based on fuzzy rolling time domain control strategy | |
JP2010134863A (en) | Control input determination means of control object | |
CN111445094B (en) | Express vehicle path optimization method and system based on time requirement | |
CN111885493B (en) | Micro-cloud deployment method based on improved cuckoo search algorithm | |
CN116187896A (en) | Green vehicle path problem solving method, device, computer equipment and medium | |
CN111476389A (en) | Method and device for pre-estimating order receiving waiting time | |
CN110097218A (en) | Unmanned commodity distribution method and system under changing environment when a kind of | |
CN110930696B (en) | AI navigation-based intelligent city traffic management operation method and system | |
CN109840625B (en) | Courier group path navigation method | |
Zhao et al. | Adaptive Swarm Intelligent Offloading Based on Digital Twin-assisted Prediction in VEC | |
CN117252318B (en) | Intelligent networking automobile group machine collaborative carpooling scheduling method and system | |
US11237008B2 (en) | System and method for controlling vehicular pollution concentration and providing maximum traffic flow throughput | |
CN114882720B (en) | Intelligent road network pushing method, intelligent road network pushing device and vehicle | |
JP6045388B2 (en) | Operation management device, operation management method, and program | |
CN115062868B (en) | Pre-polymerization type vehicle distribution path planning method and device | |
CN116777171A (en) | Point-to-point carpooling dynamic scheduling method and device for network taxi taking | |
CN108389003B (en) | Method and device for scheduling service tasks under remote health monitoring line | |
Zhang et al. | Research on taxi driver strategy game evolution with carpooling detour | |
Lai et al. | Utility-based matching of vehicles and hybrid requests on rider demand responsive systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |