CN111311905A - Particle swarm optimization wavelet neural network-based expressway travel time prediction method - Google Patents

Particle swarm optimization wavelet neural network-based expressway travel time prediction method Download PDF

Info

Publication number
CN111311905A
CN111311905A CN202010068929.3A CN202010068929A CN111311905A CN 111311905 A CN111311905 A CN 111311905A CN 202010068929 A CN202010068929 A CN 202010068929A CN 111311905 A CN111311905 A CN 111311905A
Authority
CN
China
Prior art keywords
data
layer
neural network
output
error
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
Application number
CN202010068929.3A
Other languages
Chinese (zh)
Inventor
于泉
孙瑶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202010068929.3A priority Critical patent/CN111311905A/en
Publication of CN111311905A publication Critical patent/CN111311905A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology

Abstract

The invention relates to a particle swarm optimization-based travel time prediction method for a wavelet neural network. The particle swarm optimization algorithm optimizes parameters of the wavelet neural network through continuous iteration, so that the defects of a wavelet neural network model are overcome, including slow convergence speed, easiness in falling into a local minimum value and easiness in generating an oscillation effect. Through comparison experiment analysis, the particle swarm optimization wavelet neural network model can accurately predict the change trend of the travel time and the fluctuation condition of the travel time, and the advantages of high convergence speed, high prediction precision and strong adaptability are proved.

Description

