Urban Traffic Signal Control Based on Multiobjective Joint Optimization

.is paper proposes two algorithms for signal timing optimization of single intersections, namely, microbial genetic algorithm and simulated annealing algorithm. .e basis of the optimization of these two algorithms is the original timing scheme of the SCATS, and the optimized parameters are the average delay of vehicles and the capacity. Experiments verify that these two algorithms are, respectively, improved by 67.47% and 46.88%, based on the original timing scheme.


Introduction
In order to solve the current traffic problem more effectively, this paper proposes a single intersection signal control method. e traffic signal control of a single intersection is the smallest unit of the entire city traffic signal control system. It plays an important role in improving traffic order at intersections, reducing vehicle delays, and improving traffic capacity [1]. With the continuous deepening of the concept of sustainable development of urban transportation, the goal of traffic signal control is gradually changing from a single goal to multiple goals. erefore, the multiobjective joint optimization of traffic signal control has attracted researchers' attention.
When optimizing traffic signal parameters, many factors need to be considered comprehensively, such as average delay, number of stops, and queue length [2]. Among them, the most commonly used criteria are delay and number of stops. ese service quality measures are considered to be the most important factors for users to perceive the service level [3]. Webster established a steady-state random delay model for unsaturated traffic flow [4]. Heydecker and Wey used the average delay of motor vehicle flow at the intersection as the optimization target of signal control parameters [5,6]. Robertson used a linear combination of optimization goals such as average delay and number of stops to establish an optimized model of signal control parameters [7,8].
e multiobjective problem of traffic control has become a hotspot in current research, which has aroused extensive attention of researchers. Zeng et al. established a multiobjective optimization model of signal control parameters. e cell transmission model was used to simulate the operation law of motor vehicle flow at signal control intersections. e average delay and capacity were used as optimization goals [9]. Taking the average delay, average parking rate, and capacity as the optimization goals, Yao et al. established a multiobjective optimization model but did not give a specific solution algorithm [10]. Taking the average delay of motor vehicle flow, average parking rate, and pedestrian waiting time as optimization goals, Ma et al. established a multiobjective optimization model for the period duration of a single intersection with fixed period signal control [11]. Li and Lu took the average delay and average number of vehicle stops at intersections as the optimization goals, established a multiobjective optimization model of signal control parameters for a single intersection, and used genetic algorithms for optimization [12]. Taking the average delay of motor vehicle flow, the  average number of stops, and the total passing traffic as the  optimization goals, Cao and Xu established a weighted  combination optimization model of signal control parameters  for single intersection with saturation constraints and solved the model using genetic algorithm [13].
Although there are many models for multiobjective joint optimization, there are no studies comparing the actual vehicle flow at the intersection with the average delay of vehicles and the capacity of the intersection as variables, compared with the original SCATS. erefore, the optimization algorithm in this paper combines the actual intersection road conditions to analyze the night peak traffic information, which is an innovative point.

Signal Timing Method Based on
Multiobjective Joint Optimization e basic idea of signal distribution is to allocate time resources to traffic lights at intersections according to the traffic flow in each direction, so as to conduct traffic guidance. It is very important for signal synchronization [14]. In general, the best period should be the average period. If the period is too long to significantly improve the traffic capacity, it will lead to further delays of vehicles [15]. According to the traffic distribution, geometry, intersection type, area, and relevant conditions of each intersection, each intersection has a corresponding optimal time period, so the optimal time period changes in different time periods. Traffic function and control intersection signals are the best cycle, which can provide the best length for the whole vehicle [16].
In real life, urban intersections often adopt the timing scheme, which is simple and easy to operate, but it has many limitations.
rough the transformation of the Webster formula, we can get accurate signal timing for real-time traffic and improve the efficiencies of road. is paper analyzes the traffic flow of the same intersection at the same time through two algorithms and then compares these two algorithms with the classic timing scheme of SCATS.

Notes on SCATS.
e SCATS [17], researched and developed by the New South Wales Road Traffic Authority (RTA), is one of the few advanced urban signal traffic control systems in the world. When the system is selected through a real-time solution, the real-time selection of the signal period and green ratio is based on overall demand. e corresponding green light time of the intersection shall be determined in accordance with the principle that the saturation degrees of all phases are equal or close to each other [18]. With the adjustment of the signal period, the green time of each phase also changes.

Basis of Algorithm.
is part will be introduced separately from the signal phase of the intersection, the model of the intersection, and the formula for calculating the objective functions.

