CN110930691A - Network traffic flow double-layer control method in full-automatic driving environment - Google Patents

Network traffic flow double-layer control method in full-automatic driving environment Download PDF

Info

Publication number
CN110930691A
CN110930691A CN201911011916.6A CN201911011916A CN110930691A CN 110930691 A CN110930691 A CN 110930691A CN 201911011916 A CN201911011916 A CN 201911011916A CN 110930691 A CN110930691 A CN 110930691A
Authority
CN
China
Prior art keywords
path
intersection
origin
phase
destination
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911011916.6A
Other languages
Chinese (zh)
Other versions
CN110930691B (en
Inventor
钱国敏
郭满
章立辉
王亦兵
王殿海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201911011916.6A priority Critical patent/CN110930691B/en
Publication of CN110930691A publication Critical patent/CN110930691A/en
Application granted granted Critical
Publication of CN110930691B publication Critical patent/CN110930691B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Traffic Control Systems (AREA)

Abstract

The invention discloses a double-layer control method of network traffic flow in a full-automatic driving environment. The method constructs a set of linear programming model on the upper layer to minimize the maximum time required by travel among all origin-destination points in the road network as a target, so as to control the departure time, the travel path and the signal scheme of the intersection of the automatic driving vehicle, realize that the traffic flow on all road sections in the road network presents periodic characteristics under the full automatic driving environment, and fully utilize the space-time resources of the road network as far as possible. The lower layer controls the longitudinal speed of the automatic driving vehicle to realize that all vehicles on the road section can pass through the intersection without stopping, and deduces a calculation formula of the shortest path section length required by passing through the intersection without stopping through a geometric relationship. Finally, road network traffic operation presents a stable state in a full-automatic driving environment, traffic flow on all road sections presents periodic characteristics in the state, and vehicles can pass through all signal intersections without stopping.

Description

Network traffic flow double-layer control method in full-automatic driving environment
Technical Field
The invention relates to a network traffic flow double-layer control method in a full-automatic driving environment, which is used for carrying out combined optimization on the departure time, the travel path, the longitudinal track of a road section and an intersection signal scheme of an automatic driving vehicle in a road network and belongs to the technical field of intelligent traffic.
Background
Vehicles with full-automatic driving technology have been tested in a closed test field in a large number, and some of the automatic driving automobiles of enterprises are tested on open urban roads, and some of the automatic driving automobiles even start commercial test operation. The advent of vehicle road coordination technology has further made possible the large-scale application of autonomous vehicles. The autonomous vehicles with communication capability may not only provide real-time information for road traffic status estimation, but also may be used directly as actuators for various traffic control strategies to improve traffic flow operation. Under the background that the automatic driving technology is more mature, a network traffic flow control method which can fully utilize the automatic driving vehicle technology is urgently needed to reasonably distribute and utilize road space-time resources so as to relieve traffic congestion, improve travel efficiency and save travel time of residents.
Disclosure of Invention
The invention aims to provide a double-layer control method of network traffic flow in a full-automatic driving environment. The core idea of the double-layer method is that a set of linear programming model which aims at minimizing the maximum time required by travel among all origin-destination points in the road network is constructed at the upper layer to control the departure time, the travel path and the intersection of the automatic driving vehicle, so that the traffic flow on all road sections in the road network presents periodic characteristics in the full-automatic driving environment, and the space-time resources of the road network are fully utilized as far as possible. The model mainly comprises a departure time distribution module, a trip path distribution module, a signal optimization module and a path and signal phase matching module. The lower layer controls the longitudinal speed of the automatic driving vehicle to realize that all vehicles on the road section can pass through the intersection without stopping, and deduces a calculation formula of the shortest path section length required by passing through the intersection without stopping through a geometric relationship. Finally, road network traffic operation presents a stable state in a fully automatic driving environment, traffic flow on all road sections presents periodic characteristics based on signal periods in the state, and vehicles can pass through all signal intersections without stopping.
The technical scheme adopted by the invention is as follows:
a network traffic flow double-layer control method under a full-automatic driving environment comprises the following steps:
c1, constructing a set of linear programming model aiming at minimizing the maximum time required by finishing all origin-destination opposite traveling to obtain the departure time, the traveling path and the timing scheme of each signalized intersection of the automatic driving vehicle.
c2, controlling the longitudinal speed of the automatic driving vehicle on each road section at the lower layer, and realizing that all vehicles on the road network can pass through the intersection without stopping.
c3, under the full-automatic driving environment, all vehicles travel according to the travel route obtained in the step c1 and the distributed travel time, and the network traffic flow stable state can be realized by adopting the longitudinal speed of the step c2 for driving. In this steady state, traffic flow on all road segments exhibits a periodic characteristic, and vehicles may pass through without stopping at all signal intersections. The stable state means that in one signal period, all automatic vehicles entering the road section ij from the intersection i can drive out of the road section ij through the intersection j without stopping in a certain period through longitudinal speed control. And the system is extended to the road network, the traffic flow on all road sections presents periodic characteristics, and vehicles can pass through all signal intersections without stopping.
Further, in step c1, the linear programming model aiming at minimizing the maximum time required for completing all origin-destination trips comprises a departure time allocation module for controlling a departure time of the vehicle, a trip path allocation module for controlling a trip path of the vehicle, a signal optimization module for optimizing signal timing, a path and signal phase pairing module for establishing a corresponding relationship between the trip path and an intersection signal phase, and an objective function module for minimizing the maximum time required for completing all origin-destination trips.
Further, the departure time distribution module controls the departure time of the automatic driving vehicle through the following constraint, and the situation that the vehicle travels intensively in time to cause congestion is avoided.
Figure BDA0002244450570000021
Figure BDA0002244450570000022
Figure BDA0002244450570000023
Figure BDA0002244450570000024
Figure BDA0002244450570000025
Figure BDA0002244450570000026
Figure BDA0002244450570000027
In the formula:
w represents a beginning-to-end set;
Ωwa set of paths representing origin-destination w;
i represents a certain intersection;
k represents the intersection phase number;
w represents a certain alignment, W belongs to W;
r denotes a path set ΩwIn one path, r ∈ Ωw
Figure BDA0002244450570000031
Representing the number of vehicles going out by the origin-destination w per period selected path r;
Figure BDA0002244450570000032
is a binary variable and is used for judging whether the r path of the origin-destination pair w passes through the k phase of the intersection i,
Figure BDA0002244450570000033
the r path representing the origin-destination w passes through the k phase of the intersection i, otherwise, the r path does not pass through the intersection i;
Figure BDA0002244450570000034
indicates the number of vehicles passing through the phase k at the intersection i per cycle on the r path of the origin-destination pair w,
Figure BDA0002244450570000035
when in use
Figure BDA0002244450570000036
Then
Figure BDA0002244450570000037
Otherwise
Figure BDA0002244450570000038
M is a large natural number which ensures that the number of vehicles which pass through the phase k at the intersection i in each period on the r path of the origin-destination point w does not exceed the number of vehicles which depart from the path in each period
Figure BDA0002244450570000039
nwRepresenting the number of vehicles departing per cycle of origin-destination w;
Figure BDA00022444505700000310
representing the number of lanes channelized to the k phase at the intersection i;
Figure BDA00022444505700000311
green time representing the phase k at the intersection i;
hsindicating a saturated headway.
Constraints (1) - (5) indicate that the number of vehicles arriving at the phase k of the intersection i per cycle does not exceed the allowable traffic capacity; constraint (6) indicates that the number of vehicles departing from each period of origin-destination w is equal to the sum of the number of vehicles departing from all paths, the vehicles depart at equal intervals in each period, and the interval is equal to the period duration/the number of vehicles departing from each period; the constraint (7) is an integer constraint.
Further, the travel path distribution module controls the travel path of the automatic driving vehicle through the following constraint, and the situation that the vehicles are distributed too intensively on the space to cause congestion is avoided.
Figure BDA00022444505700000312
Figure BDA00022444505700000313
Figure BDA00022444505700000314
Figure BDA00022444505700000315
Figure BDA0002244450570000041
Figure BDA0002244450570000042
Figure BDA0002244450570000043
Figure BDA0002244450570000044
Figure BDA0002244450570000045
In the formula:
n represents a road network node set, which comprises an origin-destination point set I' and an intersection set I, subscripts u and v are different nodes of a road network, and u, v belong to N and u is not equal to v;
Figure BDA0002244450570000046
indicates the order of the node v on the path r when
Figure BDA0002244450570000047
Indicating that node v is not on path r, where
Figure BDA0002244450570000048
Indicating the order of origin (o), (w) of origin (w);
Figure BDA0002244450570000049
is a binary variable which indicates whether the path r passes through the section vu or not when
Figure BDA00022444505700000410
Representing that the path r passes through the section vu, otherwise it does not pass;
Figure BDA00022444505700000411
is a binary parameter, which indicates whether the node v is the starting point of the origin-destination path r, when
Figure BDA00022444505700000412
Indicating that the node v is the starting point of the path r, otherwise, not;
Figure BDA00022444505700000413
is a binary parameter, which indicates whether the node v is the end point of the origin-destination path r, when
Figure BDA00022444505700000414
Indicating that the node v is the end point of the path r, otherwise, not;
a is a road segment set in a road network;
lvuthe length of the section vu belongs to A;
Figure BDA00022444505700000415
represents the shortest distance of origin-destination w;
β is length control parameter of origin-destination path r, β value is too large and is easy to be distributed to unreasonable path, β value is too small and affects solution space of model, and reasonable range (1,2) is recommended after repeatedly solving travel path distribution module.
Constraint (8) indicates the starting sequence number of each origin to w
Figure BDA00022444505700000416
Is set to 1; constraints (9) and (10) select the road segment through which the path r passes, the constraint (9) representing a path node number
Figure BDA00022444505700000417
Sequentially increasing by 1, the constraint (10) numbering nodes not on path r sequentially
Figure BDA00022444505700000418
Set to 0; constraints (11) and (12) prevent loops from occurring; constraints (13) and constraints (14) ensure connectivity of the paths; constraints (15) ensure the reasonability of the paths, preventing the length of the allocated paths from being much greater than the length of the shortest path; constraint (16) limiting
Figure BDA0002244450570000051
Is a binary variable.
Further, the signal optimization module performs signal timing optimization by the following constraints:
Figure BDA0002244450570000052
Figure BDA0002244450570000053
Figure BDA0002244450570000054
Figure BDA0002244450570000055
Figure BDA0002244450570000056
Cmin≤C≤Cmax(22)
Figure BDA0002244450570000057
Figure BDA0002244450570000058
in the formula:
Figure BDA0002244450570000059
respectively the green time of each phase at the intersection i, and the phases are numbered according to the NEMA phase structure;
c represents the signal period duration;
Cmin,Cmaxrespectively representing the minimum and maximum values of the signal period;
gmin,gmaxrespectively representing the minimum and maximum values of the phase duration.
Constraints (17) - (20) indicate that the green light durations of the paired phases are equal; constraint (21) means that the sum of the phase green duration of the same phase ring equals the signal period; constraints (22) and (23) limit cycle C and phase duration
Figure BDA00022444505700000510
The value range of (a); constraining (24) phase duration
Figure BDA00022444505700000511
Is an integer variable.
Further, the path and signal phase pairing module establishes a corresponding relationship between the travel path and the intersection signal phase through the following constraints.
Figure BDA00022444505700000512
Figure BDA00022444505700000513
Figure BDA0002244450570000061
Figure BDA0002244450570000062
In the formula:
a-(i, k) an entry leg representing a k phase at intersection i;
a+(i, k) represents the exit link for the k phase at i intersection.
When a certain path r of the origin-destination pair w passes through an inlet-outlet road section corresponding to the k phase of the intersection i at the same time, the vehicle which selects the path r to travel passes through the intersection i at the k phase, namely
Figure BDA0002244450570000063
Then
Figure BDA0002244450570000064
Otherwise, the path r does not pass through the k phase at the i intersection, i.e.
Figure BDA0002244450570000065
Further, the objective function module minimizes the maximum time required to complete all origin-destination trips by the following constraints:
maxrm
Figure BDA0002244450570000066
Figure BDA0002244450570000067
in the formula:
rm represents the minimum value of the reciprocal of the number of cycles required for completing the travel demand at all the origin-destination pairs;
rmwrepresenting the reciprocal of the number of cycles required for origin-destination pair w to complete its fixed travel demand;
Nwrepresenting a fixed travel demand of origin-destination w;
the objective function maxrm maximizes the minimum value of the reciprocal of the number of cycles required for all the origin-destination pairs to complete the travel demand, and is equivalent to minimizing the maximum time required for all the origin-destination pairs to complete the travel demand, the constraint (29) represents the reciprocal of the number of cycles required for the origin-destination pair w to complete the travel demand, and the constraint (30) represents the minimum value of the reciprocal of the number of cycles.
Further, in the step c2, a calculation formula of the shortest path length required for realizing the non-stop passing through the intersection through the longitudinal speed control is given, and on the section meeting the minimum path length, a reasonable speed is selected to ensure that the automatic driving vehicle on the section can pass through the intersection without stopping, and the calculation formula of the shortest path length is as follows:
Figure BDA0002244450570000068
in the formula:
Figure BDA0002244450570000071
representing the length of the shortest path segment required by controlling the longitudinal speed so as to pass through the intersection without stopping;
Vmin,Vmaxrespectively representing the minimum and maximum speeds allowed for the section vu.
The invention has the beneficial effects that:
the invention provides a network traffic flow double-layer control method in a full-automatic driving environment, which is characterized in that a set of linear programming model aiming at minimizing the maximum time required by finishing all origin-destination opposite traveling is constructed for the first time to obtain a signal scheme of a departure time, a traveling path and an intersection of an automatic driving vehicle, and the longitudinal speed of the automatic driving vehicle on each road section is controlled at the same time, so that all vehicles on a network can pass through the intersection without stopping, the network traffic flow periodically presents a stable state without stopping and passing at the intersection, and the method fully utilizes road space-time resources, minimizes the total time of all origin-destination opposite traveling, relieves traffic congestion, improves traveling efficiency and saves traveling time of residents.
Drawings
FIG. 1 traffic signal NEMA phase structure;
FIG. 2 is a sample vehicle trajectory;
FIG. 3 is a road network topology.
Detailed Description
The invention provides a network traffic flow double-layer control method in a full-automatic driving environment, which specifically comprises the following steps:
c1, constructing a set of linear programming model aiming at minimizing the maximum time required by finishing all origin-destination opposite traveling to obtain the departure time, the traveling path and the timing scheme of each signalized intersection of the automatic driving vehicle.
c2, controlling the longitudinal speed of the automatic driving vehicle on each road section at the lower layer, and realizing that all vehicles on the road network can pass through the intersection without stopping.
c3, under the full-automatic driving environment, all vehicles travel according to the travel route obtained in the step c1 and the distributed travel time, and the network traffic flow stable state can be realized by adopting the longitudinal speed of the step c2 for driving. In this steady state, traffic flow on all road segments exhibits a periodic characteristic, and vehicles may pass through without stopping at all signal intersections.
As a preferred scheme, the linear programming model includes a departure time allocation module for controlling a departure time of the vehicle, a travel path allocation module for controlling a travel path of the vehicle, a signal optimization module for signal timing optimization, a path and signal phase pairing module for establishing a corresponding relationship between the travel path and a signal phase at an intersection, and an objective function module for minimizing a maximum time required for completing all origin-destination trips.
In addition, a calculation formula for the shortest path length required by realizing the non-stop passing of the intersection through longitudinal speed control is given, and on the section meeting the minimum path length, the reasonable speed is selected to ensure that the automatic driving vehicle on the section can pass through the intersection without stopping:
suppose that a distance ij defines an allowable speed range (V)min,Vmax) In which V isminIs greater than 0. When the road section length is long enough and the number of arriving vehicles does not exceed the traffic capacity of the intersection, the vehicle passes through the position (V)min,Vmax) The proper longitudinal speed is selected to ensure that an automatic vehicle entering the section ij in the same signal cycle can drive out of the section ij through the intersection j without stopping, as shown in fig. 2.
FIG. 2 shows a diagram of the trajectories of automatic vehicles on a section ij passing through an intersection j without stopping, these automatic vehicles having different longitudinal speeds, entering the section ij from the intersection i
Figure BDA0002244450570000081
The phase is driven off the road section ij through the intersection j, wherein it is assumed
Figure BDA0002244450570000082
In order to turn the phase to the left,
Figure BDA0002244450570000083
is a straight line phase. If the first vehicle turns left in each cycle at speed vfIn that
Figure BDA0002244450570000084
After the starting moment, the vehicle leaves the intersection j, and the last left-turn automatic vehicle is driven at the speed vlIn that
Figure BDA0002244450570000085
Drive off intersection j at a time before the end, and Vf,Vl∈(Vmin,Vmax) Then all left-turn automobiles entering the road section ij in the same period on the road section ij can be distributed to proper speed so as to pass through the intersection j without stopping; similarly, for a straight-ahead vehicle on a section ij, if it can be (V)min,Vmax) Also find suitable in the scopeThe speed of the first and last straight-going automatic vehicles entering in the same period can drive out of the road section ij without stopping through the intersection j, and all the straight-going automatic vehicles entering in the same period can be distributed to proper longitudinal speed to drive out of the road section ij without stopping.
The conditions ① can be met by controlling the departure time, the travel path and the signal scheme through an upper layer, and the shortest path calculation formula required to be met by the condition ② road section is given as follows:
Figure BDA0002244450570000086
calculating the allowable speed range (V) according to the shortest path length calculation formulamin,Vmax) The shortest path length in between. When the length of the path is larger than
Figure BDA0002244450570000087
In time, an appropriate speed V (V) can be assigned according to the time when each automatic vehicle enters the section vumin≤V≤Vmax) Thereby ensuring that all automatic vehicles can drive away from the road section vu without stopping.
Taking a 4 × 4 road network as an example, as shown in fig. 3, the road network has 32 nodes, 16 intersections, 40 road segments, and 32 alignment points. The implementation steps are as follows:
step 1: solving the linear programming model to obtain the departure time, the travel path and the signal timing scheme of each intersection of each automatic vehicle, and the method comprises the following specific steps:
step 1.1: constructing a road network topology and a network traffic flow set, wherein the road network topology and the network traffic flow set comprise an intersection set I, a road network node set N, a road section set A, an origin-destination set W and a phase in-out road section matching set AP, and examples are shown in table 1:
table 1 road network topology and network traffic flow set examples
Collection name Aggregated content
Intersection set I {1,2,3,4,…,14,15,16}
Road network node set N {1,2,3,…,30,31,32}
Road section set A {(1,2),(2,1),(1,5),(5,1)…}
Origin-destination set W {17→25,25→17,19→27,27→19…}
Phase access road section matching set AP g1Phase { (10,1,4), (11,2,5) … } …
Step 1.2: investigating to obtain various model parameters including road length laOrigin-destination fixed traffic demand NwMinimum and maximum values of the period duration Cmin,CmaxMinimum and maximum phase duration gmin,gmaxNumber of channelized lanes per intersection phase k
Figure BDA0002244450570000091
Saturated headway hs
Step 1.3: calculating the shortest path length between each pair of origin-destination pairs by using Dijkstra algorithm
Figure BDA0002244450570000092
The number of paths r and path length control parameter β between each pair of origin-destination pairs, for example, origin-destination pair 17 → 25 has two paths assigned thereto, the shortest path length is 4155m, the path length control parameter β is set to 1.5, and the OD shortest path lengths are as shown in Table 2:
TABLE 2 shortest Path Length of each OD of road network
OD pair 17→25 18→26 19→27 20→28 21→29 22→30 23→31 24→32
Shortest Path Length (m) 4155 3255 3075 4215 4605 3660 3435 4155
OD pair 25→17 26→18 27→19 28→20 29→21 30→22 31→23 32→24
Shortest Path Length (m) 4155 3255 3075 4215 4605 3660 3435 4155
OD pair 17→29 18→30 19→31 20→32 25→21 26→22 27→23 28→24
Shortest Path Length (m) 3240 2790 2760 3150 3015 2715 2220 2790
OD pair 21→17 22→18 23→19 24→20 29→25 30→26 31→27 32→28
Shortest Path Length (m) 3120 2880 2745 3015 3210 3090 2325 2850
Step 1.4: and solving the linear programming model to obtain a timing scheme of each signalized intersection in the network, the number of the departure vehicles in each period and the number of the departure vehicles in each period of each path. For example, the two paths of origin-destination 19 → 27 are {19,3,2,6,5,9,13,14,27} and {19,3,2,1,5,9,10,14,27} respectively, and 24 and 4 vehicles are respectively started per cycle, and 28 vehicles are always started, and other paths and flow and intersection signal schemes are as shown in tables 3 and 4:
TABLE 3 Path and traffic distribution results
Figure BDA0002244450570000101
Figure BDA0002244450570000111
TABLE 4 intersection Signaling scheme
Figure BDA0002244450570000112
Step 1.5: all vehicles are randomly sequenced, departure time is distributed at equal intervals according to the number of departure vehicles in each period, and departure paths of the vehicles are distributed at equal intervals according to the proportion of the number of departure vehicles in each period of each path.
For example, the departure interval of origin-destination pair 19 → 27 is the period duration/number of vehicles departing per period, i.e. 120/28 equals 4.3s, the first vehicle departs at the beginning of the period, the second vehicle departs at the 4.3s of the period, the third vehicle departs at the 8.6s of the period, and so on, the departure time of each vehicle is obtained. And one of every 7 vehicles is assigned to the second path, and assuming that the first vehicle and the remaining six vehicles are assigned to the second path, the first vehicle of the next 7 vehicles is still assigned to the second path.
Step 2: according to the sequence of each vehicle entering the road section in the same-flow traffic flow, the cycle time of the road section which is driven away at the earliest at the downstream intersection is determined, and the conflict is prevented. For example, a 5 th vehicle turning left at the intersection 7 may initially leave the intersection at the 5 th s time of the left turn phase, otherwise a collision with the preceding vehicle may occur.
And step 3: selecting a suitable speed range (V) for each road section according to the length of the road sectionmin,Vmax) And determining the periodic moment of each vehicle actually driving away from the road section and distributing proper longitudinal speed to ensure that the vehicle does not stop driving away from the intersection.
For example, a 5 th vehicle turning left at the intersection 7 can initially leave the intersection at the 5 th s of the left-turn phase, but the assigned speed V > VmaxNot meeting the speed requirement, and only at this time, the speed V ismaxDriving off the intersection at a certain moment after 5s at the earliest; and if the 5 th left-turn vehicle can leave the intersection at the 5 th s moment of the left-turn phase, the speed of the road section is distributed according to the time so as to ensure that the vehicle passes through without stopping.
The invention is not the best known technology.
The above embodiments are merely illustrative of the technical ideas and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention, and not to limit the protection scope of the present invention. All equivalent changes and modifications made according to the spirit of the present invention should be covered within the protection scope of the present invention.

