CN106910350A - A kind of method for finding signalized crossing group's critical path - Google Patents

A kind of method for finding signalized crossing group's critical path Download PDF

Info

Publication number
CN106910350A
CN106910350A CN201710172728.6A CN201710172728A CN106910350A CN 106910350 A CN106910350 A CN 106910350A CN 201710172728 A CN201710172728 A CN 201710172728A CN 106910350 A CN106910350 A CN 106910350A
Authority
CN
China
Prior art keywords
path
intersection
paths
degree
group
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710172728.6A
Other languages
Chinese (zh)
Other versions
CN106910350B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201710172728.6A priority Critical patent/CN106910350B/en
Publication of CN106910350A publication Critical patent/CN106910350A/en
Application granted granted Critical
Publication of CN106910350B publication Critical patent/CN106910350B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/08Controlling traffic signals according to detected number or speed of vehicles

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of method for finding signalized crossing group's critical path, comprise the following steps:(1) every data message in all sections of intersection group and intersection is obtained;(2) the discreteness coupling index I of each path origin and destination of intersection group is calculated according to these data1;(3) the retardancy coupling index I in each section of intersection group is calculated2;(4) makeAs the total correlation index of a paths,It is all section I in the path2Average value, with the path of population algorithm search intersection group, draw the paths of I highests one as critical path.The present invention is finely adjusted to path related degree model, search degree of association highest path, all paths need not be traveled through, greatly reduce operand, so as to reduce the operation time that whole signalized crossing group time-space distribution coordinates and optimizes system, the real-time and validity of control are coordinated in enhancing, increase its manageable intersection group scale.

Description

