WO2021142917A1 - Multi-depot vehicle routing method, apparatus, computer device and storage medium - Google Patents

Multi-depot vehicle routing method, apparatus, computer device and storage medium Download PDF

Info

Publication number
WO2021142917A1
WO2021142917A1 PCT/CN2020/079881 CN2020079881W WO2021142917A1 WO 2021142917 A1 WO2021142917 A1 WO 2021142917A1 CN 2020079881 W CN2020079881 W CN 2020079881W WO 2021142917 A1 WO2021142917 A1 WO 2021142917A1
Authority
WO
WIPO (PCT)
Prior art keywords
path planning
target
individual
population
current
Prior art date
Application number
PCT/CN2020/079881
Other languages
French (fr)
Chinese (zh)
Inventor
于琪嫄
刘松柏
林秋镇
陈剑勇
Original Assignee
深圳大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 深圳大学 filed Critical 深圳大学
Publication of WO2021142917A1 publication Critical patent/WO2021142917A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083Shipping
    • G06Q10/0835Relationships between shipper or supplier and carriers
    • G06Q10/08355Routing methods
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Definitions

  • This application relates to the technical field of route planning, and in particular to a method, device, computer equipment, and storage medium for multi-yard vehicle route planning.
  • the embodiments of the present application provide a method, device, computer equipment, and storage medium for vehicle path planning in multiple depots, aiming to solve the problem that the prior art cannot quickly and accurately perform optimization when multiple distribution centers are distributed in different areas.
  • an embodiment of the present application provides a multi-park vehicle path planning method, which includes:
  • the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and each corresponding to the current number of users.
  • the user s cargo capacity and current user location information for each user corresponding to the current number of users;
  • an embodiment of the present application provides a multi-park vehicle path planning device, which includes a unit for executing the multi-park vehicle path planning method described in the first aspect.
  • an embodiment of the present application provides a computer device, which includes a memory, a processor, and a computer program stored on the memory and running on the processor, and the processor executes the computer
  • the program implements the multi-depot vehicle path planning method described in the first aspect.
  • an embodiment of the present application also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the processor executes the above-mentioned first On the one hand, the multi-depot vehicle path planning method.
  • FIG. 1 is a schematic diagram of an application scenario of a multi-park vehicle path planning method provided by an embodiment of the application;
  • FIG. 2 is a schematic flowchart of a method for multi-parking vehicle path planning provided by an embodiment of the application
  • FIG. 3 is a schematic diagram of a sub-process of a multi-park vehicle path planning method provided by an embodiment of the application;
  • FIG. 4 is a schematic diagram of another sub-flow of the method for multi-depot vehicle path planning provided by an embodiment of the application;
  • Fig. 5 is a schematic block diagram of a multi-park vehicle path planning device provided by an embodiment of the application.
  • Fig. 6 is a schematic block diagram of a computer device provided by an embodiment of the application.
  • FIG. 1 is a schematic diagram of an application scenario of a multi-park vehicle path planning method provided by an embodiment of the application
  • FIG. 2 is a schematic flowchart of a multi-park vehicle path planning method provided by an embodiment of the application.
  • the vehicle path planning method is applied to the server, and the method is executed by application software installed in the server.
  • the method includes steps S110 to S140.
  • S110 Determine whether a path planning request sent by the client is received.
  • the first is the client.
  • the client can be understood as a user terminal.
  • the user terminal can be a smart phone, tablet computer, notebook computer, desktop computer, personal digital assistant, wearable device and other electronic devices with communication functions.
  • the user terminal sends a path plan Request to the server.
  • the second is the server.
  • the server receives the path planning request sent by the client, based on the input data and constraint conditions corresponding to the path planning request, and calls the pre-stored vehicle path planning multi-objective optimization model to perform the evolutionary solution of super multi-objectives. Obtain the path optimal solution set.
  • the optimal solution set of the path is obtained from the server and sent to the client.
  • the server detects whether the path planning request sent by the client is received, and when the server receives the path planning request sent by the client, the subsequent step S120 is executed. When the server does not receive the path planning request sent by the client Upon request, step S110 is executed again after waiting for the preset delay time.
  • the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and information corresponding to the current number of users.
  • the current user location information of each user corresponding to the cargo capacity of each user and the number of current users.
  • the server receives the path planning request sent by the client, it obtains the input data and constraint conditions corresponding to the path planning request. Since the vehicle path planning multi-objective optimization model has been pre-stored in the server, it can be solved subsequently according to the input data and constraint conditions, so as to obtain the path optimal solution set.
  • the vehicle path planning multi-objective optimization model stored in the server is a multi-depot multi-vehicle path planning multi-objective optimization model.
  • the optimization objectives are all Try to achieve a satisfactory path scheduling plan.
  • the vehicle path planning multi-objective optimization model includes five optimization objective functions, which are respectively denoted as:
  • the vehicle path planning multi-objective optimization model is preset with R parking lots, each parking lot has K vehicles, the total cargo capacity of each vehicle is Q and the maximum total service time is T; R, S, K and The value of Q is a positive integer; the current number of users included in the input data corresponding to the path planning request is denoted as P; the P customer nodes corresponding to the current number of users P are denoted as node 1 to node P, and R parking lots The corresponding parking lot nodes are respectively marked as node P+1 to node P+R;
  • K r represents the r-th actual number of yard vehicles
  • RK k represents a vehicle ID of the r-yard
  • d ij represents the distance between node i to the j-th node, Represents the travel time of the k numbered vehicle of the r-th parking lot from the i-th node to the j-th node
  • s i represents the service time corresponding to the i-th node
  • s j represents the service time corresponding to the j-th node
  • p i represents the i-th node
  • the cargo capacity of the i-th customer corresponding to the node Represents the path access state of the k numbered vehicle of the rth parking lot moving from the ith node to the jth node, It means that the k-numbered vehicle of the r-th parking lot moves from the i-th node to the customer of the j-th node and is served.
  • the multi-park vehicle path planning problem can be defined as: assuming that there are R parking lots, each parking lot has K vehicles with a total cargo capacity of Q and a maximum total service time (path travel time plus customer service time) of T, and there is P A customer needs delivery, the cargo capacity of the i-th customer p i ⁇ Q, and the service time s i ⁇ T of the customer p i.
  • Each customer can be served by any vehicle but can only be served once, each vehicle can serve multiple users and can be required to return to the departure depot when the service is over.
  • a candidate solution x of the vehicle path planning multi-objective optimization model it refers to satisfying the above 5 optimization objective functions (that is, satisfying minf 1 (x), minf 2 (x), minf 3 (x), minf 4 (x) ), a path of minf 5 (x)), X represents a set containing multiple candidate solutions, and the path optimal solution set X is optimal composed of multiple path optimal solutions.
  • the model is a high-dimensional optimization model that combines the above five objective functions, namely input data and constraint conditions, to obtain the path optimal solution set X
  • the proposed optimization goals and constraints can be met to the greatest extent. More specifically, the most appropriate number of vehicle dispatches, the smaller the total delivery duration of vehicles and the total travel distance of vehicles, the appropriate maximum capacity difference of single vehicles, and the Smaller waiting time for customers.
  • the constraint conditions corresponding to the path planning request are as follows:
  • st is called subject to and means to be restricted (the expression of general constraints starts with st), for each i ⁇ 1,2,...,P ⁇ : Indicates that the values of i in the set ⁇ 1,2, whil,P ⁇ are such that Through the data and the constraint conditions, the vehicle path planning multi-objective optimization model can be solved.
  • the step S130 includes:
  • mapping set Map each individual in the normalized archive set to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to a mapping on the hypersurface Points to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max ;
  • an initial multi-target population is randomly generated under the restriction of constraint conditions.
  • the initial multi-target population is the first-generation multi-target population.
  • each target value of the ideal individual indicates that the initial multi-target population is in the corresponding objective function
  • the minimum target value of the ideal individual that is, the path output solution corresponding to the ideal individual is substituted into f 1 (x) to f 5 (x), and all correspond to the minimum target value; each target value of the worst individual represents the current population in the corresponding target function
  • the maximum target value that is, the path output solution corresponding to the worst individual corresponds to the maximum target value after substituting f 1 (x) to f 5 (x).
  • the initial multi-target population can be simulated binary crossover and polynomial mutation to obtain a subpopulation with the same total number of individuals as the initial multi-target population.
  • the generation of subpopulation is to randomly select two individuals from the current initial population to simulate binary crossover until N new individuals are crossed. Then, the N new individuals are mutated according to the mutation probability and polynomial mutation to get the updated N individuals, the updated N individuals form a subpopulation).
  • two individuals are randomly selected from the initial multi-target population to perform binary crossover in sequence until N cross-processed new individuals are generated, and the N cross-processed new individuals are subjected to polynomial mutation, and the new individual after polynomial mutation Form subpopulations.
  • N new individuals after crossover processing are obtained.
  • the process of randomly selecting two individuals for binary crossover is similar to an iterative process. Until the number of new individuals reaches the population size N, the process of multiple binary crossovers is stopped.
  • binary crossover and polynomial mutation are both conventional processing procedures, and will not be repeated here.
  • the initial multi-target population and the sub-population are combined to obtain a mixed population, and the total number of individuals included in the mixed population is twice the population size N.
  • non-dominated sorting can be performed on the individuals in the mixed population, thereby obtaining non-dominated solution sets and multi-layer solution sets.
  • the non-dominated solution set corresponding to the mixed population can be obtained through a non-dominated solution (also called Pareto solution) acquisition method.
  • a non-dominated solution also called Pareto solution
  • the definition of Pareto solution is to assume that for any two solutions S1 and S2, S1 is better than or the same as S2 for all targets, and there is at least one target, and the corresponding target value of S1 on this target is better than S2.
  • the corresponding target value on the target is called S1 dominates S2.
  • S1 is called the non-dominated solution (undominated solution), also called the Pareto solution (ie Pareto solution).
  • the obtained non-dominated solution set is denoted as Q 1 .
  • the multi-layer solution set is obtained.
  • the multi-layer solution set includes multiple solution set subsets and is respectively denoted as Q 2 to Q L , where Q 1 to The union of Q L is the mixed population, and the intersection of any two sets from Q 1 to Q L is an empty set, Q 1 ⁇ Q 2 ⁇ Q 3 ⁇ « ⁇ Q L ; where " ⁇ " indicates a dominance relationship , Q i ⁇ Q j means that Q j is dominated by the solution in Q i , and the relationship is transitive.
  • Q 1 ⁇ Q 2 means that for f 1 (x) to f 5 (x), each of Q2 The solutions are all dominated by at least one solution in Q1, and the relationship is transitive, that is, each solution in Q3 is dominated by at least one solution in Q1 or Q2, and the others are in turn.
  • solutions that exceed the population size N need to be selected at this time to form an archive set.
  • the selection mode as follows: first of all individuals in the selected Q 1, determines whether the total number of individuals beyond the current population size N, if the current total number of individuals of the population size does not exceed N, continue Q 2 in the selection of all of the individuals, the total number of individual current plus the number of Q 2 in the subject to update the current as the total number of individuals, and then determines whether the total number of individuals beyond the current The population size N, until the total number of individuals in Q 1 to Q a obtained exceeds the population size N (where a is a positive integer greater than 1 and not exceeding L) to form an archive set.
  • each individual in the archive set may be normalized according to the ideal individual and the worst individual to obtain a normalized archive set; wherein, the normalized archive set and the The normalized individual set corresponding to the non-dominated solution set is recorded as the normalized non-dominated solution set. That is, each individual in the archive set is normalized, and a normalized individual corresponding to each individual is obtained, thereby forming a normalized non-dominated solution set.
  • the purpose of the above normalization processing is to eliminate the difference between different dimensions, so as to facilitate subsequent data processing.
  • step S1308 includes:
  • NA m represents in the archive set the individual over the individual normalization
  • a represents the worst of the worst individual subject
  • a represents the individual over the individual m-th corresponding to a m.
  • a normalized archive set can be obtained, thereby eliminating the difference between different dimensions.
  • the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front can be estimated according to the normalized individuals included therein. Estimate the shape of the population in the objective function space based on all non-dominated solutions, which will help to discover the characteristics of multi-objective optimization tasks. Multi-objective optimization is usually divided into three categories, convex optimization, concave optimization, and linear optimization. In this step, the curvature values of three types of hypersurfaces will be obtained.
  • step S1309 includes:
  • each normalized non-dominated individual in the normalized non-dominated solution set Remove the normalized non-dominated individuals that are not in the standardized target space (the standardized target space is the target space formed by f 1 (x) to f 5 (x) equal to 1), and get the normalized non-dominated solution after screening set. After removing the normalized non-dominated individuals that are not in the standardized target space, the interference of these individuals is effectively eliminated, which is beneficial to the subsequent population evolution.
  • Cur represents the current curvature corresponding to the normalized non-dominated solution set.
  • the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front are obtained through the processing process of steps S13091-S13098, in order to facilitate the subsequent clustering operation, and the individuals used in the prediction process are all non-dominated individual.
  • a negative value of the distance D i means that the individual is below the target hyperplane
  • a positive value means that the individual is above the target hyperplane
  • the distance D i is 0 Indicates that the individual is on the target hyperplane.
  • the current curvature is adjusted according to the coefficient of variation cv to obtain the adjusted curvature.
  • the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front can be correspondingly obtained according to the adjusted curvature.
  • the hypersurface corresponding to the shape of the Pareto front will help to discover the characteristics of the multi-objective optimization task.
  • each individual in the normalized archive set is mapped to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to a mapping on the hypersurface Points to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max .
  • the fitness value of each mapping point needs to be obtained, specifically Obtain the fitness value of each mapping point, f i (x l ) represents the i-th target value corresponding to the l-th mapping point in the mapping set, and the value range of l is [1, L max ].
  • the fitness value of each mapping point represents the sum of all the target values of the individual corresponding mapping points. Since it tends to find a well-converged solution set, the fitness value is used as a measure of convergence.
  • a mapping point is selected from each cluster cluster in the clustering result to form a target mapping point set, so that the number of mapping points in the target mapping point set is equal to the population size N, and the mapping
  • the individual set corresponding to the point set is the population of this evolution, and it is also the parent population of the next generation. Through this selection strategy, it can ensure that the next generation population has good convergence and diversity.
  • the current iteration algebra is iterated multiple times to reach the maximum iteration algebra, and the current multi-target population output is used as a path optimal solution set.
  • each individual in the path optimal solution set is the optimal solution corresponding to the vehicle path planning multi-objective optimization model obtained according to the input data and constraint conditions.
  • the specific encoding method of each path optimal solution in the path optimal solution set (that is, the method that is finally displayed to the user) is a user-defined encoding format and saved in the server, which is not limited this time Specific encoding method.
  • the client after obtaining the path optimal solution set in the server, it can be sent to the client. Therefore, the client can assist the delivery process after determining the delivery route according to the optimal solution set of the route.
  • This method realizes that the convergence and diversity of the population are fully considered in the process of solving super-multi-objective evolution, and the effective prediction of the shape of the population is realized, and the multi-objective of vehicle path planning is obtained quickly and accurately based on input data and constraints.
  • the path optimal solution set of the optimization model is obtained quickly and accurately based on input data and constraints.
  • the embodiment of the present application also provides a multi-park vehicle path planning device, which is used to execute any embodiment of the foregoing multi-park vehicle path planning method.
  • FIG. 5 is a schematic block diagram of a multi-park vehicle path planning apparatus provided by an embodiment of the present application.
  • the multi-parking vehicle path planning device 100 may be configured in a server.
  • the multi-depot vehicle path planning device 100 includes a path planning request detection unit 110, a data condition acquisition unit 120, a path optimal solution set acquisition unit 130, and an optimal solution set transmission unit 140.
  • the path planning request detection unit 110 is used to determine whether the path planning request sent by the client is received.
  • the data condition obtaining unit 120 is configured to, if a path planning request sent by the client is received, obtain input data and constraint conditions corresponding to the path planning request; wherein the input data corresponding to the path planning request includes the current number of users , The cargo capacity of each user corresponding to the current number of users, and the current user location information of each user corresponding to the current number of users.
  • the path optimal solution set acquisition unit 130 is configured to call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and according to the constraint conditions and the control
  • the vehicle path planning multi-objective optimization model performs the evolutionary solution of super multi-objectives, and obtains the optimal solution set of the path.
  • the optimal solution set sending unit 140 is configured to send the path optimal solution set to the client.
  • the path optimal solution set obtaining unit 130 includes:
  • the initial multi-objective population generating unit is used to randomly generate an initial multi-objective population according to the constraint conditions; wherein the initial multi-objective population includes multiple individuals, and each individual corresponds to one of the vehicle path planning multi-objective optimization models Path output solution, the total number of multiple individuals included in the initial multi-target population is recorded as the population size N;
  • the first judging unit of the current iteration algebra is used to obtain the current iteration algebra, and judge whether the current iteration algebra reaches the preset maximum iteration algebra;
  • a target individual obtaining unit configured to obtain an ideal individual and a worst individual in the initial multi-target population if the current iteration algebra does not reach the maximum iteration algebra; wherein the ideal individual is input to the vehicle path planning
  • the target value obtained by the multi-objective optimization model is the smallest target value among the target values corresponding to each individual in the initial multi-objective population, and the worst individual is input to the vehicle path planning multi-objective optimization model and the target value obtained is the initial multi-objective The largest target value among the target values corresponding to each individual in the population;
  • the individual crossover mutation unit is used to simulate binary crossover and polynomial mutation on the initial multi-target population to obtain a subpopulation with the same total number of individuals as the initial multi-target population;
  • a mixed population obtaining unit configured to merge the initial multi-target population and the sub-population to obtain a mixed population
  • the non-dominated solution set obtaining unit is used to perform non-dominated sorting of the individuals in the mixed population to obtain a non-dominated solution set and a multi-layer solution set; wherein the non-dominated solution set is denoted as Q 1 , and the multi-layer
  • the solution set includes multiple solution set subsets and are respectively denoted as Q 2 to Q L , where the union of Q 1 to Q L is the mixed population, and the intersection of any two sets from Q 1 to Q L is an empty set, Q 1 ⁇ Q 2 ⁇ Q 3 ⁇ « ⁇ Q L ;
  • the archive set acquisition unit is used to sequentially merge multiple solution set subsets in the non-dominated solution set and the multi-layer solution set to obtain multiple sets until the total number of individuals exceeds the population size N to form an archive gather;
  • the normalization processing unit is configured to normalize each individual in the archive set according to the ideal individual and the worst individual to obtain a normalized archive set; wherein, the normalized archive set The set of normalized individuals corresponding to the non-dominated solution set is recorded as the normalized non-dominated solution set;
  • a hypersurface acquiring unit configured to estimate and acquire the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the normalized non-dominated solution set;
  • the individual mapping unit is used to map each individual in the normalized archive set to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to the hypersurface A mapping point above to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max ;
  • the fitness value obtaining unit is used to call a pre-stored target point fitness value algorithm to obtain the fitness value corresponding to each mapping point in the mapping set;
  • the mapping point clustering unit is used to cluster the corresponding Euclidean distance in the hyperplane between each mapping point in the mapping set and the population size N to obtain a clustering result; wherein, the The total number of clusters in the clustering result is equal to the population size N;
  • the target mapping point set acquiring unit is configured to select a mapping point in each cluster cluster of the clustering result to form a target mapping point set;
  • the initial multi-target population update unit is used to obtain the individual corresponding to each target mapping point in the target mapping point set to form a current multi-target population, and use the current multi-target population as the initial multi-target population;
  • the current iteration algebra second judging unit configured to add one to the current iteration algebra as the current iteration algebra, and return to execute the step of judging whether the current iteration algebra reaches the preset maximum iteration algebra;
  • the path optimal solution set output unit is configured to output the current multi-target population as the path optimal solution set if the current iteration algebra reaches the maximum iteration algebra.
  • the individual cross mutation unit is also used for:
  • the normalization processing unit is further used for:
  • NA m represents in the archive set the individual over the individual normalization
  • a represents the worst of the worst individual subject
  • a represents the individual over the individual m-th corresponding to a m.
  • the hypersurface acquiring unit includes:
  • the non-dominated individuals screening unit is used to remove the normalized non-dominated individuals that are not in the standardized target space among the normalized non-dominated individuals in the normalized non-dominated solution set to obtain the normalized non-dominated individuals after screening Solution set; wherein each target value corresponding to each normalized non-dominated individual located in the standardized target space does not exceed 1;
  • the non-dominated individual acquisition unit is used to acquire each non-dominated individual in the normalized non-dominated solution set after the screening, which is respectively recorded as B 1 to B n ; where the value of n is the same as the normalized non-dominated after screening.
  • the total number of normalized non-dominated individuals in the dominating solution set is the same;
  • Hyperplane from the parameter acquisition unit configured to obtain hyperplane distance D 1 to D n corresponding hyperplane distance from average and standard deviation hyperplane; wherein hyperplane distance D 1 to D n corresponding to the average distance referred hyperplane Is D avg , the standard deviation of the hyperplane distance corresponding to the hyperplane distance D 1 to D n is denoted as D std ;
  • the coefficient of variation obtaining unit is configured to perform a norm operation according to the quotient of the average value of the hyperplane distance and the standard deviation of the hyperplane distance to obtain the corresponding coefficient of variation; wherein the coefficient of variation is denoted as cv;
  • the current curvature acquisition unit is configured to acquire a first hypersurface corresponding to a curvature of 2 of the target hyperplane, and a second hypersurface corresponding to a curvature of 0.5 of the target hyperplane, according to the target hyperplane and the first hypersurface , The second hypersurface and the preset curvature determination strategy, to obtain the current curvature corresponding to the normalized non-dominated solution set; wherein, the curvature determination strategy is d(2.0) represents the distance from the peak point of the first hypersurface to the target hyperplane, and d(0.5) represents the distance from the peak point of the second hypersurface to the target hyperplane;
  • the curvature adjustment unit is configured to adjust the current curvature according to the coefficient of variation to obtain an adjusted curvature; wherein, if the coefficient of variation is less than 0.1, adjust the value of the current curvature to 1, so that the adjusted curvature The value of the curvature is 1; if the coefficient of variation is greater than or equal to 0.1, the value of the current curvature is kept unchanged, so that the adjusted curvature is equal to the current curvature;
  • the hypersurface generating unit is used to obtain the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the adjusted curvature.
  • the device realizes that the convergence and diversity of the population are fully considered in the process of solving super-multi-objective evolution, and the effective prediction of the shape of the population is realized, and the multi-objective of vehicle path planning is obtained quickly and accurately based on input data and constraint conditions.
  • the path optimal solution set of the optimization model is obtained quickly and accurately based on input data and constraint conditions.
  • the above-mentioned multi-yard vehicle path planning apparatus may be implemented in the form of a computer program, and the computer program may run on a computer device as shown in FIG. 6.
  • FIG. 6 is a schematic block diagram of a computer device according to an embodiment of the present application.
  • the computer device 500 is a server, and the server may be an independent server or a server cluster composed of multiple servers.
  • the computer device 500 includes a processor 502, a memory, and a network interface 505 connected through a system bus 501, where the memory may include a non-volatile storage medium 503 and an internal memory 504.
  • the non-volatile storage medium 503 can store an operating system 5031 and a computer program 5032.
  • the processor 502 can execute the multi-park vehicle path planning method.
  • the processor 502 is used to provide calculation and control capabilities, and support the operation of the entire computer device 500.
  • the internal memory 504 provides an environment for the operation of the computer program 5032 in the non-volatile storage medium 503.
  • the processor 502 can make the processor 502 execute the multi-park vehicle path planning method.
  • the network interface 505 is used for network communication, such as providing data information transmission.
  • the structure shown in FIG. 6 is only a block diagram of part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device 500 to which the solution of the present application is applied.
  • the specific computer device 500 may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.
  • the processor 502 is configured to run a computer program 5032 stored in a memory to implement the multi-park vehicle path planning method disclosed in the embodiment of the present application.
  • the embodiment of the computer device shown in FIG. 6 does not constitute a limitation on the specific configuration of the computer device.
  • the computer device may include more or less components than those shown in the figure. Or combine certain components, or different component arrangements.
  • the computer device may only include a memory and a processor. In such an embodiment, the structure and function of the memory and the processor are consistent with the embodiment shown in FIG. 6 and will not be repeated here.
  • the processor 502 may be a central processing unit (Central Processing Unit, CPU), and the processor 502 may also be other general-purpose processors, digital signal processors (Digital Signal Processors, DSPs), Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc.
  • the general-purpose processor may be a microprocessor or the processor may also be any conventional processor.
  • a computer-readable storage medium may be a non-volatile computer-readable storage medium.
  • the computer-readable storage medium stores a computer program, where the computer program is executed by a processor to implement the multi-parking vehicle path planning method disclosed in the embodiments of the present application.
  • the storage medium is a physical, non-transitory storage medium, such as a U disk, a mobile hard disk, a read-only memory (Read-Only Memory, ROM), a magnetic disk, or an optical disk, etc., which can store program codes. medium.
  • a physical, non-transitory storage medium such as a U disk, a mobile hard disk, a read-only memory (Read-Only Memory, ROM), a magnetic disk, or an optical disk, etc., which can store program codes. medium.

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Game Theory and Decision Science (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

A multi-depot vehicle routing method, comprising: determining whether a routing request sent by a client terminal is received (S110); if a routing request sent by a client terminal is received, acquiring input data and an end condition corresponding to the routing request, the input data corresponding to the routing request comprising a current number of users, a cargo capacity for each user corresponding to the current number of users, and current user position information of each user corresponding to the current number of users (S120); calling a pre-stored vehicle routing multi-objective optimization model, using the input data as input for the vehicle routing multi-objective optimization model, and performing many-objective evolutionary solving according to the end condition and the vehicle routing multi-objective optimization model, obtaining a route optimization solution set (S130); sending the route optimization solution set to the client terminal (S140). Also provided are a multi-depot vehicle routing apparatus, a computer device and a storage medium.

Description

多车场车辆路径规划方法、装置、计算机设备及存储介质Multi-depot vehicle path planning method, device, computer equipment and storage medium
本申请要求于2020年1月15日提交中国专利局、申请号为202010042242.2、申请名称为“多车场车辆路径规划方法、装置、计算机设备及存储介质”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。This application claims the priority of a Chinese patent application filed with the Chinese Patent Office on January 15, 2020, the application number is 202010042242.2, and the application name is "Multi-yard vehicle path planning method, device, computer equipment and storage medium", and its entire content Incorporated in this application by reference.
技术领域Technical field
本申请涉及路径规划技术领域,尤其涉及一种多车场车辆路径规划方法、装置、计算机设备及存储介质。This application relates to the technical field of route planning, and in particular to a method, device, computer equipment, and storage medium for multi-yard vehicle route planning.
背景技术Background technique
随着电子商务的蓬勃发展,物流产业的重要性越来越凸显。现代物流综合了信息、运输、仓储、库存等多种活动,运用计算机技术科学地进行物流调度仍是物流及运输产业发展的关键。With the vigorous development of e-commerce, the importance of the logistics industry has become more and more prominent. Modern logistics integrates various activities such as information, transportation, warehousing, inventory, and the use of computer technology to scientifically carry out logistics scheduling is still the key to the development of logistics and transportation industries.
对于各大物流企业,其所面临的最重要的问题之一就是如何制定科学的运输路线以高效的满足客户的配送需求,因此车辆路径的规划对整个物流系统的运输成本和效率都有极其重要的影响。随着物流运输规模的日益加大以及配送要求的不断提高,已经无法通过现有的路径规划方法在多配送中心分布的情况下进行最优路径规划。For major logistics companies, one of the most important problems they face is how to formulate scientific transportation routes to efficiently meet the distribution needs of customers. Therefore, the planning of vehicle paths is extremely important to the transportation cost and efficiency of the entire logistics system. Impact. With the increasing scale of logistics transportation and the continuous improvement of distribution requirements, it is no longer possible to use the existing path planning methods to carry out optimal path planning in the case of multiple distribution centers.
发明内容Summary of the invention
本申请实施例提供了一种多车场车辆路径规划方法、装置、计算机设备及存储介质,旨在解决现有技术中在多配送中心分布在不同区域的情况下,无法快速且准确的进行最优路径规划的问题。The embodiments of the present application provide a method, device, computer equipment, and storage medium for vehicle path planning in multiple depots, aiming to solve the problem that the prior art cannot quickly and accurately perform optimization when multiple distribution centers are distributed in different areas. The problem of path planning.
第一方面,本申请实施例提供了一种多车场车辆路径规划方法,其包括:In the first aspect, an embodiment of the present application provides a multi-park vehicle path planning method, which includes:
判断是否接收到客户端发送的路径规划请求;Determine whether the path planning request sent by the client is received;
若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息;If the path planning request sent by the client is received, the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and each corresponding to the current number of users. The user’s cargo capacity and current user location information for each user corresponding to the current number of users;
调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集;以及Call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and perform a lot of operations on the vehicle path planning multi-objective optimization model according to the constraint conditions and The evolutionary solution of the goal, the optimal solution set of the path is obtained; and
将所述路径最优解集发送至客户端。Send the path optimal solution set to the client.
第二方面,本申请实施例提供了一种多车场车辆路径规划装置,其包括用于执行上述第一方面所述的多车场车辆路径规划方法的单元。In a second aspect, an embodiment of the present application provides a multi-park vehicle path planning device, which includes a unit for executing the multi-park vehicle path planning method described in the first aspect.
第三方面,本申请实施例又提供了一种计算机设备,其包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器 执行所述计算机程序时实现上述第一方面所述的多车场车辆路径规划方法。In a third aspect, an embodiment of the present application provides a computer device, which includes a memory, a processor, and a computer program stored on the memory and running on the processor, and the processor executes the computer The program implements the multi-depot vehicle path planning method described in the first aspect.
第四方面,本申请实施例还提供了一种计算机可读存储介质,其中所述计算机可读存储介质存储有计算机程序,所述计算机程序当被处理器执行时使所述处理器执行上述第一方面所述的多车场车辆路径规划方法。In a fourth aspect, an embodiment of the present application also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the processor executes the above-mentioned first On the one hand, the multi-depot vehicle path planning method.
附图说明Description of the drawings
为了更清楚地说明本申请实施例技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the technical solutions of the embodiments of the present application more clearly, the following will briefly introduce the drawings used in the description of the embodiments. Obviously, the drawings in the following description are some embodiments of the present application. Ordinary technicians can obtain other drawings based on these drawings without creative work.
图1为本申请实施例提供的多车场车辆路径规划方法的应用场景示意图;FIG. 1 is a schematic diagram of an application scenario of a multi-park vehicle path planning method provided by an embodiment of the application;
图2为本申请实施例提供的多车场车辆路径规划方法的流程示意图;FIG. 2 is a schematic flowchart of a method for multi-parking vehicle path planning provided by an embodiment of the application;
图3为本申请实施例提供的多车场车辆路径规划方法的子流程示意图;FIG. 3 is a schematic diagram of a sub-process of a multi-park vehicle path planning method provided by an embodiment of the application;
图4为本申请实施例提供的多车场车辆路径规划方法的另一子流程示意图;FIG. 4 is a schematic diagram of another sub-flow of the method for multi-depot vehicle path planning provided by an embodiment of the application;
图5为本申请实施例提供的多车场车辆路径规划装置的示意性框图;Fig. 5 is a schematic block diagram of a multi-park vehicle path planning device provided by an embodiment of the application;
图6为本申请实施例提供的计算机设备的示意性框图。Fig. 6 is a schematic block diagram of a computer device provided by an embodiment of the application.
具体实施方式Detailed ways
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be described clearly and completely in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all of them. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of this application.
应当理解,当在本说明书和所附权利要求书中使用时,术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。It should be understood that when used in this specification and appended claims, the terms "including" and "including" indicate the existence of the described features, wholes, steps, operations, elements and/or components, but do not exclude one or The existence or addition of multiple other features, wholes, steps, operations, elements, components, and/or collections thereof.
还应当理解,在此本申请说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本申请。如在本申请说明书和所附权利要求书中所使用的那样,除非上下文清楚地指明其它情况,否则单数形式的“一”、“一个”及“该”意在包括复数形式。It should also be understood that the terms used in the specification of this application are only for the purpose of describing specific embodiments and are not intended to limit the application. As used in the specification of this application and the appended claims, unless the context clearly indicates other circumstances, the singular forms "a", "an" and "the" are intended to include plural forms.
还应当进一步理解,在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。It should be further understood that the term "and/or" used in the specification and appended claims of this application refers to any combination and all possible combinations of one or more of the associated listed items, and includes these combinations .
请参阅图1和图2,图1为本申请实施例提供的多车场车辆路径规划方法的应用场景示意图;图2为本申请实施例提供的多车场车辆路径规划方法的流程示意图,该多车场车辆路径规划方法应用于服务器中,该方法通过安装于服务器中的应用软件进行执行。Please refer to FIGS. 1 and 2. FIG. 1 is a schematic diagram of an application scenario of a multi-park vehicle path planning method provided by an embodiment of the application; FIG. 2 is a schematic flowchart of a multi-park vehicle path planning method provided by an embodiment of the application. The vehicle path planning method is applied to the server, and the method is executed by application software installed in the server.
如图2所示,该方法包括步骤S110~S140。As shown in Figure 2, the method includes steps S110 to S140.
S110、判断是否接收到客户端发送的路径规划请求。S110: Determine whether a path planning request sent by the client is received.
为了更清楚的理解本申请的技术方案,下面对所涉及到的终端进行介绍。本申请是在服务器的角度描述技术方案。In order to understand the technical solution of the present application more clearly, the terminals involved are introduced below. This application describes the technical solution from the perspective of the server.
第一是客户端,客户端可以理解为用户终端,用户终端可以是智能手机、平板电脑、笔记本电脑、台式电脑、个人数字助理和穿戴式设备等具有通信功能的电子设备,用户终端发送路径规划请求至服务器。The first is the client. The client can be understood as a user terminal. The user terminal can be a smart phone, tablet computer, notebook computer, desktop computer, personal digital assistant, wearable device and other electronic devices with communication functions. The user terminal sends a path plan Request to the server.
第二是服务器,服务器接收客户端发送的路径规划请求,根据与所述路径规划请求对应的输入数据和约束条件,及调用预先存储的车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集。服务器中得到所述路径最优解集后发送至客户端。The second is the server. The server receives the path planning request sent by the client, based on the input data and constraint conditions corresponding to the path planning request, and calls the pre-stored vehicle path planning multi-objective optimization model to perform the evolutionary solution of super multi-objectives. Obtain the path optimal solution set. The optimal solution set of the path is obtained from the server and sent to the client.
在本实施例中,通过服务器检测是否接收到客户端发送的路径规划请求,当服务器接收到客户端发送的路径规划请求时则执行后续的步骤S120,当服务器未接收到客户端发送的路径规划请求时则等待预设的延迟时间后再次执行步骤S110。In this embodiment, the server detects whether the path planning request sent by the client is received, and when the server receives the path planning request sent by the client, the subsequent step S120 is executed. When the server does not receive the path planning request sent by the client Upon request, step S110 is executed again after waiting for the preset delay time.
S120、若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息。S120. If the path planning request sent by the client is received, the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and information corresponding to the current number of users. The current user location information of each user corresponding to the cargo capacity of each user and the number of current users.
在本实施例中,若服务器接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件。由于服务器中已经预先存储了车辆路径规划多目标优化模型,后续根据所述输入数据和约束条件即可进行求解,从而得到路径最优解集。In this embodiment, if the server receives the path planning request sent by the client, it obtains the input data and constraint conditions corresponding to the path planning request. Since the vehicle path planning multi-objective optimization model has been pre-stored in the server, it can be solved subsequently according to the input data and constraint conditions, so as to obtain the path optimal solution set.
S130、调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集。S130. Invoke a pre-stored vehicle path planning multi-objective optimization model, use the input data as input to the vehicle path planning multi-objective optimization model, and perform processing on the vehicle path planning multi-objective optimization model according to the constraint conditions and The evolutionary solution of super-multi-objectives obtains the optimal solution set of the path.
在本实施例中,服务器中存储的车辆路径规划多目标优化模型,是一种多车场多车辆的路径规划多目标优化模型,通过对该车辆路径规划多目标优化模型进行求解,使得优化目标都尽可能达到满足的路径调度规划方案。In this embodiment, the vehicle path planning multi-objective optimization model stored in the server is a multi-depot multi-vehicle path planning multi-objective optimization model. By solving the vehicle path planning multi-objective optimization model, the optimization objectives are all Try to achieve a satisfactory path scheduling plan.
在一实施例中,所述车辆路径规划多目标优化模型包括5个优化目标函数,分别记为:In an embodiment, the vehicle path planning multi-objective optimization model includes five optimization objective functions, which are respectively denoted as:
调度车辆数量优化目标函数minf 1(x)、所有车辆总配送持续时间优化目标函数minf 2(x)、所有车辆总行驶距离优化目标函数minf 3(x)、单个车辆最大容量差优化目标函数minf 4(x)、客户等待时间优化目标函数min f 5(x); Dispatching the number of vehicles to optimize the objective function minf 1 (x), the total delivery duration of all vehicles to optimize the objective function minf 2 (x), the total driving distance of all vehicles to optimize the objective function minf 3 (x), the maximum capacity difference of a single vehicle to optimize the objective function minf 4 (x), customer waiting time optimization objective function min f 5 (x);
Figure PCTCN2020079881-appb-000001
Figure PCTCN2020079881-appb-000001
Figure PCTCN2020079881-appb-000002
Figure PCTCN2020079881-appb-000002
Figure PCTCN2020079881-appb-000003
Figure PCTCN2020079881-appb-000003
Figure PCTCN2020079881-appb-000004
Figure PCTCN2020079881-appb-000004
Figure PCTCN2020079881-appb-000005
Figure PCTCN2020079881-appb-000005
其中,所述车辆路径规划多目标优化模型中预先设置有R个车场、每一车场有K辆车辆,每一车辆的总货物容量为Q且最大服务总时间为T;R、S、K和Q的取值为正整数;与所述路径规划请求对应的输入数据中包括的当前用户数量记为P;当前用户数量P对应的P个客户节点分别记为节点1至节点P,R个车场对应的车场节点分别记为节点P+1至节点P+R;Wherein, the vehicle path planning multi-objective optimization model is preset with R parking lots, each parking lot has K vehicles, the total cargo capacity of each vehicle is Q and the maximum total service time is T; R, S, K and The value of Q is a positive integer; the current number of users included in the input data corresponding to the path planning request is denoted as P; the P customer nodes corresponding to the current number of users P are denoted as node 1 to node P, and R parking lots The corresponding parking lot nodes are respectively marked as node P+1 to node P+R;
K r表示第r个车场的实际使用车辆数;rk表示第r个车场的k编号的车辆;d ij表示第i个节点到第j个节点之间的距离,
Figure PCTCN2020079881-appb-000006
表示第r个车场的k编号的车辆从第i节点运动第j个节点的行驶时间,s i表示第i节点对应的服务时间,s j表示第j节点对应的服务时间,p i表示第i节点对应的第i个客户的货物容量,
Figure PCTCN2020079881-appb-000007
表示第r个车场的k编号的车辆从第i节点运动第j个节点的路径访问状态,
Figure PCTCN2020079881-appb-000008
表示第r个车场的k编号的车辆从第i节点运动第j个节点的客户被服务状态。
K r represents the r-th actual number of yard vehicles; RK k represents a vehicle ID of the r-yard; d ij represents the distance between node i to the j-th node,
Figure PCTCN2020079881-appb-000006
Represents the travel time of the k numbered vehicle of the r-th parking lot from the i-th node to the j-th node, s i represents the service time corresponding to the i-th node, s j represents the service time corresponding to the j-th node, and p i represents the i-th node The cargo capacity of the i-th customer corresponding to the node,
Figure PCTCN2020079881-appb-000007
Represents the path access state of the k numbered vehicle of the rth parking lot moving from the ith node to the jth node,
Figure PCTCN2020079881-appb-000008
It means that the k-numbered vehicle of the r-th parking lot moves from the i-th node to the customer of the j-th node and is served.
其中,当第r个车场的k编号的车辆访问了从第i节点到第j节点的路径,则
Figure PCTCN2020079881-appb-000009
为1,否则
Figure PCTCN2020079881-appb-000010
为0;同理当第r个车场的k编号的车辆服务了第i个节点(此时i为客户),则为
Figure PCTCN2020079881-appb-000011
否则
Figure PCTCN2020079881-appb-000012
为0。
Among them, when the k-numbered vehicle of the r-th parking lot has visited the path from the i-th node to the j-th node, then
Figure PCTCN2020079881-appb-000009
Is 1, otherwise
Figure PCTCN2020079881-appb-000010
Is 0; in the same way, when the k-numbered vehicle in the r-th parking lot serves the i-th node (in this case, i is a customer), then
Figure PCTCN2020079881-appb-000011
otherwise
Figure PCTCN2020079881-appb-000012
Is 0.
即多车场车辆路径规划问题可以定义为:假设拥有R个车场,每个车场有K辆总货物容量为Q且最大服务总时间(路径行驶时间加上客户服务时间)为T的车辆,有P个客户需要配送,第i个客户的货物容量p i<Q,且该客户p i的服务时间s i<T。每个客户可以被任意车辆服务但只能被服务一次,每辆车可服务多个用户且当服务结束后可以被要求返回始发车场。 That is, the multi-park vehicle path planning problem can be defined as: assuming that there are R parking lots, each parking lot has K vehicles with a total cargo capacity of Q and a maximum total service time (path travel time plus customer service time) of T, and there is P A customer needs delivery, the cargo capacity of the i-th customer p i <Q, and the service time s i <T of the customer p i. Each customer can be served by any vehicle but can only be served once, each vehicle can serve multiple users and can be required to return to the departure depot when the service is over.
对于所述车辆路径规划多目标优化模型的一个候选解x,其指满足以上5个优化目标函数(即满足minf 1(x)、minf 2(x)、minf 3(x)、minf 4(x)、minf 5(x))的一条路径,X表示包含多个候选解的集合,多个路径最优解组成的路径最优解集X 最优For a candidate solution x of the vehicle path planning multi-objective optimization model, it refers to satisfying the above 5 optimization objective functions (that is, satisfying minf 1 (x), minf 2 (x), minf 3 (x), minf 4 (x) ), a path of minf 5 (x)), X represents a set containing multiple candidate solutions, and the path optimal solution set X is optimal composed of multiple path optimal solutions.
以所述车辆路径规划多目标优化模型进行路径最优解集的获取时,该模型是一种高维优化模型,结合上述5个目标函数即输入数据和约束条件求解得到路径最优解集X 最优时,能够最大地满足所提出的优化目标和约束条件,更具体即最适当的车辆调度数量,较小的车辆总配送持续时间和车辆总行驶路程,适当的单车最大容量差,以及对客户而言较小的等待时间。 When the vehicle path planning multi-objective optimization model is used to obtain the path optimal solution set, the model is a high-dimensional optimization model that combines the above five objective functions, namely input data and constraint conditions, to obtain the path optimal solution set X When it is optimal , the proposed optimization goals and constraints can be met to the greatest extent. More specifically, the most appropriate number of vehicle dispatches, the smaller the total delivery duration of vehicles and the total travel distance of vehicles, the appropriate maximum capacity difference of single vehicles, and the Smaller waiting time for customers.
在本实施例中,与所述路径规划请求对应的约束条件如下:In this embodiment, the constraint conditions corresponding to the path planning request are as follows:
Figure PCTCN2020079881-appb-000013
Figure PCTCN2020079881-appb-000013
其中,“s.t.”全称为subject to并表示受限制于(一般约束条件的表述以s.t.开始),for each i∈{1,2,...,P}:
Figure PCTCN2020079881-appb-000014
表示对取值属于集合{1,2,……,P}中i值均 使得
Figure PCTCN2020079881-appb-000015
通过所述数据数据和所述约束条件,即可对车辆路径规划多目标优化模型进行求解。
Among them, "st" is called subject to and means to be restricted (the expression of general constraints starts with st), for each i∈{1,2,...,P}:
Figure PCTCN2020079881-appb-000014
Indicates that the values of i in the set {1,2,……,P} are such that
Figure PCTCN2020079881-appb-000015
Through the data and the constraint conditions, the vehicle path planning multi-objective optimization model can be solved.
在一实施例中,如图3所示,所述步骤S130包括:In an embodiment, as shown in FIG. 3, the step S130 includes:
S1301、根据所述约束条件随机生成初始多目标种群;其中,所述初始多目标种群中包括多个个体,每一个体对应所述车辆路径规划多目标优化模型的一个路径输出解,所述初始多目标种群中包括多个个体的总个数记为种群大小N;S1301. Randomly generate an initial multi-objective population according to the constraint conditions; wherein the initial multi-objective population includes multiple individuals, and each individual corresponds to a path output solution of the vehicle path planning multi-objective optimization model, and the initial The total number of multiple individuals in a multi-target population is recorded as the population size N;
S1302、获取当前迭代代数,判断所述当前迭代代数是否达到预设的最大迭代代数;S1302: Obtain the current iteration algebra, and determine whether the current iteration algebra reaches a preset maximum iteration algebra;
S1303、若所述当前迭代代数未达到所述最大迭代代数,获取所述初始多目标种群中的理想个体和最差个体;其中,所述理想个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最小目标值,所述最差个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最大目标值;S1303. If the current iteration algebra does not reach the maximum iteration algebra, obtain the ideal individual and the worst individual in the initial multi-objective population; wherein the ideal individual is input to the vehicle path planning multi-objective optimization model to obtain The target value of is the smallest target value among the target values corresponding to each individual in the initial multi-objective population, and the worst individual is input to the vehicle path planning multi-objective optimization model to obtain the target value for each individual in the initial multi-objective population The largest target value among the corresponding target values;
S1304、对所述初始多目标种群进行模拟二进制交叉和多项式变异,得到与所述初始多目标种群有相同个体总个数的子种群;S1304. Perform simulated binary crossover and polynomial mutation on the initial multi-target population to obtain a subpopulation with the same total number of individuals as the initial multi-target population;
S1305、将所述初始多目标种群与所述子种群进行合并,得到混合种群;S1305. Combine the initial multi-target population and the sub-population to obtain a mixed population;
S1306、将所述混合种群中的个体进行非支配排序,得到非支配解集及多层解集;所述非支配解集记为Q 1,所述多层解集中包括多个解集子集且分别记为Q 2至Q L,其中Q 1至Q L的并集为所述混合种群,Q 1至Q L中任意两个集合的交集为空集,Q 1≥Q 2≥Q 3≥……≥Q LS1306. Perform non-dominated sorting of individuals in the mixed population to obtain a non-dominated solution set and a multi-layer solution set; the non-dominated solution set is denoted as Q 1 , and the multi-layer solution set includes multiple solution set subsets And respectively denoted as Q 2 to Q L , where the union of Q 1 to Q L is the mixed population, the intersection of any two sets from Q 1 to Q L is the empty set, Q 1 ≥Q 2 ≥Q 3 ≥ ……≥Q L ;
S1307、在所述非支配解集、及多层解集中多个解集子集依序合并从而获取多个集合直至个体的总个数超出所述种群大小N,以组成存档集合;S1307. In the non-dominated solution set and the multi-layer solution set, multiple solution set subsets are sequentially merged to obtain multiple sets until the total number of individuals exceeds the population size N to form an archive set;
S1308、根据所述理想个体、所述最差个体将所述存档集合中每一个个体进行归一化处理,得到归一化存档集合;其中,所述归一化存档集合中与所述非支配解集对应的归一化个体集合记为归一化非支配解集;S1308. Normalize each individual in the archive set according to the ideal individual and the worst individual to obtain a normalized archive set; wherein, in the normalized archive set and the non-dominated The normalized individual set corresponding to the solution set is recorded as the normalized non-dominated solution set;
S1309、根据所述归一化非支配解集估计获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面;S1309: Obtain the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the normalized non-dominated solution set estimation;
S1310、将所述归一化存档集合中每一个体映射至所述超曲面上,以得到映射集合;其中,所述归一化存档集合中每一个体均对应所述超曲面上的一个映射点,以组成所述映射集合,所述映射集合中对应的映射点的个数记为L maxS1310. Map each individual in the normalized archive set to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to a mapping on the hypersurface Points to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max ;
S1311、调用预先存储的目标点适应值算法,获取所述映射集合中每一映射点对应的适应值;S1311. Invoke a pre-stored target point fitness algorithm to obtain the fitness value corresponding to each mapping point in the mapping set;
S1312、将所述映射集合中各映射点之间在所述超平面中对应的欧氏距离、及所述种群大小N进行聚类,得到聚类结果;其中,所述聚类结果中的聚类簇的总数与所述种群大小N相等;S1312. Cluster the corresponding Euclidean distance in the hyperplane between the mapping points in the mapping set and the population size N to obtain a clustering result; wherein the clustering in the clustering result The total number of clusters is equal to the population size N;
S1313、在所述聚类结果的每一聚类簇中均挑选一个映射点,以组成目标映射点集合;S1313. Select a mapping point in each cluster cluster of the clustering result to form a target mapping point set;
S1314、获取所述目标映射点集合中各目标映射点对应的个体,以组成当前多目标种群,将所述当前多目标种群作为初始多目标种群;S1314. Obtain an individual corresponding to each target mapping point in the target mapping point set to form a current multi-target population, and use the current multi-target population as an initial multi-target population;
S1315、将所述当前迭代代数加一以作为当前迭代代数,返回执行判断所述当前迭代代数是否达到预设的最大迭代代数的步骤;S1315. Add one to the current iteration algebra as the current iteration algebra, and return to execute the step of judging whether the current iteration algebra reaches the preset maximum iteration algebra;
S1316、若所述当前迭代代数达到所述最大迭代代数,将所述当前多目标种群输出作为路径最优解集。S1316: If the current iteration algebra reaches the maximum iteration algebra, output the current multi-target population as a path optimal solution set.
在本实施例中,在约束条件的限制下随机生成一个初始多目标种群,该初始多目标种群为第一代多目标种群,此时先判断当前迭代代数是否达到预设的最大迭代代数,以确定是否继续迭代执行后续步骤以获取路径最优解集。其中,当前迭代代数的初始值设置为1。若当前迭代代数达到了所述最大迭代代数,将所述当前多目标种群输出作为路径最优解集。In this embodiment, an initial multi-target population is randomly generated under the restriction of constraint conditions. The initial multi-target population is the first-generation multi-target population. At this time, it is first judged whether the current iteration algebra reaches the preset maximum iteration algebra. Determine whether to continue to iteratively execute the next steps to obtain the optimal solution set of the path. Among them, the initial value of the current iteration algebra is set to 1. If the current iteration algebra reaches the maximum iteration algebra, the current multi-target population output is used as the path optimal solution set.
若当前迭代代数未达到所述最大迭代代数,先在所述初始多目标种群寻找初始多目标种群中的理想个体和最差个体,理想个体的每个目标值表示初始多目标种群在对应目标函数的最小目标值,即该理想个体对应的路径输出解代入f 1(x)至f 5(x)后均对应最小的目标值;最差个体的每个目标值表示当前种群在对应目标函数的最大目标值,即该最差个体对应的路径输出解代入f 1(x)至f 5(x)后均对应最大的目标值。 If the current iteration algebra does not reach the maximum iteration algebra, first find the ideal individual and the worst individual in the initial multi-target population in the initial multi-target population. Each target value of the ideal individual indicates that the initial multi-target population is in the corresponding objective function The minimum target value of the ideal individual, that is, the path output solution corresponding to the ideal individual is substituted into f 1 (x) to f 5 (x), and all correspond to the minimum target value; each target value of the worst individual represents the current population in the corresponding target function The maximum target value, that is, the path output solution corresponding to the worst individual corresponds to the maximum target value after substituting f 1 (x) to f 5 (x).
之后需要根据初始多目标种群生成子种群,此过程中可以先对所述初始多目标种群进行模拟二进制交叉和多项式变异,得到与所述初始多目标种群有相同个体总个数的子种群。(子种群的生成是每次随机从当前的初始种群里选择两个个体进行模拟二进制交叉,直到交叉到N个新个体,再根据变异概率和多项式变异对N个新个体进行变异,得到更新的N个个体,这更新的N个个体组成子种群)。Afterwards, it is necessary to generate subpopulations based on the initial multi-target population. In this process, the initial multi-target population can be simulated binary crossover and polynomial mutation to obtain a subpopulation with the same total number of individuals as the initial multi-target population. (The generation of subpopulation is to randomly select two individuals from the current initial population to simulate binary crossover until N new individuals are crossed. Then, the N new individuals are mutated according to the mutation probability and polynomial mutation to get the updated N individuals, the updated N individuals form a subpopulation).
也即在所述初始多目标种群中任意挑选两个个体以依次进行二进制交叉,直到生成N个交叉处理后新个体,对N个交叉处理后新个体进行多项式变异,由多项式变异后的新个体组成子种群。That is, two individuals are randomly selected from the initial multi-target population to perform binary crossover in sequence until N cross-processed new individuals are generated, and the N cross-processed new individuals are subjected to polynomial mutation, and the new individual after polynomial mutation Form subpopulations.
在本实施例中,根据所述初始多目标种群中任意挑选两个个体进行二进制交叉处理后,得到N个交交叉处理后新个体。这里多次任意挑选两个个体进行二进制交叉的过程也类似于一种迭代过程,直到新个体数达到种群大小N,才停止上述多次二进制交叉的处理过程。另外,二进制交叉和多项式变异均为常规处理过程,此处不再赘述。In this embodiment, after two individuals are randomly selected from the initial multi-target population to perform binary crossover processing, N new individuals after crossover processing are obtained. Here, the process of randomly selecting two individuals for binary crossover is similar to an iterative process. Until the number of new individuals reaches the population size N, the process of multiple binary crossovers is stopped. In addition, binary crossover and polynomial mutation are both conventional processing procedures, and will not be repeated here.
之后将所述初始多目标种群与所述子种群进行合并,得到混合种群后,所述混合种群中所包括个体的总个数为所述种群大小N的2倍。Then, the initial multi-target population and the sub-population are combined to obtain a mixed population, and the total number of individuals included in the mixed population is twice the population size N.
此时,可对所述混合种群中各个体进行非支配排序,从而得到非支配解集和多层解集。具体对所述混合种群中各个体进行非支配排序时,可通过非支配解(也可以称为帕累托解)的获取方式,来得到与所述混合种群对应的非支配解集。其中,帕累托解的定义为假设任何二解S1及S2对所有目标而言,S1均优于或同于S2,并且存在至少一个目标,S1在该目标上对应的目标值优于S2该目标上对应的目标值,则称S1支配S2,若S1的解没有被其他解所支配,则S1称为非支配解(不受支配解),也称Pareto解(即帕累托解)。具体的,对所述混合种群中求解非支配解时,得到的非支配解集记为Q 1。所述混合种群中去掉非支配解集对应的个体之后,得到的为多层解集,所述多层解集中包括多 个解集子集且分别记为Q 2至Q L,其中Q 1至Q L的并集为所述混合种群,Q 1至Q L中任意两个集合的交集为空集,Q 1≥Q 2≥Q 3≥……≥Q L;其中,“≥”表示支配关系,Q i≥Q j表示存在Q i中的解支配Q j,该关系是具有传递性,Q 1≥Q 2表示对于f 1(x)至f 5(x)而言,Q2中的每个解都至少被Q1中的一个解所支配,该关系具有传递性,即Q3中的每个解至少被Q1或Q2中的一个解支配,其他的也依次类推。 At this time, non-dominated sorting can be performed on the individuals in the mixed population, thereby obtaining non-dominated solution sets and multi-layer solution sets. When specifically performing non-dominated sorting of the entities in the mixed population, the non-dominated solution set corresponding to the mixed population can be obtained through a non-dominated solution (also called Pareto solution) acquisition method. Among them, the definition of Pareto solution is to assume that for any two solutions S1 and S2, S1 is better than or the same as S2 for all targets, and there is at least one target, and the corresponding target value of S1 on this target is better than S2. The corresponding target value on the target is called S1 dominates S2. If the solution of S1 is not dominated by other solutions, then S1 is called the non-dominated solution (undominated solution), also called the Pareto solution (ie Pareto solution). Specifically, when solving the non-dominated solution in the mixed population, the obtained non-dominated solution set is denoted as Q 1 . After the individuals corresponding to the non-dominated solution set are removed from the mixed population, the multi-layer solution set is obtained. The multi-layer solution set includes multiple solution set subsets and is respectively denoted as Q 2 to Q L , where Q 1 to The union of Q L is the mixed population, and the intersection of any two sets from Q 1 to Q L is an empty set, Q 1 ≥Q 2 ≥Q 3 ≥……≥Q L ; where "≥" indicates a dominance relationship , Q iQ j means that Q j is dominated by the solution in Q i , and the relationship is transitive. Q 1Q 2 means that for f 1 (x) to f 5 (x), each of Q2 The solutions are all dominated by at least one solution in Q1, and the relationship is transitive, that is, each solution in Q3 is dominated by at least one solution in Q1 or Q2, and the others are in turn.
当获取了所述非支配解集、及多层解集后,此时需要挑选超出所述种群大小N的解,以组成存档集合。此时的挑选方式具体如下:先在Q 1中挑选所有的个体,判断当前个体的总个数是否超出所述种群大小N,若当前个体的总个数未超出所述种群大小N,则继续在Q 2中挑选所有的个体,将当前个体的总个数加上Q 2中所述的个体的个数以更新作为当前个体的总个数,再判断当前个体的总个数是否超出所述种群大小N,直至获取的Q 1至Q a中的个体总个数超出所述种群大小N(其中,a为大于1且未超出L的正整数),以组成存档集合。 After the non-dominated solution set and the multi-layer solution set are obtained, solutions that exceed the population size N need to be selected at this time to form an archive set. At this time, the selection mode as follows: first of all individuals in the selected Q 1, determines whether the total number of individuals beyond the current population size N, if the current total number of individuals of the population size does not exceed N, continue Q 2 in the selection of all of the individuals, the total number of individual current plus the number of Q 2 in the subject to update the current as the total number of individuals, and then determines whether the total number of individuals beyond the current The population size N, until the total number of individuals in Q 1 to Q a obtained exceeds the population size N (where a is a positive integer greater than 1 and not exceeding L) to form an archive set.
此时,可以根据所述理想个体、所述最差个体将所述存档集合中每一个个体进行归一化处理,得到归一化存档集合;其中,所述归一化存档集合中与所述非支配解集对应的归一化个体集合记为归一化非支配解集。即所述存档集合中每一个个体均进行归一化,得到了与每个个体对应的归一化个体,从而组成了归一化非支配解集。上述归一化处理的目的在于消除不同量纲之间的差异,从而便于后续的数据处理。At this time, each individual in the archive set may be normalized according to the ideal individual and the worst individual to obtain a normalized archive set; wherein, the normalized archive set and the The normalized individual set corresponding to the non-dominated solution set is recorded as the normalized non-dominated solution set. That is, each individual in the archive set is normalized, and a normalized individual corresponding to each individual is obtained, thereby forming a normalized non-dominated solution set. The purpose of the above normalization processing is to eliminate the difference between different dimensions, so as to facilitate subsequent data processing.
在一实施例中,步骤S1308包括:In an embodiment, step S1308 includes:
根据
Figure PCTCN2020079881-appb-000016
将所述存档集合中每一个个体进行归一化处理,得到与所述存档集合中每一个个体对应的归一化个体,以组成归一化存档集合;其中,NA m表示所述存档集合中第m个个体A m对应的归一化个体,A 最差个体表示所述最差个体,A 理想个体表示所述理想个体。
according to
Figure PCTCN2020079881-appb-000016
Perform normalization processing on each individual in the archive set to obtain a normalized individual corresponding to each individual in the archive set to form a normalized archive set; where NA m represents in the archive set the individual over the individual normalization, a represents the worst of the worst individual subject, a represents the individual over the individual m-th corresponding to a m.
通过上述归一化处理的模型,即可得到归一化存档集合,从而消除不同量纲之间的差异。Through the above-mentioned normalized processing model, a normalized archive set can be obtained, thereby eliminating the difference between different dimensions.
在获取了所述归一化存档集合后,可根据其中包括的归一化个体估计获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面。根据所有非支配解来估计种群在目标函数空间的形状,这将有助于发现多目标优化任务的特点,通常将多目标优化分为三类,凸优化、凹优化以及线性优化,在这一步骤中将得到代表三种类型的超曲面曲率值。After obtaining the normalized archive set, the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front can be estimated according to the normalized individuals included therein. Estimate the shape of the population in the objective function space based on all non-dominated solutions, which will help to discover the characteristics of multi-objective optimization tasks. Multi-objective optimization is usually divided into three categories, convex optimization, concave optimization, and linear optimization. In this step, the curvature values of three types of hypersurfaces will be obtained.
在一实施例中,如图4所示,步骤S1309包括:In one embodiment, as shown in FIG. 4, step S1309 includes:
S13091、将所述归一化非支配解集的各归一化非支配个体中未处于标准化目标空间的归一化非支配个体移除,得到筛选后归一化非支配解集;其中,位于所述标准化目标空间中各归一化非支配个体对应的各目标值均未超过1;S13091. Remove the normalized non-dominated individuals that are not in the standardized target space among the normalized non-dominated individuals in the normalized non-dominated solution set to obtain the normalized non-dominated solution set after screening; where Each target value corresponding to each normalized non-dominated individual in the standardized target space does not exceed 1;
S13092、获取所述筛选后归一化非支配解集中每一非支配个体,分别记为B 1至B n;其中,n的取值与所述筛选后归一化非支配解集中归一化非支配个体的总个数相同; S13092. Obtain each non-dominated individual in the normalized non-dominated solution set after the screening, which is respectively recorded as B 1 to B n ; where the value of n is normalized to the normalized non-dominated solution set after the screening The total number of non-dominant individuals is the same;
S13093、获取非支配个体B i到目标超平面对应的超平面距离D i;其中,i的取值范围为[1,n],目标超平面为f 1(x)+f 2(x)+f 3(x)+f 4(x)+f 5(x)=1; S13093. Obtain the hyperplane distance D i corresponding to the target hyperplane from the non-dominated individual B i ; where the value range of i is [1,n], and the target hyperplane is f 1 (x)+f 2 (x)+ f 3 (x)+f 4 (x)+f 5 (x)=1;
S13094、获取超平面距离D 1至D n对应的超平面距离平均值和超平面距离标准差;其中,超平面距离D 1至D n对应的超平面距离平均值记为D avg,超平面距离D 1至D n对应的超平面距离标准差记为D stdS13094, and obtaining a distance mean value hyperplane standard differential hyperplane distance from the hyperplane D 1 to D n corresponding to; wherein hyperplane distance from the hyperplane D 1 to D n corresponding to an average value referred to as D avg, from the hyperplane The standard deviation of the hyperplane distance corresponding to D 1 to D n is denoted as D std ;
S13095、根据所述超平面距离平均值与所述超平面距离标准差之商进行范数运算,得到对应的变异系数;其中,所述变异系数记为cv;S13095. Perform a norm operation according to the quotient of the hyperplane distance average value and the hyperplane distance standard deviation to obtain a corresponding coefficient of variation; wherein the coefficient of variation is denoted as cv;
S13096、获取所述目标超平面的曲率为2对应的第一超曲面,和所述目标超平面的曲率为0.5对应的第二超曲面,根据目标超平面、第一超曲面、第二超曲面及预设的曲率确定策略,获取所述归一化非支配解集对应的当前曲率;其中,所述曲率确定策略为
Figure PCTCN2020079881-appb-000017
d(2.0)表示第一超曲面的峰值点到所述目标超平面的距离,d(0.5)表示第二超曲面的峰值点到所述目标超平面的距离;
S13096. Obtain a first hypersurface corresponding to a curvature of 2 of the target hyperplane, and a second hypersurface corresponding to a curvature of 0.5 of the target hyperplane, according to the target hyperplane, the first hypersurface, and the second hypersurface And a preset curvature determination strategy to obtain the current curvature corresponding to the normalized non-dominated solution set; wherein, the curvature determination strategy is
Figure PCTCN2020079881-appb-000017
d(2.0) represents the distance from the peak point of the first hypersurface to the target hyperplane, and d(0.5) represents the distance from the peak point of the second hypersurface to the target hyperplane;
S13097、根据所述变异系数对所述当前曲率进行调整,得到调整后曲率;其中,若所述变异系数小于0.1,将所述当前曲率的取值调整为1,以使调整后曲率取值为1;若所述变异系数大于或等于0.1,将所述当前曲率的取值保持不变,以使得调整后曲率等于所述当前曲率;S13097. Adjust the current curvature according to the coefficient of variation to obtain an adjusted curvature; wherein, if the coefficient of variation is less than 0.1, adjust the value of the current curvature to 1, so that the adjusted curvature value is 1; if the coefficient of variation is greater than or equal to 0.1, keep the value of the current curvature unchanged, so that the adjusted curvature is equal to the current curvature;
S13098、根据所述调整后曲率对应获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面。S13098: Obtain the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the adjusted curvature.
在本实施例中,由于归一化非支配解集中仍可能存在未处于标准化目标空间的归一化非支配个体,此时将所述归一化非支配解集的各归一化非支配个体中未处于标准化目标空间(标准化目标空间是f 1(x)至f 5(x)均等于1而形成的目标空间)的归一化非支配个体移除,得到筛选后归一化非支配解集。移除未处于标准化目标空间的归一化非支配个体后,有效的排除了这些个体的干扰作用,有利于后续的种群进化。 In this embodiment, since there may still be normalized non-dominated individuals in the normalized non-dominated solution set that are not in the standardized target space, at this time, each normalized non-dominated individual in the normalized non-dominated solution set Remove the normalized non-dominated individuals that are not in the standardized target space (the standardized target space is the target space formed by f 1 (x) to f 5 (x) equal to 1), and get the normalized non-dominated solution after screening set. After removing the normalized non-dominated individuals that are not in the standardized target space, the interference of these individuals is effectively eliminated, which is beneficial to the subsequent population evolution.
其中,
Figure PCTCN2020079881-appb-000018
m在本申请中的取值为5,当p=2时即可获取d(2.0)的值,当p=0.5时即可获取d(0.5)的值。Cur表示所述归一化非支配解集对应的当前曲率。
in,
Figure PCTCN2020079881-appb-000018
The value of m in this application is 5. When p=2, the value of d(2.0) can be obtained, and when p=0.5, the value of d(0.5) can be obtained. Cur represents the current curvature corresponding to the normalized non-dominated solution set.
此时通过步骤S13091-S13098的处理过程获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面,是为了便于后续的聚类操作,预测过程中使用的个体为所有的非支配个体。At this time, the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front are obtained through the processing process of steps S13091-S13098, in order to facilitate the subsequent clustering operation, and the individuals used in the prediction process are all non-dominated individual.
而且,计算获取非支配个体B i到目标超平面的距离D i时,距离D i为负值表示该个体在目标超平面以下,正值表示该个体在目标超平面以上,距离D i为0表示该个体在目标超平面上。通过获取超平面距离D 1至D n对应的超平面距离平均值和超平面距离标准差,即可对应的获取后续调整曲率所需的变异系数cv。 Moreover, when calculating the distance D i from the non-dominated individual B i to the target hyperplane, a negative value of the distance D i means that the individual is below the target hyperplane, a positive value means that the individual is above the target hyperplane, and the distance D i is 0 Indicates that the individual is on the target hyperplane. By obtaining the hyperplane distance average value and the hyperplane distance standard deviation corresponding to the hyperplane distance D 1 to D n , the coefficient of variation cv required for subsequent adjustment of the curvature can be obtained correspondingly.
根据变异系数cv对所述当前曲率进行调整,得到调整后曲率,即可根据所述调整后曲率对应获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面。以帕累托前沿的形状对应的超曲面,将有助于发现多目标优化任务的特点。The current curvature is adjusted according to the coefficient of variation cv to obtain the adjusted curvature. The shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front can be correspondingly obtained according to the adjusted curvature. The hypersurface corresponding to the shape of the Pareto front will help to discover the characteristics of the multi-objective optimization task.
之后,将所述归一化存档集合中每一个体映射至所述超曲面上,以得到映射集合;其中,所述归一化存档集合中每一个体均对应所述超曲面上的一个映射点,以组成所述映射集合,所述映射集合中对应的映射点的个数记为L max。此时,需要获取每一个映射点的适应值,具体是根据
Figure PCTCN2020079881-appb-000019
获取每一个映射点的适应值,f i(x l)表示所述映射集合中第l个映射点对应的第i个目标值,l的取值范围为[1,L max]。每个映射点的适应值表示个体对应的映射点的所有目标值的总和,由于倾向于寻找到收敛良好的解集,因此将适应值作为衡量收敛的指标。
Afterwards, each individual in the normalized archive set is mapped to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to a mapping on the hypersurface Points to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max . At this time, the fitness value of each mapping point needs to be obtained, specifically
Figure PCTCN2020079881-appb-000019
Obtain the fitness value of each mapping point, f i (x l ) represents the i-th target value corresponding to the l-th mapping point in the mapping set, and the value range of l is [1, L max ]. The fitness value of each mapping point represents the sum of all the target values of the individual corresponding mapping points. Since it tends to find a well-converged solution set, the fitness value is used as a measure of convergence.
此时为了将所述映射集合中各映射点进行聚类,且使得所述聚类结果中的聚类簇的总数与所述种群大小N相等,可以根据各映射点之间在所述超平面中对应的欧氏距离进行聚类。若两个映射点之间在所述超平面中对应的欧氏距离越小,表示两个映射点对应的个体越相似。At this time, in order to cluster each mapping point in the mapping set, and make the total number of cluster clusters in the clustering result equal to the population size N, it is possible to set the distance between each mapping point on the hyperplane The corresponding Euclidean distance in clustering. If the corresponding Euclidean distance between two mapping points in the hyperplane is smaller, it means that the individuals corresponding to the two mapping points are more similar.
完成聚类后,从聚类结果中的每一聚类簇中均挑选一个映射点,以组成目标映射点集合,这样目标映射点集合中映射点的个数与所述种群大小N相等,映射点集合对应的个体集合即为本次进化的种群,同时也是下一代的父种群。通过这一选解策略,能保证下一代种群有良好的收敛性和多样性。通过多次迭代所述当前迭代代数达到所述最大迭代代数,将所述当前多目标种群输出作为路径最优解集。其中,所述路径最优解集中的每一个体即为根据所述输入数据和约束条件得到的与车辆路径规划多目标优化模型对应的最优解。After the clustering is completed, a mapping point is selected from each cluster cluster in the clustering result to form a target mapping point set, so that the number of mapping points in the target mapping point set is equal to the population size N, and the mapping The individual set corresponding to the point set is the population of this evolution, and it is also the parent population of the next generation. Through this selection strategy, it can ensure that the next generation population has good convergence and diversity. The current iteration algebra is iterated multiple times to reach the maximum iteration algebra, and the current multi-target population output is used as a path optimal solution set. Wherein, each individual in the path optimal solution set is the optimal solution corresponding to the vehicle path planning multi-objective optimization model obtained according to the input data and constraint conditions.
具体实施时,所述路径最优解集中每一路径最优解的具体编码方式(也即最终展示给用户查看的方式),是用户自定义编码格式并保存在服务器中,此次并不限定具体编码方式。In specific implementation, the specific encoding method of each path optimal solution in the path optimal solution set (that is, the method that is finally displayed to the user) is a user-defined encoding format and saved in the server, which is not limited this time Specific encoding method.
S140、将所述路径最优解集发送至客户端。S140. Send the path optimal solution set to the client.
在本实施例中,当在服务器中完成了路径最优解集的获取之后,即可发送至客户端。从而客户端可根据所述路径最优解集确定配送路径后,以辅助配送过程。In this embodiment, after obtaining the path optimal solution set in the server, it can be sent to the client. Therefore, the client can assist the delivery process after determining the delivery route according to the optimal solution set of the route.
该方法实现了在超多目标的进化求解的过程中充分考虑了种群的收敛和多样性,实现了种群形状的有效预测,实现了基于输入数据和约束条件快速且准确的获取车辆路径规划多目标优化模型的路径最优解集。This method realizes that the convergence and diversity of the population are fully considered in the process of solving super-multi-objective evolution, and the effective prediction of the shape of the population is realized, and the multi-objective of vehicle path planning is obtained quickly and accurately based on input data and constraints. The path optimal solution set of the optimization model.
本申请实施例还提供一种多车场车辆路径规划装置,该多车场车辆路径规划装置用于执行前述多车场车辆路径规划方法的任一实施例。具体地,请参阅图5,图5是本申请实施例提供的多车场车辆路径规划装置的示意性框图。该多车场车辆路径规划装置100可以被配置于服务器中。The embodiment of the present application also provides a multi-park vehicle path planning device, which is used to execute any embodiment of the foregoing multi-park vehicle path planning method. Specifically, please refer to FIG. 5, which is a schematic block diagram of a multi-park vehicle path planning apparatus provided by an embodiment of the present application. The multi-parking vehicle path planning device 100 may be configured in a server.
如图5所示,多车场车辆路径规划装置100包括路径规划请求检测单元110、数据条件获取单元120、路径最优解集获取单元130、及最优解集发送单元140。As shown in FIG. 5, the multi-depot vehicle path planning device 100 includes a path planning request detection unit 110, a data condition acquisition unit 120, a path optimal solution set acquisition unit 130, and an optimal solution set transmission unit 140.
其中,路径规划请求检测单元110,用于判断是否接收到客户端发送的路径 规划请求。Wherein, the path planning request detection unit 110 is used to determine whether the path planning request sent by the client is received.
数据条件获取单元120,用于若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息。The data condition obtaining unit 120 is configured to, if a path planning request sent by the client is received, obtain input data and constraint conditions corresponding to the path planning request; wherein the input data corresponding to the path planning request includes the current number of users , The cargo capacity of each user corresponding to the current number of users, and the current user location information of each user corresponding to the current number of users.
路径最优解集获取单元130,用于调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集。The path optimal solution set acquisition unit 130 is configured to call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and according to the constraint conditions and the control The vehicle path planning multi-objective optimization model performs the evolutionary solution of super multi-objectives, and obtains the optimal solution set of the path.
最优解集发送单元140,用于将所述路径最优解集发送至客户端。The optimal solution set sending unit 140 is configured to send the path optimal solution set to the client.
在一实施例中,所述路径最优解集获取单元130包括:In an embodiment, the path optimal solution set obtaining unit 130 includes:
初始多目标种群生成单元,用于根据所述约束条件随机生成初始多目标种群;其中,所述初始多目标种群中包括多个个体,每一个体对应所述车辆路径规划多目标优化模型的一个路径输出解,所述初始多目标种群中包括多个个体的总个数记为种群大小N;The initial multi-objective population generating unit is used to randomly generate an initial multi-objective population according to the constraint conditions; wherein the initial multi-objective population includes multiple individuals, and each individual corresponds to one of the vehicle path planning multi-objective optimization models Path output solution, the total number of multiple individuals included in the initial multi-target population is recorded as the population size N;
当前迭代代数第一判断单元,用于获取当前迭代代数,判断所述当前迭代代数是否达到预设的最大迭代代数;The first judging unit of the current iteration algebra is used to obtain the current iteration algebra, and judge whether the current iteration algebra reaches the preset maximum iteration algebra;
目标个体获取单元,用于若所述当前迭代代数未达到所述最大迭代代数,获取所述初始多目标种群中的理想个体和最差个体;其中,所述理想个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最小目标值,所述最差个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最大目标值;A target individual obtaining unit, configured to obtain an ideal individual and a worst individual in the initial multi-target population if the current iteration algebra does not reach the maximum iteration algebra; wherein the ideal individual is input to the vehicle path planning The target value obtained by the multi-objective optimization model is the smallest target value among the target values corresponding to each individual in the initial multi-objective population, and the worst individual is input to the vehicle path planning multi-objective optimization model and the target value obtained is the initial multi-objective The largest target value among the target values corresponding to each individual in the population;
个体交叉变异单元,用于对所述初始多目标种群进行模拟二进制交叉和多项式变异,得到与所述初始多目标种群有相同个体总个数的子种群;The individual crossover mutation unit is used to simulate binary crossover and polynomial mutation on the initial multi-target population to obtain a subpopulation with the same total number of individuals as the initial multi-target population;
混合种群获取单元,用于将所述初始多目标种群与所述子种群进行合并,得到混合种群;A mixed population obtaining unit, configured to merge the initial multi-target population and the sub-population to obtain a mixed population;
非支配解集获取单元,用于将所述混合种群中的个体进行非支配排序,得到非支配解集及多层解集;其中,所述非支配解集记为Q 1,所述多层解集中包括多个解集子集且分别记为Q 2至Q L,其中Q 1至Q L的并集为所述混合种群,Q 1至Q L中任意两个集合的交集为空集,Q 1≥Q 2≥Q 3≥……≥Q LThe non-dominated solution set obtaining unit is used to perform non-dominated sorting of the individuals in the mixed population to obtain a non-dominated solution set and a multi-layer solution set; wherein the non-dominated solution set is denoted as Q 1 , and the multi-layer The solution set includes multiple solution set subsets and are respectively denoted as Q 2 to Q L , where the union of Q 1 to Q L is the mixed population, and the intersection of any two sets from Q 1 to Q L is an empty set, Q 1 ≥Q 2 ≥Q 3 ≥……≥Q L ;
存档集合获取单元,用于在所述非支配解集、及多层解集中多个解集子集依序合并从而获取多个集合直至个体的总个数超出所述种群大小N,以组成存档集合;The archive set acquisition unit is used to sequentially merge multiple solution set subsets in the non-dominated solution set and the multi-layer solution set to obtain multiple sets until the total number of individuals exceeds the population size N to form an archive gather;
归一化处理单元,用于根据所述理想个体、所述最差个体将所述存档集合中每一个个体进行归一化处理,得到归一化存档集合;其中,所述归一化存档集合中与所述非支配解集对应的归一化个体集合记为归一化非支配解集;The normalization processing unit is configured to normalize each individual in the archive set according to the ideal individual and the worst individual to obtain a normalized archive set; wherein, the normalized archive set The set of normalized individuals corresponding to the non-dominated solution set is recorded as the normalized non-dominated solution set;
超曲面获取单元,用于根据所述归一化非支配解集估计获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面;A hypersurface acquiring unit, configured to estimate and acquire the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the normalized non-dominated solution set;
个体映射单元,用于将所述归一化存档集合中每一个体映射至所述超曲面 上,以得到映射集合;其中,所述归一化存档集合中每一个体均对应所述超曲面上的一个映射点,以组成所述映射集合,所述映射集合中对应的映射点的个数记为L maxThe individual mapping unit is used to map each individual in the normalized archive set to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to the hypersurface A mapping point above to form the mapping set, and the number of corresponding mapping points in the mapping set is denoted as L max ;
适应值获取单元,用于调用预先存储的目标点适应值算法,获取所述映射集合中每一映射点对应的适应值;The fitness value obtaining unit is used to call a pre-stored target point fitness value algorithm to obtain the fitness value corresponding to each mapping point in the mapping set;
映射点聚类单元,用于将所述映射集合中各映射点之间在所述超平面中对应的欧氏距离、及所述种群大小N进行聚类,得到聚类结果;其中,所述聚类结果中的聚类簇的总数与所述种群大小N相等;The mapping point clustering unit is used to cluster the corresponding Euclidean distance in the hyperplane between each mapping point in the mapping set and the population size N to obtain a clustering result; wherein, the The total number of clusters in the clustering result is equal to the population size N;
目标映射点集合获取单元,用于在所述聚类结果的每一聚类簇中均挑选一个映射点,以组成目标映射点集合;The target mapping point set acquiring unit is configured to select a mapping point in each cluster cluster of the clustering result to form a target mapping point set;
初始多目标种群更新单元,用于获取所述目标映射点集合中各目标映射点对应的个体,以组成当前多目标种群,将所述当前多目标种群作为初始多目标种群;The initial multi-target population update unit is used to obtain the individual corresponding to each target mapping point in the target mapping point set to form a current multi-target population, and use the current multi-target population as the initial multi-target population;
当前迭代代数第二判断单元,用于将所述当前迭代代数加一以作为当前迭代代数,返回执行判断所述当前迭代代数是否达到预设的最大迭代代数的步骤;The current iteration algebra second judging unit, configured to add one to the current iteration algebra as the current iteration algebra, and return to execute the step of judging whether the current iteration algebra reaches the preset maximum iteration algebra;
路径最优解集输出单元,用于若所述当前迭代代数达到所述最大迭代代数,将所述当前多目标种群输出作为路径最优解集。The path optimal solution set output unit is configured to output the current multi-target population as the path optimal solution set if the current iteration algebra reaches the maximum iteration algebra.
在一实施例中,所述个体交叉变异单元还用于:In an embodiment, the individual cross mutation unit is also used for:
在所述初始多目标种群中任意挑选两个个体以依次进行二进制交叉,直到生成N个交叉处理后新个体,对N个交叉处理后新个体进行多项式变异,由多项式变异后的新个体组成子种群。Randomly select two individuals in the initial multi-target population to perform binary crossover in sequence, until N cross-processed new individuals are generated, and the N cross-processed new individuals are subjected to polynomial mutation, and the new individuals after polynomial mutation are composed of subgroups. Population.
在一实施例中,所述归一化处理单元还用于:In an embodiment, the normalization processing unit is further used for:
根据
Figure PCTCN2020079881-appb-000020
将所述存档集合中每一个个体进行归一化处理,得到与所述存档集合中每一个个体对应的归一化个体,以组成归一化存档集合;其中,NA m表示所述存档集合中第m个个体A m对应的归一化个体,A 最差个体表示所述最差个体,A 理想个体表示所述理想个体。
according to
Figure PCTCN2020079881-appb-000020
Perform normalization processing on each individual in the archive set to obtain a normalized individual corresponding to each individual in the archive set to form a normalized archive set; where NA m represents in the archive set the individual over the individual normalization, a represents the worst of the worst individual subject, a represents the individual over the individual m-th corresponding to a m.
在一实施例中,所述超曲面获取单元包括:In an embodiment, the hypersurface acquiring unit includes:
非支配个体筛选单元,用于将所述归一化非支配解集的各归一化非支配个体中未处于标准化目标空间的归一化非支配个体移除,得到筛选后归一化非支配解集;其中,位于所述标准化目标空间中各归一化非支配个体对应的各目标值均未超过1;The non-dominated individuals screening unit is used to remove the normalized non-dominated individuals that are not in the standardized target space among the normalized non-dominated individuals in the normalized non-dominated solution set to obtain the normalized non-dominated individuals after screening Solution set; wherein each target value corresponding to each normalized non-dominated individual located in the standardized target space does not exceed 1;
非支配个体获取单元,用于获取所述筛选后归一化非支配解集中每一非支配个体,分别记为B 1至B n;其中,n的取值与所述筛选后归一化非支配解集中归一化非支配个体的总个数相同; The non-dominated individual acquisition unit is used to acquire each non-dominated individual in the normalized non-dominated solution set after the screening, which is respectively recorded as B 1 to B n ; where the value of n is the same as the normalized non-dominated after screening. The total number of normalized non-dominated individuals in the dominating solution set is the same;
超平面距离获取单元,用于获取非支配个体B i到目标超平面对应的超平面距离D i;其中,i的取值范围为[1,n],目标超平面为f 1(x)+f 2(x)+f 3(x)+f 4(x)+f 5(x)=1; Hyperplane distance acquisition unit for acquiring the individual B i to a non-dominant hyperplane hyperplane corresponding target distance D i; where i is the range [1, n], the target hyperplane f 1 (x) + f 2 (x)+f 3 (x)+f 4 (x)+f 5 (x)=1;
超平面距离参数获取单元,用于获取超平面距离D 1至D n对应的超平面距离 平均值和超平面距离标准差;其中,超平面距离D 1至D n对应的超平面距离平均值记为D avg,超平面距离D 1至D n对应的超平面距离标准差记为D stdHyperplane from the parameter acquisition unit, configured to obtain hyperplane distance D 1 to D n corresponding hyperplane distance from average and standard deviation hyperplane; wherein hyperplane distance D 1 to D n corresponding to the average distance referred hyperplane Is D avg , the standard deviation of the hyperplane distance corresponding to the hyperplane distance D 1 to D n is denoted as D std ;
变异系数获取单元,用于根据所述超平面距离平均值与所述超平面距离标准差之商进行范数运算,得到对应的变异系数;其中,所述变异系数记为cv;The coefficient of variation obtaining unit is configured to perform a norm operation according to the quotient of the average value of the hyperplane distance and the standard deviation of the hyperplane distance to obtain the corresponding coefficient of variation; wherein the coefficient of variation is denoted as cv;
当前曲率获取单元,用于获取所述目标超平面的曲率为2对应的第一超曲面,和所述目标超平面的曲率为0.5对应的第二超曲面,根据目标超平面、第一超曲面、第二超曲面及预设的曲率确定策略,获取所述归一化非支配解集对应的当前曲率;其中,所述曲率确定策略为
Figure PCTCN2020079881-appb-000021
d(2.0)表示第一超曲面的峰值点到所述目标超平面的距离,d(0.5)表示第二超曲面的峰值点到所述目标超平面的距离;
The current curvature acquisition unit is configured to acquire a first hypersurface corresponding to a curvature of 2 of the target hyperplane, and a second hypersurface corresponding to a curvature of 0.5 of the target hyperplane, according to the target hyperplane and the first hypersurface , The second hypersurface and the preset curvature determination strategy, to obtain the current curvature corresponding to the normalized non-dominated solution set; wherein, the curvature determination strategy is
Figure PCTCN2020079881-appb-000021
d(2.0) represents the distance from the peak point of the first hypersurface to the target hyperplane, and d(0.5) represents the distance from the peak point of the second hypersurface to the target hyperplane;
曲率调整单元,用于根据所述变异系数对所述当前曲率进行调整,得到调整后曲率;其中,若所述变异系数小于0.1,将所述当前曲率的取值调整为1,以使调整后曲率取值为1;若所述变异系数大于或等于0.1,将所述当前曲率的取值保持不变,以使得调整后曲率等于所述当前曲率;The curvature adjustment unit is configured to adjust the current curvature according to the coefficient of variation to obtain an adjusted curvature; wherein, if the coefficient of variation is less than 0.1, adjust the value of the current curvature to 1, so that the adjusted curvature The value of the curvature is 1; if the coefficient of variation is greater than or equal to 0.1, the value of the current curvature is kept unchanged, so that the adjusted curvature is equal to the current curvature;
超曲面生成单元,用于根据所述调整后曲率对应获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面。The hypersurface generating unit is used to obtain the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the adjusted curvature.
该装置实现了在超多目标的进化求解的过程中充分考虑了种群的收敛和多样性,实现了种群形状的有效预测,实现了基于输入数据和约束条件快速且准确的获取车辆路径规划多目标优化模型的路径最优解集。The device realizes that the convergence and diversity of the population are fully considered in the process of solving super-multi-objective evolution, and the effective prediction of the shape of the population is realized, and the multi-objective of vehicle path planning is obtained quickly and accurately based on input data and constraint conditions. The path optimal solution set of the optimization model.
上述多车场车辆路径规划装置可以实现为计算机程序的形式,该计算机程序可以在如图6所示的计算机设备上运行。The above-mentioned multi-yard vehicle path planning apparatus may be implemented in the form of a computer program, and the computer program may run on a computer device as shown in FIG. 6.
请参阅图6,图6是本申请实施例提供的计算机设备的示意性框图。该计算机设备500是服务器,服务器可以是独立的服务器,也可以是多个服务器组成的服务器集群。Please refer to FIG. 6, which is a schematic block diagram of a computer device according to an embodiment of the present application. The computer device 500 is a server, and the server may be an independent server or a server cluster composed of multiple servers.
参阅图6,该计算机设备500包括通过系统总线501连接的处理器502、存储器和网络接口505,其中,存储器可以包括非易失性存储介质503和内存储器504。Referring to FIG. 6, the computer device 500 includes a processor 502, a memory, and a network interface 505 connected through a system bus 501, where the memory may include a non-volatile storage medium 503 and an internal memory 504.
该非易失性存储介质503可存储操作系统5031和计算机程序5032。该计算机程序5032被执行时,可使得处理器502执行多车场车辆路径规划方法。The non-volatile storage medium 503 can store an operating system 5031 and a computer program 5032. When the computer program 5032 is executed, the processor 502 can execute the multi-park vehicle path planning method.
该处理器502用于提供计算和控制能力,支撑整个计算机设备500的运行。The processor 502 is used to provide calculation and control capabilities, and support the operation of the entire computer device 500.
该内存储器504为非易失性存储介质503中的计算机程序5032的运行提供环境,该计算机程序5032被处理器502执行时,可使得处理器502执行多车场车辆路径规划方法。The internal memory 504 provides an environment for the operation of the computer program 5032 in the non-volatile storage medium 503. When the computer program 5032 is executed by the processor 502, the processor 502 can make the processor 502 execute the multi-park vehicle path planning method.
该网络接口505用于进行网络通信,如提供数据信息的传输等。本领域技术人员可以理解,图6中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备500的限定,具体的 计算机设备500可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。The network interface 505 is used for network communication, such as providing data information transmission. Those skilled in the art can understand that the structure shown in FIG. 6 is only a block diagram of part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device 500 to which the solution of the present application is applied. The specific computer device 500 may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.
其中,所述处理器502用于运行存储在存储器中的计算机程序5032,以实现本申请实施例公开的多车场车辆路径规划方法。Wherein, the processor 502 is configured to run a computer program 5032 stored in a memory to implement the multi-park vehicle path planning method disclosed in the embodiment of the present application.
本领域技术人员可以理解,图6中示出的计算机设备的实施例并不构成对计算机设备具体构成的限定,在其他实施例中,计算机设备可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。例如,在一些实施例中,计算机设备可以仅包括存储器及处理器,在这样的实施例中,存储器及处理器的结构及功能与图6所示实施例一致,在此不再赘述。Those skilled in the art can understand that the embodiment of the computer device shown in FIG. 6 does not constitute a limitation on the specific configuration of the computer device. In other embodiments, the computer device may include more or less components than those shown in the figure. Or combine certain components, or different component arrangements. For example, in some embodiments, the computer device may only include a memory and a processor. In such an embodiment, the structure and function of the memory and the processor are consistent with the embodiment shown in FIG. 6 and will not be repeated here.
应当理解,在本申请实施例中,处理器502可以是中央处理单元(Central Processing Unit,CPU),该处理器502还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。其中,通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。It should be understood that, in this embodiment of the application, the processor 502 may be a central processing unit (Central Processing Unit, CPU), and the processor 502 may also be other general-purpose processors, digital signal processors (Digital Signal Processors, DSPs), Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor.
在本申请的另一实施例中提供计算机可读存储介质。该计算机可读存储介质可以为非易失性的计算机可读存储介质。该计算机可读存储介质存储有计算机程序,其中计算机程序被处理器执行时实现本申请实施例公开的多车场车辆路径规划方法。In another embodiment of the present application, a computer-readable storage medium is provided. The computer-readable storage medium may be a non-volatile computer-readable storage medium. The computer-readable storage medium stores a computer program, where the computer program is executed by a processor to implement the multi-parking vehicle path planning method disclosed in the embodiments of the present application.
所述存储介质为实体的、非瞬时性的存储介质,例如可以是U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、磁碟或者光盘等各种可以存储程序代码的实体存储介质。The storage medium is a physical, non-transitory storage medium, such as a U disk, a mobile hard disk, a read-only memory (Read-Only Memory, ROM), a magnetic disk, or an optical disk, etc., which can store program codes. medium.
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,上述描述的设备、装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。Those skilled in the art can clearly understand that, for the convenience and conciseness of description, the specific working process of the above-described equipment, device, and unit can refer to the corresponding process in the foregoing method embodiment, which will not be repeated here.
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以权利要求的保护范围为准。The above are only specific implementations of this application, but the protection scope of this application is not limited to this. Anyone familiar with the technical field can easily think of various equivalents within the technical scope disclosed in this application. Modifications or replacements, these modifications or replacements shall be covered within the protection scope of this application. Therefore, the protection scope of this application shall be subject to the protection scope of the claims.

