Bilevel Programming for Traffic Signal Coordinated Control considering Pedestrian Crossing

,


Introduction
Taking Nanjing as an example, in recent years, the growth of urban motor vehicles has been obvious.
e main artery roads of many cities are in near-saturated state all the year round.
is has brought great negative impacts on our traveling, which not only affects the efficiency of work, but also causes people emotional dissatisfaction, so it is an important mission of our traffic workers to improve traffic capacity and maximize traffic efficiency within a limited road space [1,2].
rough the coordinated control of traffic signals, we can meet people's traffic needs to the maximum extent while ensuring safety [3][4][5]. e function of an isolated intersection is mainly determined by its control parameters, and of course by other intersections. However, due to the complexity of modern transportation systems, coordinated control of traffic signals has become one of the important traffic management strategies in cities [6,7]. Mapping to reality is even more unimpeded traffic flow, fewer traffic accidents, and more stable driving speed [8].
erefore, it is particularly important to coordinate the signals at the intersections around the subway stations where congestion is increasing and traffic demand is in expansion [9]. e project will use microsimulation as a research tool to analyze the main factors that need to be considered in the coordinated control of traffic signals. e time-space conflict between the vehicle and its influence on the coordinated control of the signals proposes a traffic signal control method that is conducive to subway passengers' travel [10,11].
Here is the literature review of signal coordinated control. In terms of coordinated control of signals, as early as the 1970s, there was a MAXBAND model for coordinated control of artery roads [12]. e concept of a progression band was proposed. e two-way progression bandwidth was taken as the objective function, and the maximum objective function model was used, combining mixed integer linear programming. en there is the MULTIBAND model [13,14], which is based on different traffic requirements on different road sections, and the corresponding bandwidth is calculated to get the best coordination effect. For the regional coordinated control, an off-line optimized network signal timing of TRANSYT proposed by the British Road Traffic Research Institute has been developed to comprehensively consider the comprehensive model of cruise time, bus modeling, peak delay calculation, and fuel consumption estimation [15], as well as inductive regional traffic coordination control systems SCATS, SCOOT, and others, which can effectively improve the traffic conditions of unsaturated and supersaturated road networks [16].
After combining the bus priority strategy, Vasudevan [17], under the framework of bus priority coordination control, with the minimum green time as the constraint and the minimum per capita delay under coordinated control as the optimization goal, established the optimization of trunk coordination control based on bus priority model. Many advanced signal control systems take into account different priority transit strategies and are added to the algorithm to better and more efficiently adapt to multimode control needs. Conversely, some systems are designed to reduce per capita delays, such as the reduction of motor vehicle delays. Some studies [18,19] have added traffic priority based on traffic rules, and, for Lo [20], a series of motor vehicles are selected in advance to implement their priority control methods. In addition, some systems weigh the factors such as pedestrian delays, motor vehicle delays, bus delays [21], or by reducing the transit time of buses and the waiting time of pedestrians on downstream stations and other priority strategies for other transportation systems [22].
Recently, signal coordination control has increasingly favored the use of vehicle-coordinated data, and they have attempted to combine all of the above factors, including weighing the delays of various passers-by, taking into account various priority strategies. In the system operation, real-time traffic information needs to be obtained. Tan realized vehicle information interaction through AVL (Automated Vehicle Location) or APC (Automated Passenger Counter) [23]. However, such systems are limited to two connected main road intersections (not more), and there are more problems with multiple bus lines and multiple priority control strategies [11]. ere are still researches using machine learning methods to deal with the problems within the intersection between vehicles and pedestrians [24][25][26]; however, the results are limited and unexplainable for traffic design and control [27,28].
Although there are some researches focusing on the modeling of pedestrian crossing behavior [29][30][31], to the best of our knowledge, it is rare to combine models of coordinating signal control and pedestrian delay. So, we consider the problems of pedestrian crossing with signal control. In this paper, the situation of green light time and relative offset between two adjacent intersections is analyzed comprehensively. en, according to the average travel time of the motor vehicle, the green wave feasibility analysis is carried out to calculate the corresponding progression bandwidth and the delay of the motor vehicle flowing under the green wave (not considering the delay of the nongreen wave motor vehicle). Finally, the complete intersection green wave analysis is obtained. en, the model of pedestrian delay is analyzed. By collecting a large amount of data, the model of a different motor vehicle and pedestrian flow delays is obtained. en, the vehicle flow and pedestrian flow are simulated by Vissim to test the simulation of each model under different conditions. e degree of integration is obtained to generate a more accurate calculation model of pedestrian delay. Secondly, the algorithm is written in Python, and the above two models are merged together and applied to the actual intersection coordination model, and the optimal coordination scheme is obtained by inputting into the algorithm. Finally, the Vissim is used to simulate the actual situation to judge the rationality of the established model and the edited algorithm and draw conclusions. e remaining part of this paper is as follows: in section 2, we put forward models of progression bandwidth under different situations of cycle length, green time, and offset. Section 3 illustrates the process of establishing models of pedestrian delay. Section 4 proposes the bilevel programming fusing the aforementioned two models and an algorithm to solve it. Furthermore, it also provides a case study to test the reasonableness of them. Finally, we conclude in Section 5.