A kind of method for finding signalized crossing group's critical path
Technical field
The present invention relates to traffic signalization field, especially a kind of side for finding signalized crossing group's critical path Method.
Background technology
Intersection of the China on urban road network arterial highway is generally signalized crossing, and single point signals control is carried out respectively System.Carry out vehicle when single point signals are controlled frequently to stop at the intersection, therefore cause that road network operational efficiency is low, trip delay The traffic problems such as increase.To reduce down time of the vehicle on each intersection, intersection group time-space distribution optimization system is led to The real-time monitoring volume of traffic is crossed, control is coordinated using multiple intersections as an entirety, to reduce the congestion of intersection group.
Existing intersection group coordinated control system is divided into five parts, comprising the identification of target intersection, intersection group model Enclose delimitation, critical path retrieval, five steps of time-space distribution optimization and online real-time adjustment.Wherein, critical path search module In the model of use, the calculating of the path degree of association has one to go dimension step so that the calculation of relationship degree in each path is depended on All path retardancy degrees of association and the extreme value of the discrete type degree of association, travel through all paths in intersection unavoidably, cause the mould Block is calculated and taken considerable time, usually reach several minutes, the magnitude of traffic flow has occurred that change in this period, is weakened whole The real-time of coordinated control system, influences Adjustment effect.
The content of the invention
The technical problems to be solved by the invention are, there is provided a kind of side for finding signalized crossing group's critical path Method, can reduce operand, increase manageable intersection group scale.
In order to solve the above technical problems, the present invention provides a kind of method for finding signalized crossing group's critical path, Comprise the following steps:
(1) obtain each item data of intersection group, including intersection group scope X, intersection quantity N, crucial intersection The number of track-lines of the connected mode in K, intersection and section and each each entrance driveway of intersection, canalization scheme, queue length, car The volume of traffic of speed, each flow direction;
(2) data of all edge intersections are traveled through, the discreteness degree of association in all paths in intersection group is calculated I1, computing formula is as follows:
I1=nd/n0
nd--- the vehicle number (veh) that fleet passes through in the starting point green time of path
no--- the vehicle number (veh) that fleet passes through in path termination green time;
(3) all sections in traversal intersection group, calculate the retardancy degree of association I of its two ends intersection2, unidirectional N The retardancy degree of association of lane segment mComputing formula is:
L --- road section length (m)
--- function section length (m) in n-th of section track
--- the vehicle queue length (m) in section m nth bars track
--- deceleration distance (m)
--- perception-reaction distance (m);
(4) path of population interconnection degree I maximums in road network and then by genetic algorithm is drawn, there is the path in M bars section, its Population interconnection degree is as follows:
Preferably, in step (4), genetic algorithm is concretely comprised the following steps:
(41) for each intersection is numbered, i.e. X={ 0,1 ..., N } represents a paths, i.e. R={ x with numbered sequencen, R It is a paths, xnRepresent n-th of path intersection;The outside random search since crucial intersection, generates a two ends Path R in intersection group border, repeatedly generates multiple R, constitutes initial population
(42) existTwo intersecting paths of middle random selection:
Rand (set) --- an element is randomly choosed in set
Make X12=R1∩R2If, | X12| >=2, two intersection points of this two paths are taken at random:
x1=rand (X12),x2=rand (X12)
By R1、R2In x1、x2Between part exchange, generation filial generation path RF1、RF2;If | X12|=1, then by intersection point x ∈ X12Part afterwards exchanges, and generates RF1、RF2;The process is repeated several times, all filial generation paths that will be generated constitute progeny populationSimulate mating, the reproductive process of nature;
(43) calculateThe degree of association I in all paths, andRandomly select multiple path R:
Adaptive probability is taken, the new probability formula being selected of each path R is as follows:
pmin、pmaxIt is customized parameter
fmax——In all paths fRnMaximum
fmin——Middle fRnMinimum value
iR--- the degree of association I of path RR Descending ranking in all paths
An intersection x is randomly choosed in selected RiIf i ≠ 0 ∧ i ≠ | R | chooses xi-1Another phase Adjacent intersection x 'i
xi=rand (R)
x′i=rand (x | x ∈ X ∧ join (x, xi-1)})
Otherwise choose other edge intersection x 'i
If x 'iWith follow-up intersection xi+iIt is joined directly together, then by the x in RiReplace with x 'i;Otherwise search xi+iIt is most short Path r '={ x 'i,x′i+1,…,x′i+i′, replace the x in R with r 'i;The step simulates biological gene mutation process.
(44) recalculateThe degree of association I in all paths, screens to it, only leaves identical with parental generation population quantity A part of path, constituteR removed probability in path is:
(45) if the degree of association in the path of population is not up to the condition of convergence, also not up to predetermined maximum iteration, makesAnd return to step (42);Otherwise stop iteration and by current degree of association highest path RmaxAs critical path.
Beneficial effects of the present invention are:The present invention path related degree model is finely adjusted, per paths the degree of association with Other paths of intersection group decouple, and can individually calculate, and such that it is able to ask for approximate solution using heuritic approach, search for the degree of association Highest path, it is not necessary to travel through all paths, greatly reduce operand, so as to reduce during whole signalized crossing group Null resource coordinates and optimizes the operation time of system, and enhancing coordinates the real-time and validity of control, increases its manageable friendship Prong group's scale.
Brief description of the drawings
Fig. 1 is method of the present invention schematic flow sheet.
Specific embodiment
As shown in figure 1, a kind of method for finding signalized crossing group's critical path, comprises the following steps:
(1) investigation obtains each item data of intersection group, including the scope X of intersection group, intersection quantity N, crucial hands over The connected mode in prong K, intersection and section, the number of track-lines M of each entrance driveway, canalization scheme etc., and obtained by electronic equipment Take the real-time traffic stream informations such as queue length, speed, the volume of traffic of each flow direction and the speed of each entrance driveway.
(2) data of all edge intersections are traveled through, the discreteness degree of association in all paths in intersection group is calculated I1, computing formula is as follows:
I1=nd/n0
nd--- the vehicle number (veh) that fleet passes through in the starting point green time of path
no--- the vehicle number (veh) that fleet passes through in path termination green time
(3) all sections in traversal intersection group, calculate the retardancy degree of association I of its two ends intersection2.Unidirectional N The retardancy degree of association of lane segment mComputing formula is:
L --- road section length (m)
--- function section length (m) in n-th of section track
--- the vehicle queue length (m) in section m nth bars track
--- deceleration distance (m)
--- perception-reaction distance (m)
(4) path of population interconnection degree I maximums in road network and then by genetic algorithm is drawn.There is the path in M bars section, its Population interconnection degree is as follows:
Genetic algorithm is concretely comprised the following steps:
(40) condition that iteration terminates is set:
And maximum iteration T is set;
(41) all intersections are made for set X, is each intersection numbering, be i.e. X={ 0,1 ..., N } is represented with numbered sequence One paths.That is R={ xn, R is a paths, xnRepresent n-th of path intersection.Since crucial intersection outwards with Machine is searched for, and generation one is reached the path on intersection group border by crucial intersection:
R1={ xn|xn=rand (X) ∧ (n=0 ∨ join (xn,xn-1)}
R2={ xn|xn=rand (X) ∧ join (xn-1,xn)}
Rand (set) --- an element is randomly choosed in set
Multiple R are repeatedly generated, initial population is constituted
(42) existTwo intersecting paths of middle random selection:
Make X12=R1∩R2If, | X12| >=2, two intersection points of this two paths are taken at random:
x1=rand (X12),x2=rand (X12)
By R1、R2In x1、x2Between part exchange, generation filial generation path RF1、RF2.If | X12|=1, then by intersection point x ∈ X12Part afterwards exchanges, and generates RF1、RF2.The process is repeated several times, all filial generation paths that will be generated constitute progeny population
(43) calculateThe degree of association I in all paths, andRandomly select multiple path R:
Take adaptive probability, the selected new probability formula of path R is as follows:
pmin、pmaxIt is customized parameter
fmax——In all paths fRnMaximum
fmin——Middle fRnMinimum value
iR--- the degree of association I of path RR In descending ranking in all path degrees of association
An intersection x is randomly choosed in RiIf i ≠ 0 ∧ i ≠ | R | chooses xi-1Another Adjacent Intersections x′i
xi=rand (R)
x′i=rand (x | x ∈ X ∧ join (x, xi-1)})
Otherwise choose other edge intersection x 'i
If x 'iWith follow-up intersection xi+iIt is joined directly together, then by the x in RiReplace with x 'i;Otherwise search xi+iIt is most short Path r '={ x 'i,x′i+1,…,x′i+i′, replace the x in R with r 'i
(44) recalculateThe degree of association I in all paths, screens to it, only leaves identical with parental generation population quantity A part of path, constituteR removed probability in path is:
(45) if the degree of association in the path of population is not up to condition of convergence P, iterations is also not up to predetermined maximum and changes Generation number T, then makeAnd return to step 42);Otherwise stop iteration and by current degree of association highest path RmaxAs Critical path.
This method is intended to retrieve mould to the critical path that urban road signalized crossing group's time-space distribution coordinates system Block is targetedly improved, and improves solving speed.
Although the present invention is illustrated and has been described with regard to preferred embodiment, it is understood by those skilled in the art that Without departing from scope defined by the claims of the present invention, variations and modifications can be carried out to the present invention.