Claims (10)

  1. 一种多车场车辆路径规划方法,包括:A method for multi-depot vehicle path planning includes:
    判断是否接收到客户端发送的路径规划请求;Determine whether the path planning request sent by the client is received;
    若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息;If a path planning request sent by the client is received, the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and each corresponding to the current number of users. The user’s cargo capacity and current user location information for each user corresponding to the current number of users;
    调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集;以及Call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and perform a lot of operations on the vehicle path planning multi-objective optimization model according to the constraint conditions and The evolutionary solution of the goal, the optimal solution set of the path is obtained; and
    将所述路径最优解集发送至客户端。Send the path optimal solution set to the client.
  2. 根据权利要求1所述的多车场车辆路径规划方法,其中,所述车辆路径规划多目标优化模型包括5个优化目标函数,分别记为调度车辆数量优化目标函数minf 1(x)、所有车辆总配送持续时间优化目标函数minf 2(x)、所有车辆总行驶距离优化目标函数minf 3(x)、单个车辆最大容量差优化目标函数minf 4(x)、客户等待时间优化目标函数min f 5(x); The method of multi-depot vehicle path planning according to claim 1, wherein the vehicle path planning multi-objective optimization model includes five optimization objective functions, which are respectively recorded as the number of dispatch vehicles optimization objective function minf 1 (x), and the total number of all vehicles Delivery duration optimization objective function minf 2 (x), total travel distance optimization objective function of all vehicles minf 3 (x), single vehicle maximum capacity difference optimization objective function minf 4 (x), customer waiting time optimization objective function min f 5 ( x);
    Figure PCTCN2020079881-appb-100001
    Figure PCTCN2020079881-appb-100001
    Figure PCTCN2020079881-appb-100002
    Figure PCTCN2020079881-appb-100002
    Figure PCTCN2020079881-appb-100003
    Figure PCTCN2020079881-appb-100003
    Figure PCTCN2020079881-appb-100004
    Figure PCTCN2020079881-appb-100004
    Figure PCTCN2020079881-appb-100005
    Figure PCTCN2020079881-appb-100005
    其中,所述车辆路径规划多目标优化模型中预先设置有R个车场、每一车场有K辆车辆,每一车辆的总货物容量为Q且最大服务总时间为T;R、S、K和Q的取值为正整数;与所述路径规划请求对应的输入数据中包括的当前用户数量记为P;当前用户数量P对应的P个客户节点分别记为节点1至节点P,R个车场对应的车场节点分别记为节点P+1至节点P+R;Wherein, the vehicle path planning multi-objective optimization model is preset with R parking lots, each parking lot has K vehicles, the total cargo capacity of each vehicle is Q and the maximum total service time is T; R, S, K and The value of Q is a positive integer; the current number of users included in the input data corresponding to the path planning request is denoted as P; the P customer nodes corresponding to the current number of users P are denoted as node 1 to node P, and R parking lots The corresponding parking lot nodes are respectively marked as node P+1 to node P+R;
    K r表示第r个车场的实际使用车辆数;rk表示第r个车场的k编号的车辆;d ij表示第i个节点到第j个节点之间的距离,
    Figure PCTCN2020079881-appb-100006
    表示第r个车场的k编号的车辆从第i节点运动第j个节点的行驶时间,s i表示第i节点对应的服务时间,s j表示第j节点对应的服务时间,p i表示第i节点对应的第i个客户的货物容量,
    Figure PCTCN2020079881-appb-100007
    表示第r个车场的k编号的车辆从第i节点运动第j个节点的路径访问状态,
    Figure PCTCN2020079881-appb-100008
    表示第r个车场的k编号的车辆从第i节点运动第j个节点的客户被服务状态。
    K r represents the r-th actual number of yard vehicles; RK k represents a vehicle ID of the r-yard; d ij represents the distance between node i to the j-th node,
    Figure PCTCN2020079881-appb-100006
    Represents the travel time of the k numbered vehicle of the r-th parking lot from the i-th node to the j-th node, s i represents the service time corresponding to the i-th node, s j represents the service time corresponding to the j-th node, and p i represents the i-th node The cargo capacity of the i-th customer corresponding to the node,
    Figure PCTCN2020079881-appb-100007
    Represents the path access state of the k numbered vehicle of the rth parking lot moving from the ith node to the jth node,
    Figure PCTCN2020079881-appb-100008
    It means that the k-numbered vehicle of the r-th parking lot moves from the i-th node to the customer of the j-th node and is served.
  3. 根据权利要求2所述的多车场车辆路径规划方法,其中,所述以所述输 入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集,得到路径最优解集,包括:The method of multi-depot vehicle path planning according to claim 2, wherein the input data is used as the input of the vehicle path planning multi-objective optimization model, and the vehicle path planning is based on the constraint conditions and the multi-objective optimization model. The goal optimization model performs the evolutionary solution of super-multi-objectives, and obtains the path optimal solution set, which includes:
    根据所述约束条件随机生成初始多目标种群;其中,所述初始多目标种群中包括多个个体,每一个体对应所述车辆路径规划多目标优化模型的一个路径输出解,所述初始多目标种群中包括多个个体的总个数记为种群大小N;An initial multi-objective population is randomly generated according to the constraint conditions; wherein the initial multi-objective population includes a plurality of individuals, and each individual corresponds to a path output solution of the vehicle path planning multi-objective optimization model, and the initial multi-objective The total number of individuals in the population is recorded as the population size N;
    获取当前迭代代数,判断所述当前迭代代数是否达到预设的最大迭代代数;Acquiring the current iteration algebra, and judging whether the current iteration algebra reaches the preset maximum iteration algebra;
    若所述当前迭代代数未达到所述最大迭代代数,获取所述初始多目标种群中的理想个体和最差个体;其中,所述理想个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最小目标值,所述最差个体输入至所述车辆路径规划多目标优化模型得到的目标值为初始多目标种群中每个个体对应的目标值中最大目标值;If the current iteration algebra does not reach the maximum iteration algebra, obtain the ideal individual and the worst individual in the initial multi-objective population; wherein the ideal individual is input to the target obtained by the vehicle path planning multi-objective optimization model The value is the smallest target value among the target values corresponding to each individual in the initial multi-objective population, and the worst individual is input to the vehicle path planning multi-objective optimization model to obtain the target value corresponding to each individual in the initial multi-objective population The maximum target value among the target values;
    对所述初始多目标种群进行模拟二进制交叉和多项式变异,得到与所述初始多目标种群有相同个体总个数的子种群;Performing simulated binary crossover and polynomial mutation on the initial multi-target population to obtain a subpopulation with the same total number of individuals as the initial multi-target population;
    将所述初始多目标种群与所述子种群进行合并,得到混合种群;Combining the initial multi-target population and the sub-population to obtain a mixed population;
    将所述混合种群中的个体进行非支配排序,得到非支配解集及多层解集;其中,所述非支配解集记为Q 1,所述多层解集中包括多个解集子集且分别记为Q 2至Q L,其中Q 1至Q L的并集为所述混合种群,Q 1至Q L中任意两个集合的交集为空集,Q 1≥Q 2≥Q 3≥……≥Q LThe individuals in the mixed population are sorted in a non-dominated manner to obtain a non-dominated solution set and a multi-layer solution set; wherein the non-dominated solution set is denoted as Q 1 , and the multi-layer solution set includes a plurality of solution set subsets And respectively denoted as Q 2 to Q L , where the union of Q 1 to Q L is the mixed population, the intersection of any two sets from Q 1 to Q L is the empty set, Q 1 ≥Q 2 ≥Q 3 ≥ ……≥Q L ;
    在所述非支配解集、及多层解集中多个解集子集依序合并从而获取多个集合直至个体的总个数超出所述种群大小N,以组成存档集合;In the non-dominated solution set and the multi-layer solution set, multiple solution set subsets are sequentially merged to obtain multiple sets until the total number of individuals exceeds the population size N to form an archive set;
    根据所述理想个体、所述最差个体将所述存档集合中每一个个体进行归一化处理,得到归一化存档集合;其中,所述归一化存档集合中与所述非支配解集对应的归一化个体集合记为归一化非支配解集;According to the ideal individual and the worst individual, each individual in the archive set is normalized to obtain a normalized archive set; wherein, the normalized archive set and the non-dominated solution set The corresponding normalized individual set is recorded as the normalized non-dominated solution set;
    根据所述归一化非支配解集估计获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面;Obtaining the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the estimation of the normalized non-dominated solution set;
    将所述归一化存档集合中每一个体映射至所述超曲面上,以得到映射集合;其中,所述归一化存档集合中每一个体均对应所述超曲面上的一个映射点,以组成所述映射集合,所述映射集合中对应的映射点的个数记为L maxMapping each individual in the normalized archive set to the hypersurface to obtain a mapping set; wherein each individual in the normalized archive set corresponds to a mapping point on the hypersurface, To form the mapping set, the number of corresponding mapping points in the mapping set is denoted as L max ;
    调用预先存储的目标点适应值算法,获取所述映射集合中每一映射点对应的适应值;Calling a pre-stored target point fitness algorithm to obtain the fitness value corresponding to each mapping point in the mapping set;
    将所述映射集合中各映射点之间在所述超平面中对应的欧氏距离、及所述种群大小N进行聚类,得到聚类结果;其中,所述聚类结果中的聚类簇的总数与所述种群大小N相等;Clustering the corresponding Euclidean distance in the hyperplane between the mapping points in the mapping set and the population size N to obtain a clustering result; wherein, the clustering cluster in the clustering result The total number of is equal to the population size N;
    在所述聚类结果的每一聚类簇中均挑选一个映射点,以组成目标映射点集合;Selecting a mapping point in each cluster of the clustering result to form a target mapping point set;
    获取所述目标映射点集合中各目标映射点对应的个体,以组成当前多目标种群,将所述当前多目标种群作为初始多目标种群;Acquiring an individual corresponding to each target mapping point in the target mapping point set to form a current multi-target population, and use the current multi-target population as an initial multi-target population;
    将所述当前迭代代数加一以作为当前迭代代数,返回执行判断所述当前迭代代数是否达到预设的最大迭代代数的步骤;Adding one to the current iterative algebra as the current iterative algebra, and return to execute the step of judging whether the current iterative algebra reaches the preset maximum iterative algebra;
    若所述当前迭代代数达到所述最大迭代代数,将所述当前多目标种群输出作为路径最优解集。If the current iteration algebra reaches the maximum iteration algebra, output the current multi-target population as a path optimal solution set.
  4. 根据权利要求3所述的多车场车辆路径规划方法,其中,所述对所述初始多目标种群进行模拟二进制交叉和多项式变异,得到与所述初始多目标种群有相同个体总个数的子种群,包括:The method of multi-depot vehicle path planning according to claim 3, wherein the simulated binary crossover and polynomial mutation are performed on the initial multi-target population to obtain a sub-population with the same total number of individuals as the initial multi-target population ,include:
    在所述初始多目标种群中任意挑选两个个体以依次进行二进制交叉,直到生成N个交叉处理后新个体,对N个交叉处理后新个体进行多项式变异,由多项式变异后的新个体组成子种群。Randomly select two individuals in the initial multi-target population to perform binary crossover in sequence, until N cross-processed new individuals are generated, and the N cross-processed new individuals are subjected to polynomial mutation, and the new individuals after polynomial mutation are composed of subgroups. Population.
  5. 根据权利要求3所述的多车场车辆路径规划方法,其中,所述根据所述理想个体、所述最差个体将所述存档集合中每一个个体进行归一化处理,得到归一化存档集合,包括:The method of multi-depot vehicle path planning according to claim 3, wherein, according to the ideal individual and the worst individual, each individual in the archive set is normalized to obtain a normalized archive set ,include:
    根据
    Figure PCTCN2020079881-appb-100009
    将所述存档集合中每一个个体进行归一化处理,得到与所述存档集合中每一个个体对应的归一化个体,以组成归一化存档集合;其中,NA m表示所述存档集合中第m个个体A m对应的归一化个体,A 最差个体表示所述最差个体,A 理想个体表示所述理想个体。
    according to
    Figure PCTCN2020079881-appb-100009
    Perform normalization processing on each individual in the archive set to obtain a normalized individual corresponding to each individual in the archive set to form a normalized archive set; where NA m represents in the archive set the individual over the individual normalization, a represents the worst of the worst individual subject, a represents the individual over the individual m-th corresponding to a m.
  6. 根据权利要求3所述的多车场车辆路径规划方法,其中,所述目标点适应值算法为
    Figure PCTCN2020079881-appb-100010
    f i(x l)表示所述映射集合中第l个映射点对应的第i个目标值,l的取值范围为[1,L max]。
    The method of multi-depot vehicle path planning according to claim 3, wherein the target point fitness algorithm is
    Figure PCTCN2020079881-appb-100010
    f i (x l ) represents the i-th target value corresponding to the l-th mapping point in the mapping set, and the value range of l is [1, L max ].
  7. 根据权利要求3所述的多车场车辆路径规划方法,其中,所述根据所述归一化非支配解集估计获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面,包括:The method of multi-depot vehicle path planning according to claim 3, wherein said obtaining the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front according to the normalized non-dominated solution set estimation comprises :
    将所述归一化非支配解集的各归一化非支配个体中未处于标准化目标空间的归一化非支配个体移除,得到筛选后归一化非支配解集;其中,位于所述标准化目标空间中各归一化非支配个体对应的各目标值均未超过1;Remove the normalized non-dominated individuals that are not in the standardized target space among the normalized non-dominated individuals in the normalized non-dominated solution set, and obtain the normalized non-dominated solution set after screening; wherein, located in the Each target value corresponding to each normalized non-dominated individual in the standardized target space does not exceed 1;
    获取所述筛选后归一化非支配解集中每一非支配个体,分别记为B 1至B n;其中,n的取值与所述筛选后归一化非支配解集中归一化非支配个体的总个数相同; Obtain each non-dominated individual in the normalized non-dominated solution set after the screening, respectively denoted as B 1 to B n ; where the value of n is the same as the normalized non-dominated in the normalized non-dominated solution set after the screening The total number of individuals is the same;
    获取非支配个体B i到目标超平面对应的超平面距离D i;其中,i的取值范围为[1,n],目标超平面为f 1(x)+f 2(x)+f 3(x)+f 4(x)+f 5(x)=1; Obtain the hyperplane distance D i corresponding to the target hyperplane from the non-dominated individual B i ; where the value range of i is [1,n], and the target hyperplane is f 1 (x)+f 2 (x)+f 3 (x)+f 4 (x)+f 5 (x)=1;
    获取超平面距离D 1至D n对应的超平面距离平均值和超平面距离标准差;其中,超平面距离D 1至D n对应的超平面距离平均值记为D avg,超平面距离D 1至D n对应的超平面距离标准差记为D stdGet hyperplane distance D 1 to D n corresponding hyperplane distance from average and standard deviation hyperplane; wherein hyperplane distance from the hyperplane D 1 to D n corresponding to an average value referred to as D avg, the distance D 1 hyperplane The standard deviation of the hyperplane distance corresponding to D n is denoted as D std ;
    根据所述超平面距离平均值与所述超平面距离标准差之商进行范数运算,得到对应的变异系数;其中,所述变异系数记为cv;Perform a norm operation according to the quotient of the average value of the hyperplane distance and the standard deviation of the hyperplane distance to obtain the corresponding coefficient of variation; wherein the coefficient of variation is denoted as cv;
    获取所述目标超平面的曲率为2对应的第一超曲面,和所述目标超平面的曲率为0.5对应的第二超曲面,根据目标超平面、第一超曲面、第二超曲面及预设的曲率确定策略,获取所述归一化非支配解集对应的当前曲率;其中,所述 曲率确定策略为
    Figure PCTCN2020079881-appb-100011
    d(2.0)表示第一超曲面的峰值点到所述目标超平面的距离,d(0.5)表示第二超曲面的峰值点到所述目标超平面的距离;
    Obtain a first hypersurface corresponding to a curvature of 2 of the target hyperplane, and a second hypersurface corresponding to a curvature of 0.5 of the target hyperplane, according to the target hyperplane, the first hypersurface, the second hypersurface, and the prediction Set the curvature determination strategy to obtain the current curvature corresponding to the normalized non-dominated solution set; wherein, the curvature determination strategy is
    Figure PCTCN2020079881-appb-100011
    d(2.0) represents the distance from the peak point of the first hypersurface to the target hyperplane, and d(0.5) represents the distance from the peak point of the second hypersurface to the target hyperplane;
    根据所述变异系数对所述当前曲率进行调整,得到调整后曲率;其中,若所述变异系数小于0.1,将所述当前曲率的取值调整为1,以使调整后曲率取值为1;若所述变异系数大于或等于0.1,将所述当前曲率的取值保持不变,以使得调整后曲率等于所述当前曲率;Adjust the current curvature according to the coefficient of variation to obtain an adjusted curvature; wherein, if the coefficient of variation is less than 0.1, adjust the value of the current curvature to 1, so that the adjusted curvature takes the value of 1; If the coefficient of variation is greater than or equal to 0.1, keeping the value of the current curvature unchanged, so that the adjusted curvature is equal to the current curvature;
    根据所述调整后曲率对应获取帕累托前沿的形状,及帕累托前沿的形状对应的超曲面。According to the adjusted curvature, the shape of the Pareto front and the hypersurface corresponding to the shape of the Pareto front are correspondingly obtained.
  8. 一种多车场车辆路径规划装置,包括:A multi-parking vehicle path planning device, including:
    路径规划请求检测单元,用于判断是否接收到客户端发送的路径规划请求;The path planning request detection unit is used to determine whether the path planning request sent by the client is received;
    数据条件获取单元,用于若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息;The data condition obtaining unit is configured to obtain input data and constraint conditions corresponding to the path planning request if the path planning request sent by the client is received; wherein, the input data corresponding to the path planning request includes the current number of users, The cargo capacity of each user corresponding to the current number of users, and the current user location information of each user corresponding to the current number of users;
    路径最优解集获取单元,用于调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集;The path optimal solution set acquisition unit is used to call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and compare the constraints to the The multi-objective optimization model of vehicle path planning performs the evolutionary solution of super multi-objectives and obtains the optimal solution set of the path;
    最优解集发送单元,用于将所述路径最优解集发送至客户端。The optimal solution set sending unit is configured to send the path optimal solution set to the client.
  9. 一种计算机设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时实现以下步骤:A computer device includes a memory, a processor, and a computer program stored on the memory and capable of running on the processor, wherein the processor implements the following steps when the processor executes the computer program:
    判断是否接收到客户端发送的路径规划请求;Determine whether the path planning request sent by the client is received;
    若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息;If a path planning request sent by the client is received, the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and each corresponding to the current number of users. The user’s cargo capacity and current user location information for each user corresponding to the current number of users;
    调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集;以及Call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and perform a lot of operations on the vehicle path planning multi-objective optimization model according to the constraint conditions and The evolutionary solution of the goal, the optimal solution set of the path is obtained; and
    将所述路径最优解集发送至客户端。Send the path optimal solution set to the client.
  10. 一种计算机可读存储介质,其中,所述计算机可读存储介质存储有计算机程序,所述计算机程序当被处理器执行时使所述处理器执行以下操作:A computer-readable storage medium, wherein the computer-readable storage medium stores a computer program that, when executed by a processor, causes the processor to perform the following operations:
    判断是否接收到客户端发送的路径规划请求;Determine whether the path planning request sent by the client is received;
    若接收到客户端发送的路径规划请求,获取与所述路径规划请求对应的输入数据和约束条件;其中,与所述路径规划请求对应的输入数据包括当前用户 数量、当前用户数量对应的每一用户的货物容量、当前用户数量对应的每一用户的当前用户位置信息;If the path planning request sent by the client is received, the input data and constraint conditions corresponding to the path planning request are obtained; wherein, the input data corresponding to the path planning request includes the current number of users and each corresponding to the current number of users. The user’s cargo capacity and current user location information for each user corresponding to the current number of users;
    调用预先存储的车辆路径规划多目标优化模型,以所述输入数据为所述车辆路径规划多目标优化模型的输入,并根据所述约束条件和对所述车辆路径规划多目标优化模型进行超多目标的进化求解,得到路径最优解集;以及Call a pre-stored vehicle path planning multi-objective optimization model, use the input data as the input of the vehicle path planning multi-objective optimization model, and perform a lot of operations on the vehicle path planning multi-objective optimization model according to the constraint conditions and The evolutionary solution of the goal, the optimal solution set of the path is obtained; and
    将所述路径最优解集发送至客户端。Send the path optimal solution set to the client.