Particle swarm optimization wavelet neural network-based expressway travel time prediction method
Technical Field
The invention designs a practical model, belongs to the field of intelligent transportation, and particularly relates to a particle swarm optimization wavelet neural network-based expressway travel time prediction method capable of providing trip information with high real-time performance and reliability for travelers.
Background
The accurate prediction of the travel time provides a basis for traffic managers to make management decisions and travelers to reasonably arrange travel plans. Currently, the method of travel time prediction includes: time series models, multiple regression models, Kalman filter prediction, artificial neural networks, and the like. However, the time series model does not consider other factors that affect travel time, but predicts future travel time by analyzing similar trends in past data. The multiple regression model studies data of a target road section and an adjacent road section, and the prediction result is generally not satisfactory. The Kalman filtering model needs to continuously correct the weight in the prediction process, so that the calculation time consumption is too large, and the prediction real-time performance is poor. The artificial neural network is a common method for predicting the travel time, but has the defects of easy falling into a local optimal solution, easy generation of oscillation effect, slow convergence and the like in the application process.
Therefore, when the invention predicts the travel time, the advantages and the disadvantages of the existing method are comprehensively considered, and a new prediction model of the travel time of the expressway is established.
Disclosure of Invention
In view of the defects and research trends in the prior art, the invention aims to provide a particle swarm optimization wavelet neural network-based expressway journey time prediction method model for accurately predicting journey time, thereby providing good theoretical support and decision basis for traffic management and optimization of relevant departments.
In order to realize the purpose of the invention, the adopted technical scheme is as follows:
the method for predicting the travel time of the expressway is researched, and the research method mainly comprises the following steps:
A. processing the highway toll data, and extracting a receipt field required by research;
B. processing the extracted charging data field;
C. on the basis of data processing, predicting the travel time of the highway by using a wavelet neural network model;
D. on the basis of data processing, predicting the travel time of the expressway by utilizing a particle swarm optimization wavelet neural network model;
E. and carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE). According to the comparison result, the wavelet neural network model optimized by the particle swarm can more accurately realize the prediction of the travel time of the expressway, and shows better adaptability in the verification of actual data.
The expressway travel time prediction method is researched, wherein the specific analysis process of the step A is as follows:
A1. the highway toll in China adopts an informatization system which comprehensively covers the toll collection process, so that a large amount of toll collection data can be collected. The data fields required for the study include instantionid, time in, update, time out, as shown in table 1.
Table 1 charging data field description table
Figure BDA0002376793340000021
The expressway journey time prediction method is researched, wherein the specific analysis process in the step B is as follows:
B1. the charging data processing is carried out in two steps, first: eliminating abnormal data including missing data, error data and the like; secondly, the method comprises the following steps: screening effective data according to a quartile method, wherein the effective data refers to a data set which can effectively reflect the whole numerical condition, and the interval of the effective data is [ t ]min,tmax]。
B2. The abnormal data mainly includes the following four types:
lack of entry/exit toll station or entry/exit time information: when the charging system software is abnormal or the line transmission is wrong, the charging system cannot input complete charging data information;
same in-out toll station data: when the driver exchanges the toll card in the service area to avoid the behavior of charging, the vehicle recorded in the toll collection system has the same time at the entrance toll station and the exit toll station;
and recording abnormal time data. Although the charging system is increasingly perfect, the defect that the systems of different charging stations are not synchronous exists, so that the time of different charging stations is deviated, and the phenomenon that the entrance time is later than the exit time is caused. This typically occurs when the charging system updates the date in the morning;
and recording abnormal time data. Due to the particularity of the expressway, a service area and a parking area are generally arranged on the road section, so that the phenomenon that long-distance vehicles stay for a long time occurs.
B3. And screening effective data based on a quartile method on the basis of removing error data. The lower limit and the upper limit of the travel time valid data set are calculated according to the formula:
tmin=t25%-1.5×(t75%-t25%)
tmax=t75%+1.5×(t75%-t25%)
t25%-25% quantile
t75%- - -75% quantile
tmin-a time of flight valid data set lower limit value
tmax-a time of flight valid dataset upper limit value.
The expressway travel time prediction method is researched, wherein the specific analysis process in the step C is as follows:
C1. and (5) constructing a network. And constructing a three-layer wavelet neural network comprising an input layer, a hidden layer and an output layer. As shown in FIG. 1, wherein X1,X2,...,Xm(m represents the number of input neurons) represents the input to the wavelet neural network, Y1,Y2,...,Yn(n represents the number of output neurons), J represents the number of hidden layer nodes, ωijAnd ωjkThe weights between the input layer and the hidden layer, the hidden layer and the output layer are respectively referred to, wherein i is 1, 2. J ═ 1,2,. J; k is 1, 2.
C2. And (5) initializing the network. Randomly initializing the scaling factor a of a wavelet functionjTranslation factor bjConnection weight ωij(input layer and hidden layer) and ωjk(hidden layer and output layer) setting learning rate η of wavelet neural network1And η2
The scaling factor is adjusted as follows:
Figure BDA0002376793340000031
Figure BDA0002376793340000032
-adjusting the pre-scaling factor;
Figure BDA0002376793340000041
-adjusted scaling factor;
η2-learning rate of wavelet parameters, preferably η2=0.01。
The adjustment formula of the translation factor is as follows:
Figure BDA0002376793340000042
Figure BDA0002376793340000043
-adjusting the pre-translational factor value;
Figure BDA0002376793340000044
-adjusted translation factor value;
η2-learning rate of wavelet parameters, preferably η2=0.01。
The input layer and the hidden layer can be modified and adjusted by the following formula:
Figure BDA0002376793340000045
Figure BDA0002376793340000046
-adjusting weights between the previous input layer and the hidden layer;
Figure BDA0002376793340000047
after adjustmentInputting a weight between the layer and the implicit layer;
η1-learning rate of network weight parameters, preferably η1=0.01。
The relationship between the hidden layer and the output layer can be modified by the following formula:
Figure BDA0002376793340000048
Figure BDA0002376793340000049
-adjusting weights between the previous hidden layer and the output layer;
Figure BDA00023767933400000410
-adjusting the weights between the hidden layer and the output layer after adjustment;
η1-learning rate of network weight parameters, preferably η1=0.01。
C3. And (6) outputting the prediction. The training samples are input to the network to calculate the predicted output of the network and to calculate the network output and the expected output error e.
C4. And (4) error calculation. The error e between the predicted output and the expected output of the network is calculated. The specific calculation formula is as follows:
Figure BDA0002376793340000051
yp,k(k) -an actual output value;
tp,k-ideal output value.
C5. And (6) correcting the weight value. In order to enable the error to meet the requirement, a gradient descent method is adopted to correct the weight and the parameters of the wavelet neural network.
C6. If the training times is more than 1000 or e meets the prediction precision requirement (e is less than 0.0001, the error precision is set to be 0.0001, if the error precision is set to be too high, the network convergence speed is slow, and if the error precision is set to be too low, the accuracy of the prediction result is influenced), the training is finished, the prediction result is returned, and otherwise, the learning and the training are continued.
The process of the travel time prediction experiment based on the wavelet neural network is shown in figure 2.
The expressway travel time prediction method is researched, wherein the specific analysis process in the step D is as follows:
D1. network construction: constructing a three-layer POS-WNN neural network comprising an input layer, a hidden layer and an output layer
D2. Parameter initialization processing: the system randomly generates S particles and uses the position vectors of the particles to represent the scale factor and the translation factor a in the wavelet basis functionj,bj(ii) a And a weight value ωij(weight value between ith input and jth hidden layer) and ωj(weight value between jth hidden layer and output layer), the calculation formula is as follows. And meanwhile, setting the maximum and minimum speed, the learning rate and the maximum iteration times of the particles.
xs=(xs,1,...,xs,d,...,xs,D)=(a1,...,aJ,b1,...,bJ11,...,ω1J21,...,ω2J,...,ω291,...,ω29J1,...,ωJ)
D=m×J+J×(m-n)+(m-1)×J
m- - -number of nodes of input layer
J- - -number of hidden layer nodes
n- - -number of output layer nodes
D3. And (5) network training. And calculating the particle fitness in each iteration process based on the error between the actual output value and the ideal output value, wherein the calculation formula is as follows.
Figure BDA0002376793340000061
Q- -total number of training samples
n- - -number of network output neurons
yp,k(k) - -actual output value
tp,k-ideal output value.
D4. Taking the particle fitness as an index for judging whether the particle fitness meets the set error requirement, finishing training if the particle fitness can meet the set error requirement, namely the fitness value is basically kept unchanged, and turning to the step F; if the error requirement is still not met, the next step is carried out.
D5. Judging whether the training times reach the set maximum iteration times or not, if so, jumping out of the loop, turning to the step F, and stopping training; otherwise, the speed and the position of the particles are updated according to the following formula, and the training is continued by returning to the step C.
vs,d(i+1)=ω×vs,d(i)+c1r1[ps,d-xs,d(i)]+c2r2[pg,d-xs,d(i)]
xs,d(i+1)=xs,d(i)+vs,d(i+1)
s=1,2,..,S;d=1,2,...,D
c1,c2-acceleration factor
i- -number of current iterations
Omega-inertia factor
r1,r2-random numbers evenly distributed between 0 and 1.
D6. Network debugging: c, selecting a sample value of a specific time period of the training sample as an input value, taking an actual value of another specific time period in the training sample as an ideal output value, calculating an error according to a fitness calculation formula in the step C, and finishing operation if the set error precision is met; if the set error precision can not be met, the step C is carried out.
The experimental flow of the travel time prediction model based on the particle swarm optimization wavelet neural network is shown in fig. 3.
The expressway journey time prediction method is researched, wherein the specific analysis process of the step E is as follows:
E1. and carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE).
E2. Let tp,kRepresenting actual values of travel time,yp,k(k) The predicted travel time value is represented, and l represents the number of predicted time periods.
E3. Calculating the average absolute error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000071
mean Absolute Error (MAE) -an indicator used to assess the difference between the actual and predicted values.
E4. Calculating the average relative error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000072
mean relative error (MAPE) -an indicator of how accurate the travel time prediction is during the travel time prediction process.
E5. Calculating the mean square error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000073
mean Square Error (MSE) -the degree of change of the overall evaluation data.
And carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE). The experimental results show that: the average absolute error, the average relative error and the mean square error of the prediction result of the particle swarm optimized wavelet neural network travel time prediction model are respectively reduced by 83.36%, 82.20% and 98.15% compared with the wavelet neural network model. The particle swarm optimized wavelet neural network travel time prediction model not only has high prediction precision, but also can more accurately predict the trend and fluctuation condition of the travel time, has certain advantages in the aspect of convergence speed, and has better adaptability.
The invention discloses a particle swarm optimization wavelet neural network-based travel time prediction method. The particle swarm optimization algorithm optimizes parameters of the wavelet neural network through continuous iteration, so that the defects of a wavelet neural network model, including slow convergence speed, easiness in falling into a local minimum value and easiness in generating an oscillation effect, are overcome. Through comparative experiment analysis, the particle swarm optimization wavelet neural network model can accurately predict the change trend of the travel time and can also accurately predict the fluctuation condition of the travel time. The expressway travel time prediction experiment based on the particle swarm optimization wavelet neural network model proves that the particle swarm optimization wavelet neural network has the advantages of high convergence speed, high prediction precision and strong adaptability.
It should be understood that the present invention is not limited to the above preferred embodiments, and any other products in various forms can be obtained by the present invention, but any changes in the shape or structure thereof, which are the same or similar to the technical solution of the present invention, fall within the protection scope of the present invention.
Drawings
FIG. 1 is a block diagram of a compact wavelet neural network
FIG. 2 wavelet neural network flow chart
FIG. 3 is a flow chart of a particle swarm optimization wavelet neural network travel time prediction model experiment
The method for predicting the travel time of the expressway is researched, and the research method mainly comprises the following steps:
A. processing the highway toll data, and extracting a receipt field required by research;
B. processing the extracted charging data field;
C. on the basis of data processing, predicting the travel time of the highway by using a wavelet neural network model;
D. on the basis of data processing, predicting the travel time of the expressway by utilizing a particle swarm optimization wavelet neural network model;
E. and carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE). According to the comparison result, the wavelet neural network model optimized by the particle swarm can more accurately realize the prediction of the travel time of the expressway, and shows better adaptability in the verification of actual data.
The expressway travel time prediction method is researched, wherein the specific analysis process of the step A is as follows:
A1. the highway toll in China adopts an informatization system which comprehensively covers the toll collection process, so that a large amount of toll collection data can be collected. The data fields required for the study include instantionid, time in, update, time out, as shown in table 1.
Table 1 charging data field description table
Figure BDA0002376793340000081
Figure BDA0002376793340000091
The expressway journey time prediction method is researched, wherein the specific analysis process in the step B is as follows:
B1. the charging data processing is carried out in two steps, first: eliminating abnormal data including missing data, error data and the like; secondly, the method comprises the following steps: screening effective data according to a quartile method, wherein the effective data refers to a data set which can effectively reflect the whole numerical condition, and the interval of the effective data is [ t ]min,tmax]。
B2. The abnormal data mainly includes the following four types:
lack of entry/exit toll station or entry/exit time information: when the charging system software is abnormal or the line transmission is wrong, the charging system cannot input complete charging data information;
same in-out toll station data: when the driver exchanges the toll card in the service area to avoid the behavior of charging, the vehicle recorded in the toll collection system has the same time at the entrance toll station and the exit toll station;
and recording abnormal time data. Although the charging system is increasingly perfect, the defect that the systems of different charging stations are not synchronous exists, so that the time of different charging stations is deviated, and the phenomenon that the entrance time is later than the exit time is caused. This typically occurs when the charging system updates the date in the morning;
and recording abnormal time data. Due to the particularity of the expressway, a service area and a parking area are generally arranged on the road section, so that the phenomenon that long-distance vehicles stay for a long time occurs.
B3. And screening effective data based on a quartile method on the basis of removing error data. The lower limit and the upper limit of the travel time valid data set are calculated according to the formula:
tmin=t25%-1.5×(t75%-t25%)
tmax=t75%+1.5×(t75%-t25%)
t25%-25% quantile
t75%- - -75% quantile
tmin-a time of flight valid data set lower limit value
tmax-a time of flight valid dataset upper limit value.
The expressway travel time prediction method is researched, wherein the specific analysis process in the step C is as follows:
C1. and (5) constructing a network. And constructing a three-layer wavelet neural network comprising an input layer, a hidden layer and an output layer. As shown in FIG. 1, wherein X1,X2,...,Xm(m represents the number of input neurons) represents the input to the wavelet neural network, Y1,Y2,...,Yn(n represents the number of output neurons), J represents the number of hidden layer nodes, ωijAnd ωjkThe weights between the input layer and the hidden layer, the hidden layer and the output layer are respectively referred to, wherein i is 1, 2. J ═ 1,2,. J; k is 1, 2.
C2. And (5) initializing the network. Randomly initializing the scaling factor a of a wavelet functionjTranslation factor bjConnection weight ωij(transfusion)In-layer and hidden layer) and ωjk(hidden layer and output layer) setting learning rate η of wavelet neural network1And η2
The scaling factor is adjusted as follows:
Figure BDA0002376793340000101
Figure BDA0002376793340000102
-adjusting the pre-scaling factor;
Figure BDA0002376793340000103
-adjusted scaling factor;
η2-learning rate of wavelet parameters, preferably η2=0.01。
The adjustment formula of the translation factor is as follows:
Figure BDA0002376793340000104
Figure BDA0002376793340000105
-adjusting the pre-translational factor value;
Figure BDA0002376793340000106
-adjusted translation factor value;
η2-learning rate of wavelet parameters, preferably η2=0.01。
The input layer and the hidden layer can be modified and adjusted by the following formula:
Figure BDA0002376793340000111
Figure BDA0002376793340000112
-adjusting weights between the previous input layer and the hidden layer;
Figure BDA0002376793340000113
-the weights between the input layer and the hidden after adjustment;
η1-learning rate of network weight parameters, preferably η1=0.01。
The relationship between the hidden layer and the output layer can be modified by the following formula:
Figure BDA0002376793340000114
Figure BDA0002376793340000115
-adjusting weights between the previous hidden layer and the output layer;
Figure BDA0002376793340000116
-adjusting the weights between the hidden layer and the output layer after adjustment;
η1-learning rate of network weight parameters, preferably η1=0.01。
C3. And (6) outputting the prediction. The training samples are input to the network to calculate the predicted output of the network and to calculate the network output and the expected output error e.
C4. And (4) error calculation. The error e between the predicted output and the expected output of the network is calculated. The specific calculation formula is as follows:
Figure BDA0002376793340000117
yp,k(k) -an actual output value;
tp,k-ideal output value.
C5. And (6) correcting the weight value. In order to enable the error to meet the requirement, a gradient descent method is adopted to correct the weight and the parameters of the wavelet neural network.
C6. If the training times is more than 1000 or e meets the prediction precision requirement (e is less than 0.0001, the error precision is set to be 0.0001, if the error precision is set to be too high, the network convergence speed is slow, and if the error precision is set to be too low, the accuracy of the prediction result is influenced), the training is finished, the prediction result is returned, and otherwise, the learning and the training are continued.
The process of the travel time prediction experiment based on the wavelet neural network is shown in figure 2.
The expressway travel time prediction method is researched, wherein the specific analysis process in the step D is as follows:
D1. network construction: constructing a three-layer POS-WNN neural network comprising an input layer, a hidden layer and an output layer
D2. Parameter initialization processing: the system randomly generates S particles and uses the position vectors of the particles to represent the scale factor and the translation factor a in the wavelet basis functionj,bj(ii) a And a weight value ωij(weight value between ith input and jth hidden layer) and ωj(weight value between jth hidden layer and output layer), the calculation formula is as follows. And meanwhile, setting the maximum and minimum speed, the learning rate and the maximum iteration times of the particles.
xs=(xs,1,...,xs,d,...,xs,D)=(a1,...,aJ,b1,...,bJ11,...,ω1J21,...,ω2J,...,ω291,...,ω29J1,...,ωJ)
D=m×J+J×(m-n)+(m-1)×J
m- - -number of nodes of input layer
J- - -number of hidden layer nodes
n- - -number of output layer nodes
D3. And (5) network training. And calculating the particle fitness in each iteration process based on the error between the actual output value and the ideal output value, wherein the calculation formula is as follows.
Figure BDA0002376793340000121
Q- -total number of training samples
n- - -number of network output neurons
yp,k(k) - -actual output value
tp,k-ideal output value.
D4. Taking the particle fitness as an index for judging whether the particle fitness meets the set error requirement, finishing training if the particle fitness can meet the set error requirement, namely the fitness value is basically kept unchanged, and turning to the step F; if the error requirement is still not met, the next step is carried out.
D5. Judging whether the training times reach the set maximum iteration times or not, if so, jumping out of the loop, turning to the step F, and stopping training; otherwise, the speed and the position of the particles are updated according to the following formula, and the training is continued by returning to the step C.
vs,d(i+1)=ω×vs,d(i)+c1r1[ps,d-xs,d(i)]+c2r2[pg,d-xs,d(i)]
xs,d(i+1)=xs,d(i)+vs,d(i+1)
s=1,2,..,S;d=1,2,...,D
c1,c2-acceleration factor
i- -number of current iterations
Omega-inertia factor
r1,r2-random numbers evenly distributed between 0 and 1.
D6. Network debugging: c, selecting a sample value of a specific time period of the training sample as an input value, taking an actual value of another specific time period in the training sample as an ideal output value, calculating an error according to a fitness calculation formula in the step C, and finishing operation if the set error precision is met; if the set error precision can not be met, the step C is carried out.
The experimental flow of the travel time prediction model based on the particle swarm optimization wavelet neural network is shown in fig. 3.
The expressway journey time prediction method is researched, wherein the specific analysis process of the step E is as follows:
E1. and carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE).
E2. Let tp,kRepresenting the actual value of travel time, yp,k(k) The predicted travel time value is represented, and l represents the number of predicted time periods.
E3. Calculating the average absolute error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000131
mean Absolute Error (MAE) -an indicator used to assess the difference between the actual and predicted values.
E4. Calculating the average relative error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000132
mean relative error (MAPE) -an indicator of how accurate the travel time prediction is during the travel time prediction process.
E5. Calculating the mean square error evaluation index, wherein the specific calculation expression is as follows:
Figure BDA0002376793340000141
mean Square Error (MSE) -the degree of change of the overall evaluation data.
And carrying out comparative analysis on the precision of the 2 prediction models by using 3 evaluation indexes of Mean Absolute Error (MAE), mean relative error (MAPE) and Mean Square Error (MSE). The experimental results show that: the average absolute error, the average relative error and the mean square error of the prediction result of the particle swarm optimized wavelet neural network travel time prediction model are respectively reduced by 83.36%, 82.20% and 98.15% compared with the wavelet neural network model. The particle swarm optimized wavelet neural network travel time prediction model not only has high prediction precision, but also can more accurately predict the trend and fluctuation condition of the travel time, has certain advantages in the aspect of convergence speed, and has better adaptability.