Objective Function Calculation Formula.
TRRL method is also known as Webster method or F-B method.
is signal timing method and corresponding calculation formula were first proposed by Corbusier and Webster in the 1990s [19]. e basic idea of this model is to take the vehicle delay time as the only evaluation index when the vehicle passes the intersection and to optimize the signal timing scheme through calculation. According to the theory of signal timing at the intersection [20], the two decisive factors are signal period and green time.
e formula for calculating the objective function PI is as follows: In (1), D is the vehicle delay at the intersection, B is the capacity of the intersection, and α and β are two evaluation indicators, respectively, representing the proportion of vehicle delay and capacity (they are considered equally important in this paper, both assigned the value 0.5).
e formulas for calculating the vehicle delay D and the capacity of the intersection B are as follows: In (2) and (3), D i is the stop delay of each vehicle in phase i, q i is the flow in phase i, and B i is the capacity of the intersection in phase i. e delay at the ith phase intersection is D i , and the traffic capacity at the ith phase intersection is B i . e calculation formulas for the two of them are as follows: In (4) and (5), C is the cycle, λ i is the green ratio in phase i, y i is the traffic intensity in phase i, k i is the number of lanes in phase i, and h i is the lane saturation headway in phase i (the default value for this article is 2.5 s). e formula for period C is as follows: In (6), g i is the decision variable and there are 4 decision variables in total (g A , g B , g C , g D ), g A is the effective green time of phase A, g B is the effective green time of phase B, g C is the effective green time of phase C, and g D is the effective green time of phase D.
at is to say, C � g A + g B + g C + g D + (delay start + delay redlight) × 4; in this paper, when the green light is on, delay-start is set to 3 seconds and delay-red light is set to 2 seconds.

2
Scientific Programming λ i is the green ratio in phase i, and the calculation formula is as follows: In (7), y i is the traffic intensity in phase i, and the calculation formula is as follows: In (8), q i is the traffic flow in phase i. e calculation of the average vehicle delay D 0 and the initial intersection capacity B 0 is similar to that of D i and B i , but there is a difference in calculating the period and the green ratio; C 0 and λ i0 use the calculation method of timing period, and the calculation formula is as follows: In (9) and (10), L is the total delay time of four phases and Y is the sum of the four phases of y i . e above formula is the basis of optimizing traffic signals in this paper. e following are a detailed introduction of the process of each algorithm and the final results from different algorithm perspectives.

Signal Phase and Intersection
Model. Before determining which signal phase is used at the intersection, it is important to fully understand the condition of the vehicle as it rushes from all directions to the intersection.
In this paper, the investigation and research intersection is a three-lane intersection; after taking full account of the actual traffic situation, a four-phase traffic signal is set up. Based on the road near schools, supermarkets, residential areas, vegetable markets, and many other crowded and traffic-intensive areas, we set the signal allocation time separate for the turn (left or right), so as to alleviate the traffic problems that may face the turning vehicles, reducing the traffic pressure at the intersection. e four-phase traffic signal is shown in Figure 1, and the urban intersection model is shown in Figure 2. e numbers 1-12 in Figure 2 are the coil numbers, which are used to calculate the ratio of the actual traffic volume at the intersection to the saturated traffic volume. e directions of the arrows nos. 1-4 are the same as the four phase directions A-D in Figure 1.

Experimental Data Source.
e data source of this article is the real data of one hour (17 : 28-18 : 30) of the night peak of a certain intersection. e set of data includes the effective green light duration of the four phases, the number of vehicles passing through the set green light cycle, and its coil number. Because of the large amount of recorded traffic information, this article randomly selects 20 lists as templates.
In the comparative experiment of this article, there are four sets of timing schemes for the SCATS of the intersection of the test subject. Taking the third set of plans as reference, the green light ratio information of the four sets of plans is shown in Table 1, and the data in Table 2 are records randomly extracted from traffic information.

Experiment Procedure.
Genetic algorithm is a global random search algorithm that simulates the evolutionary process and natural selection process of natural species. e algorithm was proposed by Professor J. Holland of Michigan University in 1975. is method was based on genetic theory and natural selection process [21] and used the natural selection rule and the cross mutation process of chromosomes and genes in biological groups as simulation object [22][23][24]. In this paper, the optimization process of microbial genetic algorithm is as follows: (1) Initialize the values of 4 decision variables randomly (effective green light time for 4 phases).  Figure 3. In the figure, the abscissa represents the number of iterations, and the ordinate represents the value of the objective function at the end. e figure visualizes the entire optimization process of the microbial genetic algorithm. It is clear from the figure that the value of the corresponding objective function changes with the increase of iteration [25]. Before the objective function completes 24 iterations, the value of the objective function fluctuates with the accumulation of the iteration and the displayed state is unstable, but the value of the objective function tends to be stable, approaching 25 iterations [26]. erefore, it is not   Scientific Programming difficult to conclude from the relationship change diagram that, in terms of accuracy and training time, the microbial genetic algorithm performed better. e traffic data is derived from the traffic condition information of an-hour-long evening rush in a certain domestic city. Recording the values of vehicle delay D, traffic capacity B, and objective function PI during this period and comparing them with the SCATS at the same time, we get the results in Table 3.