PCT/CN2020/079881 2020-01-15 2020-03-18 Multi-depot vehicle routing method, apparatus, computer device and storage medium WO2021142917A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202010042242.2 2020-01-15
CN202010042242.2A CN111260129B (en) 2020-01-15 2020-01-15 Multi-yard vehicle path planning method and device, computer equipment and storage medium

Publications (1)

Publication Number Publication Date
WO2021142917A1 true WO2021142917A1 (en) 2021-07-22

Family

ID=70948839

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/079881 WO2021142917A1 (en) 2020-01-15 2020-03-18 Multi-depot vehicle routing method, apparatus, computer device and storage medium

Country Status (2)

Country Link
CN (1) CN111260129B (en)
WO (1) WO2021142917A1 (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113626954A (en) * 2021-08-17 2021-11-09 中国地质大学(武汉) Multi-target information processing method and system based on decomposition, computer equipment and terminal
CN113705921A (en) * 2021-09-03 2021-11-26 厦门闽江智慧科技有限公司 Electric vehicle dynamic path planning optimization solving method based on hybrid charging strategy
CN113780676A (en) * 2021-09-23 2021-12-10 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN113887083A (en) * 2021-10-29 2022-01-04 北京明略软件系统有限公司 Air conditioner scheduling optimization method and system, computer equipment and storage medium
CN114065512A (en) * 2021-11-16 2022-02-18 吉林大学 Induction type magnetometer coil parameter multi-objective optimization method based on NSGA-II algorithm
CN114169832A (en) * 2021-12-15 2022-03-11 中国建设银行股份有限公司 Logistics distribution scheduling method and device, storage medium and processor
CN114594744A (en) * 2022-03-14 2022-06-07 武汉理工大学 Distributed factory production and distribution integrated scheduling method and system
CN114662943A (en) * 2022-03-28 2022-06-24 湘潭大学 Strip mine truck scheduling method based on multi-target genetic algorithm
CN114781706A (en) * 2022-04-11 2022-07-22 东南大学 Logistics vehicle distribution scheduling method, vehicle, electronic equipment and storage medium
CN114819249A (en) * 2022-02-15 2022-07-29 南京邮电大学 Community material vehicle path optimization method based on improved lion group algorithm under epidemic situation
CN114997659A (en) * 2022-06-06 2022-09-02 武汉理工大学 Resource scheduling model construction method and system based on dynamic multi-objective optimization
CN115099118A (en) * 2022-06-24 2022-09-23 中国人民解放军海军工程大学 NSGA III-based high-dimensional multi-target joint parallel simulation optimization method
CN115345549A (en) * 2022-08-11 2022-11-15 上海维祥信息技术有限公司 Vehicle path adjusting method and system combined with loading scheme
CN115576343A (en) * 2022-09-08 2023-01-06 浙江财经大学 Multi-target vehicle path optimization method combining unmanned aerial vehicle distribution
CN115994459A (en) * 2023-03-24 2023-04-21 国网湖北省电力有限公司经济技术研究院 Method, system and equipment for optimizing construction time sequence of mass power grid infrastructure projects
CN116013085A (en) * 2023-02-03 2023-04-25 以萨技术股份有限公司 Data processing system for obtaining target traffic flow
CN116382099A (en) * 2023-06-02 2023-07-04 上海数字大脑科技研究院有限公司 Robot path scheduling planning method and system
CN116739467A (en) * 2023-06-16 2023-09-12 暨南大学 Circulating delivery path planning method, device and medium based on cloud edge cooperation
CN117875825A (en) * 2024-03-11 2024-04-12 万联易达物流科技有限公司 Method for matching cargo sources of transportation circulating line

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112215418B (en) * 2020-10-12 2024-04-16 苏州大学 Method, device and storage medium for automatically generating washing inspection arrangement scheme of multi-rail vehicle
CN112270047B (en) * 2020-11-16 2023-09-29 中山大学 Urban vehicle path optimization method based on data-driven group intelligent calculation
CN114493399B (en) * 2021-12-17 2024-04-09 重庆特斯联智慧科技股份有限公司 Logistics robot pickup order-oriented path planning method and system
CN115345550B (en) * 2022-08-11 2023-10-24 上海维祥信息技术有限公司 Vehicle path planning method and system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001084480A (en) * 1999-09-10 2001-03-30 Hitachi Software Eng Co Ltd Method and device for making delivery plan
JP2001283136A (en) * 2000-03-29 2001-10-12 Hitachi Software Eng Co Ltd Method and system for preparing delivery plan
CN101673382A (en) * 2009-10-21 2010-03-17 北京交通大学 Combined optimization method for agricultural chain-operation logistics delivering and loading-distribution
CN107578199A (en) * 2017-08-21 2018-01-12 南京航空航天大学 A kind of method for solving two dimension and loading constraint logistics vehicle dispatching problem
CN108510227A (en) * 2018-03-23 2018-09-07 东华大学 A kind of real-time planning system of vehicle-mounted logistics distribution based on machine learning

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106503836B (en) * 2016-10-09 2020-06-16 电子科技大学 Multi-objective optimized pure electric vehicle logistics distribution optimized scheduling method
CN106651043A (en) * 2016-12-28 2017-05-10 中山大学 Intelligent algorithm for solving a multi-objective MDVRPTW (Multi-Depot Vehicle Routing Problem with Time Window)
CN108416482B (en) * 2018-04-03 2021-08-03 浙江工业大学 Logistics distribution path planning method based on regional multiple strategies
CN109948855B (en) * 2019-03-22 2021-08-31 杭州电子科技大学 Heterogeneous hazardous chemical transportation path planning method with time window

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001084480A (en) * 1999-09-10 2001-03-30 Hitachi Software Eng Co Ltd Method and device for making delivery plan
JP2001283136A (en) * 2000-03-29 2001-10-12 Hitachi Software Eng Co Ltd Method and system for preparing delivery plan
CN101673382A (en) * 2009-10-21 2010-03-17 北京交通大学 Combined optimization method for agricultural chain-operation logistics delivering and loading-distribution
CN107578199A (en) * 2017-08-21 2018-01-12 南京航空航天大学 A kind of method for solving two dimension and loading constraint logistics vehicle dispatching problem
CN108510227A (en) * 2018-03-23 2018-09-07 东华大学 A kind of real-time planning system of vehicle-mounted logistics distribution based on machine learning

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113626954A (en) * 2021-08-17 2021-11-09 中国地质大学(武汉) Multi-target information processing method and system based on decomposition, computer equipment and terminal
CN113626954B (en) * 2021-08-17 2023-08-18 中国地质大学(武汉) Multi-target information processing method, system, computer equipment and terminal based on decomposition
CN113705921A (en) * 2021-09-03 2021-11-26 厦门闽江智慧科技有限公司 Electric vehicle dynamic path planning optimization solving method based on hybrid charging strategy
CN113705921B (en) * 2021-09-03 2024-02-27 厦门闽江智慧科技有限公司 Electric vehicle dynamic path planning optimization method based on hybrid charging strategy
CN113780676A (en) * 2021-09-23 2021-12-10 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN113780676B (en) * 2021-09-23 2023-06-23 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN113887083A (en) * 2021-10-29 2022-01-04 北京明略软件系统有限公司 Air conditioner scheduling optimization method and system, computer equipment and storage medium
CN114065512A (en) * 2021-11-16 2022-02-18 吉林大学 Induction type magnetometer coil parameter multi-objective optimization method based on NSGA-II algorithm
CN114169832A (en) * 2021-12-15 2022-03-11 中国建设银行股份有限公司 Logistics distribution scheduling method and device, storage medium and processor
CN114819249A (en) * 2022-02-15 2022-07-29 南京邮电大学 Community material vehicle path optimization method based on improved lion group algorithm under epidemic situation
CN114594744A (en) * 2022-03-14 2022-06-07 武汉理工大学 Distributed factory production and distribution integrated scheduling method and system
CN114594744B (en) * 2022-03-14 2023-11-28 武汉理工大学 Distributed factory production distribution integrated scheduling method and system
CN114662943B (en) * 2022-03-28 2024-05-31 湘潭大学 Strip mine truck scheduling method based on multi-target genetic algorithm
CN114662943A (en) * 2022-03-28 2022-06-24 湘潭大学 Strip mine truck scheduling method based on multi-target genetic algorithm
CN114781706A (en) * 2022-04-11 2022-07-22 东南大学 Logistics vehicle distribution scheduling method, vehicle, electronic equipment and storage medium
CN114781706B (en) * 2022-04-11 2024-05-03 东南大学 Logistics vehicle distribution scheduling method, vehicle, electronic equipment and storage medium
CN114997659A (en) * 2022-06-06 2022-09-02 武汉理工大学 Resource scheduling model construction method and system based on dynamic multi-objective optimization
CN115099118A (en) * 2022-06-24 2022-09-23 中国人民解放军海军工程大学 NSGA III-based high-dimensional multi-target joint parallel simulation optimization method
CN115099118B (en) * 2022-06-24 2024-04-09 中国人民解放军海军工程大学 NSGA III-based high-dimensional multi-objective joint parallel simulation optimization method
CN115345549A (en) * 2022-08-11 2022-11-15 上海维祥信息技术有限公司 Vehicle path adjusting method and system combined with loading scheme
CN115345549B (en) * 2022-08-11 2024-01-23 上海维祥信息技术有限公司 Vehicle path adjustment method and system combined with loading scheme
CN115576343A (en) * 2022-09-08 2023-01-06 浙江财经大学 Multi-target vehicle path optimization method combining unmanned aerial vehicle distribution
CN115576343B (en) * 2022-09-08 2024-06-07 浙江财经大学 Multi-target vehicle path optimization method combined with unmanned aerial vehicle distribution
CN116013085A (en) * 2023-02-03 2023-04-25 以萨技术股份有限公司 Data processing system for obtaining target traffic flow
CN116013085B (en) * 2023-02-03 2023-09-22 以萨技术股份有限公司 Data processing system for obtaining target traffic flow
CN115994459B (en) * 2023-03-24 2023-06-13 国网湖北省电力有限公司经济技术研究院 Method, system and equipment for optimizing construction time sequence of mass power grid infrastructure projects
CN115994459A (en) * 2023-03-24 2023-04-21 国网湖北省电力有限公司经济技术研究院 Method, system and equipment for optimizing construction time sequence of mass power grid infrastructure projects
CN116382099A (en) * 2023-06-02 2023-07-04 上海数字大脑科技研究院有限公司 Robot path scheduling planning method and system
CN116739467B (en) * 2023-06-16 2024-02-02 暨南大学 Circulating delivery path planning method, device and medium based on cloud edge cooperation
CN116739467A (en) * 2023-06-16 2023-09-12 暨南大学 Circulating delivery path planning method, device and medium based on cloud edge cooperation
CN117875825A (en) * 2024-03-11 2024-04-12 万联易达物流科技有限公司 Method for matching cargo sources of transportation circulating line
CN117875825B (en) * 2024-03-11 2024-06-04 万联易达物流科技有限公司 Method for matching cargo sources of transportation circulating line

Also Published As

Publication number Publication date
CN111260129A (en) 2020-06-09
CN111260129B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
WO2021142917A1 (en) Multi-depot vehicle routing method, apparatus, computer device and storage medium
US8117609B2 (en) System and method for optimizing changes of data sets
CN112288347A (en) Method, device, server and storage medium for determining route of cold chain distribution
WO2021142915A1 (en) Multi-objective flow shop scheduling method and apparatus, and computer device and storage medium
US10997540B2 (en) System and method for matching resource capacity with client resource needs
CN106910065B (en) Data processing method, device and system for calculating settlement amount based on multiple transactions
CN110099112A (en) Date storage method, device, medium and terminal device based on point to point network
CN111044062B (en) Path planning and recommending method and device
CN112751785B (en) Method and device for sending pending request, computer equipment and storage medium
CN112966054A (en) Enterprise graph node relation-based ethnic group division method and computer equipment
CN107659982B (en) Wireless network access point classification method and device
CN111858585A (en) Block chain strategy processing device, computer readable storage medium and terminal equipment
CN106681803A (en) Task scheduling method and server
CN112738213B (en) Task demand response method, device and system based on block chain and storage medium
CN115016911A (en) Task arrangement method, device, equipment and medium for large-scale federal learning
Han et al. Opportunistic coded distributed computing: An evolutionary game approach
CN113837677A (en) Method, device and equipment for determining logistics line generation strategy
CN111291792A (en) Traffic data type integrated classification method and device based on double evolution
CN115330559B (en) Power distribution network elasticity evaluation method and device based on information data time-space coordination
CN114722085B (en) Service combination method, system, equipment and medium for meeting user demand constraint
US20230054360A1 (en) Machine learning consolidation
Lu Research on the Application of Multidimensional Cluster Analysis in Customer Information
CN115422303A (en) Intelligent relationship determination method and device
CN118101749A (en) Information pushing method, device, equipment and medium
CN115883667A (en) Resource transfer routing method, device, computer equipment and storage medium

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20913883

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 11/11/2022)

122 Ep: pct application non-entry in european phase

Ref document number: 20913883

Country of ref document: EP

Kind code of ref document: A1