Claims (8)

1. A network traffic flow double-layer control method in a full-automatic driving environment is characterized by comprising the following steps:
c1, constructing a set of linear programming model aiming at minimizing the maximum time required by finishing all origin-destination opposite traveling to obtain the departure time, the traveling path and the timing scheme of each signalized intersection of the automatic driving vehicle.
c2, controlling the longitudinal speed of the automatic driving vehicle on each road section at the lower layer, and realizing that all vehicles on the road network can pass through the intersection without stopping.
c3, under the full-automatic driving environment, all vehicles travel according to the travel route obtained in the step c1 and the distributed travel time, and the network traffic flow stable state can be realized by adopting the longitudinal speed of the step c2 for driving. In this steady state, traffic flow on all road segments exhibits a periodic characteristic, and vehicles may pass through without stopping at all signal intersections.
2. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 1, characterized in that: in step c1, the linear programming model aiming at minimizing the maximum time required for completing all origin-destination trips comprises a departure time distribution module for controlling a departure time of the vehicle, a trip path distribution module for controlling a trip path of the vehicle, a signal optimization module for optimizing signal timing, a path-signal phase pairing module for establishing a corresponding relationship between the trip path and an intersection signal phase, and an objective function module for minimizing the maximum time required for completing all origin-destination trips.
3. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 2, characterized in that: the departure time allocation module controls the departure time of the autonomous vehicle by the following constraints:
Figure FDA0002244450560000011
Figure FDA0002244450560000012
Figure FDA0002244450560000013
Figure FDA0002244450560000014
Figure FDA0002244450560000015
Figure FDA0002244450560000016
Figure FDA0002244450560000017
in the formula:
w represents a beginning-to-end set;
Ωwa set of paths representing origin-destination w;
i represents a certain intersection;
k represents the intersection phase number;
w represents a certain alignment, W belongs to W;
r denotes a path set ΩwIn one path, r ∈ Ωw
Figure FDA0002244450560000021
Representing the number of vehicles going out by the origin-destination w per period selected path r;
Figure FDA0002244450560000022
is a binary variable and is used for judging whether the r path of the origin-destination pair w passes through the k phase of the intersection i,
Figure FDA0002244450560000023
the r path representing the origin-destination w passes through the k phase of the intersection i, otherwise, the r path does not pass through the intersection i;
Figure FDA0002244450560000024
indicates the number of vehicles passing through the phase k at the intersection i per cycle on the r path of the origin-destination pair w,
Figure FDA0002244450560000025
when in use
Figure FDA0002244450560000026
Then
Figure FDA0002244450560000027
Otherwise
Figure FDA0002244450560000028
M is a large natural number which ensures that the number of vehicles which pass through the phase k at the intersection i in each period on the r path of the origin-destination point w does not exceed the number of vehicles which depart from the path in each period
Figure FDA0002244450560000029
nwRepresenting the number of vehicles departing per cycle of origin-destination w;
Figure FDA00022444505600000210
representing the number of lanes channelized to the k phase at the intersection i;
Figure FDA00022444505600000211
green time representing the phase k at the intersection i;
hsindicating a saturated headway.
Constraints (1) - (5) indicate that the number of vehicles arriving at the phase k of the intersection i per cycle does not exceed the allowable traffic capacity; constraint (6) means that the number of vehicles departing from origin-destination point w per period is equal to the sum of the number of vehicles departing from all paths, and the vehicles depart at equal intervals per period; the constraint (7) is an integer constraint.
4. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 2, characterized in that: the travel path distribution module controls the travel path of the automatic driving vehicle through the following constraints:
Figure FDA00022444505600000212
Figure FDA00022444505600000213
Figure FDA00022444505600000214
Figure FDA00022444505600000215
Figure FDA0002244450560000031
Figure FDA0002244450560000032
Figure FDA0002244450560000033
Figure FDA0002244450560000034
Figure FDA0002244450560000035
in the formula:
n represents a road network node set, which comprises an origin-destination point set I' and an intersection set I, subscripts u and v are different nodes of a road network, and u, v belong to N and u is not equal to v;
Figure FDA0002244450560000036
indicates the order of the node v on the path r when
Figure FDA0002244450560000037
Indicating that node v is not on path r, where
Figure FDA0002244450560000038
Indicating the order of origin (o), (w) of origin (w);
Figure FDA0002244450560000039
is a binary variable which indicates whether the path r passes through the section vu or not when
Figure FDA00022444505600000310
Representing that the path r passes through the section vu, otherwise it does not pass;
Figure FDA00022444505600000311
is a binary parameter, which indicates whether the node v is the starting point of the origin-destination path r, when
Figure FDA00022444505600000312
Indicating that the node v is the starting point of the path r, otherwise, not;
Figure FDA00022444505600000313
is a binary parameter, which indicates whether the node v is the end point of the origin-destination path r, when
Figure FDA00022444505600000314
Indicating that the node v is the end point of the path r, otherwise, not;
a is a road segment set in a road network;
lvuthe length of the section vu belongs to A;
Figure FDA00022444505600000315
represents the shortest distance of origin-destination w;
β is the length control parameter of the origin-destination path r, with a reasonable range of (1,2) being recommended.
Constraint (8) indicates the starting sequence number of each origin to w
Figure FDA00022444505600000316
Is set to 1; constraints (9) and (10) select the road segment through which the path r passes, the constraint (9) representing a path node number
Figure FDA00022444505600000317
Sequentially increasing by 1, the constraint (10) numbering nodes not on path r sequentially
Figure FDA00022444505600000318
Set to 0; constraints (11) and (12) prevent loops from occurring(ii) a Constraints (13) and constraints (14) ensure connectivity of the paths; constraints (15) ensure the reasonability of the paths, preventing the length of the allocated paths from being much greater than the length of the shortest path; constraint (16) limiting
Figure FDA00022444505600000319
Is a binary variable.
5. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 2, characterized in that: the signal optimization module performs signal timing optimization through the following constraints:
Figure FDA0002244450560000041
Figure FDA0002244450560000042
Figure FDA0002244450560000043
Figure FDA0002244450560000044
Figure FDA0002244450560000045
Cmin≤C≤Cmax(22)
Figure FDA0002244450560000046
Figure FDA0002244450560000047
in the formula:
Figure FDA0002244450560000048
respectively the green time of each phase at the intersection i, and the phases are numbered according to the NEMA phase structure;
c represents the signal period duration;
Cmin,Cmaxrespectively representing the minimum and maximum values of the signal period;
gmin,gmaxrespectively representing the minimum and maximum values of the phase duration.
Constraints (17) - (20) indicate that the green light durations of the paired phases are equal; constraint (21) means that the sum of the phase green duration of the same phase ring equals the signal period; constraints (22) and (23) limit cycle C and phase duration
Figure FDA0002244450560000049
The value range of (a); constraining (24) phase duration
Figure FDA00022444505600000410
Is an integer variable.
6. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 2, characterized in that: the path and signal phase pairing module establishes a corresponding relation between a travel path and an intersection signal phase through the following constraints.
Figure FDA00022444505600000411
Figure FDA00022444505600000412
Figure FDA00022444505600000413
Figure FDA0002244450560000051
In the formula:
a-(i, k) an entry leg representing a k phase at intersection i;
a+(i, k) represents the exit link for the k phase at i intersection.
When a certain path r of the origin-destination pair w passes through an inlet-outlet road section corresponding to the k phase of the intersection i at the same time, the vehicle which selects the path r to travel passes through the intersection i at the k phase, namely
Figure FDA0002244450560000052
Then
Figure FDA0002244450560000053
Otherwise, the path r does not pass through the k phase at the i intersection, i.e.
Figure FDA0002244450560000054
7. The network traffic flow double-layer control method in the full-automatic driving environment according to claim 2, characterized in that: the objective function module achieves the minimum maximum time required for completing all origin-destination opposite trips through the following constraints:
max rm
Figure FDA0002244450560000055
Figure FDA0002244450560000056
in the formula:
rm represents the minimum value of the reciprocal of the number of cycles required for completing the travel demand at all the origin-destination pairs;
rmwrepresenting the reciprocal of the number of cycles required for origin-destination pair w to complete its fixed travel demand;
Nwrepresenting a fixed travel demand of origin-destination w;
the objective function maxrm maximizes the minimum value of the reciprocal of the number of cycles required for all the origin-destination pairs to complete the travel demand, and is equivalent to minimizing the maximum time required for all the origin-destination pairs to complete the travel demand, the constraint (29) represents the reciprocal of the number of cycles required for the origin-destination pair w to complete the travel demand, and the constraint (30) represents the minimum value of the reciprocal of the number of cycles.
8. The road network traffic flow double-layer control method under the full-automatic driving environment according to claim 1, characterized in that: in the step c2, a calculation formula of the shortest path length required for realizing the non-stop passing through the intersection through the longitudinal speed control is given, and on the section meeting the shortest path length, the reasonable speed is selected to ensure that the automatic driving vehicle on the section can pass through the intersection without stopping, and the calculation formula of the shortest path length is as follows:
Figure FDA0002244450560000057
in the formula:
Figure FDA0002244450560000061
representing the length of the shortest path segment required by controlling the longitudinal speed so as to pass through the intersection without stopping;
Vmin,Vmaxrespectively representing the minimum and maximum speeds allowed for the section vu.
CN201911011916.6A 2019-10-23 2019-10-23 Network traffic flow double-layer control method in full-automatic driving environment Active CN110930691B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911011916.6A CN110930691B (en) 2019-10-23 2019-10-23 Network traffic flow double-layer control method in full-automatic driving environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911011916.6A CN110930691B (en) 2019-10-23 2019-10-23 Network traffic flow double-layer control method in full-automatic driving environment

