CN110599772B - Mixed traffic flow cooperative optimization control method based on double-layer planning - Google Patents

Mixed traffic flow cooperative optimization control method based on double-layer planning Download PDF

Info

Publication number
CN110599772B
CN110599772B CN201910887658.1A CN201910887658A CN110599772B CN 110599772 B CN110599772 B CN 110599772B CN 201910887658 A CN201910887658 A CN 201910887658A CN 110599772 B CN110599772 B CN 110599772B
Authority
CN
China
Prior art keywords
vehicle
road
stage
model
time
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.)
Active
Application number
CN201910887658.1A
Other languages
Chinese (zh)
Other versions
CN110599772A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201910887658.1A priority Critical patent/CN110599772B/en
Publication of CN110599772A publication Critical patent/CN110599772A/en
Application granted granted Critical
Publication of CN110599772B publication Critical patent/CN110599772B/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
    • G08G1/0145Measuring and analyzing of parameters relative to traffic conditions for specific applications for active traffic flow control
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/075Ramp control
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

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

Abstract

The invention provides a double-layer planning-based hybrid traffic flow collaborative optimization control method, and belongs to the field of traffic engineering. The method adopts a double-layer optimization model based on dynamic programming recursion to carry out mixed traffic flow cooperative decision control, and is suitable for different traffic scenes in the mixed traffic flow, including expressway ramp vehicle convergence, intersection vehicle convergence and vehicle passing intersection; the double-layer optimization model comprises an upper layer model and a lower layer model, wherein the upper layer model is a vehicle sequencing problem solved by dynamic programming recursion, and the lower layer model is a single vehicle track optimization problem in different traffic scenes solved by dynamic programming recursion; and establishing an upper layer model and a lower layer model, and solving the models. The upper layer model and the lower layer model jointly ensure the optimal running of the system vehicles, so that the vehicle conflict of the vehicles in the process of converging or passing through the intersection is reduced under the mixed traffic flow environment, and the passing efficiency and the comfort of the vehicles are effectively improved.

Description