Model of Progression Bandwidth
Due to the high land utilization rate in large cities, the distance between adjacent intersections is very close. erefore, it is necessary to coordinate the signal scheme to maximize the traffic efficiency of the vehicle. is is called artery traffic coordination control. e so-called green wave means that the vehicle (vehicle flow) continuously passes through a certain main road and encounters less red light as much as possible.
ere are many parameters in the signal coordination control, including signal period, green signal ratio, offset, and so on. e main idea of building a model is obtained [32], which is modeled by drawing a spatiotemporal map.
is paper mainly studies the coordinated control of two adjacent intersections and calculates the delay and progression bandwidth generated by the corresponding green wave traffic by considering the situation of the vehicle queuing at the second intersection. is paper is based on and beyond the content of the former, and at the same time, corrects some errors. rough the comparison of the red and green light phase difference between the two intersections and the first intersection lipstick and green light time, the complete signal set of adjacent intersections is summarized. en, the subclass is divided by the average transit time of the motor vehicle, and then the coordination parameters of each subclass are analyzed in turn, and finally, the complete model is obtained.
Because the coordinated control is affected by multiple factors, the vehicle discretization is serious, resulting in the formation of a fleet; there are a large number of vehicles entering and exiting the road, thus destroying the formation of continuous traffic. e signal lights at the intersection are complex, with more than two phases and coordination difficulties. erefore, before establishing a model, it is necessary to make some reasonable assumptions, both not to deviate from the actual situation and to meet the conditions for establishing the model.
(1) e effects of road traffic conditions, such as weather and road geometry, are not considered.  Table 1, and the relationship of some of them is given as follows (specific symbols in models will be given correspondingly): In Figure 1, when Δ r i,i+1 > r i , Δ g i,i+1 < g i , according to the position of the last intersection corresponding to the start time of the red light at the next intersection, it is divided into three parts of four categories. e slope of red lines in Figure 2 is the basis of classification.
Category I (t i,i+1 < Δ r i,i+1 − r i (Figure 1 part ①)): in Figure 2(a), at this moment, some vehicles departing from the green light at the intersection i can pass through the i + 1th intersection, and some of them have to wait for congestion. t 2 denotes the time interval when vehicles meet with congestion at the specified speed, and t 3 denotes the congestion evacuation time (for the evacuation through the parking line, the same below). D denotes traffic delay of vehicles from the progression band, which is represented by the shadowed part in this figure. So we can derive the following equations: Category II (Δ r i,i+1 − r i < t i,i+1 < Δ r i,i+1 (Figure 1 part ②)): in Figure 2(b), the vehicles from ith intersection at red time and the link connecting ith and i + 1th intersection generate queue. So when the progression vehicles come to Offset of red starting time between I i and I i+1 (s) Δ Offset of green starting time between I i and I i+1 (s) q r i Traffic volume at I i+1 when I i is red (pcu/s) q g i Traffic volume at I i+1 when I i is green (pcu/s) S i Saturated flow rate at I i (pcu/s) Average travel time between I i and I i+1 (s) b i Progression bandwidth (s) the i + 1th intersection, the front part of them should wait. We denote t 1 as the max time interval when vehicles leave from ith intersection and wait at the next node. t 2 denotes the waiting time of the foremost vehicle from the progression band, and t 3 denotes the congestion evacuation time. D also denotes the traffic delay of vehicles from the progression band, which is represented by the shadow part in this figure. We then derive the following:  Journal of Advanced Transportation Figure 1 part ③)): in Figure 2(c), at this time, the vehicle which reaches at the start time of the next intersection at the specific speed is at the range of green time at ith intersection, and the tail of the traffic flow in the progression band must be congested, and the front of the vehicle should wait or not wait depending on evacuation rate at i + 1th node. At this time, it is divided into two categories. is one discusses that the saturated flow rate is large enough that the congestion has vanished when the head of the progression band passes. We denote t 1 as the max time interval when vehicles leave from ith intersection and wait at the next node. t 2 denotes the waiting time of the end vehicle from the progression band, and t 3 denotes the congestion evacuation time, t′ denotes threshold time between block and unblock. D also denotes the traffic delay of vehicles from the progression band, which is represented by the shadow part in this figure. We then derive the following: In Figure 2(d), when the saturated flow rate is not enough, the head of the progression band from ith intersection should wait. t 1 denotes the green time interval of the tail of traffic flow at ith intersection, t 2 denotes the green time interval of the head of traffic flow at ith intersection, t 3 denotes the congestion evacuation time, and t′ denotes threshold time between block and unblock.
In Figure 3, when Δ r i,i+1 < r i , Δ g i,i+1 < g i , according to the position of the last intersection corresponding to the start time of the red light at the next intersection, it is divided into three parts of four categories. e slope of red lines in Figure 4 is the basis of classification. (Figure 3 part ①)): in Figure 4(a), this situation is the same as category II in the interlaced one ( Figure 2(b)). So the notions are omitted, and equation (3) should be used. Figure 3 part ②)): in Figures 4(b) and 4(c), this situation is the same as category III in the interlaced one (Figures 2(c) and 2(d)), So the notions are omitted, and equation (4) and (5) should be used. Figure 3 part ③)): in Figure 4(d), under the specific speed, the vehicles from the ith intersection of the progression band must wait at the next intersection, so it can generate green wave, which means that this situation is eliminated.
In Figure 5, when Δ r i,i+1 > r i , Δ g i,i+1 > g i , according to the vehicle's position at the last node which reaches at the start Journal of Advanced Transportation time of red at the next intersection, it is divided into three parts of five categories. e slope of red lines in Figure 6 is the basis of classification.
Category Figure 5 part ①)): in Figure 6(a), at this time, the vehicle which reaches at the start time of the next intersection at the specific speed is at the range of green time at ith intersection, and the tail of the traffic flow in the progression band must be congested, and the front of the vehicle should or should not wait depending on evacuation rate at i + 1th node.
At this time, it is divided into two categories. is one discusses that the saturated flow rate is large enough that the congestion has vanished when the head of the progression band passes. We denote t 1 as the max time interval when vehicles leave from ith intersection and wait at the next node. t 2 denotes the waiting time of the end vehicle from progression band, t 3 denotes the congestion evacuation time, and t′ denotes threshold time between block and unblock. D also denotes the traffic delay of vehicles from the progression band,  Figure 3: Classification of forward overlapping spatiotemporal figures according to offset. 6 Journal of Advanced Transportation   Journal of Advanced Transportation 7 which is represented by the shadow part in this figure.
We then derive the following: In Figure 6(b), when the saturated flow rate is not enough, the head of the progression band from ith intersection should wait. t 1 denotes the green time interval of the tail of traffic flow at ith intersection, t 2 denotes the green time interval of the head of traffic flow at ith intersection, t 3 denotes the congestion evacuation time, and t′ denotes threshold time between block and unblock.
Category Figure 5 part ②)): in Figure 6(c), this situation is the same as category II in the interlaced one (Figure 2(b)). So the notions are omitted, and equation (3) should be used. Figure 5 part ③)): in Figures 6(d) and 6(e), this situation is the same as the category III in the interlaced one (Figures 2(c) and 2(d)), So the notions are omitted, and equations (4) and (5) should be used.
Finally, we have proposed progression bandwidth models considering different offset, green time, and average travel time, which is summarized in Table 2.
Above all, we get the progression bandwidth and motor vehicle delay of adjacent intersections in each case, and finally, there will be multiple progression bandwidths on different road sections. At present, the optimization methods for coordinated control of traffic signals are the maximum progression bandwidth method and the minimum delay method. e maximum progression bandwidth method is to ensure that the vehicle (vehicle flow) can occupy the longest period of one cycle through the intersection in a continuous manner. Since the cycle length of each coordinated intersection is assumed to be the same, the progression band is satisfied here. Of course, this will ensure that the most traffic flow can participate in coordinated control. For the case where the length of the progression band is different in different sections, the minimum value is taken as the inspection index. e minimum delay method is to minimize the sum of the delays of the progression motor vehicle at the intersection (the delay of the nongreen wave vehicle is not considered here), which is another reasonable solution in theory. Delays are not only related to coordinated control but also related to vehicle emissions, parking, fuel consumption, and so forth, thus affecting the corresponding driving costs.
In addition, in the equations of the aforementioned models, some models' bandwidths are negatively correlated with delays, and some models satisfy the quadratic function of the two variables, which is relatively complicated. When the progression bandwidth is the largest, the motor vehicle delay is not necessarily the smallest; otherwise, when the motor vehicle delay is the smallest, the progression bandwidth is not necessarily the largest. In addition to the pedestrian delays to be studied after this article, there are at least three evaluation indicators, as for how to assign weights, how to give comprehensive evaluation indicators reasonably, and then give them in detail in Section 4.

