# Two-Phase Model of Multistep Forecasting of Traffic State Reliability

Authors: Jufen Yang; Zhigang Liu; Guiyan Jiang; Lin Zhu
Journal: Discrete Dynamics in Nature and Society
 (2018)

---

## Abstract

Multistep prediction of traffic state is a key technology for advanced transportation information system. The research results based on the principle of multistep prediction can provide more information about the traffic operating quality in advance. Considering that prediction error increases with the increasing numbers of multistep predictions, this research proposes the concept of dynamic predictability that is related to the characteristic of historical traffic flow data used. The traffic flow is characterized by randomness, regularity, and volatility according to the traffic flow theory. Therefore, three key indexes are firstly calculated to measure the characteristics of reliability series. Then a two-phase model is established based on wavelet neural network optimized by particle swarm optimization. The upper phase is a model to estimate the number of predictable steps, and the lower phase is the multistep prediction model of reliability. Compared with that of backpropagation neural network and support vector machine, results show that the convergence time of the wavelet neural network optimized by particle swarm optimization is the lowest, which only costs 256 and 291 seconds in both two-phase models under the same conditions. The average relative error of multistep prediction reached the lowest value, 8.91% and 12.01%, respectively, for weekday and weekend data used. Moreover, the prediction performance based on weekday is better than that of weekend. The research results lay a decision-making basis for managers in determining the key parts of road network to develop future improvement measures.


---

## Body

## 1. Introduction

In recent decades, the increasing demand for road transportation has negatively affected the stability and reliability of traffic operation which caused a series of drawbacks, such as extensive waste of travel time, decreasing of environmental quality, and aggravated vehicle wear and tear. Meanwhile, the requirement for accurate prediction of traffic state is increasing. Prediction results serve not only as an important basis for traffic control and guidance but also as a decision support for travelers to adjust their travel plan.Development of Advanced Traveler Information Systems usually divides traffic state into three states: smooth, congested, and blocked. It has oversimplified the problem since traffic state parameter is a continuous variable, e.g., flow, density, and speed. This study presents a new concept to describe traffic state. It is the traffic state reliability which is defined as the degree of actual traffic flow relative to the free flow, which assumes that the reliability of free flow is highest. In order to get the traffic flow trend in the future, the forecasting of traffic state reliability is needed to quantitatively measure the reliability of each alternative in future.For example, there are five links A, B, C, D, and E for the travelers. The prediction result shows that links A, C, and D are in smooth state. However, the travelers do not know how to determine which is the best link among links A, C, and D. If the reliability prediction could provide the level of actual traffic flow relative to the most reliable free flow state of each link (that is, 0.9, 0.2, 0.92, 0.79, and 0.1). Obviously, the traveler will prioritize link C with the highest reliability value to travel.Antoniou et al. [1] present an approach for local traffic state estimation and traffic state prediction; the method exploits all available (traffic and other) information and uses data-driven computational techniques. The approach is advantageous as it can flexibly incorporate additional explanatory variables. Given that previously proposed models can outperform current state-of-the-art models, integrating them into existing traffic estimation and prediction models is valuable. Current technologies, such as global position system- (GPS-) enabled cell phones, can record vehicle trajectories and have opened a new way of collecting traffic data. Hiribarren and Herrera [2] present and assess a new method to estimate traffic states on arterials based on trajectory data. The method is based on the Lighthill–Whitham–Richards theory. Preliminary analysis based on microsimulation suggested that this method yields good traffic state estimates at congested and uncongested situations. Mannini et al. [3] explore traffic state estimation on freeways in urban areas and route-based data to properly feed a second-order traffic flow model. This model is recursively corrected by an extended Kalman filter. Considering the lack of real-time information, these authors use simulation-based data to improve the traffic state estimation accuracy.In recent years, traffic data are currently collected through various sensors, including loop detectors, probe vehicles, cell phones, Bluetooth devices, video cameras, remote sensing applications, and public transport smart cards. Nantes et al. [4] develop a new model-based methodology for real-time traffic state estimation in urban corridors from multiple sources, particularly from loop detectors and partial observations from Bluetooth and GPS devices. Zheng and Su [5] develop a novel algorithm based on compressed sensing theory to recover traffic data with Gaussian measurement noise, missing partial data, and corrupted noise using information recovered from noisy traffic data and traffic state estimation. These authors extend traffic state estimation method to handle traffic state variables of high dimensions.Traffic state estimation is a key problem with considerable implication for modern traffic management. Kong et al. [6] propose a novel approach to efficiently estimate and predict urban traffic congestion using floating car trajectory data. An innovative traffic flow prediction method using particle swarm optimization is used to calculate the traffic flow parameters for predicting traffic congestion. Then, a congestion state fuzzy division module is applied to convert the predicted flow parameters to the cognitive congestion state of citizens.Travelers should be able to select the best path to avoid traffic congestion. However, if travelers also avoid the link with high variability, then travelers will enjoy additional benefits, that is, the “reliability benefits.” To date, traffic reliability has introduced the idea of reliability into traffic research and is an important field of traffic problems analysis. Considerable research has been conducted on traffic reliability, covering from theory to practice and from model to algorithm. Frameworks for reliability analysis have also been developed.Wang et al. [7] comprehensively review the literature on traffic reliability. Basic definitions, theory, and methods are depicted accompanied with the application of traffic optimization. These authors also discuss the future development of traffic reliability research. Xiao et al. [8, 9] consider the influence of travel time variability on the congestion profile in the presence of endogenous congestion and the role of scheduling preferences. The results show that the cost of travel time variability is the same in exogenous or endogenous congestion for two classes of preferences: linear at work and constant-exponential marginal utility of time. He et al. [10] introduce road segment and network congestion indexes to, respectively, measure the congestion levels of urban road segment and network. He also carries out a traffic congestion analysis for Beijing expressway network, based on the speed performance data.Probabilistic forecasting of reliability can be used for risk-averse routing. Bezuglov and Comert [11] have studied the possible applications and accuracy levels of three gray system theory models for short-term traffic speed and travel time predictions. Gray models consistently demonstrate low prediction errors over all time series, thereby successfully improving the accuracy by approximately 50% on average in root-mean-squared errors and mean absolute percent errors. Traffic parameters can vary due to several factors, such as weather, accidents, and driving characteristics. Comert et al. [12] develop a model to predict traffic speed under abrupt change. The developed model is tested from 1-step to 45-step forecasts. The accuracy of predictions is improved until the 15-step forecast compared with nonadaptive and mean adaptive models. Although the developed model is not retrained on different data sets, the method provides better results than or close results to nonadaptive and adaptive models that are retrained on the corresponding data set. Chen and Rakha [13] develop an agent-based modeling approach to predict multistep ahead experienced travel time using real-time and historical spatiotemporal traffic data. The results show that the agent-based modeling approach produced the least prediction errors compared with other state-of-the-practice and state-of-the-art methods (such as instantaneous travel time, historical average, and k-nearest neighbor). The fast algorithm computation allows the proposed approach to be implemented in real-time applications in traffic management centers.Papathanasopoulou et al. [14] develop a methodology to realize online calibration of microscopic traffic simulation models for dynamic multistep prediction of traffic measures. The application leads to less than 10% error in speed prediction even for 10 steps into the future in all considered data sets. Rajabzadeh et al. [15] provide a two-step approach based on stochastic differential equations to improve short-term prediction. At the first step, a Hull–White model is applied to obtain a baseline prediction model from previous days. Then, the extended Vasicek model is employed to model the difference between observations and baseline predictions (residuals) during an individual day. The results show that the proposed model can accurately input the missing data in traffic data set.In summary, most studies focus on traffic state estimation and travel time reliability, which mainly emphasize the traffic quality but not the reliability of traffic quality. Advanced Traveler Information System is one of the functional areas of Intelligent Transportation Systems and aims to provide forecasted traffic information for travelers to make better decisions. There are two methods to forecast the traffic data. One is one-step prediction, which predicts the reliability in the next one step. The other is fixed multistep prediction, in which the number of prediction steps is fixed. The results of fixed multistep prediction show that it easily leads to large prediction errors and the prediction errors gradually linearly increase with the number of multisteps [14]. For example, fixed multistep prediction can obtain the reliability in the next 10 days or even 30 days; however, the prediction error may be 10% in the next 10 days and 40% in the next 30 days. If an acceptable prediction errors is less than 20% for travelers (it can be obtained through investigation), do you want an accurate 10-day forecast or an inaccurate 30-day forecast? How to determine the number of predictable steps, which can not only guarantee exceeding the acceptable prediction errors but also provide as much future reliability data as possible. If there is a certain numberS of steps to guarantee it, the numberS of predictable steps is defined as the dynamic predictability.The proposed concept of dynamic predictability needs to find the numberS of predictable steps. After determining the number of predictable steps, the multistep prediction limited to the number of predictable steps will be performed. The realization of this theoretical principle is described in the following. Firstly, the traffic state reliability is defined, and the multistep forecasting theory and basic conditions are put forward in Section 2. Section 3 describes a two-phase model which determines the number of predictable steps in the upper-phase model and performs multistep prediction in the lower-phase model. Based on this, Section 4 verifies the validity and accuracy of the proposed method by applying actual traffic data. Section 5 summarizes the research conclusions.


## 2. Methodology

The definition of traffic state reliability (TSR) is the basis for setting a method to measure reliable level of traffic state. What is the relationship between past and present reliability and the future reliability? It is the discussion focus of this part.


### 2.1. Traffic State Reliability Calculation

For a link, intersection, or road network, TSR is defined as the degree of actual traffic flow state relative to the most reliable free flow state, which assumes that the reliability of free flow state is highest. It uses the TSR index to quantify evaluation, which is the ratio between actual speed and free flow speed in this research. TRS is classified to be equal to 1, when the actual speed is larger than or equal to the free flow speed and equal toypt (see (1)) when the actual speed is less than the free flow speed. Therefore, TSR value is distributed from 0 to 1 and is calculated as (1)ypt=vptvfpwhere ypt is the TSR of link or direction linkp at the periodt; vpt is the average travel speed of link or direction linkp at the periodt; vfp is the free flow speed of link or direction linkp.


### 2.2. Multistep Prediction Principle

One-step prediction [11] means predicting the reliability index in the next one step based on current and past reliability index series. Multistep prediction [12] means predicting the reliability index in the next several steps. The reliability series yt-n+1,⋯yt-1,yt of Monday consists of TSR of the continuously (t-n+1)th Monday,  ⋯, (t-1)th Monday,  ⋯,tth Monday, also for all the other days of the week. The prediction principle of futureSth series of reliability is(2)y∧t+s=fsyt-n+1,⋯yt-1,ytwhere 1≤s≤S, 0<n≤t, andt is the current time interval; y^t+s represents the predictive value of reliability in the next(t+s)th interval; fs represents the predictive function;n is the number of time intervals.The above method can be used to offline predict the reliability of the next S steps by a kind of prediction method. The prediction results show that the average prediction errors gradually linearly increase with the number of steps. The graph shows the relationship between the average prediction errors and the number S of steps is seen in Figure1 (Data Source in Section 4.1). The acceptable prediction errors investigated are less than 20% for travelers. Ultimately, the error 20% is the error threshold; the numberS of steps corresponding to error 20% is the dynamic predictability. The prediction error is related to the characteristic of historical traffic flow data used. Traffic flow is characterized by randomness, regularity, and volatility according to traffic flow theory [8, 9]. Therefore, three key indexes should be first calculated to measure the characteristics of reliability series. Then, a model should be constructed to associate the key index with the dynamic predictability.Figure 1
The relationship between the average prediction errors and the number S of steps.Fluctuation index is used to measure the random fluctuations of TSR series data. It refers to the ratio of variance to average reliability series. The lessBpt, the better the dynamic predictability. This index is expressed as(3)Bpt=A·∑i=1nyip2-yp¯2/n-1yp¯,wheren is the number of time intervals of TSR series. Bpt is the volatility index of link or direction linkp at periodt. y¯p is the average value relative to n time interval. yip is the TSR value of link or direction linkp at theith time intervalA which is the adjustment coefficient.Tendency index is used to quantitatively measure the tendency characteristic of the continuous increasing or decreasing of TSR series data. It refers to the cumulative amplitude change between two consecutive TSR of data series. The lessQpt, the better the dynamic predictability.This index is expressed as(4)Qpt=∑i=1nyip-ytp/n-kB·kwhere Qpt is the tendency index of link or direction linkp at periodt.k is the gradient calculated based on the slope of TSR series.ytp is the reliability of linkp at the current interval.B is the scaling factor.Uncertainty characteristic index is used to quantitatively measure the disorder characteristic of the TSR series data. This index presents an inverse proportional relationship with the number of predictable steps. It can be quantitatively measured by a fuzzy entropy indicator based on fuzzy analysis theory. The lessLpt, the better the dynamic predictability. It is expressed as(5)Lpt=C·log2n+∑i=t-n+1tMiplog2Mip,whereLpt is the fuzzy entropy.C is the scaling factor. Mip is the ambiguity and is expressed as(6)Mip=yip∑i=t-n+1tyipBased on above description, it is required to findS first by using the dynamic predictability model to guarantee that the prediction error will not exceed the error threshold. Then, prediction within the maximum numberS of steps will be conducted through the multistep forecasting model. Therefore, an upper-phase model is introduced to obtain the numberS of predictable steps within the error threshold. After that, the multistep prediction limited to the number of predictable steps will be performed by the lower-phase model.


## 2.1. Traffic State Reliability Calculation

For a link, intersection, or road network, TSR is defined as the degree of actual traffic flow state relative to the most reliable free flow state, which assumes that the reliability of free flow state is highest. It uses the TSR index to quantify evaluation, which is the ratio between actual speed and free flow speed in this research. TRS is classified to be equal to 1, when the actual speed is larger than or equal to the free flow speed and equal toypt (see (1)) when the actual speed is less than the free flow speed. Therefore, TSR value is distributed from 0 to 1 and is calculated as (1)ypt=vptvfpwhere ypt is the TSR of link or direction linkp at the periodt; vpt is the average travel speed of link or direction linkp at the periodt; vfp is the free flow speed of link or direction linkp.


## 2.2. Multistep Prediction Principle

One-step prediction [11] means predicting the reliability index in the next one step based on current and past reliability index series. Multistep prediction [12] means predicting the reliability index in the next several steps. The reliability series yt-n+1,⋯yt-1,yt of Monday consists of TSR of the continuously (t-n+1)th Monday,  ⋯, (t-1)th Monday,  ⋯,tth Monday, also for all the other days of the week. The prediction principle of futureSth series of reliability is(2)y∧t+s=fsyt-n+1,⋯yt-1,ytwhere 1≤s≤S, 0<n≤t, andt is the current time interval; y^t+s represents the predictive value of reliability in the next(t+s)th interval; fs represents the predictive function;n is the number of time intervals.The above method can be used to offline predict the reliability of the next S steps by a kind of prediction method. The prediction results show that the average prediction errors gradually linearly increase with the number of steps. The graph shows the relationship between the average prediction errors and the number S of steps is seen in Figure1 (Data Source in Section 4.1). The acceptable prediction errors investigated are less than 20% for travelers. Ultimately, the error 20% is the error threshold; the numberS of steps corresponding to error 20% is the dynamic predictability. The prediction error is related to the characteristic of historical traffic flow data used. Traffic flow is characterized by randomness, regularity, and volatility according to traffic flow theory [8, 9]. Therefore, three key indexes should be first calculated to measure the characteristics of reliability series. Then, a model should be constructed to associate the key index with the dynamic predictability.Figure 1
The relationship between the average prediction errors and the number S of steps.Fluctuation index is used to measure the random fluctuations of TSR series data. It refers to the ratio of variance to average reliability series. The lessBpt, the better the dynamic predictability. This index is expressed as(3)Bpt=A·∑i=1nyip2-yp¯2/n-1yp¯,wheren is the number of time intervals of TSR series. Bpt is the volatility index of link or direction linkp at periodt. y¯p is the average value relative to n time interval. yip is the TSR value of link or direction linkp at theith time intervalA which is the adjustment coefficient.Tendency index is used to quantitatively measure the tendency characteristic of the continuous increasing or decreasing of TSR series data. It refers to the cumulative amplitude change between two consecutive TSR of data series. The lessQpt, the better the dynamic predictability.This index is expressed as(4)Qpt=∑i=1nyip-ytp/n-kB·kwhere Qpt is the tendency index of link or direction linkp at periodt.k is the gradient calculated based on the slope of TSR series.ytp is the reliability of linkp at the current interval.B is the scaling factor.Uncertainty characteristic index is used to quantitatively measure the disorder characteristic of the TSR series data. This index presents an inverse proportional relationship with the number of predictable steps. It can be quantitatively measured by a fuzzy entropy indicator based on fuzzy analysis theory. The lessLpt, the better the dynamic predictability. It is expressed as(5)Lpt=C·log2n+∑i=t-n+1tMiplog2Mip,whereLpt is the fuzzy entropy.C is the scaling factor. Mip is the ambiguity and is expressed as(6)Mip=yip∑i=t-n+1tyipBased on above description, it is required to findS first by using the dynamic predictability model to guarantee that the prediction error will not exceed the error threshold. Then, prediction within the maximum numberS of steps will be conducted through the multistep forecasting model. Therefore, an upper-phase model is introduced to obtain the numberS of predictable steps within the error threshold. After that, the multistep prediction limited to the number of predictable steps will be performed by the lower-phase model.


## 3. Study Model

Both the dynamic predictability and multistep prediction are nonlinear fitting problems. It is exactly the problem that wavelet neural network can solve. Wavelet neural network has the good capability of localization and nonlinear mapping [16]. Moreover, particle swarm optimization can further improve the convergence rate and parameter iteration process. Therefore, a two-phase model of multistep prediction is designed based on wavelet neural network optimized by particle swarm algorithm. The upper-phase model is the dynamic predictability estimation method and the lower-phase model is the multistep prediction method.


### 3.1. Two-Phase Model of Multistep Prediction Based on Wavelet Neural Network Optimized by Particle Swarm Algorithm

A three-layer wavelet neural network is typically used to fit nonlinear correlations. This study uses a three-layer wavelet neural network to design the two-phase model for multistep prediction. First, in the upper-phase model, a three-layer wavelet neural network (Figure1) is designed to estimate dynamic predictability. The inputs are the three key indexes of reliability series and the output is the numberS of predictable steps. The numberH of neurons in hidden layer is, respectively, set from 5 to 20 to achieve the fastest convergence rate. The nonlinear fitting is realized by using the Morlet function as shown in (7)S∧=∑h=1Hωhψ∑i=13ωhixi-bhahwhere ψ is the Morlet wavelet function; ωh is the weight between the output layer and thehth wavelet in the hidden layer; ωhi is the weight between thehth wavelet in the hidden layer and theith input; bh and ah are the translation factor and stretch factor of thehth wavelet in the hidden layer, respectively. The Morlet wavelet function is shown in (8)ψx=cos⁡1.75xexp⁡-x22Second, in the lower-phase model, a three-layer wavelet neural network (Figure2) is designed to carry out S steps prediction. Its input is reliability series data, and the output is S steps reliability in the future. The number of neurons in hidden layer is, respectively, set from 5 to 20 to achieve the fastest convergence rate. The hidden layer also uses the Morlet function as the excitation function, as shown in (9)yt+s∧=∑f=1Fωfsψ∑n=1Nωfnxn-bfafwhere ψ is the Morlet wavelet function; ωfs is the weight between thesth output and thefth wavelet in the hidden layer; ωfn is the weight between thefth wavelet in the hidden layer and thenth input; bf and af are the translation factor and stretch factor of thefth wavelet in the hidden layer, respectively.Figure 2
Structure of the wavelet neural network for upper phase.The iterative optimization algorithm of wavelet neural network is algorithm with descent gradient. Notably, there is a local minimization problem. Particle swarm optimization can extract the transient characteristics of target as input, which has a better convergence and prediction accuracy. The principle of particle swarm optimization is to randomly generate a certain number of particle swarm optimizations, which has three attribute indexes: position, velocity, and fitness. The position represents a possible solution to the optimization problem. The velocity is a vector, which determines the direction and size of iteration. The fitness is used to measure the good or poor position. If evaluation result is poor, then the position should be updated through iteration. Supposing the position vector of themth particle is Xm=xm,1,⋯,xm,d,⋯,xm,D, then its velocity vector is Vm=vm,1,⋯,vm,d,⋯,vm,D, where D is the number of dimensions. During every iterative process, themth particle first searches for its own best solution pbestm=pm,1,pm,2,⋯,pm,D, that is, the individual extreme value. Then, all particles find the group best solutiongbest=pg,1,pg,2,⋯,pg,D, that is, the global extreme value. When the fitness corresponding to the global extreme value satisfies the error requirement, the iteration reaches convergence and the global extreme value is obtained as the optimal solution. Otherwise, the velocity and position are iteratively updated according to the method shown in (10) and (11). The movement of particles can be appropriately restricted by setting the velocity range vmin,vmax and the position range xmin,xmax to ensure that the particles can efficiently search the global extreme value and avoid the local minimum problem.(10)vm,di+1=ω×vm,di+c1r1pm,d-xm,di+c2r2pg,d-xm,di;(11)xm,di+1=xm,di+vm,di+1,d=1,2,⋯,D,wherei is the current time of iterations. ω is the inertia factor. r1 and r2 are the random values evenly distributed between 0 and 1. c1 and c2 are the acceleration factors that are determined using the most commonly used penalty function strategy as shown in(12)c1=c1s+c1e-c1sTmax×Iter(13)c2=c2s+c2e-c2sTmax×IterThe inertia factorω controls the influences of iteration speed of previous step on the current iteration speed. Small ω is beneficial to convergence, whereas large ω is conducive to jump out of the local minimum for particle. Considering the global searching and local searching, the inertia factor ω is iterated by the linear decreasing method as shown in(14)ωIter=ωmax-ωmax-ωminTmax×Iter,whereIter is the current number of iterations, ωmax=0.95, ωmin=0.25, Tmax is the total times of iterations, C1s=2.75, C1e=1.25, C2s=0.25, andC2e=2.25.


### 3.2. Algorithm of Two-Phase Model for Multistep Prediction

For a particular place, the current reliability series isyt-n+1,⋯yt-1,yt at the timet and three key indexes are Bt,Qt,Lt.G data samples of historical series Y1,⋯Yg,⋯YG are available in database, where Yg=ygt-n+1,⋯ygt-1,ygt.Step 1.
The three key indexesT1,⋯Tg,⋯TG of data samples are calculated, where Tg=Bgt,Qgt,Lgt, and the fixedX steps (X ≫S) prediction of data samples is conducted by applying the wavelet neural network optimized by particle swarm algorithm. Then, the actual number S1,⋯Sg,⋯SG of predictable steps is determined by combining with an acceptable error threshold, in which the prediction error is continuously less than the threshold. Among the data samples in the database, two-thirds of the sample data (totalI) is used as the training data set, and the rest is used as the testing data set.Step 2.
The upper-phase model is constructed as shown in Figure2. The numberH of neurons in hidden layer is set from 5 to 20 in training the network for selecting the optimal number of neurons.Step 3.
The training parameter is optimized by particle swarm algorithm.M particles are randomly generated. The position vector of each particle is represented as the real number coded that includes the stretch factor ah, translation factor bh of initialized wavelet function, the connection weight ωih between theith input and thehth hidden layer, and the connection weight ωh between thehth hidden layer and output layers. The function is expressed as (15)Xm=xm,1,⋯,xm,d,⋯,xm,D=a1,⋯,aH,b1,⋯bH,ω11,⋯ω1H,ω21,⋯ω2H,ω31,⋯ω3H,ω1,⋯,ωHwhere D=3×H+H×1+2×H=6H.Step 4.
The upper-phase model is trained by using the data sample from the training data set. The three key indexesBit,Qit,Lit of data sample are the input and the numberS of predictable steps is the desired output. TheSAPE (average prediction error of steps) between the predicted and desired output are used as the fitness value of particle as shown in (16)SAPE=∑i=1ISi-Si∧I×Si×100%.where SAPE is the relative average prediction error of steps. Si is the actual number of predictable steps of theith data sample. S^i is the predicted number of predictable steps of theith data sample.Step 5.
The upper-phase model is tested using the data sample from the testing data set until the average relative error satisfies the requirements that means the upper-phase model fits well. Then, the predictable step numberS of the reliability series at the current periodt is forecasted through input Bt,Qt,Lt. Otherwise, Step 3 is repeated.Step 6.
After obtaining the numberSof predictable steps, a multistep prediction model is constructed to forecast the reliability in the subsequentS step based on wavelet neural network optimized by particle swarm algorithm. The actual value Y1′,⋯Yg′,⋯YG′ of reliability data sample in the futureS steps is extracted from the database, where Yg′=ygt+1,⋯ygt+S-1,ygt+S. Two-thirds of the sample data (totalI) are used as the training data set, and the rest are used as the testing data set.Step 7.
The lower-phase model is constructed as shown in Figure3. The numberF of neurons in hidden layer is set from 5 to 20 in training the network for selecting the optimal number of neurons.Figure 3
Structure of the wavelet neural network for lower phase.Step 8.
The training parameter is optimized using particle swarm algorithm andM particles are randomly generated. The position vector of each particle is represented by the real number coded that includes the stretch factor af, translation factor bf of initialized wavelet function, the connection weight ωfn between thenth input and thefth hidden layer, and the connection weight ωfs between thefth hidden layer and thesth output layer. The function is expressed as(17)Xm=xm,1,⋯,xm,d,⋯,xm,D=a1,⋯,aF,b1,⋯bF,ω11r,⋯ω1Fr,⋯,ωNFr,ω11c,⋯,ωFScwhere D=2×F+N×F+F×S=2+N+SF.Step 9.
The lower-phase model is trained by using the data sample from the training data set. The reliability series dataygt-n+1,⋯ygt-1,ygt are the input and the actual value ygt+1,⋯ygt+S-1,ygt+S in the futureS steps is the desired output. TheMAPE (average predictive error of multistep)between the predicted and desired output is used as the fitness value of particle as shown in(18)MAPE=1I×S∑s=1S∑i=1Iyis-yis∧yis×100%,where MAPE is the average predictive error of multistep prediction. yisis the actual reliability of theith data sample in the futuret+s time interval. y^is is the predicted reliability of theith data sample in the futuret+stime interval.I is the total number of data samples in training data set.S is the number of predictable steps.Step 10.
The lower-phase model is tested using the data sample from the testing data set. If the average relative errorMAPE satisfies the requirements, then the training is terminated and the reliability in futureS steps of reliability series at the current periodtis forecasted. Otherwise, Step 8 is repeated.In order to verify the effectiveness of the above model and algorithm, it is necessary to collect and analyze the actual traffic flow parameters, calculate proposed TSR, estimate the dynamic predictability, and predict the multistep reliability in future. Eventually, it is necessary to compare the research result with actual traffic state to confirm the validity of the proposed method.


## 3.1. Two-Phase Model of Multistep Prediction Based on Wavelet Neural Network Optimized by Particle Swarm Algorithm

A three-layer wavelet neural network is typically used to fit nonlinear correlations. This study uses a three-layer wavelet neural network to design the two-phase model for multistep prediction. First, in the upper-phase model, a three-layer wavelet neural network (Figure1) is designed to estimate dynamic predictability. The inputs are the three key indexes of reliability series and the output is the numberS of predictable steps. The numberH of neurons in hidden layer is, respectively, set from 5 to 20 to achieve the fastest convergence rate. The nonlinear fitting is realized by using the Morlet function as shown in (7)S∧=∑h=1Hωhψ∑i=13ωhixi-bhahwhere ψ is the Morlet wavelet function; ωh is the weight between the output layer and thehth wavelet in the hidden layer; ωhi is the weight between thehth wavelet in the hidden layer and theith input; bh and ah are the translation factor and stretch factor of thehth wavelet in the hidden layer, respectively. The Morlet wavelet function is shown in (8)ψx=cos⁡1.75xexp⁡-x22Second, in the lower-phase model, a three-layer wavelet neural network (Figure2) is designed to carry out S steps prediction. Its input is reliability series data, and the output is S steps reliability in the future. The number of neurons in hidden layer is, respectively, set from 5 to 20 to achieve the fastest convergence rate. The hidden layer also uses the Morlet function as the excitation function, as shown in (9)yt+s∧=∑f=1Fωfsψ∑n=1Nωfnxn-bfafwhere ψ is the Morlet wavelet function; ωfs is the weight between thesth output and thefth wavelet in the hidden layer; ωfn is the weight between thefth wavelet in the hidden layer and thenth input; bf and af are the translation factor and stretch factor of thefth wavelet in the hidden layer, respectively.Figure 2
Structure of the wavelet neural network for upper phase.The iterative optimization algorithm of wavelet neural network is algorithm with descent gradient. Notably, there is a local minimization problem. Particle swarm optimization can extract the transient characteristics of target as input, which has a better convergence and prediction accuracy. The principle of particle swarm optimization is to randomly generate a certain number of particle swarm optimizations, which has three attribute indexes: position, velocity, and fitness. The position represents a possible solution to the optimization problem. The velocity is a vector, which determines the direction and size of iteration. The fitness is used to measure the good or poor position. If evaluation result is poor, then the position should be updated through iteration. Supposing the position vector of themth particle is Xm=xm,1,⋯,xm,d,⋯,xm,D, then its velocity vector is Vm=vm,1,⋯,vm,d,⋯,vm,D, where D is the number of dimensions. During every iterative process, themth particle first searches for its own best solution pbestm=pm,1,pm,2,⋯,pm,D, that is, the individual extreme value. Then, all particles find the group best solutiongbest=pg,1,pg,2,⋯,pg,D, that is, the global extreme value. When the fitness corresponding to the global extreme value satisfies the error requirement, the iteration reaches convergence and the global extreme value is obtained as the optimal solution. Otherwise, the velocity and position are iteratively updated according to the method shown in (10) and (11). The movement of particles can be appropriately restricted by setting the velocity range vmin,vmax and the position range xmin,xmax to ensure that the particles can efficiently search the global extreme value and avoid the local minimum problem.(10)vm,di+1=ω×vm,di+c1r1pm,d-xm,di+c2r2pg,d-xm,di;(11)xm,di+1=xm,di+vm,di+1,d=1,2,⋯,D,wherei is the current time of iterations. ω is the inertia factor. r1 and r2 are the random values evenly distributed between 0 and 1. c1 and c2 are the acceleration factors that are determined using the most commonly used penalty function strategy as shown in(12)c1=c1s+c1e-c1sTmax×Iter(13)c2=c2s+c2e-c2sTmax×IterThe inertia factorω controls the influences of iteration speed of previous step on the current iteration speed. Small ω is beneficial to convergence, whereas large ω is conducive to jump out of the local minimum for particle. Considering the global searching and local searching, the inertia factor ω is iterated by the linear decreasing method as shown in(14)ωIter=ωmax-ωmax-ωminTmax×Iter,whereIter is the current number of iterations, ωmax=0.95, ωmin=0.25, Tmax is the total times of iterations, C1s=2.75, C1e=1.25, C2s=0.25, andC2e=2.25.


## 3.2. Algorithm of Two-Phase Model for Multistep Prediction

For a particular place, the current reliability series isyt-n+1,⋯yt-1,yt at the timet and three key indexes are Bt,Qt,Lt.G data samples of historical series Y1,⋯Yg,⋯YG are available in database, where Yg=ygt-n+1,⋯ygt-1,ygt.Step 1.
The three key indexesT1,⋯Tg,⋯TG of data samples are calculated, where Tg=Bgt,Qgt,Lgt, and the fixedX steps (X ≫S) prediction of data samples is conducted by applying the wavelet neural network optimized by particle swarm algorithm. Then, the actual number S1,⋯Sg,⋯SG of predictable steps is determined by combining with an acceptable error threshold, in which the prediction error is continuously less than the threshold. Among the data samples in the database, two-thirds of the sample data (totalI) is used as the training data set, and the rest is used as the testing data set.Step 2.
The upper-phase model is constructed as shown in Figure2. The numberH of neurons in hidden layer is set from 5 to 20 in training the network for selecting the optimal number of neurons.Step 3.
The training parameter is optimized by particle swarm algorithm.M particles are randomly generated. The position vector of each particle is represented as the real number coded that includes the stretch factor ah, translation factor bh of initialized wavelet function, the connection weight ωih between theith input and thehth hidden layer, and the connection weight ωh between thehth hidden layer and output layers. The function is expressed as (15)Xm=xm,1,⋯,xm,d,⋯,xm,D=a1,⋯,aH,b1,⋯bH,ω11,⋯ω1H,ω21,⋯ω2H,ω31,⋯ω3H,ω1,⋯,ωHwhere D=3×H+H×1+2×H=6H.Step 4.
The upper-phase model is trained by using the data sample from the training data set. The three key indexesBit,Qit,Lit of data sample are the input and the numberS of predictable steps is the desired output. TheSAPE (average prediction error of steps) between the predicted and desired output are used as the fitness value of particle as shown in (16)SAPE=∑i=1ISi-Si∧I×Si×100%.where SAPE is the relative average prediction error of steps. Si is the actual number of predictable steps of theith data sample. S^i is the predicted number of predictable steps of theith data sample.Step 5.
The upper-phase model is tested using the data sample from the testing data set until the average relative error satisfies the requirements that means the upper-phase model fits well. Then, the predictable step numberS of the reliability series at the current periodt is forecasted through input Bt,Qt,Lt. Otherwise, Step 3 is repeated.Step 6.
After obtaining the numberSof predictable steps, a multistep prediction model is constructed to forecast the reliability in the subsequentS step based on wavelet neural network optimized by particle swarm algorithm. The actual value Y1′,⋯Yg′,⋯YG′ of reliability data sample in the futureS steps is extracted from the database, where Yg′=ygt+1,⋯ygt+S-1,ygt+S. Two-thirds of the sample data (totalI) are used as the training data set, and the rest are used as the testing data set.Step 7.
The lower-phase model is constructed as shown in Figure3. The numberF of neurons in hidden layer is set from 5 to 20 in training the network for selecting the optimal number of neurons.Figure 3
Structure of the wavelet neural network for lower phase.Step 8.
The training parameter is optimized using particle swarm algorithm andM particles are randomly generated. The position vector of each particle is represented by the real number coded that includes the stretch factor af, translation factor bf of initialized wavelet function, the connection weight ωfn between thenth input and thefth hidden layer, and the connection weight ωfs between thefth hidden layer and thesth output layer. The function is expressed as(17)Xm=xm,1,⋯,xm,d,⋯,xm,D=a1,⋯,aF,b1,⋯bF,ω11r,⋯ω1Fr,⋯,ωNFr,ω11c,⋯,ωFScwhere D=2×F+N×F+F×S=2+N+SF.Step 9.
The lower-phase model is trained by using the data sample from the training data set. The reliability series dataygt-n+1,⋯ygt-1,ygt are the input and the actual value ygt+1,⋯ygt+S-1,ygt+S in the futureS steps is the desired output. TheMAPE (average predictive error of multistep)between the predicted and desired output is used as the fitness value of particle as shown in(18)MAPE=1I×S∑s=1S∑i=1Iyis-yis∧yis×100%,where MAPE is the average predictive error of multistep prediction. yisis the actual reliability of theith data sample in the futuret+s time interval. y^is is the predicted reliability of theith data sample in the futuret+stime interval.I is the total number of data samples in training data set.S is the number of predictable steps.Step 10.
The lower-phase model is tested using the data sample from the testing data set. If the average relative errorMAPE satisfies the requirements, then the training is terminated and the reliability in futureS steps of reliability series at the current periodtis forecasted. Otherwise, Step 8 is repeated.In order to verify the effectiveness of the above model and algorithm, it is necessary to collect and analyze the actual traffic flow parameters, calculate proposed TSR, estimate the dynamic predictability, and predict the multistep reliability in future. Eventually, it is necessary to compare the research result with actual traffic state to confirm the validity of the proposed method.


## 4. Case Study

The empirical data were obtained from a freeway of Shanghai. The structure of road network is shown in Figure4. The process of empirical analysis is divided into three steps: data source, experimental program, experimental results, and comparative analysis.Figure 4
Diagram of the road experimental area.


### 4.1. Data Source

Traffic data was collected by the loop detector from a freeway in Shanghai of China. The collected traffic parameters mainly include speed, traffic volume, and occupancy. The size of data set is 24 hours of every day from May to July. The time interval of collected data is 1 minute.


### 4.2. Experimental Program

The reliability of each link was calculated using (1). The data collected in May and June is used as training data set, and the data collected in July is used as testing data set. Parameters A, B, and C have been set according to analysis results of sensitivity analysis method through assigning 5 to 200, respectively. The results show that A is set as 50, B is set as 5, and C is set as 150.First, the step number of iterations is used to measure convergence effect of wavelet neural network optimized by particle swarm algorithm. Second, the mean relative error SAPE is calculated to measure performance of the upper-phase model. Third, the mean relative error MAPE is calculated to measure performance of the lower-phase model. As traffic state on weekday is usually very different from that on weekend, the number of neurons in the middle layer is separately determined, and the forecasting work is also verified separately.Liu et al. [17] conducted a prediction model of time series based on backpropagation neural network. Bao et al. [18] performed a multistep-ahead time series prediction using multiple-output support vector regression. The proposed method is verified comparatively with backpropagation neural network and support vector regression.The difference between wavelet neural network and classical artificial neural network (such as backpropagation neural network) is excitation function in middle hidden layer. The excitation function of the latter usually uses sigmoid function, which leads to nonconvergence or slow convergence. However, the wavelet neural network adopts wavelet function as excitation function. It introduces translation scale and scaling factor to extract local information. Thus, the convergence speed, approximation precision, and performance are improved compared with those of classical artificial neural network.Support vector machine is one learning method for small sample. It ignores probability measure and law of large numbers unlike existing statistical methods. It is also based on the minimization principle of structural risk. The global optimal solution is obtained through the optimization problem of convex quadratic, which has high generalization capability and an advantage in terms of generalization and classification.


### 4.3. Experimental Results and Comparative Analysis

#### 4.3.1. Convergence Analysis of Wavelet Neural Network Optimized by Particle Swarm Algorithm

The convergence rate was trained and comparatively analyzed in upper-phase and lower-phase model using the same set of reliability data series and under the same conditions. The result is shown in Figures5 and 6.Figure 5
Convergence of upper-phase model based on different methods.Figure 6
Convergence of lower-phase model based on different methods.Figures5 and 6 show that upper-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 256 seconds, whereas backpropagation neural network and support vector regression achieve convergence effect at 398 and 463 seconds, respectively. The lower-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 291 seconds, whereas the other begins convergence effect at 487 and 526 seconds, respectively. Therefore, the parameter-optimizing process based on particle swarm algorithm contributes greatly to the fast convergence speed.


#### 4.3.2. Determining the Number of Neurons, Respectively, Based on Weekday and Weekend Data

In order to determine the number of neurons in hidden layers of upper-phase and lower-phase model, the number of neurons is set from 5 to 20 to train, respectively, based on weekday and weekend data. The result is shown in Figures7 and 8.Figure 7
Relationship between MAPE/SAPE and the number of neurons on weekday.Figure 8
Relationship between MAPE/SAPE and the number of neurons on weekend.Based on weekday data, Figure7 shows that SAPE was lowest when the number of neurons in hidden layer is 12. MAPE is lowest when the number of neurons in hidden layer is 14. Meanwhile, based on weekend data, Figure 8 shows that the SAPE is lowest when the number of neurons in hidden layer was 16. Meanwhile, the MAPE is lowest when the number of neurons in the hidden layer is 18.


#### 4.3.3. Performance of Upper-Phase Model, Respectively, Based on Weekday and Weekend Data

Based on weekday data, the number of predictable steps is estimated and analyzed comparatively based on wavelet neural network optimized by particle swarm algorithm, the backpropagation neural network, and the support vector regression. The result is shown in Table1. It shows that the proposed method presents a better performance compared with referenced methods. Meanwhile, based on weekend data, the number of predictable steps is estimated and analyzed comparatively. The result is shown in Table 2, which shows that the proposed method presents a better estimation performance compared with referenced methods.Table 1
SAPE of proposed method and referenced method for upper-phase model based on weekday data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 15.41 19.87 7.92 testing sample 20.75 20.74 8.63Table 2
SAPE of proposed method and referenced method for upper-phase model based on weekend data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 16.51 22.91 10.07 testing sample 24.73 20.62 12.73


### 4.4. Performance of Lower-Phase Model, Respectively, Based on Weekday and Weekend Data

Based on weekday data, the multistep prediction is analyzed comparatively based on wavelet neural network optimized by particle swarm algorithm, the backpropagation neural network, and the support vector regression. The result is shown in Table3. It shows that MAPE of lower-phase model is less than that of referenced methods. Meanwhile, based on weekend data, the multistep prediction result is shown in Table 4. It shows that the MAPE of lower-phase model is the lowest, 12.01%.Table 3
MAPE of proposed method and referenced method for lower-phase model based on weekday data (%).
Reliability index Data set support vector regression back propagation neural network Proposed method reliability series training sample 17.92 20.7 8.91 testing sample 20.15 21.58 10.02Table 4
MAPE of proposed method and referenced method for lower-phase model based on weekend data (%).
Reliability index Data set support vector regression back propagation neural network Proposed method TSR series training sample 22.32 23.89 12.01 testing sample 26.13 23.91 13.54Moreover, 30 steps in future are predicted, respectively, based on weekday and weekend data. The result is shown in Figures9 and 10.Figure 9
Prediction results of 30 steps for weekday reliability based on different methods.Figure 10
Prediction results of 30 steps for weekend reliability based on different methods.Figures9 and 10 indicate that prediction error increases as the increase of step number which presents a positive correlation. Meanwhile, MAPE of the proposed method is less than that of referenced methods. Moreover, the prediction performance based on weekday data is better than that of weekend.Meanwhile, 24-step reliability is predicted, respectively, on Monday and Saturday for the direction link 080810 which is from south to north in the intersection. The forecasting period is from 08:00:00 to 08:23:00, and the time interval is 1 min. The result is shown in Figures11 and 12, which shows that the reliability predicted by the proposed method is closer to actual values than that of backpropagation neural network and support vector regression. The proposed method obtains the lowest average relative error. Therefore, the proposed method shows a better prediction performance.Figure 11
Multistep prediction of weekday reliability based on different methods.Figure 12
Multistep prediction of weekend reliability based on different methods.


## 4.1. Data Source

Traffic data was collected by the loop detector from a freeway in Shanghai of China. The collected traffic parameters mainly include speed, traffic volume, and occupancy. The size of data set is 24 hours of every day from May to July. The time interval of collected data is 1 minute.


## 4.2. Experimental Program

The reliability of each link was calculated using (1). The data collected in May and June is used as training data set, and the data collected in July is used as testing data set. Parameters A, B, and C have been set according to analysis results of sensitivity analysis method through assigning 5 to 200, respectively. The results show that A is set as 50, B is set as 5, and C is set as 150.First, the step number of iterations is used to measure convergence effect of wavelet neural network optimized by particle swarm algorithm. Second, the mean relative error SAPE is calculated to measure performance of the upper-phase model. Third, the mean relative error MAPE is calculated to measure performance of the lower-phase model. As traffic state on weekday is usually very different from that on weekend, the number of neurons in the middle layer is separately determined, and the forecasting work is also verified separately.Liu et al. [17] conducted a prediction model of time series based on backpropagation neural network. Bao et al. [18] performed a multistep-ahead time series prediction using multiple-output support vector regression. The proposed method is verified comparatively with backpropagation neural network and support vector regression.The difference between wavelet neural network and classical artificial neural network (such as backpropagation neural network) is excitation function in middle hidden layer. The excitation function of the latter usually uses sigmoid function, which leads to nonconvergence or slow convergence. However, the wavelet neural network adopts wavelet function as excitation function. It introduces translation scale and scaling factor to extract local information. Thus, the convergence speed, approximation precision, and performance are improved compared with those of classical artificial neural network.Support vector machine is one learning method for small sample. It ignores probability measure and law of large numbers unlike existing statistical methods. It is also based on the minimization principle of structural risk. The global optimal solution is obtained through the optimization problem of convex quadratic, which has high generalization capability and an advantage in terms of generalization and classification.


## 4.3. Experimental Results and Comparative Analysis

### 4.3.1. Convergence Analysis of Wavelet Neural Network Optimized by Particle Swarm Algorithm

The convergence rate was trained and comparatively analyzed in upper-phase and lower-phase model using the same set of reliability data series and under the same conditions. The result is shown in Figures5 and 6.Figure 5
Convergence of upper-phase model based on different methods.Figure 6
Convergence of lower-phase model based on different methods.Figures5 and 6 show that upper-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 256 seconds, whereas backpropagation neural network and support vector regression achieve convergence effect at 398 and 463 seconds, respectively. The lower-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 291 seconds, whereas the other begins convergence effect at 487 and 526 seconds, respectively. Therefore, the parameter-optimizing process based on particle swarm algorithm contributes greatly to the fast convergence speed.


### 4.3.2. Determining the Number of Neurons, Respectively, Based on Weekday and Weekend Data

In order to determine the number of neurons in hidden layers of upper-phase and lower-phase model, the number of neurons is set from 5 to 20 to train, respectively, based on weekday and weekend data. The result is shown in Figures7 and 8.Figure 7
Relationship between MAPE/SAPE and the number of neurons on weekday.Figure 8
Relationship between MAPE/SAPE and the number of neurons on weekend.Based on weekday data, Figure7 shows that SAPE was lowest when the number of neurons in hidden layer is 12. MAPE is lowest when the number of neurons in hidden layer is 14. Meanwhile, based on weekend data, Figure 8 shows that the SAPE is lowest when the number of neurons in hidden layer was 16. Meanwhile, the MAPE is lowest when the number of neurons in the hidden layer is 18.


### 4.3.3. Performance of Upper-Phase Model, Respectively, Based on Weekday and Weekend Data

Based on weekday data, the number of predictable steps is estimated and analyzed comparatively based on wavelet neural network optimized by particle swarm algorithm, the backpropagation neural network, and the support vector regression. The result is shown in Table1. It shows that the proposed method presents a better performance compared with referenced methods. Meanwhile, based on weekend data, the number of predictable steps is estimated and analyzed comparatively. The result is shown in Table 2, which shows that the proposed method presents a better estimation performance compared with referenced methods.Table 1
SAPE of proposed method and referenced method for upper-phase model based on weekday data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 15.41 19.87 7.92 testing sample 20.75 20.74 8.63Table 2
SAPE of proposed method and referenced method for upper-phase model based on weekend data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 16.51 22.91 10.07 testing sample 24.73 20.62 12.73


## 4.3.1. Convergence Analysis of Wavelet Neural Network Optimized by Particle Swarm Algorithm

The convergence rate was trained and comparatively analyzed in upper-phase and lower-phase model using the same set of reliability data series and under the same conditions. The result is shown in Figures5 and 6.Figure 5
Convergence of upper-phase model based on different methods.Figure 6
Convergence of lower-phase model based on different methods.Figures5 and 6 show that upper-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 256 seconds, whereas backpropagation neural network and support vector regression achieve convergence effect at 398 and 463 seconds, respectively. The lower-phase model based on wavelet neural network optimized by particle swarm algorithm begins to converge at 291 seconds, whereas the other begins convergence effect at 487 and 526 seconds, respectively. Therefore, the parameter-optimizing process based on particle swarm algorithm contributes greatly to the fast convergence speed.


## 4.3.2. Determining the Number of Neurons, Respectively, Based on Weekday and Weekend Data

In order to determine the number of neurons in hidden layers of upper-phase and lower-phase model, the number of neurons is set from 5 to 20 to train, respectively, based on weekday and weekend data. The result is shown in Figures7 and 8.Figure 7
Relationship between MAPE/SAPE and the number of neurons on weekday.Figure 8
Relationship between MAPE/SAPE and the number of neurons on weekend.Based on weekday data, Figure7 shows that SAPE was lowest when the number of neurons in hidden layer is 12. MAPE is lowest when the number of neurons in hidden layer is 14. Meanwhile, based on weekend data, Figure 8 shows that the SAPE is lowest when the number of neurons in hidden layer was 16. Meanwhile, the MAPE is lowest when the number of neurons in the hidden layer is 18.


## 4.3.3. Performance of Upper-Phase Model, Respectively, Based on Weekday and Weekend Data

Based on weekday data, the number of predictable steps is estimated and analyzed comparatively based on wavelet neural network optimized by particle swarm algorithm, the backpropagation neural network, and the support vector regression. The result is shown in Table1. It shows that the proposed method presents a better performance compared with referenced methods. Meanwhile, based on weekend data, the number of predictable steps is estimated and analyzed comparatively. The result is shown in Table 2, which shows that the proposed method presents a better estimation performance compared with referenced methods.Table 1
SAPE of proposed method and referenced method for upper-phase model based on weekday data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 15.41 19.87 7.92 testing sample 20.75 20.74 8.63Table 2
SAPE of proposed method and referenced method for upper-phase model based on weekend data (%).
Reliability index Calibration data set support vector regression back propagation neural network Proposed method reliability series training sample 16.51 22.91 10.07 testing sample 24.73 20.62 12.73


## 4.4. Performance of Lower-Phase Model, Respectively, Based on Weekday and Weekend Data

Based on weekday data, the multistep prediction is analyzed comparatively based on wavelet neural network optimized by particle swarm algorithm, the backpropagation neural network, and the support vector regression. The result is shown in Table3. It shows that MAPE of lower-phase model is less than that of referenced methods. Meanwhile, based on weekend data, the multistep prediction result is shown in Table 4. It shows that the MAPE of lower-phase model is the lowest, 12.01%.Table 3
MAPE of proposed method and referenced method for lower-phase model based on weekday data (%).
Reliability index Data set support vector regression back propagation neural network Proposed method reliability series training sample 17.92 20.7 8.91 testing sample 20.15 21.58 10.02Table 4
MAPE of proposed method and referenced method for lower-phase model based on weekend data (%).
Reliability index Data set support vector regression back propagation neural network Proposed method TSR series training sample 22.32 23.89 12.01 testing sample 26.13 23.91 13.54Moreover, 30 steps in future are predicted, respectively, based on weekday and weekend data. The result is shown in Figures9 and 10.Figure 9
Prediction results of 30 steps for weekday reliability based on different methods.Figure 10
Prediction results of 30 steps for weekend reliability based on different methods.Figures9 and 10 indicate that prediction error increases as the increase of step number which presents a positive correlation. Meanwhile, MAPE of the proposed method is less than that of referenced methods. Moreover, the prediction performance based on weekday data is better than that of weekend.Meanwhile, 24-step reliability is predicted, respectively, on Monday and Saturday for the direction link 080810 which is from south to north in the intersection. The forecasting period is from 08:00:00 to 08:23:00, and the time interval is 1 min. The result is shown in Figures11 and 12, which shows that the reliability predicted by the proposed method is closer to actual values than that of backpropagation neural network and support vector regression. The proposed method obtains the lowest average relative error. Therefore, the proposed method shows a better prediction performance.Figure 11
Multistep prediction of weekday reliability based on different methods.Figure 12
Multistep prediction of weekend reliability based on different methods.


## 5. Conclusions

In view of the fact that the traffic state should be a continuous variable, this study proposes the concept of TSR to measure the stability of actual traffic flow state relative to the most reliable free flow state. In order to know TSR in future, a two-phase model of multistep prediction has been designed based on wavelet neural network optimized by particle swarm algorithm. In the upper-phase model, the numberS of predictable steps within the error threshold is determined according to three key indexes that are fluctuation index, tendency index, and uncertainty index. In the lower-phase model, a three-layer wavelet neural network has been designed to carry out prediction in S steps. In two-phase model, training parameter has been optimized by particle swarm algorithm, which has a better convergence and prediction accuracy. The proposed model can guarantee not only forecasting within the acceptable prediction errors but also providing as much future reliability data as possible. Moreover, it aggregates the prevailing nonlinear mapping capability of wavelet neural network and the excellent parameter optimization effect of particle swarm algorithm. The empirical analysis shows that parameter-optimizing process based on particle swarm algorithm has fast convergence speed. The proposed method presents a better prediction performance compared with backpropagation neural network and support vector regression. Moreover, the prediction performance of weekday is better than that of weekend.



---

*Source: 7650928-2018-07-09.xml*