Modeling forMicro Traffic Flowwith the Consideration of Lateral Vehicle’s Influence

In order to make the car-following model describe the driving behavior of vehicle on urban road more accurately, existing carfollowing models are simulated using measured traffic data. According to the analysis of the simulation result, two new improved car-following models based on the optimal velocity model (OVM) are proposed in this paper. 1e lateral vehicle’s influence is introduced as the influence factor of driving behavior. By using of linear stability analysis, stability conditions of improved carfollowing models are obtained. Nonlinear analysis is carried out to investigate the traffic performances near the critical point. 1e result of numerical simulation indicates that stability of traffic flow is under the influence from lateral vehicle; the lesser the influence, the greater the stability. New cooperative car-followingmodels are verified by the traffic flow data collected in Xi’an city. It is shown that compared with the optimal velocity model, the simulation result of the second cooperative model, respectively, gets 62.89% unbiased variance reduction, 66.39% maximum absolute error reduction, and 33.4% minimum absolute error reduction.1erefore, the second cooperative model is more suitable to describe the vehicle’s actual behavior in car-following state.


Introduction
Car-following theory is one of the most important parts of modern traffic theory, which describes the dynamic behaviors of one vehicle following the other, and usually used in the study of microscopic traffic flow. e car-following theory has developed rapidly in recent decades, and an increasing number of models have been put forward, which include the early nonlinear models proposed by Pipes [1], Newell [2], the recent classic model presented by Bando et al. [3], Helbing and Tilch [4] and Jiang et al. [5], and some others in the literature [6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24]. Bando et al. [3] proposed a remarkable car-following model called the optimal velocity model (OVM) because it solves the problem of describing the motion of successive vehicles. By considering both the negative and the positive velocity differences, Jiang et al. [5] developed a full velocity difference model (FVDM). Some scholars considered more headway information of vehicles ahead and proposed some extended optimal velocity models.
Based on OVM, Wagner [6] put forward a new model called the looking forward model because it contains the location information of multiple vehicles ahead. Hasebe et al. [7] further expanded the OVM by introducing the headway information of K vehicles ahead and proved that the new model can improve the stability of traffic flow more than the OVM from the perspective of the linear response. Chen et al. [8] extended the single lane feedback model by introducing improved lane-changing rules and analyzed the influence of feedback signal on the stability of traffic flow. Li et al. [9] proposed a heterogeneous car-following model with lowand high-sensitivity vehicles by constructing the motion equation of individual vehicle. In the proposed model, all driver-vehicle units are described by the same model using two different parameter values of sensitivity. Peng et al. [10] explored the effect of the V2V communication and proposed a delayed-feedback control method for the car-following model. In our prior work [11], we developed an extended full velocity difference model by taking the short-term driving memory into account. We found that short-term driving memory can significantly affect the driver's driving behavior.
Most of the above models are based on such basic assumptions that the driver always keeps following the leading vehicle and responds immediately to the change of the driving state of the leading vehicle. In recent years, scholars found that the driver is also influenced by the driving state of the lateral adjacent vehicles. Tang et al. [12,13] put forward a new two-lane car-following model by taking the influence of adjacent-ahead vehicle on the car-following behavior into account. Ge et al. [14] believed that vehicles will be subject to the friction interference from the vehicle of neighbor lane, and the friction interference is related to the lateral distance between the vehicle and its closest-ahead vehicle in the neighbor lane. Based on this consideration, they investigated an extended two-lane traffic flow with lane-changing behaviors and derived the stability condition by using the control method.
With the consideration of influence of lateral vehicle, an improved car-following model based on the OVM is proposed in this paper. It is common to see in real traffic that for some special reasons, for example, weather condition, health conditions, driving technology, and personal habits, the driver often does not keep up with the leading vehicle but drives at a velocity similar to that of lateral vehicles to form a stable vehicle group. erefore, the lateral vehicles are in role of auxiliary reference. is phenomenon is an incomplete car-following behavior. In view of the above discussions, by using the general force framework, the influence of both the leading vehicle and the lateral adjacent vehicle is used as stimulus in the new proposed model. is paper is organized as follows. In Section 2, the driving phenomenon is described. In Section 3, an improved car-following model based on the OVM is developed by taking into account the lateral effect. Linear stability analysis is carried out in Section 4, and the stability condition of the improved model is obtained. In Section 5, the MKDV equation near the critical point is derived by using nonlinear analysis. In Section 6, a series of numerical simulations are proposed to confirm the theoretical results. Conclusions are given in Section 7.