Mixed traffic flow cooperative optimization control method based on double-layer planning
Technical Field
The invention relates to a mixed traffic flow collaborative optimization control method based on double-layer planning, and belongs to the field of traffic engineering.
Background
New technologies such as Global Positioning System (GPS), wireless communication, advanced sensing, and automatic control have prompted the rapid development of autonomous vehicles (i.e., intelligent networked vehicles, which are optimally controllable vehicles, CAVs). An autonomous vehicle is defined as a vehicle that is capable of sensing and communicating with the driving environment, and the operation of the vehicle (in part or in whole) may be performed without driver action. Autonomous vehicles have better controllability and cooperativity than the difficult cooperativity of conventional human-driven vehicles (non-optimally controlled vehicles). This may therefore provide benefits such as improved fuel/energy efficiency, traffic safety and traffic stability. Mixed traffic (i.e., a state in which human-driven vehicles and autonomous vehicles are mixed) will become a major mode of road traffic before the hundreds of popularity of autonomous vehicles. Under the environment of mixed traffic flow, under the conditions of various traffic scenes such as ramp junction convergence at a highway intersection, junction intersection or junction at a T-junction and the like, traffic conflict between an automatic driving vehicle and a human driving vehicle can occur, and certain harm is caused. The scientific theoretical framework and the modeling method are used for carrying out mixed traffic flow cooperative decision control, so that traffic conflicts are reduced or eliminated, and the vehicle track is optimized to a certain degree to become a practical problem which needs to face in future traffic for a long time.
Existing research related to autonomous vehicles is mostly based on the assumption that the permeability of the autonomous vehicle is 100%, and is mostly studied from a macroscopic perspective. The model of the microscopic level is also mainly focused on the single-vehicle track optimization of the research automatic driving vehicle, and the optimization of the system level cannot be guaranteed. From a microscopic view point, a mixed traffic flow cooperative decision control optimization model aiming at system optimization is researched. The cooperative decision control under the mixed traffic flow environment is carried out aiming at various microscopic traffic scenes, such as crossroads, T-shaped intersections, expressway ramp confluence and the like, so that traffic conflicts are eliminated to the maximum extent, and the research on traffic operation efficiency and traffic capacity is basically not available.
Disclosure of Invention
The invention aims to provide a mixed traffic flow cooperative optimization control method based on double-layer planning, aiming at the defects in the prior art, and from a microscopic view, researching a mixed traffic flow cooperative decision control optimization model aiming at system optimization. The mixed traffic flow refers to the traffic flow formed by mixing a plurality of vehicles which can be optimally controlled and vehicles which can not be optimally controlled.
The technical scheme adopted by the invention for realizing the aim is as follows:
a mixed traffic flow collaborative optimization control method based on double-layer planning is characterized in that a double-layer optimization model based on dynamic planning recursion is adopted for carrying out mixed traffic flow collaborative decision control, and the method is suitable for different traffic scenes in the mixed traffic flow, wherein the different traffic scenes comprise ramp vehicle convergence on a highway, vehicle confluence at an intersection and vehicle passing through the intersection; the double-layer optimization model comprises an upper layer model and a lower layer model; the upper model is a vehicle sequencing problem solved by dynamic programming recursion; the lower layer model is a single vehicle track optimization problem in different traffic scenes solved by dynamic programming recursion;
the method comprises the following steps:
s1, establishing the upper layer model, including:
s1-1, determining the division stage, the state variable and the decision variable of the upper layer model, which are as follows:
the upper layer model is divided: assuming that an X road and a Y road are two one-way roads with intersections, n vehicles on the X road need to sequentially converge into or pass through m +1 intervals among the vehicles on the Y road, and the behavior that each single vehicle on the X road converges into or passes through the Y road is represented as a stage, and the behavior that the k-th vehicle on the X road converges into or passes through the Y road is recorded as a k-th stage, wherein k is 1,2,3, …, n;
state variables of the upper model: the number of the vehicles which can be converged into the kth vehicle on the X road or pass through the Y road in the kth stage is skRepresents;
decision variables of the upper model: the decision made in each stage represents that the k vehicle on the X road in the k stage can merge into or pass through the s road of the Y roadkSelecting a particular xth of individual vehicle intervalskThe vehicles are alternately converged or passed;
s1-2, determining a state transition equation, a cost function and an objective function of the upper layer model, which are as follows:
the state transition equation of the upper model is as follows:
Figure BDA0002207808210000021
setting the initial condition as s0=m+1;
The state transition equation of the upper model indicates that s is equal to 1 when kk=s1The number of the vehicle intervals for the 1 st vehicle on the X road to enter or pass through the Y road in the 1 st stage is m + 1; when k is 2,3, …, n, the (k-1) th vehicle on the X road selects the (X) th vehicle at the (k-1) th stagek-1The interval of each vehicle is used as a state variable s after the vehicle is converged or passes through a Y roadkA change in (c); s0M +1 represents that the number of vehicle intervals for vehicles to merge into or pass through the road Y on the road X in the initial state is m + 1;
cost function of the upper model:
Figure BDA0002207808210000022
cost function D of the upper modelk(sk,xk) Indicating a phase index function required for making a decision at the kth phase, wherein
Figure BDA0002207808210000023
S representing that the k-th vehicle on the X road is merged into or passes through the road where the k-th vehicle can merge into or pass through the Y road under the action of the cooperative optimization control strategy in the k stagekAll possible cost costs arising from individual vehicle intervals;
Figure BDA0002207808210000024
the method is characterized in that the Y road does not directly participate in vehicle convergence or vehicle confluence or traffic flow of vehicles passing through an intersection, and the cost is consumed due to the fact that the speed of the vehicle is adjusted according to the following safety requirement of the vehicle due to the influence of vehicle convergence or passing of the vehicle on the front vehicle;
an objective function of the upper model:
Figure BDA0002207808210000031
setting the initial condition as f0(s0)=0;
Objective function f of the upper modelk(sk) Representing the cumulative cost consumption of the system vehicles from stage 1 to stage k, f0(s0) 0 denotes that the system cost is 0 in the initial state;
s2, establishing the lower layer model, including:
s2-1, determining a microscopic follow-up model, describing a follow-up state of the vehicle by using the microscopic follow-up model, and predicting an initial track of the vehicle; the following state of the vehicle comprises the speed, acceleration and position of the vehicle;
s2-2, establishing a condition constraint model for judging whether the kth vehicle on the X road at the kth stage can smoothly merge into or pass through the Y road;
s2-3, aiming at various possible situations that occur in the process that a k vehicle on a k-th stage X road converges or passes through a Y road under a mixed traffic scene, simulating a cooperative control strategy set;
s2-4, based on the vehicle initial track predicted in the step S2-1, the condition constraint model established in the step S2-2 sequentially judges that the k-th vehicle on the X-th road of the k stage is converged into or passes through the S-th road on the Y road for converging into or passing throughkThe specific situation occurring in the process of each vehicle interval is determined according to the situation that the k-th vehicle on the X-th road in the k stage can be merged into or passes through the s of the Y roadkSpecific conditions occurring in the process of each vehicle interval respectively make the cooperative control strategy corresponding to the cooperative control strategy set planned in the step S2-3;
s2-5, determining the vehicles which can be optimally controlled in the vehicles participating in the k stage as target vehicles; optimizing the running track of the target vehicle according to the cooperative control strategies made in the step S2-4 respectively, resolving the optimization problem into an optimal control problem of discrete time state constraint, and solving by using a dynamic programming idea to obtain the cooperative optimization control strategy about the target vehicle;
s2-6, calculating the k-th vehicle on the X road or the S which can be converged by the k-th vehicle or pass through the Y road under the action of the cooperative optimization control strategy in the k-th stagekAll possible cost consumptions of one vehicle interval
Figure BDA0002207808210000032
(as one input in the upper model dynamic programming recursive solution process);
s3, solving the double-layer optimization model:
s3-1, solving the upper model, and determining the decision made by each stage of the upper model when the accumulated cost consumption of the system vehicle is the lowest;
s3-2, reversely deducing the decision made at each stage of the upper layer model determined in the step S3-1 to obtain the vehicle optimized track of each single vehicle on the X road of the lower layer model converging into or passing through the Y road;
and S4, solving and obtaining a mixed traffic flow cooperative decision of the double-layer optimization model aiming at system optimization by the step S3, and acting the decision on the system vehicle to control the operation of the system vehicle.
Further, the step S2-3 specifically includes:
according to the k stage X road, the k vehicle is in the s available for the k vehicle to merge or pass through the Y roadkSelecting a particular xth of individual vehicle intervalskThe vehicles gather or pass at intervals, the k-th vehicle on the X road is marked as a vehicle k, and the X-th vehicle is positioned on the Y roadkThe front vehicle of a vehicle interval is marked as a vehicle
Figure BDA0002207808210000041
Rear vehicle as vehicle
Figure BDA0002207808210000042
Further defines a vehicle combination K, representing the combination of vehicles on the X road and the Y road directly participating in the K-th stage, and K ∈ { K }1,K2,K3Vehicle combination
Figure BDA0002207808210000043
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure BDA0002207808210000044
Vehicle k and vehicle
Figure BDA0002207808210000045
Vehicle combination
Figure BDA0002207808210000046
The vehicle participating in the kth stage comprises a vehicle k and a vehicle
Figure BDA0002207808210000047
At this time, the x-th road is positioned on the Y roadkThe interval of each vehicle has no front vehicle; vehicle combination
Figure BDA0002207808210000048
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure BDA0002207808210000049
Vehicle k, now on the Xth road of YkNo rear vehicle participates at the interval of each vehicle;
based on the vehicle track predicted by the microcosmic following model, the vehicle is driven
Figure BDA00022078082100000410
Vehicle k and vehicle
Figure BDA00022078082100000411
The relationship between the two is that the vehicle k can smoothly merge into or pass through the vehicle
Figure BDA00022078082100000412
And a vehicle
Figure BDA00022078082100000413
In the interval between the vehicles, the vehicle k can not smoothly merge into or pass through the vehicle
Figure BDA00022078082100000414
And a vehicle
Figure BDA00022078082100000415
The spacing therebetween; the vehicle k can not smoothly merge into or pass through the vehicle
Figure BDA00022078082100000416
And a vehicle
Figure BDA00022078082100000417
The interval between them is divided into four cases: the first case is denoted as R1 and represents vehicle k and vehicle
Figure BDA00022078082100000418
The distance between the two adjacent layers is too close to meet the constraint condition that the two adjacent layers can be smoothly merged or passed; the second case is denoted as R2 and represents vehicle kAnd a vehicle
Figure BDA00022078082100000419
The distance between the two adjacent layers is too close to meet the constraint condition that the two adjacent layers can be smoothly merged or passed; the third case is denoted as R3 and represents vehicle k and vehicle
Figure BDA00022078082100000420
And with vehicles
Figure BDA00022078082100000421
The process of meeting the basic spacing requirements but merging or passing in is uncomfortable; the fourth case is denoted as R4 and represents vehicle k and vehicle
Figure BDA00022078082100000422
And with vehicles
Figure BDA00022078082100000423
The constraint conditions for smooth import or passing are not met between the two groups;
based on different vehicle combinations, different vehicle type combinations and whether the vehicle k can smoothly merge into or pass through the vehicle
Figure BDA00022078082100000424
And a vehicle
Figure BDA00022078082100000425
The planning of the cooperative control strategy set under different conditions of the interval between the two strategies specifically comprises the following steps:
the collective expression for the vehicle conditions participating in the k-th phase is as follows:
Figure BDA00022078082100000426
Figure BDA00022078082100000427
or 1 or a NaN,
αkeither the number of bits is 0 or 1,
Figure BDA00022078082100000428
or 1 or a NaN,
r k0 or 1 or 2 or 3 or 4,
wherein,
Figure BDA00022078082100000429
αk
Figure BDA00022078082100000430
respectively representing vehicles
Figure BDA00022078082100000431
Vehicle k and vehicle
Figure BDA00022078082100000432
A value 0 indicates that the vehicle type is a non-optimally controllable vehicle, a value 1 indicates that the vehicle type is an optimally controllable vehicle, and a symbol NaN indicates that there is no vehicle; r iskIndicating vehicles
Figure BDA00022078082100000433
Vehicle k and vehicle
Figure BDA00022078082100000434
Relation between r k0 means that the vehicle k can smoothly merge into or pass through the vehicle
Figure BDA00022078082100000435
And a vehicle
Figure BDA00022078082100000436
Interval between rk1 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure BDA00022078082100000437
And a vehicle
Figure BDA00022078082100000438
Case of spacing between R1, R k2 means that the vehicle k is noneLaw of smooth import or passage through vehicles
Figure BDA00022078082100000439
And a vehicle
Figure BDA00022078082100000440
Case of spacing between R2, Rk3 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure BDA0002207808210000052
And a vehicle
Figure BDA0002207808210000053
Case of spacing between R3, R k4 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure BDA0002207808210000054
And a vehicle
Figure BDA0002207808210000055
The case of the interval therebetween R4;
a specific list of vehicle conditions participating in phase k is as follows:
Figure BDA0002207808210000051
and (3) aiming at the condition of the vehicle participating in the k stage, planning a corresponding cooperative control strategy:
defining decision variables of the k-th stage lower layer model at t moment
Figure BDA0002207808210000056
Reflects the acceleration of the vehicle at the moment t in the k-th stage, an
Figure BDA0002207808210000061
Wherein
Figure BDA0002207808210000062
A decision variable set of the lower layer model at the kth stage at the time t;
for vehicle k:
when in use
Figure BDA0002207808210000063
Figure BDA0002207808210000064
Figure BDA0002207808210000065
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure BDA0002207808210000066
Figure BDA0002207808210000067
when in use
Figure BDA0002207808210000068
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure BDA0002207808210000069
when in use
Figure BDA00022078082100000610
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions: v. ofk(t)+uk(t)τ≤ve
For vehicles
Figure BDA00022078082100000611
When in use
Figure BDA00022078082100000612
Figure BDA00022078082100000613
Figure BDA00022078082100000614
In time, the vehicle
Figure BDA00022078082100000615
Decision variables at time t
Figure BDA00022078082100000616
The requirements are as follows:
Figure BDA00022078082100000617
when in use
Figure BDA00022078082100000618
In time, the vehicle
Figure BDA00022078082100000619
Decision variables at time t
Figure BDA00022078082100000620
The requirements are as follows:
Figure BDA00022078082100000621
and is
Figure BDA00022078082100000622
For vehicles
Figure BDA00022078082100000623
When in use
Figure BDA00022078082100000624
Figure BDA00022078082100000625
In time, the vehicle
Figure BDA00022078082100000626
Decision variables at time t
Figure BDA00022078082100000627
The requirements are as follows:
Figure BDA00022078082100000628
when in use
Figure BDA00022078082100000629
In time, the vehicle
Figure BDA00022078082100000630
Decision variables at time t
Figure BDA00022078082100000631
The requirements are as follows:
Figure BDA00022078082100000632
and is
Figure BDA00022078082100000633
Wherein,
Figure BDA00022078082100000634
representing the vehicle condition participating in the k-th phase; τ is the reaction time of the vehicle driving; v. ofk(t)、
Figure BDA00022078082100000635
Respectively show a vehicle k and a vehicle
Figure BDA00022078082100000636
Vehicle with a steering wheel
Figure BDA00022078082100000637
Velocity at time t;
Figure BDA00022078082100000638
Figure BDA00022078082100000639
respectively a vehicle k and a vehicle
Figure BDA00022078082100000641
Vehicle with a steering wheel
Figure BDA00022078082100000640
The safe following acceleration of the vehicle at the time t is predicted according to the microcosmic following model;
Figure BDA00022078082100000642
Figure BDA00022078082100000643
respectively a vehicle k and a vehicle
Figure BDA00022078082100000645
Vehicle with a steering wheel
Figure BDA00022078082100000644
The safe following speed of the vehicle at the moment t + tau is predicted according to the microcosmic following model; v. ofeIs the desired speed; (in the above, Lk(t) represents the relative distance between the vehicle k and its following preceding vehicle at time t,
Figure BDA0002207808210000071
representing the speed of the car k before the car k follows at the time t;
Figure BDA0002207808210000072
indicating vehicle at time t
Figure BDA0002207808210000073
The relative distance between the car and the car before the car is driven,
Figure BDA0002207808210000074
indicating vehicles
Figure BDA0002207808210000075
The speed of the car before the car is followed at the time t;
Figure BDA0002207808210000076
indicating vehicle at time t
Figure BDA0002207808210000077
Relative to the car before it followsThe distance between the first and second electrodes,
Figure BDA0002207808210000078
indicating vehicles
Figure BDA0002207808210000079
The speed of the car ahead at time t. )
If the vehicle in stage k-1
Figure BDA00022078082100000710
The target vehicle is subjected to cooperative optimization control, and the vehicle is taken as the vehicle in the k stage
Figure BDA00022078082100000711
The vehicle processes that will be regarded as non-optimizable control when engaged are expressed in terms of expressions as follows:
Figure BDA00022078082100000712
Figure BDA00022078082100000713
further, the objective function of the lower model is:
Figure BDA00022078082100000714
Figure BDA00022078082100000715
setting an initial cost to
Figure BDA00022078082100000716
Setting an initial state variable to
Figure BDA00022078082100000717
Wherein, the target vehicle in the k stage is taken as the target vehiclei;
Figure BDA00022078082100000718
Indicating the time when the target vehicle i enters the control area in the k-th stage;
Figure BDA00022078082100000719
represents the time when the target vehicle i leaves the control area in the k-th stage; will be provided with
Figure BDA00022078082100000720
To
Figure BDA00022078082100000721
Is divided equally into N segments in discrete time intervals τ', i.e.
Figure BDA00022078082100000722
Defining a control decision time as
Figure BDA00022078082100000723
Objective function of the underlying model
Figure BDA00022078082100000724
Shows that after the target vehicle i is subjected to cooperative optimization control in the k stage, the
Figure BDA00022078082100000725
Cost consumption from time to time t;
Figure BDA00022078082100000726
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA00022078082100000727
reflecting the vehicle position and speed of the target vehicle i in the kth stage at the time t for the state variable of the lower layer model in the kth stage at the time t;
Figure BDA00022078082100000728
reflecting the decision variables of the lower layer model at the kth stage at the t momenti vehicle acceleration at time t;
Figure BDA00022078082100000729
a decision variable set of the lower layer model at the kth stage at the time t; the initial state variable is
Figure BDA00022078082100000730
Indicates that the target vehicle i is in the k-th stage
Figure BDA00022078082100000731
Vehicle position and speed at the moment.
Further, the decision index function of the k-th stage lower layer model at the time t
Figure BDA00022078082100000732
The method specifically comprises the following steps:
when k is equal to 1, the first step is carried out,
Figure BDA00022078082100000733
when k is 2,3, …, n,
Figure BDA00022078082100000734
Figure BDA0002207808210000081
Figure BDA0002207808210000082
Figure BDA0002207808210000083
wherein,
Figure BDA0002207808210000084
representing a decision index function of the lower layer model in the 1 st stage at the t moment;
Figure BDA0002207808210000085
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA0002207808210000086
indicating a time at which the target vehicle of the k-th stage enters the control area;
Figure BDA0002207808210000087
indicating the time when the target vehicle of the k-1 stage leaves the control area;
Figure BDA0002207808210000088
indicating the time at which the target vehicle leaves the control area at the k-th stage.
Further, the vehicle which does not directly participate in the process of vehicle confluence or vehicle passing through the intersection on the Y road in the step 1-2 is marked as a vehicle i', and the cost consumption caused by the vehicle speed adjustment of the vehicle due to the following safety requirement of the vehicle is caused by the influence of the vehicle confluence or passing of the front vehicle
Figure BDA0002207808210000089
The method specifically comprises the following steps:
Figure BDA00022078082100000810
Figure BDA00022078082100000811
wherein,
Figure BDA00022078082100000812
indicating that vehicle i' belongs to a vehicle between the rear vehicle at the k-1 th stage and the front vehicle at the k-th stage on the Y road.
Furthermore, a microscopic traffic flow simulation environment is constructed, and simulation results before and after optimization under different traffic situations are compared.
Compared with the prior art, the method has the beneficial effects that:
the invention provides a mixed traffic flow collaborative optimization control method based on double-layer planning, which is characterized in that an upper layer model of a double-layer optimization model is utilized to ensure that a system for hybrid traffic flow collaborative decision control is optimal by searching a vehicle sequence with optimal system under different traffic scenes (including ramp vehicle convergence on an expressway, vehicle convergence at an intersection and vehicle passing through the intersection), namely searching an optimal sequence with m +1 intervals existing among vehicles converged on an X road or passing through a Y road; the lower layer model of the double-layer optimization model is classified according to cooperative control strategies under different scenes, and based on a corresponding microcosmic car-following model and a condition constraint model for judging whether vehicles on an X road can smoothly merge into or pass through a Y road, vehicle tracks in the process of single vehicle convergence or vehicle confluence or vehicle passing through an intersection are subjected to optimization control, so that the vehicles on the X road can smoothly merge into or pass through the Y road, and the vehicle tracks are optimal. And the vehicle track optimization result obtained by the objective function in the lower model is used as one input in the dynamic programming recursive solving process of the upper model. The upper layer model and the lower layer model jointly ensure the optimal running of the system vehicles, so that the vehicle conflict of the vehicles in the process of converging or passing through the intersection is reduced under the mixed traffic flow environment, and the passing efficiency and the comfort of the vehicles are effectively improved.
The method can be generally applied to various microscopic traffic scenes, including highway ramp vehicle convergence, intersection vehicle convergence, vehicle passing through an intersection and the like, and can be applied only by pertinently changing the upper layer model and the lower layer model, so that the method has certain universal applicability.
A large number of simulation operation results show that the sequence of vehicle convergence or vehicle interval passing can be changed by using the double-layer optimization model, most vehicles can run at a higher speed under the action of a cooperative optimization control strategy, and an acceleration curve is smoother. The throughput of the confluence region segment increases by approximately 18% when the permeability of the autonomous vehicle (i.e. the vehicle that can be optimally controlled) is high. Under the cooperative optimization control mechanism of the mixed traffic flow, the traffic capacity of the road section can be further improved by 10-15%.
The present invention will be described in further detail with reference to the following detailed description and the accompanying drawings, which are not intended to limit the scope of the invention.
Drawings
FIG. 1 is a diagram of the basic architecture of a two-layer optimization model according to an embodiment of the present invention.
Fig. 2 is a schematic view of vehicle convergence of a highway ramp in a mixed traffic flow scene according to an embodiment of the invention.
Fig. 3 is a basic diagram of the permeability of the autonomous vehicle and the road section passing capacity under the GIPPS following model according to the embodiment of the invention.
Fig. 4 is a microscopic trace diagram of 23 vehicles on the ramp passing through the confluent region section when the permeability of the autonomous vehicle is 90% and is not controlled by the double-layer optimization model according to the embodiment of the invention.
Fig. 5 is a micro-motion trace diagram of 23 vehicles on a ramp passing through a confluence region segment when the penetration rate of the autonomous vehicle is 90% and controlled by a double-layer optimization model according to an embodiment of the invention.
Detailed Description
The following description of the embodiments of the present invention will be made with reference to the accompanying drawings.
A mixed traffic flow collaborative optimization control method based on double-layer planning is characterized in that a double-layer optimization model based on dynamic planning recursion is adopted for carrying out mixed traffic flow collaborative decision control, and the method is suitable for different traffic scenes in the mixed traffic flow, wherein the different traffic scenes comprise ramp vehicle convergence on a highway, vehicle confluence at an intersection and vehicle passing through the intersection; the double-layer optimization model comprises an upper layer model and a lower layer model; the upper model is a vehicle sequencing problem solved by dynamic programming recursion; the lower model is a single vehicle trajectory optimization problem in different traffic scenarios solved recursively with dynamic programming. Fig. 1 is a diagram showing a basic architecture of a two-layer optimization model according to an embodiment of the present invention.
Example one
The present embodiment describes the method of the present invention by taking highway ramp vehicle convergence as an example.
As shown in fig. 2, the schematic view of vehicle convergence on a ramp of an expressway under a mixed traffic flow scene is shown, where a main road (i.e., a Y road) and a ramp (i.e., an X road) are both one-way roads, a plurality of traffic flows randomly arranged between an automatically-driven vehicle a (i.e., an optimally-controlled vehicle) and a human-driven vehicle H (i.e., a non-optimally-controlled vehicle) are provided on the main road, and a plurality of traffic flows randomly arranged between the automatically-driven vehicle a and the human-driven vehicle H are provided on the ramp and need to converge between the traffic flows on the main road through a convergence region. Assume that the autonomous vehicle a and the human-driven vehicle H on the main road and the ramp both follow the microscopic follow-up model.
Description of the upper layer model: the process of converging vehicles on a ramp of an expressway, namely vehicles on the ramp are sequentially converged into a main road, starting time of an initial stage from the moment when a first ramp vehicle is ready to be converged, and then convergence of a next ramp vehicle is the next stage. The first ramp vehicle faces the selection of a plurality of spaces into which vehicles can be converged on the main road, and the number of the spaces into which vehicles can be converged on the main road is correspondingly changed along with the successful convergence of the vehicles on each ramp. The upper dynamic planning model is mainly used for finding the optimal sequence of the interval of vehicles on the ramp and merging into the main road.
The lower model describes: in the junction area of the ramps of the expressway, when vehicles normally run, the traffic flow on the same lane follows the traffic flow according to the microcosmic car following model, namely, the vehicles on the main road
Figure BDA0002207808210000101
And the vehicles k on the ramp follow the microscopic follow-up model to run. Suppose that when vehicle k enters the merge area on a ramp, vehicle k may observe a conflicting vehicle on the arterial road
Figure BDA0002207808210000102
And
Figure BDA0002207808210000103
at this time, according to the vehicle k and the vehicle
Figure BDA0002207808210000104
Whether the relative distance between the two meets the condition of conflux or not; if yes, the vehicle k successfully enters the main road; if not, judging the reason why the vehicle k cannot be converged successfully, optimizing the vehicle running track by applying a corresponding cooperative control strategy and combining with dynamic planning, and calculating corresponding consumption.
The method comprises the following steps:
s1, establishing the upper layer model, including:
s1-1, determining the division stage, the state variable and the decision variable of the upper layer model, which are as follows:
the upper layer model is divided: assuming that the main road and the ramp are both one-way roads, n vehicles on the ramp need to sequentially merge into m +1 intervals existing among the vehicles on the main road, the behavior that each single vehicle on the ramp merges into the main road is represented as a stage, and the behavior that the kth vehicle on the ramp merges into the main road is recorded as the kth stage, wherein k is 1,2,3, …, n.
State variables of the upper model: the number of vehicle intervals for leading the kth vehicle on the ramp to merge into the main road in the kth stage is skAnd (4) showing.
Decision variables of the upper model: the decision made at each stage represents that the k-th vehicle on the ramp at the k-th stage can merge into the main road skSelecting a particular xth of individual vehicle intervalskThe vehicles are alternately merged.
S1-2, determining a state transition equation, a cost function and an objective function of the upper layer model, which are as follows:
the state transition equation of the upper model is as follows:
Figure BDA0002207808210000105
setting the initial condition as s0=m+1;
The state transition equation of the upper model indicates that s is equal to 1 when kk=s1Is shown inIn the stage 1, the number of vehicle intervals for the 1 st vehicle on the ramp to flow into the main road is m + 1; when k is 2,3, …, n, the x-th vehicle on the ramp is selected at the k-1 stagek-1The interval of each vehicle is used as a state variable s after being converged into the main roadkA change in (c); s0M +1 represents that the number of vehicle intervals for vehicles on the ramp to merge into the main road in the initial state is m + 1.
Cost function of the upper model:
Figure BDA0002207808210000111
cost function D of the upper modelk(sk,xk) Indicating a phase index function required for making a decision at the kth phase, wherein
Figure BDA0002207808210000112
S for showing that the k-th vehicle on the ramp is merged into the main road for the k-th vehicle to merge into the main road under the action of the cooperative optimization control strategy in the k-th stagekAll possible cost costs arising from individual vehicle intervals;
Figure BDA0002207808210000113
the method is characterized in that the main road does not directly participate in the traffic flow in the vehicle converging process, and the cost is consumed due to the fact that the speed of the vehicle is adjusted according to the following safety requirement of the vehicle because the front vehicle is influenced by the vehicle converging;
as described above
Figure BDA0002207808210000114
The mathematical expression of (a) is as follows:
Figure BDA0002207808210000115
Figure BDA0002207808210000116
wherein,
Figure BDA0002207808210000117
it indicates that the vehicle i' belongs to a vehicle between the rear vehicle of the k-1 stage and the front vehicle of the k stage on the main road.
An objective function of the upper model:
Figure BDA0002207808210000118
setting the initial condition as f0(s0)=0;
Objective function f of the upper modelk(sk) Representing the cumulative cost consumption of the system vehicles from stage 1 to stage k, f0(s0) 0 indicates that the system cost is 0 in the initial state.
S2, establishing the lower layer model, including:
s2-1, determining a microscopic follow-up model, describing a follow-up state of the vehicle by using the microscopic follow-up model, and predicting an initial track of the vehicle; the following state of the vehicle includes a speed, an acceleration, and a position of the vehicle.
The microcosmic car-following model can be selected from specific microcosmic car-following models, such as a GIPPS car-following model, an IDM/EIDM car-following model and the like, and the specific model can be selected to be judged according to the simulation condition.
The GIPPS following model was selected for this example as follows:
v(t+τ)=min(ve,v(t)+aτ,vsafe),
Figure BDA0002207808210000119
l(t+τ)=l(t)-v(t)τ-0.5u(t)τ2
the velocity of the vehicle at time t + τ is denoted by v (t + τ), where veV (t) + a τ represents the velocity obtained by acceleration at a constant acceleration a, v (t) + a τ being the desired velocitysafeIndicating a safe speed. τ is the reaction time of the vehicle driving; b is a constant deceleration; v (t) represents the speed of the vehicle at time t; u (t) represents the acceleration of the vehicle at time t; l (t) watchShowing the position of the vehicle at time t; v. oflead(t) represents the speed of the preceding vehicle followed by the vehicle at time t; llead(t) represents a position of a preceding vehicle followed by the vehicle at time t; laIndicating the vehicle body length; l is0Indicating a minimum safe following distance between the vehicle and a preceding vehicle to which the vehicle is following.
S2-2, establishing a condition constraint model for judging whether the kth vehicle on the kth stage ramp can smoothly merge into the main road.
Suppose that the k-th vehicle on the ramp (denoted as vehicle k) is in the merge region (i.e., z)1<lk(t)≤z0) Which will merge into two vehicles of the continuous flow of the main road
Figure BDA0002207808210000121
And
Figure BDA0002207808210000122
wherein the vehicle
Figure BDA0002207808210000123
Is a front vehicle or a vehicle
Figure BDA0002207808210000124
For the rear vehicle, the following conditional constraint model (i.e., the convergence model) is used to judge whether the vehicle k can smoothly converge into the main road:
Figure BDA0002207808210000125
Figure BDA0002207808210000126
Figure BDA0002207808210000127
u abovek(t) represents the effect of the confluence, reflects the comfort level during the confluence, and is the distance between vehicles during the confluence, the confluent vehicle k on the ramp during the confluence and the rear vehicle of the main road
Figure BDA0002207808210000128
Is calibrated. Wherein,
Figure BDA0002207808210000129
what is shown is the utility of the bus action when unconstrained. laThe length of the vehicle body is long,
Figure BDA00022078082100001210
for a minimum safe distance of autonomous vehicle convergence,
Figure BDA00022078082100001211
minimum safe distance for human-driven vehicles to converge. During confluence, a confluent vehicle k on the ramp actually follows the front vehicle on the main road
Figure BDA00022078082100001212
Rear cars running on the main road
Figure BDA00022078082100001213
The actual following bus vehicles k run, and their accelerations can be calculated from the vehicle-following model.
Figure BDA00022078082100001214
An absolute value representing the acceleration of the vehicle k;
Figure BDA00022078082100001215
rear vehicle on main road
Figure BDA00022078082100001216
Absolute value of acceleration of (a); bsafeIndicating the maximum allowed deceleration. PhiAFor autonomous vehicle set, phiHSet of vehicles for human driving η1And η2Respectively representing a safety factor and a polite factor, the safety factor η1Is a constant, polite coefficient η2In a piecewise continuous form, vthIs a given speed threshold, veTo desired speed, β1And β2Is a constant. By means of Ik(t + τ) represents a convergence decision, a value of 0 indicates that the vehicle k cannot smoothly complete convergence at time t + τ, and a value of 1 indicates that the vehicle k can smoothly complete convergence at time t + τ. When making convergence decision IkThe value of (t + tau) is 0, namely when the vehicle k cannot smoothly complete confluence, the vehicle k on the ramp can use the end of the ramp as a stopped virtual front vehicle, continuously decelerates or even stops to wait by following a microscopic following model until a vehicle interval meeting confluence appears on the main road, and the vehicle k is converged into the main road.
S2-3, aiming at various possible situations that occur in the process that the kth vehicle on the kth stage ramp is converged into the main road under the mixed traffic scene, simulating a cooperative control strategy set. The method specifically comprises the following steps:
according to the k vehicle on the k stage ramp, the k vehicle can be converged into the main roadkSelecting a particular xth of individual vehicle intervalskThe vehicles are converged at intervals, the kth vehicle on the ramp is marked as a vehicle k, and the xth vehicle is positioned on the main roadkThe front vehicle of a vehicle interval is marked as a vehicle
Figure BDA0002207808210000131
Rear vehicle as vehicle
Figure BDA0002207808210000132
Further defining a vehicle combination K representing a combination of vehicles on the ramp and the main road directly participating in the K-th stage, and K ∈ { K }1,K2,k3Vehicle combination
Figure BDA0002207808210000133
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure BDA0002207808210000134
Vehicle k and vehicle
Figure BDA0002207808210000135
Vehicle combination
Figure BDA0002207808210000136
The vehicle participating in the kth stage comprises a vehicle k and a vehicle
Figure BDA0002207808210000137
Now located on the x-th main roadkThe interval of each vehicle has no front vehicle; vehicle combination
Figure BDA0002207808210000138
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure BDA0002207808210000139
Vehicle k, now on the x-th roadkNo rear vehicle participates at the interval of each vehicle;
based on the vehicle track predicted by the microcosmic following model, the vehicle is driven
Figure BDA00022078082100001310
Vehicle k and vehicle
Figure BDA00022078082100001311
The relationship between the two is that the vehicle k can smoothly merge into the vehicle
Figure BDA00022078082100001312
And a vehicle
Figure BDA00022078082100001313
The vehicle k can not smoothly merge into the vehicle
Figure BDA00022078082100001314
And a vehicle
Figure BDA00022078082100001315
The spacing therebetween; the vehicle k can not smoothly merge into the vehicle
Figure BDA00022078082100001316
And a vehicle
Figure BDA00022078082100001317
BetweenThe interval of (c) is divided into four cases: the first case is denoted as R1 and represents vehicle k and vehicle
Figure BDA00022078082100001318
The distance between the two adjacent beams is too close to meet the constraint condition of smooth import; the second case is denoted as R2 and represents vehicle k and vehicle
Figure BDA00022078082100001319
The distance between the two adjacent beams is too close to meet the constraint condition of smooth import; the third case is denoted as R3 and represents vehicle k and vehicle
Figure BDA00022078082100001320
And with vehicles
Figure BDA00022078082100001321
The process of meeting the basic spacing requirement but merging is not comfortable; the fourth case is denoted as R4 and represents vehicle k and vehicle
Figure BDA00022078082100001322
And with vehicles
Figure BDA00022078082100001323
The constraint conditions which can be smoothly merged are not met;
based on different vehicle combinations, different vehicle type combinations and whether the vehicle k can smoothly merge into the vehicle
Figure BDA00022078082100001324
And a vehicle
Figure BDA00022078082100001325
The planning of the cooperative control strategy set under different conditions of the interval between the two strategies specifically comprises the following steps:
the collective expression for the vehicle conditions participating in the k-th phase is as follows:
Figure BDA00022078082100001326
Figure BDA00022078082100001327
or 1 or a NaN,
αkeither the number of bits is 0 or 1,
Figure BDA00022078082100001328
or 1 or a NaN,
r k0 or 1 or 2 or 3 or 4,
wherein,
Figure BDA00022078082100001329
αk
Figure BDA00022078082100001330
respectively representing vehicles
Figure BDA00022078082100001331
Vehicle k and vehicle
Figure BDA00022078082100001332
A value
0 indicates that the vehicle type is a human-driven vehicle, a value 1 indicates that the vehicle type is an autonomous vehicle, and a symbol NaN indicates that there is no vehicle; r iskIndicating vehicles
Figure BDA00022078082100001333
Vehicle k and vehicle
Figure BDA00022078082100001334
Relation between r k0 means that the vehicle k can smoothly merge into the vehicle
Figure BDA00022078082100001335
And a vehicle
Figure BDA00022078082100001336
Interval between rk1 means that the vehicle k cannot smoothly merge into the vehicle
Figure BDA00022078082100001337
And a vehicle
Figure BDA00022078082100001338
Case of spacing between R1, R k2 means that the vehicle k cannot smoothly merge into the vehicle
Figure BDA0002207808210000142
And a vehicle
Figure BDA0002207808210000143
Case of spacing between R2, Rk3 means that the vehicle k cannot smoothly merge into the vehicle
Figure BDA0002207808210000144
And a vehicle
Figure BDA0002207808210000145
Case of spacing between R3, R k4 means that the vehicle k cannot smoothly merge into the vehicle
Figure BDA0002207808210000146
And a vehicle
Figure BDA0002207808210000147
The case of the interval therebetween R4;
a specific list of vehicle conditions participating in phase k is as follows:
Figure BDA0002207808210000141
and (3) aiming at the condition of the vehicle participating in the k stage, planning a corresponding cooperative control strategy:
defining decision variables of the k-th stage lower layer model at t moment
Figure BDA0002207808210000148
Reflects the acceleration of the vehicle at the moment t in the k-th stage, an
Figure BDA0002207808210000151
Wherein
Figure BDA0002207808210000152
A decision variable set of the lower layer model at the kth stage at the time t;
for vehicle k:
when in use
Figure BDA0002207808210000153
Figure BDA0002207808210000154
Figure BDA0002207808210000155
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure BDA0002207808210000156
Figure BDA0002207808210000157
when in use
Figure BDA0002207808210000158
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure BDA0002207808210000159
when in use
Figure BDA00022078082100001510
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions: v. ofk(t)+uk(t)τ≤ve
For vehicles
Figure BDA00022078082100001511
When in use
Figure BDA00022078082100001512
Figure BDA00022078082100001513
Figure BDA00022078082100001514
In time, the vehicle
Figure BDA00022078082100001515
Decision variables at time t
Figure BDA00022078082100001516
The requirements are as follows:
Figure BDA00022078082100001517
when in use
Figure BDA00022078082100001518
In time, the vehicle
Figure BDA00022078082100001520
Decision variables at time t
Figure BDA00022078082100001519
The requirements are as follows:
Figure BDA00022078082100001521
and is
Figure BDA00022078082100001522
For vehicles
Figure BDA00022078082100001523
When in use
Figure BDA00022078082100001524
Figure BDA00022078082100001525
In time, the vehicle
Figure BDA00022078082100001526
Decision variables at time t
Figure BDA00022078082100001527
The requirements are as follows:
Figure BDA00022078082100001528
when in use
Figure BDA00022078082100001529
In time, the vehicle
Figure BDA00022078082100001530
Decision variables at time t
Figure BDA00022078082100001531
The requirements are as follows:
Figure BDA00022078082100001532
and is
Figure BDA00022078082100001533
Wherein,
Figure BDA00022078082100001534
representing the vehicle condition participating in the k-th phase; τ is the reaction time of the vehicle driving; v. ofk(t)、
Figure BDA00022078082100001535
Respectively show a vehicle k and a vehicle
Figure BDA00022078082100001536
Vehicle with a steering wheel
Figure BDA00022078082100001537
Velocity at time t;
Figure BDA00022078082100001538
Figure BDA00022078082100001539
are respectivelyVehicle k and vehicle
Figure BDA00022078082100001541
Vehicle with a steering wheel
Figure BDA00022078082100001540
The safe following acceleration of the vehicle at the time t is predicted according to the microcosmic following model;
Figure BDA00022078082100001542
Figure BDA00022078082100001543
respectively a vehicle k and a vehicle
Figure BDA00022078082100001545
Vehicle with a steering wheel
Figure BDA00022078082100001544
The safe following speed of the vehicle at the moment t + tau is predicted according to the microcosmic following model; v. ofeIs the desired speed.
If the vehicle in stage k-1
Figure BDA0002207808210000161
The target vehicle is subjected to cooperative optimization control, and the vehicle is taken as the vehicle in the k stage
Figure BDA0002207808210000162
The process of taking part in the process of being treated as human driving a vehicle is expressed by the following expression:
Figure BDA0002207808210000163
Figure BDA0002207808210000164
s2-4, based on the vehicle initial track predicted in the step S2-1, sequentially judging that the k-th vehicle on the ramp of the k-th stage merges into the main road through the condition constraint model established in the step S2-2S on which the kth vehicle can merge into the trunk roadkThe specific situation of each vehicle interval is determined according to the fact that the k-th vehicle on the kth stage ramp can be converged into the main road on the main road for the k-th vehicle to be converged into the main roadkThe specific situation occurring in the process of each vehicle interval respectively makes the cooperative control strategy corresponding to the cooperative control strategy set planned in the step S2-3.
S2-5, determining the vehicles which can be optimally controlled in the vehicles participating in the k stage as target vehicles; and optimizing the running track of the target vehicle according to the cooperative control strategies made in the step S2-4 respectively, resolving the optimization problem into an optimal control problem of discrete time state constraint, and solving by using a dynamic planning idea to obtain the cooperative optimization control strategy about the target vehicle.
S2-6, calculating S for the k-th vehicle on the ramp to merge into the main road under the action of the cooperative optimization control strategy in the k-th stagekAll possible cost consumptions of one vehicle interval
Figure BDA0002207808210000165
The objective function of the lower model is:
Figure BDA0002207808210000166
Figure BDA0002207808210000167
setting an initial cost to
Figure BDA0002207808210000168
Setting an initial state variable to
Figure BDA0002207808210000169
Wherein, the target vehicle in the k stage is recorded as a target vehicle i;
Figure BDA00022078082100001610
indicating the time when the target vehicle i enters the control area in the k-th stage;
Figure BDA00022078082100001611
represents the time when the target vehicle i leaves the control area in the k-th stage; will be provided with
Figure BDA00022078082100001612
To
Figure BDA00022078082100001613
Is divided equally into N segments in discrete time intervals τ', i.e.
Figure BDA00022078082100001614
Defining a control decision time as
Figure BDA00022078082100001615
Objective function of the underlying model
Figure BDA00022078082100001616
Shows that after the target vehicle i is subjected to cooperative optimization control in the k stage, the
Figure BDA00022078082100001617
Cost consumption from time to time t;
Figure BDA00022078082100001618
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA00022078082100001619
reflecting the vehicle position and speed of the target vehicle i in the kth stage at the time t for the state variable of the lower layer model in the kth stage at the time t;
Figure BDA00022078082100001620
reflecting the vehicle acceleration of the target vehicle i at the moment t in the kth stage for the decision variable of the lower layer model at the moment t in the kth stage;
Figure BDA00022078082100001621
a decision variable set of the lower layer model at the kth stage at the time t; the initial state variable is
Figure BDA00022078082100001622
Indicates that the target vehicle i is in the k-th stage
Figure BDA00022078082100001623
Vehicle position and speed at the moment.
Decision index function of the k-th stage lower layer model at t moment
Figure BDA00022078082100001624
The method specifically comprises the following steps:
when k is equal to 1, the first step is carried out,
Figure BDA0002207808210000172
when k is 2,3, …, n,
Figure BDA0002207808210000173
Figure BDA0002207808210000174
Figure BDA0002207808210000175
Figure BDA0002207808210000176
wherein,
Figure BDA0002207808210000177
representing a decision index function of the lower layer model in the 1 st stage at the t moment;
Figure BDA0002207808210000178
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA0002207808210000179
indicating a time at which the target vehicle of the k-th stage enters the control area;
Figure BDA00022078082100001710
indicating the time when the target vehicle of the k-1 stage leaves the control area;
Figure BDA00022078082100001711
indicating the time at which the target vehicle leaves the control area at the k-th stage.
S3, solving the double-layer optimization model:
s3-1, solving the upper model, and determining the decision made by each stage of the upper model when the accumulated cost consumption of the system vehicle is the lowest;
s3-2, and reversely deducing the decision made at each stage of the upper layer model determined in the step S3-1 to obtain the vehicle optimized track of each single vehicle on the lower layer model ramp merging into the main road.
And S4, solving and obtaining a mixed traffic flow cooperative decision of the double-layer optimization model aiming at system optimization by the step S3, and acting the decision on the system vehicle to control the operation of the system vehicle.
The method establishes a microscopic traffic flow simulation environment through programming, and the parameters and values of the simulation experiment environment are realized by computer programming as shown in the following table:
Figure BDA0002207808210000171
Figure BDA0002207808210000181
fig. 3 is a basic graph of autonomous vehicle permeability versus road segment traffic capacity for the present embodiment under the GIPPS-followed model, showing the flow-density relationship for different autonomous vehicle permeabilities observed just downstream of the confluence region segment. Wherein graph a represents a comparison of the flow-density relationship for a 100% autonomous vehicle and a 100% human-driven vehicle; the graph B, C, D reflects the flow-density relationship for autonomous vehicle permeabilities of 30%, 50%, and 70%, respectively. From fig. 3 it can be found that: when the permeability of the automatic driving vehicle is 0, namely all vehicles are human driving vehicles, the road section traffic capacity is about 2244 veh/h; when the permeability of the automatic driving vehicle is 50%, namely the ratio of the human driving vehicle to the automatic driving vehicle is 1:1, compared with the permeability of the automatic driving vehicle is 0, the road section traffic capacity is improved by about 11.8%, and reaches about 2508 veh/h; when the permeability of the automatic driving vehicle is 100 percent, namely all the automatic driving vehicles are automatic driving vehicles, compared with the permeability of the automatic driving vehicles of 0, the road section traffic capacity is improved by about 17.7 percent and reaches 2640 veh/h. Experiments using other microscopic follow models also showed similar results. It follows from this that: if the permeability of the autonomous vehicle increases, the traffic capacity of the road segment is improved.
Fig. 4 is a microscopic movement trace diagram of 23 vehicles on the ramp passing through the confluence region segment when the penetration rate of the autonomous vehicle is 90% and is not controlled by the double-layer optimization model in the embodiment. Wherein, the graph (a) is a velocity-time relation graph, and the graph (b) is an acceleration-time relation graph.
Fig. 5 is a microscopic movement trace diagram of 23 vehicles on the ramp passing through the confluence region segment when the penetration rate of the automatic driving vehicle is 90% and controlled by the double-layer optimization model. Wherein, the graph (a) is a velocity-time relation graph, and the graph (b) is an acceleration-time relation graph.
By comparing fig. 4 and fig. 5, it can be found that: by applying the double-layer optimization model, the vehicle convergence sequence can be changed, so that the global optimum is ensured. Under the action of a cooperative optimization control strategy, most vehicles can run at a higher speed, and an acceleration curve is smoother. When the permeability of the automatic driving vehicle is higher, the traffic capacity of the road section can be further improved by 10-15% under the cooperative optimization control mechanism of the mixed traffic flow.
Example two
This embodiment will explain the method of the present invention by taking intersection vehicle merging as an example.
Merging vehicles at the intersection: the X road and the Y road are two crossed one-way roads, a plurality of automatic driving vehicles (vehicles which can be optimally controlled) and human driving vehicles (vehicles which cannot be optimally controlled) are arranged on the Y road to form irregularly arranged traffic flows, and a plurality of automatic driving vehicles and human driving vehicles which form irregularly arranged traffic flows need to turn at the intersection and converge into intervals among the traffic flows on the Y road. It is assumed that the autonomous vehicle and the human-driven vehicle on both the X-road and the Y-road follow the microscopic follow-up model.
Description of the upper layer model: the intersection vehicle confluence, namely the process that vehicles to be confluent turn and converge into straight traffic flow at the intersection in sequence, is started when the first vehicle to be confluent is ready to converge into straight traffic flow, and then each vehicle to be confluent turns and converges into straight traffic flow in one stage. The first vehicle to be merged is selected to face a plurality of intervals of vehicles capable of merging, and the number of the intervals of the vehicles capable of merging is correspondingly changed along with successful merging of the vehicles to be merged. The upper dynamic planning model is mainly used for searching an optimal sequence of vehicle intervals for converging vehicles to be converged into a straight traffic flow.
The lower model describes: under the mixed traffic environment of urban intersections, when vehicles normally run, the traffic flow on the same lane follows the vehicle according to the microcosmic car following model, namely, the vehicles on the Y road
Figure BDA0002207808210000191
And the vehicle k on the X road runs according to the microcosmic following model. Suppose that when vehicle k on the X road enters the collision area, vehicle k can observe a colliding vehicle on the Y road
Figure BDA0002207808210000192
And
Figure BDA0002207808210000193
at this time, according to the vehicle k and the vehicle
Figure BDA0002207808210000194
Whether the relative distance between the two meets the condition of confluent flow or not; if the vehicle k meets the requirement, the vehicle k successfully enters the straight-going traffic flow; if not, judging the reason why the vehicle k cannot successfully converge into the straight traffic flow, optimizing the running track of the vehicle by applying a corresponding cooperative control strategy and combining with dynamic planning, and calculating corresponding consumption.
The method comprises the following steps:
s1, establishing the upper layer model, including:
s1-1, determining the division stage, the state variable and the decision variable of the upper layer model, which are as follows:
the upper layer model is divided: assuming that an X road and a Y road are two one-way roads with intersections, n vehicles on the X road need to sequentially converge into m +1 intervals existing among the vehicles on the Y road, the behavior that each one-way vehicle on the X road converges into the Y road is represented as a stage, and the behavior that the k-th vehicle on the X road converges into the Y road is represented as a k-th stage, wherein k is 1,2,3, …, n.
State variables of the upper model: the number of the vehicle intervals for the k-th vehicle on the X road to merge into the Y road in the k stage is skAnd (4) showing.
Decision variables of the upper model: the decision made in each stage represents that the k vehicle on the X road in the k stage can be converged into the s road of the Y roadkSelecting a particular xth of individual vehicle intervalskThe vehicles are alternately merged.
S1-2, determining a state transition equation, a cost function and an objective function of the upper layer model, which are as follows:
the state transition equation of the upper model is as follows:
Figure BDA0002207808210000201
setting the initial condition as s0=m+1;
The state transition equation of the upper model indicates that s is equal to 1 when kk=s1The number of vehicle intervals for enabling the 1 st vehicle on the X road to merge into the Y road in the 1 st stage is m + 1; when k is 2,3, …, n, the (k-1) th vehicle on the X road selects the (X) th vehicle at the (k-1) th stagek-1The interval of each vehicle is used as a state variable s after converging into the Y roadkA change in (c); s0M +1 indicates that the number of vehicle intervals at which vehicles on the X road can merge into the Y road in the initial state is m + 1.
Cost function of the upper model:
Figure BDA0002207808210000202
cost function D of the upper modelk(sk,xk) Indicating a phase index function required for making a decision at the kth phase, wherein
Figure BDA0002207808210000203
S for showing that the k-th vehicle on the X road is converged into the road available for converging the k-th vehicle into the Y road under the action of the cooperative optimization control strategy in the k stagekAll possible cost costs arising from individual vehicle intervals;
Figure BDA0002207808210000204
the method is characterized in that the vehicle flow which directly participates in the vehicle confluence process on the Y road is not generated, and the cost is consumed due to the fact that the speed of the vehicle is adjusted by the vehicle due to the fact that the vehicle ahead is influenced by vehicle confluence;
as described above
Figure BDA0002207808210000205
The mathematical expression of (a) is as follows:
Figure BDA0002207808210000206
Figure BDA0002207808210000207
wherein,
Figure BDA0002207808210000208
indicating vehicle iAnd the vehicles belong to the vehicles between the rear vehicle at the k-1 stage and the front vehicle at the k stage on the Y road.
An objective function of the upper model:
Figure BDA0002207808210000209
setting the initial condition as f0(s0)=0;
Objective function f of the upper modelk(sk) Representing the cumulative cost consumption of the system vehicles from stage 1 to stage k, f0(s0) 0 indicates that the system cost is 0 in the initial state.
S2, establishing the lower layer model, including:
s2-1, determining a microscopic follow-up model, describing a follow-up state of the vehicle by using the microscopic follow-up model, and predicting an initial track of the vehicle; the following state of the vehicle includes a speed, an acceleration, and a position of the vehicle.
The microcosmic car-following model can be selected from specific microcosmic car-following models, such as a GIPPS car-following model, an IDM/EIDM car-following model and the like, and the specific model can be selected to be judged according to the simulation condition.
And S2-2, establishing a condition constraint model for judging whether the k vehicle on the X road at the k stage can smoothly merge into the Y road.
Based on a microcosmic car-following model, acceleration constraint, distance constraint and safety constraint are added to establish a condition constraint model which accords with the vehicle confluence characteristic of the intersection by combining the road geometric characteristics of the intersection.
S2-3, aiming at various possible situations that occur in the process that the k vehicle on the X road at the k stage under the mixed traffic scene converges into the Y road, simulating a cooperative control strategy set.
S2-4, based on the vehicle initial track predicted in the step S2-1, sequentially judging that the k-th vehicle on the X road at the k stage converges to the Y road for converging the k-th vehicle on the Y road by the condition constraint model established in the step S2-2S into Y roadkThe specific situation of each vehicle interval is determined according to the situation that the k-th vehicle on the X-th road in the k stage can be supplied to the Y road on the Y road for the k-th vehicle to enter the Y roadkThe specific situation occurring in the process of each vehicle interval respectively makes the cooperative control strategy corresponding to the cooperative control strategy set planned in the step S2-3.
S2-5, determining the vehicles which can be optimally controlled in the vehicles participating in the k stage as target vehicles; and optimizing the running track of the target vehicle according to the cooperative control strategies made in the step S2-4 respectively, resolving the optimization problem into an optimal control problem of discrete time state constraint, and solving by using a dynamic planning idea to obtain the cooperative optimization control strategy about the target vehicle.
S2-6, calculating S for the k-th vehicle on the X road to merge into the Y road under the action of the cooperative optimization control strategy in the k stagekAll possible cost consumptions of one vehicle interval
Figure BDA0002207808210000211
The objective function of the lower model is:
Figure BDA0002207808210000212
Figure BDA0002207808210000213
setting an initial cost to
Figure BDA0002207808210000214
Setting an initial state variable to
Figure BDA0002207808210000215
Wherein, the target vehicle in the k stage is recorded as a target vehicle i;
Figure BDA0002207808210000216
indicating the time when the target vehicle i enters the control area in the k-th stage;
Figure BDA0002207808210000217
represents the time when the target vehicle i leaves the control area in the k-th stage; will be provided with
Figure BDA0002207808210000218
To
Figure BDA0002207808210000219
Is divided equally into N segments in discrete time intervals τ', i.e.
Figure BDA00022078082100002110
Defining a control decision time as
Figure BDA00022078082100002111
Objective function of the underlying model
Figure BDA00022078082100002112
Shows that after the target vehicle i is subjected to cooperative optimization control in the k stage, the
Figure BDA00022078082100002113
Cost consumption from time to time t;
Figure BDA00022078082100002114
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA00022078082100002115
reflecting the vehicle position and speed of the target vehicle i in the kth stage at the time t for the state variable of the lower layer model in the kth stage at the time t;
Figure BDA00022078082100002116
reflecting the vehicle acceleration of the target vehicle i at the moment t in the kth stage for the decision variable of the lower layer model at the moment t in the kth stage;
Figure BDA00022078082100002117
a decision variable set of the lower layer model at the kth stage at the time t; the initial state variable is
Figure BDA00022078082100002118
Indicates that the target vehicle i is in the k-th stage
Figure BDA00022078082100002119
Vehicle position and speed at the moment.
Decision index function of the k-th stage lower layer model at t moment
Figure BDA0002207808210000221
The method specifically comprises the following steps:
when k is equal to 1, the first step is carried out,
Figure BDA0002207808210000222
when k is 2,3, …, n,
Figure BDA0002207808210000223
Figure BDA0002207808210000224
Figure BDA0002207808210000225
Figure BDA0002207808210000226
wherein,
Figure BDA0002207808210000227
representing a decision index function of the lower layer model in the 1 st stage at the t moment;
Figure BDA0002207808210000228
is shown asA decision index function of the lower layer model at the moment t in the stage k;
Figure BDA0002207808210000229
indicating a time at which the target vehicle of the k-th stage enters the control area;
Figure BDA00022078082100002210
indicating the time when the target vehicle of the k-1 stage leaves the control area;
Figure BDA00022078082100002211
indicating the time at which the target vehicle leaves the control area at the k-th stage.
S3, solving the double-layer optimization model:
s3-1, solving the upper model, and determining the decision made by each stage of the upper model when the accumulated cost consumption of the system vehicle is the lowest;
s3-2, and reversely deducing the decision made at each stage of the upper layer model determined in the step S3-1 to obtain the vehicle optimized track of each single vehicle on the X road of the lower layer model converging into the Y road.
And S4, solving and obtaining a mixed traffic flow cooperative decision of the double-layer optimization model aiming at system optimization by the step S3, and acting the decision on the system vehicle to control the operation of the system vehicle.
EXAMPLE III
The present embodiment describes the method of the present invention by taking the vehicle passing through an intersection as an example.
Vehicle passing through the intersection: the X road and the Y road are two crossed one-way roads, a plurality of automatic driving vehicles (vehicles capable of being controlled in an optimized mode) and human driving vehicles (vehicles incapable of being controlled in an optimized mode) are arranged on the Y road to form irregularly arranged traffic flows, and the plurality of automatic driving vehicles and the human driving vehicles which form the irregularly arranged traffic flows need to pass through the intersection and sequentially pass through intervals among the traffic flows on the Y road. It is assumed that the autonomous vehicle and the human-driven vehicle on both the X-road and the Y-road follow the microscopic follow-up model.
Description of the upper layer model: the process that vehicles passing through the intersection, namely vehicles to pass through the intersection on an X road pass through vehicle intervals on a Y road in sequence at the intersection is started when a first vehicle to pass through is ready to pass through, and then each vehicle to pass through the intersection is a stage, wherein the first vehicle to pass through faces the selection of a plurality of passing vehicle intervals on the Y road, and the number of the passing vehicle intervals on the Y road is correspondingly changed along with the successful passing of each vehicle to pass through the intersection. The upper dynamic planning model is mainly used for finding the optimal sequence of the vehicle intervals of the vehicles to pass through the intersection on the Y road.
The lower model describes: under the mixed traffic environment of urban intersections, when vehicles normally run, the traffic flow on the same lane follows the vehicle according to the microcosmic car following model, namely, the vehicles on the Y road
Figure BDA0002207808210000231
And the vehicle k on the X road runs according to the microcosmic following model. Suppose that when vehicle k on the X road enters the collision area, vehicle k can observe a colliding vehicle on the Y road
Figure BDA0002207808210000232
And
Figure BDA0002207808210000233
at this time, according to the vehicle k and the vehicle
Figure BDA0002207808210000234
Whether the relative distance between the two meets the condition of passing through the intersection; if the vehicle k meets the requirement, the vehicle k successfully passes through the intersection; if not, judging the reason why the vehicle k cannot successfully pass through the intersection, optimizing the running track of the vehicle by applying a corresponding cooperative control strategy and combining with dynamic planning, and calculating corresponding consumption.
The method comprises the following steps:
s1, establishing the upper layer model, including:
s1-1, determining the division stage, the state variable and the decision variable of the upper layer model, which are as follows:
the upper layer model is divided: assuming that an X road and a Y road are two one-way roads with intersections, n vehicles on the X road need to sequentially pass through m +1 intervals existing between the vehicles on the Y road, and a behavior of each single vehicle on the X road passing through the Y road is represented as a stage, and a behavior of a k-th vehicle on the X road passing through the Y road is represented as a k-th stage, where k is 1,2,3, …, n.
State variables of the upper model: the number of the vehicle intervals for the k-th vehicle on the X road to pass through the Y road in the k stage is skAnd (4) showing.
Decision variables of the upper model: the decision made in each stage indicates that the k-th vehicle on the X-th road in the k-th stage can pass through the s-th road of the Y-th roadkSelecting a particular xth of individual vehicle intervalskThe vehicles pass through at intervals.
S1-2, determining a state transition equation, a cost function and an objective function of the upper layer model, which are as follows:
the state transition equation of the upper model is as follows:
Figure BDA0002207808210000235
setting the initial condition as s0=m+1;
The state transition equation of the upper model indicates that s is equal to 1 when kk=s1The number of vehicle intervals for allowing the 1 st vehicle on the X road to pass through the Y road in the 1 st stage is m + 1; when k is 2,3, …, n, the (k-1) th vehicle on the X road selects the (X) th vehicle at the (k-1) th stagek-1Individual vehicle interval as state variable s after Y road passingkA change in (c); s0M +1 indicates that the number of vehicle intervals at which vehicles on the X road can pass through the Y road in the initial state is m + 1.
Cost function of the upper model:
Figure BDA0002207808210000236
cost function D of the upper modelk(sk,xk) Indicating a phase index function required for making a decision at the kth phase, wherein
Figure BDA0002207808210000241
S representing that the k-th vehicle on the X road passes through the Y road for the k-th vehicle to pass through under the action of the cooperative optimization control strategy in the k stagekAll possible cost costs arising from individual vehicle intervals;
Figure BDA0002207808210000242
the method is characterized in that the Y road does not directly participate in the traffic flow of the vehicles passing through the intersection, and the cost is consumed due to the fact that the speed of the vehicles is adjusted according to the following safety requirement of the vehicles of the Y road due to the fact that the vehicles in front of the Y road are influenced by the passing of the vehicles;
as described above
Figure BDA0002207808210000243
The mathematical expression of (a) is as follows:
Figure BDA0002207808210000244
Figure BDA0002207808210000245
wherein,
Figure BDA0002207808210000246
indicating that vehicle i' belongs to a vehicle between the rear vehicle at the k-1 th stage and the front vehicle at the k-th stage on the Y road.
An objective function of the upper model:
Figure BDA0002207808210000247
setting the initial condition as f0(s0)=0;
Objective function f of the upper modelk(sk) To representCumulative cost consumption of system vehicles from stage 1 to stage k, f0(s0) 0 indicates that the system cost is 0 in the initial state.
S2, establishing the lower layer model, including:
s2-1, determining a microscopic follow-up model, describing a follow-up state of the vehicle by using the microscopic follow-up model, and predicting an initial track of the vehicle; the following state of the vehicle includes a speed, an acceleration, and a position of the vehicle.
The microcosmic car-following model can be selected from specific microcosmic car-following models, such as a GIPPS car-following model, an IDM/EIDM car-following model and the like, and the specific model can be selected to be judged according to the simulation condition.
And S2-2, establishing a condition constraint model for judging whether the k vehicle on the k-th road in the k stage X can smoothly pass through the Y road.
Based on the microcosmic car-following model, the geometric characteristics of the road of the intersection are combined, and the acceleration constraint, the distance constraint and the safety constraint are added to establish a condition constraint model which accords with the characteristics of the vehicle passing the intersection.
S2-3, aiming at various possible situations that occur in the process that the k vehicle on the X road at the k stage passes through the Y road under the mixed traffic scene, a cooperative control strategy set is prepared.
S2-4, based on the vehicle initial track predicted in the step S2-1, the condition constraint model established in the step S2-2 sequentially judges that the k-th vehicle on the X-th road in the k stage passes through the S for the k-th vehicle to pass through the Y road on the Y roadkThe specific situation of each vehicle interval is determined according to the condition that the k-th vehicle on the X-th road in the k stage can pass through the s of the Y road for the k-th vehicle to pass through the Y roadkThe specific situation occurring in the process of each vehicle interval respectively makes the cooperative control strategy corresponding to the cooperative control strategy set planned in the step S2-3.
S2-5, determining the vehicles which can be optimally controlled in the vehicles participating in the k stage as target vehicles; and optimizing the running track of the target vehicle according to the cooperative control strategies made in the step S2-4 respectively, resolving the optimization problem into an optimal control problem of discrete time state constraint, and solving by using a dynamic planning idea to obtain the cooperative optimization control strategy about the target vehicle.
S2-6, calculating S for the k-th vehicle on the X road to pass through the road through which the k-th vehicle can pass through the Y road under the action of the cooperative optimization control strategy in the k-th stagekAll possible cost consumptions of one vehicle interval
Figure BDA0002207808210000251
The objective function of the lower model is:
Figure BDA0002207808210000252
Figure BDA0002207808210000253
setting an initial cost to
Figure BDA0002207808210000254
Setting an initial state variable to
Figure BDA0002207808210000255
Wherein, the target vehicle in the k stage is recorded as a target vehicle i;
Figure BDA0002207808210000256
indicating the time when the target vehicle i enters the control area in the k-th stage;
Figure BDA0002207808210000257
represents the time when the target vehicle i leaves the control area in the k-th stage; will be provided with
Figure BDA0002207808210000258
To
Figure BDA0002207808210000259
Is divided into N segments equally according to discrete time interval tauI.e. by
Figure BDA00022078082100002510
Defining a control decision time as
Figure BDA00022078082100002511
Objective function of the underlying model
Figure BDA00022078082100002512
Shows that after the target vehicle i is subjected to cooperative optimization control in the k stage, the
Figure BDA00022078082100002513
Cost consumption from time to time t;
Figure BDA00022078082100002514
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA00022078082100002515
reflecting the vehicle position and speed of the target vehicle i in the kth stage at the time t for the state variable of the lower layer model in the kth stage at the time t;
Figure BDA00022078082100002516
reflecting the vehicle acceleration of the target vehicle i at the moment t in the kth stage for the decision variable of the lower layer model at the moment t in the kth stage;
Figure BDA00022078082100002517
a decision variable set of the lower layer model at the kth stage at the time t; the initial state variable is
Figure BDA00022078082100002518
Indicates that the target vehicle i is in the k-th stage
Figure BDA00022078082100002519
Vehicle position and speed at the moment.
Decision index function of the k-th stage lower layer model at t moment
Figure BDA00022078082100002520
The method specifically comprises the following steps:
when k is equal to 1, the first step is carried out,
Figure BDA00022078082100002521
when k is 2,3, …, n,
Figure BDA00022078082100002522
Figure BDA00022078082100002523
Figure BDA00022078082100002524
Figure BDA0002207808210000261
wherein,
Figure BDA0002207808210000262
representing a decision index function of the lower layer model in the 1 st stage at the t moment;
Figure BDA0002207808210000263
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure BDA0002207808210000264
indicating a time at which the target vehicle of the k-th stage enters the control area;
Figure BDA0002207808210000265
indicating the time when the target vehicle of the k-1 stage leaves the control area;
Figure BDA0002207808210000266
indicating departure of the target vehicle from the control zone in phase kThe time of day.
S3, solving the double-layer optimization model:
s3-1, solving the upper model, and determining the decision made by each stage of the upper model when the accumulated cost consumption of the system vehicle is the lowest;
s3-2, and the vehicle optimized track of each single vehicle on the X road of the lower model through the Y road is obtained by reversely deducing the decision made at each stage of the upper model determined in the step S3-1.
And S4, solving and obtaining a mixed traffic flow cooperative decision of the double-layer optimization model aiming at system optimization by the step S3, and acting the decision on the system vehicle to control the operation of the system vehicle.
While the present invention has been described above by way of example with reference to the accompanying drawings, it is to be understood that the invention is not limited to the specific embodiments shown herein.