Model of Pedestrian Delay
Since this study not only considers the coordinated control of traffic signals, but also considers a large number of pedestrian flows, it is necessary to consider not only the coordination and optimization of the progression bands at each intersection, but also the minimization of pedestrian crossing delays. In fact, signal control interacts with pedestrians crossing the street. Pedestrians will affect the pedestrian clearing time at the intersection, thus affecting the green time of the opposite vehicle. Too many vehicles will 8 Journal of Advanced Transportation aggravate their conflict with pedestrians, resulting in increasing delays. In addition, there have been many studies on the delays of pedestrians at the intersection of signal timing, so we should consider the degree of adaptation of different models in different specific situations. e specific process of this chapter is as follows: firstly, we introduce several pedestrian delay models, analyze their advantages and disadvantages, applicable conditions, and then use Vissim microsimulation to obtain three cases of motor vehicles through cluster analysis, simulating pedestrian flow from less to more. en the simulation results are compared with the calculation results, and the errors are calculated to determine the adaptability of each model under different flow rates, which lays a theoretical foundation for the actual analysis of pedestrian flow in the following section.

Different Models of Pedestrian Delay
where C denotes cycle length of intersection(s), G denotes green time of a cycle(s). e model is based on the equal pedestrian arrival rate, the same signal period, and there is no conflict between pedestrians and motor vehicles. It is an ideal model which is applicable to intersections where motor vehicles and pedestrian traffic are small, which is not suitable for the background of this project. Nevertheless, a series of improved models based on this model have been promoted even more.