Publications (2)

Publication Number Publication Date
CN110930691A true CN110930691A (en) 2020-03-27
CN110930691B CN110930691B (en) 2021-02-23

Family

ID=69849194

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911011916.6A Active CN110930691B (en) 2019-10-23 2019-10-23 Network traffic flow double-layer control method in full-automatic driving environment

Country Status (1)

Country Link
CN (1) CN110930691B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113936452A (en) * 2020-06-29 2022-01-14 阿里巴巴集团控股有限公司 Traffic control method and device and electronic equipment
CN115083177A (en) * 2022-06-28 2022-09-20 浙江大学 Control method for realizing intelligent vehicle series arrangement on three lanes and above

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB432900A (en) * 1934-02-05 1935-08-06 Automatic Electric Co Ltd Improvements in or relating to traffic control signals
DE3126481C2 (en) * 1981-07-04 1989-10-19 Volkswagen Ag
JPH11241622A (en) * 1998-02-24 1999-09-07 Honda Motor Co Ltd Engine driving force controller
DE19938266A1 (en) * 1999-08-12 2001-02-15 Volkswagen Ag Method and equipment for visual recognition of traffic signs, electronically identifies obscure traffic signs by comparison with known signs
CN1526596A (en) * 2003-01-30 2004-09-08 Rail transportation system of small vehicle
US20040225434A1 (en) * 2003-05-07 2004-11-11 Gotfried Bradley L. Vehicle navigation and safety systems
CN102136104A (en) * 2011-03-22 2011-07-27 西安电子科技大学 Load balance and Lin-Kernighan (LK) algorithm based vehicle route planning method
CN103236177A (en) * 2013-03-26 2013-08-07 清华大学 Intelligent interactive system with vehicular network multi-system fusion, and control method thereof
CN103383812A (en) * 2013-06-27 2013-11-06 姜廷顺 System and method for releasing travel time by being accurate to second
CN104021675A (en) * 2014-06-25 2014-09-03 北京易华录信息技术股份有限公司 System and method for predicating travel time needed by express way in future time
CN205721998U (en) * 2016-02-18 2016-11-23 华南理工大学 Turn pike net traffic information collection based on path identifying system and inducible system
CN106251016A (en) * 2016-08-01 2016-12-21 南通大学 A kind of parking system paths planning method based on dynamic time windows
CN106809223A (en) * 2017-02-16 2017-06-09 常文作 A kind of intelligent track traffic system
CN108647475A (en) * 2018-05-25 2018-10-12 东南大学 City discrete network design problem R language implementation methods based on load balancing
CN108711300A (en) * 2018-05-25 2018-10-26 深圳市易成自动驾驶技术有限公司 Swarm intelligence is parked method, system and computer readable storage medium
CN108804801A (en) * 2018-05-25 2018-11-13 东南大学 City discrete network design problem R language implementation methods based on target flow
CN108961784A (en) * 2018-06-11 2018-12-07 江苏济通智慧科技有限公司 Speed and intersection signal bidirection control method under a kind of automatic Pilot

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB432900A (en) * 1934-02-05 1935-08-06 Automatic Electric Co Ltd Improvements in or relating to traffic control signals
DE3126481C2 (en) * 1981-07-04 1989-10-19 Volkswagen Ag
JPH11241622A (en) * 1998-02-24 1999-09-07 Honda Motor Co Ltd Engine driving force controller
DE19938266A1 (en) * 1999-08-12 2001-02-15 Volkswagen Ag Method and equipment for visual recognition of traffic signs, electronically identifies obscure traffic signs by comparison with known signs
CN1526596A (en) * 2003-01-30 2004-09-08 Rail transportation system of small vehicle
US20040225434A1 (en) * 2003-05-07 2004-11-11 Gotfried Bradley L. Vehicle navigation and safety systems
CN102136104A (en) * 2011-03-22 2011-07-27 西安电子科技大学 Load balance and Lin-Kernighan (LK) algorithm based vehicle route planning method
CN103236177A (en) * 2013-03-26 2013-08-07 清华大学 Intelligent interactive system with vehicular network multi-system fusion, and control method thereof
CN103383812A (en) * 2013-06-27 2013-11-06 姜廷顺 System and method for releasing travel time by being accurate to second
CN104021675A (en) * 2014-06-25 2014-09-03 北京易华录信息技术股份有限公司 System and method for predicating travel time needed by express way in future time
CN205721998U (en) * 2016-02-18 2016-11-23 华南理工大学 Turn pike net traffic information collection based on path identifying system and inducible system
CN106251016A (en) * 2016-08-01 2016-12-21 南通大学 A kind of parking system paths planning method based on dynamic time windows
CN106809223A (en) * 2017-02-16 2017-06-09 常文作 A kind of intelligent track traffic system
CN108647475A (en) * 2018-05-25 2018-10-12 东南大学 City discrete network design problem R language implementation methods based on load balancing
CN108711300A (en) * 2018-05-25 2018-10-26 深圳市易成自动驾驶技术有限公司 Swarm intelligence is parked method, system and computer readable storage medium
CN108804801A (en) * 2018-05-25 2018-11-13 东南大学 City discrete network design problem R language implementation methods based on target flow
CN108961784A (en) * 2018-06-11 2018-12-07 江苏济通智慧科技有限公司 Speed and intersection signal bidirection control method under a kind of automatic Pilot

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
SUNEET KUMAR SINGH ET AL.: "Impact Analysis of Start-Up Lost Time at Major Intersections on Sathorn Road Using a Synchro Optimization and a Microscopic SUMO Traffic Simulation", 《IEEE》 *
TAKASHI NAGATANI: "Effect of speed fluctuations on a green-light path in a 2d traffic network controlled by signals", 《EL SEVIER》 *
YIMING BIE ET AL.: "Modelling correlation degree between two adjacent signalised intersections for dynamic subarea partition", 《IEEE》 *
保丽霞: "基于信息集成的城市交通流诱导与交通控制协同的关键理论及技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *
吴伟 等: "车路协同环境下基于路径的信号协调优化模型", 《吉林大学学报(工学版)》 *
张得志 等: "在时变速度下的车辆出行时刻和路径联合优化", 《铁道科学与工程学报》 *
赵小方: "基于博弈论的交通控制和动态网络均衡组合优化研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113936452A (en) * 2020-06-29 2022-01-14 阿里巴巴集团控股有限公司 Traffic control method and device and electronic equipment
CN115083177A (en) * 2022-06-28 2022-09-20 浙江大学 Control method for realizing intelligent vehicle series arrangement on three lanes and above
CN115083177B (en) * 2022-06-28 2024-05-10 浙江大学 Control method for realizing intelligent vehicle series arrangement on three lanes and above

