CN102074112A - Time sequence multiple linear regression-based virtual speed sensor design method - Google Patents
Time sequence multiple linear regression-based virtual speed sensor design method Download PDFInfo
- Publication number
- CN102074112A CN102074112A CN 201110033437 CN201110033437A CN102074112A CN 102074112 A CN102074112 A CN 102074112A CN 201110033437 CN201110033437 CN 201110033437 CN 201110033437 A CN201110033437 A CN 201110033437A CN 102074112 A CN102074112 A CN 102074112A
- Authority
- CN
- China
- Prior art keywords
- pseudo
- speed
- node
- traffic flow
- velocity sensor
- 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
Links
Images
Landscapes
- Traffic Control Systems (AREA)
Abstract
The invention discloses a time sequence multiple linear regression-based virtual speed sensor design method, which belongs to the technical field of traffic state acquisition and relates to data processing analysis and data fusion of a traffic flow sensor. The method comprises the following steps of: first, setting a virtual speed sensor between adjacent sensor nodes; then, mapping speed data on a time sequence to a space sequence and training a weight matrix which correlates to the speed of the virtual speed sensor and the speed data obtained by the traffic flow sensor by a least square method; and finally, estimating the speed of each virtual speed sensor through multiple linear regression by utilizing a weight coefficient and the speed data of the traffic flow sensor. The method can acquire the spatial distribution of the speed on a road section, can acquire more accurate traffic information and provides the basis for traffic management and control.
Description
Technical field
The invention belongs to traffic behavior and obtain technical field, be specifically related to the Data Management Analysis of traffic flow sensor and a kind of pseudo-velocity sensor design method of data fusion based on the time series multiple linear regression.
Background technology
Transport information has the time-space domain characteristic, and traffic flow all has continuity in the time-space domain.Transport information is not only relevant with the current and historical traffic behavior in somewhere, and is simultaneously also relevant with the traffic behavior of upstream and downstream.To more reasonably determine transport information, need determine on the highway section continually varying traffic behavior on time and space.But owing to do not have measurement mechanism between the adjacent sensor node, the traffic data of sparse spatial sequence is difficult to directly obtain transport information accurately.
If the traffic behavior of each locus, sensor node place can be determined on the highway section, then when sensor abundant (trending towards infinite), can obtain traffic behavior accurately.Therefore we have defined virtual-sensor node (not actual traffic current detection functionality), if the traffic behavior of each locus, virtual-sensor node place can determine that the traffic behavior that then can obtain the highway section distributes.
Traffic flow all has the characteristic of stream on time and space, so the continuity of the traffic behavior of the spatial sequence middle and lower reaches historical traffic behavior that is upstreams is flowed.Therefore we utilize the speed data on the known sensor node time sequence, and the speed of the pseudo-velocity sensor node of diverse location on the spatial sequence is estimated.Utilize least square method to obtain weight matrix with pseudo-velocity sensor node speed and traffic flow sensor velocity correlation, estimate the speed of pseudo-velocity sensor node, finally can obtain the space distribution of highway section speed by multiple linear regression.
Summary of the invention
In order to overcome the deficiency of the existing pick-up unit quantity of road, the invention provides a kind of pseudo-velocity sensor design method based on the time series multiple linear regression, it is characterized in that, may further comprise the steps:
(1) step that the pseudo-velocity sensor node is demarcated: adopt approach based on linear interpolation between adjacent traffic flow sensor, to set the pseudo-velocity sensor;
(2) according to the speed data of calibrated pseudo-velocity sensor node and traffic flow sensor, the step of training weight matrix: the seasonal effect in time series speed data is mapped to spatial sequence, the average velocity between the adjacent node is adopted to average is similar to, utilizes least square method that weight matrix is trained:
(3) estimate the step of pseudo-velocity sensor node speed: utilize the data on the sensor node time series and the weight matrix of training, estimate the speed of pseudo-velocity sensor node by multiple linear regression.
Node in the described adjacent node comprises traffic flow sensor node and pseudo-velocity sensor node.
The described pseudo-velocity sensor node is demarcated is to suppose that the distance between traffic flow sensor node i and the i+1 is d
i, the distance setting between the adjacent node (node comprises traffic flow sensor and pseudo-velocity sensor) is a
i, then the number of the pseudo-velocity sensor node between traffic flow sensor node i and the i+1 is m
i=[d
i/ a
i].
In the described step (2) the seasonal effect in time series speed data being mapped to spatial sequence is: at t
NConstantly, (i, speed k) is defining virtual speed pickup node
Pseudo-velocity sensor node number between traffic flow sensor node i and the i+1 is m
i, then by
Wherein, n
iNumber for the speed data of the traffic flow sensor i of needs; w
KjBe weight coefficient,
w
Kj≤ 1; w
Kj〉=0; t
NBe current time; Δ t is the time interval of traffic flow sensor acquisition data; v
i(t
N-q Δ t) (q=1,2...n
i) be traffic flow sensor i current time t
NThe q time preceding historical data.
Beneficial effect of the present invention:
The present invention is based on the method for time series multiple linear regression, by between the adjacent sensors node, setting the pseudo-velocity sensor and weight matrix being trained based on least square method, utilize the speed data of real sensor node time sequence, obtain the speed of pseudo-velocity sensor node, realize design the pseudo-velocity sensor node.By the result of this design, can draw the space distribution of traffic behavior, can estimate transport information more accurately, for traveler provides real-time and effective information, help them to carry out routing better, realizing route is induced.
Description of drawings
Fig. 1 pseudo-velocity sensor node is demarcated synoptic diagram.
The time series speed data synoptic diagram of Fig. 2 sensor node i
Embodiment
The invention provides a kind of pseudo-velocity sensor design method based on the time series multiple linear regression, overcome the deficiency of the existing pick-up unit quantity of road, more reasonable aspect traffic flow theory, obtain the speed of pseudo-velocity sensor by the multiple linear regression on the time series, can obtain the space distribution of speed, thereby estimate transport information more accurately, and provide foundation for traffic administration and control.
A kind of pseudo-velocity sensor design method based on the time series multiple linear regression, contain following steps:
(1) the pseudo-velocity sensor node is demarcated.Suppose that the distance between traffic flow sensor node i and the i+1 is d
i(unit, rice), the distance setting between the adjacent node is a
i, then the number of the pseudo-velocity sensor node between traffic flow sensor node i and the i+1 is m
i=[d
i/ a
i]; Node in the described adjacent node comprises pseudo-velocity sensor node and traffic flow sensor node; Be labeled as (i, 1), (i, 2) respectively ... (i, m
i), as shown in Figure 1.
(2), train the step of weight matrix according to the speed data of calibrated pseudo-velocity sensor node and traffic flow sensor
Step 1: the seasonal effect in time series speed data is mapped to spatial sequence
Traffic flow all has the characteristic of stream on time and space, so the continuity of the traffic behavior of the spatial sequence middle and lower reaches historical traffic behavior that is upstreams is flowed.Therefore can utilize the speed data on the known traffic flow sensor node time series, the speed of the pseudo-velocity sensor node of diverse location on the spatial sequence is estimated.
Traffic flow sensor node i (i=1,2 ... n) speed can directly be obtained, so pseudo-velocity sensor node (i, k) (k=1,2 ... m
i) speed can be by the seasonal effect in time series data of sensor node i by calculating.(i, speed k) is the pseudo-velocity sensor node
Wherein, n
iNumber for the speed data of the traffic flow sensor i of needs; w
KjBe weight coefficient,
w
Kj≤ 1; w
Kj〉=0; t
NBe current time; Δ t is the time interval of traffic flow sensor acquisition data; v
i(t
N-q Δ t) (q=1,2 ... n
i) be traffic flow sensor i current time t
NThe q time preceding historical data.The signal of traffic flow sensor node i seasonal effect in time series speed data as shown in Figure 2.
Weight coefficient w
KjCan constitute a weight matrix θ,
θ is relevant with the speed data of pseudo-velocity sensor node speed and traffic flow sensor.
On road conditions and the similar highway section of traffic flow situation, weight matrix is constant.Pseudo-velocity sensor node number between sensor node i and the sensor node i+1 is m
i, can get matrix by (1):
Step 2: adopt the method for averaging to be similar to the average velocity between the adjacent sensors
If enough pseudo-velocity sensors are arranged on the highway section, then the distance between the adjacent node (node comprises traffic flow sensor and pseudo-velocity sensor) will be very short.So, can adopt the method for averaging to be similar to for the average velocity between the adjacent traffic flow sensor, obtain following relation:
Wherein, T
i(t
N) for the start time be t
NSensor node i and the real hourage between the sensor node i+1.
Step 3: utilize least square method that weight matrix is trained
Can get by (2) (3),
Then following formula one has m
i* n
iIndividual known variables is utilized abundant data source, lists the overdetermination matrix equation that constitutes shape such as A θ=b, wherein,
Wherein, l 〉=m
i* n
i
Adopt least square method to estimate weight matrix, make the quadratic sum of error
Minimum, thus draw the estimated value of weight matrix θ.
(3), estimate the step of pseudo-velocity sensor node speed according to weight matrix after the training and the speed data on the traffic flow sensor time sequence
Utilize the data on the sensor node time series and the weight matrix of training, estimate the speed of the pseudo-velocity sensor node of different time by multiple linear regression.
Embodiment 2:
Adopt the Xizhimenwai Dajie to carry out case verification.The sensor placement in this highway section and between distance as shown in Figure 4.Wherein, 41003,41004 is microwave detector, and speed data is provided; 67,68,69 is detecting device hourage, and the journey time of each car process adjacent hourage of detecting device is provided.
Distance between 41003 to 41004 is 1.2020km, the number of setting the pseudo-velocity sensor node is [1202.0/100]=12, distance between the adjacent node (comprising real sensor node and pseudo-velocity sensor node) is 1202.0/12=100.2m, from 41003 to 41004 demarcate successively to (41003,1), (41003,2) ... (41003,12).
The pseudo-velocity number of sensors is 12, need 12 historical speed data to calculate the speed of pseudo-velocity sensor, then per two adjacent nodes have 12 parameters between (comprising real sensor node and pseudo-velocity sensor node), and one has 144 parameters and need determine on 41003 to the 41004 whole highway sections.
The journey time of all vehicles in the same data collection cycle is on average obtained this collection period hourage of the zero hour.With 67 to 68 and 68 to 69 different corresponding constantly hourages the ratio by distance be converted into hourage between this moment speed pickup 41003 to 41004, be designated as T respectively
41003(t
1), T
41003(t
2) ... T
41003(t
144).With these hourages as real hourage.
Can get by (6),
Wherein, m
3=12, v
41003(t
N), v
41004(t
N) (N=1,2 ... 144) be respectively t
N41003 and 41004 speed datas that detect constantly.Extract 2008.08.13 to the data of the 9:00-12:00 of 2008.08.15 and calculate matrix b.
A whenever definite b (i) (i=1,2 ... 144), all in the tabulation of 41003 historical speed from t
NThe zero hour get 12 speed datas forward, be designated as A (i, 1) respectively, A (i, 2) ... A (i, 12) obtains preceding 12 of A matrix and is listed as.Data collection cycle Δ t is 120s.By Algorithm Analysis as can be known, A should be one 144 * 144 matrix, and per 12 row all are the repetitions of preceding 12 row.Thereby determine the A matrix.
Utilize least square method that parameter is trained, the parameter training model is as follows:
Min(Aθ-b)
2
Thereby estimate parameter matrix θ.
Claims (4)
1. the pseudo-velocity sensor design method based on the time series multiple linear regression is characterized in that, may further comprise the steps:
(1) step that the pseudo-velocity sensor node is demarcated: adopt approach based on linear interpolation between adjacent traffic flow sensor, to set the pseudo-velocity sensor;
(2) according to the speed data of calibrated pseudo-velocity sensor node and traffic flow sensor, the step of training weight matrix: the seasonal effect in time series speed data is mapped to spatial sequence, the average velocity between the adjacent node is adopted to average is similar to, utilizes least square method that weight matrix is trained:
(3) estimate the step of pseudo-velocity sensor node speed: utilize the data on the sensor node time series and the weight matrix of training, estimate the speed of pseudo-velocity sensor node by multiple linear regression.
2. the pseudo-velocity sensor design method based on the time series multiple linear regression according to claim 1 is characterized in that the node in the described adjacent node comprises traffic flow sensor node and pseudo-velocity sensor node.
3. the pseudo-velocity sensor design method based on the time series multiple linear regression according to claim 1 is characterized in that, the described pseudo-velocity sensor node is demarcated is to suppose that the distance between traffic flow sensor node i and the i+1 is d
i, the distance setting between the adjacent node (node comprises traffic flow sensor and pseudo-velocity sensor) is a
i, then the number of the pseudo-velocity sensor node between traffic flow sensor node i and the i+1 is m
i=[d
i/ a
i].
4. the pseudo-velocity sensor design method based on the time series multiple linear regression according to claim 1 is characterized in that, in the described step (2) the seasonal effect in time series speed data is mapped to spatial sequence and is: at t
NConstantly, (i, speed k) is defining virtual speed pickup node
Wherein, n
iNumber for the speed data of the traffic flow sensor i of needs; K represents the pseudo-velocity sensor node; w
KjBe weight coefficient,
w
Kj≤ 1; w
Kj〉=0; t
NBe current time; Δ t is the time interval of traffic flow sensor acquisition data; v
i(t
N-q Δ t) (q=1,2...n
i) be traffic flow sensor i current time t
NThe q time preceding historical data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110033437 CN102074112B (en) | 2011-01-30 | 2011-01-30 | Time sequence multiple linear regression-based virtual speed sensor design method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110033437 CN102074112B (en) | 2011-01-30 | 2011-01-30 | Time sequence multiple linear regression-based virtual speed sensor design method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102074112A true CN102074112A (en) | 2011-05-25 |
CN102074112B CN102074112B (en) | 2012-12-26 |
Family
ID=44032636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110033437 Expired - Fee Related CN102074112B (en) | 2011-01-30 | 2011-01-30 | Time sequence multiple linear regression-based virtual speed sensor design method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102074112B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102226932A (en) * | 2011-05-31 | 2011-10-26 | 东北大学 | Vehicle speed measuring method based on radio frequency identification technology |
CN102271348A (en) * | 2011-07-08 | 2011-12-07 | 电子科技大学 | Link quality estimation system and method for cyber physical system |
CN104751630A (en) * | 2014-12-31 | 2015-07-01 | 浙江工业大学 | Road traffic state acquisition method based on Kernel-KNN matching |
CN105632193A (en) * | 2015-12-25 | 2016-06-01 | 银江股份有限公司 | Data missing road speed calculation method based on time-space relevance |
CN105894809A (en) * | 2014-12-25 | 2016-08-24 | 杭州远眺科技有限公司 | Sectional type urban road traffic state estimation method |
CN114283590A (en) * | 2021-09-02 | 2022-04-05 | 青岛海信网络科技股份有限公司 | Traffic flow peak prediction method and device and electronic equipment |
-
2011
- 2011-01-30 CN CN 201110033437 patent/CN102074112B/en not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
《第27届中国控制会议论文集》 20080701 王旭 基于混合PETRI网的列车速度调整的分析 全文 1-4 , * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102226932A (en) * | 2011-05-31 | 2011-10-26 | 东北大学 | Vehicle speed measuring method based on radio frequency identification technology |
CN102271348A (en) * | 2011-07-08 | 2011-12-07 | 电子科技大学 | Link quality estimation system and method for cyber physical system |
CN102271348B (en) * | 2011-07-08 | 2014-03-05 | 电子科技大学 | Link quality estimation system and method for cyber physical system |
CN105894809A (en) * | 2014-12-25 | 2016-08-24 | 杭州远眺科技有限公司 | Sectional type urban road traffic state estimation method |
CN104751630A (en) * | 2014-12-31 | 2015-07-01 | 浙江工业大学 | Road traffic state acquisition method based on Kernel-KNN matching |
CN105632193A (en) * | 2015-12-25 | 2016-06-01 | 银江股份有限公司 | Data missing road speed calculation method based on time-space relevance |
CN105632193B (en) * | 2015-12-25 | 2017-12-22 | 银江股份有限公司 | A kind of shortage of data section speed calculation method based on space-time relationship |
CN114283590A (en) * | 2021-09-02 | 2022-04-05 | 青岛海信网络科技股份有限公司 | Traffic flow peak prediction method and device and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN102074112B (en) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105788252B (en) | Arterial street track of vehicle reconstructing method based on fixed point detector and signal timing dial data fusion | |
CN102074112B (en) | Time sequence multiple linear regression-based virtual speed sensor design method | |
CN104134349B (en) | A kind of public transport road conditions disposal system based on traffic multisource data fusion and method | |
Qiu et al. | Estimation of freeway traffic density with loop detector and probe vehicle data | |
CN102509454B (en) | Road state merging method based on floating car data (FCD) and earth magnetism detector | |
CN104200687B (en) | A kind of driver's speed control behavior monitoring device and monitoring method | |
CN106710215B (en) | Bottleneck upstream lane grade traffic status prediction system and implementation method | |
CN104680789B (en) | Rapid road congestion index estimation and prediction method | |
CN105023434B (en) | Method for obtaining congestion index of motorway | |
CN108492562A (en) | Intersection vehicles trajectory reconstruction method based on fixed point detection with the alert data fusion of electricity | |
Yang et al. | Arterial roadway energy/emissions estimation using modal-based trajectory reconstruction | |
CN104197945A (en) | Global voting map matching method based on low-sampling-rate floating vehicle data | |
CN104778834A (en) | Urban road traffic jam judging method based on vehicle GPS data | |
CN101894461A (en) | Method for predicting travel time on urban ground level roads | |
CN109272746A (en) | A kind of MFD estimating and measuring method based on BP neural network data fusion | |
CN103942957A (en) | Method for calculating signalized intersection vehicle queuing length under saturation condition | |
CN107437339A (en) | Variable information advices plate control method for coordinating and system under a kind of information guidance | |
CN105930614A (en) | Cell transmission model parameter calibration and verification method specific to variable speed limit control | |
Zhu et al. | An automated vehicle fuel economy benefits evaluation framework using real-world travel and traffic data | |
CN105632193A (en) | Data missing road speed calculation method based on time-space relevance | |
CN106355882A (en) | Traffic state estimation method based on in-road detector | |
CN103632537B (en) | A kind of urban road AADT method of estimation based on Floating Car | |
CN103093625A (en) | City road traffic condition real-time estimation method based on reliability verification | |
An et al. | Lane-based traffic arrival pattern estimation using license plate recognition data | |
CN103927437B (en) | The method measuring space headway in non-rectilinear section |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121226 Termination date: 20220130 |