CN102867424A - Area coordinating traffic control method - Google Patents

Abstract

The invention relates to the field of traffic control, in particular to an area coordinating traffic control method. The area in the area coordinating traffic control method refers to a five-intersection area. The control method comprises the steps of collecting data through a detector, storing data through a database and processing data through a computer, and finally carrying out dynamic coordinating control on a traffic signal area, that is to say, lightening moments of the red and green lights of traffic lights are output and controlled. The control scheme has the advantages of optimized traffic time, high traffic efficiency and good application effect.

Description

A kind of regional coordination traffic control method

Technical field

The present invention relates to the traffic control field, relate in particular to the regional coordination traffic control method.

Background technology

After entering the nineties, development along with computer technology and automatic control technology, and the development of traffic flow theory is perfect, urban traffic control begins to informationization, the intelligent direction development, namely adopt computing machine (crossing computing machine, zone principal computer and control center's computing machine) networking control, Real-Time Traffic Volume according to the crossing, traffic model and software by development are determined crossing traffic lights timing scheme, realize the timing optimization of whole traffic network, formed at present multiple municipal intelligent traffic control system, such as the SCOOT system of in practice widespread use, the SCATS system.

By the research to China's urban transportation present situation, it is many-sided that discovery causes the reason of traffic congestion, the reason that roadnet itself is arranged, the reason that the aspects such as traffic administration and control and urban land use exploitation are also arranged, and the crossing is as " throat " of traffic capacity in the road net, and traffic behavior is more complicated, subjects to be subject to the impact of traffic environment, the stream of people, wagon flow, be the multiple ground of traffic jam and accident, become " bottleneck " that affect Traffic Capacity of Urban Road.We need a kind of region control method of base unit that is take the crossing, coordinate to control the traffic lights of each crossing.

Also produced a lot of control methods in the traffic control field at present, but mostly all too simple and mechanical.Such as; by under the perfect condition loop that is open to traffic being retrained; the constraint of this loop is only processed the computer data that ideal is open to traffic under the state; there are not prerequisite and restriction; often can cause loop periphery and downstream wagon flow to be open to traffic extremely unreasonable; therefore constraint lacks overall treatment to loop, also is a large reason that causes not having efficient available area traffic control method.Regional traffic control should be dynamic and static combination, do specific calculation for the variety classes zone, the zone is very many in the actual life, can not and bring good effect by all complicated regional traffic control problems of the disposable solution of a kind of universal method.Therefore being necessary territorial classification, such as a kind of zone of five crossings as shown in Figure 1, is exactly to have very widely zone in the real road.This zone connects five crossings by many tracks and forms, each crossing is the cross junction that is made of two mutual square crossings in track, wherein cross junction is regularly arranged, existing traffic control field is to the control program piece very act of Single Intersection, but the control program to regional traffic is considerably less, especially for the zone of above-mentioned these class five crossings.

In addition, existing traffic control method all is to calculate green light interval by Computer basically, then transport in the traffic lights, consider again rationality, be that green time can not selected to be less than 10 seconds, can be greater than 60 seconds, in any case therefore control yet, considering aspect the present rationality factor that the green time scope is generally at 10-60 between second.Therefore traffic control method is exactly generally further to retrain in the scope of 10-60 between second in the green time scope, makes it reasonable.The span of existing green time is generally optimized at the mean value that 10-60 calculated between second, this mean value immobilizes, not with vehicle flowrate, left-hand rotation rate lamp factors vary, this control method Consideration is few, it is static control method, do not consider the factor that is open to traffic of dynamic change, namely do not have the reasonable utilization data that are open to traffic dynamically to do further to optimize, the efficient that causes being open to traffic is not high.

Summary of the invention

Technical matters to be solved by this invention provides a kind of regional coordination traffic control method, and this traffic control method is open to traffic the time by Traffic signal control and optimization, improves the efficient that is open to traffic in zone.

The present invention solves the problems of the technologies described above, by the following technical solutions:

A kind of regional coordination traffic control method, described zone connects five crossings by many tracks and forms, each crossing is the cross junction that is made of two mutual square crossings in track, wherein four cross junctions are regularly arranged, and the track connects the first control subregion of four right-angled intersection interruption-forming groined types; Another crossing is the straight line extension in a track therein, orthogonal two crossings connect respectively two crossings on two groups of parallel tracks, upstream in this crossing, and another crossing links to each other by two crossings on the two groups of parallel tracks in track and its upstream and forms the second control subregion;

With A, C, B and D representative, the intersection center of cross junction represents with O four bearing of trends of cross junction respectively; Two tracks of mutual square crossing all are back and forth two way zones in the described cross junction, wherein a track is to be kept straight on through O by A to keep straight on to the back and forth two way zone of A through O to B or by B, and another perpendicular track is to be kept straight on through O by C to keep straight on to the back and forth two way zone of C through O to D or by D;

Effective controlling party that there are 8 Vehicle Driving Cycles each crossing to, respectively:

With numeral 1 represent B drive towards controlling party that O place left-hand rotation drives towards D to,

With numeral 2 represent A drives towards B through the craspedodrome of O place controlling party to,

With numeral 3 represent D drive towards controlling party that O place left-hand rotation drives towards A to,

With numeral 4 represent C drives towards D through the craspedodrome of O place controlling party to,

With numeral 5 represent A drive towards controlling party that O place left-hand rotation drives towards C to,

With numeral 6 represent B drives towards A through the craspedodrome of O place controlling party to,

With numeral 7 represent C drive towards controlling party that O place craspedodrome drives towards B to,

With numeral 8 represent D drives towards C through the craspedodrome of O place controlling party to,

Above-mentioned 8 controlling parties are to distinguishing corresponding 8 phase places at control field, the 1st phase place, the 2nd phase place, the 3rd phase place, the 4th phase place, the 5th phase place, the 6th phase place, the 7th phase place and the 8th phase place are corresponding with above-mentioned 1,2,3,4,5,6,7 and 8 respectively;

The required hardware of this coordination traffic control method comprises a plurality of detecting devices, many traffic signal controlling machines, regional coordinations control database server, regional coordination control computing machine, traffic lights; This control step of coordinating traffic control method is as follows:

(1) a plurality of detecting devices is installed respectively in above-mentioned five crossings, be installed in the detecting device of crossing to the collection of day part traffic data, the traffic data that collects is sent to traffic signal controlling machine, and traffic signal controlling machine connects the traffic lights of each crossing;

(2) the traffic signal controlling machine regional coordination that traffic data uploaded to control center by Ethernet or the GPRS network server of controlling database;

(3) traffic data that the regional coordination control computing machine that is arranged in control center extracts database server is processed and is predicted;

According to the traffic flow data that gathers each time period, the transport need of calculating each phase place, i.e. flow rate, the calculation procedure of flow rate is as follows:

1. detect first the time headway at stop line place by the detecting device of each crossing;

The traffic flow data that 2. will gather each time period is processed, and calculates time headway, adopts h to represent average headway;

3. adopt v to represent transport need, i.e. flow rate calculates flow rate by following formula:

$v=3600\frac{1}{h},$

Above-mentioned time headway in 1. refers to when start the time interval between the adjacent two car headstocks; V wherein _iImplication be i ∈ 1,2 ..., the transport need of 8} phase place, the i.e. flow rate of i phase place;

4. detect and calculate the saturation volume rate of each crossing phase place, saturation volume rate represents with s;

The data that step ()～(three) gather and calculate comprise: the time gap that the route time data between each phase data, the Adjacent Intersections, positive negative direction green wave band center begin to each phase place green light, flow rate and the saturation volume rate data of each phase place;

(4) flow rate and the saturation volume rate data with each phase place of calculating in (three) input to regional coordination control computing machine, regional coordination control computing machine is processed each dynamic flow rate and saturation volume rate data constantly again, export the maximum green time of each crossing, maximum green time g ^MaxExpression, the maximum green time g of each crossing ^MaxThe maximum green time that specifically comprises the 1st phase place and the 5th phase place, the 2nd phase place and the 6th phase place, the 3rd phase place and the 7th phase place and the 4th phase place and the 8th phase place;

(5) above-mentioned maximum green time is retrained, constraint condition is:

If 10≤g ^Max≤ 60, value g then ^MaxIf g ^MaxLess than 10 seconds, then got 10 seconds; If g ^MaxGreater than 60 seconds, then got 60 seconds;

(6) definite area wherein a crossing be with reference to the crossing, the data that detect and calculate gained in the step (three) are inputed to regional coordination control computing machine carry out the data processing, draw in the subregion each crossing and with reference to the phase differential between the crossing, according to the phase differential that calculates, the cycle zero-time of phase place is coordinated in the crossing in the definite area;

Cycle start time information and each maximum green time information of each crossing are delivered in the traffic signal controlling machine at crossing, traffic signal controlling machine is according to cycle zero-time and maximum green time, carry out dynamic traffic signals regional coordination control, namely export and control the traffic lights bright light of each traffic lights constantly.

Further, each maximum green time draws by following method processing in the described step (four):

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\Σ}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

Wherein:

X _dBe the crucial v/c ratio of expectation, X _dValue is 0.9;

L is the total losses time in the one-period, equals to multiply by lost time the phase place number, and be 4 seconds lost time, and L equals to multiply by 4 phase places lost time, equals 16 seconds;

C is Cycle Length;

Saturated flow rate represents with s in the step (four), wherein s _iRefer to i ∈ in any one crossing 1,2 ..., 8} phase place saturation volume rate;

Use g ^MaxRepresent maximum green time, wherein:

The maximum green time that represents the 1st phase place and the 5th phase place;

The maximum green time that represents the 2nd phase place and the 6th phase place;

The maximum green time that represents the 3rd phase place and the 7th phase place;

The maximum green time that represents the 4th phase place and the 8th phase place.

As preferably, each maximum green time is processed by following method and is drawn in the described step (four):

(1) determines phase sequence according to the left-hand rotation rate of each direction: according to left-hand rotation rate and threshold value k ₂Relation, in two phase place sequence, three phase sequence a, three phase sequence b and four phase sequences, select and determine a phase sequence;

(2) calculate maximum green time: according to flow rate and the saturation volume rate of each phase place, and (1) is determined to select and definite phase sequence, just adopt following corresponding formula to calculate the maximum green time of each phase place, wherein the two phase place sequence adopts following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+2+5+6}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence a adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_1,v_5)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_2,v_6)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_4,v_8)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence b adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+2+5+6}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Four phase sequences adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

Wherein, X _dBe the crucial v/c ratio of expectation, value is 0.9,

v _iBe i ∈ 1,2 ..., the transport need of 8} phase place;

s _iBe i ∈ 1,2 ..., the saturation volume rate of 8} phase place;

The total losses time in the L one-period equals to multiply by the phase place number lost time, gets lost time 4 seconds, and the phase place number is 4, here L=4*4=16 second;

The C Cycle Length;

The maximum green time of phase place 1+2+5+6;

The maximum green time of phase place 3+4+7+8;

The maximum green time of phase place 1+5;

The maximum green time of phase place 2+6;

The maximum green time of phase place 3+7;

The maximum green time of phase place 4+8.

As preferably, the track that described many continuous first places join consists of main line, and described zone connects five crossings by 4 main lines and forms, and five crossings are respectively the 1st crossing, the 2nd crossing, the 3rd crossing, the 4th crossing and the 5th crossing; Article 4, main line is respectively: the 1st main line that connects successively the 1st crossing, the 2nd crossing, the 3rd crossing; The 2nd main line that connects successively the 4th crossing, the 5th crossing, the 3rd crossing; The 3rd main line that connects the 4th crossing, the 1st crossing; The 4th main line that connects the 5th crossing, the 2nd crossing, arrange by positive dirction the crossing on every main line;

Wherein interconnective each track consists of in the first control subregion between the 1st crossing, the 2nd crossing, the 4th crossing and the 5th crossing, and the track that connects the 2nd crossing, the 5th crossing and the 3rd crossing consists of the second control subregion;

In the first control subregion:

The track that connects the 4th crossing with the 1st crossing is defined as highway section 14;

The track that connects the 5th crossing with the 4th crossing is defined as highway section 45;

The track that connects the 2nd crossing with the 5th crossing is defined as highway section 52;

The track that connects the 1st crossing with the 2nd crossing is defined as highway section 21;

Highway section 14, highway section 45, highway section 52 and highway section 21 form the first loop;

In the second control subregion:

The track that connects the 5th crossing with the 2nd crossing is defined as highway section 25;

The track that connects the 3rd crossing with the 5th crossing is defined as highway section 53;

The track that connects the 2nd crossing with the 3rd crossing is defined as highway section 32;

Highway section 25, highway section 53 and highway section 32 form the second loop;

Regional coordination control computing machine carries out the data processing in the described step (six), draw in the subregion each crossing and with reference to the phase differential between the crossing, according to the phase differential that calculates, the cycle zero-time of crossing coordination phase place is to be undertaken by following two loop constrained procedures in the definite area:

The loop constrained procedure of described the first loop is as follows:

(1) determine that constraint condition is as follows:

Will $\begin{array}{c}{\mathrm{\&phi;}}_{(a,2),(b.2)}+{\mathrm{\&phi;}}_{(b,2),(b,4)}+{\mathrm{\&phi;}}_{(b,4),(c,4)}+{\mathrm{\&phi;}}_{(c,4),(c,2)}\\ +{\mathrm{\&phi;}}_{(c,2),(d,2)}+{\mathrm{\&phi;}}_{(d,2),(d,4)}+{\mathrm{\&phi;}}_{(d,4),(a,4)}+{\mathrm{\&phi;}}_{(a,4),(a,2)}=I,\end{array}$ This formula is defined as the A formula, wherein:

φ _{(a, 2), (b, 2)}Be that the 2nd phase place of a crossing is to the phase differential of the 2nd phase place of b crossing;

φ _{(b, 2), (b, 4)}Be that the 2nd phase place of b crossing is to the phase differential of the 4th phase place of b crossing;

φ _{(b, 4), (c, 4)}Be that the 4th phase place of b crossing is to the phase differential of the 4th phase place of c crossing;

φ _{(c, 4), (c, 2)}Be that the 4th phase place of c crossing is to the phase differential of the 2nd phase place of c crossing;

φ _{(c, 2), (d, 2)}Be that the 2nd phase place of c crossing is to the phase differential of the 2nd phase place of d crossing;

φ _{(d, 2), (d, 4}) be that the 2nd phase place of d crossing is to the phase differential of the 4th phase place of d crossing;

φ _{(d, 4), (a, 4)}Be that the 4th phase place of d crossing is to the phase differential of the 4th phase place of a crossing;

φ _{(a, 4), (a, 2)}Be that the 4th phase place of a crossing is to the phase differential of the 2nd phase place of a crossing;

The value of above-mentioned I is integer;

Wherein a, b, c and d get respectively one of them numeral in 1,2,4 and 5, but will guarantee to exist between a and b, b and c, c and d and d and a upstream and downstream relation, namely all are adjacent crossings in the first loop highway section;

φ _{(a,2)，(b，2)}=(t _ab+0.5γ _a,2L _a-0.5γ _b,2L _b)z+w _a,2-w _b,2+0.5(γ _b,2-γ _a,2)-I ₁,

φ _{(b，4)，(c，4)}=(t _bc+0.5γ _b,4L _b-0.5γ _c,4L _c)z+w _b,4-w _c,4+0.5(γ _c,4-γ _b,4)-I ₂,

φ _{(c，2),(d，2)}=(t _cd+0.5γ _b,2L _d-0.5γ _c,2L _c)z+v _d，2-v _c，2+0.5(γ _d,2-γ _c,2)-I ₃,

φ _{(d, 4), (a, 4)}=(t _Da+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4})-I ₄, this formula group is B formula group, wherein:

Z refers to the inverse in cycle, is the unknown quantity that needs restraint and calculate in above-mentioned formula;

t _AbThat a crossing is to the route time of b crossing;

t _BcThat the b crossing is to the route time of c crossing;

t _CdThat the c crossing is to the route time of d crossing;

t _DaThat the d crossing is to the route time of a crossing;