Also Published As

Publication number Publication date
CN110930691B (en) 2021-02-23

Similar Documents

Publication Publication Date Title
CN106875710B (en) A kind of intersection self-organization control method towards net connection automatic driving vehicle
CN105096622B (en) A kind of dynamic route guidance method communicated based on bus or train route
CN106548646B (en) Road information service system and method when being blocked up based on the city that mist calculates
CN104766484B (en) Traffic Control and Guidance system and method based on Evolutionary multiobjective optimization and ant group algorithm
CN105938572B (en) A kind of more automatic guided vehicle paths planning methods of the pre- anti-interference of logistics storage system
CN109410619B (en) Automatic traffic control method and system under one-way road network
CN107016857B (en) Signal control intersection left-turn traffic combination design optimization method
CN103000027B (en) Intelligent traffic guidance method based on floating car under congestion condition
CN101593419B (en) Intelligent coordination control method for traffic flow of urban road network with public traffic priority
CN103578281B (en) A kind of main line of communication signal lamp optimal control method and device
CN110930691B (en) Network traffic flow double-layer control method in full-automatic driving environment
CN102307373B (en) VANET clustering method taking regard of vehicle traffic characteristic
CN104916142B (en) A kind of crossing self-adapting traffic signal control method of giving priority to trunk roads
CN107507430A (en) A kind of urban road crossing traffic control method and system
CN108647802B (en) Anti-congestion method based on double-layer traffic network model
CN107389076A (en) A kind of real-time dynamic path planning method of energy-conservation suitable for intelligent network connection automobile
CN104080056B (en) The message distributing method of the vehicular ad hoc network perceived based on degree of communication probability
CN109887289A (en) A kind of network vehicle flowrate maximization approach of urban traffic network model
CN110136456A (en) Traffic lights anti-clogging control method and system based on deeply study
CN104637328B (en) Distributed real-time navigation method based on RSU in a kind of vehicle self-organizing network
CN102982688A (en) Regional traffic signal control method based on arterial road coordination priority
CN107490384A (en) A kind of optimal static path system of selection based on city road network
CN108629455A (en) A kind of real-time route planing method based on vehicle self-organizing net
CN107343024A (en) A kind of centralized car networking MAC layer merges prediction of collision and avoiding method
CN103236169A (en) Bidirectional segmental green-wave signal setting method for urban main line bus stops

Legal Events

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