BR Model [33]
where F denotes the proportion of pedestrians obeying traffic rules. is model takes into account pedestrians who do not follow traffic rules and believes they will not be delayed. Especially suitable for the situation in developing areas, people generally have a weak sense of traffic. However, the model does not consider the conflict between pedestrians and motor vehicles, and there is no clear behavioral assumption. e applicable conditions are subject to great constraints. [34].

LI Model
where d G denotes average delays of pedestrians(s), k NU denotes adjustment factor of uneven arrival rate, R E denotes effective red time(s), P w denotes the proportion of violations, and q denotes average arrival rate of vehicles (pcu/s). is model assumes that the pedestrian arrival rate is not fixed throughout a cycle and that pedestrians arriving in the green time and pedestrians who do not obey the traffic rules will still have delays. Furthermore, the average delay of pedestrians varies with the entire cycle, and these changes are certain. At the green time, the average delay of pedestrians is fixed. At the beginning of the red time, the value is the largest, then linearly decreases, and periodically changes. is model is not suitable for areas where pedestrian flow is high, and the calculated value will be small. [35].

Z Model
where s denotes saturated pedestrian flow (ped/h), q denotes average arrival rate of pedestrians (ped/h), λ 1 denotes arrival Table 2: Summary of progression bandwidth models.