Description of Driving Phenomenon.
In actual traffic, one can often observe a phenomenon that although the headway between vehicle and its leader vehicle is larger, the driver still keeps a similar velocity to the surrounding vehicles and does not accelerate to shorten the distance from the leader vehicle. For example, the two investigated vehicles in Figure 1 are, respectively, marked as car 1 and car 2. Obviously, the driving state of car 2 is in accordance with the driving rules that can be described by the classic car-following model. However, despite a greater distance from the leader vehicle, car 1 still keeps the same velocity as its adjacent car on the right side. is phenomenon cannot be described by carfollowing models mentioned above because car 1 did not keep up with its leading car according to the rules describing by those car-following models. Figure 1 is intercepted from video captured by a camera. An overpass with 7.6 meter high in the South Second Ring Road of Xi'an City is selected for the field data collection. e location of the camera is identified as coordinate origin, and the distance between the coordinate origin and the lower edge of Figure 1 is 17.5 meters. For ensuring the measurement accuracy, measurement range is set from 17.5 meters to 89.5 meters because the larger shooting distance will reduce the measurement accuracy. e video is intercepted a set of images by frame/0.2 seconds. With the method of geometric transformation, one-to-one relationship can be established between image coordinates and vehicle's actual position. Data of the vehicles in Figure 1 can be seen in Table 1. Since the data obtained from images are discrete data, average velocity of vehicle in every image can be calculated. Velocity data are fitted by polynomial, and the function curves of velocity-time are obtained in Figure 2.

Analysis of Driving Phenomenon.
Figures 2(a) and 2(b) are, respectively, the headway-time curves and velocity-time curves of vehicles. From Figure 2(a), the headway of car 1 is larger than that of car 2. Most of the time, the headway of car 1 is larger than 15 meters, while the headway of car 2 is smaller than 5 meters. From Figure 2(b), it can be found that the two investigated vehicles always have similar velocity and that the trends of velocity changes are roughly the same.
Two classical car-following models, optimal velocity model (OVM) and full velocity difference model (FVDM), are used to simulate the driving state of car 1 in Figure 1: (2) Equations (1) and (2) are, respectively, the expressions of OVM and FVDM, where a denotes sensitivity; x j (t) and v j (t) are, respectively, the position and velocity of the vehicle j at time t; Δx j (t) and Δv j (t) are the headway and the velocity difference between the preceding vehicle and the following vehicle, respectively; V(Δx j (t)) presents the optimal velocity function; and λ denotes the sensitivity of velocity difference Δv j (t). Before the simulation, the parameters in OVM and FVDM need to be fitted. Parameters of OVM are the sensitivity coefficient a, maximum velocity of vehicle v max , and safety distance hc. Parameters of FVDM are, respectively, a, hc, and the sensitivity of velocity difference λ. e measured data of car 2 in Figure 1 can be used for parameter fitting because it is in accordance with most of vehicles' driving rules which are described by classic car-following models.  Journal of Advanced Transportation e results of fitting are a � 0.0877, v max � 16.7, and hc � 6.9781 for OVM, and a � 0.0593, hc � 13.0123, and λ � 0.4332 for FVDM. Substituting above parameters into OVM and FVDM, respectively, to simulate the velocity of car 1 and comparing the simulation result with real velocity, we can obtain Figure 3.
From Figure 3, we can see that there is a large error between the simulation results of the two classic car-following models and this set of traffic data for car 1. e simulation velocities of OVM continuously increase with the increase of headway. e difference between simulation velocity of FVDM and the velocity curve of car 1 keeps a larger value. Moreover, the variation trend of simulation velocity is opposite to that of real velocity after t � 3.8 seconds, which means an increasing error.
From the perspective of model structure, OVM only considers the influence of headway, and FVDM takes into account the influence of both headway and velocity difference between the following vehicle and its preceding vehicle. e previous studies indicate that besides the influence of the leading vehicle, the driving behavior also can be affected by the lateral adjacent vehicle. However, both OVM and FVDM did not consider the factors of lateral influence. With above consideration, new improved carfollowing models accounting for the lateral influence are proposed in the next section. In Section 6, new improved models are validated by multiple sets of measured traffic data.