Claims (2)

1. it is a kind of find signalized crossing group critical path method, it is characterised in that comprise the following steps:
(1) obtain each item data of intersection group, including intersection group scope X, intersection quantity N, key intersection K, friendship It is the connected mode in prong and section and the number of track-lines of each each entrance driveway of intersection, canalization scheme, queue length, speed, each The volume of traffic of flow direction;
(2) data of all edge intersections are traveled through, the discreteness degree of association I in all paths in intersection group is calculated1, meter Calculate formula as follows:
I1=nd/n0
nd--- the vehicle number (veh) that fleet passes through in the starting point green time of path
no--- the vehicle number (veh) that fleet passes through in path termination green time;
(3) all sections in traversal intersection group, calculate the retardancy degree of association I of its two ends intersection2, unidirectional N tracks road The retardancy degree of association of section mComputing formula is:
I 2 m = m a x ( D 1 m , D 2 m , ... , D n m , ... , D N m ) L
D n m = d 1 n m + d 2 n m + d 3 n m
L --- road section length (m)
--- function section length (m) in n-th of section track
--- the vehicle queue length (m) in section m nth bars track
--- deceleration distance (m)
--- perception-reaction distance (m);
(4) path of population interconnection degree I maximums in road network and then by genetic algorithm is drawn, has the path in M bars section, its totality The degree of association is as follows:
I = I 1 + I 2 m ‾
I 2 m ‾ = Σ m = 1 M I 2 m M .
2. the method for finding signalized crossing group's critical path as claimed in claim 1, it is characterised in that step (4) In, genetic algorithm is concretely comprised the following steps:
(41) for each intersection is numbered, i.e. X={ 0,1 ..., N } represents a paths, i.e. R={ x with numbered sequencen, R is one Paths, xnRepresent n-th of path intersection;The outside random search since crucial intersection, one two ends of generation are in The path R on intersection group border, repeatedly generates multiple R, constitutes initial population
(42) existTwo intersecting paths of middle random selection:
Rand (set) --- an element is randomly choosed in set
Make X12=R1∩R2If, | X12| >=2, two intersection points of this two paths are taken at random:
x1=rand (X12),x2=rand (X12)
By R1、R2In x1、x2Between part exchange, generation filial generation path RF1、RF2;If | X12|=1, then by intersection point x ∈ X12Afterwards Part exchange, generate RF1、RF2;The process is repeated several times, all filial generation paths that will be generated constitute progeny populationSimulate mating, the reproductive process of nature;
(43) calculateThe degree of association I in all paths, andRandomly select multiple path R:
Adaptive probability is taken, the new probability formula being selected of each path R is as follows:
p ( R ) = p min , I x &GreaterEqual; I &OverBar; p max - ( p m a x - p min ) ( f R - f m i n ) f m a x - f m i n , I x < I &OverBar;
pmin、pmaxIt is customized parameter
fmax——In all paths fRnMaximum
fmin——Middle fRnMinimum value
iR--- the degree of association I of path RR Descending ranking in all paths
An intersection x is randomly choosed in selected RiIf i ≠ 0 ∧ i ≠ | R | chooses xi-1Another adjacent intersection Mouth x 'i
xi=rand (R)
x′i=rand (x | x ∈ X ∧ join (x, xi-1)})
Otherwise choose other edge intersection x 'i
If x 'iWith follow-up intersection xi+iIt is joined directly together, then by the x in RiReplace with x 'i;Otherwise search xi+iMost short path R '={ x 'i,x′i+1,…,x′i+i′, replace the x in R with r 'i;The step simulates biological gene mutation process.
(44) recalculateThe degree of association I in all paths, screens to it, only leaves and parental generation population quantity identical one Part path, is constitutedR removed probability in path is:
(45) if the degree of association in the path of population is not up to the condition of convergence, also not up to predetermined maximum iteration, makesAnd return to step (42);Otherwise stop iteration and by current degree of association highest path RmaxAs critical path.
CN201710172728.6A 2017-03-22 2017-03-22 A method of finding signalized crossing group critical path Active CN106910350B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710172728.6A CN106910350B (en) 2017-03-22 2017-03-22 A method of finding signalized crossing group critical path

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710172728.6A CN106910350B (en) 2017-03-22 2017-03-22 A method of finding signalized crossing group critical path

