Nothing Special   »   [go: up one dir, main page]

CN117132011A - Intercity travel vehicle route determination method, system, electronic device and medium - Google Patents

Intercity travel vehicle route determination method, system, electronic device and medium Download PDF

Info

Publication number
CN117132011A
CN117132011A CN202311209584.9A CN202311209584A CN117132011A CN 117132011 A CN117132011 A CN 117132011A CN 202311209584 A CN202311209584 A CN 202311209584A CN 117132011 A CN117132011 A CN 117132011A
Authority
CN
China
Prior art keywords
travel
vehicle
order
updated
solution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202311209584.9A
Other languages
Chinese (zh)
Inventor
蔡奕侨
林子帆
王依婷
柳培忠
喻小光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
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 Huaqiao University filed Critical Huaqiao University
Priority to CN202311209584.9A priority Critical patent/CN117132011A/en
Publication of CN117132011A publication Critical patent/CN117132011A/en
Pending legal-status Critical Current

Links

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
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • 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

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method, a system, electronic equipment and a medium for determining an inter-city travel vehicle path, which relate to the fields of intelligent optimization algorithms and vehicle scheduling, wherein the method for determining the path comprises the following steps: constructing an inter-city travel vehicle path model by taking the maximized average number of passengers per trip, the minimized number of vehicles, the minimized total travel distance of the vehicles and the minimized total waiting time of the passengers as objective functions and taking the vehicle passenger capacity, the service quality, the time constraint and the safety constraint as constraint conditions; acquiring an inter-city order to be travelled; and determining the travel vehicles corresponding to each travel order according to the inter-city travel vehicle path model. The invention can provide a path planning scheme for meeting a plurality of requirements for the travel of the inter-city network about vehicles.

Description

一种城际出行车辆路径确定方法、系统、电子设备及介质Method, system, electronic device and medium for determining intercity travel vehicle path

技术领域Technical Field

本发明涉及智能优化算法和车辆调度领域,特别是涉及一种城际出行车辆路径确定方法、系统、电子设备及介质。The present invention relates to the field of intelligent optimization algorithms and vehicle scheduling, and in particular to a method, system, electronic equipment and medium for determining a vehicle path for inter-city travel.

背景技术Background Art

随着智能手机和应用程序的普及,越来越多的公司开发搭建由专业驾驶员和商用车辆组成的乘车共享平台,以互联网技术为基础,通过整合供需信息,合理匹配乘客和司机,高效安全规划出行路线,实现非传统出租车预订服务的商业活动。这项新的移动服务经历了快速增长,并以指数级增长,为跨越66个国家的大多数大都市地区提供服务。由于这些优势,网约车在中国迅速普及。城际出行可以充分利用车辆空位资源来满足无车用户的出行需求、避免拥挤、节约出行时间,同时减少了私家车车主的车辆出行成本、提高了车辆利用率、缓解了交通拥堵、节约资源以及保护环境,因此城际出行具有重要的现实意义。With the popularity of smartphones and applications, more and more companies are developing ride-sharing platforms consisting of professional drivers and commercial vehicles. Based on Internet technology, they integrate supply and demand information, reasonably match passengers and drivers, and plan travel routes efficiently and safely to realize commercial activities of non-traditional taxi booking services. This new mobile service has experienced rapid growth and has grown exponentially to provide services to most metropolitan areas across 66 countries. Due to these advantages, online ride-hailing has quickly become popular in China. Intercity travel can make full use of vehicle vacancy resources to meet the travel needs of car-free users, avoid congestion, and save travel time. At the same time, it reduces the vehicle travel costs of private car owners, improves vehicle utilization, alleviates traffic congestion, saves resources, and protects the environment. Therefore, intercity travel has important practical significance.

而当前面向城际出行的车辆路径问题却很少人研究,目前对于拼车的应用场景大多为城市内的拼车出行,市内出行有四大特点:出行需求高、出行频率高、出行距离短、出行方便。这些特征使得市内出行通常是实时进行的。平台为了保证服务质量需要在很短的时间内响应乘客的请求。而与市内出行不同,城际出行具有以下特点:出行需求低、出行频率低、出行距离长、出行选择少。这些特点使得城际出行通常是高度计划的,尤其是非自驾出行。由于城际出行是长途出行,大部分出行时间都花在两个城市之间的高速公路上。乘客不会对市内出行时间的微小变化敏感。因此这一特征不允许平台以实时方式工作。与优化司机和乘客匹配以及定价的市内网约车平台不同,平台主要优化日常运营中的车辆路线。在现实世界中,乘客需要提前几个小时向平台提供出行信息。However, few people have studied the vehicle routing problem for intercity travel. At present, the application scenarios of carpooling are mostly carpooling travel within the city. Intracity travel has four characteristics: high travel demand, high travel frequency, short travel distance, and convenient travel. These characteristics make intracity travel usually carried out in real time. In order to ensure the quality of service, the platform needs to respond to passengers' requests in a very short time. Unlike intracity travel, intercity travel has the following characteristics: low travel demand, low travel frequency, long travel distance, and few travel options. These characteristics make intercity travel usually highly planned, especially non-self-driving travel. Since intercity travel is a long-distance trip, most of the travel time is spent on the highway between two cities. Passengers are not sensitive to small changes in intracity travel time. Therefore, this feature does not allow the platform to work in real time. Unlike intracity online car-hailing platforms that optimize driver and passenger matching and pricing, the platform mainly optimizes vehicle routes in daily operations. In the real world, passengers need to provide travel information to the platform several hours in advance.

因此由于当前研究中并没有模型能充分代表面向城际出行的多目标车辆路径问题,鉴于面向城际出行的车辆路径问题本质上是一个多目标优化问题,而多目标优化问题是优化问题中的核心问题之一,许多学者对其进行研究并提出了多种元启发式方法的变种用于解决不同情况下的多目标优化问题。而对于大多数现有方法来说,有效地解决多个冲突目标的优化问题仍然是一个巨大的挑战。同时,在该问题中,不同的目标在物理意义、规模和优化难度上具有不同的特性,因此,针对一些用常规邻域搜索操作难以优化的目标,需要设计针对目标的邻域搜索操作。故如何设计针对城际出行车辆路径问题的高效算法框架具有较高的理论意义和应用价值。Therefore, since there is no model in the current research that can fully represent the multi-objective vehicle routing problem for intercity travel, and given that the vehicle routing problem for intercity travel is essentially a multi-objective optimization problem, and the multi-objective optimization problem is one of the core problems in optimization problems, many scholars have studied it and proposed a variety of meta-heuristic method variants to solve multi-objective optimization problems in different situations. For most existing methods, effectively solving the optimization problem of multiple conflicting objectives is still a huge challenge. At the same time, in this problem, different objectives have different characteristics in physical meaning, scale and optimization difficulty. Therefore, for some objectives that are difficult to optimize using conventional neighborhood search operations, it is necessary to design neighborhood search operations for the objectives. Therefore, how to design an efficient algorithm framework for the intercity travel vehicle routing problem has high theoretical significance and application value.

发明内容Summary of the invention

本发明的目的是提供一种城际出行车辆路径确定方法、系统、电子设备及介质,能够为城际网约车出行提供满足多个需求的路径规划方案。The purpose of the present invention is to provide a method, system, electronic device and medium for determining a vehicle path for inter-city travel, which can provide a path planning solution that meets multiple needs for inter-city online car-hailing travel.

为实现上述目的,本发明提供了如下方案:To achieve the above object, the present invention provides the following solutions:

一种城际出行车辆路径确定方法,所述确定方法包括:A method for determining a path for an inter-city travel vehicle, the method comprising:

以最大化平均每趟出行的乘客数、最小化车辆数目、最小化车辆总行驶距离和最小化乘客总等待时间为目标函数,以车辆载客量、服务质量、时间约束和安全约束为约束条件,构建城际出行车辆路径模型;Taking maximizing the average number of passengers per trip, minimizing the number of vehicles, minimizing the total distance traveled by vehicles and minimizing the total waiting time of passengers as the objective function, and taking vehicle capacity, service quality, time constraints and safety constraints as constraints, an inter-city travel vehicle path model is constructed.

获取待出行城际订单;所述待出行城际订单包括多个出行订单;各所述出行订单包括乘车人数、出发时间、出发地点和目的地点;Obtaining intercity travel orders to be traveled; the intercity travel orders to be traveled include multiple travel orders; each of the travel orders includes the number of passengers, departure time, departure location and destination point;

根据所述城际出行车辆路径模型,确定各所述出行订单对应的出行车辆;其中,一个所述出行车辆和由所述出行车辆完成的多个所述出行订单构成一个出行解;多个所述出行解构成出行外部存档;将所述出行外部存档作为所述待出行城际订单的路径规划方案。According to the inter-city travel vehicle path model, the travel vehicles corresponding to each travel order are determined; wherein, one travel vehicle and multiple travel orders completed by the travel vehicle constitute a travel solution; multiple travel solutions constitute a travel external archive; and the travel external archive is used as the path planning solution for the inter-city order to be traveled.

可选地,所述目标函数为:Optionally, the objective function is:

f2=|Mk|f 2 = |M k |

其中,目标f1表示最大化平均每趟行程的乘客数,Mt、Mk分别表示行程t、车辆k的数量,Qt表示行程t的上座率,目标f2表示最小化车辆数目,目标f3表示最小化车辆总行驶距离,Mt表示行程t的数量,lt表示行程t的总行驶距离,用来判断行程t是否安排至车辆k,目标f4表示最小化乘客总等待时间。Among them, the goal f1 represents maximizing the average number of passengers per trip, Mt and Mk represent the number of trips t and vehicles k respectively, Qt represents the occupancy rate of trip t, the goal f2 represents minimizing the number of vehicles, the goal f3 represents minimizing the total driving distance of vehicles, Mt represents the number of trips t, lt represents the total driving distance of trip t, It is used to determine whether trip t is arranged to vehicle k. Objective f 4 represents minimizing the total waiting time of passengers.

可选地,根据所述城际出行车辆路径模型,确定各所述出行订单对应的出行车辆,具体包括:Optionally, determining the travel vehicles corresponding to each of the travel orders according to the inter-city travel vehicle path model specifically includes:

根据所述约束条件和所述目标函数,确定与各所述出行订单对应的多个待出行行程;其中,各所述待出行行程包括至少一个所述出行订单且各所述待出行行程满足所述约束条件和所述目标函数;According to the constraint conditions and the objective function, a plurality of itineraries to be traveled corresponding to each of the travel orders are determined; wherein each of the itineraries to be traveled includes at least one of the travel orders and each of the itineraries to be traveled satisfies the constraint conditions and the objective function;

根据所述约束条件和所述目标函数,确定与各所述待出行行程对应的车辆;其中,一个所述车辆和由所述车辆完成的多个所述待出行行程构成一个解;多个所述解构成外部存档;According to the constraint conditions and the objective function, determine the vehicle corresponding to each of the to-be-traveled trips; wherein one of the vehicles and a plurality of the to-be-traveled trips completed by the vehicle constitute a solution; and a plurality of the solutions constitute an external archive;

对各所述解应用邻域操作算子,得到各所述解对应的搜索结果;Applying a neighborhood operation operator to each of the solutions to obtain search results corresponding to each of the solutions;

对各所述解和各所述解对应的搜索结果应用ε占优存档更新策略,得到出行外部存档。The ε-dominant archive update strategy is applied to each of the solutions and the search results corresponding to each of the solutions to obtain a travel external archive.

可选地,根据所述约束条件和所述目标函数,确定与各所述出行订单对应的多个待出行行程,具体包括:Optionally, determining a plurality of to-be-traveled trips corresponding to each of the travel orders according to the constraint conditions and the objective function specifically includes:

将预处理订单集中出发时间最早的出行订单作为初始行程的基准订单;所述初始行程至少包括一个所述出行订单;所述预处理订单集包括出发地点相同的多个所述出行订单;The travel order with the earliest departure time in the pre-processed order set is used as the benchmark order for the initial trip; the initial trip includes at least one travel order; the pre-processed order set includes multiple travel orders with the same departure location;

从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列;所述待处理订单集为所述预处理订单集去掉所述基准订单后的剩余订单;所述订单队列为从所述待处理订单集中去掉一个所述预设预处理订单后的剩余订单;A preset pre-processed order is selected from the pending order set and added to the initial itinerary to obtain an updated initial itinerary and an order queue; the pending order set is the remaining orders after the pre-processed order set is removed from the benchmark order; the order queue is the remaining orders after one of the preset pre-processed orders is removed from the pending order set;

判断所述更新后的初始行程是否满足所述约束条件;Determining whether the updated initial itinerary satisfies the constraint condition;

当所述更新后的初始行程不满足所述约束条件时,将所述预设预处理订单从所述更新后的初始行程中删除,并返回步骤“从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列”继续执行;When the updated initial itinerary does not satisfy the constraint condition, the preset pre-processed order is deleted from the updated initial itinerary, and the process returns to the step of "selecting a preset pre-processed order from the pending order set and adding it to the initial itinerary to obtain an updated initial itinerary and order queue" for further execution;

当所述更新后的初始行程满足所述约束条件时,判断所述更新后的初始行程的乘车人数是否达到所述车辆载客量;When the updated initial trip satisfies the constraint condition, determining whether the number of passengers in the updated initial trip reaches the vehicle passenger capacity;

当所述更新后的初始行程的乘车人数达到所述车辆载客量时,所述更新后的初始行程作为所述最终的初始行程,并将所述最终的初始行程中的出行订单从所述预处理订单集中删除,得到更新后的预处理订单集;When the number of passengers in the updated initial trip reaches the passenger capacity of the vehicle, the updated initial trip is used as the final initial trip, and the travel orders in the final initial trip are deleted from the pre-processed order set to obtain an updated pre-processed order set;