Model
e second-order car-following model belongs to the stimulus response model, which can be expressed as follows: Equation (3) expresses that the jth vehicle's acceleration depends on speed of the jth vehicle, the headway, and relative velocity between the jth vehicle and the (j + 1)th vehicle. f sti represents response of the jth vehicle to the (j + 1)th vehicle. Real driving behavior is determined mainly by the following two motivations: achieving an expectation speed and maintaining an expected distance [22].
By using the general force framework, we, respectively, put forward two improved car-following models, named model 1 and model 2, and compared to choose the optimum model. e new improved model with the general force framework can be expressed as follows: where f + j (t) and f − j (t) represent the influence of the leading car and the lateral vehicle's effect, respectively. Because of mainly depending on Δx j (t), f + j (t) can be expressed as OVM. f − j (t) is determined by the impact of lateral adjacent vehicle and plays a role only when headway Δx j (t) is in a certain range. When hc > Δx j (t), in order to avoid the danger of collision, the driver often neglects the influence of lateral vehicle and changes the velocity to increase the headway. Only when hc < Δx j (t), the driver will consider the lateral effect.
where x j ′ (t) represents the position of the jth vehicle's lateral adjacent vehicle. Obviously, lateral vehicles can play an auxiliary reference role only under certain conditions; for example, the relative speed between the investigated car and the lateral car is larger, or the distance is far. erefore, we set the parameter c � const|Δv is the headway between the jth vehicle's lateral adjacent vehicle and its leading vehicle.
respectively, represent the headway and velocity difference between the jth vehicle and its lateral adjacent vehicle. 0 < ζ ≤ 0.1 is a constant. us, the first improved model (model 1) is given as follows: According to Reference [25], the second improved model (model 2) can be defined as where 0 ≤ p ≤ 1 is a coefficient which represents the weight of lateral effect. When p � 0, equation (7) is reduced to the OVM. It means that the driver is completely unaffected by the lateral vehicle, and the driving state is determined only by the driving state of the leading vehicle. When p � 1, the driver is only affected by lateral vehicle, and the vehicle exits the state of car-following.
Next, the two improved models are verified by the measured data to select the optimal model. In this paper, the least square method is used for parameter fitting. Result of fitting is as follows: a � 0.0877, hc � 6.9781, p � 0.7993, and c � 1.0246. e fitted parameters are taken into the models, and velocity of the vehicle is simulated by the two improved models. Simulation results are compared with the actual data in Figure 4. Figure 4 shows the comparison between simulation results of improved models and a set of traffic data for car 1 in Figure 1. e solid line is the curve of real velocity data, and the dashed line is the simulation result obtained by using the improved models. From Figure 4, we can find that for the data of car 1 in Figure 1, the curve of simulated velocity by using the first improved model has an increasing deviation from the real velocity after 3.6 seconds, while the curve obtained by the second improved model is close to the measured velocity curve during most of the time.

The Linear Stability Analysis
e linear stability theory is applied to derive linear stability condition of the improved model by the linear stability theory. e uniform traffic flow is defined by such a state that all vehicles on lane l move with the optimal velocity V l (Δx l,n (t)) and the identical headway h l .
e solution x (0) l,j (t) is given by where l � 1, 2, h l � D/N l , D is length of road section, and N l represents the total number of vehicles in lane l. Assuming x l,j (t) � x (0) l,j (t) + y l,j (t), where x (0) l,j (t) � h l + V(h l )t, and y l,j (t) be a small deviation from the steady state x (0) l,j (t), then we can rewrite the two improved models (equations (6) and (7)) as where Δy l,j (t) � y l,j+1 (t) − y l,j (t), y l′,j (t) � y 3− l,j (t). For convenience of calculation, we set y 3− l,j (t) � μy l,j , μ � h 3− l / h l .Expanding y l,j (t) in the Fourier-modes y l,j (t) ∝ exp (ik l j + z l t), we have Substituting z l � z 1l (ik) + z 2l (ik) 2 + · · · into equations (10) and (11), we obtain the first-and second-order terms of coefficients in the expression of z l as follows: According to the long wave theory, the uniform steady state will become unstable when z 2l is negative. us, the neutral stability curves for model 1 and model 2 are given by For model 1, the uniform traffic flow will be stable if For model 2, the uniform traffic flow will be stable if When c � 0 and P � 0, the stability condition (equations (15) and (16)) is reduced to the stability condition in the literature [3]. Figure 5 shows the neutral stability line in the headwaysensitivity space (Δx, a l ) with μ � 0.5 and different sets of parameters of lateral influence, where (a) and (b) are, respectively, for model 1 and model 2. e apex of each curve represents the critical point (hc, a lc ). It can be seen clearly from Figure 5(a) that for model 1, the critical points decrease with increasing c. It means that with decrease of lateral influence, traffic flow is more stable. On the contrary, from Figure 5(b), we can find that critical points rise with increasing p. It is in accordance with common sense that lateral interference can affect driving.

