Research Article  Open Access
Shaopei Chen, Ji Yang, Yong Li, Jingfeng Yang, "Multiconstrained Network Intensive Vehicle Routing Adaptive Ant Colony Algorithm in the Context of Neural Network Analysis", Complexity, vol. 2017, Article ID 8594792, 9 pages, 2017. https://doi.org/10.1155/2017/8594792
Multiconstrained Network Intensive Vehicle Routing Adaptive Ant Colony Algorithm in the Context of Neural Network Analysis
Abstract
Neural network models have recently made significant achievements in solving vehicle scheduling problems. Adaptive ant colony algorithm provides a new idea for neural networks to solve complex system problems of multiconstrained network intensive vehicle routing models. The pheromone in the path is changed by adjusting the volatile factors in the operation process adaptively. It effectively overcomes the tendency of the traditional ant colony algorithm to fall easily into the local optimal solution and slow convergence speed to search for the global optimal solution. The multiconstrained network intensive vehicle routing algorithm based on adaptive ant colony algorithm in this paper refers to the interaction between groups. Adaptive transfer and pheromone update strategies are introduced based on the traditional ant colony algorithm to optimize the selection, update, and coordination mechanisms of the algorithm further. Thus, the search task of the objective function for a feasible solution is completed by the search ants. Through the division and collaboration of different kinds of ants, pheromone adaptive strategy is combined with polymorphic ant colony algorithm. It can effectively overcome some disadvantages, such as premature stagnation, and has a theoretical significance to the study of largescale multiconstrained vehicle routing problems in complex traffic network systems.
1. Introduction
Network intensive vehicle service, which is an important part of the complex urban public transport network system, is an innovation representative of the shared economy model in the era of Internet +. The Internet service platform has been attracting increasing attention in the shared economy era by creating a new commercial and management mode to adapt the Internet economy. It is also perceived as a representative of the development direction of the new economy and business model. Moreover, optimizing the network intensive vehicle scheduling model and improving the scheduling efficiency, response speed, and cost savings are important measures to improve the quality of network intensive vehicle services. Thus, studies on supplyanddemand information matching and of the path optimization algorithm for developing network intensive vehicle services are of great significance. However, the vehicle scheduling problem on network intensive car rental is different from that of the classic vehicle scheduling problem. In the classic vehicle scheduling problem, the control center is assumed to know all the information related to optimal scheduling before the optimal scheduling instruction is executed, and the information does not change over time. In the era of Internet +, passenger demand information, such as location, time of demand, quantity demanded, travel time, and service time, is unknown and even dynamic before the path optimization. Furthermore, new information may have arrived or existing information may have changed. Therefore, dynamic vehicle routing problems (DVRPs) are more prominent in the network intensive vehicle scheduling problem. Moreover, the supplyanddemand information of network intensive vehicles is matched in real time seamlessly and dynamically and does not require manual docking and transmission. This overcomes a problem that network booking taxi is forbidden to parade showmanship on the road under the existing management system. Therefore, the significant difference between the network intensive and classic vehicle scheduling problems is the continuous development of and change in the demand type and quantity in the former based on the realtime interaction through the Internet. In addition, realtime dynamic resource scheduling is performed according to the change in demand. However, the first principle of passenger transport in network intensive vehicle scheduling is to minimize the response time and relatively lessen the cost factors. This is also an important feature that distinguishes it from the traditional transportation problem and the vehicle routing problem (VRP). Network intensive vehicle scheduling problem is a multiconstrained path planning problem that integrates the dynamic network flow problem with multiconditioned and multistartandstop points.
Studies on vehicle scheduling problem at home and abroad have gradually deepened in recent years. Studies on the target selection of modeling and consideration of constraints are mostly based on the simple extension of the traditional transportation problem and VRP. The single starting point and demand problem is extended to multiple starting points and demand problems. Moreover, the goal is extended from the single focus on cost minimization to the multiple objective optimizations of efficiency and cost [1–3]. Presently, no literature on the multiconstrained path problem of network intensive vehicles has been found, although some scholars have recently started to introduce some particular constraints into vehicle scheduling modeling. For instance, some scholars have established many kinds of dynamic vehicle scheduling models and designed concise and efficient heuristic algorithms for goods distribution or emergency material transportation vehicle scheduling problems [4–7]. Emergency material transportation vehicle scheduling models consider multitransport ways, multimaterials, multisupply points, and multidemand points, as they suppose that the number of vehicles, supply amount of materials and goods, and other parameters are dynamic. However, the goal of emergency material transportation is to minimize the total quantity of the goods that are not satisfied in the transport cycle. Although constraint is comprehensive in emergency material transportation vehicle scheduling models, vehicle scheduling under a dynamic demand still adopts the optimization method in a single cycle and limits the consideration of global optimization in multiple cycles [6, 7]. Other scholars similarly consider the vehicle dynamic scheduling problem under a vehicle shortage hypothesis and introduce minimum transportation time and cost into the goal. Nevertheless, the improvement of the vehicle scheduling efficiency is still limited to the optimal one in a single stage. The global consideration of dynamic vehicle scheduling problem is inadequate [8, 9]. Thus, many unsolved problems are elucidated in the research on network intensive vehicle scheduling problem. Therefore, considering the characteristics of DVRP with multiple constraints (DVRPMC), the construction and exploration of a simple intelligent algorithm with excellent performance are important for the construction of a network intensive vehicle scheduling service system.
2. DVRPMC Description and Mathematical Model
2.1. Problem Description
A network intensive vehicle scheduling problem is oriented to the personal travel mode. Therefore, it should involve objective and subjective conditions and is a typical vehicle routing problem with multiple constraints, involving time, vehicle, traffic, safety, and so on. As a result, a network intensive vehicle schedule problem, that is, DVRPMC, should consider passenger priority in the order of booking time, demand of vehicle pattern, billing method, traffic conditions, time window, capacitance, safe surveillance, and other conditions. This is to ensure that network intensive scheduling is universal and practical. The detailed problem description is as follows. (1) A car is assigned to the passenger (i.e., customer point) according to the nearby principle. (2) The nearest dispatch point (i.e., parking spot of booking car) is selected from a total of callable vehicles based on the time window limit to reach the customer point. If the vehicle arrives in advance, it must wait to serve the customer. (3) Different priorities are assigned according to passenger booking time. The highest grade is 1, and the lowest grade is R. High priority means that the passenger has significant weight. Thus, all aspects of priority should be satisfied. (4) Vehicles only consider passenger demands that are less than or equal to their carrying capacities. (5) Each passenger can select only one type of vehicle and charging method pattern and is only serviced by the selected car. (6) The passenger arrives at the target point and guides the vehicle to the nearest scheduled point. (7) For a transport service, the proper vehicle type, dispatch time, and path can minimize the time cost and obtain the highest response efficiency.
Information in DVRPMC changes over time. DVRPMC should consider time factors. Time axis is introduced here. The whole scheduling cycle (such as a working day) in a dynamic scheduling environment is constructed as the time axis. The time each new demand triggers a cycle is . At , demand information is divided into four kinds according to the state of the vehicle: (1) vehicles that have completed the task, (2) vehicles in service or on the way to serve, (3) vehicles that have responded to the service demand but do not proceed to the demand spots, and (4) vehicles waiting for new demands. A series of key points can be constructed on the time axis to represent the dynamic problem, and a vehicle scheduling system can arrange a scheduling plan according to the identified key points. The following scheduling strategies are adopted. (1) The vehicle from a certain scheduling point returns to the nearest scheduling point at the end of the delivery. (2) The fairness principle is followed to prevent customer points from obtaining any vehicle service response. Each customer service response shall not be lower than the level of its own demand. (3) The goal of the network intensive vehicle transport service is to minimize the response time, that is, from the initial start of the demand to the time when the vehicle reaches the demand point.
Compared with traditional VRP, the essence of DVRPMC is to create a vehicle scheduling plan for each customer point in the case of a known customer point distribution and the realtime development of any customer point. Moreover, many different vehicle models should be considered in DVRPMC. Thus, dynamic model selection and path arrangement are performed according to the passenger capacities of the models, the effect on the path, and the model requirements of passengers.
2.2. Mathematical Model
The mathematical model needs to satisfy the dynamic scheduling description. Therefore, the analysis mechanics can be described as follows: Firstly, this model is set to be , which refers to the set of all key points at (Class 4). It is used to represent the dynamic problem by constructing a series of key points on the time axis. refers to the completed transport task, where the vehicle returned to the nearest scheduling point (parking point). refers to a passenger being served or a vehicle on the way to the service. refers to the vehicle that responded to the service demand but is not yet at the customer point. represents a vehicle in waiting state. Class 2 in the implementation of the task cannot be changed. The other three classes, as schedulable network intensive vehicles, are recorded as vehicles scheduled from the scheduling point (parking point) to customers at . Then, the scheduling task can be represented by a weighted graph at , where refers to a set of nodes, 0 refers to the scheduling point, nodes refer to customer points, and refers to a set of paths from to . Last, considering road conditions and various factors, a coefficient of road condition is used to represent the influence of road conditions on the scheduled vehicles. For the standard path, . If is better than the standard path, ; otherwise, . The coefficient of the road condition is multiplied by the actual length of the path between distribution points. is equal to the equivalent path length in consideration of the path influence. refers to a feasible path. refers to the cost corresponding to this path, that is, objective function. refers to the quantity booked by the passenger. refers to the maximum load of vehicle . The importance of the priority of a customer point is represented by weight . The time window of passenger is . The model demand of passenger is . The billing method of passenger is . refers to the travel time of vehicle from the scheduling point to the customer point. The time vehicle served at the starting point of passenger is . is a decision variable that represents the moment vehicle reaches customer point . is also a decision variable and can be represented as follows:
DVRPMC aims to determine a set of paths. The travel time of each path should be minimal, and no vehicle should exceed its carrying capacity . A feasible solution under a different priority, model demand, different billing method, and road conditions is sought for highly efficient vehicle scheduling at minimum cost.
First, considering the cost of vehicle scheduling, refers to the fixed cost of using vehicle . refers to the operating cost of vehicle in the unit distance. refers to the travel speed of vehicle under standard road conditions. Scheduling cost mainly includes the fixed cost of vehicle and operating cost . , and . Second, the cost for waiting for customer is considered, and . refers to the wait cost in the unit of time. Finally, extra cost is incurred if vehicle arrives beyond the time window of the passenger. To set linear increasewhere refers to the wait cost of vehicle in customer point with priority per unit time and refers to the penalty per unit time if vehicle arrives after the time window of the customer with priority .
Therefore, the minimum objective function can be represented as
Constraint conditions are as follows:
Equation (3) refers to the objective function, which represents the minimum scheduling cost. Equation (4) refers to the sum of load capacities on a feasible path, which is not more than the maximum load of the vehicle. Equation (5) refers to a vehicle (vehicle model and billing method that satisfies the requirements) that is distributed once in each customer point. Equations (6)–(8) refer to the subpath of each vehicle. A vehicle starts from the scheduling point, arrives at the customer point to provide transport service, brings passengers to the objective point, and finally returns to the scheduling point. Equation (9) indicates that if a vehicle travels directly from node to , the arrival time at node is equal to the sum of the arrival time node , service time, and travel time from node to . Equation (10) indicates that path distance is equal to the product of standard speed and travel time under the influence of road conditions. Equation (11) refers to the decision variable , which is .
3. Design of Optimization Algorithm of DVRPMC Based on Adaptive Ant Colony Algorithm
Intelligent algorithms have made great achievements in solving vehicle scheduling problems, especially for adaptive ant colony algorithm. Adaptive ant colony algorithm provides not only a new idea for solving the complex combination of optimization problems, but also a scientific perspective to deal with VRPs with multiple constraints. The pheromone in the path is changed by adaptively adjusting the volatile factors in the operation process. It effectively overcomes the tendency of the traditional ant colony algorithm to fall easily into the local optimal solution and slow convergence speed to search for the global optimal solution [10, 11]. Considering DVRPMC, this paper proposes an adaptive ant colony algorithm to introduce adaptive transfer and pheromone update strategies based on a traditional ant colony algorithm to optimize the selection, update, and coordination mechanisms of the algorithm further. Thus, considering the diversity of the ant colony, a task with many constraint conditions is provided to reconnaissance ants. Therefore, the search task for the feasible solution for the objective function is completed by the search ants. Through the division and collaboration of different kinds of ants, pheromone adaptive strategy is combined with polymorphic ant colony algorithm to effectively overcome some disadvantages, including long computing time and susceptibility to premature stagnation.
3.1. Principle and Implementation of Basic Ant Colony Algorithm
Ant colony algorithm is proposed based on the research on the real ant colony behavior in the natural world. It is a kind of simulated evolutionary algorithm based on population and belongs to the random search algorithm. Dorigo et al. [12, 13] first proposed this method and made full use of the similarities between the food search process of an ant colony and the famous traveling salesman problem (TSP) to determine the shortest path from the ant nest to the food source through information exchange and cooperation among individuals to solve the TSP [14–18]. The principle and the method of the basic algorithm are as follows.
ants are placed on nodes selected at random. An ant selects the next node or cycle it has not yet visited based on a criterion, , which is the concentration of residual information on the path from node to node at . This is the information provided by the algorithm itself, and is the initial information from node transferred to node . This initial information is provided with the problem to be solved. refers to an a priori value of node to . Thus, the probability that ant at node selects node as the objective node at issupposing that .
is the relative importance of the residual information and is the relative importance of the expected value.
refers to all possible objective nodes. Nodes mean no access. Each ant maintains a list , which records all cities it has visited up to the present, to avoid multiple visits to the same node. refers to the probability that ant transfers from node to node .
After each ant has accessed all nodes (i.e., a cycle), residual information must be updated and old information must be weakened to prevent residual information from inundating inspired information caused by excessive residual information. New information on the ant access path must be added to .
refers to the reserve part of residual information. refers to the weakened part of the residual information. must be less than 1 to prevent the unlimited accumulation of information. refers to the residual information content ant left in the path from to in the access period to .
GambardeIIa and Dorigo et al. provided the AntQ (Quantity) algorithm model based on the basic ant colony algorithm [12]. AntQ algorithm is a reinforcement learning algorithm based on the Cellular Computing paradigm. AntQ’s results on the vehicle routing problem, which is conceptually similar to fuel reload, are better than other ant colony optimization algorithms, such as AntC (Cycle) and AntD (Density) or the genetic algorithm (GA).
Otherwise, if ant selects the path () in the period from to (), is a constant, and refers to the total path length ant selected in this cycle. If it does not select this path,
According to the AntQ algorithm concept, each step (i.e., from to ()) requires the update of the concentration of residual information, but not the update of the residual information concentration after waiting for all ants to complete the access to all cities. Therefore, in the AntQ algorithm model, , where refers to the distance from node to node and is the concentration of residual information; that is, the concentration of residual information increases with the decrease of the distance to a city.
Therefore, the basic ant colony algorithm is mainly based on the combination of the principle of positive information feedback and a certain heuristic algorithm. This algorithm employs the random selection strategy in the process of constructing the solution. This selection strategy slows down the evolution. The positive feedback principle is designed to enhance the performance of the solution but is prone to stagnation. This is the root of the deficiency of the ant colony algorithm. Therefore, improving the selection strategy is necessary; that is, the selection strategy should be adapted based on the combination of deterministic and random selections. The evolution to a certain algebra dynamically adjusts the probability of deterministic selection in the search process. Thus, the evolution direction is basically determined. Dynamic adjustment at a particular time is based on the amount of information on the path [5]. Therefore, adaptive ant colony algorithm has increasingly attracted attention. This algorithm shrinks the gap in the quantity of information between the best and worst paths and appropriately enlarges the probability of random selection to conduct a complete search of the solution space with less than to overcome the shortcomings of the basic ant colony algorithm effectively.
3.2. Principle and Design of Adaptive Ant Colony Algorithm
Drawing lessons from the polymorphic ant colony method of Yang et al. [1], this paper divides ants in the ant colony algorithm into two types, namely, reconnaissance and search ants, because of the multiconstrained and dynamic characteristics of network intensive vehicle routing problems. Reconnaissance ants complete local reconnaissance and search ants complete global search.
Reconnaissance ants are endowed with multiconstrained tasks by the polymorphic ant colony method. Each customer point is regarded as the center to make local reconnaissance. Furthermore, a reconnaissance element is used to mark the reconnaissance result, in order to provide supplementary information for a search ant, and select the next customer point after arriving at this point. The local search method of the reconnaissance ant colony places reconnaissance ants in passenger demand points, and each reconnaissance ant regards its location as the center from which to scout the feasibility of points. The greater the feasibility, the higher the reconnaissance element on the path. Moreover, the reconnaissance element consists of three parts: (1) capacityconstrained reconnaissance element with a weight coefficient of , (2) time window matching reconnaissance element with a weight coefficient of , and (3) reconnaissance element in several points nearest to reconnaissance center (i.e., 20 points are taken from 100 demand points) with a weight coefficient of , and , where . If the capacity constraint is satisfactory, the contribution to the reconnaissance element is recorded as . If it is unsatisfactory, the contribution to the reconnaissance element is 0. The matching factor of the time window in the routing starting point determines the reconnaissance element contributed by the time window matching degree . Matching factor refers to the matching degree of the time window, and . refers to the length of the overlapping parts of the two time windows. Under a fixed value of [], is high if the overlapping part is large. That is, point , which is the precursor of , is reasonable from the angle of the time window. Passenger priority factors are considered in the search ant pheromone; that is, the factors are considered in the objective function but are not considered in the reconnaissance element. Moreover, to reduce the scout scope of reconnaissance ants, the statistic result of Quan and Wen [14] shows that the next node of a point is selected in several points nearest to this point in the path optimization solution without calculating all the remaining viable nodes. Thus, the optimization process can be accelerated and the solution is quite reliable, which can converge to the optimal solution with the probability of approaching 1. The reconnaissance result in ascending order is combined with the existing prior knowledge (integrated with max(PC) [14]) to generate another reconnaissance element, which is recorded as , on the path from point to . is represented as follows:where refers to the nearest distance from city as the center to other cities . To synthesize the total reconnaissance element of the above reconnaissance ants,
Considering the need for a certain pheromone on the path at the initial moment, the amount of information of each path at the initial moment is set aswhere and refer to the minimum and maximum distances from point , as the center, to other points . refers to the concentration of pheromone on each path at the initial moment and is a constant. Through the reconnaissance element trail marked by the reconnaissance ants, research ants can make a directional search with the assistance of the pheromone. This can improve efficiency and determine the optimal solution.
The global search of the search ant colony draws lessons from the concept of the AntQ algorithm of Dorigo et al. [13] where the selection strategy is combined with the deterministic and random selection [1, 4]: (1) is the number of search ants, refers to the visibility of side (), and , which reflects the heuristic degree transfer from node to and is different from in the basic ant colony algorithm model. The objective function is placed in the global scope given the best visibility. It reflects the bigpicture thinking of ants and does not necessarily have to select the closest point in a large probability as the next point. However, each step of the ant must consider the objective function of overall optimization. (2) refers to the pheromone trail intensity of side (), refers to state mobility probability of ant from node to node , and refers to the node that is not yet visited. When each ant selects the next node, two points should be considered while selecting the next service point under the premise of adhering to the vehicle capacity and time window constraints. One is visibility and the amount of information access to the next service point. Second is priority conditions priorityofsmall time window and priorityofshortest waiting time.
The search ant colony is tasked to conduct a global search. At each customer point, reconnaissance element and the pheromone at each side select the next service point until they find and mark the best route. The formula for the probability of the search ant colony that transfers from to at during the movement is as follows:where refers to the set of the service points and destinations that ant is allowed to select for the next step; α and β parameters reflect the relative importance the pheromone ants accumulated during the movement and heuristic information provided that the ants select the path, respectively; and are weight coefficients that satisfy , , , and , which refers to the taboo list () each ant creates, where ( refers to the sum of search ant colony) for recording the nodes that ant has visited at and the ants are forbidden to visit again in this cycle. The tabu list is cleared at the end of the cycle.
3.3. Adaptive Information Update Strategy
Traveling ants often cause blockage and stagnation. Ant colony algorithm easily leads to premature and local convergences. To solve this problem and improve the global convergence of ant colony algorithm and search speed, many studies proposed different strategies to update existing information [19–25]; for example, in the standard ant colony algorithm, when information is updated and as long as the ant is traversed, selecting the path can update the information on the path. This enhances the information on the path for optimal fitness and weakens the information on other paths; the other algorithms based on level change enable the ant to set several paths with relatively good fitness, and the degree of excellence of its solution determines the magnitude of information. These algorithms presented above are different. They mainly update the amount of information using the increasing or decreasing proportion of the fixed amount of information and ignore the distribution features of the solution. They improve the characteristics of the ant colony algorithm to a certain extent and are used to treat smallscale problems. This paper proposes a new adaptive information update strategy for solving the largescale problems, which starts from the distribution state of the solution.
In largescale problems, the existence of a volatile coefficient reduces the amount of information on the path that has never been searched to be close to 0, which lowers the search ability of the algorithm on these paths. If the amount of information on a path is large, the amount of information on these paths is increased, and the chance that the path that has been searched is selected again is high. Moreover, the global search ability of the algorithm is influenced. A fixed changing volatile coefficient can improve the global search ability but reduces the convergence speed of the algorithm. Therefore, an adaptive method for changing is proposed to update the pheromone. The formula is as follows [11]:while ,while ,where is a function proportional to the number of convergences . If the number of convergences is higher, is larger, as denoted by the equation below:where is a constant. The distribution of the solution adaptively updates the amount of information to dynamically adjust the intensity of information on each path, causing the moderate concentration or scattering of ants to avoid premature convergence and local convergence and improving global search ability. The adaptive information update strategy adaptively updates the amount of information based on the distribution of the solution. Thus, it dynamically adjusts the intensity of information on each path, increases the diversity of the solution space, improves the global search ability, and prevents premature convergence and local convergence.
4. Case Analysis
The DVRPMC optimal model and algorithm based on ant colony algorithm have many constrained conditions and parameters. To verify the performance of the algorithm, the customer point priority, service time of customer point, coefficient of road condition, vehicle wait cost, and other constrained conditions have not yet been considered in this case analysis. The objective function in the model is also simplified as the path distance. In this case analysis, 15 customer points were distributed in a square region in the edge length of 8 km. Table 1 illustrates the customer points’ coordinates and related demand for passenger capacity, respectively. The region has two scheduling points (parking points). The coordinates are () and (), which own five and eight vehicles, respectively. The maximum passenger capacity of each vehicle is four.