当所述更新后的初始行程的乘车人数未达到所述车辆载客量时,判断所述订单队列中出行订单的数量是否为0;When the number of passengers in the updated initial trip does not reach the passenger capacity of the vehicle, determining whether the number of travel orders in the order queue is 0;

当所述订单队列中出行订单的数量不为0时,将所述订单队列作为待处理订单集,并返回步骤“从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列”;When the number of travel orders in the order queue is not 0, the order queue is used as a pending order set, and the process returns to the step of "selecting a preset pre-processed order from the pending order set and adding it to the initial trip to obtain an updated initial trip and order queue";

当所述订单队列中出行订单的数量为0时,返回步骤“所述更新后的初始行程作为所述最终的初始行程,并将所述最终的初始行程中的出行订单从所述预处理订单集中删除,得到更新后的预处理订单集”;When the number of travel orders in the order queue is 0, return to the step of "taking the updated initial itinerary as the final initial itinerary, and deleting the travel orders in the final initial itinerary from the pre-processed order set, to obtain an updated pre-processed order set";

当所述更新后的预处理订单集中出行订单的数量不为0时,将所述更新后的预处理订单集作为预处理订单集,并返回步骤“将预处理订单集中出发时间最早的出行订单作为初始行程的基准订单”继续执行;When the number of travel orders in the updated pre-processed order set is not 0, the updated pre-processed order set is used as the pre-processed order set, and the process returns to the step of "using the travel order with the earliest departure time in the pre-processed order set as the benchmark order for the initial trip" to continue execution;

当所述更新后的预处理订单集中出行订单的数量为0时,确定所述最终的初始行程为与各所述出行订单对应的多个待出行行程。When the number of travel orders in the updated pre-processed order set is 0, the final initial itinerary is determined to be a plurality of to-be-traveled itineraries corresponding to each of the travel orders.

可选地,根据所述约束条件和所述目标函数,确定与各所述待出行行程对应的车辆,具体包括:Optionally, determining the vehicle corresponding to each of the to-be-traveled trips according to the constraint condition and the objective function specifically includes:

将预处理行程集中出发时间最早的待出行行程作为车辆的基准行程;所述车辆至少包括一个所述待出行行程;所述预处理行程集包括出发地点相同的多个所述待出行行程;The trip to be traveled with the earliest departure time in the pre-processed trip set is used as the reference trip of the vehicle; the vehicle includes at least one trip to be traveled; the pre-processed trip set includes multiple trips to be traveled with the same departure point;

从待处理行程集中选择一个预设预处理行程加入到所述初始行程中,得到更新后的车辆和行程队列;所述待处理行程集为所述预处理行程集去掉所述基准行程后的剩余行程;所述行程队列为从所述待处理行程集中去掉一个所述预设预处理行程后的剩余行程;A preset pre-processed trip is selected from the to-be-processed trip set and added to the initial trip to obtain an updated vehicle and trip queue; the to-be-processed trip set is the remaining trip after the pre-processed trip set is removed from the reference trip; the trip queue is the remaining trip after one of the preset pre-processed trips is removed from the to-be-processed trip set;

判断所述更新后的车辆是否满足所述约束条件;Determining whether the updated vehicle satisfies the constraint condition;

当所述更新后的车辆不满足所述约束条件时,将所述预设预处理行程从所述更新后的车辆中删除,并返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述车辆中,得到更新后的车辆和行程队列”继续执行;When the updated vehicle does not satisfy the constraint condition, the preset pre-processed trip is deleted from the updated vehicle, and the process returns to the step of "selecting a preset pre-processed trip from the set of trips to be processed and adding it to the vehicle to obtain an updated vehicle and trip queue" to continue the execution;

当所述更新后的车辆满足所述约束条件时,所述更新后的车辆作为所述最终的车辆,并将所述最终的车辆中的待出行行程从所述预处理行程集中删除,得到更新后的预处理行程集;When the updated vehicle satisfies the constraint condition, the updated vehicle is used as the final vehicle, and the to-be-traveled trips in the final vehicle are deleted from the pre-processed trip set to obtain an updated pre-processed trip set;

判断所述行程队列中待出行行程的数量是否为0;Determine whether the number of pending trips in the trip queue is 0;

当所述行程队列中待出行行程的数量不为0时,将所述行程队列作为待处理行程集,并返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述车辆中,得到更新后的车辆和行程队列”;When the number of pending trips in the trip queue is not 0, the trip queue is used as a pending trip set, and the process returns to the step of "selecting a preset pre-processed trip from the pending trip set and adding it to the vehicle to obtain an updated vehicle and trip queue";

当所述行程队列中待出行行程的数量为0时,返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述初始行程中,得到更新后的车辆和行程队列”;When the number of pending trips in the trip queue is 0, return to the step of "selecting a preset pre-processed trip from the pending trip set to add to the initial trip to obtain an updated vehicle and trip queue";

判断所述更新后的预处理行程集中待出行行程的数量是否为0;Determine whether the number of pending trips in the updated pre-processed trip set is 0;

当所述更新后的预处理行程集中待出行行程的数量不为0时,将所述更新后的预处理行程集作为预处理行程集,并返回步骤“将预处理行程集中出发时间最早的待出行行程作为车辆的基准行程”继续执行;When the number of to-be-traveled trips in the updated pre-processed trip set is not 0, the updated pre-processed trip set is used as the pre-processed trip set, and the process returns to the step of "using the to-be-traveled trip with the earliest departure time in the pre-processed trip set as the reference trip of the vehicle" to continue the execution;

当所述更新后的预处理行程集中待出行行程的数量为0时,确定所述最终的车辆为与各所述待出行行程对应的多个待出行行程。When the number of the to-be-traveled trips in the updated pre-processed trip set is 0, the final vehicle is determined to be a plurality of to-be-traveled trips corresponding to the to-be-traveled trips.

可选地,对各所述解应用邻域操作算子,得到各所述解对应的搜索结果,具体包括:Optionally, applying a neighborhood operation operator to each of the solutions to obtain search results corresponding to each of the solutions specifically includes:

对各所述解应用多个邻域操作算子进行搜索,得到各所述解对应的多个初始搜索结果;Applying a plurality of neighborhood operation operators to search for each of the solutions to obtain a plurality of initial search results corresponding to each of the solutions;

根据多个所述初始搜索结果,计算各所述解对应的提升度;Calculating the lift corresponding to each of the solutions according to the plurality of initial search results;

当所述搜索次数是否小于预设次数时,返回步骤“对所述解中预设解应用邻域操作算子,得到初始搜索结果”继续执行;When the number of searches is less than the preset number, returning to the step of "applying a neighborhood operation operator to the preset solution in the solution to obtain an initial search result" to continue execution;

当所述搜索次数大于等于预设次数时,根据各所述解和各所述解对应的所述初始搜索结果,计算各所述解的改进拥挤距离;When the number of searches is greater than or equal to a preset number of times, calculating an improved crowding distance of each solution according to each solution and the initial search result corresponding to each solution;

将各所述解的改进拥挤距离进行归一化,得到各解的解选择概率;Normalizing the improved crowding distance of each solution to obtain the solution selection probability of each solution;

根据各所述解对应的提升度,确定各邻域操作算子的操作选择概率;Determining the operation selection probability of each neighborhood operation operator according to the lifting degree corresponding to each of the solutions;

根据解选择概率,从多个所述解中选择一个待处理的解;Selecting a solution to be processed from a plurality of solutions according to a solution selection probability;

根据操作选择概率,从多个所述邻域操作算子中确定目标邻域操作;Determining a target neighborhood operation from a plurality of neighborhood operation operators according to the operation selection probability;

应用所述目标邻域操作,对所述待处理的解进行局部搜索,得到所述待处理的解对应的最终搜索结果;Applying the target neighborhood operation to perform a local search on the solution to be processed, and obtaining a final search result corresponding to the solution to be processed;

根据所述最终搜索结果对所述初始搜索结果进行更新,得到各所述解对应的搜索结果。The initial search results are updated according to the final search results to obtain search results corresponding to each of the solutions.

一种城际出行车辆路径确定系统,应用上述的城际出行车辆路径确定方法,所述确定系统包括:A system for determining a path for an inter-city travel vehicle, applying the above-mentioned method for determining a path for an inter-city travel vehicle, the determination system comprising:

构建模块,用于以最大化平均每趟出行的乘客数、最小化车辆数目、最小化车辆总行驶距离和最小化乘客总等待时间为目标函数,以车辆载客量、服务质量、时间约束和安全约束为约束条件,构建城际出行车辆路径模型;A construction module is used to construct an inter-city travel vehicle path model with the objective function of maximizing the average number of passengers per trip, minimizing the number of vehicles, minimizing the total vehicle travel distance and minimizing the total passenger waiting time, and with the vehicle passenger capacity, service quality, time constraint and safety constraint as constraints;

获取模块,用于获取待出行城际订单;所述待出行城际订单包括多个出行订单;各所述出行订单包括乘车人数、出发时间、出发地点和目的地点;An acquisition module is used to acquire inter-city travel orders to be traveled; the inter-city travel orders to be traveled include multiple travel orders; each of the travel orders includes the number of passengers, departure time, departure location and destination point;

方案确定模块,用于根据所述城际出行车辆路径模型,确定各所述出行订单对应的出行车辆;其中,一个所述出行车辆和由所述出行车辆完成的多个所述出行订单构成一个出行解;多个所述出行解构成出行外部存档;将所述出行外部存档作为所述待出行城际订单的路径规划方案。A solution determination module is used to determine the travel vehicles corresponding to each travel order according to the inter-city travel vehicle path model; wherein, one travel vehicle and multiple travel orders completed by the travel vehicle constitute a travel solution; multiple travel solutions constitute a travel external archive; and the travel external archive is used as the path planning solution for the inter-city order to be traveled.

一种电子设备,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行所述计算机程序以使所述电子设备执行上述的城际出行车辆路径确定方法。An electronic device includes a memory and a processor, wherein the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the above-mentioned method for determining a vehicle path for inter-city travel.

一种计算机可读存储介质,其存储有计算机程序,所述计算机程序被处理器执行时实现上述的城际出行车辆路径确定方法。A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned method for determining a vehicle path for inter-city travel.

根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

本发明将城际网约车路径优化问题定义为一个包含四个目标的多目标问题,以最大化平均每趟出行的乘客数、最小化车辆数目、最小化车辆总行驶距离和最小化乘客总等待时间为目标函数,以车辆载客量、服务质量、时间约束和安全约束为约束条件,构建城际出行车辆路径模型;并根据获取待出行城际订单对城际出行车辆路径模型进行求解,确定各出行订单对应的出行车辆;本发明能够更加全面真实地反映城际出行问题的本质,为城际网约车出行提供满足多个需求的路径规划方案。The present invention defines the inter-city online car-hailing route optimization problem as a multi-objective problem with four objectives, with maximizing the average number of passengers per trip, minimizing the number of vehicles, minimizing the total vehicle driving distance and minimizing the total waiting time of passengers as the objective function, and with vehicle passenger capacity, service quality, time constraints and safety constraints as constraints, to construct an inter-city travel vehicle route model; and solve the inter-city travel vehicle route model based on the inter-city orders to be traveled, and determine the travel vehicles corresponding to each travel order; the present invention can more comprehensively and truly reflect the essence of the inter-city travel problem, and provide a route planning solution that meets multiple needs for inter-city online car-hailing travel.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

图1为本发明城际出行车辆路径确定方法流程图;FIG1 is a flow chart of a method for determining a vehicle path for inter-city travel according to the present invention;

图2为本发明求解多目标城际出行车辆路径问题的自适应局部搜索方法的流程图。FIG2 is a flow chart of an adaptive local search method for solving a multi-objective intercity travel vehicle routing problem according to the present invention.

图3为本发明基于时间序列的初始化构造方法的流程图。FIG3 is a flow chart of the initialization construction method based on time series of the present invention.

具体实施方式DETAILED DESCRIPTION

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

本发明的目的是提供一种城际出行车辆路径确定方法、系统、电子设备及介质,通过能够为城际网约车出行提供满足多个需求的路径规划方案。The purpose of the present invention is to provide a method, system, electronic device and medium for determining a vehicle path for inter-city travel, which can provide a path planning solution that meets multiple needs for inter-city online car-hailing travel.

本发明在选定优化目标以及驾驶安全,疲劳驾驶限制等方面均有特殊性,因此拟分析城际出行在实际应用中的各种约束,结合城际出行的特性,构建具有一定通用性的路径优化模型。致力于针对城际出行车辆路径问题的多目标特性和现实应用,建立一个考虑全面、贴合实际的通用多目标路径优化模型,并基于真实世界数据,构建城际出行案例的测试算例。并设计针对各目标的邻域搜索操作,最后结合基于学习的自适应选择框架来解决城际出行车辆路径问题,从而最终高效智能地规划城际网约车的出行路径,给出行客户带来舒适的出行体验,还能节约运输成本,给公司和司机带来更多收益。在城际网约车行业中,对于出行路径方案的规划是网约车服务中的关键问题之一。The present invention has particularities in terms of selecting optimization objectives, driving safety, and fatigue driving restrictions. Therefore, it is intended to analyze the various constraints of intercity travel in practical applications, and to build a path optimization model with certain universality in combination with the characteristics of intercity travel. It is committed to establishing a comprehensive and practical universal multi-objective path optimization model for the multi-objective characteristics and practical applications of intercity travel vehicle path problems, and to construct a test case for intercity travel cases based on real-world data. And design neighborhood search operations for each goal, and finally combine a learning-based adaptive selection framework to solve the intercity travel vehicle path problem, so as to ultimately plan the travel path of intercity online car-hailing efficiently and intelligently, bringing a comfortable travel experience to travelers, saving transportation costs, and bringing more benefits to companies and drivers. In the intercity online car-hailing industry, the planning of travel path plans is one of the key issues in online car-hailing services.