Nonlinear Analysis
In this section, nonlinear analysis is carried out to investigate the traffic performances near the critical point (h c , a c ) of improved car-following model 2.
e mKdV equation is deduced to describe the traffic congestion in the unstable area of traffic flow. For extracting slow scales with the space variable n and the time variable t, the slow variables X and T are defined as follows: Journal of Advanced Transportation 5 where b is a constant to be determined. We set the headway as We rewrite equation (7) as follows: For convenience of calculation, we assume that the vehicles around the investigated vehicle are in normal driving process and have similar driving state. erefore, the driving condition of the lateral vehicle can be approximately expressed by that of other vehicles around it. In view of above considerations, we set (d 2 Δx j ′ (t)/dt 2 ) � λ 1 (d 2 Δx j+1 (t)/dt 2 ) + λ 2 (d 2 Δx j− 1 (t)/dt 2 ), where λ 1 and λ 2 are constant, and λ 1 + λ 2 � 1. Substituting equations (17) and (18) into equation (19) and making the Taylor expansion to the fifth order of ε lead to the following expression: Near the critical point (h c , a c ), taking a c /a � (1 + ε 2 ), b � V ′ (h c ), and eliminating both the second-order and the third-order term of ε, equation (14) can be simplified as where Making such transformations as T � (1/g 1 )T m and R � ����� g 1 /g 2 R m , we obtain the modified KDV equation with a o(ε) correction term: If we ignore o(ε) term, this is just the modified KDV equation, and the solution is as follows: where B � (5g 2 g 3 /2g 2 g 4 − 3g 1 g 5 ). us, we obtain the kink solution of the headway:

Model Validation with Measured Data.
In this subsection, the OVM and improved models are, respectively, verified by using real traffic data. e traffic data can be seen in Tables 2-5. e model parameters need to be fitted by using the least square method. Due to the difference of driving behavior, the model parameters obtained from different sets of data are different. Hence, the data of each group need to be fitted, respectively, and the fitting results are shown in Table 6. e skewness test and kurtosis test of the fitted parameters are carried out. U 1 and U 2 are, respectively, denote sample skewness and sample kurtosis of the fitted parameters. Calculation results for skewness and kurtosis are shown in Table 7.
We set the significant level as α � 0.05, so the reject region is (− ∞, 1.96) ∪ (1.96, +∞). From Table 2, it can be seen that the calculation results of U 1 and U 2 are all outside the reject region. erefore, the fitted parameters obey normal distribution.
In order to compare the performance of the models, we substitute the fitted parameters into OVM and two improved models, respectively, carry out simulation by using the three models, and then obtain the curves of velocity. Comparing the simulation curves with the measured velocity curve, we obtain Figure 6.
In Figure 6, real line is the measured velocity of vehicle, "." represents the curve of OVM, and "+" and " * " denote the curve of improved model 1 and model 2, respectively. e simulation results show that the velocity curves of the two improved models are more close to the real velocity than that of the OVM and that in the comparison of the two improved models, and the velocity curve of the improved model 2 is closer to the actual speed.
In this paper, the following three performance indexes are used to evaluate the simulation results: the mean square error, the maximum absolute error, and the minimum absolute error. Table 8 shows the comparison of three performance evaluation indices for each group of data. From Table 8, we can find that for mean square deviation and maximum absolute error, numerical values of OVM are the largest in four sets of data, while the numerical values of improved model 2 are the smallest. Moreover, least absolute errors of improved model 1 are larger than that of OVM in the 1 st and 3 rd sets of data, while improved model 2 has the smallest least absolute errors in all four sets of data. Compared with the OVM, the simulation result of improved model 2, respectively, gets 62.89% unbiased variance reduction, 66.39% maximum absolute error reduction, and 33.4% minimum absolute error reduction, while that of improved model 1 only gets 27.53% unbiased variance reduction, 28.4%       maximum absolute error reduction, and 19.97% minimum absolute error reduction, respectively. e results indicate that improved model 2 is more suitable than OVM and improved model 1 for describing the traffic phenomena mentioned in this article.