Experimental Results. Time dependent objective function is shown in
As the data in Table 3 show, we set the same external environmental factors, taking the same vehicle information in the same time period for comparison. e microbial genetic algorithm has improved by 67.465% in optimizing the objective function [27]. Obviously, the microbial genetic algorithm proposed in this paper is superior to SCATS in traffic signal optimization. e SCATS is the basis of the algorithm optimization in this paper, with a period of 230 s. e time given in Table 4 is the effective green light time when the microbial genetic algorithm obtains the best solution.

Algorithm Research History and Main Announcement.
American physicist N. Metropolis and colleagues published an article in 1953 that studied complex systems and calculated the energy distribution in them. ey used Monte Carlo simulation to calculate the energy distribution of molecules in multimolecular systems [28]. is is a term often referred to in the simulated annealing algorithm as the Metropolis criterion.
In 1983, physicists S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi of the IBM Corporation in the United States borrowed the method of Metropolis and others to explore a rotating glass state system [29]; they discovered the energy of the physical system and the cost of some optimal combination problems. e function is quite similar: seeking the lowest cost is like seeking the lowest energy. As a result, they developed an algorithm based on the metropolitan process to find the best solution. e optimal solution (the maximum value) of the problem is similar to the minimum energy of the system; when the system cools down, the energy also gradually decreases, and in the same sense, the solution to the problem falls down to its maximum. e calculation formula for the Metropolis criterion is as follows: In (11) and (12), p is the probability of cooling, K is the constant, and exp is the natural index. e cooling calculation formula of the simulated annealing algorithm is as follows:     Scientific Programming In (13), T 0 is the initial temperature.
In the process, the maximum tempering temperature gradually decreases with each annealing process, and the algorithm operation is stopped until the tempering end temperature is reached. e tempering process during the inner loop of the algorithm increases the algorithm's ability to accept nonoptimal solutions, which helps the algorithm to jump out of the local optimal solution and obtain the global optimal solution.

Experimental Results.
In this simulation algorithm, the relationship between the temperature and the iterative algorithm in the driving process, and between the objective function and the external iterative algorithm is as follows. Figure 4 shows the linear change in temperature as the number of iterations increases. Meanwhile, in another aspect, what is presented in Figure 5 is the change corresponding to the objective function when the number of iterations is accumulated. During the accumulation of the number of iterations, the temperature decreases accordingly. When the temperature reaches 0.4°C, the functional relationship is stable and convergent. In the graph of the relationship between the objective function and the number of outer iterations, the abscissa is each outer loop, and the ordinate is the objective function value. It can be seen that with the decrease of iteration depth and temperature, the value of the objective function gradually becomes smaller [30]. e optimization effect of the simulated annealing algorithm is getting better and better. Meanwhile, it can be seen from the figure that, at the 500th cycle, the SA model has approached convergence and the target value is around 7.7.
Similar to the microbial genetic algorithm, the values of vehicle delay D, traffic capacity B, and objective function PI are also recorded. e recorded data is shown in Table 5. Compared with the SCATS timing scheme, the objective function is improved by 46.848%.
By setting the same variables and taking the same dataset, the simulated annealing algorithm has been partially improved compared to the timing of the SCATS. e time given in Table 6 is the effective green light time when the simulated annealing algorithm obtains the best solution.

Conclusions
is paper uses two algorithms to optimize the signal timing scheme of intersections. e innovation of this paper lies in the use of mutation genetic algorithm and simulated annealing algorithm to optimize the timing scheme of SCATS. rough the real traffic flow data, the effectiveness of the two optimization algorithms in improving the objective function is demonstrated. e advantage of the microbial genetic algorithm is that it will select the one with poor fitness from the two individuals as the mutation object; therefore, the quality of the optimal solution can be improved. At the same time, the microbial genetic algorithm is not highly dependent on the initial population, and two individuals can be randomly selected to calculate their objective function value and individual fitness. e global search ability of the algorithm is robust, and it will not fall into the trap of local optimal solution. e advantage of the simulated annealing algorithm is that its entire operation process is very simple. e algorithm has strong robustness and applicability. Simultaneously, in order to ensure that the global optimal solution is not missed during the solution process, a higher initial temperature is set.
Both algorithms require a lot of iterations, so there are drawbacks in solving real-time problems. On the one hand, the signal optimization in this paper is only for a specific   intersection, which lacks general applicability. On the other hand, when selecting the objective function, other traffic factors are not fully considered, and there are certain errors here. Future research on signal control of intersections may consider more parameters, and the vehicle condition information of adjacent intersections will also affect that of the others. Incorporating more variables into the reference model, extending the research on individual intersections to intersection groups, and comparing the original SCATS timing scheme with actual traffic flow data can be the directions of hard research in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest regarding the publication of this paper.