在本发明中,旨在完成将订单分配至行程,将行程分配至车辆(出行任务)的两阶段分配。通过优化车辆路线来最小化车辆总行驶距离、最小化乘客总等待时间、最大化出行上座率和最小化车辆数目。因此,提出了一种基于时间序列的初始化构造方法来生成初始化规划方案。在该方法中,基于时间顺序的随机选择,因此可以保证得到规划方案的同时而不让规划方案过早收敛,在达到解质量相对较优的同时具有多样性。同时注意,由于本发明研究模型为城际出行场景,故存在两种类型的订单集(即从城市a到城市b的订单集合Ord1和从城市b到城市a的订单集合Ord2,在此基础上得到的行程集Trips1/Trips2和车辆集Car1/Car2也为两种类型)。In the present invention, it is intended to complete the two-stage allocation of allocating orders to trips and allocating trips to vehicles (travel tasks). By optimizing the vehicle route, the total vehicle driving distance, the total passenger waiting time, the maximum travel occupancy rate and the number of vehicles are minimized. Therefore, an initialization construction method based on time series is proposed to generate an initialization planning scheme. In this method, based on the random selection of time sequence, it is possible to ensure that the planning scheme is obtained without allowing the planning scheme to converge prematurely, and it has diversity while achieving relatively good solution quality. At the same time, it should be noted that since the research model of the present invention is an inter-city travel scenario, there are two types of order sets (i.e., the order set Ord 1 from city a to city b and the order set Ord 2 from city b to city a, and the trip set Trips 1 /Trips 2 and vehicle set Car 1 /Car 2 obtained on this basis are also two types).

为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

实施例一Embodiment 1

如图1所示,本发明提供了一种城际出行车辆路径确定方法,所述确定方法包括:As shown in FIG1 , the present invention provides a method for determining a path of an inter-city travel vehicle, the method comprising:

步骤S1:以最大化平均每趟出行的乘客数、最小化车辆数目、最小化车辆总行驶距离和最小化乘客总等待时间为目标函数,以车辆载客量、服务质量、时间约束和安全约束为约束条件,构建城际出行车辆路径模型。Step S1: Taking maximizing the average number of passengers per trip, minimizing the number of vehicles, minimizing the total vehicle travel distance and minimizing the total passenger waiting time as the objective function, and taking vehicle passenger capacity, service quality, time constraints and safety constraints as constraints, an inter-city travel vehicle path model is constructed.

具体地,所述目标函数为:Specifically, the objective function is:

其中,目标f1表示最大化平均每趟行程的乘客数,Mt、Mk分别表示行程t、车辆k的数量,Qt表示行程t的上座率,目标f2表示最小化车辆数目,目标f3表示最小化车辆总行驶距离,Mt表示行程t的数量,lt表示行程t的总行驶距离,用来判断行程t是否安排至车辆k,目标f4表示最小化乘客总等待时间。Among them, the goal f1 represents maximizing the average number of passengers per trip, Mt and Mk represent the number of trips t and vehicles k respectively, Qt represents the occupancy rate of trip t, the goal f2 represents minimizing the number of vehicles, the goal f3 represents minimizing the total driving distance of vehicles, Mt represents the number of trips t, lt represents the total driving distance of trip t, It is used to determine whether trip t is arranged to vehicle k. Objective f 4 represents minimizing the total waiting time of passengers.

城际网约车的路径规划问题可定义为一家私人公司已经开发了一个平台,为两个城市之间的乘客提供城际乘车共享(IRS)服务。两个城市之间的乘客需要提前几个小时向平台提供他们的旅行信息。具体描述为:一组订单(记作C={1,2,…,n})由网约车车辆座位规模相同的公司提供服务,每个订单都有一个乘车人数即需求qi和一个出发时间窗口[bi,ei],并只能由一辆车服务,且每辆车都属于车辆起点和目的地所在的城市。该车根据每个订单的上下车点规划接送路线,并依次将每个订单的乘客送至对应下车点。且要求任何车辆k∈K在接下一次的新乘客之前,需要让每次旅行的所有乘客下车。出于城际出行特点的考虑,本发明乘客订单的时间信息为软时间窗,即晚于最晚出发时间的差值定义为乘客的等待时间。同时为了确保驾驶安全,每个司机必须在国家规定的连续行驶时间4小时内进行最少20分钟的休息。同时,每个订单的乘客数不能超过该车辆最大载客量且同一条行程分配的所有订单的总乘客数不能超过该车辆的最大载客数。The path planning problem of intercity online ride-hailing can be defined as a private company has developed a platform to provide intercity ride-sharing (IRS) services for passengers between two cities. Passengers between the two cities need to provide their travel information to the platform a few hours in advance. The specific description is: a group of orders (denoted as C = {1, 2, ..., n}) are served by companies with the same seat size of online ride-hailing vehicles. Each order has a number of passengers, i.e., demand q i and a departure time window [b i , e i ], and can only be served by one vehicle, and each vehicle belongs to the city where the vehicle's starting point and destination are located. The vehicle plans the pick-up route according to the boarding and alighting points of each order, and sends the passengers of each order to the corresponding boarding and alighting points in turn. And it is required that any vehicle k∈K needs to let all passengers of each trip get off before picking up the next new passenger. Considering the characteristics of intercity travel, the time information of the passenger order of the present invention is a soft time window, that is, the difference later than the latest departure time is defined as the waiting time of the passenger. At the same time, in order to ensure driving safety, each driver must take a minimum of 20 minutes of rest within the 4 hours of continuous driving time stipulated by the state. At the same time, the number of passengers for each order cannot exceed the maximum passenger capacity of the vehicle, and the total number of passengers for all orders assigned to the same trip cannot exceed the maximum passenger capacity of the vehicle.

在城际网约车路径规划问题的模型中,包含多个约束条件,具体定义如下:The model of the inter-city online car-hailing route planning problem contains multiple constraints, which are defined as follows:

1)容量约束:所有车辆服务的乘客数量不超过最大载客数,即每个订单的乘客数不能超过该车辆最大载客数且同一条行程分配的所有订单的总乘客数不能超过该车辆的最大载客数。即满足:1) Capacity constraint: The number of passengers served by all vehicles does not exceed the maximum passenger capacity, that is, the number of passengers in each order cannot exceed the maximum passenger capacity of the vehicle and the total number of passengers in all orders assigned to the same trip cannot exceed the maximum passenger capacity of the vehicle. That is, it satisfies:

其中,表示行程t的第p个订单的乘客数量,Mt表示该车辆行程t的数量。in, represents the number of passengers in the pth order of trip t, and M t represents the number of trips t of the vehicle.

2)服务约束一:为了保证网约车的服务质量,任何车辆k∈K在接下一批的新乘客之前,需要让这次路径的所有乘客下车。即满足:2) Service Constraint 1: In order to ensure the service quality of online ride-hailing, any vehicle k∈K needs to get off all passengers on this route before picking up the next batch of new passengers. That is, it satisfies:

其中,表示车辆k服务完接送j后的当前乘客人数,是0-1变量,用于判断接送点i是否被安排至车辆k,是0-1变量,用于判断行程t是否安排至车辆k,是0-1变量,判断出发点i是否为行程t的最后一个接送点。可知该式成立当且仅当即此时车辆k服务的接送j是该行程的最后一个接送点,故此时表示车辆k服务完最后一个接送点j后乘客人数为0。in, represents the current number of passengers after vehicle k has completed its service to pick up j, is a 0-1 variable used to determine whether the pick-up point i is assigned to vehicle k. is a 0-1 variable used to determine whether trip t is scheduled to vehicle k. is a 0-1 variable that determines whether the starting point i is the last pick-up point of trip t. It can be seen that this formula is valid if and only if That is, the pick-up point j served by vehicle k is the last pick-up point of the trip, so at this time It means that the number of passengers after vehicle k serves the last pick-up point j is 0.

3)服务约束二:每个订单由一辆车提供服务。3) Service constraint 2: Each order is served by one vehicle.

其中,是0-1变量,用于判断订单p是否被安排至车辆k。in, is a 0-1 variable used to determine whether order p is assigned to vehicle k.

4)服务约束三:在城际出行中,任意一个订单的上下车接送点仅会出现在一条且同一条路径中,并且一个订单在出发城市的上车点必须出现在该订单在目的城市的下车点之前。4) Service Constraint Three: In inter-city travel, the pick-up and drop-off points of any order will only appear on one and the same route, and the pick-up point of an order in the departure city must appear before the drop-off point of the order in the destination city.

其中,是0-1变量,用于判断接送点i,j是否被安排至订单p。是0-1变量,用于判断订单p是否被安排至车辆k。表示车辆k到达接送点i的时间。in, It is a 0-1 variable used to determine whether the pick-up point i, j is assigned to order p. is a 0-1 variable used to determine whether order p is assigned to vehicle k. represents the time when vehicle k arrives at pick-up point i.

5)服务约束四:所有车辆一次行程里服务一个地点的次数不超过一次。5) Service constraint 4: All vehicles shall serve a location no more than once in one trip.

其中,是0-1变量,用于判断接送点i是否被安排至行程t。in, is a 0-1 variable used to determine whether pick-up point i is scheduled for trip t.

6)服务约束五:表示所有车辆若服务于该乘客的出发点,那么必须服务于该乘客的目的地。6) Service constraint five: means that all vehicles that serve the passenger’s departure point must also serve the passenger’s destination.

其中,是0-1变量,用于判断接送点i,j是否被安排至订单p。是0-1变量,用于判断订单p是否被安排至车辆k,是0-1变量,用于判断接送点i是否被安排至车辆k。in, It is a 0-1 variable used to determine whether the pick-up point i, j is assigned to order p. is a 0-1 variable used to determine whether order p is assigned to vehicle k. It is a 0-1 variable used to determine whether the pick-up point i is assigned to vehicle k.

7)时间约束:为了保证城际网约车的服务质量,开始服务时间必须满足服务时间窗口。即满足:7) Time constraints: In order to ensure the service quality of inter-city online ride-hailing services, the start time of service must meet the service time window. That is, it must meet:

其中,是0-1变量,用于判断订单p是否被安排至车辆k。eP表示订单p的最迟开始服务时间,表示车辆k到达接送点i的时间。bP表示订单p的最早开始服务时间。in, is a 0-1 variable used to determine whether order p is assigned to vehicle k. e P represents the latest service start time of order p. represents the time when vehicle k arrives at pick-up point i. b P represents the earliest service start time of order p.

8)安全约束:考虑到城际出行一次行程的行驶时间相比城市内出行更久,为了保证城际网约车出行时乘客和驾驶员的安全以及保证驾驶员的出行效率,因此要求所有车辆在达到国家规定驾驶员最长连续行驶时间后必须进行国家规定的最少休息时间。即满足:8) Safety constraints: Considering that the driving time of an inter-city trip is longer than that of an intra-city trip, in order to ensure the safety of passengers and drivers during inter-city online car-hailing trips and to ensure the travel efficiency of drivers, all vehicles are required to take the minimum rest time prescribed by the state after reaching the maximum continuous driving time prescribed by the state. That is, it meets:

其中,表示车辆k到达接送点i的时间,ti,j表示车辆从接送i到接送j的行驶时间,gi表示服务时间,TR表示国家规定超过连续行驶时间的最小休息时间,是0-1决策变量,用于判断若车辆k到达接送点i时的连续行驶时间是否超过规定的最大连续驾驶时间,在中国,规定的最大连续驾驶时间为4小时。in, represents the time when vehicle k arrives at pick-up point i, ti ,j represents the driving time of the vehicle from pick-up point i to pick-up point j, gi represents the service time, and TR represents the minimum rest time required by the state for continuous driving time. is a 0-1 decision variable used to determine whether the continuous driving time of vehicle k when it arrives at pick-up point i exceeds the prescribed maximum continuous driving time. In China, the prescribed maximum continuous driving time is 4 hours.

步骤S2:获取待出行城际订单;所述待出行城际订单包括多个出行订单;各所述出行订单包括乘车人数、出发时间、出发地点和目的地点。Step S2: Obtain intercity travel orders to be traveled; the intercity travel orders to be traveled include multiple travel orders; each of the travel orders includes the number of passengers, departure time, departure location and destination point.

步骤S3:根据所述城际出行车辆路径模型,确定各所述出行订单对应的出行车辆;其中,一个所述出行车辆和由所述出行车辆完成的多个所述出行订单构成一个出行解;多个所述出行解构成出行外部存档;将所述出行外部存档作为所述待出行城际订单的路径规划方案。Step S3: Determine the travel vehicles corresponding to each travel order according to the inter-city travel vehicle path model; wherein, one travel vehicle and multiple travel orders completed by the travel vehicle constitute a travel solution; multiple travel solutions constitute a travel external archive; and use the travel external archive as the path planning solution for the inter-city order to be traveled.

S3具体包括:S3 specifically includes:

步骤S31:根据所述约束条件和所述目标函数,确定与各所述出行订单对应的多个待出行行程;其中,各所述待出行行程包括至少一个所述出行订单且各所述待出行行程满足所述约束条件和所述目标函数。Step S31: determining a plurality of itineraries to be traveled corresponding to each of the travel orders according to the constraints and the objective function; wherein each of the itineraries to be traveled includes at least one of the travel orders and each of the itineraries to be traveled satisfies the constraints and the objective function.

S31具体包括:S31 specifically includes:

步骤S3100:将预处理订单集中出发时间最早的出行订单作为初始行程的基准订单;所述初始行程至少包括一个所述出行订单;所述预处理订单集包括出发地点相同的多个所述出行订单。Step S3100: taking the travel order with the earliest departure time in the pre-processed order set as the benchmark order for the initial trip; the initial trip includes at least one travel order; the pre-processed order set includes multiple travel orders with the same departure point.