Publications (2)

Publication Number Publication Date
CN106910350A true CN106910350A (en) 2017-06-30
CN106910350B CN106910350B (en) 2019-06-18

Family

ID=59195140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710172728.6A Active CN106910350B (en) 2017-03-22 2017-03-22 A method of finding signalized crossing group critical path

Country Status (1)

Country Link
CN (1) CN106910350B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918731A (en) * 2019-01-31 2019-06-21 上海卫星工程研究所 Satellite task planning simulation analysis method and system based on critical path
CN110097757A (en) * 2019-05-22 2019-08-06 东南大学 A kind of intersection group critical path recognition methods based on depth-first search
CN112466116A (en) * 2019-12-27 2021-03-09 北京航空航天大学 Intersection group distributed control method based on edge calculation
CN113112084A (en) * 2020-07-31 2021-07-13 中国海洋大学 Training plane rear body research and development flow optimization method and device
CN113223293A (en) * 2021-05-06 2021-08-06 杭州海康威视数字技术股份有限公司 Road network simulation model construction method and device and electronic equipment
CN115100857A (en) * 2022-06-17 2022-09-23 广州运星科技有限公司 Dijkstra algorithm-based road network subregion key path identification method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006194603A (en) * 2005-01-11 2006-07-27 Sanyo Electric Co Ltd Route search device and method using genetic algorithm
CN103794065A (en) * 2014-01-24 2014-05-14 东南大学 Active urban road area signal timing parameter collaborative optimization method
CN103824446A (en) * 2013-12-13 2014-05-28 华南理工大学 Sub-area multi-intersection group decision-making control method
CN105809958A (en) * 2016-03-29 2016-07-27 中国科学院深圳先进技术研究院 Traffic control method and system based on intersection group

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006194603A (en) * 2005-01-11 2006-07-27 Sanyo Electric Co Ltd Route search device and method using genetic algorithm
CN103824446A (en) * 2013-12-13 2014-05-28 华南理工大学 Sub-area multi-intersection group decision-making control method
CN103794065A (en) * 2014-01-24 2014-05-14 东南大学 Active urban road area signal timing parameter collaborative optimization method
CN105809958A (en) * 2016-03-29 2016-07-27 中国科学院深圳先进技术研究院 Traffic control method and system based on intersection group

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李岩等: ""基于小波变换的关联交叉口群关键路径识别方法"", 《中国公路学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918731A (en) * 2019-01-31 2019-06-21 上海卫星工程研究所 Satellite task planning simulation analysis method and system based on critical path
CN109918731B (en) * 2019-01-31 2023-04-07 上海卫星工程研究所 Satellite task planning simulation analysis method and system based on critical path
CN110097757A (en) * 2019-05-22 2019-08-06 东南大学 A kind of intersection group critical path recognition methods based on depth-first search
CN110097757B (en) * 2019-05-22 2021-08-06 东南大学 Intersection group critical path identification method based on depth-first search
CN112466116A (en) * 2019-12-27 2021-03-09 北京航空航天大学 Intersection group distributed control method based on edge calculation
CN113112084A (en) * 2020-07-31 2021-07-13 中国海洋大学 Training plane rear body research and development flow optimization method and device
CN113112084B (en) * 2020-07-31 2023-01-31 中国海洋大学 Training plane rear body research and development flow optimization method and device
CN113223293A (en) * 2021-05-06 2021-08-06 杭州海康威视数字技术股份有限公司 Road network simulation model construction method and device and electronic equipment
CN113223293B (en) * 2021-05-06 2023-08-04 杭州海康威视数字技术股份有限公司 Road network simulation model construction method and device and electronic equipment
CN115100857A (en) * 2022-06-17 2022-09-23 广州运星科技有限公司 Dijkstra algorithm-based road network subregion key path identification method
CN115100857B (en) * 2022-06-17 2024-03-19 广州运星科技有限公司 Road network subzone key path identification method based on Dijkstra algorithm