Claims (6)

1. A mixed traffic flow collaborative optimization control method based on double-layer planning is characterized in that a double-layer optimization model based on dynamic programming recursion is adopted to carry out mixed traffic flow collaborative decision control, and the method is suitable for different traffic scenes in the mixed traffic flow, wherein the different traffic scenes comprise ramp vehicle convergence on a highway, vehicle confluence at an intersection and vehicle passing through the intersection; the double-layer optimization model comprises an upper layer model and a lower layer model; the upper model is a vehicle sequencing problem solved by dynamic programming recursion; the lower layer model is a single vehicle track optimization problem in different traffic scenes solved by dynamic programming recursion;
the method comprises the following steps:
s1, establishing the upper layer model, including:
s1-1, determining the division stage, the state variable and the decision variable of the upper layer model, which are as follows:
the upper layer model is divided: assuming that an X road and a Y road are two one-way roads with intersections, n vehicles on the X road need to sequentially converge into or pass through m +1 intervals among the vehicles on the Y road, the behavior that each single vehicle on the X road converges into or passes through the Y road is represented as a stage, and the behavior that a k-th vehicle on the X road converges into or passes through the Y road is recorded as a k-th stage, wherein k is 1,2,3,..., n;
state variables of the upper model: the number of the vehicles which can be converged into the kth vehicle on the X road or pass through the Y road in the kth stage is skRepresents;
decision variables of the upper model: the decision made in each stage represents that the k vehicle on the X road in the k stage can merge into or pass through the s road of the Y roadkSelecting a particular xth of individual vehicle intervalskThe vehicles are alternately converged or passed;
s1-2, determining a state transition equation, a cost function and an objective function of the upper layer model, which are as follows:
the state transition equation of the upper model is as follows:
Figure FDA0002555080180000011
setting the initial condition as s0=m+1;
The state transition equation of the upper model indicates that s is equal to 1 when kk=s1The number of the vehicle intervals for the 1 st vehicle on the X road to enter or pass through the Y road in the 1 st stage is m + 1; when k is 2, 3.., n, the (k-1) th vehicle on the X road selects the (X) th vehicle at the (k-1) th stagek-1The interval of each vehicle is used as a state variable s after the vehicle is converged or passes through a Y roadkA change in (c); s0M +1 represents that the number of vehicle intervals for vehicles to merge into or pass through the road Y on the road X in the initial state is m + 1;
cost function of the upper model:
Figure FDA0002555080180000012
cost function D of the upper modelk(sk,xk) Indicating a phase index function required for making a decision at the kth phase, wherein
Figure FDA0002555080180000013
S representing that the k-th vehicle on the X road is merged into or passes through the road where the k-th vehicle can merge into or pass through the Y road under the action of the cooperative optimization control strategy in the k stagekAll possible cost costs arising from individual vehicle intervals;
Figure FDA0002555080180000021
the vehicle speed regulation system is characterized in that the vehicle does not directly participate in vehicle convergence or vehicle confluence or vehicle crossing intersection process on a Y road, and the cost is consumed due to the fact that the vehicle speed regulation of the vehicle is required for following safety because a front vehicle is influenced by the vehicle convergence or passing;
an objective function of the upper model:
Figure FDA0002555080180000022
setting the initial condition as f0(s0)=0;
Objective function f of the upper modelk(sk) Representing the cumulative cost consumption of the system vehicles from stage 1 to stage k, f0(s0) 0 denotes that the system cost is 0 in the initial state;
s2, establishing the lower layer model, including:
s2-1, determining a microscopic follow-up model, describing a follow-up state of the vehicle by using the microscopic follow-up model, and predicting an initial track of the vehicle; the following state of the vehicle comprises the speed, acceleration and position of the vehicle;
s2-2, establishing a condition constraint model for judging whether the kth vehicle on the X road at the kth stage can smoothly merge into or pass through the Y road;
s2-3, aiming at various possible situations that occur in the process that a k vehicle on a k-th stage X road converges or passes through a Y road under a mixed traffic scene, simulating a cooperative control strategy set;
s2-4, based on the vehicle initial track predicted at step S2-1The condition constraint model established in the step S2-2 sequentially judges S that the k-th vehicle on the X road of the k stage is converged into or passes through the Y road and can be converged into or pass through the Y roadkThe specific situation occurring in the process of each vehicle interval is determined according to the situation that the k-th vehicle on the X-th road in the k stage can be merged into or passes through the s of the Y roadkSpecific conditions occurring in the process of each vehicle interval respectively make the cooperative control strategy corresponding to the cooperative control strategy set planned in the step S2-3;
s2-5, determining the vehicles which can be optimally controlled in the vehicles participating in the k stage as target vehicles; optimizing the running track of the target vehicle according to the cooperative control strategies made in the step S2-4 respectively, resolving the optimization problem into an optimal control problem of discrete time state constraint, and solving by using a dynamic programming idea to obtain the cooperative optimization control strategy about the target vehicle;
s2-6, calculating the k-th vehicle on the X road or the S which can be converged by the k-th vehicle or pass through the Y road under the action of the cooperative optimization control strategy in the k-th stagekAll possible cost consumptions of one vehicle interval
Figure FDA0002555080180000023
S3, solving the double-layer optimization model:
s3-1, solving the upper model, and determining the decision made by each stage of the upper model when the accumulated cost consumption of the system vehicle is the lowest;
s3-2, reversely deducing the decision made at each stage of the upper layer model determined in the step S3-1 to obtain the vehicle optimized track of each single vehicle on the X road of the lower layer model converging into or passing through the Y road;
and S4, solving and obtaining a mixed traffic flow cooperative decision of the double-layer optimization model aiming at system optimization by the step S3, and acting the decision on the system vehicle to control the operation of the system vehicle.
2. The mixed traffic flow cooperative optimization control method based on double-layer planning according to claim 1, wherein the step S2-3 is specifically:
according to the k stage X road, the k vehicle is in the s available for the k vehicle to merge or pass through the Y roadkSelecting a particular xth of individual vehicle intervalskThe vehicles gather or pass at intervals, the k-th vehicle on the X road is marked as a vehicle k, and the X-th vehicle is positioned on the Y roadkThe front vehicle of a vehicle interval is marked as a vehicle
Figure FDA0002555080180000031
Rear vehicle as vehicle
Figure FDA0002555080180000032
Further defines a vehicle combination K, representing the combination of vehicles on the X road and the Y road directly participating in the K-th stage, and K ∈ { K }1,K2,K3Vehicle combination
Figure FDA0002555080180000033
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure FDA0002555080180000034
Vehicle k and vehicle
Figure FDA0002555080180000035
Vehicle combination
Figure FDA0002555080180000036
The vehicle participating in the kth stage comprises a vehicle k and a vehicle
Figure FDA0002555080180000037
At this time, the x-th road is positioned on the Y roadkThe interval of each vehicle has no front vehicle; vehicle combination
Figure FDA0002555080180000038
Indicating that the vehicle participating in the k-th stage has a vehicle
Figure FDA0002555080180000039
Vehicle k, now on the Xth road of YkNo rear vehicle participates at the interval of each vehicle;
based on the vehicle track predicted by the microcosmic following model, the vehicle is driven
Figure FDA00025550801800000310
Vehicle k and vehicle
Figure FDA00025550801800000311
The relationship between the two is that the vehicle k can smoothly merge into or pass through the vehicle
Figure FDA00025550801800000312
And a vehicle
Figure FDA00025550801800000313
In the interval between the vehicles, the vehicle k can not smoothly merge into or pass through the vehicle
Figure FDA00025550801800000314
And a vehicle
Figure FDA00025550801800000315
The spacing therebetween; the vehicle k can not smoothly merge into or pass through the vehicle
Figure FDA00025550801800000316
And a vehicle
Figure FDA00025550801800000317
The interval between them is divided into four cases: the first case is denoted as R1 and represents vehicle k and vehicle
Figure FDA00025550801800000318
The distance between the two adjacent layers is too close to meet the constraint condition that the two adjacent layers can be smoothly merged or passed; the second case is denoted as R2 and represents vehicle k and vehicle
Figure FDA00025550801800000319
The distance between the two adjacent layers is too close to meet the constraint condition that the two adjacent layers can be smoothly merged or passed; the third case is denoted as R3 and represents vehicle k and vehicle
Figure FDA00025550801800000320
And with vehicles
Figure FDA00025550801800000321
The process of meeting the basic spacing requirements but merging or passing in is uncomfortable; the fourth case is denoted as R4 and represents vehicle k and vehicle
Figure FDA00025550801800000322
And with vehicles
Figure FDA00025550801800000323
The constraint conditions for smooth import or passing are not met between the two groups;
based on different vehicle combinations, different vehicle type combinations and whether the vehicle k can smoothly merge into or pass through the vehicle
Figure FDA00025550801800000324
And a vehicle
Figure FDA00025550801800000325
The planning of the cooperative control strategy set under different conditions of the interval between the two strategies specifically comprises the following steps:
the collective expression for the vehicle conditions participating in the k-th phase is as follows:
Figure FDA00025550801800000326
Figure FDA00025550801800000327
or 1 or a NaN,
αkeither the number of bits is 0 or 1,
Figure FDA00025550801800000328
or 1 or a NaN,
rk0 or 1 or 2 or 3 or 4,
wherein,
Figure FDA0002555080180000041
αk
Figure FDA0002555080180000042
respectively representing vehicles
Figure FDA0002555080180000043
Vehicle k and vehicle
Figure FDA0002555080180000044
A value 0 indicates that the vehicle type is a non-optimally controllable vehicle, a value 1 indicates that the vehicle type is an optimally controllable vehicle, and a symbol NaN indicates that there is no vehicle; r iskIndicating vehicles
Figure FDA0002555080180000045
Vehicle k and vehicle
Figure FDA0002555080180000046
Relation between rk0 means that the vehicle k can smoothly merge into or pass through the vehicle
Figure FDA0002555080180000047
And a vehicle
Figure FDA0002555080180000048
Interval between rk1 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure FDA0002555080180000049
And a vehicle
Figure FDA00025550801800000410
Case of spacing between R1, Rk2 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure FDA00025550801800000411
And a vehicle
Figure FDA00025550801800000412
Case of spacing between R2, Rk3 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure FDA00025550801800000413
And a vehicle
Figure FDA00025550801800000414
Case of spacing between R3, Rk4 means that the vehicle k cannot smoothly merge into or pass through the vehicle
Figure FDA00025550801800000415
And a vehicle
Figure FDA00025550801800000416
The case of the interval therebetween R4;
a specific list of vehicle conditions participating in phase k is as follows:
Figure FDA00025550801800000417
Figure FDA0002555080180000051
and (3) aiming at the condition of the vehicle participating in the k stage, planning a corresponding cooperative control strategy:
defining decision variables of the k-th stage lower layer model at t moment
Figure FDA0002555080180000052
Reflects the acceleration of the vehicle at the moment t in the k-th stage, an
Figure FDA0002555080180000053
Wherein
Figure FDA0002555080180000054
A decision variable set of the lower layer model at the kth stage at the time t;
for vehicle k:
when in use
Figure FDA0002555080180000055
Figure FDA0002555080180000056
Figure FDA0002555080180000057
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure FDA0002555080180000058
Figure FDA0002555080180000059
when in use
Figure FDA00025550801800000510
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions:
Figure FDA00025550801800000511
when in use
Figure FDA00025550801800000512
Time t, decision variable u of vehicle k at time tk(t) satisfies the following conditions: v. ofk(t)+uk(t)τ≤ve
For vehicles
Figure FDA00025550801800000513
When in use
Figure FDA00025550801800000515
Figure FDA00025550801800000516
In time, the vehicle
Figure FDA00025550801800000517
Decision variables at time t
Figure FDA00025550801800000518
The requirements are as follows:
Figure FDA00025550801800000519
when in use
Figure FDA00025550801800000520
In time, the vehicle
Figure FDA00025550801800000521
Decision variables at time t
Figure FDA00025550801800000522
The requirements are as follows:
Figure FDA00025550801800000523
and is
Figure FDA00025550801800000524
For vehicles
Figure FDA00025550801800000525
When in use
Figure FDA00025550801800000526
Figure FDA00025550801800000527
In time, the vehicle
Figure FDA00025550801800000528
Decision variables at time t
Figure FDA00025550801800000529
The requirements are as follows:
Figure FDA00025550801800000530
when in use
Figure FDA00025550801800000531
In time, the vehicle
Figure FDA00025550801800000532
Decision variables at time t
Figure FDA00025550801800000533
The requirements are as follows:
Figure FDA00025550801800000534
and is
Figure FDA00025550801800000535
Wherein,
Figure FDA0002555080180000061
representing the vehicle condition participating in the k-th phase; τ is the reaction time of the vehicle driving; v. ofk(t)、
Figure FDA0002555080180000062
Respectively show a vehicle k and a vehicle
Figure FDA0002555080180000063
Vehicle with a steering wheel
Figure FDA0002555080180000064
Velocity at time t;
Figure FDA0002555080180000065
Figure FDA0002555080180000066
respectively a vehicle k and a vehicle
Figure FDA0002555080180000067
Vehicle with a steering wheel
Figure FDA0002555080180000068
The safe following acceleration of the vehicle at the time t is predicted according to the microcosmic following model;
Figure FDA0002555080180000069
Figure FDA00025550801800000610
respectively a vehicle k and a vehicle
Figure FDA00025550801800000611
Vehicle with a steering wheel
Figure FDA00025550801800000612
The safe following speed of the vehicle at the moment t + tau is predicted according to the microcosmic following model; v. ofeIs the desired speed;
if the vehicle in stage k-1
Figure FDA00025550801800000613
The target vehicle is subjected to cooperative optimization control, and the vehicle is taken as the vehicle in the k stage
Figure FDA00025550801800000614
The vehicle processes that will be regarded as non-optimizable control when engaged are expressed in terms of expressions as follows:
Figure FDA00025550801800000637
Figure FDA00025550801800000616
3. the double-layer planning-based mixed traffic flow cooperative optimization control method according to claim 1, wherein the objective function of the lower layer model is:
Figure FDA00025550801800000617
Figure FDA00025550801800000618
setting an initial cost to
Figure FDA00025550801800000619
Setting an initial state variable to
Figure FDA00025550801800000620
Wherein, the target vehicle in the k stage is recorded as a target vehicle i;
Figure FDA00025550801800000621
indicating the time when the target vehicle i enters the control area in the k-th stage;
Figure FDA00025550801800000622
represents the time when the target vehicle i leaves the control area in the k-th stage; will be provided with
Figure FDA00025550801800000623
To
Figure FDA00025550801800000624
Is divided equally into N segments in discrete time intervals τ', i.e.
Figure FDA00025550801800000625
Defining a control decision time as
Figure FDA00025550801800000626
Objective function of the underlying model
Figure FDA00025550801800000627
Shows that after the target vehicle i is subjected to cooperative optimization control in the k stage, the
Figure FDA00025550801800000628
Cost consumption from time to time t;
Figure FDA00025550801800000629
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure FDA00025550801800000630
reflecting the vehicle position and speed of the target vehicle i in the kth stage at the time t for the state variable of the lower layer model in the kth stage at the time t;
Figure FDA00025550801800000631
reflecting the vehicle acceleration of the target vehicle i at the moment t in the kth stage for the decision variable of the lower layer model at the moment t in the kth stage;
Figure FDA00025550801800000632
a decision variable set of the lower layer model at the kth stage at the time t; the initial state variable is
Figure FDA00025550801800000633
Indicates that the target vehicle i is in the k-th stage
Figure FDA00025550801800000634
Vehicle position and speed at the moment.
4. The double-layer planning-based mixed traffic flow cooperative optimization control method according to claim 3, wherein the k-th stage lower layer model is a decision index function at the time t
Figure FDA00025550801800000635
The method specifically comprises the following steps:
when k is equal to 1, the first step is carried out,
Figure FDA00025550801800000636
when k is 2,3, …, n,
Figure FDA0002555080180000071
Figure FDA0002555080180000072
Figure FDA0002555080180000073
Figure FDA0002555080180000074
wherein,
Figure FDA0002555080180000075
representing a decision index function of the lower layer model in the 1 st stage at the t moment;
Figure FDA0002555080180000076
representing a decision index function of the lower layer model at the kth stage at the time t;
Figure FDA0002555080180000077
indicating a time at which the target vehicle of the k-th stage enters the control area;
Figure FDA0002555080180000078
indicating the time when the target vehicle of the k-1 stage leaves the control area;
Figure FDA0002555080180000079
indicating the time at which the target vehicle leaves the control area at the k-th stage.
5. The double-deck planning-based mixed traffic flow cooperative optimization control method according to claim 1, wherein the vehicle on the Y road in the step S1-2, which does not directly participate in the vehicle confluence or vehicle crossing process, is marked as vehicle i', and the cost is consumed due to the vehicle speed adjustment of the own vehicle for the follow-up safety requirement because the front vehicle is influenced by the vehicle confluence or passing
Figure FDA00025550801800000710
The method specifically comprises the following steps:
Figure FDA00025550801800000711
Figure FDA00025550801800000712
wherein,
Figure FDA00025550801800000713
indicating that vehicle i' belongs to a vehicle between the rear vehicle at the k-1 th stage and the front vehicle at the k-th stage on the Y road.
6. The double-layer planning-based mixed traffic flow cooperative optimization control method according to any one of claims 1 to 5, characterized in that a microscopic traffic flow simulation environment is constructed, and simulation results before and after optimization under different traffic situations are compared.
CN201910887658.1A 2019-09-19 2019-09-19 Mixed traffic flow cooperative optimization control method based on double-layer planning Active CN110599772B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910887658.1A CN110599772B (en) 2019-09-19 2019-09-19 Mixed traffic flow cooperative optimization control method based on double-layer planning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910887658.1A CN110599772B (en) 2019-09-19 2019-09-19 Mixed traffic flow cooperative optimization control method based on double-layer planning