步骤S3101:从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列;所述待处理订单集为所述预处理订单集去掉所述基准订单后的剩余订单;所述订单队列为从所述待处理订单集中去掉一个所述预设预处理订单后的剩余订单。Step S3101: Select a preset preprocessed order from the pending order set and add it to the initial itinerary to obtain an updated initial itinerary and order queue; the pending order set is the remaining orders after removing the benchmark order from the preprocessed order set; the order queue is the remaining orders after removing one of the preset preprocessed orders from the pending order set.

步骤S3102:判断所述更新后的初始行程是否满足所述约束条件。Step S3102: Determine whether the updated initial itinerary satisfies the constraint condition.

步骤S3103:当所述更新后的初始行程不满足所述约束条件时,将所述预设预处理订单从所述更新后的初始行程中删除,并返回步骤“从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列”继续执行。Step S3103: When the updated initial itinerary does not satisfy the constraint condition, the preset pre-processed order is deleted from the updated initial itinerary, and the process returns to the step of "selecting a preset pre-processed order from the set of pending orders and adding it to the initial itinerary to obtain an updated initial itinerary and order queue" to continue execution.

步骤S3104:当所述更新后的初始行程满足所述约束条件时,判断所述更新后的初始行程的乘车人数是否达到所述车辆载客量。Step S3104: When the updated initial trip satisfies the constraint condition, it is determined whether the number of passengers in the updated initial trip reaches the vehicle passenger capacity.

步骤S3105:当所述更新后的初始行程的乘车人数达到所述车辆载客量时,所述更新后的初始行程作为所述最终的初始行程,并将所述最终的初始行程中的出行订单从所述预处理订单集中删除,得到更新后的预处理订单集。Step S3105: When the number of passengers in the updated initial trip reaches the passenger capacity of the vehicle, the updated initial trip is used as the final initial trip, and the travel orders in the final initial trip are deleted from the pre-processed order set to obtain an updated pre-processed order set.

步骤S3106:当所述更新后的初始行程的乘车人数未达到所述车辆载客量时,判断所述订单队列中出行订单的数量是否为0。Step S3106: When the number of passengers in the updated initial trip does not reach the passenger capacity of the vehicle, determine whether the number of travel orders in the order queue is 0.

步骤S3107:当所述订单队列中出行订单的数量不为0时,将所述订单队列作为待处理订单集,并返回步骤“从待处理订单集中选择一个预设预处理订单加入到所述初始行程中,得到更新后的初始行程和订单队列”。Step S3107: When the number of travel orders in the order queue is not 0, the order queue is used as a pending order set, and the process returns to the step of "selecting a preset pre-processed order from the pending order set and adding it to the initial trip to obtain an updated initial trip and order queue".

步骤S3108:当所述订单队列中出行订单的数量为0时,返回步骤“所述更新后的初始行程作为所述最终的初始行程,并将所述最终的初始行程中的出行订单从所述预处理订单集中删除,得到更新后的预处理订单集”。Step S3108: When the number of travel orders in the order queue is 0, return to the step of "taking the updated initial itinerary as the final initial itinerary, and deleting the travel orders in the final initial itinerary from the pre-processed order set to obtain an updated pre-processed order set".

步骤S3109:当所述更新后的预处理订单集中出行订单的数量不为0时,将所述更新后的预处理订单集作为预处理订单集,并返回步骤“将预处理订单集中出发时间最早的出行订单作为初始行程的基准订单”继续执行。Step S3109: When the number of travel orders in the updated pre-processed order set is not 0, the updated pre-processed order set is used as the pre-processed order set, and the process returns to the step of "using the travel order with the earliest departure time in the pre-processed order set as the benchmark order for the initial trip" to continue execution.

步骤S3110:当所述更新后的预处理订单集中出行订单的数量为0时,确定所述最终的初始行程为与各所述出行订单对应的多个待出行行程。Step S3110: When the number of travel orders in the updated pre-processed order set is 0, determining that the final initial itinerary is a plurality of to-be-traveled itineraries corresponding to each of the travel orders.

在实际应用中,如图3所示,确定与各所述出行订单对应的多个待出行行程是基于时间序列的初始化构造方法,包括以下七个步骤:In practical applications, as shown in FIG3 , determining a plurality of to-be-traveled trips corresponding to each of the travel orders is based on an initialization construction method of a time series, and includes the following seven steps:

步骤一,对两种订单集Ord1/Ord2按订单最早出发时间升序排序,选择待处理订单集Ord1/Ord2的第一个订单o1为基准生成一个行程节点向量t。将订单o1从待处理订单集Ord1/Ord2移除。将待处理订单集Ord1/Ord2中的订单添加到订单队列S,并将所有订单的状态设为“未处理”。Step 1: Sort the two order sets Ord 1 /Ord 2 in ascending order according to the earliest departure time of the order, and select the first order o 1 in the pending order set Ord 1 /Ord 2 as the benchmark to generate a trip node vector t. Remove order o 1 from the pending order set Ord 1 /Ord 2. Add the orders in the pending order set Ord 1 /Ord 2 to the order queue S, and set the status of all orders to "unprocessed".

在实际应用中,一个路径规划方案X是有k个路线所组成的集合O={o1,..oi...,ok}来表示的,其中是一条由包含Ni个订单,2Ni个接送点的访问序列构成的路径,表示第i条路径的第j个客户点。由于每一个订单包含该客户在出发城市的上车点和目标城市的下车点,因此,在每一条路径中,每一个订单表示成两个接送点,即该订单在出发城市的上车点和目标城市的下车点。在一个路径规划方案中,任意一个订单的两个接送点仅会出现在一条且同一条路径中。In practical applications, a path planning scheme X is represented by a set of k routes O = {o1,..oi...,ok}, where is a path consisting of a visit sequence containing Ni orders and 2Ni pick-up and drop-off points. represents the jth customer point of the i-th path. Since each order contains the customer's pick-up point in the departure city and the drop-off point in the target city, each order is represented as two pick-up points in each path, namely the pick-up point in the departure city and the drop-off point in the target city. In a path planning scheme, the two pick-up points of any order will only appear in one and the same path.

步骤二,从当前待处理订单集Ord1/Ord2随机选择订单o加入t,将订单o状态设为“已处理”。同时检查订单o加入后的行程t是否满足模型约束,若满足,则进入步骤三,若不满足,则将订单o从行程t中移除,返回步骤二。Step 2: Randomly select order o from the current pending order set Ord 1 /Ord 2 and add it to t, and set the status of order o to "processed". At the same time, check whether the trip t after order o is added meets the model constraints. If so, proceed to step 3. If not, remove order o from trip t and return to step 2.

步骤三,检查当前添加到行程t的订单乘客数量之和是否达到车辆载客量,若满足,进入步骤四,若不满足,则进入步骤六。Step 3: Check whether the sum of the number of passengers currently added to the trip t reaches the vehicle capacity. If so, proceed to step 4; if not, proceed to step 6.

步骤四,将t添加到行程集Trips1/Trips2中,并将行程t的订单从待处理订单集Ord1/Ord2移除,同时清空t。Step 4: Add t to the itinerary set Trips 1 /Trips 2 , remove the order of itinerary t from the pending order set Ord 1 /Ord 2 , and clear t at the same time.

步骤五,判断待处理订单集Ord1/Ord2是否为空,若满足,则进入步骤七,若不满足,则返回步骤一。Step 5: Determine whether the pending order set Ord 1 /Ord 2 is empty. If so, proceed to step 7. If not, return to step 1.

步骤六,判断订单队列S中是否存在订单状态为“未处理”,若满足,则返回步骤二,若不满足,则返回步骤四。Step 6, determine whether there is an order in the order queue S with an order status of "unprocessed". If so, return to step 2; if not, return to step 4.

步骤七,初始化行程集Trips1/Trips2成功,Step 7: Initialize the trip set Trips 1 /Trips 2 successfully.

步骤S32:根据所述约束条件和所述目标函数,确定与各所述待出行行程对应的车辆;其中,一个所述车辆和由所述车辆完成的多个所述待出行行程构成一个解;多个所述解构成外部存档。Step S32: Determine the vehicle corresponding to each of the to-be-traveled trips according to the constraint conditions and the objective function; wherein one of the vehicles and a plurality of the to-be-traveled trips completed by the vehicle constitute a solution; and a plurality of the solutions constitute an external archive.

S32具体包括:S32 specifically includes:

步骤S3201:将预处理行程集中出发时间最早的待出行行程作为车辆的基准行程;所述车辆至少包括一个所述待出行行程;所述预处理行程集包括出发地点相同的多个所述待出行行程。Step S3201: taking the trip to be traveled with the earliest departure time in the pre-processed trip set as the reference trip of the vehicle; the vehicle includes at least one trip to be traveled; the pre-processed trip set includes multiple trips to be traveled with the same departure point.

步骤S3202:从待处理行程集中选择一个预设预处理行程加入到所述初始行程中,得到更新后的车辆和行程队列;所述待处理行程集为所述预处理行程集去掉所述基准行程后的剩余行程;所述行程队列为从所述待处理行程集中去掉一个所述预设预处理行程后的剩余行程。Step S3202: Select a preset pre-processed trip from the set of trips to be processed and add it to the initial trip to obtain an updated vehicle and trip queue; the set of trips to be processed is the remaining trips after removing the reference trip from the pre-processed trip set; the trip queue is the remaining trips after removing one of the preset pre-processed trips from the set of trips to be processed.

步骤S3203:判断所述更新后的车辆是否满足所述约束条件。Step S3203: Determine whether the updated vehicle satisfies the constraint condition.

步骤S3204:当所述更新后的车辆不满足所述约束条件时,将所述预设预处理行程从所述更新后的车辆中删除,并返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述车辆中,得到更新后的车辆和行程队列”继续执行。Step S3204: When the updated vehicle does not satisfy the constraint condition, the preset pre-processed trip is deleted from the updated vehicle, and the process returns to the step of "selecting a preset pre-processed trip from the set of trips to be processed and adding it to the vehicle to obtain an updated vehicle and trip queue" to continue execution.

步骤S3205:当所述更新后的车辆满足所述约束条件时,所述更新后的车辆作为所述最终的车辆,并将所述最终的车辆中的待出行行程从所述预处理行程集中删除,得到更新后的预处理行程集。Step S3205: When the updated vehicle satisfies the constraint condition, the updated vehicle is used as the final vehicle, and the to-be-traveled trips in the final vehicle are deleted from the pre-processed trip set to obtain an updated pre-processed trip set.

步骤S3206:判断所述行程队列中待出行行程的数量是否为0。Step S3206: Determine whether the number of pending trips in the trip queue is 0.

步骤S3207:当所述行程队列中待出行行程的数量不为0时,将所述行程队列作为待处理行程集,并返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述车辆中,得到更新后的车辆和行程队列”。Step S3207: When the number of pending trips in the trip queue is not 0, the trip queue is used as a pending trip set, and the process returns to the step of "selecting a preset pre-processed trip from the pending trip set and adding it to the vehicle to obtain an updated vehicle and trip queue".

步骤S3208:当所述行程队列中待出行行程的数量为0时,返回步骤“从待处理行程集中选择一个预设预处理行程加入到所述初始行程中,得到更新后的车辆和行程队列”。Step S3208: When the number of pending trips in the trip queue is 0, return to the step of "selecting a preset pre-processed trip from the pending trip set and adding it to the initial trip to obtain an updated vehicle and trip queue".

步骤S3209:判断所述更新后的预处理行程集中待出行行程的数量是否为0。Step S3209: Determine whether the number of pending trips in the updated pre-processed trip set is 0.

步骤S3210:当所述更新后的预处理行程集中待出行行程的数量不为0时,将所述更新后的预处理行程集作为预处理行程集,并返回步骤“将预处理行程集中出发时间最早的待出行行程作为车辆的基准行程”继续执行。Step S3210: When the number of to-be-traveled trips in the updated pre-processed trip set is not 0, the updated pre-processed trip set is used as the pre-processed trip set, and the process returns to the step of "using the to-be-traveled trip with the earliest departure time in the pre-processed trip set as the vehicle's reference trip" to continue execution.

步骤S3211:当所述更新后的预处理行程集中待出行行程的数量为0时,确定所述最终的车辆为与各所述待出行行程对应的多个待出行行程。Step S3211: when the number of the to-be-traveled trips in the updated pre-processed trip set is 0, determining that the final vehicle is a plurality of to-be-traveled trips corresponding to each of the to-be-traveled trips.

在实际应用中,如图3所示,确定与各所述待出行行程对应的车辆的过程具体包括:In practical applications, as shown in FIG3 , the process of determining the vehicle corresponding to each of the to-be-traveled trips specifically includes:

步骤一,并对两种行程集Trips1/Trips2按行程最早开始时间升序排序。Step 1: Sort the two trip sets Trips 1 /Trips 2 in ascending order according to the earliest start time of the trip.

步骤二,选择待处理行程集Trips1/Trips2的第一个行程t1为基准生成一个车辆节点向量c。将行程t1从待处理行程集Trips1/Trips2移除。将待处理行程集Trips1/Trips2中的行程添加到行程队列T,并将所有行程的状态设为“未处理”。Step 2: Select the first trip t 1 of the pending trip set Trips 1 /Trips 2 as a reference to generate a vehicle node vector c. Remove trip t 1 from the pending trip set Trips 1 /Trips 2. Add the trips in the pending trip set Trips 1 /Trips 2 to the trip queue T, and set the status of all trips to "unprocessed".

步骤三,从当前待处理行程集Trips1/Trips2随机选择行程t加入c,将行程t状态设为“已处理”,同时检查行程t加入后的车辆c是否满足模型约束,若满足,则进入步骤四,若不满足,则将该行程t从车辆c中移除,返回步骤三。Step 3: Randomly select trip t from the current set of pending trips Trips 1 /Trips 2 and add it to c, set the status of trip t to "processed", and check whether the vehicle c after adding trip t meets the model constraints. If yes, go to step 4; if not, remove trip t from vehicle c and return to step 3.