Density Wave.
In order to check the validity of the above theoretical results, we discuss the space-time evolution of the headway for various values of p by using the Euler method in this section. In two adjacent lanes which are named lane 1 and lane 2, all vehicles drive in the same direction with the constant headway h � 4.0 m. Under the periodic boundary condition, the following initial conditions in the two lanes are given as follows: Δx 2,j (0) � Δx 2,j (1) � 4.0, for n ≠ 50, 51, Equations (26) and (27) are, respectively, the initial conditions in lane 1 and lane 2. e total number of vehicles in each lane is n � 100. e safety distance is h c � 4.0. Other input parameters for the simulation are a � 2, v max � 3. In order to simplify the experiment, we consider the vehicles in lane 2 are in normal driving process; therefore, the driving state of those vehicles can be described by the classic carfollowing model. e vehicles in lane 1 are investigated. en, space-time evolution of the headway in lane 1 after a sufficiently long time and the corresponding hysteresis loops will be obtained in Figures 7 and 8.   Figure 7 shows the space-time evolution of the headway for different sets of p after t � 10000 s. e patterns (a)-(d) are, respectively, corresponding to p � 0, 0.2, 0.4, and 0.6. Muramatsu and Nagatani [26] present the solution of MKDV equation can be exhibited as kink-antikink wave. From Figure 7, because of the lateral influence from the vehicles in lane 2, disturbance generates in the uniform flow in lane 1. Small disturbance added to the uniform traffic flow is amplified, and the kink-antikink wave appears as traffic jam. Only parts of the vehicles in lane1 are affected by the lateral vehicle, so the amplitudes of the density wave are not uniform. In addition, it can be seen clearly from Figure 7 that the amplitudes of the density wave increase with increasing p. It indicates that the greater the lateral influence, the more unstable the traffic flow is. Figure 8 exhibits the hysteresis loops after t � 10000 s with different sets of p. e solid line and dotted line are, respectively, obtained with p � 0 and p � 0.6. Hysteresis loop has a larger size means that the stability of traffic flow is lower, and the stability of traffic flow is improved with the decrease of the size of the loop. From Figure 8, we can see that the hysteresis loop obtained with p � 0.6 is much larger than that with p � 0. e result indicates that by taking lateral influence into account, the stability of the traffic flow decreases, which is in accordance with the linear stability analysis and nonlinear analysis in Sections 4 and 5, respectively.

Conclusion
In this paper, we describe a kind of traffic phenomenon that the driver often does not keep up with the leading vehicle but drives at a velocity similar to that of the surrounding vehicles. Because of special reasons such as weather condition, health conditions, driving technology, and personal habits, the driver does not always keep up with the leading vehicle but maintains a larger distance from the leading vehicle to leave more reaction time. is phenomenon is an incomplete car-following behavior. With the consideration of the lateral vehicle's influence, we put forward two improved models with different structures on the base of the OV model in this paper.
e new improved models are verified by multiple sets of measured data. Linear stability analysis and nonlinear analysis are carried out to study the traffic characteristics.
e results of linear stability analysis and numerical simulation show that with increase of lateral influence, traffic flow is more unstable, which is in accordance with the traffic phenomenon that lateral interference can affect driving.
Data Availability e real traffic data used to support the findings of this study are included within the article. ere are no restrictions on data access.

Conflicts of Interest
e authors declare that they have no conflicts of interest.