The VRP standard case library of Solomon and Desrosiers [26] is used for detection. The Solomon and Desrosiers case library is divided into three categories: Class C, Class R, and Class RC. The customer points in Class C are distributed in cluster types. The customer points in Class R are distributed at random. The customer points in Class RC are the mixedhalf cluster type. This case analysis is conducted in the hardware environment of Intel Pentium IV 3.0 with 1 G memory. VC language is used for the programming of the algorithm and to calculate R and C. The result is compared with the optimal solution from several existing algorithms. The parameters of the algorithm are α = 1, β = 2, and ρ = 0.8. The iteration number is 50. The results of 10 times the calculation are shown in Table 2.

Table 2 shows that the optimal solutions of RC and R are mostly in the twenties. The calculation time is not more than 16 s in the hardware environment of Intel Pentium IV 3.0 with 1 G memory. Thus, the optimal solution convergence is faster than that of the standard ant colony algorithm. To compare the calculation results, standard genetic, adaptive genetic, and ant colony optimal adaptive algorithms are used to stimulate the case study. The obtained parameters in the case of 10 search times for the Class RC are shown in Table 3.

To compare the calculation results, the average number of vehicles in Table 3 up to two digits after the decimal point is taken. The calculation result shows that, in Class RC, the DVRPMC calculation results of the average number of vehicles and traveling distances are the best. Thus, the calculation result of the adaptive multiconstrained ant colony algorithm is stable with short traveling distance and high calculation efficiency.
5. Conclusion
In the context of neural network analysis, this paper analyzes the existing ant colony algorithm and proposes a multiconstrained network intensive vehicle routing optimal model and algorithm based on the adaptive ant colony algorithm. Drawing lessons from the adaptive polymorphic ant colony algorithm, the adaptive multiconstrained network intensive vehicle ant colony algorithm in this paper is used to collaboratively solve the problems through different types of ants. It overcomes the shortcomings of the traditional ant colony algorithm of having one ant colony with low algorithm efficiency. The ant colony can exchange planning information during the problemsolving process and maintain diversity in the search process because of the interaction among populations. Therefore, through the further use of the adaptive information update strategy, the algorithm of introducing the interaction among populations performs better than the traditional algorithm. Finally, this paper conducts a simulation experiment using specific examples. The experimental result shows that the adaptive multiconstrained ant colony algorithm has better solving efficiency and results than the standard ant colony algorithm. Therefore, for largescale DVRPMC in network intensive vehicle scheduling, the adaptive multiconstrained and colony algorithm has significant superiority. Through the division and collaboration of different kinds of ants, pheromone adaptive strategy is combined with polymorphic ant colony algorithm; the algorithm given in this paper can effectively overcome some disadvantages, such as premature stagnation, and has a theoretical significance to the study of largescale multiconstrained vehicle routing problems in complex traffic network systems. The established general mathematical model can be used in future studies, involving verifying further its result and efficiency in the case study of simulation under different passenger priorities, different road conditions, spatiotemporal traffic networks, and other cases’ study of the optimal routing search method of realtime network intensive vehicle scheduling in a dynamic environment intensively. Furthermore, the algorithm provided in this paper should be significant for the further study on routing problem in the dynamic environment of a spatial and temporal complexity traffic network.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by the Guangdong Natural Science Fund (2015A030313626), Special Fund of Frontier and Key Technology Innovation of Guangdong Province (Provincial Major Science and Technology Project) (2014B010112008, 2015B010106001, 2015B010129003, and 2016B010109007), Pearl River S&T Nova Program of Guangzhou (201610010034), Implementation of Innovation Driven Development Capacity Building Special Funds of Guangdong Academy of Sciences (2017GDASCX0101, 2017GDASCX0601), Guangdong Science and Technology Program (2013B01010201, 2013B010102018, and 2013B090600152), and Guangzhou Science and Technology Project (2014Y200044).
References
 C. Yang, Z. Li, and J. Li, “Trajectory planning and optimized adaptive control for a class of wheeled inverted pendulum vehicle models,” IEEE Transactions on Cybernetics, vol. 43, no. 1, pp. 24–36, 2013. View at: Publisher Site  Google Scholar
 A. GarciaNajera and J. A. Bullinaria, “An improved multiobjective evolutionary algorithm for the vehicle routing problem with time windows,” Computers & Operations Research, vol. 38, no. 1, pp. 287–300, 2010. View at: Publisher Site  Google Scholar
 S. F. Ghannadpour, S. Noori, R. TavakkoliMoghaddam, and K. Ghoseiri, “A multiobjective dynamic vehicle routing problem with fuzzy time windows: model, solution and application,” Applied Soft Computing Journal, vol. 14, no. 1, pp. 504–527, 2014. View at: Publisher Site  Google Scholar
 M. Chen, Z. Zhang, and J. Shi, “Vehicle routing problem with multiple constraints using adaptive and polymorphic ant colony algorithm,” Journal of Southeast University (Natural Science Edition), vol. 38, no. 1, pp. 37–42, 2008. View at: Google Scholar
 C. Yang, Z. Li, R. Cui, and B. Xu, “Neural Networkbased Motion Control of An Underactuated Wheeled Inverted Pendulum Model,” IEEE Transactions on Neural Networks and Learning Systems, vol. 25, no. 11, pp. 2004–2016, 2014. View at: Google Scholar
 H. Wang, B. Li, and K. Liu, “Demand allocation and network flow assignment under emergency rescue circumstance,” Systems EngineeringTheory & Practice, vol. 35, no. 6, pp. 1457–1464, 2015. View at: Google Scholar
 H. Xiao, Z. Li, C. Yang et al., “Robust stabilization of a wheeled mobile robot using model predictive control based on neurodynamics optimization,” IEEE Transactions on Industrial Electronics, vol. 64, no. 1, pp. 505–516, 2017. View at: Publisher Site  Google Scholar
 K. Li, L. Tang, and S. Chen, “Modeling and optimizing for the integrated problem with container storage allocation and truck scheduling,” System EngineeringTheory & Practice, vol. 34, no. 1, pp. 115–121, 2014. View at: Google Scholar
 Y. Chen and J. Wang, “Improved Ant Colony Algorithm for Vehicle Routing Problem with Time Windows,” Journal of Computer Engineering and Applications, vol. 43, no. 29, pp. 218219, 2006. View at: Google Scholar
 L. Shao, “Ant Colony Optimizing Adaptive Gene Algorism Resolving Logistics Vehicle Schedule,” Computer Measurement & Control, vol. 20, no. 5, pp. 1423–1425, 2012. View at: Google Scholar
 Y. Wang and J. Xie, “An adaptive ant colony optimization algorithm and simulation,” Journal of System Simulation, vol. 14, no. 1, pp. 31–33, 2014. View at: Google Scholar
 C. Yang, Y. Jiang, Z. Li, W. He, and C.Y. Su, “Neural control of bimanual robots with guaranteed global stability and motion precision,” IEEE Transactions on Industrial Informatics, vol. 13, no. 3, pp. 1162–1171, 2017. View at: Publisher Site  Google Scholar
 M. Dorigo and M. Luca, “A study of some properties of AntO,” Tech. Rep. TR/IRIDIA/19964, IRIDIA, Universit Libre de BruxeIIes, 1996. View at: Google Scholar
 H. Quan and G. Wen, “Subspacegenetical algorithm for TSP,” Mathematical Theory and Applications, vol. 22, no. 1, pp. 36–39, 2002. View at: Google Scholar
 A. G. Qureshi, E. Taniguchi, and T. Yamada, “An exact solution approach for vehicle routing and scheduling problems with soft time windows,” Transportation Research Part E: Logistics & Transportation Review, vol. 45, no. 6, pp. 960–977, 2009. View at: Publisher Site  Google Scholar
 I. H. Osman, “Metastrategy simulated annealing and tabu search algorithms for the vehicle routing problem,” Annals of Operations Research, vol. 41, no. 4, pp. 421–451, 1993. View at: Publisher Site  Google Scholar
 J. Schulze and T. Fahle, “A parallel algorithm for the vehicle routing problem with time window constraints,” Annals of Operations Research, vol. 86, pp. 585–607, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 A. S. Kenyon and D. P. Morton, “Stochastic vehicle routing with random travel times,” Transportation Science, vol. 37, no. 1, pp. 69–82, 2003. View at: Publisher Site  Google Scholar
 Q. Cao, X. Liu, and X. Ren, “Vehicle schedule problem based on plant growth simulation algorithm,” System EngineeringTheory & Practice, vol. 35, no. 6, pp. 1449–1456, 2015. View at: Google Scholar
 M. Wei and W. Jin, “Discrete particle swarm optimization algorithm for vehicle scheduling problem,” Computer Science, vol. 37, no. 4, pp. 187–191, 2010. View at: Google Scholar
 Z. Zhao, Y. Liu, and F. Luo, “Output feedback boundary control of an axially moving system with input saturation constraint,” ISA Transactions, vol. 68, pp. 22–32, 2017. View at: Publisher Site  Google Scholar
 Z. Zhao, Y. Liu, F. Guo, and Y. Fu, “Vibration control and boundary tension constraint of an axially moving string system,” Nonlinear Dynamics, vol. 89, no. 4, pp. 2431–2440, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Zhao, Y. Liu, and F. Luo, “Boundary Control for a Vibrating String System with Bounded Input,” Asian Journal of Control, In press. View at: Publisher Site  Google Scholar
 Z. Zhao, Y. Liu, W. He, and F. Luo, “Adaptive boundary control of an axially moving belt system with high acceleration/deceleration,” IET Control Theory & Applications, vol. 10, no. 11, pp. 1299–1306, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 Z.J. Zhao, Y. Liu, F. Guo, and Y. Fu, “Modelling and control for a class of axially moving nonuniform system,” International Journal of Systems Science. Principles and Applications of Systems and Integration, vol. 48, no. 4, pp. 849–861, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 M. M. Solomon and J. Desrosiers, “Time window constrained routing and scheduling problems,” Operations Research Society of America. Transportation Science Section. Transportation Science, vol. 22, no. 1, pp. 1–11, 1988. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Shaopei Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.