步骤四,判断行程队列T中是否存在行程状态为“未处理”,若满足,则返回步骤三,若不满足,则进入步骤五。Step 4: Determine whether there is a trip in the trip queue T with a status of "unprocessed". If so, return to step 3; if not, proceed to step 5.

步骤五,将c添加到车辆集Car1/Car2中,并将车辆c的行程从待处理行程集Trips1/Trips2移除,同时清空c。判断待处理行程集Trips1/Trips2是否为空,若满足,则进入步骤六,若不满足,则返回步骤二。Step 5: add c to the vehicle set Car 1 /Car 2 , remove the trip of vehicle c from the pending trip set Trips 1 /Trips 2 , and clear c. Determine whether the pending trip set Trips 1 /Trips 2 is empty. If so, proceed to step 6. If not, return to step 2.

步骤六,初始化车辆集Car1/Car2成功(两种车辆集的所有车辆即构成一种解决方案)。Step 6: Initialize the vehicle set Car 1 /Car 2 successfully (all vehicles in the two vehicle sets constitute a solution).

步骤S33:对各所述解应用邻域操作算子,得到各所述解对应的搜索结果。Step S33: applying a neighborhood operation operator to each of the solutions to obtain search results corresponding to each of the solutions.

具体地,本发明提供的邻域操作算子为八邻域操作算子,分别为:Specifically, the neighborhood operation operators provided by the present invention are eight neighborhood operation operators, which are:

邻域操作一,针对单个客户点的删除再插入操作。邻域操作二,针对单个客户点的交换操作。邻域操作三、针对多个客户点的删除再插入操作。邻域操作四,针对多个交换的删除再插入操作。邻域操作五,针对单个行程的删除再插入操作。邻域操作六,针对单个行程的交换操作。邻域操作七,针对多个行程的删除再插入操作。邻域操作八,针对多个行程的交换操作。Neighborhood operation 1: delete and then insert operation for a single customer point. Neighborhood operation 2: exchange operation for a single customer point. Neighborhood operation 3: delete and then insert operation for multiple customer points. Neighborhood operation 4: delete and then insert operation for multiple exchanges. Neighborhood operation 5: delete and then insert operation for a single trip. Neighborhood operation 6: exchange operation for a single trip. Neighborhood operation 7: delete and then insert operation for multiple trips. Neighborhood operation 8: exchange operation for multiple trips.

进一步地,本发明设计的八种针对城际出行问题特点的邻域操作算子来进行对解的局部搜索优化。鉴于面向城际出行的车辆路径问题本质上是一个多目标优化问题,而多目标优化问题是优化问题中的核心问题之一,在该问题中,不同的目标在物理意义、规模和优化难度上具有不同的特性,因此,针对一些用常规邻域搜索操作难以优化的目标,需要设计针对城际出行特性的邻域搜索操作。具体定义如下:Furthermore, the present invention designs eight neighborhood operation operators targeting the characteristics of intercity travel problems to perform local search optimization for solutions. Given that the vehicle routing problem for intercity travel is essentially a multi-objective optimization problem, and the multi-objective optimization problem is one of the core problems in optimization problems, in which different objectives have different characteristics in terms of physical meaning, scale, and optimization difficulty. Therefore, for some objectives that are difficult to optimize using conventional neighborhood search operations, it is necessary to design neighborhood search operations targeting the characteristics of intercity travel. The specific definitions are as follows:

邻域操作一:删除一个客户点重新插入最优位置,随机选择一条路径上的一个客户点进行删除,并将其客户点重新插入到最佳位置。其中路径下车点顺序应用极角排序的方法确定,具体为:根据计算下车点的极角并排序后,重新安排客户点的下车顺序(以下邻域操作后均会对客户点下车路线重新调整,采用的方法与其一致)。Neighborhood operation 1: Delete a customer point and reinsert it to the optimal position. Randomly select a customer point on a path to delete it and reinsert it to the optimal position. The order of the drop-off points on the path is determined by the polar angle sorting method. Specifically, the order of the drop-off points is rearranged after calculating the polar angles of the drop-off points and sorting them (the following neighborhood operations will readjust the drop-off routes of the customer points, and the methods used are the same).

邻域操作二:交换单个客户点,首先,随机选择一条路径上的一个客户点准备进行交换,第二步,将该客户点依次尝试和其他路径上的客户点进行交换。记录delta值,为优化后该目标的值。第三步,比较得到min_delta值,将该客户点与此时得到该目标最小值的情况下的客户点进行交换。Neighborhood operation 2: Exchange a single customer point. First, randomly select a customer point on a path to exchange. Second, try to exchange this customer point with customer points on other paths in turn. Record the delta value, which is the value of the target after optimization. Third, compare and get the min_delta value, and exchange this customer point with the customer point that gets the minimum value of the target at this time.

邻域操作三:删除多个客户点重新插入最优位置,随机选择一条路径上的多个客户点进行删除,并将这些客户点重新插入到最佳位置。Neighborhood operation three: Delete multiple customer points and reinsert them into the optimal position. Randomly select multiple customer points on a path to delete and reinsert them into the optimal position.

邻域操作四:交换多个客户点(two-opt操作),首先,随机选择一条路径上的一个客户点作为起点准备进行交换,第二步,轮流选择其他路径上的一个客户点作为第二条交换的起点,将该客户点为起点的路径尝试和交换路径上的客户点为起点的路径进行交换。记录delta值,为优化后该目标的值。第三步,比较得到min_delta值,将以该客户点为起点的路径与此时得到该目标最小值的情况下的客户点进行交换。Neighborhood operation 4: Exchange multiple customer points (two-opt operation). First, randomly select a customer point on a path as the starting point for exchange. Second, select a customer point on another path as the starting point for the second exchange in turn, and try to exchange the path starting from the customer point on the exchange path with the path starting from the customer point on the exchange path. Record the delta value, which is the value of the target after optimization. Third, compare and get the min_delta value, and exchange the path starting from the customer point with the customer point that gets the minimum value of the target at this time.

邻域操作五:删除一条路径并将路径上的客户点插入其他路径(针对单个行程删除再插入),第一步,随机选择一条路径(假设为路径trip_ind)。第二步,将路径上的每个客户点尝试插入其他路径,1:该路径上的客户点不可以全部删除插入成功,则还原路径返回false,2:若该路径上的客户点全部删除插入成功,则此时路径trip_ind为空路径。此时分为两种情况继续进行考虑:Neighborhood operation 5: Delete a path and insert the customer points on the path into other paths (delete and then insert for a single trip). The first step is to randomly select a path (assuming it is the path trip_ind). The second step is to try to insert each customer point on the path into other paths. 1: If all customer points on the path cannot be deleted and inserted successfully, the restore path returns false. 2: If all customer points on the path are deleted and inserted successfully, the path trip_ind is now an empty path. There are two cases to consider:

2-1:路径trip_ind是该车辆的最后一条路径,则可以直接进行删除路径trip_ind。返回true。2-1: If the path trip_ind is the last path of the vehicle, then the path trip_ind can be deleted directly. Return true.

2-2:路径trip_ind不是该车辆的最后一条路径,则需要检查路径trip_ind的后一条路径trip_next是否也可按照同样方法变成空路径,即将路径trip_next上的每个客户点尝试插入其他路径,若路径trip_next上的客户点不能全部删除插入成功,则还原路径trip_ind和trip_next,返回false。若路径trip_next上的客户点全部删除插入成功,则此时路径trip_next为空路径。则可以删除路径trip_ind和路径trip_next,返回true。2-2: If the path trip_ind is not the last path of the vehicle, it is necessary to check whether the path trip_next after the path trip_ind can also be turned into an empty path in the same way, that is, try to insert each customer point on the path trip_next into other paths. If the customer points on the path trip_next cannot be deleted and inserted successfully, restore the paths trip_ind and trip_next and return false. If the customer points on the path trip_next are all deleted and inserted successfully, the path trip_next is now an empty path. Then the paths trip_ind and trip_next can be deleted and true is returned.

邻域操作六:交换单个行程,首先,随机选择两辆车辆上的两条不同路径准备进行交换,第二步,将这两条路径进行交换。记录delta值,为优化后该目标的值。第三步,若delta值比原先的目标值小,则将两条路径进行交换。Neighborhood operation six: swapping a single trip. First, randomly select two different paths on two vehicles to swap. Second, swap the two paths. Record the delta value, which is the value of the target after optimization. Third, if the delta value is smaller than the original target value, swap the two paths.

邻域操作七:删除一辆车辆并将其车辆的回程插入其他车辆(针对多个行程删除再插入),第一步,选择当前路径数最少的一辆车辆car_ind。第二步,将车辆上的路径以回程为单位或特殊的单个行程(该行程没有下一条行程)尝试遍历插入其他车辆,(1)若该车辆上的路径不可以全部删除插入成功,则还原车辆返回false,(2)若该车辆上的路径全部删除插入成功,则此时车辆car_ind为空路径。则可以删除车辆car_ind,返回true。Neighborhood operation seven: Delete a vehicle and insert its return trip into other vehicles (delete and then insert for multiple trips). The first step is to select a vehicle car_ind with the least number of current paths. The second step is to try to traverse and insert other vehicles into the paths on the vehicle in units of return trips or special single trips (this trip has no next trip). (1) If the paths on the vehicle cannot be deleted and inserted successfully, the restoration vehicle returns false. (2) If the paths on the vehicle are deleted and inserted successfully, the vehicle car_ind is an empty path. Then the vehicle car_ind can be deleted and true is returned.

邻域操作八:交换单个回程(交换多个行程),首先,随机选择两辆车辆上的两条不同回程准备进行交换,第二步,将这两条回程进行交换。记录delta值,为优化后该目标的值。第三步,若delta值比原先的目标值小,则将两条回程进行交换。Neighborhood operation eight: exchange a single return trip (exchange multiple trips). First, randomly select two different return trips on two vehicles to exchange. Second, exchange the two return trips. Record the delta value, which is the value of the target after optimization. Third, if the delta value is smaller than the original target value, exchange the two return trips.

在本发明中,关于局部搜索的应用,采用以下策略:概率+概率策略:每次根据不同解的潜力值,概率从外部存档解集EP选择一个解,并根据各邻域操作的贡献度,概率地选择一种邻域操作方式进行一次深度局部搜索。其中关于邻域操作的贡献度,更应该关系它对收敛性的贡献,即对于各目标值的提升情况。所以,每个邻域操作的选择概率可以定义为L窗口(近L次搜索)内对于外部存档在收敛性上的贡献,其中,L定义为各邻域操作的各自的被选择轮数,同时设置一个矩阵来维护概率列表。关于外部存档解的潜力值,由于解集中的都是非占优集,其潜力值反应的是对于多样性的潜力,所以,每个解的选择概率可以定义为拥挤距离。即本发明采用自适应局部搜索策略用于产生新的路径规划方案。具体定义如下:In the present invention, regarding the application of local search, the following strategies are adopted: probability + probability strategy: each time according to the potential values of different solutions, a solution is selected from the external archive solution set EP with probability, and according to the contribution of each neighborhood operation, a neighborhood operation mode is selected with probability to perform a deep local search. Among them, regarding the contribution of the neighborhood operation, it should be more related to its contribution to convergence, that is, the improvement of each target value. Therefore, the selection probability of each neighborhood operation can be defined as the contribution to the convergence of the external archive within the L window (nearly L searches), where L is defined as the number of selected rounds of each neighborhood operation, and a matrix is set to maintain the probability list. Regarding the potential value of the external archive solution, since the solution set is a non-dominant set, its potential value reflects the potential for diversity, so the selection probability of each solution can be defined as the crowding distance. That is, the present invention adopts an adaptive local search strategy to generate a new path planning scheme. The specific definitions are as follows:

1.解的潜力值定义:关于外部存档解的潜力值,由于解集中的都是非占优集,其潜力值反应的是对于多样性的潜力,所以,每个解的潜力值可以定义为拥挤距离。而因为如今大多采用的计算方法为NSGA-II中拥挤距离计算方法,其具有简单快捷的特点,但是由于按照一个目标函数排序计算拥挤距离,导致最终保留的解集仅在该目标函数上是均匀分布的,而在其他函数上存在不均匀的情况。因此在其基础上,选择了改进后的拥挤距离计算方法,其计算公式如下:1. Definition of the potential value of the solution: Regarding the potential value of the external archived solution, since the solution set is a non-dominant set, its potential value reflects the potential for diversity, so the potential value of each solution can be defined as the crowding distance. Because the calculation method most commonly used today is the crowding distance calculation method in NSGA-II, it is simple and fast, but because the crowding distance is calculated by sorting according to an objective function, the final retained solution set is only uniformly distributed on the objective function, and there is an uneven situation on other functions. Therefore, on this basis, the improved crowding distance calculation method is selected, and its calculation formula is as follows:

其中:表示按第m维目标函数值从小到大排序后,第i个例子的下一个粒子的第m维目标函数值;表示按第m维目标函数值从小到大排序后,第i个粒子的前一个粒子的第m维目标函数值。改进后的拥挤距离计算方法不但依然仅考虑前后两个粒子,而且综合考略了每一维目标函数,更具合理性。此外还需注意,当i=1或i=N时,此时disi=inf,若考虑这两个极值解,在后面的计算公式中解的被选择概率可视为100%(必选),且由于实际意义中这两个极值解可视为精英解,因此直接将这两个极值解进入外部存档保存,不参与接下来的解选择。in: It represents the m-th dimension objective function value of the next particle of the ith example after sorting the m-th dimension objective function values from small to large; It represents the m-th dimension objective function value of the particle before the i-th particle after the m-th dimension objective function value is sorted from small to large. The improved crowding distance calculation method not only still only considers the two particles before and after, but also comprehensively considers each dimension of objective function, which is more reasonable. In addition, it should be noted that when i=1 or i=N, dis i =inf. If these two extreme value solutions are considered, the probability of solution selection in the subsequent calculation formula can be regarded as 100% (mandatory), and because these two extreme value solutions can be regarded as elite solutions in practical significance, these two extreme value solutions are directly saved in the external archive and do not participate in the subsequent solution selection.