Also Published As

Publication number Publication date
CN106910350B (en) 2019-06-18

Similar Documents

Publication Publication Date Title
CN106910350A (en) A kind of method for finding signalized crossing group&#39;s critical path
CN108197739B (en) Urban rail transit passenger flow prediction method
CN103593535B (en) Urban traffic complex self-adaptive network parallel simulation system and method based on multi-scale integration
Heung et al. Coordinated road-junction traffic control by dynamic programming
CN103996289B (en) A kind of flow-speeds match model and Travel Time Estimation Method and system
CN108010316A (en) A kind of road traffic multisource data fusion processing method based on road net model
CN104658297B (en) A kind of center type dynamic route guidance method based on Sarsa study
CN106408124B (en) Moving path hybrid prediction method oriented to data sparse environment
CN102693633A (en) Short-term traffic flow weighted combination prediction method
CN103680158B (en) Based on the control work zone method for dynamically partitioning of C-average fuzzy cluster analysis
CN105070042A (en) Modeling method of traffic prediction
CN111710162B (en) Urban road network traffic operation condition monitoring method and system
CN114202120A (en) Urban traffic travel time prediction method aiming at multi-source heterogeneous data
CN112489426A (en) Urban traffic flow space-time prediction scheme based on graph convolution neural network
CN100533475C (en) Traffic signal off-line time distribution optimizing method basedon particle group operation method
CN104112366B (en) Method for traffic signal optimization based on latent semantic model
CN107591004A (en) A kind of intelligent traffic guidance method based on bus or train route collaboration
CN114495500B (en) Traffic prediction method based on dual dynamic space-time diagram convolution
CN105096614A (en) Newly established crossing traffic flow prediction method based on generating type deep belief network
CN113516277B (en) Internet intelligent traffic path planning method based on road network dynamic pricing
CN109269516A (en) A kind of dynamic route guidance method based on multiple target Sarsa study
CN108682147A (en) A kind of highway traffic congestion dredges decision-making technique
CN112884229B (en) Large-scale public place people flow guiding path planning method based on differential evolution algorithm
CN107153896A (en) Traffic network path prediction method and system based on node pair entropy
Shen et al. ST-AFN: a spatial-temporal attention based fusion network for lane-level traffic flow prediction

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