CN101930666A - Testing method of one-way road section traffic volume composition - Google Patents

Testing method of one-way road section traffic volume composition Download PDF

Info

Publication number
CN101930666A
CN101930666A CN 201010139939 CN201010139939A CN101930666A CN 101930666 A CN101930666 A CN 101930666A CN 201010139939 CN201010139939 CN 201010139939 CN 201010139939 A CN201010139939 A CN 201010139939A CN 101930666 A CN101930666 A CN 101930666A
Authority
CN
China
Prior art keywords
node
traffic
impedance
highway section
moment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 201010139939
Other languages
Chinese (zh)
Other versions
CN101930666B (en
Inventor
段进
翁芳玲
徐中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN2010101399398A priority Critical patent/CN101930666B/en
Publication of CN101930666A publication Critical patent/CN101930666A/en
Application granted granted Critical
Publication of CN101930666B publication Critical patent/CN101930666B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a testing method of one-way road section traffic volume composition comprising the following steps of initializing traffic demand; initializing city road network; initializing impedance, traffic volume and traffic volume composition of each one-way road section; initializing cycle index i which is equal to 1; judging whether i is not more than n; if the i is not more than n, executing an impedance matrix generation algorithm, calculating the current impedance matrix by using the impedance matrix generation algorithm and preparing for searching the shortest path in traffic allocation by using any one node as a source point and a preposing node relative to the other node, and if the i is more than 1 ending to obtain the traffic volume composition on each road section; allocating the traffic demand ODi matrix to a previous directed network Gi-1(V, E) by using a traffic allocation algorithm to form a current directed network Gi (V, E); recording the traffic volume composition on each one-way road section of the directed network Gi (V, E), wherein the traffic volume composition comprises a starting point, an ending point and volume; supposing i=i+1, continuously judging whether the i is not more than n, and circularly executing until all OD matrixes are allocated; accumulating all traffic volumes on the one-way road section according to the starting point and the ending point to form final one-way road section traffic volume composition which is represented in a form of table.

Description

The assay method that one-way road section traffic volume constitutes
Technical field
The present invention relates to the technical field of traffic forecast and evaluation method, relate in particular to the assay method that a kind of road section traffic volume constitutes.
Background technology
This method is the quadravalence section to traditional traffic forecast " quadravalence section " method---the improvement that traffic distributes.Traffic forecast " quadravalence section " method distributes four parts to constitute by plot generation and traffic attraction, OD distribution, mode of transportation division and traffic.This method is by the improvement to the conventional traffic assigning process, and the volume of traffic of having analyzed each highway section constitutes, and each volume of traffic by which plot produces, attracts and amount, represents with the form of form at last.
Summary of the invention
1, technical matters:
The invention provides the assay method that a kind of one-way road section traffic volume constitutes, the advantage that can provide each volume of traffic in the unidirectional highway section that origin and destination and amount thereof take place is provided.
2, technical scheme:
The assay method that a kind of one-way road section traffic volume constitutes:
Step 1: the initialization transport need on average is split as ∏ five equilibrium: OD with the OD matrix of expressing transport need 1Matrix, OD 2Matrix ..., OD Matrix, ∏ are the positive integer greater than 1, the OD matrix after each splits comprise Γ starting and terminal point to transport need, promptly
q i(o 1, d 1), q i(o 2, d 2) ..., q i(o n, d n) ..., q i(o Γ, d Γ), q wherein i(o n, d n) expression OD iStarting point is a node in the matrix
Figure GSA00000073706000011
Terminal point is a node
Figure GSA00000073706000012
Transport need,
Step 2: urban road network is set at resistive directed networks G, and (V, E), (V, node E) form set of node V, are designated as: V={v as described resistive directed networks G with the terminal in crossing in the urban road network and unidirectional highway section 1, v 2..., v τ, τ is the number of node, v 1, v 2..., v τBe the node of directed networks, (the limit collection is designated as: E={e for V, directed edge E) as described resistive directed networks G with the unidirectional highway section in the urban road network 1, e 2..., e σ..., e η, e 1, e 2..., e σ..., e ηThe unidirectional highway section of the directed edge representative of expression directed networks, σ is the subscript in a certain unidirectional highway section, and η is the number in unidirectional highway section in the directed networks, and (V E) carries out initialization, obtains initialization directed networks G to directed networks G 0(V, E),
Step 3: initialization impedance, the volume of traffic and the volume of traffic constitute, and make directed networks G 0(V, E) in the accumulative total volume of traffic on all unidirectional highway sections
Figure GSA00000073706000013
σ=1,2 ..., η, Represent unidirectional highway section e σOn the accumulative total volume of traffic; Make the volume of traffic on all unidirectional highway sections constitute
Figure GSA00000073706000015
σ=1,2 ..., η,
Figure GSA00000073706000016
Represent unidirectional highway section e σGo up by node v oTo v dThe volume of traffic, the impedance in unidirectional highway section is:
w e σ = w e σ 0 ( 1 + α ( x e σ c e σ ) β )
Wherein
Figure GSA00000073706000018
For when on the unidirectional highway section during without any vehicle, the time cost that vehicle passes through, it is that the geometrical length in known unidirectional highway section is divided by the known design speed of a motor vehicle;
Figure GSA00000073706000019
Be the one-way road section traffic volume amount,
Figure GSA000000737060000110
Be the maximum traffic capacity in unidirectional highway section,
Figure GSA000000737060000111
Be unidirectional highway section saturation degree; α and β are empirical parameter, α=0.15, and β=4.0,
Step 4: loop initialization pointer i=1,
Step 5: judge i≤∏, if execution in step 6 then; Otherwise, finishing, the volume of traffic that obtains on every highway section constitutes
Figure GSA000000737060000112
σ=1,2 ..., η, 1≤o≤τ, 1≤d≤τ,
Step 6: calculate current impedance matrix and be source point and with any node respectively with the impedance matrix generating algorithm, search shortest path for traffic in distributing and prepare with respect to the preposition node of another node arbitrarily,
Step 7: use the traffic allocation algorithm with transport need OD iMatrix allocation arrives last time directed networks G I-1(V, E) on, form current directed networks G i(V E), and records to network G i(V, E) volume of traffic on every the unidirectional highway section on constitutes,
Step 8: make i=i+1, return step 5,
Described impedance matrix generating algorithm is:
Step 6.1: loop initialization pointer j=1,
Step 6.2: if j≤τ, then execution in step 6.3, otherwise impedance matrix generates and to finish, and obtains the impedance between each node, and to obtain respectively with any node be source point and with respect to the preposition node of another node arbitrarily,
Step 6.3: adopt single source impedance generating algorithm to calculate arbitrary node v jArrive all node v as source point kImpedance r j[v k], k=1,2 ... τ, and obtain preposition node with respect to arbitrary node of source point,
Step 6.4: make j=j+1, return step 6.2,
Wherein, the described single source impedance generating algorithm of step 6.3 is:
Step 6.3.1: initialization p j[] array, r j[] array, S set, set Q,
p j[v k] represent with node v jNode v for source point kPreposition node, k=1,2 ... τ, the order:
p j[v k]=φ
r j[v k] represent with source point v jBe starting point, node v kBe the path impedance of terminal point, k=1,2 ... τ, order:
r j [ v k ] = 0 k = j ∝ k ≠ j
S set is used to deposit the node of having handled, and initialization S makes S=φ; Set Q is used to deposit untreated node, and initialization Q makes Q=V,
Step 6.3.2: judge that whether the Q set is empty, if be empty, then obtains node v jArrive all node v as source point kImpedance r j[v k], k=1,2 ..., τ, and obtain preposition node with respect to arbitrary node of source point; Otherwise, execution in step 6.3.3,
Step 6.3.3: establishing among the current set Q has Ψ node, and Ψ≤τ according to node subscript series arrangement from small to large, and corresponds respectively to l with tactic node subscript to the node among the current Q 1, l 2..., l Ψ, obtain
Figure GSA00000073706000022
1≤l ξ≤ τ,
Step 6.3.4: check
Figure GSA00000073706000023
Therefrom find out minimum And with node
Figure GSA00000073706000025
As node u,
Step 6.3.5: determine source point v respectively jTo the impedance of each adjacent node of node u, and definite source point v jThe preposition node of adjacent node of node u,
Step 6.3.6: node u is moved to the S set by the Q set, execution in step 6.3.2,
Wherein, the described definite respectively source point v of step 6.3.5 jTo the impedance of each adjacent node of node u, and definite source point v jThe method of preposition node of adjacent node of node u be:
Step 6.3.5.1: establishing node u has Ω adjacent node, and Ω≤τ according to subscript series arrangement from small to large, and corresponds respectively to m with tactic node subscript to the adjacent node of node u 1, m 2..., m Ω, obtain the adjacent node of tactic node u
v m 1 , v m 2 , . . . , v m ρ , . . . , v m Ω , 1 ≤ , m ρ ≤ τ ,
Step 6.3.5.2: establish cycle index g=1,
Step 6.3.5.3: if g≤Ω, execution in step 6.3.5.4 then, otherwise, current source point v obtained jTo the impedance of each adjacent node of node u, and source point v jThe preposition node of adjacent node of node u,
Step 6.3.5.4: if
Figure GSA00000073706000027
Execution in step 6.3.5.5 then; Otherwise, skip to step 6.3.5.7, wherein,
Figure GSA00000073706000028
Expression by node u to adjacent node
Figure GSA00000073706000029
The impedance of unidirectional highway section, if u is extremely
Figure GSA000000737060000210
Between many unidirectional highway sections are arranged, then
Figure GSA000000737060000211
Refer to the impedance on that unidirectional highway section of impedance minimum,
Step 6.3.5.5: make source point v jTo adjacent node
Figure GSA000000737060000212
Impedance equal source point v jTo the impedance of u add u with
Figure GSA000000737060000213
Impedance on the unidirectional highway section, even
Figure GSA000000737060000214
Step 6.3.5.6: order With respect to v jFor the preposition node of starting point is u, even
Figure GSA000000737060000216
Step 6.3.5.7: make g=g+1, return step 6.3.5.3;
Step 7 is described " uses the traffic allocation algorithm with transport need OD iMatrix allocation to directed networks G last time (V, E) on, form current directed networks G (V, E), and directed networks G before the record (V, E) volume of traffic on every the unidirectional highway section on constitutes " method be:
Step 7.1: establish transport need OD iComprise in the matrix Γ starting and terminal point to transport need, promptly
q i(o 1, d 1), q i(o 2, d 2) ..., q i(o n, d n) ..., q i(o Γ, d Γ), q wherein i(o n, d n) represent that starting point is a node in the ODi matrix
Figure GSA00000073706000031
Terminal point is a node Transport need,
Step 7.2: establish cycle index n=1,
Step 7.3: if n≤Γ, then execution in step 7.4; Otherwise, obtain current one-way road section traffic volume amount and road section traffic volume and form,
Step 7.4: with starting point-terminal point to allocation algorithm distribute the starting point that obtains by step 1
Figure GSA00000073706000033
Terminal point
Figure GSA00000073706000034
Between transport need q i(o n, d n),
Step 7.5: make n=n+1, return step 7.3,
Wherein, the starting point-terminal point described in the step 7.4 to allocation algorithm be:
Step 7.4.1: establish pointer s=d n, d nBe the node subscript of the terminal point of point-to-point transmission transport need,
Step 7.4.2: if pointer s=o n, then obtain transport need q i(o n, d n) in starting point be
Figure GSA00000073706000035
Terminal point is
Figure GSA00000073706000036
Shortest path on allocation result, otherwise, execution in step 7.4.3,
Step 7.4.3: find source point to be
Figure GSA00000073706000037
Node v sPreposition node
Figure GSA00000073706000038
Node v then sWith source point be Node v sPreposition node v ' sBetween the unidirectional highway section of unidirectional highway section impedance minimum be e a,
Step 7.4.4: at the unidirectional highway section e of unidirectional highway section impedance minimum aLast distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e a = x e a + q i ( o n , d n ) ,
Step 7.4.5: record is to the unidirectional highway section e of highway section impedance minimum aOn the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e a ( o n , d n ) = T e a ( o n , d n ) + q i ( o n , d n ) ,
Step 7.4.6: make pointer s=s ', return step 7.4.2 then.
3, technique effect:
(1) each volume of traffic starting point that list of missing articles distributes to the highway section in the conventional allocation method, the shortcoming of endpoint information have been overcome, the assay method that one-way road section traffic volume constitutes has write down the starting point of each volume of traffic that distributes, the point-to-point transmission of terminal point in assigning process, the traffic that forms the highway section constitutes.
(2) in the assigning process, this method has adopted the impedance matrix generating algorithm, and is more superior than traditional floyd-warshall algorithm.It is, under the prerequisite of given adjacency list, is τ with the time complexity of impedance matrix generating algorithm 0(η+τ log τ), and the time complexity of floyd-warshall algorithm is τ 3, the impedance matrix generating algorithm is more excellent.
(3) assay method that constitutes of one-way road section traffic volume provides a kind of system of new analysis road traffic for the urban planner.In fact, soil utilization and urban transportation are the relations of presenting mutually.Traditional traffic distributes only provides the total amount of one-way road section traffic volume, and the specifying information of these traffic can not be provided.This method provides starting point, the terminal point of each traffic in the unidirectional highway section, and the volume of traffic summation that can add up identical starting point and terminal point, by the record in the assigning process, three complete information that traffic constitutes in the unidirectional highway section are provided, make the urban planner can pass through these information, adjust generation, the traffic attraction in corresponding plot targetedly.For the urban planner has solved in land use planning, only provided the traffic distribution result in the past, and do not have adjust according to the difficult problem of goal directness.
Description of drawings
Fig. 1 is urban road network's exemplary plot.
Fig. 2 is main process figure.
Fig. 3 is the impedance matrix process flow diagram.
Fig. 4 is single source impedance product process figure.
Fig. 5 is the node processing algorithm flow chart.
Fig. 6 is traffic allocation flow figure.
Fig. 7 is a point-to-point transmission allocation algorithm process flow diagram.
Fig. 8 is urban road network's exemplary plot.
Fig. 9 is the figure of initialization urban road network.
Embodiment
The assay method main body that road section traffic volume constitutes: the assay method that road section traffic volume constitutes is based upon on follow-on method of traffic assignment.This method mainly is made of four steps, comprises that traffic generation and attraction, Traffic Distribution, mode of transportation are divided and traffic distributes, and it is mainly to the quadravalence section---and the traffic distribution improves.
General thought is: " initialization demand (fractionation) " → " initialization road network " → " the initialization volume of traffic " → " handling each part transport need successively ".To each part transport need: at first set up impedance matrix according to the accumulative total volume of traffic on road network and current each highway section; Carry out traffic again and distribute, with result and the stack of the accumulative total volume of traffic; Write down road section traffic volume simultaneously and constitute, constitute stack with the accumulative total road section traffic volume, as Fig. 2.
Embodiment 1
Step 1: the initialization transport need on average is split as ∏ five equilibrium: OD with the OD matrix of expressing transport need 1Matrix, OD 2Matrix ..., OD Matrix, ∏ are the positive integer greater than 1,
Step 2: the initialization road network, urban road network is set at resistive directed networks G, and (V, E), (V, node E) form set of node (V), are designated as: V={v as described resistive directed networks G with the terminal in crossing in the urban road network and highway section 1, v 2... v τ, τ is the number of node, v 1, v 2... v τBe the node of directed networks, (the limit collection is designated as: E={e for V, directed edge E) as described resistive directed networks G with the unidirectional highway section in the urban road network 1, e 2..., e a... e η, e 1, e 2..., e a... e ηThe limit of expression directed networks, a is the subscript in a certain highway section, and η is the number in highway section in the directed networks, and as Fig. 1, the calculating formula of impedance is:
w e a = w e a 0 ( 1 + α ( x e a c e a ) β )
Wherein For when on the highway section during without any vehicle, the time cost that vehicle passes through, it is relevant with geometrical length, the design speed in known highway section;
Figure GSA00000073706000043
Be the road section traffic volume amount,
Figure GSA00000073706000044
Be the maximum traffic capacity in highway section,
Figure GSA00000073706000045
Be the highway section saturation degree; α and β are empirical parameter, α=0.15, and β=4.0,
Step 3: the initialization volume of traffic and the volume of traffic constitute, and make the accumulative total volume of traffic on all highway sections
Figure GSA00000073706000046
Expression highway section e aOn the accumulative total volume of traffic; Make the volume of traffic on all highway sections constitute Expression highway section e aGo up by node v oTo v dThe volume of traffic,
Step 4: loop initialization pointer i=1,
Step 5: judge i≤∏, if execution in step 6 then; Otherwise, finish,
Step 6: calculate current impedance matrix with the impedance matrix generating algorithm, and arbitrary node point prepares for searching shortest path in the traffic distribution with respect to the preposition node of arbitrary node,
Step 7: at first, find out shortest path, use the traffic allocation algorithm with transport need OD by the result of step 6 iMatrix allocation is to road network, and the volume of traffic that writes down on every highway section constitutes,
Step 8: make i=i+1, return step 5.
The impedance matrix generating algorithm
One, the generating algorithm of impedance matrix
At first, be source point with each node, calculate the minimum impedance that they arrive all nodes; Again result of calculation is compiled and formed impedance matrix.As Fig. 3.
It is as follows that impedance matrix generates step:
Step 6.1: loop initialization pointer j=1,
Step 6.2: if j≤τ, then execution in step 6.3, otherwise, impedance matrix generates and finishes, obtain the impedance between each node, and to obtain each node being the preposition node of starting point with respect to arbitrary node, is the preposition node of starting point with respect to arbitrary node by obtaining each node, we can obtain internodal shortest path at the traffic allocated phase
Step 6.3: adopt single source impedance generating algorithm to calculate arbitrary node v jArrive all node v as source point kImpedance r j[v k], k=1,2 ... τ, and to obtain arbitrary node be the preposition node of starting point with respect to source point, prepares for later traffic distributes,
Step 6.4: make j=j+1, return step 6.2.
Two, single source impedance generating algorithm
What algorithm solved is the shortest route problem of individual node in the transportation network (source point) to all nodes.
The input of algorithm has comprised a resistive transportation network G, and one of them comes source node v j, j=1,2, ..., τ.We represent the set of all nodes among the G with V.Limit among each figure all is that two nodes are formed has order elements right.We represent the set on all limits with E, and the impedance on limit then by impedance function w definition, spends in the traffic assignment problem middle impedance w representative time.The definition of w is Bureau of Public Roads model (BPR), wherein w 0For when on the highway section during without any vehicle, the time cost that vehicle passes through, it is relevant with geometrical length, the design speed in known highway section;
Figure GSA00000073706000051
Be link flow,
Figure GSA00000073706000052
Be the highway section capacity,
Figure GSA00000073706000053
Be the highway section saturation degree; α and β are empirical parameter, α=0.15, β=4.0.
w ( e a ) = w 0 ( e a ) ( 1 + α ( x e a c e a ) β )
The impedance on limit is two time costs between the node.Appointing the impedance of inter-two-point path, is exactly the impedance summation in all highway sections on this path.Node v is arranged among the known V j, algorithm can find from a node v jTake office the what shortest path of his node.
This algorithm is by being each node v kKeep found so far from v jTo v kShortest path come work.When initial, source point v jThe path value be 0 (r by tax j[v j]=0), simultaneously the path r of every other node j[v k] be made as infinity, represent that promptly we do not know that any path of these nodes of leading to is (for all node v among the V kRemove v jOuter r j[v k]=∞); The father node of all nodes is composed and is empty (p simultaneously j[v k]=φ).When algorithm finishes, array r jWhat store in [] is from v jTo v kThe impedance of shortest path, if perhaps the non-existent words in path are infinitely great.
Concrete steps are as follows:
Step 6.3.1: initialization p j[] array, r j[] array, S set, set Q,
p j[v k] expression v kWith respect to v jBe the preposition node of starting point, k=1,2 ... τ, order:
p j[v k]=φ
r j[v k] the expression starting point is v j, terminal point is v kThe impedance of shortest path, k=1,2 ... τ, order:
r j [ v k ] = 0 k = j ∝ k ≠ j
S set is used to deposit the node of having handled, and initialization S makes S=φ; Set Q is used to deposit untreated node, and initialization Q makes Q=V,
Step 6.3.2: judge that whether the Q set is empty, if be empty, then finishes; Otherwise, execution in step 6.6.3,
Step 6.3.3: establishing among the set Q has Ψ node, Ψ≤τ.Node among the current Q is arranged from small to large according to subscript:
Figure GSA00000073706000056
L wherein ξBe the subscript of node v, 1≤l ξ≤ τ, l ξBe the little node subscript of node subscript ξ among the Q,
Step 6.3.4: check Therefrom find out minimum
Figure GSA00000073706000058
Order
Figure GSA00000073706000059
Step 6.3.5: the adjacent node of determining node u is with respect to source point v jPreposition node and source point v for starting point jTo the impedance of the adjacent node of u, prepare for obtaining shortest path in the traffic allocation algorithm hereinafter,
Step 6.3.6: node u is moved to S set, execution in step 6.3.2 by the Q set.
Wherein, the adjacent node of described definite node u is with respect to source point v jPreposition node and source point v for starting point jImpedance algorithms to the adjacent node of u is:
Step 6.3.5.1: establishing u has Ω adjacent node, Ω≤τ.Adjacent node to u is arranged from small to large according to subscript:
Figure GSA000000737060000510
M wherein gBe the subscript of node v, 1≤m g≤ τ, m gBe the little node subscript of subscript g in the adjacent node of u,
Step 6.3.5.2: establish cycle index g=1,
Step 6.3.5.3: if g≤Ω, execution in step 6.3.5.4 then, otherwise node processing finishes,
Step 6.3.5.4: if
Figure GSA000000737060000511
Execution in step 6.3.5.5 then; Otherwise, skip to step 6.3.5.7, wherein,
Figure GSA000000737060000512
Expression by node u to adjacent node
Figure GSA000000737060000513
The highway section impedance, if u is extremely
Figure GSA000000737060000514
Between many highway sections are arranged, then
Figure GSA000000737060000515
Refer to the impedance on that highway section of impedance minimum,
Step 6.3.5.5: make source point v jTo adjacent node
Figure GSA000000737060000516
Impedance equal source point v jTo the impedance of u add u with
Figure GSA000000737060000517
Impedance on the highway section, even
r j [ v m g ] = r j [ u ] + w u → v m g ,
Step 6.3.5.6: order
Figure GSA000000737060000519
With respect to v jFor the preposition node of starting point is u, even
Figure GSA000000737060000520
Step 6.3.5.7: make g=g+1, return step 6.3.5.3.
The traffic allocation algorithm
One, mathematical model
The apportion model commonly used of urban road network has: entirely by completely without model (All or Nothing), user equilibrium model (UserEquilibrium), polynary Probit model.Here we use the user equilibrium model to come the distribute traffic amount.The mathematic(al) representation of user equilibrium model is:
min z ( x ) = Σ a = 1 η ∫ 0 x e a w e a ( t ) dt
Two, arthmetic statement
Finding the solution of user equilibrium model usually is the Frank-Wolfe method of finding the solution the linear restriction quadratic programming by a kind of.Can obtain exact solution by this method, but the operation efficiency of this method is lower, is not suitable for engineering calculation.Therefore we adopt the increment apportion design to solve the user equilibrium problem in program.The increment distribution is meant OD demand average mark is slit into some parts, is loaded on the road network successively.The highway section can become along with the increase of the volume of traffic and block up, and the highway section also can increase by the time, and original time shortest path will change.Distribute when finishing, which bar road no matter the user of identical starting and terminal point pass through, and institute's time spent is identical, has just reached " user equilibrium ".
Concrete steps are as follows:
Step 7.1: establish transport need OD iComprise in the matrix Γ starting and terminal point to transport need, promptly
q i(o 1, d 1), q i(o 2, d 2) ..., q i(o n, d n) ..., q i(o Γ, d Γ), q wherein i(o n, d n) expression OD iStarting point is in the matrix Terminal point is Transport need,
Step 7.2: establish cycle index n=1,
Step 7.3: if n≤Γ, then execution in step 7.4; Otherwise, distribute and finish,
Step 7.4: handle with the point-to-point transmission allocation algorithm
Figure GSA00000073706000064
Between transport need q i(o n, d n),
Step 7.5: make n=n+1, return step 7.3.
Wherein, the point-to-point transmission allocation algorithm is:
Step 7.4.1: establish pointer s=d n, d nBe the node subscript of the terminal point of point-to-point transmission transport need,
Step 7.4.2: if the node subscript s=o of pointed starting point n, then finish, otherwise, execution in step 7.4.3,
Step 7.4.3: find node v sPreposition node
Figure GSA00000073706000065
And node v sWith its preposition node v ' sBetween highway section e a,
Step 7.4.4: at highway section e aLast distribute traffic amount, order accumulative total road section traffic volume amount
Figure GSA00000073706000066
Step 7.4.5: record highway section e aOn the volume of traffic constitute, the order accumulative total volume of traffic constitutes
Figure GSA00000073706000067
Step 7.4.6: make the subscript of its preposition node of pointer s sensing, even s=s ' returns step 7.4.2 then.
Embodiment 2
The 1st step: suppose that transport need OD matrix is
Table 1.OD matrix
O\D? v 1 v 2 v 3 v 4
v 1 0? 0? 100? 0?
v 2 0? 0? 0? 50?
v 3 0? 0? 0? 0?
v 4 0? 0? 0? 0?
It is split as 2 equal OD matrix: OD 1, OD 2,
OD matrix after table 2. splits
O\D? v 1 v 2 v 3 v 4
v 1 0? 0? 50? 0?
v 2 0? 0? 0? 25?
v 3 0? 0? 0? 0?
v 4 0? 0? 0? 0?
OD then 1={ q 1(1,3), q 1(2,4) }, q wherein 1(1,3)=50, q 1(2,4)=25,
OD 2={ q 2(1,3), q 2(2,4) }, q wherein 2(1,3)=50, q 2(2,4)=25,
The 2nd step: suppose urban road network such as Fig. 8,
Set of node V={v 1, v 2, v 3, v 4, node number τ=4, limit collection E={e 1, e 2, e 3, e 4Represent the unidirectional highway section of each bar in the road network, and limit number η=8, the one-way road section traffic volume capacity of establishing every limit is 60.
The 3rd step: make the one-way road section traffic volume amount
Figure GSA00000073706000071
One-way road section traffic volume is formed
Suppose that the zero flow impedance on the unidirectional highway section is
Figure 267186DEST_PATH_GSB00000294871700074
Figure 476768DEST_PATH_GSB00000294871700076
Because the volume of traffic on the unidirectional highway section is 0, so
Figure 112280DEST_PATH_GSB00000294871700077
Figure 379313DEST_PATH_GSB00000294871700078
Figure 578213DEST_PATH_GSB00000294871700079
Figure 302773DEST_PATH_GSB000002948717000711
Figure 57102DEST_PATH_GSB000002948717000712
Figure 794114DEST_PATH_GSB000002948717000713
Figure 130548DEST_PATH_GSB000002948717000714
As Fig. 9,
The 4th step: loop initialization pointer i=1,
The 5th step: judge i≤∏, this moment i=1, ∏=2, the judgement formula is 1≤2, correct judgment continues to carry out,
The 6th step: generate impedance matrix, loop initialization pointer j=1 at first,
The 7th step: judge j≤τ, this moment j=1, τ=4, the judgement formula is 1≤4, correct judgment continues to carry out,
The 8th step: use single source impedance method of formation to calculate the impedance r of v1 to all nodes 1[v k] and with respect to v 1Preposition node,
The 9th step: initialization p 1[] array, r 1[] array, S set, set Q,
The 10th step: p 1[v 1]=φ, p 1[v 2]=φ, p 1[v 3]=φ, p 1[v 4]=φ,
The 11st step: r 1[v 1]=0, r 1[v 2]=∞, r 1[v 3]=∞, r 1[v 4]=∞,
The 12nd step: S=φ, Q={v 1, v 2, v 3, v 4,
The 13rd step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 14th step: check r 1[v 1], r 1[v 2], r 1[v 3], r 1[v 4], be respectively 0.0, ∞, ∞, so ∞ is l λ=1, make u=v 1,
The 15th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 16th step: establish cycle index g=1,
The 17th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 18th step: judge
Figure GSA000000737060000715
This moment j=1, g=1, mg=2, r j[m g]=∞, u=v 1, r j[u]=0.0,
Figure GSA000000737060000716
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 19th step: order
Figure GSA000000737060000717
Be r 1[v 2]=0.0+1.0=1.0,
The 20th step: order
Figure GSA000000737060000718
Be p 1[v 2]=v 1,
The 21st step: make g=g+1, thus g=2,
The 22nd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 23rd step: judge
Figure GSA000000737060000719
This moment j=1, g=2, mg=4, r j[m g]=∞, u=v 1, r j[u]=0.0, Therefore, the judgement formula is ∞>2.0, and correct judgment then continues to carry out,
The 24th step: order
Figure GSA000000737060000721
Be r 1[v 4]=0.0+2.0=2.0,
The 25th step: order
Figure GSA000000737060000722
Be p 1[v 4]=v 1,
The 26th step: make g=g+1, thus g=3,
The 27th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 28th step: u is moved to S set, u=v by the Q set 1,
The 29th step: judge whether current set Q is empty, at this moment Q={v 2, v 3, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=2, l 2=3, l 3=4,
The 30th step: check r 1[v 2], r 1[v 3], r 1[v 4], be respectively 1.0, ∞, 2.0, so l λ=2, make u=v 2,
The 31st step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 32nd step: establish cycle index g=1,
The 33rd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 34th step: judge
Figure GSA00000073706000081
This moment j=1, g=1, m g=1, r j[m g]=0.0, u=v 2, r j[u]=1.0,
Figure GSA00000073706000082
Therefore, the judgement formula is 0.0>2.0, and next adjacent node is judged in misjudgment
The 35th step: make g=g+1, thus g=2,
The 36th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 37th step: judge
Figure GSA00000073706000083
This moment j=1, g=2, m g=3, r j[m g]=∞, u=v 2, r j[u]=1.0,
Figure GSA00000073706000084
Therefore, the judgement formula is ∞>4.0, and correct judgment then continues to carry out,
The 38th step: order
Figure GSA00000073706000085
Be r 1[v 3]=0.0+3.0=3.0,
The 39th step: order Be p 1[v 3]=v 2,
The 40th step: make g=g+1, thus g=3,
The 41st step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 42nd step: u is moved to S set, u=v by the Q set 2,
The 43rd step: judge whether current set Q is empty, at this moment Q={v 3, v 4, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=3, l 2=4,
The 44th step: check r 1[v 3], r 1[v 4], be respectively 4.0,2.0, so l λ=4, make u=v 4,
The 45th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 46th step: establish cycle index g=1,
The 47th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 48th step: judge
Figure GSA00000073706000087
This moment j=1, g=1, m g=1, r j[m g]=0.0, u=v 4, r j[u]=2.0, Therefore, the judgement formula is 0.0>4.0, and next adjacent node is judged in misjudgment
The 49th step: make g=g+1, thus g=2,
The 50th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 51st step: judge
Figure GSA00000073706000089
This moment j=1, g=2, m g=3, r j[m g]=4.0, u=v 4, r j[u]=2.0,
Figure GSA000000737060000810
Therefore, the judgement formula is 4.0>3.0, and correct judgment then continues to carry out,
The 52nd step: order
Figure GSA000000737060000811
Be r 1[v 3]=0.0+1.0=1.0,
The 53rd step: order
Figure GSA000000737060000812
Be p 1[v 3]=v 4,
The 54th step: make g=g+1, thus g=3,
The 55th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 56th step: u is moved to S set, u=v by the Q set 4,
The 57th step: judge whether current set Q is empty, at this moment Q={v 3, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=3,
The 58th step: check r 1[v 3], be respectively 3.0, so l λ=3, make u=v 3,
The 59th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 60th step: establish cycle index g=1,
The 61st step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 62nd step: judge
Figure GSA00000073706000091
This moment j=1, g=1, m g=2, r j[m g]=1.0, u=v 3, r j[u]=3.0,
Figure GSA00000073706000092
Therefore, the judgement formula is 1.0>6.0, and next adjacent node is judged in misjudgment
The 63rd step: make g=g+1, thus g=2,
The 64th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 65th step: judge
Figure GSA00000073706000093
This moment j=1, g=2, m g=4, r j[m g]=2.0, u=v 3, r j[u]=3.0,
Figure GSA00000073706000094
Therefore, the judgement formula is 2.0>4.0, and next adjacent node is judged in misjudgment
The 66th step: make g=g+1, thus g=3,
The 67th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 68th step: u is moved to S set, u=v by the Q set 3,
The 69th step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 70th step: judge j≤τ, this moment j=2, τ=4, the judgement formula is 2≤4, correct judgment continues to carry out,
The 71st step: use single source impedance method of formation to calculate v 2Impedance r to all nodes 2[v k] and with respect to v 2Preposition node,
The 72nd step: initialization p 2[] array, r 2[] array, S set, set Q,
The 73rd step: p 2[v 1]=φ, p 2[v 2]=φ, p 2[v 3]=φ, p 2[v 4]=φ,
The 74th step: r 2[v 1]=∞, r 2[v 2]=0, r 2[v 3]=∞, r 2[v 4]=∞,
The 75th step: S=φ, Q=[v 1, v 2, v 3, v 4,
The 76th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 77th step: check r 2[v 1], r 2[v 2], r 2[v 3], r 2[v 4], be respectively ∞, 0.0, ∞, so ∞ is l λ=2, make u=v 2,
The 78th step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript from little not big series arrangement, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 79th step: establish cycle index g=1,
The 80th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 81st step: judge
Figure GSA00000073706000095
This moment j=2, g=1, m g=1, r j[m g]=∞, u=v 2, r j[u]=0.0,
Figure GSA00000073706000096
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 82nd step: order
Figure GSA00000073706000097
Be r 2[v 1]=0.0+1.0=1.0,
The 83rd step: order
Figure GSA00000073706000098
Be p 2[v 1]=v 2,
The 84th step: make g=g+1, thus g=2,
The 85th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 86th step: judge
Figure GSA00000073706000099
This moment j=2, g=2, m g=3, r j[m g]=∞, u=v 2, r j[u]=0.0,
Figure GSA000000737060000910
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 87th step: order
Figure GSA000000737060000911
Be r 2[v 3]=0.0+3.0=3.0,
The 88th step: order
Figure GSA000000737060000912
Be p 2[v 3]=v 2,
The 89th step: make g=g+1, thus g=3,
The 90th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 91st step: u is moved to S set, u=v by the Q set 2,
The 92nd step: judge whether current set Q is empty, at this moment Q={v 1, v 3, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=3, l 3=4,
The 93rd step: check r 2[v 1], r 2[v 3], r 2[v 4], being respectively 1.0,3.0, so ∞ is l λ=1, make u=v 1,
The 94th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 95th step: establish cycle index g=1,
The 96th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 97th step: judge
Figure GSA00000073706000101
This moment j=2, g=1, m g=2, r j[m g]=0.0, u=v 1, r j[u]=1.0,
Figure GSA00000073706000102
Therefore, the judgement formula is 0.0>2.0, and next adjacent node is judged in misjudgment
The 98th step: make g=g+1, thus g=2,
The 99th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 100th step: judge
Figure GSA00000073706000103
This moment j=2, g=2, m g=4, r j[m g]=∞, u=v 1, r j[u]=1.0,
Figure GSA00000073706000104
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 101st step: order
Figure GSA00000073706000105
Be r 2[v 4]=0.0+2.0=2.0,
The 102nd step: order
Figure GSA00000073706000106
Be p 2[v 4]=v 1,
The 103rd step: make g=g+1, thus g=3,
The 104th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then looked into the distance out in misjudgment,
The 105th step: u is moved to S set, u=v by the Q set 1,
The 106th step: judge whether current set Q is empty, at this moment Q={v 3, v 4, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=3, l 2=4,
The 107th step: check r 2[v 3], r 2[v 4], be respectively 3.0,3.0, so l λ=3, make u=v 3,
The 108th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 109th step: establish cycle index g=1,
The 110th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 111st step: judge This moment j=2, g=1, m g=2, r j[m g]=0.0, u=v 3, r j[u]=3.0,
Figure GSA00000073706000108
Therefore, the judgement formula is 0.0>6.0, and next adjacent node is judged in misjudgment
The 112nd step: make g=g+1, thus g=2,
The 113rd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 114th step: judge
Figure GSA00000073706000109
This moment j=2, g=2, m g=4, r j[m g]=3.0, u=v 3, r j[u]=3.0,
Figure GSA000000737060001010
Therefore, the judgement formula is 3.0>4.0, and next adjacent node is judged in misjudgment
The 115th step: make g=g+1, thus g=3,
The 116th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 117th step: u is moved to S set, u=v by the Q set 3,
The 118th step: judge whether current set Q is empty, at this moment Q={v 4, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=4,
The 119th step: check r 2[v 4], be respectively 3.0, so l λ=4, make u=v 4,
The 120th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 121st step: establish cycle index g=1,
The 122nd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 123rd step: judge
Figure GSA000000737060001011
This moment j=2, g=1, m g=1, r j[m g]=1.0, u=v 4, r j[u]=3.0,
Figure GSA00000073706000111
Therefore, the judgement formula is 1.0>5.0, and next adjacent node is judged in misjudgment
The 124th step: make g=g+1, thus g=2,
The 125th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 126th step: judge
Figure GSA00000073706000112
This moment j=2, g=2, m g=3, r j[m g]=3.0, u=v 4, r j[u]=3.0, Therefore, the judgement formula is 3.0>4.0, and next adjacent node is judged in misjudgment
The 127th step: make g=g+1, thus g=3,
The 128th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 129th step: u is moved to S set, u=v by the Q set 4,
The 130th step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 131st step: judge j≤τ, this moment j=3, τ=4, the judgement formula is 3≤4, correct judgment continues to carry out,
The 132nd step: use single source impedance method of formation to calculate v 3Impedance r to all nodes 3[v k] and with respect to v 3Preposition node,
The 133rd step: initialization p 3[] array, r 3[] array, S set, set Q,
The 134th step: p 3[v 1]=φ, p 3[v 2] a=φ, p 3[v 3]=φ, p 3[v 4]=φ,
The 135th step: r 3[v 1]=∞, r 3[v 2]=∞, r 3[v 3]=0, r 3[v 4]=∞,
The 136th step: S=φ, Q={v 1, v 2, v 3, v 4,
The 137th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 138th step: check r 3[v 1], r 3[v 2], r 3[v 3], r 3[v 4], be respectively ∞, ∞, 0.0, so ∞ is l λ=3, make u=v 3,
The 139th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 140th step: establish cycle index g=1,
The 141st step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 142nd step: judge
Figure GSA00000073706000114
This moment j=3, g=1, m g=2, r j[m g]=∞, u=v 3, r j[u]=0.0,
Figure GSA00000073706000115
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 143rd step: order
Figure GSA00000073706000116
Be r 3[v 2]=0.0+3.0=3.0,
The 144th step: order
Figure GSA00000073706000117
Be p 3[v 2]=v 3,
The 145th step: make g=g+1, thus g=2,
The 146th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 147th step: judge
Figure GSA00000073706000118
This moment j=3, g=2, m g=4, r j[m g]=∞, u=v 3, r j[u]=0.0,
Figure GSA00000073706000119
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 148th step: order
Figure GSA000000737060001110
Be r 3[v 4]=0.0+1.0=1.0,
The 149th step: order Be p 3[v 4]=v 3,
The 150th step: make g=g+1, thus g=3,
The 151st step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 152nd step: u is moved to S set, u=v by the Q set 3,
The 153rd step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=4,
The 154th step: check r 3[v 1], r 3[v 2], r 3[v 4], be respectively ∞, 3.0,1.0, so l λ=4, make u=v 4,
The 155th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 156th step: establish cycle index g=1,
The 157th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 158th step: judge This moment j=3, g=1, m g=1, r j[m g]=∞, u=v 4, r j[u]=1.0, Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 159th step: order
Figure GSA00000073706000123
Be r 3[v 1]=0.0+2.0=2.0,
The 160th step: order
Figure GSA00000073706000124
Be p 3[v 1]=v 4,
The 161st step: make g=g+1, thus g=2,
The 162nd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 163rd step: judge This moment j=3, g=2, m g=3, r j[m g]=0.0, u=v 4, r j[u]=1.0,
Figure GSA00000073706000126
Therefore, the judgement formula is 0.0>2.0, and next adjacent node is judged in misjudgment
The 164th step: make g=g+1, thus g=3,
The 165th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 166th step: u is moved to S set, u=v by the Q set 4,
The 167th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2,
The 168th step: check r 3[v 1], r 3[v 2], be respectively 3.0,3.0, so l λ=1, make u=v 1,
The 169th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 170th step: establish cycle index g=1,
The 171st step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 172nd step: judge
Figure GSA00000073706000127
This moment j=3, g=1, m g=2, r j[m g]=3.0, u=v 1, r j[u]=3.0, Therefore, the judgement formula is 3.0>4.0, and next adjacent node is judged in misjudgment
The 173rd step: make g=g+1, thus g=2,
The 174th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 175th step: judge
Figure GSA00000073706000129
This moment j=3, g=2, m g=4, r j[m g]=1.0, u=v 1, r j[u]=3.0,
Figure GSA000000737060001210
Therefore, the judgement formula is 1.0>5.0, and next adjacent node is judged in misjudgment
The 176th step: make g=g+1, thus g=3,
The 177th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 178th step: u is moved to S set, u=v by the Q set 1,
The 179th step: judge whether current set Q is empty, at this moment Q={v 2, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=2,
The 180th step: check r 3[v 2], be respectively 3.0, so l λ=2, make u=u 2,
The 181st step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 182nd step: establish cycle index g=1,
The 183rd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 184th step: judge
Figure GSA000000737060001211
This moment j=3, g=1, m g=1, r j[m g]=3.0, u=v 2, r j[u]=3.0, Therefore, the judgement formula is 3.0>4.0, and next adjacent node is judged in misjudgment
The 185th step: make g=g+1, thus g=2,
The 186th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 187th step: judge
Figure GSA00000073706000131
This moment j=3, g=2, m g=3, r j[m g]=0.0, u=v 2, r j[u]=3.0, Therefore, the judgement formula is 0.0>6.0, and next adjacent node is judged in misjudgment
The 188th step: make g=g+1, thus g=3,
The 189th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 190th step: u is moved to S set, u=v by the Q set 2,
The 191st step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 192nd step: judge j≤τ, this moment j=4, τ=4, the judgement formula is 4≤4, correct judgment continues to carry out,
The 193rd step: use single source impedance method of formation to calculate v 4Impedance r to all nodes 4[v k] and with respect to v 4Preposition node,
The 194th step: initialization p 4[] array, r 4[] array, S set, set Q,
The 195th step: p 4[v 1]=φ, p 4[v 2]=φ, p 4[v 3]=φ, p 4[v 4]=φ,
The 196th step: r 4[v 1]=∞, r 4[v 2]=∞, r 4[v 3]=∞, r 4[v 4]=0,
The 197th step: S=φ, Q={v 1, v 2, v 3, v 4,
The 198th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 199th step: check r 4[v 1], r 4[v 2], r 4[v 3], r 4[v 4], be respectively ∞, ∞, ∞, 0.0, so l λ=4, make u=v 4,
The 200th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 201st step: establish cycle index g=1,
The 202nd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 203rd step: judge
Figure GSA00000073706000133
This moment j=4, g=1, m g=1, r j[m g]=∞, u=v 4, r j[u]=0.0,
Figure GSA00000073706000134
Therefore, the judgement formula is ∞>2.0, and correct judgment then continues to carry out,
The 204th step: order
Figure GSA00000073706000135
Be r 4[v 1]=0.0+2.0=2.0,
The 205th step: order
Figure GSA00000073706000136
Be p 4[v 1]=v 4,
The 206th step: make g=g+1, thus g=2,
The 207th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 208th step: judge
Figure GSA00000073706000137
This moment j=4, g=2, m g=3, r j[m g]=∞, u=v 4, r j[u]=0.0,
Figure GSA00000073706000138
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 209th step: order
Figure GSA00000073706000139
Be r 4[v 3]=0.0+1.0=1.0,
The 210th step: order
Figure GSA000000737060001310
Be p 4[v 3]=v 4,
The 211st step: make g=g+1, thus g=3,
The 212nd step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 213rd step: u is moved to S set, u=v by the Q set 4,
The 214th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3,
The 215th step: check r 4[v 1], r 4[v 2] a, r 4[v 3], be respectively 2.0, ∞, 1.0, so l λ=3, make u=v 3,
The 216th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 217th step: establish cycle index g=1,
The 218th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 219th step: judge
Figure GSA00000073706000141
This moment j=4, g=1, m g=2, r j[m g]=∞, u=v 3, r j[u]=1.0,
Figure GSA00000073706000142
Therefore, the judgement formula is ∞>4.0, and correct judgment then continues to carry out,
The 220th step: order
Figure GSA00000073706000143
Be r 4[v 2]=0.0+3.0=3.0,
The 221st step: order
Figure GSA00000073706000144
Be p 4[v 2]=v 3,
The 222nd step: make g=g+1, thus g=2,
The 223rd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 224th step: judge
Figure GSA00000073706000145
This moment j=4, g=2, m g=4, r j[m g]=0.0, u=v 3, r j[u]=1.0,
Figure GSA00000073706000146
Therefore, the judgement formula is 0.0>2.0, and next adjacent node is judged in misjudgment
The 225th step: make g=g+1, thus g=3,
The 226th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 227th step: u is moved to S set, u=v by the Q set 3,
The 228th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2,
The 229th step: check r 4[v 1], r 4[v 2], be respectively 2.0,4.0, so l λ=1, make u=v 1,
The 230th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 231st step: establish cycle index g=1,
The 232nd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 233rd step: judge
Figure GSA00000073706000147
This moment j=4, g=1, m g=2, r j[m g]=4.0, u=v 1, r j[u]=2.0,
Figure GSA00000073706000148
Therefore, the judgement formula is 4.0>3.0, and correct judgment then continues to carry out,
The 234th step: order
Figure GSA00000073706000149
Be r 4[v 2]=0.0+1.0=1.0,
The 235th step: order Be p 4[v 2]=v 1,
The 236th step: make g=g+1, thus g=2,
The 237th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 238th step: judge
Figure GSA000000737060001411
This moment j=4, g=2, m g=4, r j[m g]=0.0, u=v 1, r j[u]=2.0,
Figure GSA000000737060001412
Therefore, the judgement formula is 0.0>4.0, and next adjacent node is judged in misjudgment
The 239th step: make g=g+1, thus g=3,
The 240th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 241st step: u is moved to S set, u=v by the Q set 1,
The 242nd step: judge whether current set Q is empty, at this moment Q={v 2, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=2,
The 243rd step: check r 4[v 2], be respectively 3.0, so l λ=2, make u=v 2,
The 244th step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 245th step: establish cycle index g=1,
The 246th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 247th step: judge
Figure GSA000000737060001413
This moment j=4, g=1, m g=1, r j[m g]=2.0, u=v 2, r j[u]=3.0,
Figure GSA00000073706000151
Therefore, the judgement formula is 2.0>4.0, and next adjacent node is judged in misjudgment
The 248th step: make g=g+1, thus g=2,
The 249th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 250th step: judge
Figure GSA00000073706000152
This moment j=4, g=2, m g=3, r j[m g]=1.0, u=v 2, r j[u]=3.0,
Figure GSA00000073706000153
Therefore, the judgement formula is 1.0>6.0, and next adjacent node is judged in misjudgment
The 251st step: make g=g+1, thus g=3,
The 252nd step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 253rd step: u is moved to S set, u=v by the Q set 2,
The 254th step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 255th step: establish cycle index n=1,
The 256th step: judge n≤Γ, this moment n=1, Γ=2, correct judgment then continues to carry out,
The 257th step: establish pointer s=d n, n=1, d n=3, so s=3
The 258th step: judge s=o n, this moment n=1, s=3, o n=1, misjudgment then continues to carry out,
The 259th step: finding source point is v 1Node v 3Preposition node p 1[v 3]=v 4, node v then 3With node v 4Between the unidirectional highway section of unidirectional highway section impedance minimum be e 6,
The 260th step: at the unidirectional highway section e of unidirectional highway section impedance minimum 6Last distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e 6 = x e 6 + 100.00 = 100.00 ,
The 261st step: record is to the unidirectional highway section e of highway section impedance minimum 6On the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e 6 ( 1,3 ) = T e 6 ( 1,3 ) + 100.00 = 100.00 ,
The 262nd step: make pointer s=s '=4,
The 263rd step: judge s=o n, this moment n=1, s=4, o n=1, misjudgment then continues to carry out,
The 264th step: finding source point is v 1Node v 4Preposition node p 1[v 4]=v 1, node v then 4With node v 1Between the unidirectional highway section of unidirectional highway section impedance minimum be e 8,
The 265th step: at the unidirectional highway section e of unidirectional highway section impedance minimum 8Last distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e 8 = x e 8 + 100.00 = 100.00 ,
The 266th step: record is to the unidirectional highway section e of highway section impedance minimum 8On the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e 8 ( 1,3 ) = T e 8 ( 1,3 ) + 100.00 = 100.00 ,
The 267th step: make pointer s=s '=1,
The 268th step: judge s=o n, this moment n=1, s=1, o n=1, correct judgment is then jumped out circulation,
The 269th step: make n=n+1, thus n=2,
The 270th step: judge n≤Γ, this moment n=2, Γ=2, correct judgment then continues to carry out,
The 271st step: establish pointer s=d n, n=2, d n=4, so s=4
The 272nd step: judge s=o n, this moment n=2, s=4, o n=2, misjudgment then continues to carry out,
The 273rd step: finding source point is v 2Node v 4Preposition node p 2[v 4]=v 1, node v then 4With node v 1Between the unidirectional highway section of unidirectional highway section impedance minimum be e 8,
The 274th step: at the unidirectional highway section e of unidirectional highway section impedance minimum 8Last distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e 8 = x e 8 + 50.00 = 150.00 ,
The 275th step: record is to the unidirectional highway section e of highway section impedance minimum 8On the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e 8 ( 2,4 ) = T e 8 ( 2,4 ) + 50.00 = 50.00 ,
The 276th step: make pointer s=s '=1,
The 277th step: judge s=o n, this moment n=2, s=1, o n=2, misjudgment then continues to carry out,
The 278th step: finding source point is v 2Node v 1Preposition node p 2[v 1]=v 2, node v then 1With node v 2Between the unidirectional highway section of unidirectional highway section impedance minimum be e 2,
The 279th step: at the unidirectional highway section e of unidirectional highway section impedance minimum 2Last distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e 2 = x e 2 + 50.00 = 50.00 ,
The 280th step: record is to the unidirectional highway section e of highway section impedance minimum 2On the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e 2 ( 2,4 ) = T e 2 ( 2,4 ) + 50.00 = 50.00 ,
The 281st step: make pointer s=s '=2,
The 282nd step: judge s=o n, this moment n=2, s=2, o n=2, correct judgment is then jumped out circulation,
The 283rd step: make n=n+1, thus n=3,
The 284th step: judge n≤Γ, this moment n=3, Γ=2, circulation is then jumped out in misjudgment,
The 285th step: judge i≤∏, this moment i=2, ∏=2, the judgement formula is 2≤2, correct judgment continues to carry out,
The 286th step: generate impedance matrix, loop initialization pointer j=1 at first,
The 287th step: judge j≤τ, this moment j=1, τ=4, the judgement formula is 1≤4, correct judgment continues to carry out,
The 288th step: use single source impedance method of formation to calculate v 1Impedance v to all nodes 1[v k] and with respect to v 1Preposition node,
The 289th step: initialization p 1[] array, r 1[] array, S set, set Q,
The 290th step: p 1[v 1]=φ, p 1[v 2]=φ, p 1[v 3]=φ, p 1[v 4]=φ,
The 291st step: r 1[v 1]=0, r 1[v 2]=∞, r 1[v 3]=∞, r 1[v 4]=∞,
The 292nd step: S=φ, Q={v 1, v 2, v 3, v 4,
The 293rd step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 294th step: check r 1[v 1], r 1[v 2], r 1[v 3], r 1[v 4], be respectively 0.0, ∞, ∞, so ∞ is l λ=1, make u=v 1,
The 295th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 296th step: establish cycle index g=1,
The 297th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 298th step: judge
Figure GSA00000073706000163
This moment j=1, g=1, m g=2, r j[m g]=∞, u=v 1, r j[u]=0.0,
Figure GSA00000073706000164
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 299th step: order
Figure GSA00000073706000165
Be r 1[v 2]=0.0+1.0=1.0,
The 300th step: order
Figure GSA00000073706000166
Be p 1[v 2]=v 1,
The 301st step: make g=g+1, thus g=2,
The 302nd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 303rd step: judge
Figure GSA00000073706000167
This moment j=1, g=2, m g=4, r j[m g]=∞, u=v 1, r j[u]=0.0,
Figure GSA00000073706000168
Therefore, the judgement formula is ∞>13.7, and correct judgment then continues to carry out,
The 304th step: order
Figure GSA00000073706000169
Be r 1[v 4]=0.0+13.7=13.7,
The 305th step: order
Figure GSA000000737060001610
Be p 1[v 4]=v 1,
The 306th step: make g=g+1, thus g=3,
The 307th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 308th step: u is moved to S set, u=v by the Q set 1,
The 309th step: judge whether current set Q is empty, at this moment Q={v 2, v 3, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=2, l 2=3, l 3=4,
The 310th step: check r 1[v 2], r 1[v 3], r 1[v 4], be respectively 1.0, ∞, 13.7, so l λ=2, make u=v 2,
The 311st step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 312nd step: establish cycle index g=1,
The 313rd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 314th step: judge
Figure GSA00000073706000171
This moment j=1, g=1, m g=1, r j[m g]=0.0, u=v 2, r j[u]=1.0,
Figure GSA00000073706000172
Therefore, the judgement formula is 0.0>2.1, and next adjacent node is judged in misjudgment
The 315th step: make g=g+1, thus g=2,
The 316th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 317th step: judge This moment j=1, g=2, m g=3, r j[m g]=∞, u=v 2, r j[u]=1.0,
Figure GSA00000073706000174
Therefore, the judgement formula is ∞>4.0, and correct judgment then continues to carry out,
The 318th step: order
Figure GSA00000073706000175
Be r 1[v 3]=0.0+3.0=3.0,
The 319th step: order Be p 1[v 3]=v 2,
The 320th step: make g=g+1, thus g=3,
The 321st step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 322nd step: u is moved to S set, u=v by the Q set 2,
The 323rd step: judge whether current set Q is empty, at this moment Q={v 3, v 4, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=3, l 2=4,
The 324th step: check r 1[v 3], r 1[v 4], be respectively 4.0,13.7, so l λ=3, make u=v 3,
The 325th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 326th step: establish cycle index g=1,
The 327th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 328th step: judge
Figure GSA00000073706000177
This moment j=1, g=1, m g=2, r j[m g]=1.0, u=v 3, r j[u]=4.0,
Figure GSA00000073706000178
Therefore, the judgement formula is 1.0>7.0, and next adjacent node is judged in misjudgment
The 329th step: make g=g+1, thus g=2,
The 330th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 331st step: judge This moment j=1, g=2, m g=4, r j[m g]=13.7, u=v 3, r j[u]=4.0,
Figure GSA000000737060001710
Therefore, the judgement formula is 13.7>5.0, and correct judgment then continues to carry out,
The 332nd step: order
Figure GSA000000737060001711
Be r 1[v 4]=0.0+1.0=1.0,
The 333rd step: order
Figure GSA000000737060001712
Be p 1[v 4]=v 3,
The 334th step: make g=g+1, thus g=3,
The 335th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 336th step: u is moved to S set, u=v by the Q set 3,
The 337th step: judge whether current set Q is empty, at this moment Q={v 4, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=4,
The 338th step: check r 1[v 4], be respectively 5.0, so l λ=4, make u=v 4,
The 339th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 340th step: establish cycle index g=1,
The 341st step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 342nd step: judge
Figure GSA000000737060001713
This moment j=1, g=1, m g=1, r j[m g]=0.0, u=v 4, v j[u]=5.0,
Figure GSA00000073706000181
Therefore, the judgement formula is 0.0>7.0, and next adjacent node is judged in misjudgment
The 343rd step: make g=g+1, thus g=2,
The 344th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 345th step: judge
Figure GSA00000073706000182
This moment j=1, g=2, m g=3, r j[m g]=4.0, u=v 4, r j[u]=5.0, Therefore, the judgement formula is 4.0>7.2, and next adjacent node is judged in misjudgment
The 346th step: make g=g+1, thus g=3,
The 347th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 348th step: u is moved to S set, u=v by the Q set 4,
The 349th step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 350th step: judge j≤τ, this moment j=2, τ=4, the judgement formula is 2≤4, correct judgment continues to carry out,
The 351st step: use single source impedance method of formation to calculate v 2Impedance r to all nodes 2[v k] and with respect to v 2Preposition node,
The 352nd step: initialization p 2[] array, r 2[] array, S set, set Q,
The 353rd step: p 2[v 1]=φ, p 2[v 2]=φ, p 2[v 3]=φ, p 2[v 4]=φ,
The 354th step: r 2[v 1]=∞, r 2[v 2]=0, r 2[v 3]=∞, r 2[v 4]=∞,
The 355th step: S=φ, Q={v 1, v 2, v 3, v 4,
The 356th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 357th step: check r 2[v 1], r 2[v 2], r 2[v 3], r 2[v 4], be respectively ∞, 0.0, ∞, so ∞ is l λ=2, make u=v 2,
The 358th step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 359th step: establish cycle index g=1,
The 360th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 361st step: judge
Figure GSA00000073706000184
This moment j=2, g=1, m g=1, r j[m g]=∞, u=v 2, r j[u]=0.0,
Figure GSA00000073706000185
Therefore, the judgement formula is ∞>1.1, and correct judgment then continues to carry out,
The 362nd step: order
Figure GSA00000073706000186
Be r 2[v 1]=0.0+1.1=1.1,
The 363rd step: order
Figure GSA00000073706000187
Be p 2[v 1]=v 2,
The 364th step: make g=g+1, thus g=2,
The 365th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 366th step: judge This moment j=2, g=2, m g=3, r j[m g]=∞, u=v 2, r j[u]=0.0,
Figure GSA00000073706000189
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 367th step: order
Figure GSA000000737060001810
Be r 2[v 3]=0.0+3.0=3.0,
The 368th step: order
Figure GSA000000737060001811
Be p 2[v 3]=v 2,
The 369th step: make g=g+1, thus g=3,
The 370th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 371st step: u is moved to S set, u=v by the Q set 2,
The 372nd step: judge whether current set Q is empty, at this moment Q={v 1, v 3, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=3, l 3=4,
The 373rd step: check r 2[v 1], r 2[v 3], r 2[v 4], being respectively 1.1,3.0, so ∞ is l λ=1, make u=v 1,
The 374th step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 375th step: establish cycle index g=1,
The 376th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 377th step: judge
Figure GSA00000073706000191
This moment j=2, g=1, m g=2, r j[m g]=0.0, u=v 1, r j[u]=1.1, Therefore, the judgement formula is 0.0>2.1, and next adjacent node is judged in misjudgment
The 378th step: make g=g+1, thus g=2,
The 379th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 380th step: judge
Figure GSA00000073706000193
This moment j=2, g=2, m g=4, r j[m g]=∞, u=v 1, r j[u]=1.1,
Figure GSA00000073706000194
Therefore, the judgement formula is ∞>14.8, and correct judgment then continues to carry out,
The 381st step: order
Figure GSA00000073706000195
Be r 2[v 4]=0.0+13.7=13.7,
The 382nd step: order
Figure GSA00000073706000196
Be p 2[v 4]=v 1,
The 383rd step: make g=g+1, thus g=3,
The 384th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 385th step: u is moved to S set, u=v by the Q set 1,
The 386th step: judge whether current set Q is empty, at this moment Q={v 3, v 4, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=3, l 2=4,
The 387th step: check r 2[v 3], r 2[v 4], be respectively 3.0,14.8, so l λ=3, make u=v 3,
The 388th step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 389th step: establish cycle index g=1,
The 390th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 391st step: judge
Figure GSA00000073706000197
This moment j=2, g=1, m g=2, r j[m g]=0.0, u=v 3, r j[u]=3.0,
Figure GSA00000073706000198
Therefore, the judgement formula is 0.0>6.0, and next adjacent node is judged in misjudgment
The 392nd step: make g=g+1, thus g=2,
The 393rd step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 394th step: judge
Figure GSA00000073706000199
This moment j=2, g=2, m g=4, r j[m g]=14.8, u=v 3, r j[u]=3.0,
Figure GSA000000737060001910
Therefore, the judgement formula is 14.8>4.0, and correct judgment then continues to carry out,
The 395th step: order Be r 2[v 4]=0.0+1.0=1.0,
The 396th step: order
Figure GSA000000737060001912
Be p 2[v 4]=v 3,
The 397th step: make g=g+1, thus g=3,
The 398th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 399th step: u is moved to S set, u=v by the Q set 3,
The 400th step: judge whether current set Q is empty, at this moment Q={v 4, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=4,
The 401st step: check r 2[v 4], be respectively 4.0, so l λ=4, make u=v 4,
The 402nd step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 403rd step: establish cycle index g=1,
The 404th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 405th step: judge
Figure GSA00000073706000201
This moment j=2, g=1, m g=1, r j[m g]=1.1, u=v 4, r j[u]=4.0,
Figure GSA00000073706000202
Therefore, the judgement formula is 1.1>6.0, and next adjacent node is judged in misjudgment
The 406th step: make g=g+1, thus g=2,
The 407th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 408th step: judge
Figure GSA00000073706000203
This moment j=2, g=2, m g=3, r j[m g]=3.0, u=v 4, r j[u]=4.0, Therefore, the judgement formula is 3.0>6.2, and next adjacent node is judged in misjudgment
The 409th step: make g=g+1, thus g=3,
The 410th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 411st step: u is moved to S set, u=v by the Q set 4,
The 412nd step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 413rd step: judge j≤τ, this moment j=3, τ=4, the judgement formula is 3≤4, correct judgment continues to carry out,
The 414th step: use single source impedance method of formation to calculate v 3Impedance r to all nodes 3[v k] and with respect to v 3Preposition node,
The 415th step: initialization p 3[] array, r 3[] array, S set, set Q,
The 416th step: p 3[v 1]=φ, p 3[v 2]=φ, p 3[v 3]=φ, p 3[v 4]=φ,
The 417th step: r 3[v 1]=∞, r 3[v 2]=∞, r 3[v 3]=0, r 3[v 4]=∞,
The 418th step: S=φ, Q={v 1, v 2, v 3, v 4,
The 419th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 420th step: check r 3[v 1], r 3[v 2], r 3[v 3], r 3[v 4], be respectively ∞, ∞, 0.0, so ∞ is l λ=3, make u=v 3,
The 421st step: establishing node u has Ω adjacent node, u=v 3, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 422nd step: establish cycle index g=1,
The 423rd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 424th step: judge
Figure GSA00000073706000205
This moment j=3, g=1, m g=2, r j[m g]=∞, u=v 3, r j[u]=0.0,
Figure GSA00000073706000206
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 425th step: order
Figure GSA00000073706000207
Be r 3[v 2]=0.0+3.0=3.0,
The 426th step: order
Figure GSA00000073706000208
Be p 3[v 2]=v 3,
The 427th step: make g=g+1, thus g=2,
The 428th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 429th step: judge
Figure GSA00000073706000209
This moment j=3, g=2, m g=4, r j[m g]=∞, u=v 3, r j[u]=0.0,
Figure GSA000000737060002010
Therefore, the judgement formula is ∞>1.0, and correct judgment then continues to carry out,
The 430th step: order
Figure GSA000000737060002011
Be r 3[v 4]=0.0+1.0=1.0,
The 431st step: order
Figure GSA000000737060002012
Be p 3[v 4]=v 3,
The 432nd step: make g=g+1, thus g=3,
The 433rd step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 434th step: u is moved to S set, u=v by the Q set 3,
The 435th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 4, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=4,
The 436th step: check r 3[v 1], r 3[v 2], r 3[v 4], be respectively ∞, 3.0,1.0, so l λ=4, make u=v 4,
The 437th step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 438th step: establish cycle index g=1,
The 439th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 440th step: judge This moment j=3, g=1, m g=1, r j[m g]=∞, u=v 4, r j[u]=1.0,
Figure GSA00000073706000212
Therefore, the judgement formula is ∞>3.0, and correct judgment then continues to carry out,
The 441st step: order
Figure GSA00000073706000213
Be r 3[v 1]=0.0+2.0=2.0,
The 442nd step: order
Figure GSA00000073706000214
Be p 3[v 1]=v 4,
The 443rd step: make g=g+1, thus g=2,
The 444th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 445th step: judge
Figure GSA00000073706000215
This moment j=3, g=2, m g=3, r j[m g]=0.0, u=v 4, r j[u]=1.0,
Figure GSA00000073706000216
Therefore, the judgement formula is 0.0>3.2, and next adjacent node is judged in misjudgment
The 446th step: make g=g+1, thus g=3,
The 447th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 448th step: u is moved to S set, u=v by the Q set 4,
The 449th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, 2 nodes are arranged among the Q, Ψ=2 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2,
The 450th step: check r 3[v 1], r 3[v 2], be respectively 3.0,3.0, so l λ=1, make u=v 1,
The 451st step: establishing node u has Ω adjacent node, u=v 1, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=2, m 2=4,
The 452nd step: establish cycle index g=1,
The 453rd step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 454th step: judge This moment j=3, g=1, m g=2, r j[m g]=3.0, u=v 1, r j[u]=3.0,
Figure GSA00000073706000218
Therefore, the judgement formula is 3.0>4.0, and next adjacent node is judged in misjudgment
The 455th step: make g=g+1, thus g=2,
The 456th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 457th step: judge This moment j=3, g=2, m g=4, r j[m g]=1.0, u=v 1, r j[u]=3.0,
Figure GSA000000737060002110
Therefore, the judgement formula is 1.0>16.7, and next adjacent node is judged in misjudgment
The 458th step: make g=g+1, thus g=3,
The 459th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 460th step: u is moved to S set, u=v by the Q set 1,
The 461st step: judge whether current set Q is empty, at this moment Q={v 2, 1 node is arranged among the Q, Ψ=1 according to node subscript series arrangement from small to large, obtains l to the node among the current Q 1=2,
The 462nd step: check r 3[v 2], be respectively 3.0, so l λ=2, make u=v 2,
The 463rd step: establishing node u has Ω adjacent node, u=v 2, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 464th step: establish cycle index g=1,
The 465th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 466th step: judge
Figure GSA000000737060002111
This moment j=3, g=1, m g=1, r j[m g]=3.0, u=v 2, r j[u]=3.0,
Figure GSA00000073706000221
Therefore, the judgement formula is 3.0>4.1, and next adjacent node is judged in misjudgment
The 467th step: make g=g+1, thus g=2,
The 468th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 469th step: judge
Figure GSA00000073706000222
This moment j=3, g=2, m g=3, r j[m g]=0.0, u=v 2, r j[u]=3.0,
Figure GSA00000073706000223
Therefore, the judgement formula is 0.0>6.0, and next adjacent node is judged in misjudgment
The 470th step: make g=g+1, thus g=3,
The 471st step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 472nd step: u is moved to S set, u=v by the Q set 2,
The 473rd step: judge whether current set Q is empty, and this moment, Q=φ then jumped out circulation,
The 474th step: judge j≤τ, this moment j=4, τ=4, the judgement formula is 4≤4, correct judgment continues to carry out,
The 475th step: use single source impedance method of formation to calculate v 4Impedance r to all nodes 4[u k] and with respect to v 4Preposition node,
The 476th step: initialization p 4[] array, r 4[] array, S set, set Q,
The 477th step: p 4[v 1]=φ, p 4[v 2]=φ, p 4[v 3]=φ, p 4[v 4]=φ,
The 478th step: r 4[v 1]=∞, r 4[v 2]=∞, r 4[v 3]=∞, r 4[v 4]=0,
The 479th step: S=φ, Q={v 1, v 2, v 3, v 4,
The 480th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, v 4, 4 nodes are arranged among the Q, Ψ=4 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3, l 4=4,
The 481st step: check r 4[v 1], r 4[v 2], r 4[v 3], r 4[v 4], be respectively ∞, ∞, ∞, 0.0, so l λ=4, make u=v 4,
The 482nd step: establishing node u has Ω adjacent node, u=v 4, then u has 2 adjacent nodes, Ω=2, to the adjacent node of node u according to subscript series arrangement from small to large, and will order corresponding to m 1, m 2, obtain m 1=1, m 2=3,
The 483rd step: establish cycle index g=1,
The 484th step: judge g≤Ω, this moment g=1, Ω=2, correct judgment then continues to carry out,
The 485th step: judge
Figure GSA00000073706000224
This moment j=4, g=1, m g=1, r j[m g]=∞, u=v 4, r j[u]=0.0,
Figure GSA00000073706000225
Therefore, the judgement formula is ∞>2.0, and correct judgment then continues to carry out,
The 486th step: order Be r 4[v 1]=0.0+2.0=2.0,
The 487th step: order
Figure GSA00000073706000227
Be p 4[v 1]=v 4,
The 488th step: make g=g+1, thus g=2,
The 489th step: judge g≤Ω, this moment g=2, Ω=2, correct judgment then continues to carry out,
The 490th step: judge
Figure GSA00000073706000228
This moment j=4, g=2, m g=3, r j[m g]=∞, u=v 4, r j[u]=0.0, Therefore, the judgement formula is ∞>2.2, and correct judgment then continues to carry out,
The 491st step: order
Figure GSA000000737060002210
Be r 4[v 3]=0.0+2.2=2.2,
The 492nd step: order
Figure GSA000000737060002211
Be p 4[v 3]=v 4,
The 493rd step: make g=g+1, thus g=3,
The 494th step: judge g≤Ω, this moment g=3, Ω=2, circulation is then jumped out in misjudgment,
The 495th step: u is moved to S set, u=v by the Q set 4,
The 496th step: judge whether current set Q is empty, at this moment Q={v 1, v 2, v 3, 3 nodes are arranged among the Q, Ψ=3 according to node subscript series arrangement from small to large, obtain l to the node among the current Q 1=1, l 2=2, l 3=3,
The 497th step: check r 4[v 1], r 4[v 2], r 4[v 3], be respectively 2.0, ∞, 2.2, so l λ=1, make u=v 1,
The 558th step: finding source point is v 2Node v 3Preposition node p 2[v 3]=v 2, node v then 3With node v 2Between the unidirectional highway section of unidirectional highway section impedance minimum be e 3,
The 559th step: at the unidirectional highway section e of unidirectional highway section impedance minimum 3Last distribute traffic amount, order accumulative total one-way road section traffic volume amount
x e 3 = x e 3 + 50.00 = 150.00 ,
The 560th step: record is to the unidirectional highway section e of highway section impedance minimum 3On the volume of traffic constitute, the order accumulative total volume of traffic constitutes
T e 3 ( 2,4 ) = T e 3 ( 2,4 ) + 50.00 = 50.00 ,
The 561st step: make pointer s=s '=2,
The 562nd step: judge s=o n, this moment n=2, s=2, o n=2, correct judgment is then jumped out circulation,
The 563rd step: make n=n+1, thus n=3,
The 564th step: judge n≤Γ, this moment n=3, Γ=2, circulation is then jumped out in misjudgment,
The 565th step: judge i≤∏, this moment i=3, ∏=2, the judgement formula is 3≤2, misjudgment, then whole process finishes, and obtains one-way road section traffic volume and forms T.

Claims (1)

1. the assay method that constitutes of an one-way road section traffic volume is characterized in that:
Step 1: the initialization transport need on average is split as ∏ five equilibrium: OD with the OD matrix of expressing transport need 1Matrix, OD 2Matrix ..., OD Matrix, ∏ are the positive integer greater than 1, the OD matrix after each splits comprise Γ starting and terminal point to transport need, i.e. q i(o 1, d 1), q i(o 2, d 2) ..., q i(o n, d n) ..., q i(o Γ, d Γ), q wherein i(o n, d n) expression OD iStarting point is a node in the matrix
Figure FSA00000073705900011
, terminal point is node
Figure FSA00000073705900012
Transport need,
Step 2: with urban road network be set at resistive directed networks G (V, E), with the terminal in crossing in the urban road network and unidirectional highway section as described resistive directed networks G (V, node E) form set of node V, are designated as:
V={ υ 1, υ 2..., υ τ, τ is the number of node, υ 1, υ 2, υ τBe the node of directed networks, with the unidirectional highway section in the urban road network as described resistive directed networks G (the limit collection is designated as for V, directed edge E):
E={e 1, e 2..., e σ..., e η, e 1, e 2..., e σ..., e ηThe unidirectional highway section of the directed edge representative of expression directed networks, σ is the subscript in a certain unidirectional highway section, and η is the number in unidirectional highway section in the directed networks, and (V E) carries out initialization, obtains initialization directed networks G to directed networks G 0(V, E),
Step 3: initialization impedance, the volume of traffic and the volume of traffic constitute, and make directed networks G 0(V, E) in the accumulative total volume of traffic on all unidirectional highway sections
Figure FSA00000073705900013
σ=1,2 ..., η,
Figure FSA00000073705900014
Represent unidirectional highway section e σOn the accumulative total volume of traffic; Make the volume of traffic on all unidirectional highway sections constitute
Figure FSA00000073705900015
σ=1,2 ..., η,
Figure FSA00000073705900016
Represent unidirectional highway section e σGo up by node υ oTo υ dThe volume of traffic, the impedance in unidirectional highway section is:
w e σ = w e σ 0 ( 1 + α ( x e σ c e σ ) β )
Wherein
Figure FSA00000073705900018
For when on the unidirectional highway section during without any vehicle, the time cost that vehicle passes through, it is that the geometrical length in known unidirectional highway section is divided by the known design speed of a motor vehicle; Be the one-way road section traffic volume amount, Be the maximum traffic capacity in unidirectional highway section, Be unidirectional highway section saturation degree; α and β are empirical parameter, α=0.15, and β=4.0,
Step 4: loop initialization pointer i=1,
Step 5: judge i≤∏, if execution in step 6 then; Otherwise, finishing, the volume of traffic that obtains on every highway section constitutes σ=1,2 ..., η, 1≤o≤τ, 1≤d≤τ,
Step 6: calculate current impedance matrix and be source point and with any node respectively with the impedance matrix generating algorithm, search shortest path for traffic in distributing and prepare with respect to the preposition node of another node arbitrarily,
Step 7: use the traffic allocation algorithm with transport need OD iMatrix allocation arrives last time directed networks G I-1(V, E) on, form current directed networks G i(V E), and records to network G i(V, E) volume of traffic on every the unidirectional highway section on constitutes,
Step 8: make i=i+1, return step 5,
Described impedance matrix generating algorithm is:
Step 6.1: loop initialization pointer j=1,
Step 6.2: if j≤τ, then execution in step 6.3, otherwise impedance matrix generates and to finish, and obtains the impedance between each node, and to obtain respectively with any node be source point and with respect to the preposition node of another node arbitrarily,
Step 6.3: adopt single source impedance generating algorithm to calculate arbitrary node υ jArrive all node υ as source point kImpedance r jk], k=1,2 ... τ, and obtain preposition node with respect to arbitrary node of source point,
Step 6.4: make j=j+1, return step 6.2,
Step 6.3.1: initialization p j[] array, r j[] array, S set, set Q, p jk] represent with node υ jNode υ for source point kPreposition node, k=1,2 ... τ, the order:
p jk]=φ
r jk] represent with source point υ jBe starting point, node υ kBe the path impedance of terminal point, k=1,2 ... τ, order:
r j [ υ k ] = 0 k = j ∞ k ≠ j
S set is used to deposit the node of having handled, initialization S, order
Figure FSA00000073705900023
Set Q is used to deposit untreated node, and initialization Q makes Q=V,
Step 6.3.2: judge that whether the Q set is empty, if be empty, then obtains node υ jArrive all node υ as source point kImpedance r jk], k=1,2 ..., τ, and obtain preposition node with respect to arbitrary node of source point; Otherwise, execution in step 6.3.3,
Step 6.3.3: establishing among the current set Q has ψ node, and ψ≤τ according to node subscript series arrangement from small to large, and corresponds respectively to l with tactic node subscript to the node among the current Q 1, l 2..., l ψ, obtain
Figure FSA00000073705900024
1≤l ζ≤ τ,
Step 6.3.4: check
Figure FSA00000073705900025
Therefrom find out minimum
Figure FSA00000073705900026
And with node
Figure FSA00000073705900027
As node u,
Step 6.3.5: determine source point υ respectively jTo the impedance of each adjacent node of node u, and definite source point υ jThe preposition node of adjacent node of node u,
Step 6.3.6: node u is moved to the S set by the Q set, execution in step 6.3.2,
Wherein, the described definite respectively source point υ of step 6.3.5 jTo the impedance of each adjacent node of node u, and definite source point υ jThe method of preposition node of adjacent node of node u be:
Step 6.3.5.1: establishing node u has Ω adjacent node, and Ω≤τ according to subscript series arrangement from small to large, and corresponds respectively to m with tactic node subscript to the adjacent node of node u 1, m 2..., m Ω, obtain the adjacent node of tactic node u
Figure FSA00000073705900031
1≤m ρ≤ τ,
Step 6.3.5.2: establish cycle index g=1,
Step 6.3.5.3: if g≤Ω, execution in step 6.3.5.4 then, otherwise, current source point υ obtained jTo the impedance of each adjacent node of node u, and source point υ jThe preposition node of adjacent node of node u,
Step 6.3.5.4: if
Figure FSA00000073705900032
Execution in step 6.3.5.5 then; Otherwise, skip to step
6.3.5.7, wherein, Expression by node u to adjacent node
Figure FSA00000073705900034
The impedance of unidirectional highway section, if u is extremely
Figure FSA00000073705900035
Between many unidirectional highway sections are arranged, then
Figure FSA00000073705900036
Refer to the impedance on that unidirectional highway section of impedance minimum,
Step 6.3.5.5: make source point υ jTo adjacent node
Figure FSA00000073705900037
Impedance equal source point υ jTo the impedance of u add u with
Figure FSA00000073705900038
Impedance on the unidirectional highway section, even
Figure FSA00000073705900039
Step 6.3.5.6: order
Figure FSA000000737059000310
With respect to υ jFor the preposition node of starting point is υ, even
Figure FSA000000737059000311
Step 6.3.5.7: make g=g+1, return step 6.3.5.3,
Step 7 is described " uses the traffic allocation algorithm with transport need OD iMatrix allocation to directed networks G last time (V, E) on, form current directed networks G (V, E), and directed networks G before the record (V, E) volume of traffic on every the unidirectional highway section on constitutes " method be:
Step 7.1: establish transport need OD iComprise in the matrix Γ starting and terminal point to transport need, promptly
q i(o 1, d 1), q i(o 2, d 2) ..., q i(o n, d n) ..., q i(o Γ, d Γ), q wherein i(o n, d n) expression OD iStarting point is a node in the matrix
Figure FSA000000737059000312
Terminal point is a node
Figure FSA000000737059000313
Transport need,
Step 7.2: establish cycle index n=1,
Step 7.3: if n≤Γ, then execution in step 7.4; Otherwise, obtain current one-way road section traffic volume amount and road section traffic volume and form,
Step 7.4: with starting point-terminal point to allocation algorithm distribute the starting point that obtains by step 1
Figure FSA000000737059000314
Terminal point
Figure FSA000000737059000315
Between transport need q i(o n, d n),
Step 7.5: make n=n+1, return step 7.3,
Wherein, the starting point-terminal point described in the step 7.4 to allocation algorithm be:
Step 7.4.1: establish pointer s=d n, d nBe the node subscript of the terminal point of point-to-point transmission transport need,
Step 7.4.2: if pointer s=o n, then obtain transport need q i(o n, d n) in starting point be Terminal point is
Figure FSA000000737059000317
Shortest path on allocation result, otherwise, execution in step 7.4.3,
Step 7.4.3: find source point to be
Figure FSA00000073705900041
Node υ sPreposition node Node υ then sWith source point be
Figure FSA00000073705900043
Node υ sPreposition node υ ' sBetween the unidirectional highway section of unidirectional highway section impedance minimum be e a,
Step 7.4.4: at the unidirectional highway section e of unidirectional highway section impedance minimum aLast distribute traffic amount, order accumulative total one-way road section traffic volume amount x e a = x e a + q i ( o n , d n ) ,
Step 7.4.5: record is to the unidirectional highway section e of highway section impedance minimum aOn the volume of traffic constitute, the order accumulative total volume of traffic constitutes T e a ( o n , d n ) = T e a ( o n , d n ) + q i ( o n , d n ) ,
Step 7.4.6: make pointer s=s ', return step 7.4.2 then.
CN2010101399398A 2010-04-02 2010-04-02 Testing method of one-way road section traffic volume composition Expired - Fee Related CN101930666B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101399398A CN101930666B (en) 2010-04-02 2010-04-02 Testing method of one-way road section traffic volume composition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101399398A CN101930666B (en) 2010-04-02 2010-04-02 Testing method of one-way road section traffic volume composition

Publications (2)

Publication Number Publication Date
CN101930666A true CN101930666A (en) 2010-12-29
CN101930666B CN101930666B (en) 2012-04-11

Family

ID=43369821

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101399398A Expired - Fee Related CN101930666B (en) 2010-04-02 2010-04-02 Testing method of one-way road section traffic volume composition

Country Status (1)

Country Link
CN (1) CN101930666B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005182383A (en) * 2003-12-18 2005-07-07 Ryuichi Kitamura Od data estimation method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005182383A (en) * 2003-12-18 2005-07-07 Ryuichi Kitamura Od data estimation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《清华大学学报(自然科学版)》 20001231 段进宇 《由路段交通流量反估出行od矩阵技术的应用》 123-126页 1 , 2 *
《第五届中国智能交通年会暨第六届国际节能与新能源汽车创新发展论坛优秀论文集(上册)--智能交通》 20091231 刘斌; 杨晓光 《OD估计技术在交通预报中的应用研究》 全文 1 , 2 *

Also Published As

Publication number Publication date
CN101930666B (en) 2012-04-11

Similar Documents

Publication Publication Date Title
CN101930666B (en) Testing method of one-way road section traffic volume composition
CN107190622B (en) Permanent deformation for asphalt pavement method for quick predicting based on full temperature range Temperature Distribution
CN102012356B (en) Quick test method for aggregate grade of asphalt concrete pavement
CN104499396A (en) Determination of optimal blending ratio of regenerative old materials of bituminous pavement and method thereof
CN1082060C (en) Atactic copolymers of propylene with ethylene
CN102863182A (en) Ultrathin rubber asphalt wearing layer mixing material for preventive conservation
CN109166601A (en) It is a kind of to judge phylogenetic method between species using three concepts
CN204492085U (en) A kind of ground heating earth construction
CN102558215A (en) Cracking process of organosilicon high-boiling components
Ravasi et al. San Giovanni in Laterano 2 Project (SGL2)
Bai et al. Decision analysis of slope ecological restoration based on AHP
CN206092406U (en) 20 40 two screw air compressor's of horsepower rotor combination
CN102682596B (en) Traffic information predicating method based on chaos theory
CN208502160U (en) A kind of ground ruggedized construction
Ding et al. Multi-objective disassembly line balancing via modified ant colony optimization algorithm
Jain A Decision Making in Selection of Bricks Using Multiple Attribute Decision Making Methods
CN204921366U (en) 7. 5 -10 horsepower twin -screw air compressor's rotor combination
CN109811982A (en) The aqueous half prefabricated high-environmental athletic floor of one kind and its construction method
CN206092415U (en) 100 150 two screw air compressor's of horsepower rotor combination
Crumley et al. Improving snow modeling by assimilating observational data collected by citizen scientists
ICHITANI et al. Spatial Variation of Slip Deficit Rate at the Nankai Trough, Southwest Japan Inferred from Three-Dimensional GPS Crustal Velocity Fields—Repeated Geodetic Inversion Analyses for the Shifted Target Area—
Ciuna et al. Parametric measurement of partial damage in building
RU2007126457A (en) METHOD FOR BUILDING FASTENING IN WINTER
CN206111551U (en) 125 175 two screw air compressor's of horsepower rotor combination
Imran et al. Development of Paving Blocks to Reduce Environmental Heat by Using Mixed Materials

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120411

Termination date: 20160402

CF01 Termination of patent right due to non-payment of annual fee