γ _{A, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of a crossing and a crossing;

γ _{B, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of b crossing and b crossing;

γ _{B, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of b crossing and b crossing;

γ _{C, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of c crossing and c crossing;

γ _{D, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of d crossing and d crossing;

γ _{C, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of c crossing and c crossing;

γ _{A, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of a crossing and a crossing;

γ _{D, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of d crossing and d crossing;

L _a, L _b, L _cAnd L _dIn the one-period, for a crossing, b crossing, c crossing and all phase places of d crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _a, L _b, L _cAnd L _dAll equate and value 16 seconds;

w _{A, 2}Be the forward green wave band position of the 2nd phase place of a crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 2}Be the forward green wave band position of the 2nd phase place of b crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 4}Be the forward green wave band position of the 4th phase place of b crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

w _{C, 4}Be the forward green wave band position of the 4th phase place of c crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{D, 2}Refer to the reverse green wave band position of the 2nd phase place of d crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{C, 2}Refer to the reverse green wave band position of the 2nd phase place of c crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{A, 4}Refer to the reverse green wave band position of the 4th phase place of a crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{D, 4}Refer to the reverse green wave band position of the 4th phase place of d crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₁Refer to the most approaching but less than (t _Ab+ 0.5 γ _{A, 2}L _a-0.5 γ _{B, 2}L _b) z+w _{A, 2}-w _{B, 2}+ 0.5 (γ _{B, 2}-γ _{A, 2}) value integer;

I ₂Refer to the most approaching but less than (t _Bc+ 0.5 γ _{B, 4}L _b-0.5 γ _{C, 4}L _c) z+w _{B, 4}-w _{C, 4}+ 0.5 (γ _{C, 4}-γ _{B, 4}) value integer;

I ₃Refer to the most approaching but less than (t _Cd+ 0.5 γ _{D, 2}L _d-0.5 γ _{C, 2}L _c) z+v _{D, 2}-v _{C, 2}+ 0.5 (γ _{D, 2}-γ _{C, 2}) value integer;

I ₄Refer to the most approaching but less than (t _Da+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4}) value integer;

(3) determine equation

${\mathrm{\&phi;}}_{(a,4)(a,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)({\mathrm{\&gamma;}}_{a,1}-{\mathrm{\&gamma;}}_{a,1}{L}_{a}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,2}{L}_{a}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,4}+{\mathrm{\&gamma;}}_{a,4}{L}_{a}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)({\mathrm{\&gamma;}}_{b,3}-{\mathrm{\&gamma;}}_{b,3}{L}_{b}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,2}{L}_{b}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,4}+{\mathrm{\&gamma;}}_{b,4}{L}_{b}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)({\mathrm{\&gamma;}}_{c,1}-{\mathrm{\&gamma;}}_{c,1}{L}_{c}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,2}{L}_{c}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,4}+{\mathrm{\&gamma;}}_{c,4}{L}_{c}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(d,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)({\mathrm{\&gamma;}}_{d,3}-{\mathrm{\&gamma;}}_{d,3}{L}_{d}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">d</figure-callout>}2}+{\mathrm{\&gamma;}}_{d,2}{L}_{d}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{d,4}+{\mathrm{\&gamma;}}_{d,4}{L}_{d}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(a,4)(a,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(d,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

This formula group is defined as C formula group, wherein:

φ _{(a, 2), (a, 4)}Be that the 2nd phase place of a crossing is to the phase differential of the 4th phase place of a crossing;

t _LRefer to the lost time of a phase place, got here 4 seconds;

γ _{A, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of a crossing and a crossing;

γ _{B, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of b crossing and d crossing;

γ _{C, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of c crossing and c crossing;

γ _{D, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of d crossing and c crossing;

(4) with B formula group and C formula group substitution A formula, draw the loop constraint formulations:

(4t _L+t _ab-γ _b,2L _b+t _bc-γ _c,4L _c+t _cd-γ _c,2L _c+t _da-γ _d,4L _d

+λ(γ _a,1)(-γ _a,1L _a+t _L)+λ(γ _b,3)(-γ _b,3L _b+t _L)

+λ(γ _c,1)(-γ _c,1L _c+t _L)+λ(γ _d,3)(-γ _d,3L _b+t _L))z

+w _a,2-w _b,2+w _b,4-w _c,4-v _c,2+v _d,2-v _d,4+v _a,4-I

=-γ _b,2-γ _c,4-γ _c,2-γ _d,4-λ(γ _a,1)γ _a,1-λ(γ _b,3)γ _b,3-λ(γ _c,1)γ _c,1-λ(γ _d,3)γ _d，3，

(5) data input area detected in the data that the described step of claim 1 ()～(three) gather and calculate is coordinated the control computing machine, recycling loop constraint formulations, calculate in the first control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing;

The loop constrained procedure of described the second loop is as follows:

(I). will $\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}+{\mathrm{\&phi;}}_{(f,4),(f,2)}+{\mathrm{\&phi;}}_{(f,2),(g,2)}+{\mathrm{\&phi;}}_{(g,2),(g,4)}\\ +{\mathrm{\&phi;}}_{(g,4),(e,4)}+{\mathrm{\&phi;}}_{(e,4),(e,2)}=I,\end{array}$ Formula is defined as the E formula, wherein:

φ _{(e, 2), (f, 2)}Be that the 2nd phase place of e crossing is to the phase differential of the 2nd phase place of f crossing;

φ _{(f, 4), (f, 2)}Be that the 4th phase place of f crossing is to the phase differential of the 2nd phase place of f crossing;

φ _{(f, 2), (g, 2)}Be that the 2nd phase place of f crossing is to the phase differential of the 2nd phase place of g crossing;

φ _{(g, 2), (g, 4)}Be that the 2nd phase place of g crossing is to the phase differential of the 4th phase place of g crossing;

φ _{(g, 4), (e, 4)}Be that the 4th phase place of g crossing is to the phase differential of the 4th phase place of e crossing;

φ _{(e, 4), (e, 2)}Be that the 4th phase place of e crossing is to the phase differential of the 2nd phase place of e crossing;

(II) can get according to the A formula:

$\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}=(t_{\mathrm{ef}}+0.5{\mathrm{\&gamma;}}_{e,2}{L}_e-0.5{\mathrm{\&gamma;}}_{f,4}{L}_f)z+w_{e,2}-w_{f,4}+0.5({\mathrm{\&gamma;}}_{f,4}-{\mathrm{\&gamma;}}_{e,2})-I_5,\\ {\mathrm{\&phi;}}_{(f,2),(g,2)}=(t_{\mathrm{fg}}+0.5{\mathrm{\&gamma;}}_{g,2}{L}_g-0.5{\mathrm{\&gamma;}}_{f,2}{L}_f)z+v_{g,2}-v_{f,2}+0.5({\mathrm{\&gamma;}}_{g,2}-{\mathrm{\&gamma;}}_{f,2})-I_6,\\ {\mathrm{\&phi;}}_{(g,4),(e,4)}=(t_{\mathrm{ge}}+0.5{\mathrm{\&gamma;}}_{e,4}{L}_e-0.5{\mathrm{\&gamma;}}_{g,4}{L}_g)z+v_{e,4}-v_{g,4}+0.5({\mathrm{\&gamma;}}_{g,4}-{\mathrm{\&gamma;}}_{e,4})-I_7,\end{array}$ This formula group is F formula group, wherein:

t _EfBe the e crossing to the route time of f crossing, unit is second;

t _FgBe the f crossing to the route time of g crossing, unit is second;

t _GeBe the g crossing to the route time of e crossing, unit is second;

γ _{E, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of e crossing and e crossing;

γ _{F, 4}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing;

γ _{G, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of g crossing and g crossing;

γ _{F, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing;

γ _{E, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing;

γ _{G, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of g crossing and g crossing;

L _e, L _fAnd L _dIn the one-period, for e crossing, f crossing and all phase places of g crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _e, L _fAnd L _dAll equal and value is 16 seconds;

w _{E, 2}Be the forward green wave band position of the 2nd phase place of e crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{F, 4}Be the forward green wave band position of the 4th phase place of f crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{G, 2}Be the reverse green wave band position of the 2nd phase place of g crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights begin;

v _{F, 2}Be the reverse green wave band position of the 4th phase place of f crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights begin;

v _{E, 4}Refer to the reverse green wave band position of the 4th phase place of e crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{G, 4}Refer to the reverse green wave band position of the 4th phase place of g crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₅Refer to the most approaching but less than (t _Ef+ 0.5 γ _{E, 2}L _e-0.5 γ _{F, 4}L _f) z+w _{E, 2}-w _{F, 4}+ 0.5 (γ _{F, 4}-γ _{E, 2}) value integer;

I ₆Refer to the most approaching but less than (t _Fg+ 0.5 γ _{G, 2}L _g-0.5 γ _{F, 2}L _f) z+v _{G, 2}-v _{F, 2}+ 0.5 (γ _{G, 2}-γ _{F, 2}) value integer;

I ₇Refer to the most approaching but less than (t _Ge+ 0.5 γ _{E, 4}L _e-0.5 γ _{G, 4}L _g) z+v _{E, 4}-v _{G, 4}+ 0.5 (γ _{G, 4}-γ _{E, 4}) value integer;

(III) determines equation

${\mathrm{\&phi;}}_{(e,4)(e,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)({\mathrm{\&gamma;}}_{e,1}-{\mathrm{\&gamma;}}_{e,1}{L}_{e}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,2}{L}_{e}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,4}+{\mathrm{\&gamma;}}_{e,4}{L}_{e}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

${\mathrm{\&phi;}}_{(f,4)(f,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)({\mathrm{\&gamma;}}_{f,1}-{\mathrm{\&gamma;}}_{f,1}{L}_{f}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,2}{L}_{f}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,4}+{\mathrm{\&gamma;}}_{f,4}{L}_{f}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)({\mathrm{\&gamma;}}_{g,3}-{\mathrm{\&gamma;}}_{g,3}{L}_{g}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,2}{L}_{g}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,4}+{\mathrm{\&gamma;}}_{g,4}{L}_{g}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(f,4)(f,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

${\mathrm{\&phi;}}_{(e,4)(e,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

This formula group is defined as G formula group, wherein:

γ _{E, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of e crossing and e crossing;

γ _{E, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing;

γ _{E, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of e crossing and e crossing;

γ _{F, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of f crossing and e crossing;

γ _{G, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of g crossing and g crossing;

(IV) draws the loop constraint formulations with F formula group and G formula group substitution E formula:

（3t _L+t _ef-γ _f,2L _f+t _fg-γ _g,4L _g-γ _g,2L _g

+λ(γ _e,1)(-γ _e,1L _e+t _L)+λ(γ _f,3)(-γ _f,3L _f+t _L)

+λ(γ _g,1)(-γ _g,1L _g+t _L))z

+w _e,2-w _f,2+w _f,4-w _g,4-v _g,2+v _e,4-I

=-γ _f,2-γ _g,4-γ _g,2-λ(γ _e,1)γ _e,1-λ(γ _f,3)γ _f,3-λ(γ _g,1)γ _g,1;

(V) data input area that step described in the claim 1 (1) is detected is coordinated the control computing machine, the loop constraint formulations of recycling step (IV), calculate in the second control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing in the second control area.

Adopt a kind of regional coordination traffic control method of such scheme design, it is the coordination traffic control method of making for the zone, five crossings described in the scheme, this traffic control method according to reasonable green time of the prior art as restriction, so that maximum green time is no less than 10 seconds, be not more than 60 seconds, further calculate maximum green time by dynamic detection data at 10-60 between second, the restriction green time.Simultaneously, according to dynamic detection data, the cycle zero-time of phase place is coordinated in the crossing in the calculative determination zone, dynamic calculation maximum green time information and cycle start time information out inputed to the control program that computer integrated draws optimization, finally control the traffic lights bright light of each traffic lights constantly.This traffic control method and existing traffic control method, difference be for five crossings of extensive existence but with, be as the basis take dynamic, as to consider more detection factors detection data with the phase data of crossing, process by computing information, the signal lamp bright light data that output is optimized are finally controlled constantly to the traffic lights bright light of outside each traffic lights.Therefore the advantage of this control program is: optimize be open to traffic time, the efficient that is open to traffic height, effect is good.

Description of drawings

Fig. 1 be the present invention based on zone, 5 crossing exemplary construction synoptic diagram;

Fig. 2 is the synoptic diagram of phase place definition of the systematic nomenclature of National Electrical Manufacturers Association (National Electrical Manufacturers Association, the NEMA) signal phase formulating and publish;

Fig. 3 is two phase place sequence synoptic diagram;

Fig. 4 is three phase sequence a synoptic diagram;

Fig. 5 is three phase sequence b synoptic diagram;

Fig. 6 is four phase sequence synoptic diagram;

Fig. 7 is loop constraint synoptic diagram;

Fig. 8 is from the phase differential synoptic diagram;

Fig. 9 is cycle initial time differential intention;

Figure 10 is phase differential and the poor synoptic diagram that concerns of phase place green light initial time;

Figure 11 is the process flow diagram of determining phase sequence according to the left-hand rotation rate of each direction.

Embodiment

The below is described in further detail technical scheme of the present invention in conjunction with above-mentioned accompanying drawing:

Embodiment 1:

A kind of regional coordination traffic control method, as shown in Figure 1 this traffic control method for the zone connect five crossings by many tracks and form, each crossing is the cross junction that is made of two mutual square crossings in track, and wherein cross junction is regularly arranged.Many the track formation main lines that continuous first place joins namely should be comprised of five crossings of 4 main lines connections in the zone.Five crossings are respectively the 1st crossing, the 2nd crossing, the 3rd crossing, the 4th crossing and the 5th crossing, the track connects the first control subregion of four right-angled intersection interruption-forming groined types, and namely interconnective each track consists of in first control subregion between wherein the 1st crossing, the 2nd crossing, the 4th crossing and the 5th crossing.Another crossing is the straight line extension in a track therein, orthogonal two crossings connect respectively two crossings on two groups of parallel tracks, upstream in this crossing, another crossing links to each other by two crossings on the two groups of parallel tracks in track and its upstream and forms the second control subregion, and the track that namely connects the 2nd crossing, the 5th crossing and the 3rd crossing consists of the second control subregion.

In above-mentioned the first control subregion:

Take the 1st crossing as starting point, the track that connects the 4th crossing is defined as highway section 14;

Take the 4th crossing as starting point, the track that connects the 5th crossing is defined as highway section 45;

Take the 5th crossing as starting point, the track that connects the 2nd crossing is defined as highway section 52;

Take the 2nd crossing as starting point, the track that connects the 1st crossing is defined as highway section 21;

Highway section 14, highway section 45, highway section 52 and highway section 21 form the first loop.

In above-mentioned the second control subregion:

Take the 2nd crossing as starting point, the track that connects the 5th crossing is defined as highway section 25;

Take the 5th crossing as starting point, the track that connects the 3rd crossing is defined as highway section 53;

Take the 3rd crossing as starting point, the track that connects the 2nd crossing is defined as highway section 32;

Highway section 25, highway section 53 and highway section 32 form the second loop.

Above-mentioned 4 main lines are respectively: the 1st main line that connects successively the 1st crossing, the 2nd crossing, the 3rd crossing; The 2nd main line that connects successively the 4th crossing, the 5th crossing, the 3rd crossing; The 3rd main line that connects the 4th crossing, the 1st crossing; The 4th main line that connects the 5th crossing, the 2nd crossing, arrange by positive dirction the crossing on every main line.

Four of each cross junction bearing of trends are respectively with A, C, B and D representative as shown in Figure 2, and the intersection center of cross junction represents with O.Two tracks of mutual square crossing all are back and forth two way zones in the cross junction, wherein a track is to be kept straight on through O by A to keep straight on to the back and forth two way zone of A through O to B or by B, and another perpendicular track is to be kept straight on through O by C to keep straight on to the back and forth two way zone of C through O to D or by D.

Effective controlling party that there are 8 Vehicle Driving Cycles each crossing to, respectively:

With numeral 1 represent B drive towards controlling party that O place left-hand rotation drives towards D to,

With numeral 2 represent A drives towards B through the craspedodrome of O place controlling party to,

With numeral 3 represent D drive towards controlling party that O place left-hand rotation drives towards A to,

With numeral 4 represent C drives towards D through the craspedodrome of O place controlling party to,

With numeral 5 represent A drive towards controlling party that O place left-hand rotation drives towards C to,

With numeral 6 represent B drives towards A through the craspedodrome of O place controlling party to,

With numeral 7 represent C drive towards controlling party that O place craspedodrome drives towards B to,

With numeral 8 represent D drives towards C through the craspedodrome of O place controlling party to,

Above-mentioned 8 controlling parties are to distinguish corresponding 8 phase places at control field, the 1st phase place, the 2nd phase place, the 3rd phase place, the 4th phase place, the 5th phase place, the 6th phase place, the 7th phase place and the 8th phase place are corresponding with above-mentioned 1,2,3,4,5,6,7 and 8 respectively, define above-mentioned letter and parameter, be convenient to statement and the calculating of following control.

The required hardware of this coordination traffic control method comprises a plurality of detecting devices, many traffic signal controlling machines, regional coordinations control database server, regional coordination control computing machine, traffic lights (being traffic lights).This control step of coordinating traffic control method is as follows:

(1) a plurality of detecting devices is installed respectively in above-mentioned five crossings, be installed in the detecting device of crossing to the collection of day part traffic data, the traffic data that collects is sent to traffic signal controlling machine, and traffic signal controlling machine connects the traffic lights of each crossing.Namely generally be that detecting device by being installed in upstream crossing is to the collection of day part traffic data, the traffic data that collects is sent to traffic signal controlling machine, when installing and using first, also must measure the road section length between 2 crossings, be stored in regional coordination and control database in the server.The crossing, upstream refers to the crossing of direction to the car upstream, current control crossing, and the data result of calculation of its collection is commonly used to control current crossing.In case of necessity, the meeting of the detecting device of upstream crossing and current crossing are jointly detected data and are gathered.

(2) traffic signal controlling machine is that it comprises with the embedded control chip of STM32F103ZE as main control unit: system initialization module, control processing module, Communications Processor Module, system's detection module, data acquisition module, data processing module, demonstration and change control parameter module, system's primary module.

The regional coordination that traffic signal controlling machine uploads to control center by Ethernet or GPRS network with the traffic data server of controlling database.Consider the difference in the different cities infrastructure construction, Ethernet has not been buried in some city underground, can upload the data to the regional coordination server of controlling database by the GPRS module, and the GPRS module of employing is emerging MG2639 among the ZTE.

(3) the regional coordination control computing machine that is arranged in control center extracts the control database traffic data of server of regional coordination and processes and predict.

According to the traffic flow data that gathers each time period, the transport need of calculating each phase place, i.e. flow rate, the calculation procedure of flow rate is as follows:

1. detect first the time headway at stop line place by the detecting device of each crossing;

The traffic flow data that 2. will gather each time period is processed, and calculates time headway, adopts h to represent average headway;

3. adopt v to represent transport need, i.e. flow rate calculates flow rate by following formula:

$v=3600\frac{1}{h},$

For example, the average headway of wagon flow is 2 seconds, and the flow rate of this wagon flow is 1800/hour (3600 seconds/hour * 0.5/second) so.

Above-mentioned time headway in 1. refers to when start the time interval between the adjacent two car headstocks; V wherein _iImplication be i ∈ 1,2 ..., the transport need of 8} phase place, the i.e. flow rate of i phase place.

4. detect and calculate the saturation volume rate of each crossing phase place, saturation volume rate represents with s, specifically records first the saturation headway at stop line place, then utilizes saturation headway to converse saturation volume rate.

The data that step ()～(three) gather and calculate comprise: the time gap that the route time data between each phase data, the Adjacent Intersections, positive negative direction green wave band center begin to each phase place green light, flow rate and the saturation volume rate data of each phase place.

(4) flow rate and the saturation volume rate data with each phase place of calculating in (three) input to regional coordination control computing machine, regional coordination control computing machine is processed each dynamic flow rate and saturation volume rate data constantly again, export the maximum green time of each crossing, maximum green time is processed by following method and is drawn:

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max ({\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

Wherein:

X _dBe the crucial v/c ratio of expectation, X _dValue is 0.9;

L is the total losses time in the one-period, equals to multiply by lost time the phase place number, and be 4 seconds lost time, and L equals to multiply by 4 phase places lost time, equals 16 seconds;

C is Cycle Length;

Saturated flow rate represents with s in the step (three), wherein s _iRefer to i ∈ in any one crossing 1,2 ..., 8} phase place saturation volume rate;

Use g ^MaxRepresent maximum green time, specifically:

The maximum green time that represents the 1st phase place and the 5th phase place;

The maximum green time that represents the 2nd phase place and the 6th phase place;

The maximum green time that represents the 3rd phase place and the 7th phase place;

The maximum green time that represents the 4th phase place and the 8th phase place.

Wherein above-mentioned v/c ratio also is called critical v/c ratio, reflects that whole signalized intersections to the load capacity of current transport need, can calculate with following formula

$v/c=\frac{v}{s(1-\frac{L}{C})}$

Wherein, v/c is the v/c ratio,

V is current transport need, and namely the vehicle arrival rate of whole crossing can be used the following formula approximate treatment

$v=\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4}),$

v _iBe the flow rate of i track group,

The saturation volume rate of the whole crossing of s can calculate with following formula

$s=\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (s_i,s_{i+4}),$

s _iBe the saturation volume rate of i track group,

The L total losses time,

C current period length.

(5) above-mentioned maximum green time is retrained, constraint condition is:

If 10≤g ^Max≤ 60, value g then ^MaxIf g ^MaxLess than 10 seconds, then got 10 seconds; If g ^MaxGreater than 60 seconds, then got 60 seconds.

(6) definite area wherein a crossing be with reference to the crossing, the data that detect and calculate gained in the step (three) are inputed to regional coordination control computing machine carry out the data processing, draw in the subregion each crossing and with reference to the phase differential between the crossing, according to the phase differential that calculates, the cycle zero-time of phase place is coordinated in the crossing in the definite area.The cycle zero-time of crossing coordination phase place is specifically undertaken by following two loop constrained procedures in the definite area:

The loop constrained procedure of the first loop is as follows:

(1) determine that constraint condition is as follows:

Will $\begin{array}{c}{\mathrm{\&phi;}}_{(a,2),(b.2)}+{\mathrm{\&phi;}}_{(b,2),(b,4)}+{\mathrm{\&phi;}}_{(b,4),(c,4)}+{\mathrm{\&phi;}}_{(c,4),(c,2)}\\ +{\mathrm{\&phi;}}_{(c,2),(d,2)}+{\mathrm{\&phi;}}_{(d,2),(d,4)}+{\mathrm{\&phi;}}_{(d,4),(a,4)}+{\mathrm{\&phi;}}_{(a,4),(a,2)}=I,\end{array}$ This formula is defined as the A formula, wherein:

φ _{(a, 2), (b, 2)}Be that the 2nd phase place of a crossing is to the phase differential of the 2nd phase place of b crossing;

φ _{(b, 2), (b, 4)}Be that the 2nd phase place of b crossing is to the phase differential of the 4th phase place of b crossing;

φ _{(b, 4), (c, 4)}Be that the 4th phase place of b crossing is to the phase differential of the 4th phase place of c crossing;

φ _{(c, 4), (c, 2)}Be that the 4th phase place of c crossing is to the phase differential of the 2nd phase place of c crossing;

φ _{(c, 2), (d, 2)}Be that the 2nd phase place of c crossing is to the phase differential of the 2nd phase place of d crossing;

φ _{(d, 2), (d, 4)}Be that the 2nd phase place of d crossing is to the phase differential of the 4th phase place of d crossing;

φ _{(d, 4), (a, 4)}Be that the 4th phase place of d crossing is to the phase differential of the 4th phase place of a crossing;

φ _{(a, 4), (a, 2)}Be that the 4th phase place of a crossing is to the phase differential of the 2nd phase place of a crossing;

The value of above-mentioned I is integer;

Wherein a, b, c and d get respectively one of them numeral in 1,2,4 and 5, but will guarantee to exist between a and b, b and c, c and d and d and a upstream and downstream relation, namely all are adjacent crossings in the first loop highway section.As shown in Figure 7, got by the A formula:

φ _{(a，2)，(b，2)}=(t _ab+0.5γ _a,2L _a-0.5γ _b,2L _b)z+w _a,2-w _b,2+0.5(γ _b,2-γ _a,2)-I ₁,

φ _{(b，4)，(c，4)}=(t _bc+0.5γ _b,4L _b-0.5γ _c,4L _c)z+w _b,4-w _c,4+0.5(γ _c,4-γ _b,4)-I ₂,

φ _{(c，2)，(d，2)}=(t _cd+0.5γ _d,2L _d-0.5γ _c,2L _c)z+v _d，2-v _c，2+0.5(γ _d,2-γ _c,2)-I ₃,

φ _{(d, 4), (a, 4)}=(t _Da+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4})-I ₄, this formula group is B formula group, wherein:

Z refers to the inverse in cycle, is the unknown quantity that needs restraint and calculate in above-mentioned formula;

t _AbThat a crossing is to the route time of b crossing;

t _BcThat the b crossing is to the route time of c crossing;

t _CdThat the c crossing is to the route time of d crossing;

t _DaThat the d crossing is to the route time of a crossing;

γ _{A, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of a crossing and a crossing;

γ _{B, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of b crossing and b crossing;

γ _{B, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of b crossing and b crossing;

γ _{C, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of c crossing and c crossing;

γ _{D, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of d crossing and d crossing;

γ _{C, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of c crossing and c crossing;

γ _{A, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of a crossing and a crossing;

γ _{D, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of d crossing and d crossing;

L _a, L _b, L _cAnd L _dIn the one-period, for a crossing, b crossing, c crossing and all phase places of d crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _a, L _b, L _cAnd L _dAll equal and value is 4*4=16 second;

w _{A, 2}Be the forward green wave band position of the 2nd phase place of a crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 2}Be the forward green wave band position of the 2nd phase place of b crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 4}Be the forward green wave band position of the 4th phase place of b crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

w _{C, 4}Be the forward green wave band position of the 4th phase place of c crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{D, 2}Refer to the reverse green wave band position of the 2nd phase place of d crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{C, 2}Refer to the reverse green wave band position of the 2nd phase place of c crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{A, 4}Refer to the reverse green wave band position of the 4th phase place of a crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{D, 4}Refer to the reverse green wave band position of the 4th phase place of d crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₁Refer to the most approaching but less than (t _Ab+ 0.5 γ _{A, 2}L _a-0.5 γ _{B, 2}L _b) z+w _{A, 2}-w _{B, 2}+ 0.5 (γ _{B, 2}-γ _{A, 2}) value integer;

I ₂Refer to the most approaching but less than (t _Bc+ 0.5 γ _{B, 4}L _b-0.5 γ _{C, 4}L _c) z+w _{B, 4}-w _{C, 4}+ 0.5 (γ _{C, 4}-γ _{B, 4}) value integer;

I ₃Refer to the most approaching but less than (t _Cd+ 0.5 γ _{D, 2}L _d-0.5 γ _{C, 2}L _c) z+v _{D, 2}-v _{C, 2}+ 0.5 (γ _{D, 2}-γ _{C, 2}) value integer;

I ₄Refer to the most approaching but less than (t _Da+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4}) value integer.

(3) with reference to figure 8, determine equation

${\mathrm{\&phi;}}_{(a,4)(a,2)}=1+t_{L}z+\mathrm{\&lambda;}\left(g_{a,1}\right)(g_{a,1}+t_{L}z)-\frac{1}{2}{r}_{a,2}-\frac{1}{2}{r}_{a,4}$

$=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)({\mathrm{\&gamma;}}_{a,1}-{\mathrm{\&gamma;}}_{a,1}{L}_{a}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,2}{L}_{a}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,4}+{\mathrm{\&gamma;}}_{a,4}{L}_{a}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=1+t_{L}z+\mathrm{\&lambda;}\left(g_{b,3}\right)(g_{b,3}+t_{L}z)-\frac{1}{2}{r}_{b,2}-\frac{1}{2}{r}_{b,4}$

$=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)({\mathrm{\&gamma;}}_{b,3}-{\mathrm{\&gamma;}}_{b,3}{L}_{b}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,2}{L}_{b}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,4}+{\mathrm{\&gamma;}}_{b,4}{L}_{b}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=1+t_{L}z+\mathrm{\&lambda;}\left(g_{c,1}\right)(g_{c,1}+t_{L}z)-\frac{1}{2}{r}_{c,2}-\frac{1}{2}{r}_{c,4}$

$=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)({\mathrm{\&gamma;}}_{c,1}-{\mathrm{\&gamma;}}_{c,1}{L}_{c}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,2}{L}_{c}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,4}+{\mathrm{\&gamma;}}_{c,4}{L}_{c}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(=d,4)}=1+t_{L}z+\mathrm{\&lambda;}\left(g_{d,3}\right)(g_{d,3}+t_{L}z)-\frac{1}{2}{r}_{d,2}-\frac{1}{2}{r}_{d,4}$

$=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)({\mathrm{\&gamma;}}_{d,3}-{\mathrm{\&gamma;}}_{d,3}{L}_{d}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">d</figure-callout>}2}+{\mathrm{\&gamma;}}_{d,2}{L}_{d}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{d,4}+{\mathrm{\&gamma;}}_{d,4}{L}_{d}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(a,4)(a,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(d,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

This formula group is defined as C formula group, wherein:

φ _{(a, 2), (a, 4)}Be that the 2nd phase place of a crossing is to the phase differential of the 4th phase place of a crossing;

t _LRefer to the lost time of a phase place, got here 4 seconds;

γ _{A, 1}Be the ratio of all phase place flow rate sums of the flow rate of the 1st phase place of a crossing and a crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{B, 3}Be the ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of b crossing and d crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{C, 1}Be the ratio of all phase place flow rate sums of the flow rate of the 1st phase place of c crossing and c crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{D, 3}Be the ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of d crossing and c crossing, obtained by the historical data statistics of above-mentioned detection.

(4) with B formula group and C formula group substitution A formula, draw the loop constraint formulations:

(4t _L+t _ab-γ _b,2L _b+t _bc-γ _c,4L _c+t _cd-γ _c,2L _c+t _da-γ _d,4L _d

+λ(γ _a,1)(-γ _a,1L _a+t _L)+λ(γ _b,3)(-γ _b,3L _b+t _L)

+λ(γ _c,1)(-γ _c,1L _c+t _L)+λ(γ _d,3)(-γ _d,3L _b+t _L))z

+w _a,2-w _b,2+w _b,4-w _c,4-v _c,2+v _d,2-v _d,4+v _a,4-I

=-γ _b,2-γ _c,4-γ _c,2-γ _d,4-λ(γ _a,1)γ _a,1-λ(γ _b,3)γ _b,3-λ(γ _c,1)γ _c,1-λ(γ _d,3)γ _d，3。

(5) detected data input area is coordinated the control computing machine in the step (), recycling loop constraint formulations, calculate in the first control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing.

In like manner, the loop constrained procedure of the second loop is as follows:

(I). will $\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}+{\mathrm{\&phi;}}_{(f,4),(f,2)}+{\mathrm{\&phi;}}_{(f,2),(g,2)}+{\mathrm{\&phi;}}_{(g,2),(g,4)}\\ +{\mathrm{\&phi;}}_{(g,4),(e,4)}+{\mathrm{\&phi;}}_{(e,4),(e,2)}=I,\end{array}$ Formula is defined as the E formula, wherein:

φ _{(e, 2), (f, 2)}Be that the 2nd phase place of e crossing is to the phase differential of the 2nd phase place of f crossing;

φ _{(f, 4), (f, 2)}Be that the 4th phase place of f crossing is to the phase differential of the 4th phase place of f crossing;

φ _{(f, 2), (g, 2)}Be that the 2nd phase place of f crossing is to the phase differential of the 2nd phase place of g crossing;

φ _{(g, 2), (g, 4)}Be that the 2nd phase place of g crossing is to the phase differential of the 4th phase place of g crossing;

φ _{(g, 4), (e, 4)}Be that the 4th phase place of g crossing is to the phase differential of the 4th phase place of e crossing;

φ _{(e, 4), (e, 2)}Be that the 4th phase place of e crossing is to the phase differential of the 2nd phase place of e crossing;

(II) can get according to the A formula:

$\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}=(t_{\mathrm{ef}}+0.5{\mathrm{\&gamma;}}_{e,2}{L}_e-0.5{\mathrm{\&gamma;}}_{f,4}{L}_f)z+w_{e,2}-w_{f,4}+0.5({\mathrm{\&gamma;}}_{f,4}-{\mathrm{\&gamma;}}_{e,2})-I_5,\\ {\mathrm{\&phi;}}_{(f,2),(g,2)}=(t_{\mathrm{fg}}+0.5{\mathrm{\&gamma;}}_{g,2}{L}_g-0.5{\mathrm{\&gamma;}}_{f,2}{L}_f)z+v_{g,2}-v_{f,2}+0.5({\mathrm{\&gamma;}}_{g,2}-{\mathrm{\&gamma;}}_{f,2})-I_6,\\ {\mathrm{\&phi;}}_{(g,4),(e,4)}=(t_{\mathrm{ge}}+0.5{\mathrm{\&gamma;}}_{e,4}{L}_e-0.5{\mathrm{\&gamma;}}_{g,4}{L}_g)z+v_{e,4}-v_{g,4}+0.5({\mathrm{\&gamma;}}_{g,4}-{\mathrm{\&gamma;}}_{e,4})-I_7,\end{array}$ This formula group is F formula group, wherein:

t _EfBe the e crossing to the route time of f crossing, unit is second;

t _FgBe the f crossing to the route time of g crossing, unit is second;

t _GeBe the g crossing to the route time of e crossing, unit is second;

γ _{E, 2}Be the ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of e crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{F, 4}Be the ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{G, 2}Be the ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of g crossing and g crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{F, 2}Be the ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{E, 4}Be the ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{G, 4}Be the ratio of all phase place flow rate sums of the flow rate of the 4th phase place of g crossing and g crossing, obtained by the historical data statistics of above-mentioned detection;

L _e, L _fAnd L _dIn the one-period, for e crossing, f crossing and all phase places of g crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _e, L _fAnd L _dAll equal and value is 16 seconds;

w _{E, 2}Be the forward green wave band position of the 2nd phase place of e crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{F, 4}Be the forward green wave band position of the 4th phase place of f crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{G, 2}Be the reverse green wave band position of the 2nd phase place of g crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights begin;

v _{F, 2}Be the reverse green wave band position of the 4th phase place of c crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights begin;

v _{E, 4}Refer to the reverse green wave band position of the 4th phase place of e crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{G, 4}Refer to the reverse green wave band position of the 4th phase place of g crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₅Refer to the most approaching but less than (t _Ef+ 0.5 γ _{E, 2}L _e-0.5 γ _{F, 4}L _f) z+w _{E, 2}-w _{F, 4}+ 0.5 (γ _{F, 4}-γ _{E, 2}) value integer;

I ₆Refer to the most approaching but less than (t _Fg+ 0.5 γ _{G, 2}L _g-0.5 γ _{F, 2}L _f) z+v _{G, 2}-v _{F, 2}+ 0.5 (γ _{G, 2}-γ _{F, 2}) value integer;

I ₇Refer to the most approaching but less than (t _Ge+ 0.5 γ _{E, 4}L _e-0.5 γ _{G, 4}L _g) z+v _{E, 4}-v _{G, 4}+ 0.5 (γ _{G, 4}-γ _{E, 4}) value integer;

(III) determines equation

${\mathrm{\&phi;}}_{(e,4)(e,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)({\mathrm{\&gamma;}}_{e,1}-{\mathrm{\&gamma;}}_{e,1}{L}_{e}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,2}{L}_{e}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,4}+{\mathrm{\&gamma;}}_{e,4}{L}_{e}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

${\mathrm{\&phi;}}_{(f,4)(f,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)({\mathrm{\&gamma;}}_{f,1}-{\mathrm{\&gamma;}}_{f,1}{L}_{f}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,2}{L}_{f}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,4}+{\mathrm{\&gamma;}}_{f,4}{L}_{f}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)({\mathrm{\&gamma;}}_{g,3}-{\mathrm{\&gamma;}}_{g,3}{L}_{g}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,2}{L}_{g}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,4}+{\mathrm{\&gamma;}}_{g,4}{L}_{g}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(f,4)(f,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

${\mathrm{\&phi;}}_{(e,4)(e,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

This formula group is defined as G formula group, wherein:

γ _{E, 3}Be the ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of e crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{E, 4}Be the ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{E, 1}Be the ratio of all phase place flow rate sums of the flow rate of the 1st phase place of e crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{F, 1}Be the ratio of all phase place flow rate sums of the flow rate of the 1st phase place of f crossing and e crossing, obtained by the historical data statistics of above-mentioned detection;

γ _{G, 3}Be the ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of g crossing and g crossing, obtained by the historical data statistics of above-mentioned detection;

(IV) draws the loop constraint formulations with F formula group and G formula group substitution E formula:

（3t _L+t _ef-γ _f,2L _f+t _fg-γ _g,4L _g-γ _g,2L _g

+λ(γ _e,1)(-γ _e,1L _e+t _L)+λ(γ _f,3)(-γ _f,3L _f+t _L)

+λ(γ _g,1)(-γ _g,1L _g+t _L))z

+w _e,2-w _f,2+w _f,4-w _g,4-v _g,2+v _e,4-I

=-γ _f,2-γ _g,4-γ _g,2-λ(γ _e,1)γ _e,1-λ(γ _f,3)γ _f,3-λ(γ _g,1)γ _g,1;

(V) data input area that step described in the claim 1 (1) is detected is coordinated the control computing machine, the loop constraint formulations of recycling step (IV), calculate in the second control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing in the second control area.

Particularly, as shown in Figure 1 the green wave band maximization problems of control area is described below:

$\mathrm{<figure-callout\; id="1"\; label="max"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">max</figure-callout>}\frac{1}{12}(b_{12}+b_{21}+b_{23}+b_{32}+b_{45}+b_{54}+b_{53}+b_{35}+b_{14}+b_{41}+b_{25}+b_{52})$

s.t.:

Article 1, the 1st of main line the crossing is the 1st crossing:

0.5b ₁₂-w _1，2≤0,

${\mathrm{\&gamma;}}_{\mathrm{1,2}}{n}_p^1{t}_{l}z+0.5b_{12}+{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{1,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{1,2}},$

0.5b ₂₁-v _1,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{1,2}}{n}_p^1{t}_{l}z+0.5b_{21}+{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{1,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{1,2}},$

Article 1, the 2nd of main line the crossing is the 2nd crossing:

0.5b ₁₂-w _2，2≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_{l}z+0.5b_{12}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,2}},$

0.5b ₂₁-v _2,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_{l}z+0.5b_{21}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,2}},$

0.5b ₂₃-w _2,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_{l}z+0.5b_{23}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,2}},$

0.5b ₃₂-v _2,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_{l}z+0.5b_{32}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,2}},$

Article 1, the 3rd of main line the crossing is the 3rd crossing:

0.5b ₂₃-w _3，2≤0,

${\mathrm{\&gamma;}}_{\mathrm{3,2}}{n}_p^3{t}_{l}z+0.5b_{23}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{3,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{3,2}},$

0.5b ₃₂-v _3,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{3,2}}{n}_p^3{t}_{l}z+0.5b_{23}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{3,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{3,2}},$

Article 2, the 1st of main line the crossing is the 4th crossing:

0.5b ₄₅-w _4,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{4,2}}{n}_p^4{t}_{L}z+0.5b_{45}+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{4,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{4,2}},$

0.5b ₅₄-v _4,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{4,2}}{n}_p^4{t}_{L}z+0.5b_{54}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{4,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{4,2}},$

Article 2, the 2nd of main line the crossing is the 5th crossing:

0.5b ₄₅-w _5，2≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_{L}z+0.5b_{45}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,2}},$

0.5b ₅₄-v _5,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_{L}z+0.5b_{54}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{5,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,2}},$

0.5b ₅₃-w _5,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_{L}z+0.5b_{53}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,2}},$

0.5b ₃₅-v _5,2≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_{L}z+0.5b_{35}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{5,2}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,2}},$

Article 2, the 3rd of main line the crossing is the 3rd crossing:

0.5b ₅₃-w _3，4≤0,

${\mathrm{\&gamma;}}_{\mathrm{3,4}}{n}_p^3{t}_{L}z+0.5b_{53}+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{3,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{3,4}},$

0.5b ₃₅-v _3,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{3,4}}{n}_p^3{t}_{L}z+0.5b_{35}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{3,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{3,4}},$

Article 3, the 1st of main line the crossing is the 4th crossing:

0.5b ₄₁-w _4,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{4,4}}{n}_p^4{t}_{L}z+0.5b_{41}+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{4,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{4,4}}$

0.5b ₁₄-v _4,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{4,4}}{n}_p^4{t}_{L}z+0.5b_{14}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{4,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{4,4}},$

Article 3, the 2nd of main line the crossing is the 1st crossing:

0.5b ₄₁-w _1,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{1,4}}{n}_p^1{t}_{L}z+0.5b_{41}+{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{1,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{1,4}},$

0.5b ₁₄-v _1,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{1,4}}{n}_p^1{t}_{L}z+0.5b_{14}+{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{1,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{1,4}},$

Article 4, the 1st of main line the crossing is the 5th crossing:

0.5b ₅₂-w _5，4≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,4}}{n}_p^5{t}_{L}z+0.5b_{52}+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,4}},$

0.5b ₂₅-v _5,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{5,4}}{n}_p^5{t}_{L}z+0.5b_{25}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{5,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{5,4}},$

Article 4, the 2nd of main line the crossing is the 2nd crossing:

0.5b ₅₂-w _2,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,4}}{n}_p^2{t}_{L}z+0.5b_{52}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{2,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,4}}$

0.5b ₂₅-v _2,4≤0,

${\mathrm{\&gamma;}}_{\mathrm{2,4}}{n}_p^2{t}_{L}z+0.5b_{25}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{2,4}}\&le;{\mathrm{\&gamma;}}_{\mathrm{2,4}},$

Article 1, the 1st of main line the highway section is highway section 12:

$(t_{12}+t_{21}+{\mathrm{\&gamma;}}_{\mathrm{1,2}}{n}_p^1{t}_L-{\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_L)z+{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{1,2}}+{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{1,2}}-{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{2,2}}-{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{2,2}}-I_{12}={\mathrm{\&gamma;}}_{\mathrm{1,2}}-{\mathrm{\&gamma;}}_{\mathrm{2,2}},$

Article 1, the 2nd of main line the highway section is highway section 23:

$(t_{23}+t_{32}+{\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_L-{\mathrm{\&gamma;}}_{\mathrm{3,2}}{n}_p^3{t}_L)z+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{2,2}}-w_{\mathrm{3,2}}-v_{\mathrm{3,2}}-I_{23}={\mathrm{\&gamma;}}_{\mathrm{2,2}}-{\mathrm{\&gamma;}}_{\mathrm{3,2}},$

Article 2, the 1st of main line the highway section is highway section 45:

$(t_{45}+t_{54}+{\mathrm{\&gamma;}}_{\mathrm{4,2}}{n}_p^4{t}_L-{\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_L)z+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{4,2}}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{4,2}}-w_{\mathrm{5,2}}-v_{\mathrm{5,2}}-I_{45}={\mathrm{\&gamma;}}_{\mathrm{4,2}}-{\mathrm{\&gamma;}}_{\mathrm{5,2}},$

Article 2, the 2nd of main line the highway section is highway section 53:

$(t_{53}+t_{35}+{\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_L-{\mathrm{\&gamma;}}_{\mathrm{3,4}}{n}_p^3{t}_L)z+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,2}}+v_{\mathrm{5,2}}-w_{\mathrm{3,4}}-v_{\mathrm{3,4}}-I_{53}={\mathrm{\&gamma;}}_{\mathrm{5,2}}-{\mathrm{\&gamma;}}_{\mathrm{3,4}},$

Article 3, the 1st of main line the highway section is highway section 41:

$(t_{41}+t_{14}+{\mathrm{\&gamma;}}_{\mathrm{4,4}}{n}_p^4{t}_L-{\mathrm{\&gamma;}}_{\mathrm{1,4}}{n}_p^1{t}_L)z+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{4,4}}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{4,4}}-{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{1,4}}-{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{1,4}}-I_{41}={\mathrm{\&gamma;}}_{\mathrm{4,4}}-{\mathrm{\&gamma;}}_{\mathrm{1,4}},$

Article 4, the 1st of main line the highway section is highway section 52:

$({\mathrm{<figure-callout\; id="52"\; label="t"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">t</figure-callout>}}_{52}+t_{25}+{\mathrm{\&gamma;}}_{\mathrm{5,4}}{n}_p^5{t}_L-{\mathrm{\&gamma;}}_{\mathrm{2,4}}{n}_p^2{t}_L)z+{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,4}}+v_{\mathrm{5,4}}-{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{2,4}}-{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{2,4}}-{\mathrm{<figure-callout\; id="52"\; label="I"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">I</figure-callout>}}_{52}={\mathrm{\&gamma;}}_{\mathrm{5,4}}-{\mathrm{\&gamma;}}_{\mathrm{2,4}},$

The first loop:

$(t_{14}-{\mathrm{\&gamma;}}_{\mathrm{1,4}}{n}_p^1{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{1,3}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{1,3}}{n}_p^1{t}_L)+t_L$

$+t_{45}-{\mathrm{\&gamma;}}_{\mathrm{5,2}}{n}_p^5{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{4,1}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{4,1}}{n}_p^4{t}_L)+t_L$

$+t_{52}-{\mathrm{\&gamma;}}_{\mathrm{2,4}}{n}_p^2{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{5,3}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{5,3}}{n}_p^5{t}_L)+t_L$

$+t_{21}-{\mathrm{\&gamma;}}_{\mathrm{2,2}}{n}_p^2{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{2,1}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{2,1}}{n}_p^2{t}_L)+t_L)z$

$-{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{1,4}}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{4,4}}-{\mathrm{<figure-callout\; id="4"\; label="w"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>\; </figure-callout>}}_{\mathrm{4,2}}+{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{5,2}}+w_{\mathrm{5,4}}-{\mathrm{<figure-callout\; id="2"\; label="w"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">w</figure-callout>}}_{\mathrm{2,4}}-{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{2,2}}+{\mathrm{<figure-callout\; id="1"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{1,2}}-I_{1452}$

$=-{\mathrm{\&gamma;}}_{\mathrm{1,4}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{1,3}}\right){\mathrm{\&gamma;}}_{\mathrm{1,3}}-{\mathrm{\&gamma;}}_{\mathrm{5,2}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{4,1}}\right){\mathrm{\&gamma;}}_{\mathrm{4,1}}$

$-{\mathrm{\&gamma;}}_{\mathrm{2,4}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{5,3}}\right){\mathrm{\&gamma;}}_{\mathrm{5,3}}-{\mathrm{\&gamma;}}_{\mathrm{2,2}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{2,1}}\right){\mathrm{\&gamma;}}_{\mathrm{2,1}},$

The second loop:

$(t_{25}-{\mathrm{\&gamma;}}_{\mathrm{2,4}}{n}_p^1{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{2,3}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{2,3}}{n}_p^2{t}_L)+t_L$

$+t_{53}-{\mathrm{\&gamma;}}_{\mathrm{3,4}}{n}_p^3{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{5,1}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{5,1}}{n}_p^5{t}_L)+t_L$

$+t_{32}-{\mathrm{\&gamma;}}_{\mathrm{3,2}}{n}_p^3{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{3,1}}\right)(t_L-{\mathrm{\&gamma;}}_{\mathrm{3,1}}{n}_p^3{t}_L)+t_L)z$

$-{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">\; <figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>\; </figure-callout>}}_{\mathrm{2,4}}+{\mathrm{<figure-callout\; id="4"\; label="v"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{5,4}}+w_{\mathrm{5,2}}-w_{\mathrm{3,4}}-{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{3,2}}+{\mathrm{<figure-callout\; id="2"\; label="v"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">v</figure-callout>}}_{\mathrm{2,2}}-I_{253}$

$=-{\mathrm{\&gamma;}}_{\mathrm{2,4}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{2,3}}\right){\mathrm{\&gamma;}}_{\mathrm{2,3}}-{\mathrm{\&gamma;}}_{\mathrm{3,4}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{5,1}}\right){\mathrm{\&gamma;}}_{\mathrm{5,1}}-{\mathrm{\&gamma;}}_{\mathrm{3,2}}-\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{\mathrm{3,1}}\right){\mathrm{\&gamma;}}_{\mathrm{3,1}},$

$\frac{1}{C_{\max }}\&le;z\&le;\frac{1}{C_{\min }}$

According to above-mentioned loop constraint condition, namely can calculate in the zone each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing.

Cycle start time information and each maximum green time information of above-mentioned each crossing are delivered in the traffic signal controlling machine at crossing, traffic signal controlling machine is according to cycle zero-time and maximum green time, carry out dynamic traffic signals regional coordination control, the traffic lights of namely exporting and control each traffic lights are the bright light moment of red light and the green light of traffic lights.

Below be that general formula is shifted in one group of loop constraint onto:

$\max \left(\underset{j=1}{\overset{n_{\mathrm{arte}}}{\mathrm{\&Sigma;}}}\underset{i=1}{\overset{n_{\mathrm{inte}}^j-1}{\mathrm{\&Sigma;}}}(k_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}{b}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}+k_{{\mathrm{\&alpha;}}_{i+1,j}{,\mathrm{\&alpha;}}_{\mathrm{ij}}}{b}_{{\mathrm{\&alpha;}}_{i+1,j},{\mathrm{\&alpha;}}_{\mathrm{ij}}})\right)$

s.t.:

$j=1,...,n_{\mathrm{arte}}:\left\{\begin{array}{c}\left\{\begin{array}{c}0.5b_{{\mathrm{\&alpha;}}_{1j},{\mathrm{\&alpha;}}_{2j}}-w_{{\mathrm{\&alpha;}}_{1j},p_{1j}}\&le;0\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{1j},p_{1j}}{n}_p^{{\mathrm{\&alpha;}}_{1j}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{1j},{\mathrm{\&alpha;}}_{2j}}+w_{{\mathrm{\&alpha;}}_{1j},p_{1j}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{1j},p_{1j}},\\ 0.5b_{{\mathrm{\&alpha;}}_{2j},{\mathrm{\&alpha;}}_{1j}}-v_{{\mathrm{\&alpha;}}_{1j},p_{1j}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{1j},p_{1j}}{n}_p^{{\mathrm{\&alpha;}}_{1j}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{2j},{\mathrm{\&alpha;}}_{1j}}+v_{{\mathrm{\&alpha;}}_{1j},p_{1j}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{1j},p_{1j}},\end{array}\right.\\ i=2,...,n_{\mathrm{inte}}^j-1:\left\{\begin{array}{c}0.5b_{{\mathrm{\&alpha;}}_{i-1,j},{\mathrm{\&alpha;}}_{\mathrm{ij}}}-w_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;0,\\ {\mathrm{\&gamma;}}_{a_{\mathrm{ij}},p_{\mathrm{ij}}}{n}_p^{{\mathrm{\&alpha;}}_{\mathrm{ij}}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{i-1,j}{\mathrm{\&alpha;}}_{\mathrm{ij}}}+w_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}},\\ 0.5b_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i-1,j}}-v_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}{n}_p^{{\mathrm{\&alpha;}}_{\mathrm{ij}}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i-1,j}}+v_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}},\\ 0.5b_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}-w_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}{n}_p^{{\mathrm{\&alpha;}}_{\mathrm{ij}}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}+w_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}},\\ 0.5b_{{\mathrm{\&alpha;}}_{i+1,j},{\mathrm{\&alpha;}}_{\mathrm{ij}}}-v_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}{n}_p^{{\mathrm{\&alpha;}}_{\mathrm{ij}}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{i+1,j},{\mathrm{\&alpha;}}_{\mathrm{ij}}}+v_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}},\end{array}\right.\\ \left\{\begin{array}{c}0.5b_{{\mathrm{\&alpha;}}_{(n_{\mathrm{inte}}^j-1),j},{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}}-w_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},p_{n_{\mathrm{inte}}^j,j}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}{p}_{n_{\mathrm{inte}}^j,j}}{n}_p^{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}}{t}_{L}z+0.5b_{{\mathrm{\&alpha;}}_{(n_{\mathrm{inte}}^j-1),j},{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}}+w_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}{p}_{n_{\mathrm{inte}}^j,j}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},p_{n_{\mathrm{inte}}^j,j}},\\ 0.5b_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},{\mathrm{\&alpha;}}_{(n_{\mathrm{inte}}^j-1),j}}-v_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},p_{n_{\mathrm{inte}}^j,j}}\&le;0,\\ {\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,h},p_{n_{\mathrm{inte}}^j,j}}{n}_p^{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j}}{t}_{L}z0.5b_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},{\mathrm{\&alpha;}}_{(n_{\mathrm{inte}}^j-1),j}}+v_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},p_{n_{\mathrm{inte}}^j,j}}\&le;{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{n_{\mathrm{inte}}^j,j},p_{n_{\mathrm{inte}}^j,j}},\end{array}\right.\end{array}\right.$

$j=1,...,n_{\mathrm{arte}}:\{i=1,...,n_{\mathrm{inte}}^j-1:\left\{\begin{array}{c}(t_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}+t_{{\mathrm{\&alpha;}}_{i+1,j}{\mathrm{\&alpha;}}_{\mathrm{ij}}}+{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}{n}_p^{{\mathrm{\&alpha;}}_{\mathrm{ij}}}{t}_L-{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{i+1,j}{p}_{i+1,j}}{n}_p^{{\mathrm{\&alpha;}}_{i+1,j}}{t}_L)z\\ +w_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}+v_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}-w_{{\mathrm{\&alpha;}}_{i+1,j},p_{i+1,j}}-v_{{\mathrm{\&alpha;}}_{i+1,j},p_{i+1,j}}\\ -I_{{\mathrm{\&alpha;}}_{\mathrm{ij}},{\mathrm{\&alpha;}}_{i+1,j}}={\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{\mathrm{ij}},p_{\mathrm{ij}}}-{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{i+1,j},p_{i+1,j}},\end{array}\right.$

$h=1,...,n_{\mathrm{loop}}:\left\{\begin{array}{c}\underset{k=1}{\overset{n_{\mathrm{inte}\_l}^h}{\mathrm{\&Sigma;}}}(t_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^e,},{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^e}}-{\mathrm{\&gamma;}}_{{\stackrel{\&OverBar;}{\mathrm{\&alpha;}}}_{\mathrm{kh}},{\stackrel{\&OverBar;}{p}}_{\mathrm{kh}}}{n}_p^{{\stackrel{\&OverBar;}{\mathrm{\&alpha;}}}_{\mathrm{kh}}}{t}_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s},(p_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s}-1)}\right)(-{\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s},(p_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s}-1)}{n}_p^{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s}}{t}_L+t_L)+t_L)z\\ +\underset{k=1}{\overset{n_{\mathrm{inte}\_l}^h}{\mathrm{\&Sigma;}}}{\mathrm{\&eta;}}_{\mathrm{kh}}({\mathrm{\&mu;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s},p_{j_{\mathrm{kh}}^a,j_{\mathrm{kh}}^s}}-{\mathrm{\&mu;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^e},p_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^e}})-I_{{\mathrm{\&alpha;}}_{j_{1h}^a,i_{1h}^s},{\mathrm{\&alpha;}}_{j_{2h}^a,i_{2h}^s},...,{\mathrm{\&alpha;}}_{j_{\left(n_{\mathrm{inte}\_l}^h\right),h}^a{,i}_{\left(n_{\mathrm{inte}\_l}^h\right),h}^s}}\\ =-\underset{k=1}{\overset{n_{\mathrm{inte}\_l}^h}{\mathrm{\&Sigma;}}}({\mathrm{\&gamma;}}_{{\stackrel{\&OverBar;}{\mathrm{\&alpha;}}}_{\mathrm{kh}},{\stackrel{\&OverBar;}{p}}_{\mathrm{kh}}}+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s},(p_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s}-1)}\right){\mathrm{\&gamma;}}_{{\mathrm{\&alpha;}}_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s},(p_{j_{\mathrm{kh}}^a,i_{\mathrm{kh}}^s}-1)}),\\ {\mathrm{\&eta;}}_{\mathrm{kh}}=\mathrm{sgn}(i_{\mathrm{kh}}^s-i_{\mathrm{kh}}^e),\\ \mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0.\end{array}\right.\end{array}\right.$

$\frac{1}{C_{\max }}\&le;z\&le;\frac{1}{C_{\min }}.$

Table one: the implication of parameters is described in the general formula

Table two: the symbol of the basic variable of signalized intersections place traffic flow, definition and unit

When implementing signal controlling, the life cycle initial time is poor convenient.The initial time in cycle is defined as: the initial time of the 1st phase place green time.Crossing a is to the poor ψ of cycle initial time of crossing b _AbAs shown in Figure 9.If the green time of the 1st phase place is zero, then the initial time in cycle is the initial time of the 2nd phase place.

In order to provide the poor ψ of cycle initial time _AbWith phase difference _{(a, i), (b, j)}, i=1,2,3,4; J=1,2,3, the relation between 4, we provide first phase difference _{(a, i), (b, j)}Poor with phase place green light initial time

Between relation.

As shown in figure 10, phase difference _{(a, i), (b, j)}Poor with phase place green light initial time

Between the pass be

The poor ψ of cycle initial time _AbPoor with phase place green light initial time

Between the pass be

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0.\end{array}\right.$

Wherein, i=1,2,3,4, j=1,2,3,4,

If the ψ that is calculated by above-mentioned publicity _AbGreater than 1, then only get fraction part as ψ _Ab

According to the public crossing between each sub regions, determine to belong in the whole control area crossing of other subregions and with reference to the phase differential between the crossing, the cycle zero-time of each crossing of calculative determination.

Embodiment 2:

All the other steps are all consistent with embodiment 1, only in step (four), added definite phase sequence step, determine first phase sequence, calculate maximum green time by the phase sequence formula of optimizing again, be used for the computing method of maximum green time in the alternative embodiment 1.

Specifically determine phase sequence according to the left-hand rotation rate of each direction, determine that flow process as shown in figure 11.When the left-hand rotation rate greater than threshold value k ₂The time, then adopt special left turn phase, k ₂Value is 2.

The as shown in figure 11 transport need of each phase place and saturation volume rate, and phase sequence all determines, just can utilize following formula to calculate the maximum green time of each phase place.(max (sltr (t), nltr (t))＞k ₂) * 2+ (max (eltr (t), wltr (t))〉k ₂) the formula explanation: if sltr (t), the higher value among the nltr (t) is greater than k ₂, then (max (sltr (t), nltr (t))〉k ₂) value be 1, otherwise be 0; If eltr (t), the higher value among the wltr (t) is greater than k ₂, (max (eltr (t), wltr (t))＞k then ₂) value be 1, otherwise be 0; Formula (max (sltr (t), nltr (t))〉k ₂) * 2+ (max (eltr (t), wltr (t))＞k ₂) result of calculation following several value is arranged: 0,1,2,3.If result of calculation is 0, adopt the two phase place sequence; Result of calculation is 1, adopts three phase sequence a; Result of calculation is 2, adopts three phase sequence b; Result of calculation is 3, adopts four phase sequences.

According to the definition of the phase place among Fig. 2, two phase place sequence, three phase sequence a, three phase sequence b and four phase sequences are respectively such as Fig. 3, Fig. 4, Fig. 5 and shown in Figure 6.Wherein, the two phase place sequence adopts following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+2+5+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence a adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence b adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+2+5+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Four phase sequences adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

${\mathrm{<figure-callout\; id="1"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="2"\; label="g"\; filenames="HDA00002193493100051.png,HDA00002193493100053.png"\; state="\{\{state\}\}">g</figure-callout>}}_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

${\mathrm{<figure-callout\; id="4"\; label="g"\; filenames="HDA00002193493100051.png"\; state="\{\{state\}\}">g</figure-callout>}}_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L)$

Wherein, X _dFor the crucial v/c ratio of expectation, get 0.9 here,

v _iBe i ∈ 1,2 ..., the transport need of 8} phase place;

s _iBe i ∈ 1,2 ..., the saturation volume rate of 8} phase place;

The total losses time in the L one-period equals to multiply by the phase place number lost time, gets lost time 4 seconds, here L=4*4=16 second;

The C Cycle Length;

The maximum green time of phase place 1+2+5+6;

The maximum green time of phase place 3+4+7+8;

The maximum green time of phase place 1+5;

The maximum green time of phase place 2+6;

The maximum green time of phase place 3+7;

The maximum green time of phase place 4+8.

Test example:

Coordinate the variation of control front and back road network traffic efficiency in order to verify contrast, adopt VISSIM software to carry out traffic simulation the control method of above-described embodiment 1 and embodiment 2, this software carries out respectively integral body and partial simulation to each main line in five crossing compositing area road networks and the road network, simulation result shows, Regional Road Network vehicle pass-through efficient has on average improved 21%, and each Trunk Road Network vehicle on average improves 26% by efficient.

Above-mentioned traffic control method is applied to the zone that a plurality of five crossings, Longwan District new city zone, Wenzhou City form, and Regional Road Network vehicle pass-through efficient has on average improved 17%.

Claims (4)

1. regional coordination traffic control method is characterized in that:

Described zone connects five crossings by many tracks and forms, each crossing is the cross junction that is made of two mutual square crossings in track, wherein four cross junctions are regularly arranged, and the track connects the first control subregion of four right-angled intersection interruption-forming groined types; Another crossing is the straight line extension in a track therein, orthogonal two crossings connect respectively two crossings on two groups of parallel tracks, upstream in this crossing, and another crossing links to each other by two crossings on the two groups of parallel tracks in track and its upstream and forms the second control subregion;

With A, C, B and D representative, the intersection center of cross junction represents with O four bearing of trends of cross junction respectively; Two tracks of mutual square crossing all are back and forth two way zones in the described cross junction, wherein a track is to be kept straight on through O by A to keep straight on to the back and forth two way zone of A through O to B or by B, and another perpendicular track is to be kept straight on through O by C to keep straight on to the back and forth two way zone of C through O to D or by D;

Effective controlling party that there are 8 Vehicle Driving Cycles each crossing to, respectively:

With numeral 1 represent B drive towards controlling party that O place left-hand rotation drives towards D to,

With numeral 2 represent A drives towards B through the craspedodrome of O place controlling party to,

With numeral 3 represent D drive towards controlling party that O place left-hand rotation drives towards A to,

With numeral 4 represent C drives towards D through the craspedodrome of O place controlling party to,

With numeral 5 represent A drive towards controlling party that O place left-hand rotation drives towards C to,

With numeral 6 represent B drives towards A through the craspedodrome of O place controlling party to,

With numeral 7 represent C drive towards controlling party that O place craspedodrome drives towards B to,

With numeral 8 represent D drives towards C through the craspedodrome of O place controlling party to,

Above-mentioned 8 controlling parties are to distinguishing corresponding 8 phase places at control field, the 1st phase place, the 2nd phase place, the 3rd phase place, the 4th phase place, the 5th phase place, the 6th phase place, the 7th phase place and the 8th phase place are corresponding with above-mentioned 1,2,3,4,5,6,7 and 8 respectively;

The required hardware of this coordination traffic control method comprises a plurality of detecting devices, many traffic signal controlling machines, regional coordinations control database server, regional coordination control computing machine, traffic lights; This control step of coordinating traffic control method is as follows:

(1) a plurality of detecting devices is installed respectively in above-mentioned five crossings, be installed in the detecting device of crossing to the collection of day part traffic data, the traffic data that collects is sent to traffic signal controlling machine, and traffic signal controlling machine connects the traffic lights of each crossing;

(2) the traffic signal controlling machine regional coordination that traffic data uploaded to control center by Ethernet or the GPRS network server of controlling database;

(3) traffic data that the regional coordination control computing machine that is arranged in control center extracts database server is processed and is predicted;

According to the traffic flow data that gathers each time period, the transport need of calculating each phase place, i.e. flow rate, the calculation procedure of flow rate is as follows:

1. detect first the time headway at stop line place by the detecting device of each crossing;

The traffic flow data that 2. will gather each time period is processed, and calculates time headway, adopts h to represent average headway;

3. adopt v to represent transport need, i.e. flow rate calculates flow rate by following formula:

$v=3600\frac{1}{h},$

Above-mentioned time headway in 1. refers to when start the time interval between the adjacent two car headstocks; V wherein _iImplication be i ∈ 1,2 ..., the transport need of 8} phase place, the i.e. flow rate of i phase place;

4. detect and calculate the saturation volume rate of each crossing phase place, saturation volume rate represents with s;

The data that step ()～(three) gather and calculate comprise: the time gap that the route time data between each phase data, the Adjacent Intersections, positive negative direction green wave band center begin to each phase place green light, flow rate and the saturation volume rate data of each phase place;

(4) flow rate and the saturation volume rate data with each phase place of calculating in (three) input to regional coordination control computing machine, regional coordination control computing machine is processed each dynamic flow rate and saturation volume rate data constantly again, export the maximum green time of each crossing, maximum green time g ^MaxExpression, the maximum green time g of each crossing ^MaxThe maximum green time that specifically comprises the 1st phase place and the 5th phase place, the 2nd phase place and the 6th phase place, the 3rd phase place and the 7th phase place and the 4th phase place and the 8th phase place;

(5) above-mentioned maximum green time is retrained, constraint condition is:

If 10≤g ^Max≤ 60, value g then ^MaxIf g ^MaxLess than 10 seconds, then got 10 seconds; If g ^MaxGreater than 60 seconds, then got 60 seconds;

(6) definite area wherein a crossing be with reference to the crossing, the data that detect and calculate gained in the step (three) are inputed to regional coordination control computing machine carry out the data processing, draw in the subregion each crossing and with reference to the phase differential between the crossing, according to the phase differential that calculates, the cycle zero-time of phase place is coordinated in the crossing in the definite area;

Cycle start time information and each maximum green time information of each crossing are delivered in the traffic signal controlling machine at crossing, traffic signal controlling machine is according to cycle zero-time and maximum green time, carry out dynamic traffic signals regional coordination control, namely export and control the traffic lights bright light of each traffic lights constantly.

2. a kind of regional coordination traffic control method according to claim 1 is characterized in that:

Each maximum green time is processed by following method and is drawn in the described step (four):

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

$g_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

Wherein:

X _dBe the crucial v/c ratio of expectation, X _dValue is 0.9;

L is the total losses time in the one-period, equals to multiply by lost time the phase place number, and be 4 seconds lost time, and L equals to multiply by 4 phase places lost time, equals 16 seconds;

C is Cycle Length;

Saturated flow rate represents with s in the step (four), wherein s _iRefer to i ∈ in any one crossing 1,2 ..., 8} phase place saturation volume rate;

Use g ^MaxRepresent maximum green time, wherein:

The maximum green time that represents the 1st phase place and the 5th phase place;

The maximum green time that represents the 2nd phase place and the 6th phase place;

The maximum green time that represents the 3rd phase place and the 7th phase place;

The maximum green time that represents the 4th phase place and the 8th phase place.

3. a kind of regional coordination traffic control method according to claim 1 is characterized in that:

Each maximum green time is processed by following method and is drawn in the described step (four):

(1) determines phase sequence according to the left-hand rotation rate of each direction: according to left-hand rotation rate and threshold value k ₂Relation, in two phase place sequence, three phase sequence a, three phase sequence b and four phase sequences, select and determine a phase sequence;

(2) calculate maximum green time: according to flow rate and the saturation volume rate of each phase place, and (1) is determined to select and definite phase sequence, just adopt following corresponding formula to calculate the maximum green time of each phase place, wherein the two phase place sequence adopts following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

$g_{1+2+5+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence a adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

$g_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+4+7+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{1,2,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Three phase sequence b adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

$g_{1+2+5+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

$g_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i\&Element;\left\{\mathrm{2,3,4}\right\}}{\mathrm{\&Sigma;}}\max (v_i,v_{i+4})}(C-L),$

Four phase sequences adopt following computing formula:

$C=\frac{X_{d}L}{X_d-\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (\frac{v_i}{s_i},\frac{v_{i+4}}{s_{i+4}})},$

$g_{1+5}^{\max }=\frac{\max (v_1,v_5)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{2+6}^{\max }=\frac{\max (v_2,v_6)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{3+7}^{\max }=\frac{\max (v_3,v_7)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

$g_{4+8}^{\max }=\frac{\max (v_4,v_8)}{\underset{i=1}{\overset{4}{\mathrm{\&Sigma;}}}\max (v_i,v_{i+4})}(C-L),$

Wherein, X _dBe the crucial v/c ratio of expectation, value is 0.9,

v _iBe i ∈ 1,2 ..., the transport need of 8} phase place;

s _iBe i ∈ 1,2 ..., the saturation volume rate of 8} phase place;

The total losses time in the L one-period equals to multiply by the phase place number lost time, gets lost time 4 seconds, and the phase place number is 4, here L=4*4=16 second;

The C Cycle Length;

The maximum green time of phase place 1+2+5+6;

The maximum green time of phase place 3+4+7+8;

The maximum green time of phase place 1+5;

The maximum green time of phase place 2+6;

The maximum green time of phase place 3+7;

The maximum green time of phase place 4+8.

4. a kind of regional coordination traffic control method according to claim 1 is characterized in that:

The track that described many continuous first places join consists of main line, and described zone connects five crossings by 4 main lines and forms, and five crossings are respectively the 1st crossing, the 2nd crossing, the 3rd crossing, the 4th crossing and the 5th crossing; Article 4, main line is respectively: the 1st main line that connects successively the 1st crossing, the 2nd crossing, the 3rd crossing; The 2nd main line that connects successively the 4th crossing, the 5th crossing, the 3rd crossing; The 3rd main line that connects the 4th crossing, the 1st crossing; The 4th main line that connects the 5th crossing, the 2nd crossing, arrange by positive dirction the crossing on every main line;

Wherein interconnective each track consists of in the first control subregion between the 1st crossing, the 2nd crossing, the 4th crossing and the 5th crossing, and the track that connects the 2nd crossing, the 5th crossing and the 3rd crossing consists of the second control subregion;

In the first control subregion:

The track that connects the 4th crossing with the 1st crossing is defined as highway section 14;

The track that connects the 5th crossing with the 4th crossing is defined as highway section 45;

The track that connects the 2nd crossing with the 5th crossing is defined as highway section 52;

The track that connects the 1st crossing with the 2nd crossing is defined as highway section 21;

Highway section 14, highway section 45, highway section 52 and highway section 21 form the first loop;

In the second control subregion:

The track that connects the 5th crossing with the 2nd crossing is defined as highway section 25;

The track that connects the 3rd crossing with the 5th crossing is defined as highway section 53;

The track that connects the 2nd crossing with the 3rd crossing is defined as highway section 32;

Highway section 25, highway section 53 and highway section 32 form the second loop;

Regional coordination control computing machine carries out the data processing in the described step (six), draw in the subregion each crossing and with reference to the phase differential between the crossing, according to the phase differential that calculates, the cycle zero-time of crossing coordination phase place is to be undertaken by following two loop constrained procedures in the definite area:

The loop constrained procedure of described the first loop is as follows:

(1) determine that constraint condition is as follows:

Will $\begin{array}{c}{\mathrm{\&phi;}}_{(a,2),(b.2)}+{\mathrm{\&phi;}}_{(b,2),(b,4)}+{\mathrm{\&phi;}}_{(b,4),(c,4)}+{\mathrm{\&phi;}}_{(c,4),(c,2)}\\ +{\mathrm{\&phi;}}_{(c,2),(d,2)}+{\mathrm{\&phi;}}_{(d,2),(d,4)}+{\mathrm{\&phi;}}_{(d,4),(a,4)}+{\mathrm{\&phi;}}_{(a,4),(a,2)}=I,\end{array}$ This formula is defined as the A formula, wherein:

φ _{(a, 2), (b, 2)}Be that the 2nd phase place of a crossing is to the phase differential of the 2nd phase place of b crossing;

φ _{(b, 2), (b, 4)}Be that the 2nd phase place of b crossing is to the phase differential of the 4th phase place of b crossing;

φ _{(b, 4), (c, 4)}Be that the 4th phase place of b crossing is to the phase differential of the 4th phase place of c crossing;

φ _{(c, 4), (c, 2)}Be that the 4th phase place of c crossing is to the phase differential of the 2nd phase place of c crossing;

φ _{(c, 2), (d, 2)}Be that the 2nd phase place of c crossing is to the phase differential of the 2nd phase place of d crossing;

φ _{(d, 2), (d, 4)}Be that the 2nd phase place of d crossing is to the phase differential of the 4th phase place of d crossing;

φ _{(d, 4), (a, 4)}Be that the 4th phase place of d crossing is to the phase differential of the 4th phase place of a crossing;

φ _{(a, 4), (a, 2)}Be that the 4th phase place of a crossing is to the phase differential of the 2nd phase place of a crossing;

The value of above-mentioned I is integer;

Wherein a, b, c and d get respectively one of them numeral in 1,2,4 and 5, but will guarantee to exist between a and b, b and c, c and d and d and a upstream and downstream relation, namely all are adjacent crossings in the first loop highway section;

φ _(a,2),(b,2)=(t _ab+0.5γ _a,2L _a-0.5γ _b,2L _b)z+w _a,2-w _b,2+0.5(γ _b,2-γ _a,2)-I ₁,

φ _{(b，4)，(c，4)}=(t _bc+0.5γ _b,4L _b-0.5γ _c,4L _c)z+w _b,4-w _c,4+0.5(γ _c,4-γ _,4)-I ₂,

φ _{(c，2)，(d，2)}=(t _cd+0.5γ _d,2L _d-0.5γ _c,2L _c)z+v _d，2-v _c，2+0.5(γ _d,2-γ _c,2)-I ₃,

φ _{(d, 4), (a, 4)}=(t _Da+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4})-I ₄, this formula group is B formula group, wherein:

Z refers to the inverse in cycle, is the unknown quantity that needs restraint and calculate in above-mentioned formula;

t _AbThat a crossing is to the route time of b crossing;

t _BcThat the b crossing is to the route time of c crossing;

t _CdThat the c crossing is to the route time of d crossing;

t _DaThat the d crossing is to the route time of a crossing;

γ _{A, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of a crossing and a crossing;

γ _{B, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of b crossing and b crossing;

γ _{B, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of b crossing and b crossing;

γ _{C, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of c crossing and c crossing;

γ _{D, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of d crossing and d crossing;

γ _{C, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of c crossing and c crossing;

γ _{A, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of a crossing and a crossing;

γ _{D, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of d crossing and d crossing;

L _a, L _b, L _cAnd L _dIn the one-period, for a crossing, b crossing, c crossing and all phase places of d crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _a, L _b, L _cAnd L _dAll equate and value 16 seconds;

w _{A, 2}Be the forward green wave band position of the 2nd phase place of a crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 2}Be the forward green wave band position of the 2nd phase place of b crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{B, 4}Be the forward green wave band position of the 4th phase place of b crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

w _{C, 4}Be the forward green wave band position of the 4th phase place of c crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{D, 2}Refer to the reverse green wave band position of the 2nd phase place of d crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{C, 2}Refer to the reverse green wave band position of the 2nd phase place of c crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights finish;

v _{A, 4}Refer to the reverse green wave band position of the 4th phase place of a crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{D, 4}Refer to the reverse green wave band position of the 4th phase place of d crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₁Refer to the most approaching but less than (t _Ab+ 0.5 γ _{A, 2}L _a-0.5 γ _{B, 2}L _b) z+w _{A, 2}-w _{B, 2}+ 0.5 (γ _{B, 2}-γ _{A, 2}) value integer;

I ₂Refer to the most approaching but less than (t _Bc+ 0.5 γ _{B, 4}L _b-0.5 γ _{C, 4}L _c) z+w _{B, 4}-w _{C, 4}+ 0.5 (γ _{C, 4}-γ _{B, 4}) value integer;

I ₃Refer to the most approaching but less than (t _Cd+0.5 γ _{D, 2}L _d-0.5 γ _{C, 2}L _c) z+v _{D, 2}-v _{C, 2}+ 0.5 (γ _{D, 2}-γ _{C, 2}) value integer;

I ₄Refer to the most approaching but less than (t _Dc+ 0.5 γ _{A, 4}L _a-0.5 γ _{D, 4}L _d) z+v _{A, 4}-v _{D, 4}+ 0.5 (γ _{D, 4}-γ _{A, 4}) value integer;

(3) determine equation

${\mathrm{\&phi;}}_{(a,4)(a,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)({\mathrm{\&gamma;}}_{a,1}-{\mathrm{\&gamma;}}_{a,1}{L}_{a}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,2}{L}_{a}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{a,4}+{\mathrm{\&gamma;}}_{a,4}{L}_{a}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)({\mathrm{\&gamma;}}_{b,3}-{\mathrm{\&gamma;}}_{b,3}{L}_{b}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,2}{L}_{b}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{b,4}+{\mathrm{\&gamma;}}_{b,4}{L}_{b}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)({\mathrm{\&gamma;}}_{c,1}-{\mathrm{\&gamma;}}_{c,1}{L}_{c}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,2}{L}_{c}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{c,4}+{\mathrm{\&gamma;}}_{c,4}{L}_{c}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(d,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)({\mathrm{\&gamma;}}_{d,3}-{\mathrm{\&gamma;}}_{d,3}{L}_{d}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{d2}+{\mathrm{\&gamma;}}_{d,2}{L}_{d}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{d,4}+{\mathrm{\&gamma;}}_{d,4}{L}_{d}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(a,4)(a,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right)(-{\mathrm{\&gamma;}}_{a,1}{L}_a+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{a,2}{L}_a-\frac{1}{2}{\mathrm{\&gamma;}}_{a,4}{L}_a)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{a,1}\right){\mathrm{\&gamma;}}_{a,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{a,2}+{\mathrm{\&gamma;}}_{a,4}),$

${\mathrm{\&phi;}}_{(b,2)(b,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right)(-{\mathrm{\&gamma;}}_{b,3}{L}_b+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{b,2}{L}_b-\frac{1}{2}{\mathrm{\&gamma;}}_{b,4}{L}_b)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{b,3}\right){\mathrm{\&gamma;}}_{b,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{b,2}+{\mathrm{\&gamma;}}_{b,4}),$

${\mathrm{\&phi;}}_{(c,4)(c,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right)(-{\mathrm{\&gamma;}}_{c,1}{L}_c+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{c,2}{L}_c-\frac{1}{2}{\mathrm{\&gamma;}}_{c,4}{L}_c)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{c,1}\right){\mathrm{\&gamma;}}_{c,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{c,2}+{\mathrm{\&gamma;}}_{c,4}),$

${\mathrm{\&phi;}}_{(d,2)(d,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right)(-{\mathrm{\&gamma;}}_{d,3}{L}_d+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{d,2}{L}_d-\frac{1}{2}{\mathrm{\&gamma;}}_{d,4}{L}_d)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{d,3}\right){\mathrm{\&gamma;}}_{d,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{d,2}+{\mathrm{\&gamma;}}_{d,4}),$

This formula group is defined as C formula group, wherein:

φ _{(a, 2), (a, 4)}Be that the 2nd phase place of a crossing is to the phase differential of the 4th phase place of a crossing;

t _LRefer to the lost time of a phase place, got here 4 seconds;

γ _{A, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of a crossing and a crossing;

γ _{B, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of b crossing and d crossing;

γ _{C, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of c crossing and c crossing;

γ _{D, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of d crossing and c crossing;

(4) with B formula group and C formula group substitution A formula, draw the loop constraint formulations:

(4t _L+t _ab-γ _b,2L _b+t _bc-γ _c,4L _c+t _cd-γ _c,2L _c+t _da-γ _d,4L _d

+λ(γ _a,1)(-γ _a,1L _a+t _L)+λ(γ _b,3)(-γ _b,3L _b+t _L)

+λ(γ _c,1)(-γ _c,1L _c+t _L)+λ(γ _d,3)(-γ _d,3L _b+t _L))z

+w _a,2-w _b,2+w _b,4-w _c,4-v _c,2+v _d,2-v _d,4+v _a,4-I

=-γ _b,2-γ _c,4-γ _c,2-γ _d,4-λ(γ _a,1)γ _a,1-λ(γ _b,3)γ _b,3-λ(γ _c,1)γ _c,1-λ(γ _d,3)γ _d，3，

(5) data input area detected in the data that the described step of claim 1 ()～(three) gather and calculate is coordinated the control computing machine, recycling loop constraint formulations, calculate in the first control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing;

The loop constrained procedure of described the second loop is as follows:

(I). will $\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}+{\mathrm{\&phi;}}_{(f,4),(f,2)}+{\mathrm{\&phi;}}_{(f,2),(g,2)}+{\mathrm{\&phi;}}_{(g,2),(g,4)}\\ +{\mathrm{\&phi;}}_{(g,4),(e,4)}+{\mathrm{\&phi;}}_{(e,4),(e,2)}=I,\end{array}$ Formula is defined as the E formula, wherein:

φ _{(e, 2), (f, 2)}Be that the 2nd phase place of e crossing is to the phase differential of the 2nd phase place of f crossing;

φ _{(f, 4), (f, 2)}Be that the 4th phase place of f crossing is to the phase differential of the 2nd phase place of f crossing;

φ _{(f, 2), (g, 2)}Be that the 2nd phase place of f crossing is to the phase differential of the 2nd phase place of g crossing;

φ _{(g, 2), (g, 4)}Be that the 2nd phase place of g crossing is to the phase differential of the 4th phase place of g crossing;

φ _{(g, 4), (e, 4)}Be that the 4th phase place of g crossing is to the phase differential of the 4th phase place of e crossing;

φ _{(e, 4), (e, 2)}Be that the 4th phase place of e crossing is to the phase differential of the 2nd phase place of e crossing;

(II) can get according to the A formula:

$\begin{array}{c}{\mathrm{\&phi;}}_{(e,2),(f,4)}=(t_{\mathrm{ef}}+0.5{\mathrm{\&gamma;}}_{e,2}{L}_e-0.5{\mathrm{\&gamma;}}_{f,4}{L}_f)z+w_{e,2}-w_{f,4}+0.5({\mathrm{\&gamma;}}_{f,4}-{\mathrm{\&gamma;}}_{e,2})-I_5,\\ {\mathrm{\&phi;}}_{(f,2),(g,2)}=(t_{\mathrm{fg}}+0.5{\mathrm{\&gamma;}}_{g,2}{L}_g-0.5{\mathrm{\&gamma;}}_{f,2}{L}_f)z+v_{g,2}-v_{f,2}+0.5({\mathrm{\&gamma;}}_{g,2}-{\mathrm{\&gamma;}}_{f,2})-I_6,\\ {\mathrm{\&phi;}}_{(g,4),(e,4)}=(t_{\mathrm{ge}}+0.5{\mathrm{\&gamma;}}_{e,4}{L}_e-0.5{\mathrm{\&gamma;}}_{g,4}{L}_g)z+v_{e,4}-v_{g,4}+0.5({\mathrm{\&gamma;}}_{g,4}-{\mathrm{\&gamma;}}_{e,4})-I_7;\end{array}$ This formula group is F formula group, wherein:

t _EfBe the e crossing to the route time of f crossing, unit is second;

t _FgBe the f crossing to the route time of g crossing, unit is second;

t _GeBe the g crossing to the route time of e crossing, unit is second;

γ _{E, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of e crossing and e crossing;

γ _{F, 4}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing;

γ _{G, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of g crossing and g crossing;

γ _{F, 2}The ratio of all phase place flow rate sums of the flow rate of the 2nd phase place of f crossing and f crossing;

γ _{E, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing;

γ _{G, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of g crossing and g crossing;

L _e, L _fAnd L _dIn the one-period, for e crossing, f crossing and all phase places of g crossing, because Phase-switching, and that time that causes the crossing not used by any direction wagon flow, namely sum lost time of each phase place; L _e, L _fAnd L _dAll equal and value is 16 seconds;

w _{E, 2}Be the forward green wave band position of the 2nd phase place of e crossing, be defined as the time gap that forward green wave band center to the 2 phase place green lights begin;

w _{F, 4}Be the forward green wave band position of the 4th phase place of f crossing, be defined as the time gap that forward green wave band center to the 4 phase place green lights begin;

v _{G, 2}Be the reverse green wave band position of the 2nd phase place of g crossing, be defined as the time gap that reverse green wave band center to the 2 phase place green lights begin;

v _{F, 2}Be the reverse green wave band position of the 4th phase place of f crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights begin;

v _{E, 4}Refer to the reverse green wave band position of the 4th phase place of e crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

v _{G, 4}Refer to the reverse green wave band position of the 4th phase place of g crossing, be defined as the time gap that reverse green wave band center to the 4 phase place green lights finish;

I ₅Refer to the most approaching but less than (t _Ef+ 0.5 γ _{E, 2}L _e-0.5 γ _{F, 4}L _f) z+w _{E, 2}-w _{F, 4}+ 0.5 (γ _{F, 4}-γ _{E, 2}) value integer;

I ₆Refer to the most approaching but less than (t _Fg+ 0.5 γ _{G, 2}L _g-0.5 γ _{F, 2}L _f) z+v _{G, 2}-v _{F, 2}+ 0.5 (γ _{G, 2}-γ _{F, 2}) value integer;

I ₇Refer to the most approaching but less than (t _Ge+ 0.5 γ _{E, 4}L _e-0.5 γ _{G, 4}L _g) z+v _{E, 4}-v _{G, 4}+ 0.5 (γ _{G, 4}-γ _{E, 4}) value integer;

(III) determines equation

${\mathrm{\&phi;}}_{(e,4)(e,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)({\mathrm{\&gamma;}}_{e,1}-{\mathrm{\&gamma;}}_{e,1}{L}_{e}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,2}{L}_{e}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{e,4}+{\mathrm{\&gamma;}}_{e,4}{L}_{e}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

${\mathrm{\&phi;}}_{(f,4)(f,2)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)({\mathrm{\&gamma;}}_{f,1}-{\mathrm{\&gamma;}}_{f,1}{L}_{f}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,2}{L}_{f}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{f,4}+{\mathrm{\&gamma;}}_{f,4}{L}_{f}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=1+t_{L}z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)({\mathrm{\&gamma;}}_{g,3}-{\mathrm{\&gamma;}}_{g,3}{L}_{g}z+t_{L}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,2}{L}_{g}z)-\frac{1}{2}(1-{\mathrm{\&gamma;}}_{g,4}+{\mathrm{\&gamma;}}_{g,4}{L}_{g}z)$

$=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

$\mathrm{\&lambda;}\left(x\right)=\left\{\begin{array}{cc}0& \mathrm{ifx}=0,\\ 1& \mathrm{ifx}\&NotEqual;0,\end{array}\right.$

That is:

${\mathrm{\&phi;}}_{(f,4)(f,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right)(-{\mathrm{\&gamma;}}_{f,1}{L}_f+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{f,2}{L}_f-\frac{1}{2}{\mathrm{\&gamma;}}_{f,4}{L}_f)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{f,1}\right){\mathrm{\&gamma;}}_{f,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{f,2}+{\mathrm{\&gamma;}}_{f,4}),$

${\mathrm{\&phi;}}_{(g,2)(g,4)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right)(-{\mathrm{\&gamma;}}_{g,3}{L}_g+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{g,2}{L}_g-\frac{1}{2}{\mathrm{\&gamma;}}_{g,4}{L}_g)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{g,3}\right){\mathrm{\&gamma;}}_{g,3}+\frac{1}{2}({\mathrm{\&gamma;}}_{g,2}+{\mathrm{\&gamma;}}_{g,4}),$

${\mathrm{\&phi;}}_{(e,4)(e,2)}=(t_L+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right)(-{\mathrm{\&gamma;}}_{e,1}{L}_e+t_L)-\frac{1}{2}{\mathrm{\&gamma;}}_{e,2}{L}_e-\frac{1}{2}{\mathrm{\&gamma;}}_{e,4}{L}_e)z+\mathrm{\&lambda;}\left({\mathrm{\&gamma;}}_{e,1}\right){\mathrm{\&gamma;}}_{e,1}+\frac{1}{2}({\mathrm{\&gamma;}}_{e,2}+{\mathrm{\&gamma;}}_{e,4}),$

This formula group is defined as G formula group, wherein:

γ _{E, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of e crossing and e crossing;

γ _{E, 4}The ratio of all phase place flow rate sums of the flow rate of the 4th phase place of e crossing and e crossing;

γ _{E, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of e crossing and e crossing;

γ _{F, 1}The ratio of all phase place flow rate sums of the flow rate of the 1st phase place of f crossing and e crossing;

γ _{G, 3}The ratio of all phase place flow rate sums of the flow rate of the 3rd phase place of g crossing and g crossing;

(IV) draws the loop constraint formulations with F formula group and G formula group substitution E formula:

(3t _L+t _ef-γ _f,2L _f+t _fg-γ _g,4L _g-γ _g,2L _g

+λ(γ _e,1)(-γ _e，1L _e+t _L)+λ(γ _f,3)(-γ _f,3L _f+t _L)

+λ(γ _g,1)(-γ _g,1L _g+t _L))z

+w _e,2-w _f,2+w _f,4-w _g,4-v _g,2+v _e,4-I

=-γ _f,2-γ _g,4-γ _g,2-λ(γ _e,1)γ _e,1-λ(γ _f,3)γ _f,3-λ(γ _g,1)γ _g,1;

(V) data input area that step described in the claim 1 (1) is detected is coordinated the control computing machine, the loop constraint formulations of recycling step (IV), calculate in the second control area each crossing and with reference to the phase differential between the crossing, determine the cycle zero-time of each crossing in the second control area.

Images