从而在disi的基础上得到PVi(EP中第i个解的被选择概率):Thus, based on dis i, we can get PV i (the probability of selecting the i-th solution in EP):

其中,n表示当前EP中解的数量,ε为设置的最小被选择概率。确保各个解在整个局部搜索过程中都有机会被使用。Where n is the number of solutions in the current EP, and ε is the minimum probability of being selected, ensuring that each solution has a chance to be used during the entire local search process.

2.邻域操作的贡献度定义:关于邻域操作搜索的贡献度,更应该关系它对收敛性的贡献,所以,第k个邻域操作的选择概率可以定义为L窗口内对于外部存档在收敛性上的贡献,即对于各目标值的提升情况DOIk(Degree of improvement):2. Definition of the contribution of neighborhood operations: Regarding the contribution of neighborhood operation search, it should be more about its contribution to convergence. Therefore, the selection probability of the kth neighborhood operation can be defined as the contribution to the convergence of the external archive within the L window, that is, the improvement of each target value DOI k (Degree of improvement):

其中,Obj'i表示归一化后该邻域操作后的解的第i个目标值,Obji表示归一化后邻域操作前的解的第i个目标值。num表示该模型中的目标数量,这里应为4。同时注意该计算公式当且仅当邻域操作优化后的解进入外部存档解集EP更新才使用,其余情况DOIk为0。其中对目标值的归一化处理采用离差标准化公式,是对原始数据的线性变换,使结果落到[0,1]区间,转换函数如下:Among them, Obj'i represents the i-th target value of the solution after the neighborhood operation after normalization, and Obj i represents the i-th target value of the solution before the neighborhood operation after normalization. Num represents the number of targets in the model, which should be 4 here. At the same time, it should be noted that this calculation formula is used only when the solution after the neighborhood operation optimization enters the external archive solution set EP update, and DOI k is 0 in other cases. The normalization of the target value adopts the deviation standardization formula, which is a linear transformation of the original data so that the result falls into the [0,1] interval. The conversion function is as follows:

其中,fi表示该邻域操作前解的第i个目标值,maxi和mini分别表示EP中所有解中第i个目标中最大的值、EP中所有解中第i个目标中最小的值。Among them, fi represents the i-th objective value of the solution before the neighborhood operation, max i and min i represent the maximum value of the i-th objective among all solutions in EP and the minimum value of the i-th objective among all solutions in EP, respectively.

从而在DOIk基础上得到LSk(第k个邻域操作的被选择概率):Thus, based on DOI k, we get LS k (the probability of selection of the kth neighborhood operation):

其中,gen表示当前该邻域操作被选择次数,ε为设置的最小被选择概率。确保各个邻域操作算子在整个局部搜索过程中都有机会被使用。Among them, gen represents the number of times the neighborhood operation is currently selected, and ε is the set minimum selection probability. It ensures that each neighborhood operation operator has the opportunity to be used during the entire local search process.

S33具体包括:S33 specifically includes:

步骤S3301:对各所述解应用多个邻域操作算子进行搜索,得到各所述解对应的多个初始搜索结果。Step S3301: Apply multiple neighborhood operation operators to search for each of the solutions to obtain multiple initial search results corresponding to each of the solutions.

步骤S3302:根据多个所述初始搜索结果,计算各所述解对应的提升度。Step S3302: Calculate the lift corresponding to each of the solutions based on the multiple initial search results.

步骤S3303:当所述搜索次数是否小于预设次数时,返回步骤“对所述解中预设解应用邻域操作算子,得到初始搜索结果”继续执行。Step S3303: When the number of searches is less than the preset number, return to the step of "applying a neighborhood operation operator to the preset solution in the solution to obtain an initial search result" to continue execution.

步骤S3304:当所述搜索次数大于等于预设次数时,根据各所述解和各所述解对应的所述初始搜索结果,计算各所述解的改进拥挤距离。Step S3304: when the number of searches is greater than or equal to a preset number of searches, the improved crowding distance of each solution is calculated according to each solution and the initial search result corresponding to each solution.

步骤S3305:将各所述解的改进拥挤距离进行归一化,得到各解的解选择概率。Step S3305: normalize the improved crowding distance of each solution to obtain the solution selection probability of each solution.

步骤S3306:根据各所述解对应的提升度,确定各邻域操作算子的操作选择概率。Step S3306: Determine the operation selection probability of each neighborhood operation operator according to the lifting degree corresponding to each solution.

步骤S3307:根据解选择概率,从多个所述解中选择一个待处理的解。Step S3307: Select a solution to be processed from the multiple solutions according to the solution selection probability.

步骤S3308:根据操作选择概率,从多个所述邻域操作算子中确定目标邻域操作。Step S3308: Determine a target neighborhood operation from the plurality of neighborhood operation operators according to the operation selection probability.

步骤S3309:应用所述目标邻域操作,对所述待处理的解进行局部搜索,得到所述待处理的解对应的最终搜索结果。Step S3309: Apply the target neighborhood operation to perform a local search on the solution to be processed to obtain a final search result corresponding to the solution to be processed.

步骤S3310:根据所述最终搜索结果对所述初始搜索结果进行更新,得到各所述解对应的搜索结果。Step S3310: updating the initial search results according to the final search results to obtain search results corresponding to each of the solutions.

在实际应用中,如图2所示,得到各所述解对应的搜索结果的具体过程如下所述:In practical applications, as shown in FIG2 , the specific process of obtaining the search results corresponding to each solution is as follows:

步骤一,随机选择外部存档EP中的一个解,按平均(1/8)的概率随机选择邻域操作算子LSn(n=1...8)进行局部搜索,同时利用ε占优更新外部存档解集EP。Step 1: randomly select a solution in the external archive EP, randomly select a neighborhood operator LSn (n=1...8) with an average probability of (1/8) to perform local search, and simultaneously update the external archive solution set EP using ε dominance.

步骤二,邻域操作LSn的搜索次数mn进行加一,并计算邻域操作的提升度并将结果保存到学习矩阵中。Step 2: The number of searches mn of the neighborhood operation LSn is increased by one, and the improvement of the neighborhood operation is calculated and the result is saved in the learning matrix.

步骤三,判断是否存在邻域操作的搜索次数mn小于m次,若满足,则返回步骤一,如不满足,则进入步骤四。在实际应用中,m是预先设定的值;可选地,m的值可以根据经验值或预实验决定。Step 3: determine whether the number of searches for the neighborhood operation m n is less than m times. If so, return to step 1; if not, proceed to step 4. In practical applications, m is a preset value; optionally, the value of m can be determined based on experience or preliminary experiments.

步骤四,对外部存档解集EP的每个解进行改进拥挤距离计算,并归一化处理得到每个解的选择概率,同时基于学习矩阵按LSk的计算公式得到每个邻域操作的选择概率。Step 4: Perform improved crowding distance calculation on each solution in the external archive solution set EP, and normalize it to obtain the selection probability of each solution. At the same time, the selection probability of each neighborhood operation is obtained based on the learning matrix according to the calculation formula of LS k .

步骤五,对于外部存档解集EP,按解选择概率从其中选出一个解,再按操作概率选择一个邻域操作进行一次深度局部搜索。Step 5: For the external archive solution set EP, a solution is selected from it according to the solution selection probability, and then a neighborhood operation is selected according to the operation probability to perform a deep local search.

步骤六,利用ε占优更新外部存档解集EP并更新学习矩阵。返回步骤四。Step 6: Use ε-dominance to update the external archive solution set EP and update the learning matrix. Return to step 4.

具体地,所述的深度局部搜索操作过程,包含如下步骤:Specifically, the deep local search operation process includes the following steps:

步骤一,设置已搜索次数m=0,最大搜索深度为D。其中,D为预设的值;可选地,D的值可以根据经验值或预实验决定。Step 1: Set the number of searches m=0, and the maximum search depth to D. Wherein D is a preset value; optionally, the value of D can be determined based on experience or preliminary experiments.

步骤二,对所选择的解进行邻域操作,判断该邻域搜索操作后的解是否能进入外部存档EP,若满足,则返回步骤一,若不满足,进入步骤三。Step 2: Perform neighborhood operation on the selected solution to determine whether the solution after the neighborhood search operation can enter the external archive EP. If it is satisfied, return to step 1; if not, proceed to step 3.

步骤三,已搜索次数m进行加一操作,判断已搜索次数m是否等于D,若满足,则结束,若不满足,则返回步骤二。Step 3, add 1 to the number of searches m to determine whether the number of searches m is equal to D. If so, end the process; if not, return to step 2.

步骤S34:对各所述解和各所述解对应的搜索结果应用ε占优存档更新策略,得到出行外部存档。Step S34: applying the ε-dominant archive update strategy to each of the solutions and the search results corresponding to each of the solutions to obtain a travel external archive.

在实际应用中,所述ε占优存档更新策略包括:In practical applications, the ε-dominant archive update strategy includes:

如果Archive为空,则将产生的路径规划方案Xnew加入到Archive中。其中,Archive是多目标优化算法理论中已定义的词,表示存档的意思,即用于在算法优化过程中储存每一代产生的优秀个体,即非支配解。If Archive is empty, the generated path planning solution X new is added to Archive. Archive is a term defined in the theory of multi-objective optimization algorithms, which means archiving, that is, it is used to store the excellent individuals generated in each generation during the algorithm optimization process, that is, non-dominated solutions.

如果Archive非空,则将产生的路径规划方案Xnew与已有的路径规划方案进行ε占优比较;如果存在已有方案占优Xnew或者与Xnew相同,则将Xnew丢弃;如果Xnew占优已有的路径规划方案,则将被占优的方案全部删除,并将Xnew加入到Archive中;如果Xnew与所有的路径规划方案互不占优,则将Xnew加入到Archive中。If Archive is not empty, the generated path planning scheme X new is compared with the existing path planning schemes for ε-dominance; if there is an existing scheme that is superior to X new or is the same as X new , X new is discarded; if X new is superior to the existing path planning schemes, all the superior schemes are deleted and X new is added to Archive; if X new is not superior to all the path planning schemes, X new is added to Archive.

其中,ε占优用于存档更新。具体来说,存档中的每个非占优解都被分配了一个关联向量B={B1,B2,…,B4},4表示4个优化目标数,其中参数变量ε用于控制存档的规模。Specifically, each non-dominant solution in the archive is assigned an association vector B = {B 1 ,B 2 ,…,B 4 }, where 4 represents the number of optimization objectives. The parameter variable ε is used to control the size of the archive.

存档更新策略:在多目标城际网约车路径规划问题中,路径规划方案之间的比较是通过多目标占优关系来进行的。本发明所涉及的占优关系的定义如下:对于路径规划方案X和Y,如果Archive update strategy: In the multi-objective intercity online car-hailing route planning problem, the comparison between route planning schemes is carried out through the multi-objective dominance relationship. The definition of the dominance relationship involved in the present invention is as follows: For route planning schemes X and Y, if

1)对于所有的目标值,fj(X)≤fj(Y),j=1,2。1) For all target values, f j (X) ≤ f j (Y), j = 1, 2.

2)至少存在一个j,使得fj(X)<fj(Y)。2) There exists at least one j such that f j (X) < f j (Y).

同时满足以上两个条件,则称X占优Y;否则,则称X和Y互不占优,X和Y是非占优解。If the above two conditions are met at the same time, then X is said to dominate Y; otherwise, X and Y are said to be non-dominant solutions.

此外,在多目标占优关系的基础上,加入ε占优策略,用于限定搜索过程中外部存档的规模。在ε占优的存档中,存档中每个非占优解对应一个关联向量B={B1,B2,B3,B4},其中,Bi=log(fi+1)/log(1+ε),每个超立方体中存储一个非占优的解。因此,基于ε占优的存档不仅使非占优解均匀分布,而且限定了搜索过程中外部存档的规模。In addition, based on the multi-objective dominance relationship, an ε-dominance strategy is added to limit the size of the external archive during the search process. In the ε-dominant archive, each non-dominant solution in the archive corresponds to an association vector B = {B 1 , B 2 , B 3 , B 4 }, where Bi = log( fi +1)/log(1+ε), and each hypercube stores a non-dominant solution. Therefore, the archive based on ε-dominance not only makes the non-dominant solutions evenly distributed, but also limits the size of the external archive during the search process.

根据上述占优关系的定义,本发明所解决的多目标城际网约车路径规划问题的存档更新策略如下:如果Archive为空,则将产生的路径规划方案Xnew加入到Archive中。According to the definition of the above dominant relationship, the archive update strategy for the multi-objective intercity online car-hailing path planning problem solved by the present invention is as follows: if Archive is empty, the generated path planning solution X new is added to Archive.

如果Archive非空,则将产生的路径规划方案Xnew与已有的路径规划方案进行ε占优比较;如果存在已有方案占优Xnew或者与Xnew相同,则将Xnew丢弃;如果Xnew占优已有的路径规划方案,则将被占优的方案全部删除,并将Xnew加入到Archive中;如果Xnew与所有的路径规划方案互不占优,则将Xnew加入到Archive中。If Archive is not empty, the generated path planning scheme X new is compared with the existing path planning schemes for ε-dominance; if there is an existing scheme that is superior to X new or is the same as X new , X new is discarded; if X new is superior to the existing path planning schemes, all the superior schemes are deleted and X new is added to Archive; if X new is not superior to all the path planning schemes, X new is added to Archive.

在实际应用本发明时,获取待出行城际订单后,执行下述操作:When the present invention is actually applied, after obtaining the intercity travel order to be traveled, the following operations are performed:

(1)将两种订单集按订单最早出发时间升序排序,按方法生成满足约束的行程,完成两种行程集的初始化并对其按行程最早开始时间升序排序,按方法生成满足约束的车辆(出行任务,下同),完成两种车辆集的分配,生成一种初始可行的规划方案,并将该方案更新到外部存档解集EP中。(1) Sort the two order sets in ascending order according to the earliest departure time of the orders, generate a trip that satisfies the constraints according to the method, complete the initialization of the two trip sets and sort them in ascending order according to the earliest start time of the trip, generate vehicles (travel tasks, the same below) that meet the constraints according to the method, complete the allocation of the two vehicle sets, generate an initial feasible planning scheme, and update the scheme to the external archive solution set EP.

(2)判断初始化次数n是否是否达到设置的最大循环次数nmax,若达到,则进入步骤(3);否则,返回步骤(1)。其中,最大循环次数nmax为预先设定的值;可选地,最大循环次数nmax可以根据经验值或预实验决定。(2) Determine whether the initialization number n reaches the set maximum number of cycles n max . If so, proceed to step (3); otherwise, return to step (1). The maximum number of cycles n max is a preset value; optionally, the maximum number of cycles n max can be determined based on experience or preliminary experiments.

(3)对于外部存档解集EP中的路径规划方案,基于城际出行特点设计的八种邻域操作算子,进行自适应局部搜索策略更新EP。(3) For the path planning solutions in the external archive solution set EP, eight neighborhood operation operators designed based on the characteristics of inter-city travel are used to perform an adaptive local search strategy to update EP.

(4)判断运行时间h是否达到设置的最大运行时间hmax,若达到,则进入步骤(5),否则,返回步骤(3)。(4) Determine whether the running time h reaches the set maximum running time h max . If so, proceed to step (5); otherwise, return to step (3).

(5)按决策者的需求(根据四种目标的实际需要)选择一种外部存档解集EP中的路径规划方案,将方案中的每一条路径序列分配给对应该路径序列的车辆(司机)。(5) According to the decision maker's needs (based on the actual needs of the four objectives), a path planning solution in the external archive solution set EP is selected, and each path sequence in the solution is assigned to the vehicle (driver) corresponding to the path sequence.

由于当前被广泛使用的算例在用于城际出行路径优化问题存在两点缺陷:第一,数据集中客户点之间的行驶距离和行驶时间过于理想化,其依赖于欧几里得距离且假设满足三角不等式关系,但在现实世界中往往不成正比例关系且不一定满足三角不等式。第二,在该算例中,其目标之间呈现弱相关性,因此不能应用在多目标优化的研究中。由于缺少适合本研究面向城际出行的车辆路径问题测试算例,同时为了进一步验证本研究提出的模型和算法的有效性。因此本发明设计了适合本研究面向城际出行的车辆路径问题测试算例:基于某城际网约车平台的真实订单数据进行测试。其中测试算例具体定义如下:Because the currently widely used examples have two defects in the inter-city travel path optimization problem: first, the driving distance and driving time between customer points in the data set are too idealized. It relies on the Euclidean distance and assumes that the triangle inequality relationship is satisfied, but in the real world, it is often not in direct proportion and does not necessarily satisfy the triangle inequality. Second, in this example, there is a weak correlation between its objectives, so it cannot be applied in the study of multi-objective optimization. Due to the lack of test examples suitable for the vehicle path problem for inter-city travel in this study, and in order to further verify the effectiveness of the model and algorithm proposed in this study. Therefore, the present invention designs a test example for the vehicle path problem for inter-city travel suitable for this study: the test is based on the real order data of a certain inter-city online car-hailing platform. The test example is specifically defined as follows:

基于厦门某网约车平台实际场景的数据,选择厦门和漳州之间的实际订单数据设计了测试算例,客户点和车库的坐标基于实际位置生成,同时任意两点之间的行驶时间和行驶距离均基于实际路网计算得到,实际路网可以是百度地图,可以得到更加符合实际出行规律的数据。距离精确到米,时间精确到分钟。同时为了尽可能保证数据的真实性和随机性,订单的出发时间理应满足早晚高峰的实际情形,故按比例生成对应时间窗口内的随机出发时间,订单的乘车人数则由单纯随机的方式进行生成。此外,通过查阅文献和调研,拟通过变化订单数量、出行需求比例和时间窗类型定义了三类规模的场景,每类规模又由于出行需求比例的不同分成五种数量的客户点,同时根据三种时间窗类型的不同进一步将场景组合起来,可得到共45个测试实例。其最大客户点数量可达600个,故为大规模的测试算例。Based on the data of the actual scene of a certain online car-hailing platform in Xiamen, the actual order data between Xiamen and Zhangzhou was selected to design a test case. The coordinates of the customer point and the garage were generated based on the actual location. At the same time, the driving time and driving distance between any two points were calculated based on the actual road network. The actual road network can be Baidu Map, which can obtain data that is more in line with the actual travel rules. The distance is accurate to meters and the time is accurate to minutes. At the same time, in order to ensure the authenticity and randomness of the data as much as possible, the departure time of the order should meet the actual situation of the morning and evening peaks. Therefore, the random departure time within the corresponding time window is generated in proportion, and the number of passengers in the order is generated by a simple random method. In addition, by consulting literature and conducting research, it is proposed to define three types of scale scenarios by changing the number of orders, travel demand ratio and time window type. Each type of scale is divided into five numbers of customer points due to the different travel demand ratios. At the same time, the scenarios are further combined according to the different types of three time windows, and a total of 45 test instances can be obtained. The maximum number of customer points can reach 600, so it is a large-scale test case.

通过该测试算例进行动态地模拟,并与该公司专业调度客服进行人工挑战的路径规划方案进行对比,本发明提出的一种求解多目标城际出行车辆路径问题的自适应局部搜索方法在4个目标方面均有明显下降,并且在上座率方面和乘客评价反馈方面也有显著的提升。同时从其中取出算例的所有解决方案,将其得到的目标值进行两两对比的权衡图(即四个目标、两两对比的权衡比较,例如目标f1作为x轴,目标f2作为y轴),从结果图来看,各个目标之间关系均大致呈反比例函数,说明本模型选择的4个目标之间均存在负相关关系,验证了本发明所建立的多目标城际出行车辆路径问题模型的目标选择的正确性。除此之外,为了进一步验证提出的方法的有效性,将所提出的方法与LSMOVRPTW算法、MOEA/D算法、NSGA2算法的对比实验结果进行统计分析。并采用了多目标优化算法广泛使用的两个指标,即:IGD和HV。它们被用来从不同的角度表明算法的收敛性与多样性,其中IGD和HV是一元指标,结果表明,本发明方法在IGD和HV指标方面都明显优于LSMOVRPTW算法。具体来说,在在IGD指标方面,根据Wilcoxon检验的单问题分析,本发明方法在40个算例中显著优于LSMOVRPTW算法,1个算例中劣于LSMOVRPTW算法。在HV指标方面,本发明方法在37个算例中显著优于LSMOVRPTW算法,在3个算例中则较差一些。基于Wilcoxon检验的多问题分析,本发明方法在IGD和HV指标方面均获得高于R-的R+值。说明在所有rwMOVRPTW算例中,本发明方法总体上优于对比算法。同时为了直观地展示本发明方法和LSMOVRPTW算法的收敛性和多样性,在几个代表性算例中得到的非支配解在第二个目标与第三个目标(f2-f3)、第二个目标与第四个目标(f2-f4)上的映射。相较于对比算法,本发明方法的非支配解分布更加广泛,并且也更接近帕累托前沿。同时,针对本发明提出的八个邻域操作进行模块对比测试,分别省略了每个局部搜索算子,并运行算法30次。省略一个邻域操作算子后,每个实例的平均目标的减少。可以观察到,由于缺乏任何邻域操作算子,特别是当实例的节点数量增加时,本发明的方法的性能在目标方面变得更差。这意味着所有的邻域操作算子对于本发明所提出的算法是非常重要的。利用IGD、HV指标、Wilcoxon检验和Frediman测试,可以观察到,缺少的每个邻域操作算子后,都劣于缺失前的结果,其中邻域操作符L4、L5、L6和L8的影响更为显著。综上所述,本发明提出的方法能够高效智能地处理城际网约车出行的路径规划问题。Through dynamic simulation of the test case and comparison with the path planning scheme of the company's professional dispatch customer service for manual challenge, the adaptive local search method for solving the multi-objective intercity travel vehicle path problem proposed in the present invention has a significant decrease in the four objectives, and has also significantly improved in the occupancy rate and passenger evaluation feedback. At the same time, all the solutions of the example are taken out, and the obtained target values are compared in pairs (i.e., four targets, two-by-two comparisons, for example, target f1 as the x-axis, target f2 as the y-axis). From the result diagram, the relationship between each target is roughly an inverse proportional function, indicating that there is a negative correlation between the four targets selected by this model, which verifies the correctness of the target selection of the multi-objective intercity travel vehicle path problem model established by the present invention. In addition, in order to further verify the effectiveness of the proposed method, the comparative experimental results of the proposed method with the LSMOVRPTW algorithm, MOEA/D algorithm, and NSGA2 algorithm are statistically analyzed. And two indicators widely used in multi-objective optimization algorithms are adopted, namely: IGD and HV. They are used to demonstrate the convergence and diversity of the algorithm from different perspectives, among which IGD and HV are univariate indicators. The results show that the method of the present invention is significantly better than the LSMOVRPTW algorithm in terms of both IGD and HV indicators. Specifically, in terms of the IGD indicator, according to the single-question analysis of the Wilcoxon test, the method of the present invention is significantly better than the LSMOVRPTW algorithm in 40 examples, and worse than the LSMOVRPTW algorithm in 1 example. In terms of the HV indicator, the method of the present invention is significantly better than the LSMOVRPTW algorithm in 37 examples, and worse in 3 examples. Based on the multi-question analysis of the Wilcoxon test, the method of the present invention obtains R+ values higher than R- in both IGD and HV indicators. It shows that in all rwMOVRPTW examples, the method of the present invention is generally better than the comparison algorithm. At the same time, in order to intuitively demonstrate the convergence and diversity of the method of the present invention and the LSMOVRPTW algorithm, the mapping of the non-dominated solutions obtained in several representative examples on the second objective and the third objective (f 2 -f 3 ), the second objective and the fourth objective (f 2 -f 4 ) is shown. Compared with the comparison algorithm, the non-dominated solutions of the method of the present invention are more widely distributed and closer to the Pareto front. At the same time, a module comparison test is performed on the eight neighborhood operations proposed by the present invention, each local search operator is omitted respectively, and the algorithm is run 30 times. After omitting a neighborhood operation operator, the average target of each instance is reduced. It can be observed that due to the lack of any neighborhood operation operator, especially when the number of nodes of the instance increases, the performance of the method of the present invention becomes worse in terms of the target. This means that all neighborhood operation operators are very important for the algorithm proposed by the present invention. Using IGD, HV indicators, Wilcoxon test and Frediman test, it can be observed that the results after each missing neighborhood operation operator are inferior to the results before the missing, among which the influence of neighborhood operators L4, L5, L6 and L8 is more significant. In summary, the method proposed in the present invention can efficiently and intelligently handle the path planning problem of inter-city online car-hailing travel.

本发明将城际网约车路径优化问题定义为一个包含四个目标的多目标问题,更加全面真实地反映城际出行问题的本质;通过基于时间序列的初始化构造方法构造初始非占优的车辆路径方案;然后基于城际出行特点设计的八种邻域操作算子,利用局部搜索的收敛性和解方案的多样性特性,制定解的选择策略和局部搜索选择策略,利用自适应局部搜索策略对路径规划方案进行迭代优化,实现利用多个邻域操作之间潜在的协同作用进行有效配合,同时,实现多目标解的质量的进一步提升。这些机制的有效结合,不仅能够高效地规划城际网约车的出行路径,而且可以利用多目标优化方法为城际网约车服务提供不同需求的高质量的规划方案集合。The present invention defines the inter-city online car-hailing route optimization problem as a multi-objective problem with four objectives, which more comprehensively and truly reflects the essence of the inter-city travel problem; constructs an initial non-dominant vehicle route plan through a time series-based initialization construction method; then, eight neighborhood operation operators designed based on the characteristics of inter-city travel, utilize the convergence of local search and the diversity of solution plans, formulate solution selection strategies and local search selection strategies, and use adaptive local search strategies to iteratively optimize the route planning plan, so as to achieve effective coordination using the potential synergy between multiple neighborhood operations, and at the same time, further improve the quality of multi-objective solutions. The effective combination of these mechanisms can not only efficiently plan the travel routes of inter-city online car-hailing, but also use multi-objective optimization methods to provide a high-quality set of planning solutions for different needs for inter-city online car-hailing services.

实施例二Embodiment 2

为了执行上述实施例一对应的方法,以实现相应的功能和技术效果,下面提供一种城际出行车辆路径确定系统,所述确定系统包括:In order to execute the method corresponding to the above embodiment 1 to achieve corresponding functions and technical effects, a system for determining a path for an inter-city travel vehicle is provided below, and the determination system includes:

构建模块,用于以最大化平均每趟出行的乘客数、最小化车辆数目、最小化车辆总行驶距离和最小化乘客总等待时间为目标函数,以车辆载客量、服务质量、时间约束和安全约束为约束条件,构建城际出行车辆路径模型。The construction module is used to construct an inter-city travel vehicle path model with the objective function of maximizing the average number of passengers per trip, minimizing the number of vehicles, minimizing the total vehicle travel distance and minimizing the total passenger waiting time, and with the vehicle passenger capacity, service quality, time constraints and safety constraints as constraints.

获取模块,用于获取待出行城际订单;所述待出行城际订单包括多个出行订单;各所述出行订单包括乘车人数、出发时间、出发地点和目的地点。The acquisition module is used to obtain intercity travel orders to be traveled; the intercity travel orders to be traveled include multiple travel orders; each of the travel orders includes the number of passengers, departure time, departure location and destination point.