Claims (5)

1. A particle swarm optimization wavelet neural network-based expressway travel time prediction method is characterized by comprising the following steps:
A. processing the highway toll data, and extracting a receipt field required by research;
B. processing the extracted charging data field;
C. on the basis of data processing, predicting the travel time of the highway by using a wavelet neural network model;
D. on the basis of data processing, a particle swarm optimization wavelet neural network model is utilized to predict the travel time of the expressway.
2. The method for predicting the travel time of the highway based on the particle swarm optimization wavelet neural network according to claim 1, wherein the specific analysis process of the step A is as follows:
A1. the highway toll in China adopts an informatization system which comprehensively covers the toll collection process, so that a large amount of toll collection data can be collected; the data fields required for the study include instantionid, time in, EXITSTATION, time out.
3. The method for predicting the travel time of the highway based on the particle swarm optimization wavelet neural network according to claim 1, wherein the specific analysis process of the step B is as follows:
B1. the charging data processing is carried out in two steps, first: eliminating abnormal data including missing data, error data and the like; secondly, the method comprises the following steps: screening effective data according to a quartile method, wherein the effective data refers to a data set reflecting the overall numerical condition, and the interval of the effective data is [ t ]min,tmax];
B2. The anomaly data contains the following four types:
lack of entry/exit toll station or entry/exit time information: when the charging system software is abnormal or the line transmission is wrong, the charging system cannot input complete charging data information;
same in-out toll station data: when the driver exchanges the toll card in the service area to avoid the behavior of charging, the vehicle recorded in the toll collection system has the same time at the entrance toll station and the exit toll station;
recording abnormal time data; although the charging system is increasingly perfect, the defect that the systems of different charging stations are not synchronous exists, so that the time of different charging stations is deviated, and the phenomenon that the entrance time is later than the exit time is caused;
recording abnormal time data; due to the particularity of the expressway, a service area and a parking area are arranged on a road section, so that the phenomenon that long-distance vehicles stay for a long time occurs;
B3. screening effective data based on a quartile method on the basis of removing error data; the lower limit and the upper limit of the travel time valid data set are calculated according to the formula:
tmin=t25%-1.5×(t75%-t25%)
tmax=t75%+1.5×(t75%-t25%)
t25%-25% quantile
t75%- - -75% quantile
tmin-a time of flight valid data set lower limit value
tmax-a time of flight valid dataset upper limit value.
4. The particle swarm optimization wavelet neural network-based expressway travel time prediction method according to claim 1, wherein the specific analysis process of the step C is as follows:
C1. constructing a network; constructing a three-layer wavelet neural network comprising an input layer, a hidden layer and an output layer; as shown in FIG. 1, wherein X1,X2,...,XmRepresenting the input of the wavelet neural network, wherein m represents the number of input neurons; y is1,Y2,...,YnWhere n represents the number of output neurons, J represents the number of hidden layer nodes, ωijAnd ωjkThe weights between the input layer and the hidden layer, the hidden layer and the output layer are respectively referred to, wherein i is 1, 2. J ═ 1,2,. J; k is 1,2,. n;
C2. initializing a network; randomly initializing the scaling factor a of a wavelet functionjTranslation factor bjConnection weight ω between input layer and hidden layerijAnd the connection weight ω between the hidden layer and the output layerjkSetting learning rate η of wavelet neural network1And η2
The scaling factor is adjusted as follows:
Figure FDA0002376793330000021
Figure FDA0002376793330000022
-adjusting the pre-scaling factor;
Figure FDA0002376793330000023
-adjusted scaling factor;
η2-learning rate of wavelet parameters, preferably η2=0.01;
The adjustment formula of the translation factor is as follows:
Figure FDA0002376793330000031
Figure FDA0002376793330000032
-adjusting the pre-translational factor value;
Figure FDA0002376793330000033
-adjusted translation factor value;
η2-learning rate of wavelet parameters, preferably η2=0.01;
The input layer and the hidden layer are modified and adjusted by the following formula:
Figure FDA0002376793330000034
Figure FDA0002376793330000035
-adjusting weights between the previous input layer and the hidden layer;
Figure FDA0002376793330000036
-the weights between the input layer and the hidden after adjustment;
η1-learning rate of network weight parameters, preferably η1=0.01;
The relationship between the hidden layer and the output layer can be modified by the following formula:
Figure FDA0002376793330000037
Figure FDA0002376793330000038
-adjusting weights between the previous hidden layer and the output layer;
Figure FDA0002376793330000039
-adjusting the weights between the hidden layer and the output layer after adjustment;
η1-learning rate of network weight parameters, preferably η1=0.01;
C3. Predicting and outputting; inputting the training samples into the network, thereby calculating the predicted output of the network and calculating the error e between the network output and the expected output;
C4. calculating an error; calculating an error e between the predicted output and the expected output of the network; the specific calculation formula is as follows:
Figure FDA00023767933300000310
yp,k(k) -an actual output value;
tp,k-an ideal output value;
C5. correcting the weight value; correcting the weight and parameters of the wavelet neural network by adopting a gradient descent method;
C6. if the training times are more than 1000 or e meets the prediction precision requirement; e is less than 0.0001, the error precision is set to be 0.0001, the prediction result is returned after the training is finished, and otherwise, the learning and the training are continued.
5. The expressway travel time prediction method based on the particle swarm optimization wavelet neural network according to claim 1,
the specific analysis process of the step D is as follows:
D1. network construction: constructing a three-layer POS-WNN neural network comprising an input layer, a hidden layer and an output layer
D2. Parameter initialization processing: the system randomly generates S particles and uses the position vectors of the particles to represent the scale factor and the translation factor a in the wavelet basis functionj,bj(ii) a And a weight value omega between the ith input and the jth hidden layerijAnd weight value omega between jth hidden layer and output layerjThe calculation formula is as follows; setting the maximum and minimum speed, learning rate and maximum iteration times of the particles;
xs=(xs,1,...,xs,d,...,xs,D)=(a1,...,aJ,b1,...,bJ11,...,ω1J21,...,ω2J,...,ω291,...,ω29J1,...,ωJ)
D=m×J+J×(m-n)+(m-1)×J
m- - -number of nodes of input layer
J- - -number of hidden layer nodes
n- - -number of output layer nodes
D3. Network training; calculating the particle fitness in each iteration process based on the error between the actual output value and the ideal output value, wherein the calculation formula is as follows;
Figure FDA0002376793330000041
q- -total number of training samples
n- - -number of network output neurons
yp,k(k) - -actual output value
tp,k-an ideal output value;
D4. taking the particle fitness as an index for judging whether the particle fitness meets the set error requirement, finishing training if the particle fitness can meet the set error requirement, namely the fitness value is basically kept unchanged, and turning to the step F; if the error requirement is still not met, the next step is carried out;
D5. judging whether the training times reach the set maximum iteration times or not, if so, jumping out of the loop, turning to the step F, and stopping training; otherwise, updating the speed and the position of the particles according to the following formula, and returning to the step C to continue training;
vs,d(i+1)=ω×vs,d(i)+c1r1[ps,d-xs,d(i)]+c2r2[pg,d-xs,d(i)]
xs,d(i+1)=xs,d(i)+vs,d(i+1)
s=1,2,..,S;d=1,2,...,D
c1,c2-acceleration factor
i- -number of current iterations
Omega-inertia factor
r1,r2-random numbers uniformly distributed between 0 and 1;
D6. network debugging: c, selecting a sample value of a specific time period of the training sample as an input value, taking an actual value of another specific time period in the training sample as an ideal output value, calculating an error according to a fitness calculation formula in the step C, and finishing operation if the set error precision is met; if the set error precision can not be met, the step C is carried out.
CN202010068929.3A 2020-01-21 2020-01-21 Particle swarm optimization wavelet neural network-based expressway travel time prediction method Pending CN111311905A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010068929.3A CN111311905A (en) 2020-01-21 2020-01-21 Particle swarm optimization wavelet neural network-based expressway travel time prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010068929.3A CN111311905A (en) 2020-01-21 2020-01-21 Particle swarm optimization wavelet neural network-based expressway travel time prediction method

Publications (1)

Publication Number Publication Date
CN111311905A true CN111311905A (en) 2020-06-19

Family

ID=71144901

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010068929.3A Pending CN111311905A (en) 2020-01-21 2020-01-21 Particle swarm optimization wavelet neural network-based expressway travel time prediction method

Country Status (1)

Country Link
CN (1) CN111311905A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112633592A (en) * 2020-12-30 2021-04-09 鱼快创领智能科技(南京)有限公司 Vehicle frequent-running route calculation method and system based on machine learning clustering algorithm
CN112906993A (en) * 2021-01-12 2021-06-04 西安石油大学 Expressway green traffic station-passing inspection time prediction method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3099934B2 (en) * 1994-09-08 2000-10-16 株式会社東芝 Travel time prediction device
JP3212241B2 (en) * 1995-08-25 2001-09-25 株式会社東芝 Traffic flow prediction device
CN103077610A (en) * 2012-12-31 2013-05-01 清华大学 Road trip time estimating method and system
CN105355049A (en) * 2015-11-05 2016-02-24 北京航空航天大学 Highway running state evaluation method based on macroscopic fundamental diagram
CN105761488A (en) * 2016-03-30 2016-07-13 湖南大学 Real-time limit learning machine short-time traffic flow prediction method based on fusion
CN105913661A (en) * 2016-06-15 2016-08-31 北京航空航天大学 Highway road section traffic state discrimination method based on charging data
CN106372722A (en) * 2016-09-18 2017-02-01 中国科学院遥感与数字地球研究所 Subway short-time flow prediction method and apparatus
CN106971536A (en) * 2017-04-17 2017-07-21 武汉云砥信息科技有限公司 It is a kind of to merge Floating Car and the highway congestion status detection method of trackside video

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3099934B2 (en) * 1994-09-08 2000-10-16 株式会社東芝 Travel time prediction device
JP3212241B2 (en) * 1995-08-25 2001-09-25 株式会社東芝 Traffic flow prediction device
CN103077610A (en) * 2012-12-31 2013-05-01 清华大学 Road trip time estimating method and system
CN105355049A (en) * 2015-11-05 2016-02-24 北京航空航天大学 Highway running state evaluation method based on macroscopic fundamental diagram
CN105761488A (en) * 2016-03-30 2016-07-13 湖南大学 Real-time limit learning machine short-time traffic flow prediction method based on fusion
CN105913661A (en) * 2016-06-15 2016-08-31 北京航空航天大学 Highway road section traffic state discrimination method based on charging data
CN106372722A (en) * 2016-09-18 2017-02-01 中国科学院遥感与数字地球研究所 Subway short-time flow prediction method and apparatus
CN106971536A (en) * 2017-04-17 2017-07-21 武汉云砥信息科技有限公司 It is a kind of to merge Floating Car and the highway congestion status detection method of trackside video

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
姚佳辰: "基于浮动车GPS数据的城市路况判别与短时预测", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
李少华等: "基于改进的PSO算法的小波神经网络", 《科研新视野》 *
陈娇娜: "大数据驱动下的高速公路交通运行状态评价与分析", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112633592A (en) * 2020-12-30 2021-04-09 鱼快创领智能科技(南京)有限公司 Vehicle frequent-running route calculation method and system based on machine learning clustering algorithm
CN112633592B (en) * 2020-12-30 2023-07-18 鱼快创领智能科技(南京)有限公司 Vehicle constant running route calculation method and system based on machine learning clustering algorithm
CN112906993A (en) * 2021-01-12 2021-06-04 西安石油大学 Expressway green traffic station-passing inspection time prediction method

Similar Documents

Publication Publication Date Title
Shu et al. A short-term traffic flow prediction model based on an improved gate recurrent unit neural network
CN109243172B (en) Traffic flow prediction method for optimizing LSTM neural network based on genetic algorithm
CN110070713B (en) Traffic flow prediction method based on bidirectional nested LSTM neural network
CN106981198B (en) Deep learning network model for travel time prediction and establishing method thereof
CN110458336B (en) Online appointment vehicle supply and demand prediction method based on deep learning
WO2023056696A1 (en) Urban rail transit short-term passenger flow forecasting method based on recurrent neural network
CN113096388B (en) Short-term traffic flow prediction method based on gradient lifting decision tree
CN110570651A (en) Road network traffic situation prediction method and system based on deep learning
CN111951553B (en) Prediction method based on traffic big data platform and mesoscopic simulation model
Ma et al. Short-term traffic flow forecasting by selecting appropriate predictions based on pattern matching
Sun et al. Spatio-temporal segmented traffic flow prediction with ANPRS data based on improved XGBoost
CN113538910B (en) Self-adaptive full-chain urban area network signal control optimization method
CN110210648B (en) Gray long-short term memory network-based control airspace strategic flow prediction method
CN110675632B (en) Vehicle short-time trajectory prediction control method aiming at multi-feature space and data sparseness
CN111179592B (en) Urban traffic prediction method and system based on spatio-temporal data flow fusion analysis
CN111027662A (en) SD-LSSVR short-time traffic flow prediction method based on chaos quantum particle swarm optimization
Shao et al. The traffic flow prediction method using the incremental learning-based CNN-LTSM model: the solution of mobile application
CN110991776A (en) Method and system for realizing water level prediction based on GRU network
CN111311905A (en) Particle swarm optimization wavelet neural network-based expressway travel time prediction method
CN113674524A (en) LSTM-GASVR-based multi-scale short-time traffic flow prediction modeling and prediction method and system
CN114510778A (en) Track irregularity prediction method based on hybrid intelligent optimization LSTM
Zou et al. The method of dynamic identification of the maximum speed limit of expressway based on electronic toll collection data
CN113327417B (en) Traffic flow prediction method based on 3D dynamic space-time residual convolution associated network
CN112098869B (en) Self-adaptive electric vehicle SOC estimation method based on big data
CN115438453B (en) Method for constructing road network facility topological structure by using observation data

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200619

WD01 Invention patent application deemed withdrawn after publication