CN104183119B - Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD - Google Patents

Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD Download PDF

Info

Publication number
CN104183119B
CN104183119B CN201410410008.5A CN201410410008A CN104183119B CN 104183119 B CN104183119 B CN 104183119B CN 201410410008 A CN201410410008 A CN 201410410008A CN 104183119 B CN104183119 B CN 104183119B
Authority
CN
China
Prior art keywords
traffic
od
time
section
volume
Prior art date
Application number
CN201410410008.5A
Other languages
Chinese (zh)
Other versions
CN104183119A (en
Inventor
周素红
柳林
杨靖芸
郝新华
陈建国
肖露子
Original Assignee
中山大学
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 中山大学 filed Critical 中山大学
Priority to CN201410410008.5A priority Critical patent/CN104183119B/en
Publication of CN104183119A publication Critical patent/CN104183119A/en
Application granted granted Critical
Publication of CN104183119B publication Critical patent/CN104183119B/en

Links

Abstract

The invention discloses a kind of based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, for predicting that real-time section is gone on a journey, attracted the volume of traffic, and carry out arithmetic for real-time traffic flow distribution.It is an advantage of the current invention that: the division of current scattered traffic zone, the distribution of traffic zone traffic volume forecast, OD Prediction of Traffic Flow in Real-time, link flow have been carried out the integration of system, build traffic zone based on road network division methods and the prognoses system of real-time traffic flow distribution, the knowledge and technology of trip origin and destination and Traffic growth rate is carried out organic combination, has been realized from software respective.

Description

Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD

Technical field

The present invention relates to technical field of intelligent traffic, especially relate to generation and the traffic attraction of a kind of road grid traffic Real-time predicting method.

Background technology

Build and extensively applying and real time information system in management in smart city along with intelligent transportation system Growing, traffic administration gradually forms intelligent, Statistical, the feature of real-time.Vehicle supervision department Obtain substantial amounts of real time traffic data with research institution, and be applied to the field of traffic control such as real-time road. Substantial amounts of traffic information data builds historical data base, provides important data for traffic administration person and research worker Support.In order to carry out effective traffic administration and provide Traffic Information accurately for traveler, need base Trip generation forecast and traffic attraction prediction accurately is carried out in Real-time Traffic Information.Carry out in short-term based on real time traffic data Traffic trip and traffic attraction prediction contribute to trip distribution modeling, this is by the important of traffic administration and planning Part.

Passing through traffic study, according to personnel's structure of each traffic zone traffic trip and the prediction of traffic attraction at present more The factors such as one-tenth, motor vehicles retention calculate travel amount and the traffic attraction of traffic zone.First, the friendship of large sample amount Logical investigation needs substantial amounts of time and fund input;Secondly, divide the scope of traffic zone and can make according to difference Become the inaccurate of prediction, and the prediction of traffic zone cannot carry out Proper Match with traffic network, it was predicted that result Disagree with practical situation;3rd, based on static investigation and prediction be difficult to reflect real-time traffic stream mode and Carry out accurate short-time traffic flow forecast.

Summary of the invention

In order to overcome above-mentioned the deficiencies in the prior art, the invention provides a kind of based on section OD counter push away real-time Traffic flow distribution forecasting method, is traffic trip and the suction dividing the real-time OD distribution carried out based on traffic zone The Forecasting Methodology of the amount of drawing, for improving the accuracy of trip forecast of distribution.

To achieve these goals, the technical solution adopted in the present invention is:

The present invention be to " traffic zone based on section divide, the volume of traffic based on road section traffic volume community the most pre- Survey, the distribution of OD Prediction of Traffic Flow in Real-time, link flow " carry out the integration of system, build based on section Traffic zone division methods and the Forecasting Methodology of real-time traffic flow distribution, it is achieved to travelling OD between hand in real time The prediction of flux, and further according to shortest path distribution method, by OD between real-time estimate traffic Amount is assigned in each section of road network, obtains real-time estimate traffic flow distribution situation.

A kind of based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, comprise the following steps:

1.1, every section in transportation network is considered as traffic zone, to transport information (Floating Car GPS gathered Data) carry out OD extraction and section coupling;

1.2, according to prediction time, this is pre-to choose a few sky in the Floating Car OD traffic data storehouse of history The OD data of each half an hour before and after the survey moment, and the calculating that is averaged, based on result of calculation OD matrix, substitutes in OD matrix;

1.3, according to prediction time, the OD traffic data in the first two moment prediction time on the same day is chosen, right Its Time Series Method once moved, obtains prediction time each traffic zone (section) O, D Traffic volume forecast value, and as in the aggregation of forecasts value substitution OD matrix of each O, D;

1.4, push away travelling OD matrix is counter according to balanced growth Y-factor method Y, obtain OD under predicted time The volume of traffic between to;

1.5, according to the algorithm of shortest path, the prediction OD volume of traffic is carried out traffic flow distribution, obtains in real time Road traffic delay prediction.

Compared with prior art, the invention has the beneficial effects as follows: section is considered as traffic zone, mesh can be solved The problem that go on a journey in front traffic zone, the attraction volume of traffic is bad with road network, meanwhile, can reduce traffic zone The cost of investigation, obtains the volume of traffic gone on a journey between traffic zone, attract by real time traffic data, solves real Time traffic assignation forecasting problem.

Accompanying drawing explanation

Fig. 1 is the flow chart of the open method of the present invention.

Fig. 2 is to use the traffic zone division methods of the present invention to interrupt the road network of Panyu District of Guangzhou City to carry with OD Take schematic diagram.

Fig. 3 is that 0-1 distributes method flow chart.

Fig. 4 is that increment distributes method flow chart schematic diagram.

Fig. 5 is straight average method flow chart.

Detailed description of the invention

The present invention is further described below in conjunction with the accompanying drawings.

The techniqueflow of the present invention as shown in Figure 1, including division, the data base of traffic zone based on section Structure, real-time OD traffic volume forecast, arithmetic for real-time traffic flow distribution four steps.

1.1, traffic zone based on section divides:

The unit of traditional traffic four-phase model is traffic zone, the present invention section is replaced traffic zone as The unit of four-phase model, the section to whole region, according to the principle having crossing just to interrupt by the road in region Net interrupts the section into more junior unit.Traffic zone is instead based on section, anti-with the OD that this serves below Push away and traffic assignation.Specific algorithm flow process is as follows:

Choose the basic road network of research, interrupt according to crossing and section road axis angle of turning is less than The principle that the corner of 150 ° is interrupted, is broken into the section of less yardstick by Regional Road Network, the section after interrupting Traffic zone will be replaced, it is numbered and adds up the minimum single of OD travel amount as four-phase model Unit.Traditional OD matrix is all with traffic analysis community as analytic unit, so hands over carrying out Path selection Need the traffic flow of each traffic zone to match on the road network of correspondence when of through-flow distribution, this process meeting Introduce error.The most directly use section as the analytic unit of OD matrix, can avoid at this Systematic error is introduced during one.

1.2, the structure of data base:

The structure of data base be based on gps data is carried out real-time OD extraction, section coupling and based on section OD collect statistics three parts.

1. O, D extract in real time: be first ranked up gps data according to car number and date, time, Then O point or D point are judged according to the change of passenger carrying status field in gps data, such as: carrying shape In state, 1 is load-carrying, and 0 is unloaded, then when the passenger carrying status of same car is become 0 from 1, for D;When same The passenger carrying status of one car from 0 become 1 time, for O;When finding O point, the most down find corresponding D Point;When finding D point, then corresponding O point is found in backtracking, mates OD point with this right.

2. section coupling: O, D of extracting point is converted to shapefile dot format, uses nearest section Method coupling section, i.e. search out the section that GPS point is nearest, and section numbering be assigned to GPS point.

3. OD based on section collects statistics: to matching the GPS OD point in section, segment number of satisfying the need is carried out Collect statistics, passenger carrying status is carried out counting statistics, obtain O based on section, D matrix.

1.3, real-time OD traffic volume forecast:

In real time OD traffic volume forecast is divided into three steps: first, determines basis OD matrix according to historical data;The Two, determine OD Matrix prediction value according to real time data;3rd, OD are counter estimates the OD friendship calculating prediction time Flux.

1. basis OD matrix is determined according to historical data

Historical data based on predicted time point and the proportionate relationship of history Road network traffic, be predicted.As wanted In prediction morning 8:00 road network each OD between the volume of traffic, by extracting each OD of road network to 8: The historical traffic amount of 00, and it is averaged value as historical data, bring in OD matrix.

2. OD Matrix prediction value is determined according to real time data

Forecasting Methodology based on lineary system theory specifically includes that linear statistical recurrence, historical trend method, time Sequence prediction method, Kalman prediction method, Smoothing Forecasting and adaptive weighting forecast model Method etc..The present embodiment carries out short-term traffic flow as a example by Time Series Forecasting Methods Single moving average method Prediction.

If x1, x2, x3..., xnFor specific O point, the life with t as time interval, in every period Become the volume of traffic, or specific D point, the attraction volume of traffic with t as time interval, in every period.

Therefore, it is a time series, and sample size is n.

The formula of moving average method is as follows:

M T ( 1 ) = x T + x T - 1 + ... + x T - N + 1 N , T ≥ N

Wherein, N is the data amount check that each rolling average comprises;It it is the Single moving average of T phase Value.

The effect of rolling average is smoothing data, eliminates some random disturbances, makes long-term trend reveal, Thus can be used for trend analysis and prediction.

Generally, if time series does not has obvious mechanical periodicity and Long-term change trend, the available T phase Single moving average value is as the predictive value of T+1 phase, it may be assumed that

x ^ T + 1 = M T ( 1 ) = x T + x T - 1 + ... + x T - N + 1 N

Abbreviation obtains recurrence formula:

x ^ T + 1 = M T ( 1 ) = M T - 1 ( 1 ) + x T - x T - N N = x ^ T + x T - x T - N N

When n is large, amount of calculation can be greatly reduced with recurrence formula.

Therefore, according to OD between historical data, the volume of traffic in following time interval is predicted.

Generation traffic volume forecast value and the attraction of each D obtaining each O in OD matrix may finally be added up Traffic volume forecast value, by historical data base, it was predicted that around the moment, the traffic flow data in the range of certain time is pressed It is assigned in OD matrix according to origin and destination, as the basis matrix of OD matrix.Meanwhile, according to time series side Method and the traffic flow data in several moment before prediction time, extrapolate the origin and destination traffic volume forecast of prediction time Value, as the predictive value in OD matrix.

3. OD is counter estimates the real-time OD volume of traffic

This research uses the distribution of growth factor methed prognosis traffic volume.

The algorithm steps of growth factor methed is as follows:

The first step, makes calculation times m=0

Second step, in given present situation OD tableU in Tm and in the future OD tablei、Vj

WhereinRepresent the m time when calculating, the volume of traffic between traffic zone i to traffic zone j;

Representing the m time when calculating, there is the volume of traffic in i-th traffic zone;

Represent the m time when calculating, the attraction volume of traffic of jth traffic zone;

UiRepresent the generation volume of traffic in OD table in future;

VjRepresent the attraction volume of traffic in OD table in future;

3rd step, obtains occurring and the rate of increase attracting the volume of traffic of each community

F O i m = U i / O i m

F D j m = V j / D j m

4th step, seeks the approximation of the m+1 time distribution volume of traffic

q i j m + 1 = q i j m × f ( F O i m , F D j m )

5th step, convergence differentiates, it determines formula is:

O i m + 1 = Σ j q i j m + 1

D j m + 1 = Σ j q i j m + 1

1 - &epsiv; < F O i m + 1 = U i / O i m + 1 < 1 + &epsiv;

1 - &epsiv; < F D j m + 1 = V j / D j m + 1 < 1 + &epsiv;

In formula, UiRepresent the generation volume of traffic in OD table in future;

VjRepresent the attraction volume of traffic in OD table in future;

Represent that the m time calculating of i community occurs growth factor;

Represent that the m time calculating of j community attracts growth factor;

It is the iterative manner of growth factor, relevant with every kind of method;

ε represents any given error constant, typically takes 3%.

According to functionKind different, growth factor methed can be divided into normal growth factor methed (Unique Growth Factor Method), balanced growth Y-factor method Y (Average Growth Factor Method), the end special Law (Detroit Method), Follett, Mary Parker method (Fratar Method) and Buddhist Nice method (Furness Method).

The present invention uses balanced growth Y-factor method Y, and the method formula is simple and clear in calculating, and shortcoming is that convergence is slow, repeatedly In generation, computational accuracy was low often.Balanced growth coefficient is expressed as follows:

F O i 0 = U i / O i 0

F D j 0 = V j / D j 0

First approximation:

(4) arithmetic for real-time traffic flow distribution:

Generally traffic flow distribution method is divided into equilibrium assignmen and the big class of non-equilibrium distribution two in the world.For completely Meet the poised state of Wardrop principle definition, referred to as EqP method;For use heuristic or its The distribution model of his approximation method, the most non-equilibrium distribution method.

Obtained after real-time OD matrix by OD estimation model, then be respectively adopted 0-1 distribution method, increment divides Join method, straight average method etc. three kinds non-equilibrium allocation algorithm development behavior traffic flow distribution model, by new real-time OD Traffic growth rate, to road network, obtains the arithmetic for real-time traffic flow distribution situation of road network.

1. 0-1 distributes method

0-1 distribution method also referred to as has entirely completely without distribution method, is the most most basic Traffic growth rate method, in the U.S. In the traffic programme of city, Chicago, obtain application first.This algorithm has 2 features: one is to think the out of shape of section Time is the most relevant with the length in section, is a definite value, unrelated with the size of the volume of traffic on section, Bu Sui road The degree of crowding on road and be varied from;Two is to think that all drivers of OD of same group select identical Route.The main processes of calculation of this distribution method is to find shortest path.

0-1 distribution method is simple and efficiency is high.It is applicable in the rare remote districts of road, and general city road network Traffic growth rate in be not suitable for using this distribution method.Calculation procedure can be summarized as follows:

The first step, calculate each OD in network between shortest path;

Second step, is fully allocated to the volume of traffic between OD on corresponding shortest path.

Owing to 0-1 distribution method can not reflect crowded effect, it is mainly used for some uncongested road network, this distribution method For there is no the situation of network that the traffic capacity limits.Therefore, it is suggested that range is: road between intercity The unrestricted area of the traffic capacity can use;The traffic flow distribution of general urban road network should not use the party Method.In practice due to its simple and practical characteristic, generally as the basis of other various distribution techniques, increasing Reusability in the amount method such as distribution method and EqP method.

2. increment distribution method

Increment distribution method is the EqP method of a kind of approximation.It is that the OD volume of traffic is divided equally into some equal portions, Cyclically by the OD Traffic growth rate of each equal portions to network.The OD of cycle assignment one equal portions hands over each time Flux, on corresponding shortest path, often distributes the travel time the most just recalculating and updating each section, so The shortest path between each OD of network is recalculated afterwards by the travel time after updating.Subsequent cycle is pressed after updating Shortest path distribute the OD volume of traffic of next equal portions.Its calculation procedure is as follows:

The first step, initializes.It is divided equally into N equal portions, order by often organizing the OD volume of traffic

Second step, updates the travel time in section. Calculate shortest path between OD.

3rd step, increment distributes.Gained is calculated by second stepBy 0-1 distribution method by the OD traffic of 1/N AmountIt is assigned in network.So obtain one group of additional traffic amount

4th step, the volume of traffic adds up.Order

5th step, it is determined that.If n=N, stop calculating.Current link counting is i.e. last solution;As Really n < N, makes n=n+1, returns second step.

Increment distribution method obtains complexity and solves accuracy all between 0-1 distribution method and EqP method. The result distributing method as N=1 with 0-1 is consistent;As N → ∞, its solution is consistent with the solution of EqP method. The method is simple, is all provided with using in many ripe business softwares, and degree of accuracy can be according to N's Size regulates, and is often applied in actual road network Traffic growth rate.But this distribution method is still It is a kind of approximation method, sometimes by section the least for too much Traffic growth rate to some capacity.

3. straight average method

Straight average method is a kind of round-robin technique between increment distribution method and EqP method.It is Constantly adjustment has been assigned to the volume of traffic on each section and has moved closer to or arrive equilibrium assignmen.Often walking circulation In, carry out a 0-1 distribution according to the volume of traffic being assigned on each section and obtain the additional of one group of each section The volume of traffic.Then the additional traffic amount obtained in circulating with this with the distribution volume of traffic in section each in last circulation It is weighted the section distribution volume of traffic that is average and that obtain in this circulation.When the distribution in adjacent two circulations is handed over Flux equal to definite value time, can stop calculate, this equivalence is a value the least, approximates 1, i.e. it is believed that The distribution volume of traffic in adjacent two circulations is sufficiently close to.The distribution volume of traffic obtained in last circulation It it is i.e. the final volume of traffic.Its main algorithm step is as follows:

The first step, initializes.Carry out a 0-1 distribution according to the free travel time in each section, obtain each road The distribution volume of traffic of section

Second step, according to the distribution volume of traffic in current each sectionCalculate each Link Travel Time, and calculate OD Between shortest path.

3rd step, the Link Travel Time calculated according to second step and the OD volume of traffic carry out 0-1 distribution, Obtain additional traffic amount F in each sectiona

4th step, calculates the current amount of traffic in each section by the method for weighted mean methodAs follows:

x a n = ( 1 - &theta; ) x a n - 1 + &theta;F a , 0 &le; &theta; &le; 1

5th step, ifThen loop termination, otherwise re-starts 0-1 distribution, Perform circulation.

In the 4th step and the 5th step, weight coefficient θ and parameter K need to be defined by reckoner oneself.Parameter Number of times and the operation efficiency of circulation that can determine circulation is set.

Straight average method is the most simply to be suitable for, again closest to a kind of distribution method of EqP method;Often step is followed In ring, the value of weight coefficient θ is in strict accordance with mathematical programming model value, can be balanced the solution of distribution.

By above method, the service level of road network can be obtained with Real-time and Dynamic, for special traffic events (as interior Flood, vehicle accident etc.) traffic flow is dredged under situation offer decision support.

Above-mentioned based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, combine traffic zone and draw Point, short-term traffic flow prediction, OD Prediction of Traffic Flow in Real-time, the aspect such as link flow distribution, it is achieved that to going out Row OD between the prediction of real-time traffic amount, and further according to the method for shortest path distribution, by OD pair Between real-time estimate Traffic growth rate in each section of road network, obtain real-time estimate traffic flow distribution feelings Condition.

Compared with conventional method, the invention has the beneficial effects as follows: one, section is considered as traffic zone by the present invention, The problem that go on a journey in current traffic zone, the attraction volume of traffic is bad with road network can be solved;Two, the present invention can With the cost of reduction traffic zone investigation, obtain going on a journey between traffic zone by real time traffic data, attract The volume of traffic, solves real-time traffic distribution forecast problem.

Based on above feature, the present invention announce based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD Can play a great role at aspects such as traffic flow real-time estimate, road condition assessment, also can be digital city and wisdom The development in city increases power-assisted.

Embodiment described above only have expressed the embodiment that the present invention is possible, and it describes more concrete and detailed, But therefore can not be interpreted as the restriction to the scope of the claims of the present invention.It should be pointed out that, for this area For those of ordinary skill, without departing from the inventive concept of the premise, it is also possible to make some deformation and improvement, These broadly fall into protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be with claims It is as the criterion.

Claims (16)

1. one kind based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, it is characterised in that include Following steps:
S1. every section in transportation network is considered as traffic zone, i.e. O and D in OD matrix;
S2. according to the GPS data from taxi passenger carrying status field obtained, it is judged that OD pair of trip, go out Row OD point extracts, and with traffic zone, the OD point of extraction is carried out section mates, and then collects and obtain each Traffic zone OD between the volume of traffic;
S3. the traffic zone OD obtained according to step S2 between the volume of traffic collect, build history OD and hand over Flux data storehouse;
S4. according to prediction time, choose in history OD traffic data storehouse around this prediction time of a few sky A range of OD data, and the calculating that is averaged, by OD matrix based on result of calculation, Substitute in OD matrix;
S5. according to prediction time, the OD traffic data in several moment before choosing prediction time on the same day, according to It is moved by Time Series Method, obtains the traffic volume forecast value of prediction time each traffic zone O, D, And the aggregation of forecasts value as each O, D substitutes in OD matrix;
S6. according to step S4 and S5 calculated basis OD matrix and O, D aggregation of forecasts value, use Balanced growth Y-factor method Y pushes away OD matrix is counter, obtain each OD of prediction time between the volume of traffic Distribution;
S7. use shortest path distribution algorithm, by each OD between real-time estimate Traffic growth rate arrive On each section of road network, obtain the real-time traffic amount distribution forecast in each section in this prediction time road network.
The most according to claim 1 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its Being characterised by, step S1 is to be interrupted less than 150 ° according to crossing and section angle in each section in road network, And each section after interrupting is as traffic zone.
The most according to claim 2 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its Being characterised by, step S2 judges that the method for OD pair of trip is: to gps data according to car number and day Phase, time are ranked up, then according to the change of passenger carrying status field in gps data judge O point or D point, when the passenger carrying status of same car is become zero load from load-carrying, for D;Passenger carrying status when same car From zero load become load-carrying time, for O;When finding O point, the most down find corresponding D point;When finding D During point, then corresponding O point is found in backtracking, mates OD point with this right.
The most according to claim 3 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its It is characterised by, the OD point of extraction is carried out section with traffic zone by step S2 and mates, collect and obtain each Traffic zone OD between the detailed process of the volume of traffic be: O, D of extracting point is converted to Shapefile dot format, uses the method coupling section in nearest section, i.e. searches out the section that GPS point is nearest, And section numbering is assigned to GPS point: to matching the GPS OD point in section, segment number of satisfying the need is converged Total statistics, carries out counting statistics to passenger carrying status, obtains O based on section, D matrix.
The most according to claim 4 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its It is characterised by, according to prediction time in step S5, chooses the OD traffic in the first two moment prediction time on the same day Amount data, the Time Series Method that it is once moved, obtain prediction time each traffic zone O, D Traffic volume forecast value.
The most according to claim 5 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its It is characterised by, the Time Series Method that it is once moved, specifically uses Single moving average method, tool Body is:
If x1, x2, x3..., xnFor specific O point, the life with t as time interval, in every period Become the volume of traffic, or x1, x2, x3..., xnFor specific D point, with t as time interval, every period The interior attraction volume of traffic;
Therefore, it is a time series, and sample size is n;
The formula of moving average method is as follows:
M T ( 1 ) = x T + x T - 1 + ... + x T - N + 1 N , T &GreaterEqual; N
Wherein, N is the data amount check that each rolling average comprises;It it is the Single moving average of T phase Value;
When time series does not has obvious mechanical periodicity and Long-term change trend, make by the Single moving average value of T phase It it is the predictive value of T+1 phaseThat is:
x ^ T + 1 = M T ( 1 ) = x T + x T - 1 + ... + x T - N + 1 N
Abbreviation obtains recurrence formula:
x ^ T + 1 = M T ( 1 ) = M T - 1 ( 1 ) + x T - x T - N N = x ^ T + x T - x T - N N .
The most according to claim 6 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its It is characterised by, step S6 uses the balanced growth Y-factor method Y detailed process that push away counter to OD matrix be:
The first step, makes calculation times m=0;
Second step, in given present situation OD tableAnd the in the future U in OD tablei、Vj
WhereinRepresent the m time when calculating, the volume of traffic between traffic zone i to traffic zone j;
Representing the m time when calculating, there is the volume of traffic in i-th traffic zone;
Represent the m time when calculating, the attraction volume of traffic of jth traffic zone;
UiRepresent the generation volume of traffic in OD table in future;
VjRepresent the attraction volume of traffic in OD table in future;
3rd step, obtains the rate of increase that the volume of traffic occurs of each communityWith the rate of increase attracting the volume of traffic
F O i m = U i / O i m
F D j m = V j / D j m
4th step, seeks the approximation of the m+1 time distribution volume of traffic
q i j m + 1 = q i j m &times; f ( F O i m , F D j m )
5th step, convergence differentiates, it determines formula is:
O i m + 1 = &Sigma; j q i j m + 1
D j m + 1 = &Sigma; j q i j m + 1
1 - &epsiv; < F O i m + 1 = U i / O i m + 1 < 1 + &epsiv;
1 - &epsiv; < F D j m + 1 = V j / D j m + 1 < 1 + &epsiv;
In formula,
Represent that the m time calculating of i-th community occurs growth factor;
Represent that the m time calculating of jth community attracts growth factor;
It is the iterative manner of growth factor, relevant with every kind of method;
ε represents any given error constant;
According to functionKind different, growth factor methed is divided into normal growth factor methed, balanced growth Y-factor method Y, Detroit method, Follett, Mary Parker method or Buddhist Nice method.
The most according to claim 7 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its Being characterised by, described growth factor methed uses balanced growth Y-factor method Y, its balanced growth coefficientTable Show as follows:
F O i 0 = U i / O i 0
F D j 0 = V j / D j 0
First approximation:
The most according to claim 8 based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD, its Being characterised by, step S7 is to use 0-1 distribution method, increment distribution method, straight average method three kinds non-equilibrium point Join algorithm development dynamic traffic assignment model, by new real-time OD Traffic growth rate to road network, its concrete mistake Cheng Wei:
Wherein 0-1 distributes method, be calculate each OD in network between shortest path, by the friendship between OD Flux is fully allocated on corresponding shortest path;
Increment distribution method, is that the OD volume of traffic is divided equally into some equal portions, is cyclically handed over by the OD of each equal portions Flux is assigned in network;Every time the OD volume of traffic of cycle assignment one equal portions is on corresponding shortest path, often The travel time in each section is the most just recalculated and updates in distribution, then the travel time after updating again based on Calculate the shortest path between each OD of network;Subsequent cycle distributes next equal portions by the shortest path after updating The OD volume of traffic;
Straight average method, continuous adjustment has been assigned to the volume of traffic on each section and has moved closer to or arrive balance point Join, often walking in circulation, carry out a 0-1 distribution according to the volume of traffic being assigned on each section and obtain one Organize the additional traffic amount in each section;Then obtain with in the distribution volume of traffic in section each in last circulation and this circulation To additional traffic amount be weighted the section distribution volume of traffic that is average and that obtain in this circulation;When adjacent two When the distribution volume of traffic in circulation trends towards certain definite value time (this definite value can determine according to actual state), i.e. Can stop calculating;The distribution volume of traffic obtained in last circulation is i.e. the final volume of traffic.
CN201410410008.5A 2014-08-19 2014-08-19 Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD CN104183119B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410410008.5A CN104183119B (en) 2014-08-19 2014-08-19 Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410410008.5A CN104183119B (en) 2014-08-19 2014-08-19 Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD

Publications (2)

Publication Number Publication Date
CN104183119A CN104183119A (en) 2014-12-03
CN104183119B true CN104183119B (en) 2016-08-24

Family

ID=51964122

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410410008.5A CN104183119B (en) 2014-08-19 2014-08-19 Based on the anti-arithmetic for real-time traffic flow distribution forecasting method pushed away of section OD

Country Status (1)

Country Link
CN (1) CN104183119B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104731963B (en) * 2015-04-03 2018-06-15 重庆邮电大学 Recommend method and system in a kind of gridding path based on car networking
CN105006148B (en) * 2015-07-06 2017-10-31 同济大学 A kind of crossing turns to vehicle number estimation method and system
CN106504535B (en) * 2016-11-30 2018-10-12 东南大学 A kind of trip distribution modeling method of combination Gravity Models and Fratar models
CN108257378A (en) * 2016-12-29 2018-07-06 杭州海康威视数字技术股份有限公司 traffic flow forecasting method and device
CN106951549A (en) * 2017-03-27 2017-07-14 重庆邮电大学 A kind of passenger's traffic path recognition methods based on track IC-card and mobile phone signaling data
CN107038862A (en) * 2017-05-04 2017-08-11 东南大学 A kind of traffic assignation modification method under the conditions of tail number restricted driving
CN107294794A (en) * 2017-07-31 2017-10-24 国网辽宁省电力有限公司 A kind of large-scale ip communication service matrix estimation method and system
CN109493449A (en) * 2018-11-23 2019-03-19 北京航空航天大学 A kind of lorry loading method for estimating state based on lorry GPS track data and high speed transaction data

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3414843B2 (en) * 1993-06-22 2003-06-09 三菱電機株式会社 Transportation control device
DE10022812A1 (en) * 2000-05-10 2001-11-22 Daimler Chrysler Ag A method for determining the traffic situation on the basis of reporting vehicle data for a traffic network with traffic-controlled network nodes
CN101436345B (en) * 2008-12-19 2010-08-18 天津市市政工程设计研究院 System for forecasting harbor district road traffic requirement based on TransCAD macroscopic artificial platform
CN101930668B (en) * 2009-04-29 2012-08-22 上海电器科学研究所(集团)有限公司 Processing method of road traffic OD (Optical Density) information collection system for license plate recognition
US8798897B2 (en) * 2010-11-01 2014-08-05 International Business Machines Corporation Real-time traffic analysis through integration of road traffic prediction and traffic microsimulation models
CN102682598B (en) * 2012-05-31 2014-03-26 山东交通学院 Traffic information distribution method capable of eliminating oscillation

Also Published As

Publication number Publication date
CN104183119A (en) 2014-12-03

Similar Documents

Publication Publication Date Title
Ghosh et al. Multivariate short-term traffic flow forecasting using time-series analysis
US8755991B2 (en) Method and structure for vehicular traffic prediction with link interactions and missing real-time data
Ben-Akiva et al. Dynamic model of peak period traffic congestion with elastic arrival rates
Chen et al. Freeway performance measurement system: mining loop detector data
Hofleitner et al. Arterial travel time forecast with streaming data: A hybrid approach of flow modeling and machine learning
Yuan et al. T-finder: A recommender system for finding passengers and vacant taxis
US9424515B2 (en) Predicting taxi utilization information
US20060178806A1 (en) Method and apparatus for predicting future travel times over a transportation network
Wey Model formulation and solution algorithm of traffic signal control in an urban network
De Palma et al. METROPOLIS: Modular system for dynamic traffic simulation
CN101373559A (en) Method for evaluating city road net traffic state based on floating vehicle data
CN101038700A (en) Mixed controlling method of single dot signal controlling crossing
CN104157139B (en) A kind of traffic congestion Forecasting Methodology and method for visualizing
Osorio et al. A surrogate model for traffic optimization of congested networks: an analytic queueing network approach
Christofa et al. Person-based traffic responsive signal control optimization
CN102081846A (en) Expressway charge data track matching based traffic state recognition method
El Esawey et al. Development of daily adjustment factors for bicycle traffic
Chen et al. Design of limited-stop bus service with capacity constraint and stochastic travel time
CN103413443B (en) Short-term traffic flow forecasting method based on hidden Markov model
CN102044149B (en) City bus operation coordinating method and device based on time variant passenger flows
CN103730008A (en) Bus congestion degree analysis method based on real-time data of bus GPS (Global Position System) and IC (Integrated Circuit) cards
CN103245347A (en) Intelligent navigation method and system based on road condition prediction
Moreira-Matias et al. On predicting the taxi-passenger demand: A real-time approach
CN103761430A (en) Method for identifying peak periods of road networks on basis of floating cars
CN103489335B (en) A kind of Solving Flight Conflicts method of multiple-objection optimization

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