方案确定模块,用于根据所述城际出行车辆路径模型,确定各所述出行订单对应的出行车辆;其中,一个所述出行车辆和由所述出行车辆完成的多个所述出行订单构成一个出行解;多个所述出行解构成出行外部存档;将所述出行外部存档作为所述待出行城际订单的路径规划方案。A solution determination module is used to determine the travel vehicles corresponding to each travel order according to the inter-city travel vehicle path model; wherein, one travel vehicle and multiple travel orders completed by the travel vehicle constitute a travel solution; multiple travel solutions constitute a travel external archive; and the travel external archive is used as the path planning solution for the inter-city order to be traveled.

实施例三Embodiment 3

本发明实施例提供一种电子设备,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行计算机程序以使电子设备执行实施例一的城际出行车辆路径确定方法。An embodiment of the present invention provides an electronic device, including a memory and a processor, wherein the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the inter-city travel vehicle path determination method of embodiment 1.

可选地,上述电子设备可以是服务器。Optionally, the above-mentioned electronic device may be a server.

另外,本发明实施例还提供一种计算机可读存储介质,其存储有计算机程序,该计算机程序被处理器执行时实现实施例一的城际出行车辆路径确定方法。In addition, an embodiment of the present invention further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the inter-city travel vehicle path determination method of embodiment 1.

本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。This article uses specific examples to illustrate the principles and implementation methods of the present invention. The above examples are only used to help understand the method and core ideas of the present invention. At the same time, for those skilled in the art, according to the ideas of the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims (9)

1. An inter-city travel vehicle path determining method, comprising:
Constructing an inter-city travel vehicle path model by taking the maximized average number of passengers per trip, the minimized number of vehicles, the minimized total travel distance of the vehicles and the minimized total waiting time of the passengers as objective functions and taking the vehicle passenger capacity, the service quality, the time constraint and the safety constraint as constraint conditions;
acquiring an inter-city order to be travelled; the inter-city orders to be travelled comprise a plurality of travelling orders; each travel order comprises the number of passengers, departure time, departure place and destination place;
determining travel vehicles corresponding to the travel orders according to the inter-city travel vehicle path model; wherein one travel vehicle and a plurality of travel orders completed by the travel vehicle form a travel solution; a plurality of travel solutions form a travel external archive; and taking the trip external archive as a path planning scheme of the inter-city order to be tripped.
2. The inter-city travel vehicle path determining method of claim 1, wherein the objective function is:
f 2 =|M k |
wherein the object f 1 Representing the maximum average number of passengers per trip, M t 、M k Respectively represent the number of the travel t and the vehicle k, Q t Indicating the seating rate of the stroke t, the target f 2 Representing a minimized number of vehicles, target f 3 Representing a minimum total travel distance of the vehicle, M t Indicating the number of strokes t, l t The total travel distance of the journey t is indicated,for determining whether the journey t is to be routed to the vehicle k, the target f 4 Meaning minimizing the total waiting time of the passengers.
3. The inter-city travel vehicle path determining method according to claim 1, wherein determining a travel vehicle corresponding to each travel order according to the inter-city travel vehicle path model specifically comprises:
determining a plurality of travel waiting trips corresponding to each travel order according to the constraint conditions and the objective function; wherein each travel to be traveled includes at least one of the travel orders and each travel to be traveled satisfies the constraint condition and the objective function;
determining vehicles corresponding to the travel routes to be traveled according to the constraint conditions and the objective function; wherein one of the vehicles and a plurality of the travel routes to be traveled completed by the vehicle form a solution; a plurality of said solutions constitute an external archive;
applying a neighborhood operator to each solution to obtain a search result corresponding to each solution;
and applying epsilon-dominant archive update strategies to each solution and the search results corresponding to each solution to obtain the travel external archive.
4. The inter-city travel vehicle path determining method according to claim 3, wherein determining a plurality of travel routes to be traveled corresponding to each of the travel orders according to the constraint condition and the objective function, specifically comprises:
taking the travel order with the earliest departure time in the preprocessing order set as a reference order of an initial journey; the initial journey at least comprises one travel order; the preprocessing order set comprises a plurality of travel orders with the same departure place;
selecting a preset pre-processing order from the to-be-processed order set, and adding the preset pre-processing order into the initial journey to obtain an updated initial journey and an order queue; the to-be-processed order set is the rest order after the reference order is removed from the pre-processing order set; the order queue is the rest order after one preset pre-processing order is removed from the to-be-processed order set;
judging whether the updated initial travel meets the constraint condition or not;
when the updated initial journey does not meet the constraint condition, deleting the preset pre-processing order from the updated initial journey, and returning to the step of selecting one preset pre-processing order from the to-be-processed order set to be added into the initial journey to obtain an updated initial journey and an order queue, and continuing to execute;
When the updated initial journey meets the constraint condition, judging whether the number of passengers in the updated initial journey reaches the passenger capacity of the vehicle;
when the number of passengers in the updated initial journey reaches the passenger capacity of the vehicle, the updated initial journey is used as a final initial journey, and travel orders in the final initial journey are deleted from the preprocessing order set, so that an updated preprocessing order set is obtained;
when the number of passengers in the updated initial journey does not reach the passenger capacity of the vehicle, judging whether the number of travel orders in the order queue is 0;
when the number of the travel orders in the order queue is not 0, the order queue is used as an order set to be processed, and a step of selecting a preset pre-processing order from the order set to be processed and adding the preset pre-processing order into the initial journey is returned to obtain an updated initial journey and the order queue;
when the number of the travel orders in the order queue is 0, returning to the step of taking the updated initial journey as the final initial journey, and deleting the travel orders in the final initial journey from the preprocessing order set to obtain an updated preprocessing order set;
When the number of the travel orders in the updated pretreatment order set is not 0, the updated pretreatment order set is used as a pretreatment order set, and the step of continuously executing the travel order with the earliest departure time in the pretreatment order set as a reference order of an initial journey is returned;
and when the number of the updated pre-processing order sets for the travel orders is 0, determining the final initial travel as a plurality of travel waiting travels corresponding to the travel orders.
5. The inter-city travel vehicle path determining method according to claim 3, wherein determining vehicles corresponding to each of the travel routes to be traveled according to the constraint condition and the objective function, specifically comprises:
taking the travel route to be traveled with the earliest departure time in the preprocessing route set as a reference route of the vehicle; the vehicle at least comprises one travel route to be traveled; the pretreatment travel set comprises a plurality of travel to be performed travel routes with the same departure place;
selecting a preset pretreatment travel from the travel set to be treated, and adding the preset pretreatment travel set to the vehicle to obtain an updated vehicle and travel queue; the stroke set to be processed is the rest stroke of the preprocessing stroke set after the reference stroke is removed; the stroke queue is the residual stroke after one preset pretreatment stroke is removed from the to-be-treated stroke set;
Judging whether the updated vehicle meets the constraint condition or not;
when the updated vehicle does not meet the constraint condition, deleting the preset pretreatment travel from the updated vehicle, and returning to the step of selecting one preset pretreatment travel from the to-be-treated travel set to be added into the vehicle to obtain an updated vehicle and a travel queue, and continuously executing;
when the updated vehicle meets the constraint condition, the updated vehicle is used as a final vehicle, and the travel route to be traveled in the final vehicle is deleted from the pretreatment route set, so that an updated pretreatment route set is obtained;
judging whether the number of travel strokes to be performed in the stroke queue is 0;
when the number of travel routes to be traveled in the route queue is not 0, the route queue is used as a route set to be processed, and a step of selecting a preset pretreatment route from the route set to be processed and adding the preset pretreatment route into the vehicle is returned to obtain an updated vehicle and route queue;
when the number of travel routes to be traveled in the route queue is 0, returning to the step of selecting a preset pretreatment route from the travel route set to be treated and adding the preset pretreatment route into the initial route to obtain an updated vehicle and route queue;
Judging whether the number of travel routes to be traveled in the updated pretreatment route set is 0;
when the number of travel routes to be traveled in the updated pretreatment route set is not 0, the updated pretreatment route set is used as a pretreatment route set, and the step of continuously executing the travel route to be traveled with the earliest departure time in the pretreatment route set as a reference route of the vehicle is returned;
and when the number of the travel routes to be traveled in the updated pretreatment route set is 0, determining that the final vehicle is a plurality of travel routes to be traveled corresponding to each travel route to be traveled.
6. The method for determining an inter-city travel vehicle path according to claim 1, wherein applying a neighborhood operator to each solution to obtain a search result corresponding to each solution, comprises:
searching each solution by applying a plurality of neighborhood operators to obtain a plurality of initial search results corresponding to each solution;
calculating the lifting degree corresponding to each solution according to a plurality of initial search results;
when the searching times are smaller than the preset times, returning to the step of applying a neighborhood operator to the preset solutions in the solutions to obtain initial searching results, and continuing to execute the steps;
When the searching times are greater than or equal to preset times, calculating an improved crowding distance of each solution according to each solution and the initial searching result corresponding to each solution;
normalizing the improved crowding distance of each solution to obtain the solution selection probability of each solution;
determining the operation selection probability of each neighborhood operator according to the lifting degree corresponding to each solution;
selecting a solution to be processed from a plurality of solutions according to the solution selection probability;
determining target neighborhood operation from a plurality of neighborhood operation operators according to the operation selection probability;
applying the target neighborhood operation to perform local search on the solution to be processed to obtain a final search result corresponding to the solution to be processed;
and updating the initial search result according to the final search result to obtain the search result corresponding to each solution.
7. An intercity travel vehicle path determination system, the determination system comprising:
the construction module is used for constructing an inter-city travel vehicle path model by taking the maximized average number of passengers per trip, the minimized number of vehicles, the minimized total running distance of the vehicles and the minimized total waiting time of the passengers as objective functions and taking the vehicle passenger capacity, the service quality, the time constraint and the safety constraint as constraint conditions;
The acquisition module is used for acquiring an inter-city order to be travelled; the inter-city orders to be travelled comprise a plurality of travelling orders; each travel order comprises the number of passengers, departure time, departure place and destination place;
the scheme determining module is used for determining the travel vehicles corresponding to the travel orders according to the inter-city travel vehicle path model; wherein one travel vehicle and a plurality of travel orders completed by the travel vehicle form a travel solution; a plurality of travel solutions form a travel external archive; and taking the trip external archive as a path planning scheme of the inter-city order to be tripped.
8. An electronic device comprising a memory for storing a computer program and a processor that runs the computer program to cause the electronic device to perform the inter-urban travel vehicle path determination method according to any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that it stores a computer program which, when executed by a processor, implements the inter-urban travel vehicle path determination method according to any one of claims 1 to 6.
CN202311209584.9A 2023-09-19 2023-09-19 Intercity travel vehicle route determination method, system, electronic device and medium Pending CN117132011A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311209584.9A CN117132011A (en) 2023-09-19 2023-09-19 Intercity travel vehicle route determination method, system, electronic device and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311209584.9A CN117132011A (en) 2023-09-19 2023-09-19 Intercity travel vehicle route determination method, system, electronic device and medium

Publications (1)

Publication Number Publication Date
CN117132011A true CN117132011A (en) 2023-11-28

Family

ID=88852699

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311209584.9A Pending CN117132011A (en) 2023-09-19 2023-09-19 Intercity travel vehicle route determination method, system, electronic device and medium

Country Status (1)

Country Link
CN (1) CN117132011A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117808273A (en) * 2024-02-29 2024-04-02 华侨大学 Inter-city carpooling scheduling method and device for passenger departure time cooperation and stage feedback

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117808273A (en) * 2024-02-29 2024-04-02 华侨大学 Inter-city carpooling scheduling method and device for passenger departure time cooperation and stage feedback
CN117808273B (en) * 2024-02-29 2024-05-31 华侨大学 Inter-city carpooling scheduling method and device for passenger departure time cooperation and stage feedback

Similar Documents

Publication Publication Date Title
CN108764777B (en) Electric logistics vehicle scheduling method and system with time window
CN105070044B (en) Dynamic scheduling method for customized buses and car pooling based on passenger appointments
CN109948855A (en) A transport route planning method for heterogeneous hazardous chemicals with time window
CN112466122B (en) Method and device for generating alternative line set and planning line of public traffic line network
CN109214756B (en) Vehicle logistics scheduling method and device, storage medium and terminal
CN113485429B (en) Route optimization method and device for air-ground cooperative traffic inspection
CN111144618B (en) Two-stage demand response type customized public transportation network planning method
Grzybowska et al. Decision support system for real-time urban freight management
CN109948854B (en) An order allocation method for intercity car-hailing based on multi-objective optimization
CN112308389A (en) An orderly charging scheduling system and method for electric vehicles based on cloud computing
Hou et al. TASeT: Improving the efficiency of electric taxis with transfer-allowed rideshare
CN113642761A (en) Robotaxi automatic driving shared network car booking resource allocation method
Lowalekar et al. Zone path construction (zac) based approaches for effective real-time ridesharing
CN109934380A (en) Optimization method of shared electric vehicle vehicle and personnel scheduling based on two-level programming
CN117132011A (en) Intercity travel vehicle route determination method, system, electronic device and medium
CN116090931A (en) A terminal distribution method and device based on customer classification
CN111860957B (en) Multi-vehicle-type vehicle path planning method considering secondary distribution and balancing
Sun et al. Risk-aware operation modeling for ride-hailing fleet in order grabbing mode: A distributional reinforcement learning approach
CN106845727A (en) Consider the highway quick charge station heuristic programming algorithm of distribution trend constraint
CN112149906B (en) Comprehensive optimization method for travel line of electric vehicle considering charging time
He et al. The charging station and swapping station site selection with many-objective evolutionary algorithm
CN111651899A (en) Robust location and capacity determination method and system for battery swap stations considering user selection behavior
CN113761397B (en) Recommendation method, system, equipment and storage medium for customizing passenger transport route
Xiao et al. A reinforcement-learning-based bus scheduling model
Viana et al. Optimization of a demand responsive transport service using multi-objective evolutionary algorithms

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination