CN109741599A - Traffic circulation evaluation method - Google Patents
Traffic circulation evaluation method Download PDFInfo
- Publication number
- CN109741599A CN109741599A CN201811621961.9A CN201811621961A CN109741599A CN 109741599 A CN109741599 A CN 109741599A CN 201811621961 A CN201811621961 A CN 201811621961A CN 109741599 A CN109741599 A CN 109741599A
- Authority
- CN
- China
- Prior art keywords
- road
- index
- congestion
- crossing
- street
- 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.)
- Pending
Links
Landscapes
- Traffic Control Systems (AREA)
Abstract
A kind of traffic circulation evaluation method, acquire road information, information of vehicles, evaluation data are calculated according to road information, information of vehicles, Traffic Evaluation is carried out according to evaluation data, the beneficial effect is that: on the one hand traffic circulation can be evaluated and accurate characteristic evidences are provided, on the other hand there is stronger adaptability, can satisfy the traffic flow evaluation under complicated road network and the larger situation of traffic density.
Description
Technical field
The present invention relates to traffic circulations to evaluate field, especially a kind of traffic circulation evaluation method.
Background technique
In field of traffic, it includes handing in traffic flow that traffic flow, which is that traffic programme, traffic control etc. provide important foundation,
Flux, traffic density and traffic speed in existing method, mainly moor model etc. using Greenshields model, Green.
Existing method cannot be built upon on the basis of ideal road, not be able to satisfy the demand of actual traffic state,
But also the assessment of the traffic flow in the lesser situation of traffic density is required, bad adaptability, when road network is complicated, traffic density is huge
In the case where, existing method then cannot carry out accurate evaluation to traffic flow modes.
Summary of the invention
The purpose of the present invention is to solve the above problems, devise a kind of traffic circulation evaluation method.Specific design side
Case are as follows:
A kind of traffic circulation evaluation method acquires road information, information of vehicles, is calculated according to road information, information of vehicles
Data are evaluated out, carry out Traffic Evaluation according to evaluation data, the road information includes crossing length, and the information of vehicles includes
Car data is crossed, the evaluation data include region congestion index, region delay index, region unbalance index, operation vehicle delay
The calculating of index, the evaluation data includes being based on average travel speedDetermining road section traffic volume Operation class, section are comprehensive
Mileage ratio Ls, based on area road rank, road section traffic volume Operation class, the comprehensive mileage ratio L in sectionsDetermining region congestion index,
Region postpones index, region unbalance index, commerial vehicle index, and the Traffic Evaluation includes data normalization processing step, meter
Calculate related data matrix step, the characteristic value for calculating related data matrix and feature vector step, calculate main composition contribution rate and
Contribution rate of accumulative total step is converted into comprehensive score step.
In the calculating of the evaluation data, the average travel speed of the region congestion indexCalculation formula are as follows:
Wherein, A, B are respectively to press two successive crossings of driving direction；
L is the distance between two crossings A, B；
N is the vehicles number for driving towards the crossing B in statistics time interval from the crossing A；
ti APass through the time at the crossing A for i-th vehicle；
ti BPass through the time at the crossing B for i-th vehicle,
The road section traffic volume Operation class include through street jam level, trunk roads jam level, secondary distributor road jam level,
Branch jam level,
The road section traffic volume Operation class includes unobstructed, substantially unobstructed, slight congestion, moderate congestion, heavy congestion,
The unobstructed grade in through street is average travel speed
The substantially unobstructed grade in through street is average travel speed
Jam level that through street is slight is average travel speed
Through street moderate jam level is average travel speed
Through street heavy congestion grade is average travel speed
The unobstructed grade in mobile jib road is average travel speed
The unobstructed grade of mobile jib roadbed sheet is average travel speed
Jam level that mobile jib road is slight is average travel speed
Mobile jib road moderate jam level is average travel speed
Mobile jib road heavy congestion grade is average travel speed
The secondary unobstructed grade in pole line is average travel speed
The substantially unobstructed grade in secondary pole line is average travel speed
The secondary slight jam level in pole line is average travel speed
Secondary pole line moderate jam level is average travel speed
Secondary pole line heavy congestion grade is average travel speed
The unobstructed grade of branch is average travel speed
The substantially unobstructed grade of branch is average travel speed
Jam level that branch is slight is average travel speed
Branch moderate jam level is average travel speed
Branch heavy congestion grade is average travel speed
The comprehensive mileage ratio L in the sectionsIn calculating, including
Through street integrates mileage ratioLength is Lf；
Trunk roads integrate mileage ratioLength is Lm；
Secondary distributor road integrates mileage ratioLength is Lc；
Branch integrates mileage ratioLength is Lz；
The same section total length is L,
The through street, secondary distributor road, includes unobstructed Operation class road section length L in branch at trunk roads1, it is substantially unobstructed
Operation class road section length L2, slight congestion Operation class road section length L3, moderate congestion Operation class road section length L4, it is serious
Congestion Operation class road section length L5,
I.e. through street includes:
The unobstructed Operation class road section length in through street
The substantially unobstructed Operation class road section length in through street
Congestion Operation class road section length that through street is slight
Through street moderate congestion Operation class road section length
Through street heavy congestion Operation class road section length
Trunk roads include:
The unobstructed Operation class road section length of trunk roads
The substantially unobstructed Operation class road section length of trunk roads
Congestion Operation class road section length that trunk roads are slight
Trunk roads moderate congestion Operation class road section length
Trunk roads heavy congestion Operation class road section length
Secondary distributor road includes:
The unobstructed Operation class road section length of secondary distributor road
The substantially unobstructed Operation class road section length of secondary distributor road
Congestion Operation class road section length that secondary distributor road is slight
Secondary distributor road moderate congestion Operation class road section length
Secondary distributor road heavy congestion Operation class road section length
Branch includes:
The unobstructed Operation class road section length of branch
The substantially unobstructed Operation class road section length of branch
Congestion Operation class road section length that branch is slight
Branch moderate congestion Operation class road section length
Branch heavy congestion Operation class road section length
Unobstructed Operation class section mileage ratio is Lu, including the unobstructed Operation class section mileage ratio in through streetTrunk roads
Unobstructed Operation class section mileageThe unobstructed Operation class section mileage of secondary distributor roadThe unobstructed Operation class section mileage of branch
Substantially unobstructed Operation class section mileage ratio is Lb, including the substantially unobstructed Operation class section mileage ratio in through streetThe substantially unobstructed Operation class section mileage of trunk roadsThe substantially unobstructed Operation class section mileage of secondary distributor roadZhi Luji
This unobstructed Operation class section mileage
Slight congestion Operation class section mileage ratio is Ll, including through street slight congestion Operation class section mileage ratioTrunk roads are slight congestion Operation class section mileageSecondary distributor road is slight congestion Operation class section mileageBranch is light
Spend congestion Operation class section mileage
Moderate congestion Operation class section mileage ratio is Lm, including through street moderate congestion Operation class section mileage ratioTrunk roads moderate congestion Operation class section mileageSecondary distributor road moderate congestion Operation class section mileageIn branch
Spend congestion Operation class section mileage
Heavy congestion Operation class section mileage ratio is Lg, including through street heavy congestion Operation class section mileage ratioTrunk roads heavy congestion Operation class section mileageSecondary distributor road heavy congestion Operation class section mileageBranch
Heavy congestion Operation class section mileage
Section integrates mileage ratio Ls=0.1*Lu+0.3*Lb+0.5*Ll+0.7*Lm+0.9*Lg；
Wherein,
Through street integrates mileage ratio
Trunk roads integrate mileage ratio
Secondary distributor road integrates mileage ratio
Secondary distributor road integrates mileage ratio
In the calculating of the evaluation data, road weight VKT, the weight VKT=L*n of every road are calculated, wherein L
For road section length, n is the equivalent narrow-gauge locomotive quantity for passing through stop line in statistics time interval, then through street road weight
Wherein k is the quantity of through street
Similarly, trunk roads weight is
Secondary distributor road weight is
Branch weight is
Calculate through street weight proportion are as follows:
Trunk roads weight proportion are as follows:
Secondary distributor road weight proportion are as follows:
Secondary distributor road weight proportion are as follows:
Calculate road network congestion mileage ratio:
TPI=10VKT '
Road network congestion mileage ratio be VKT ' [0,20%) when, road grid traffic run index TPI be [0,2), road grid traffic fortune
Row grade is usual；
Road network congestion mileage ratio be VKT ' [20%, 40%) when being, road grid traffic run index TPI be [2,4), road network is handed over
Logical Operation class is substantially usually；
Road network congestion mileage ratio be VKT ' [40%, 60%) when being, road grid traffic run index TPI be [4,6), road network is handed over
Logical Operation class is slight congestion；
Road network congestion mileage ratio be VKT ' [60%, 80%) when being, road grid traffic run index TPI be [6,8), road network is handed over
Logical Operation class is moderate congestion；
When road network congestion mileage ratio VKT' >=80%, road grid traffic run index TPI be [8,10), road grid traffic operation etc.
Grade is moderate congestion.
During the congestion delay index calculates, the through street congestion delay formula of index are as follows:
WhereinIt goes on a journey average hourage for resident is practical,It is (daily for the average hourage under freestream conditions
1:00 AM to 4 points), start crossing A and terminate crossing B, spacing is L in statistics time interval, is known based on the alert bayonet of electricity
It is n that other data analysis vehicle, which drives towards B quantity from A, and the i-th vehicle is by the time at the crossing AWhat it is at the crossing B is by the time
I-th vehicle is being by the time at the crossing A under stream modeWhat it is at the crossing B is by the time
Vehicle is not generally distinguished to section delay index calculating.But when any in calculating commerial vehicle rule track
When the delay index in section, the hourage of selection must be the hourage generated by target commerial vehicle,
Through street congestion delay index D is calculated separately according to above-mentioned formulaf, trunk roads congestion be delayed index Dm, secondary distributor road gathers around
Stifled delay index Dc, branch congestion be delayed index Dz,
There are average value the gathering around as the rank section for when a plurality of same level section, calculating Multiple Sections congestion delay index
Stifled delay index.
During the congestion delay index calculates,
Region road network congestion delay index calculation formula be
D=Mf*Df+Mm*Dm+Mc*Dc+Mz*Dz
Wherein, MfThe weight coefficient for being through street in urban road network；
MmThe weight coefficient for being through street in urban road network；
McThe weight coefficient for being through street in urban road network；
MzThe weight coefficient for being through street in urban road network；
Its calculation formula calculation formula are as follows:
Wherein, QfFor the intraday effect in through street 1 year,
QmFor the intraday effect in through street 1 year；
QcFor the intraday effect in through street 1 year；
QzFor the intraday effect in through street 1 year.
In the calculating of the region unbalance index, the TPI index of crossing input section x, in same crossing, the road are calculated
The unbalance index of mouth is the standard deviation for one group of number that the traffic circulation index for inputting section by every is constituted, its calculation formula is:
Wherein i is i-th of crossing, and n is n-th of input section at i-th of crossing, and k is the input section in i-th of crossing
Quantity,
City area-traffic unbalance index calculation formula are as follows:
Wherein d is the crossing quantity in region, and Mi is weight coefficient of i-th of crossing in urban road network.
In the Traffic Evaluation step,
In the data normalization processing step, original matrix are as follows:
Wherein, n is sample size, and p is index quantity,
In order to realize the standardization of sample data, the mean value and variance of sample data should be asked.The standardization of sample number is base
It is carried out in the mean value and variance of data, because in practical applications, the dimension that often there is index is different, so in calculating
The preceding influence that need to first eliminate dimension, and initial data standardization is made into standardized transformation to each index component, transformation for mula
Are as follows:
Wherein:For sample average, SjFor sample standard deviation, calculation formula are as follows:
Obtain canonical matrix are as follows:
It calculates in related data matrix step,
For given n sample (region), the related coefficient in (region) between sample is sought.It is each in correlation matrix
A element is expressed as by corresponding related coefficient:
Wherein, rij(i, j=1,2 ..., p) is primal variable xiWith xjRelated coefficient, its calculation formula is:
Because R is real symmetric matrix (i.e. rij=rji), triangle element thereon or lower triangle element need to be only calculated i.e.
It can.
In the eigen vector step for calculating correlation matrix,
Solution characteristic equation | λ I-R |=0 (E is unit matrix) finds out eigenvalue λi(i=1,2 ..., p), and make it by big
Small sequence arrangement, i.e. λ1≥λ2>=..., >=λp>=0, corresponding standardized feature vector γ '1、γ'2.......γ'pMake respectively
For coefficient vector.S1=γ '1X,S2=γ '2X,.....Sp=γ 'pX is referred to as the first principal component of random vector X, second
Principal component ..., pth principal component
It calculates in principal component contributor rate and accumulative contribution process,
Principal component YiContribution rate:Contribution rate of accumulative total:
Contribution rate of accumulative total is generally taken to reach the eigenvalue λ of 85-95%1,λ2,…λm, corresponding the first, the second ... ..., m
(m≤p) a principal component.
It calculates in comprehensive principal component scores index step:
It is converted into comprehensive score step,
Score y=12.78*x+55 is evaluated in traffic flow
Wherein x is comprehensive principal component scores index.
The traffic circulation evaluation method that above-mentioned technical proposal through the invention obtains, the beneficial effect is that:
On the one hand traffic circulation can be evaluated and accurate characteristic evidences is provided, on the other hand there is stronger adaptability,
It can satisfy the traffic flow evaluation under complicated road network and the larger situation of traffic density.
Detailed description of the invention
Fig. 1 is the flow chart of traffic circulation evaluation method of the present invention.
Specific embodiment
The present invention is specifically described with reference to the accompanying drawing.
When original congestion index calculation method is smaller at the regional level in industry, it may appear that congestion index is in maximum value and minimum
Hopping phenomenon between value.This calculation method optimizes former congestion index calculation method, compensates for and gathers around when regional scope position is smaller
Not the shortcomings that stifled index meter is not calculated accurately really and is jumped.
At present in transportation industry to unbalance index without clear calculation method, this unbalance index algorithm is general using mathematical standard difference
It reads, science innovation road network unbalance calculation method.
Embodiment 1
Step1 calculates road-section average and forms speed, and such as two crossings trunk roads A, B, spacing is 1000 meters, based on the alert card of electricity
Mouthful identification data analysis there are 2 vehicles to drive towards B from A, the crossing A by the time be respectively 10 points 30 minutes and 10 points 31 minutes, at the crossing B
By the time be respectively 10 points 32 minutes and 10 points 32 minutes, then A drives towards the average travel speed of B road in the time interval are as follows:
Step2 runs equivalent partition table according to road section traffic volume and show that the road section traffic volume state grade is unimpeded
Step3 tells about method according to Step1 to Step2, determines that all sections are trunk roads, traffic circulation etc. in region
Grade is unimpeded
Trunk roads are in the section mileage ratio L of unimpeded Operation class in the region Step4u=1, in substantially unimpeded operation etc.
The section mileage ratio L of gradeb=0, the section mileage ratio L in slight congestion Operation classl=0, it is in moderate congestion Operation class
Section mileage ratio Lm=0, the section mileage ratio L in heavy congestion Operation classg=0.
Trunk roads integrate mileage ratio
Step5 calculates each grade road weight VKT (truck kilometer number indicates the aggregate-value of VMT Vehicle-Miles of Travel) and VKT ratio
Example, because simplifying all trunk roads in region, then trunk roads road weight is that other road weights are 0
Step6 calculates road network congestion mileage ratio
Step7 show that road traffic operation refers to according to road network congestion mileage ratio table corresponding with road grid traffic operation index conversion
Number (i.e. congestion index) is 1
Step8 show that the Regional Road Network operation level is unimpeded according to road grid traffic Operation class division table.
Embodiment 2
Step1 calculates road-section average travel speed, such as two crossings trunk roads C, D, spacing
It is 500 meters, there are 2 vehicles to drive towards D from C based on the alert bayonet identification data analysis of electricity, be respectively by the time at the crossing C
11 points 30 minutes and 11 points 31 minutes, the crossing D by the time be respectively 11 points 32 minutes and 11 points 32 minutes, then A is sailed in the time interval
To the average travel speed of B road are as follows:
Step2 runs equivalent partition table according to road section traffic volume and show that the road section traffic volume state grade is slight congestion
Step3 tells about method according to Step1 to Step2, determines that all sections are trunk roads, traffic circulation etc. in region
Grade is slight congestion
Trunk roads are in the section mileage ratio L of unimpeded Operation class in the region Step4u=0, in substantially unimpeded operation etc.
The section mileage ratio L of gradeb=0, the section mileage ratio L in slight congestion Operation classl=1, it is in moderate congestion Operation class
Section mileage ratio Lm=0, the section mileage ratio L in heavy congestion Operation classg=0.
Trunk roads integrate mileage ratio
Step5 calculates each grade road weight VKT (truck kilometer number indicates the aggregate-value of VMT Vehicle-Miles of Travel) and VKT ratio
Example, because simplifying all trunk roads in region, then trunk roads road weight is 1, other road weights are 0
Step6 calculates road network congestion mileage ratio
Step7 show that road grid traffic operation refers to according to road network congestion mileage ratio table corresponding with road grid traffic operation index conversion
Number (i.e. congestion index) TPI is 5
Step8 show that the Regional Road Network operation level is slight congestion according to road grid traffic Operation class division table
Embodiment 3
Step1 calculates road-section average travel speed, and such as two crossings trunk roads E, F, spacing is 200 meters, based on the alert card of electricity
Mouthful identification data analysis there are 2 vehicles to drive towards F from E, the crossing E by the time be respectively 11 points 30 minutes and 11 points 31 minutes, at the crossing F
By the time be respectively 11 points 32 minutes and 11 points 32 minutes, then A drives towards the average travel speed of B road in the time interval are as follows:
Step2 runs equivalent partition table according to road section traffic volume and show that the road section traffic volume state grade is heavy congestion
Step3 tells about method according to Step1 to Step2, determines that all sections are trunk roads, traffic circulation etc. in region
Grade is heavy congestion
Trunk roads are in the section mileage ratio L of unimpeded Operation class in the region Step4u=0, in substantially unimpeded operation etc.
The section mileage ratio L of gradeb=0, the section mileage ratio L in slight congestion Operation classl=0, it is in moderate congestion Operation class
Section mileage ratio Lm=0, the section mileage ratio L in heavy congestion Operation classg=1.
Trunk roads integrate mileage ratio
Step5 calculates each grade road weight VKT (truck kilometer number indicates the aggregate-value of VMT Vehicle-Miles of Travel) and VKT ratio
Example, because simplifying all trunk roads in region, then trunk roads road weight is 1, other road weights are 0
Step6 calculates road network congestion mileage ratio
Step7 show that road grid traffic operation refers to according to road network congestion mileage ratio table corresponding with road grid traffic operation index conversion
Number (i.e. congestion index) TPI is 9
Step8 show that the Regional Road Network operation level is heavy congestion according to road grid traffic Operation class division table.
Embodiment 4
Step1 standardizes area data
Input data such as following table
| Input matrix | Index 1 | Index 2 | Index 3 | Index 4 | Index 5 |
| Region 1 | -10 | -10 | -10 | -10 | 0 |
| Region 2 | -5 | -5.5 | -5 | -5.5 | 40 |
| Region 3 | 0 | -1 | 0 | -1 | 80 |
| Region 4 | -8 | -7 | -8 | -6 | 20 |
| Region 5 | -4 | -3.5 | -3 | -2.5 | 45 |
| Region 6 | -0.5 | -1.8 | -1 | -1.5 | 60 |
Standardization is done to matrix and obtains new matrix
| Normalized matrix | Index 1 | Index 2 | Index 3 | Index 4 | Index 5 |
| Region 1 | -1.49 | -1.68 | -1.53 | -1.78 | -1.58 |
| Region 2 | -0.11 | -0.23 | -0.14 | -0.35 | -0.03 |
| Region 3 | 1.26 | 1.23 | 1.25 | 1.09 | 1.51 |
| Region 4 | -0.94 | -0.71 | -0.97 | -0.51 | -0.80 |
| Region 5 | 0.16 | 0.42 | 0.42 | 0.61 | 0.16 |
| Region 6 | 1.12 | 0.97 | 0.97 | 0.93 | 0.74 |
The covariance matrix of Step2 normalized matrix
| 1.2 | 1.18 | 1.19 | 1.14 | 1.18 |
| 1.18 | 1.2 | 1.19 | 1.19 | 1.17 |
| 1.19 | 1.19 | 1.2 | 1.16 | 1.18 |
| 1.14 | 1.19 | 1.16 | 1.2 | 1.14 |
| 1.18 | 1.17 | 1.18 | 1.14 | 1.2 |
Step3 calculates eigenvalue matrix and eigenvectors matrix
Eigenvalue matrix
| 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0000 |
| 0.0000 | 0.0067 | 0.0000 | 0.0000 | 0.0000 |
| 0.0000 | 0.0000 | 0.0276 | 0.0000 | 0.0000 |
| 0.0000 | 0.0000 | 0.0000 | 0.0806 | 0.0000 |
| 0.0000 | 0.0000 | 0.0000 | 0.0000 | 5.8851 |
Eigenvectors matrix
| -0.22 | 0.57 | -0.53 | 0.39 | 0.45 |
| 0.86 | 0.13 | 0.01 | -0.22 | 0.45 |
| -0.08 | -0.81 | -0.35 | 0.15 | 0.45 |
| -0.44 | 0.10 | 0.11 | -0.76 | 0.44 |
| -0.13 | 0.01 | 0.77 | 0.44 | 0.45 |
Step4 show that principal component number is 1 according to principal component contributor rate and contribution rate of accumulative total
Principal component characteristic value is that 5.8851 corresponding feature vectors are last column of feature vector
Step5 calculates six regional complex results by principal component
-3.61、-0.38、2.84、-1.76、0.79、2.12
Step6 calculates six region final score results
8.93、50.10、91.27、32.49、65.12、82.10。
Above-mentioned technical proposal only embodies the optimal technical scheme of technical solution of the present invention, those skilled in the art
The principle of the present invention is embodied to some variations that some of them part may be made, belongs to the scope of protection of the present invention it
It is interior.
Claims (8)
1. a kind of traffic circulation evaluation method acquires road information, information of vehicles, is calculated according to road information, information of vehicles
Data are evaluated, carry out Traffic Evaluation according to evaluation data, which is characterized in that the road information includes crossing length, the vehicle
Information included car data, and the evaluation data include region congestion index, region delay index, region unbalance index, fortune
The calculating of driving delay index, the evaluation data includes being based on average travel speedDetermining road section traffic volume operation etc.
Grade, the comprehensive mileage ratio L in sections, based on area road rank, road section traffic volume Operation class, the comprehensive mileage ratio L in sectionsDetermining area
Domain congestion index, region postpone index, region unbalance index, commerial vehicle index, and the Traffic Evaluation includes data normalization
Processing step, calculate related data matrix step, calculate related data matrix characteristic value and feature vector step, calculate it is main at
Part contribution rate and contribution rate of accumulative total step are converted into comprehensive score step.
2. traffic circulation evaluation method according to claim 1, which is characterized in that in the calculating of the evaluation data,
The average travel speed of the region congestion indexCalculation formula are as follows:
Wherein, A, B are respectively to press two successive crossings of driving direction；
L is the distance between two crossings A, B；
N is the vehicles number for driving towards the crossing B in statistics time interval from the crossing A；
ti APass through the time at the crossing A for i-th vehicle；
ti BPass through the time at the crossing B for i-th vehicle,
The road section traffic volume Operation class includes through street jam level, trunk roads jam level, secondary distributor road jam level, branch
Jam level,
The road section traffic volume Operation class includes unobstructed, substantially unobstructed, slight congestion, moderate congestion, heavy congestion.
3. traffic circulation evaluation method according to claim 2, which is characterized in that the comprehensive mileage ratio L in the sectionsMeter
In calculation, including
Through street integrates mileage ratioLength is Lf；
Trunk roads integrate mileage ratioLength is Lm；
Secondary distributor road integrates mileage ratioLength is Lc；
Branch integrates mileage ratioLength is Lz；
The same section total length is L,
The through street, secondary distributor road, includes unobstructed Operation class road section length L in branch at trunk roads1, substantially unobstructed operation
Grade road section length L2, slight congestion Operation class road section length L3, moderate congestion Operation class road section length L4, heavy congestion
Operation class road section length L5, the comprehensive mileage ratio L in sections=0.1*Lu+0.3*Lb+0.5*Ll+0.7*Lm+0.9*Lg。
4. traffic circulation evaluation method according to claim 3, which is characterized in that the calculating of the evaluation data
In, road weight VKT, the weight VKT=L*n of every road are calculated, wherein L is road section length, and n is to lead in statistics time interval
The equivalent narrow-gauge locomotive quantity of stop line is crossed, then through street road weight
Wherein k is the quantity of through street
Similarly, trunk roads weight is
Secondary distributor road weight is
Branch weight is
Calculate through street weight proportion are as follows:
Trunk roads weight proportion are as follows:
Secondary distributor road weight proportion are as follows:
Secondary distributor road weight proportion are as follows:
Calculate road network congestion mileage ratio:
TPI=10VKT '.
5. traffic circulation evaluation method according to claim 3, which is characterized in that the congestion delay index calculates
In, the through street congestion delay formula of index are as follows:
WhereinFor resident's realityIt goes on a journey average hourage,For the average travelling under freestream conditions
Time starts crossing A and terminates crossing B, spacing is L in statistics time interval, based on the alert bayonet identification data analysis of electricity
It is n that vehicle, which drives towards B quantity from A, and the i-th vehicle is by the time at the crossing AWhat it is at the crossing B is by the timeI-th vehicle is flowing automatically
It is by the time at the crossing A under stateWhat it is at the crossing B is by the time
Through street congestion delay index D is calculated separately according to above-mentioned formulaf, trunk roads congestion be delayed index Dm, secondary distributor road congestion prolongs
When index Dc, branch congestion be delayed index Dz。
6. traffic circulation evaluation method according to claim 5, which is characterized in that the congestion delay index calculates
In,
Region road network congestion delay index calculation formula be
D=Mf*Df+Mm*Dm+Mc*Dc+Mz*Dz
Wherein, MfThe weight coefficient for being through street in urban road network；
MmThe weight coefficient for being through street in urban road network；
McThe weight coefficient for being through street in urban road network；
MzThe weight coefficient for being through street in urban road network；
Its calculation formula calculation formula are as follows:
Wherein, QfFor the intraday effect in through street 1 year,
QmFor the intraday effect in through street 1 year；
QcFor the intraday effect in through street 1 year；
QzFor the intraday effect in through street 1 year.
7. traffic circulation evaluation method according to claim 6, which is characterized in that the calculating of the region unbalance index
In, the TPI index of crossing input section x is calculated, in same crossing, the unbalance index at the crossing is to input section by every
The standard deviation for one group of number that traffic circulation index is constituted, its calculation formula is:
Wherein i is i-th of crossing, and n is n-th of input section at i-th of crossing, and k is the input section number in i-th of crossing
Amount,
City area-traffic unbalance index calculation formula are as follows:
Wherein d is the crossing quantity in region, and Mi is weight coefficient of i-th of crossing in urban road network.
8. traffic circulation evaluation method according to claim 1, which is characterized in that in the Traffic Evaluation step,
In the data normalization processing step, original matrix are as follows:
Wherein, n is sample size, and p is index quantity,
In order to realize the standardization of sample data, the mean value and variance of sample data should be asked.The standardization of sample number is based on number
According to mean value and variance carry out because in practical applications, the dimension that often there is index is different, so needing before the computation
The influence of dimension is first eliminated, and initial data is standardized
Standardized transformation, transformation for mula are made to each index component are as follows:
Wherein:For sample average, SjFor sample standard deviation, calculation formula are as follows:
Obtain canonical matrix are as follows:
It calculates in related data matrix step,
For n given sample, the related coefficient of sample interval is sought.Each of correlation matrix element is by corresponding
Related coefficient is expressed as:
Wherein, rij(i, j=1,2 ..., p) is primal variable xiWith xjRelated coefficient, its calculation formula is:
In the eigen vector step for calculating correlation matrix,
Solution characteristic equation | λ I-R |=0 (E is unit matrix) finds out eigenvalue λi(i=1,2 ..., p), and make its sequence by size
Arrangement, i.e. λ1≥λ2>=..., >=λp>=0, corresponding standardized feature vector γ '1、γ'2.......γ'pRespectively as coefficient
Vector.S1=γ '1X,S2=γ '2X,.....Sp=γ 'pX be referred to as the first principal component of random vector X, second it is main at
Point ..., pth principal component
It calculates in principal component contributor rate and accumulative contribution process,
Principal component YiContribution rate:Contribution rate of accumulative total:
Contribution rate of accumulative total is generally taken to reach the eigenvalue λ of 85-95%1,λ2,…λm, corresponding the first, the second ... ..., m (m≤
P) a principal component.
It calculates in comprehensive principal component scores index step:
It is converted into comprehensive score step,
Score y=12.78*x+55 is evaluated in traffic flow
Wherein x is comprehensive principal component scores index.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811621961.9A CN109741599A (en) | 2018-12-28 | 2018-12-28 | Traffic circulation evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811621961.9A CN109741599A (en) | 2018-12-28 | 2018-12-28 | Traffic circulation evaluation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109741599A true CN109741599A (en) | 2019-05-10 |
Family
ID=66361771
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811621961.9A Pending CN109741599A (en) | 2018-12-28 | 2018-12-28 | Traffic circulation evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109741599A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110223514A (en) * | 2019-06-06 | 2019-09-10 | 北京交通发展研究院 | Urban transportation running state analysis method, apparatus and electronic equipment |
| CN110363989A (en) * | 2019-07-11 | 2019-10-22 | 汉王科技股份有限公司 | Magnitude of traffic flow detection method, device, electronic equipment and storage medium |
| CN110648530A (en) * | 2019-09-18 | 2020-01-03 | 南京信息工程大学 | Car traffic control management method based on working unit zone bit characteristics |
| CN111599171A (en) * | 2020-04-24 | 2020-08-28 | 重庆科技学院 | Intelligent control system, method and storage medium for applying big data to traffic |
| CN114898556A (en) * | 2022-05-11 | 2022-08-12 | 华南理工大学 | Digital road network traffic state calculation method based on multi-scale calculation |
| CN115662143A (en) * | 2022-11-21 | 2023-01-31 | 吉林大学 | Dynamic prediction system and method for operation safety situation of public transport enterprise |
| CN117315934A (en) * | 2023-09-25 | 2023-12-29 | 阜阳交通能源投资有限公司 | Expressway traffic flow real-time monitoring and congestion prediction system based on unmanned aerial vehicle |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102637357A (en) * | 2012-03-27 | 2012-08-15 | 山东大学 | Regional traffic state assessment method |
| CN102819955A (en) * | 2012-09-06 | 2012-12-12 | 北京交通发展研究中心 | Road network operation evaluation method based on vehicle travel data |
| WO2013011796A1 (en) * | 2011-07-20 | 2013-01-24 | 住友電気工業株式会社 | Traffic evaluation device, computer program and traffic evaluation method |
| CN105139645A (en) * | 2015-07-23 | 2015-12-09 | 合肥革绿信息科技有限公司 | Urban regional road network operation index assessment method based on floating car technology |
-
2018
- 2018-12-28 CN CN201811621961.9A patent/CN109741599A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013011796A1 (en) * | 2011-07-20 | 2013-01-24 | 住友電気工業株式会社 | Traffic evaluation device, computer program and traffic evaluation method |
| CN102637357A (en) * | 2012-03-27 | 2012-08-15 | 山东大学 | Regional traffic state assessment method |
| CN102819955A (en) * | 2012-09-06 | 2012-12-12 | 北京交通发展研究中心 | Road network operation evaluation method based on vehicle travel data |
| CN105139645A (en) * | 2015-07-23 | 2015-12-09 | 合肥革绿信息科技有限公司 | Urban regional road network operation index assessment method based on floating car technology |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110223514A (en) * | 2019-06-06 | 2019-09-10 | 北京交通发展研究院 | Urban transportation running state analysis method, apparatus and electronic equipment |
| CN110223514B (en) * | 2019-06-06 | 2020-05-19 | 北京交通发展研究院 | Urban traffic running state analysis method and device and electronic equipment |
| CN110363989A (en) * | 2019-07-11 | 2019-10-22 | 汉王科技股份有限公司 | Magnitude of traffic flow detection method, device, electronic equipment and storage medium |
| CN110648530A (en) * | 2019-09-18 | 2020-01-03 | 南京信息工程大学 | Car traffic control management method based on working unit zone bit characteristics |
| CN111599171A (en) * | 2020-04-24 | 2020-08-28 | 重庆科技学院 | Intelligent control system, method and storage medium for applying big data to traffic |
| CN114898556A (en) * | 2022-05-11 | 2022-08-12 | 华南理工大学 | Digital road network traffic state calculation method based on multi-scale calculation |
| CN114898556B (en) * | 2022-05-11 | 2023-09-29 | 华南理工大学 | Digital road network traffic state calculating method based on multi-scale calculation |
| CN115662143A (en) * | 2022-11-21 | 2023-01-31 | 吉林大学 | Dynamic prediction system and method for operation safety situation of public transport enterprise |
| CN117315934A (en) * | 2023-09-25 | 2023-12-29 | 阜阳交通能源投资有限公司 | Expressway traffic flow real-time monitoring and congestion prediction system based on unmanned aerial vehicle |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109741599A (en) | Traffic circulation evaluation method | |
| CN103366557B (en) | Traffic congestion evaluation method based on congestion index | |
| CN104658252B (en) | Method for evaluating traffic operational conditions of highway based on multisource data fusion | |
| CN105070056A (en) | Intersection traffic congestion index calculation method based on floating car | |
| CN103895649B (en) | A kind of driver safety driving warning method | |
| CN107145989A (en) | Real-road Driving Cycle construction method based on people's car traffic | |
| CN104575050B (en) | A kind of fast road ramp intellectual inducing method and device based on Floating Car | |
| CN105913661A (en) | Highway road section traffic state discrimination method based on charging data | |
| CN107330217A (en) | A kind of middle sight oil consumption Forecasting Methodology based on RBFNN | |
| CN103761589A (en) | Distribution method for urban rail transit | |
| CN109035784A (en) | Dynamic wagon flow OD estimation method based on multi-source heterogeneous data | |
| CN101739823A (en) | Road-section average travel time measuring method suitable for low-frequency sampling | |
| CN106297285A (en) | Freeway traffic running status fuzzy synthetic appraisement method based on changeable weight | |
| CN110363990A (en) | A kind of public transport is passed unimpeded index acquisition methods, system and device | |
| CN108665191A (en) | Subway station outside association area calculation method based on aggregation effect | |
| CN107798867A (en) | A kind of congested traffic network equalization methods based on electric automobile and diesel locomotive mixed flow | |
| CN105551250A (en) | Method for discriminating urban road intersection operation state on the basis of interval clustering | |
| CN104361460A (en) | Road service level evaluation method adopting fuzzy synthetic evaluation method | |
| CN107833459A (en) | A kind of city bus operation conditions evaluation method based on gps data | |
| CN107067729A (en) | A kind of urban road traffic safety state evaluating method | |
| CN105185107A (en) | GPS-based traffic running tendency prediction method | |
| Tao et al. | Evaluation of effectiveness of coordinated signal control in reducing vehicle emissions during peak hours versus nonpeak hours | |
| CN105547310A (en) | Apparatus and method for route planning based on PM2.5 healthy trip | |
| CN110766940A (en) | Method for evaluating running condition of road signalized intersection | |
| CN103258440A (en) | Algorithm for restoring wheel path based on road attributes and real-time road conditions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190510 |