Conditions
Analysis position (ith node) Congestion situation Equation Green-green Red-green Green-red-green Mid Head Tail Head and tail Red-green Green-red-green Red-green-red Head Tail Head and tail All flow Green-red-green Red-green Green-red-green Tail Head and tail Head Tail Head and tail rate of left-turn vehicles (pcu/s), λ 2 denotes arrival rate of right-turn vehicles (pcu/s), τ denotes safety gap(s), W denotes width of motorway (m), V p denotes average speed of pedestrians (m/s), R denotes reaction time(s), and l denotes time through vehicle(s). is model considers the delays of pedestrians waiting at the red light and also takes the delays caused by left-or right-turning vehicles as they cross the crosswalk into account. [36].

MV Model
where α 1 denotes adjustment factor of uneven arrival rate, α 2 denotes pedestrian ratio of noncompliance with traffic regulations, Δ t denotes ideal time of crossing(s), c denotes adjustment factor of pedestrian crossing, D VIT denotes conflict delay of pedestrian and vehicle(s). e model actually has three components, including pedestrian waiting time delay, transit time delay, and conflict delay with the motor vehicle. e pedestrian waiting time delay is obtained by correcting HCM delay formula. e transit time delay is obtained by the speed expansion theory. e third part is obtained through a binary logit model, which was obtained by directly looking up the table during the project. e model is suitable for complex intersections and assumes that the pedestrian arrival rate is uneven; pedestrians will violate traffic rules, pedestrian speed changes, and pedestrian-vehicle conflicts.

Adaptive Analysis.
We have obtained five different models from the literature review, and we have two variables: vehicles flow and pedestrian flow. To make it simple, we discrete them, that is, take a value every 50 and 100, respectively, in Vissim simulation. However, we also need a dimensionality reduction. In this study, we use cluster analysis. For more information, please see the appendix. Here we directly show the results. e relationship between pedestrian delay and motor vehicle flow is not significant, while significant in pedestrian flow. According to the flow rate of the motor vehicle, it is divided into three categories, that is, q � 600 pcu/h when the motor flow is small; q � 1500 pcu/h when the motor flow is moderate; q � 2400 pcu/h when the motor flow is large.
Here are the details in Figure 7. ere are six subplots, in the left column are the pedestrian delay(s) as the function of pedestrian flow (ped/h), in the right column are the relative error as the function of pedestrian flow (ped/h). In the first row, the two figures depict the situation under small vehicle flow (q � 600 pcu/h). In the second row, the two figures depict the situation under moderate vehicle flow (q � 1500 pcu/h). In the third row, the two figures depict the situation under large vehicle flow (q � 2400 pcu/h). e gray (only in the left column) lines denote delay in the simulation, the yellow lines denote delay and relative error of HCM model, the red lines denote delay and relative error of BR model, the green lines denote delay and relative error of LI model, the blue lines denote delay and relative error of Z model, and the purple lines denote delay and relative error of MV model. In Figures 7(a) and 7(b), we find that when the arrival rate of pedestrians is less than 600 ped/h; HCM model is more in line with this situation. When the arrival rate of pedestrians is between 600 ped/h and 1200 ped/h, LI model is more accurate because as the pedestrian flow increases, there are more people who do not obey the rules, and these people can also generate delays. When the flow rate exceeds 1200 ped/h, the MV model fits well, because when pedestrians increase, the conflict with the motor vehicle increases gradually (although the motor vehicle flow is small). At this time, both the MV model and the Z model consider the conflict between pedestrians and motor vehicles. However, since the Z model is limited by the maximum flow of pedestrians, the error near the maximum flow is significantly large, and the conflict between pedestrians and motor vehicles in the MV model should be obtained from many curves of surveys, so the fault tolerance rate is higher.
In Figures 7(c) and 7(d), when the pedestrian flow is less than 1200 ped/h, the LI model has a good fitting degree. When the pedestrian flow is greater than 1200 ped/h, the MV model has a high degree of fitting. While other models are not suitable under these conditions, especially unlike before, the HCM model has a low fit at low flow rates. is is because the delay model just only considers cycle length and green time as the vehicle flow increases. ere are also a number of correction factors to consider, including pedestrian behavioral factors, delays caused by pedestrians and motor vehicle conflicts. ese two models take the above effects into account.
In Figures 7(e) and 7(f ), when the pedestrian flow is less than 900 ped/h, the Z model has a higher degree of fitting. When the pedestrian flow is between 900-1600 ped/h, the LI model fits to a higher degree. When the pedestrian flow is greater than 1600 ped/h, only the MV model can be used. is is because when the pedestrian flow is small, the MV model overconsiders the delay caused by the conflict between the motor vehicle and the pedestrians, so the error is relatively large; on the contrary, when the pedestrian flow is large, the collision between the pedestrian and the car is gradually intensified. is should be taken into consideration, and the Z model is nonlinearly increasing due to the limitation of the pedestrian flow, and the growth rate is gradually increasing. erefore, the Z model is no longer applicable when approaching the upper limit of pedestrian flow. At the same time, the LI model takes care of both. 10 Journal of Advanced Transportation Finally, we summarize the equations of calculating pedestrian delay as follows:

Algorithm and Case Study
From Sections 2 and 3, we have established coordinated control progression band models and pedestrian delay models according to different actual situations. e most crucial point of this project is to combine the two models. For example, in the vicinity of the subway station, it is necessary to consider not only the optimization of coordinated control, but also the large number of pedestrian flows periodically entering and leaving the subway station. Firstly, we propose an algorithm to deal with this bilevel problem and then use a case to test the rationality of models and this algorithm.

Bilevel Programming and Algorithm.
From the above analysis, the problem can be transformed into the following: Journal of Advanced Transportation min f C, g 1 , g 2 , . . . , g n , q s.t. G 1 C, g 1 , g 2 , . . . , g n , q ⊂ U, max f 1 C, g 1 , g 2 , . . . , g n , p min f 2 C, g 1 , g 2 , . . . , g n , p Since the study focuses on the delay of pedestrian crossing, the minimization of pedestrian delay is defined as the objective function of the upper planning, where q is the other variable related to pedestrian delay except for the cycle length and the green time, which is represented by a letter and it represents a difference from the lower level planning decision variables; the constraints can be expressed as a general formula, specifically the constraints in the actual problem, including the green light time being less than the signal cycle length and within a certain range such as (C/4, 3C/4). e lower level planning is a multiobjective nonlinear programming. One is that the progression bandwidth is the largest, and the other is that the total delay of the vehicle at the intersection is the smallest. As with the upper plan, a parameter p is used here to represent the decision variables specific to the underlying plan. e constraints of the lower-level planning are more complicated, including not only the limitation of the green time, but also the constraints of various situations in the model.
It is apparent from the equation that correlating the models of the upper and lower layers by the green time and cycle length has several advantages in the analysis of the problem. ere are conflicts and dependencies between the upper and lower layers. Since the road right in the road is certain, the more travel time is given to the pedestrian, the corresponding motor vehicle delay will increase, and the green wave will be difficult to form. Conversely, the smoother the vehicle travels, the less time it takes for pedestrians to cross the street. However, the decision variables in their objective function are repeated and therefore dependent on each other. e decision of the lower layer not only affects the value of its own objective function, but also affects the decision of the upper layer. However, it does not completely unconditionally conform to the upper level plan, and it retains certain autonomy.
ere are many bilevel programming solutions, including fuzzy algorithm, branch and bound method, KT method, and pole algorithm. However, the settings are problem-specific under this scenario, like the cycle length is mainly between 60 and 180's as an integer and the green time is positive.If the range of the specific optimal value is not known in advance, it cannot be obtained by branch and bound. And the n-dimensional model can be established only counting the green time at n intersections, and then the offset of adjacent intersections is calculated. e final model must exceed n-dimensional, so the general gradient algorithm and iterative step size algorithm are very difficult to execute. Considering the particularity of the problem, this study turns the double goal into a single objective function as follows: Comparing the five objective functions [37] of traffic signals and cost-effect function proposed by Dr. Peñabaena-Nibbles [38], we put forward this final objective function, where i denotes ith intersection, j denotes direction of entrance of intersection, b ij denotes width of progression band, d v ij denotes delays of vehicles, d p ij denotes delays of pedestrians, and q ij denotes flow of vehicles. e method needs to list the combination of different green times in the same cycle, calculating the corresponding delay and progression bandwidth and the specific computational structure of each layer is shown in Figure 8. Progression bandwidth and motor vehicle delays require three columns as subsystems for the green time and offset between the two intersections and each intersection requires only green time (and other factors which do not change with cycle length) when calculating pedestrian delays. So, each green light column is listed as a subsystem. We then think of it as solving the shortest path problem in the calculation, using the Dijkstra algorithm here. Each subsystem is regarded as the path in the shortest path problem, and the calculated value of the subsystem is the weight. erefore, the two layers of the problem can be regarded as solving the shortest path problem with n nodes. e pseudocode of the algorithm is in Algorithm 1.
In Algorithm 1, s, i, j denotes node (intersection) and d denotes the total distance of each node. We then suppose an origin and a destination, finding the shortest path between these two nodes. And the corresponding cycle length and offset between the adjacent intersection can be calculated, which means the optimization scheme will be obtained.