Publications (2)

Publication Number Publication Date
CN110599772A CN110599772A (en) 2019-12-20
CN110599772B true CN110599772B (en) 2020-08-28

Family

ID=68861227

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910887658.1A Active CN110599772B (en) 2019-09-19 2019-09-19 Mixed traffic flow cooperative optimization control method based on double-layer planning

Country Status (1)

Country Link
CN (1) CN110599772B (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111785088B (en) * 2020-06-23 2021-08-20 大连理工大学 Double-layer collaborative optimization method for merging network vehicle ramps
CN111724603B (en) * 2020-07-01 2021-05-18 中南大学 CAV state determination method, device, equipment and medium based on traffic track data
CN114201847B (en) * 2020-10-12 2024-06-21 西南交通大学 Simplifying method for multi-vehicle track control under mixed traffic flow
CN112735126B (en) * 2020-12-24 2023-02-03 成都格林希尔德交通科技有限公司 Mixed traffic flow cooperative optimization control method based on model predictive control
CN113232651B (en) * 2021-02-08 2022-07-22 杭州电子科技大学 Fuel automobile ecological self-adaptive cruise control method based on ADHDP
CN113055474B (en) * 2021-03-12 2022-11-11 成都格林希尔德交通科技有限公司 Micro road right transaction system
CN113033883A (en) * 2021-03-12 2021-06-25 西南交通大学 Optimization control and profit allocation method in mixed traffic flow environment
CN113034047B (en) * 2021-04-21 2023-06-30 河南工业职业技术学院 Flexible manufacturing workshop optimal scheduling method and system
CN113362619B (en) * 2021-06-04 2022-08-12 东南大学 Intelligent network vehicle ramp cooperative confluence control method and device under mixed traffic environment
CN113489793B (en) * 2021-07-07 2022-04-22 重庆大学 Expressway double-lane cooperative control method in mixed traffic scene
CN113781770B (en) * 2021-08-13 2023-03-31 同济大学 Simulation double-clutch type confluence area vehicle-road cooperative control system
CN114330984A (en) * 2021-11-12 2022-04-12 中智行(上海)交通科技有限公司 Multi-vehicle cooperative vehicle merging method and system based on vehicle-road cooperation
CN114495527B (en) * 2022-02-23 2023-10-27 长安大学 Internet-connected intersection vehicle road collaborative optimization method and system in mixed traffic environment
CN116009549B (en) * 2023-01-06 2023-12-05 清华大学 Vehicle queue dynamic reorganization method and device for non-networked tandem vehicles
CN116504069B (en) * 2023-06-26 2023-09-05 中国市政工程西南设计研究总院有限公司 Urban road network capacity optimization method, device and equipment and readable storage medium

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140210646A1 (en) * 2012-12-28 2014-07-31 Balu Subramanya Advanced parking and intersection management system
CN105809958A (en) * 2016-03-29 2016-07-27 中国科学院深圳先进技术研究院 Traffic control method and system based on intersection group
CN107730883B (en) * 2017-09-11 2020-02-04 北方工业大学 Intersection area vehicle scheduling method in Internet of vehicles environment
CN108647802B (en) * 2018-03-26 2021-06-25 复旦大学 Anti-congestion method based on double-layer traffic network model
CN108806252B (en) * 2018-06-19 2019-10-01 西南交通大学 A kind of Mixed Freeway Traffic Flows collaboration optimal control method

Also Published As

Publication number Publication date
CN110599772A (en) 2019-12-20

Similar Documents

Publication Publication Date Title
CN110599772B (en) Mixed traffic flow cooperative optimization control method based on double-layer planning
CN110570049B (en) Expressway mixed traffic flow convergence collaborative optimization bottom layer control method
CN110473404B (en) Cooperative optimization bottom layer control method for mixed traffic flow at urban crossroad
CN112735126B (en) Mixed traffic flow cooperative optimization control method based on model predictive control
CN108806252B (en) A kind of Mixed Freeway Traffic Flows collaboration optimal control method
CN112907946B (en) Traffic control method and system for automatically driving vehicle and other vehicles to run in mixed mode
CN114067559B (en) Confluence optimization control method for merging special lane for automatic vehicle into common lane
CN111267831A (en) Hybrid vehicle intelligent time-domain-variable model prediction energy management method
CN109410561B (en) Uniform and heterogeneous formation driving control method for vehicles on highway
CN111325967B (en) Intelligent networking automobile formation control method and device based on cooperative assignment
CN112258875A (en) Intersection intelligent control method based on vehicle-road cooperation and automatic driving
Wei et al. Game theoretic merging behavior control for autonomous vehicle at highway on-ramp
CN112965476A (en) High-speed unmanned vehicle trajectory planning system and method based on multi-window sampling
CN114013443A (en) Automatic driving vehicle lane change decision control method based on hierarchical reinforcement learning
CN114023108B (en) Mixed traffic lane changing model and lane changing simulation method
Lu et al. A game-theoretic approach on conflict resolution of autonomous vehicles at unsignalized intersections
Liu et al. Eco-Friendly On-Ramp Merging Strategy for Connected and Automated Vehicles in Heterogeneous Traffic
CN113823076B (en) Instant-stop and instant-walking road section blockage relieving method based on networked vehicle coordination control
CN115050183A (en) Method for generating simulated traffic flow
CN113838305A (en) Control method for motorcade to converge into intelligent networking dedicated channel
Dabiri et al. Optimized speed trajectories for cyclists, based on personal preferences and traffic light information-a stochastic dynamic programming approach
Liu et al. Integrated optimization of traffic signal timings and vehicle trajectories considering mandatory lane-changing at isolated intersections
CN117360537A (en) Vehicle control method, electronic device and storage medium
CN115830908A (en) Collaborative lane changing method and system of unmanned vehicle queue in mixed traffic flow
Ning et al. Research on Multi-Objective Optimization Models for Intersection Crossing of Connected Autonomous Vehicles With Traffic Signals

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