Case Study.
is study analyzes the road between Jiming Temple, Fuqiao, and Daxinggong Station of Nanjing Metro Line 3, as shown in Figure 9, namely, Taiping North Road (North-South), Beijing East Road (East-West) from the north and Zhongshan East Road to the south. (Eastwest), with a total length of approximately 1.9 km, with several intersections, including eight signal-controlled intersections. e section is located in the Xuanwu Lake District of Nanjing and is one of the north-south traffic skeleton roads. Judging from the spatial structure of the whole city, Taiping North Road (Beijing East Road-Zhongshan East Road) is an important channel for communication between the north and the south of Xuanwu District, which mitigates the main flow of three subway stations to various functional areas; from the surrounding areas in terms of function, the area of Taiping North Road is an important urban block integrating office, residence, study, shopping, and other functions. ere are Southeast University Sipailou Campus, Huahai 3C Plaza, Sipai Building Business Center and Presidential Palace and major learning, shopping, and office spaces such as libraries. erefore, the traffic design of this main road is not only related to the efficiency of the entire urban transportation system, but also closely related to the work and life of the surrounding residents. It has important practical significance.
In addition, there are many places for students to study and conduct on this main road, including Southeast University and Chengxian Street Primary School. e walking demand is very large, which is in line with the background of pedestrians crossing the street. Furthermore, many intersections, such as Shipopo-Taiping North Road, Wendui Bridge-Taiping North Road, and Sipai Building-Taiping North Road are T-shaped intersections, and Beijing East Road-Taiping North Road is only eastbound. In addition, some north-south roads are forbidden to turn left. All the conditions are similar to our assumption of coordinated control, that is, " e road between intersections does not contain a large number of vehicles that are transferred in or out" is similar. erefore, this study selects this main road as the simulation object.
We collected corresponding data (flow, speed) by field investigation. en according to the aforementioned algorithm, we can finally get the results of traffic signal coordinated control in Table 3. Here we should note that the cycle length among the eight intersections is identical, that is, 120 s, which is one of the assumptions in this study. en we   Table 4. We get classification of vehicle and pedestrian service levels to evaluate the traffic situation.
It can be seen that in addition to Zhongshan Road-Taiping North Road and Zhujiang Road-Taiping North Road, the delays of other coordinated intersections have been improved, and some service levels have even improved by one level. e average queue length is 4 to 8 vehicles in accordance with the average headway distance of 6 m, which is acceptable. However, motor vehicle delays at the two complex intersections increased after the improvement, with an average queue of more than 15 vehicles. In terms of pedestrians, the average delay of the 8 intersections is within 30 s, the service level was more in the B and C sections, and the effect on pedestrian control reached expectations.
In short, although the delay of the two separate intersections increased slightly, the control effect of the motor vehicles and pedestrians on the whole road was improved due to the formation of green waves, which basically achieves the purpose of the project, as well as motor vehicle progression band model, pedestrian delay model, and the rationality of the algorithm. To further improve the result, we will analysis and connect findings between branch-and-bound search [39] and algorithm in this paper.

Conclusion and Discussions
In this paper, a progression band model is established for the progression band and delay generated by the motor vehicle passing through adjacent intersections. On the basis of the existing model, the conditions are further refined, so that the model covers different intersection cycle length and green time, and the calculation formulas of progression bandwidth and delay under different conditions are given. In order to obtain a more accurate model for calculating pedestrian crossing delays, different models and formulas established in several documents were integrated, and their advantages and disadvantages and applicable conditions are analyzed. Furthermore, three kinds of motor vehicle flow speed were established with Vissim. In this case, different pedestrian flows are set from small to large, and the simulation delay is compared with the calculation delay. Finally, we get the  In order to fuse the above two models, a bilevel programming problem is established. en, we analyze the characteristic of this study to simplify the objective function and use the shortest path searching to get the optimization results. In order to verify the rationality of the algorithm and the corresponding models, we find artery roads in line with the background of the project, conduct research and analysis, and obtain relevant data. After integration, on the one hand, the numerical classification will be carried into the algorithm to obtain the results of coordinated intersection, that is, cycle length, green time and offset, and the progression bandwidth based on the above models; on the other hand, the results are input into Vissim for simulation, in an environment with certain behavioral assumptions. e operation is carried out  to obtain the effect and evaluation index after the improvement under the actual situation. Finally, the service level analysis is carried out, and the certain rationality of the project model and algorithm can be guaranteed. e biggest contribution of this paper is to put forward models of progression bandwidth and pedestrian delay and establish a correlation algorithm with intersections as subsystems to solve the corresponding signal coordination parameters.
Here is the future direction of research in my mind. In the establishment of the progression bandwidth model, it is assumed that the upstream vehicles are input with continuous and uniform green wave velocity in the green time. When multiple intersections are considered simultaneously, the situation becomes more and more complicated. For a driver there are two choices at the intersection, waiting or passing. erefore, for a main road with n intersections, there are 2 n− 1 possibilities for a car through the main road. In this paper, only the progression band of two adjacent intersections is simply analyzed, which is deviated from the actual situation. However, direct analysis is very complicated. erefore, in the progression band model, how to efficiently analyze the green wave between multiple intersections becomes a research direction that needs in-depth study.

Cluster Analysis of the Relationship among Variables
To get the dimensionality reduction and better know the relationship between pedestrian flow and two forms of flow, we use K-means clustering. Firstly, k points of clustering number are randomly selected as the starting center, then the remaining points and the points closest to the starting center are classified into one category, and the center of gravity of each category is calculated, then the center of gravity is used as the new center to continue iteration until the convergence conditions are met. In addition, this method needs a given number of clusters and is sensitive to abnormal data and data noise.
System clustering regards each object as an independent class, merges the two closest samples to form a new point, and stops when all the data are combined into a cluster or reaches some convergence rule. We simulate the relationship in Vissim and get the results in Figure 10.
From the figure, we can preliminarily find that the horizontal variation of delay is not obvious; that is, the relationship between pedestrian delay and vehicle flow is not significant. However, the vertical variation is obvious, indicating that the relationship between pedestrian delay and pedestrian flow is significant. In order to better explain the relationship between them, we will use cluster analysis to explain. In addition, according to the figure, we believe that clustering is carried out according to three categories. Table 5 shows the clustering center of vehicle flow and pedestrian delay. Table 6 shows clustering variance of vehicle flow and pedestrian delay. erefore, we can get the result that q � 600 pcu/h when the motor flow is small; q � 1500 pcu/h when the motor flow is moderate; q � 2400 pcu/h when the motor flow is large.

Data Availability
e traffic flow data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.