# Modeling for Micro Traffic Flow with the Consideration of Lateral Vehicle’s Influence

Authors: Da-wei Liu; Zhong-ke Shi; Wen-Huan Ai
Journal: Journal of Advanced Transportation
 (2020)

---

## Abstract

In order to make the car-following model describe the driving behavior of vehicle on urban road more accurately, existing car-following 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. The lateral vehicle’s influence is introduced as the influence factor of driving behavior. By using of linear stability analysis, stability conditions of improved car-following models are obtained. Nonlinear analysis is carried out to investigate the traffic performances near the critical point. The 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-following models 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. Therefore, the second cooperative model is more suitable to describe the vehicle’s actual behavior in car-following state.


---

## Body

## 1. 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. The 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–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 low- and 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. Therefore, the lateral vehicles are in role of auxiliary reference. This 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. This paper is organized as follows. In Section2, 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.


## 2. Description and Analysis of Driving Phenomenon

### 2.1. 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 Figure1 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. This phenomenon cannot be described by car-following 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
Real traffic.


### 2.2. Analysis of Driving Phenomenon

Figure1 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. The 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. The 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 Figure1 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.Table 1
Data of the vehicle in Figure1.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayVelocity of leading carVelocity of following carPosition of leading carPosition of following carHeadwayVelocity of leading carVelocity of following car037.0917.7519.3411.189.3622.5517.445.1110.3610.670.239.5519.8619.6910.189.1024.9119.895.029.9210.410.441.3221.6119.719.488.9126.7921.884.919.6010.220.642.7823.0519.739.028.7928.2123.544.679.3610.080.84525.1119.898.748.7230.5725.844.739.199.991.046.5426.5120.038.598.6932.1127.484.639.079.921.248.3728.5819.798.528.7034.1829.94.288.979.861.450.1430.4119.738.498.7435.8531.913.948.899.801.651.631.7719.838.478.7937.3333.3348.809.731.853.8333.9719.868.438.8539.4235.763.668.719.652.055.0235.6719.358.358.9241.1837.863.328.609.542.256.5137.0719.448.218.9842.6339.233.48.469.402.458.6239.2319.398.019.0344.5141.293.228.319.232.659.7340.6919.047.759.0646.0142.923.098.139.032.861.4842.9218.567.439.0647.9944.973.027.938.813.063.0144.6918.327.069.0449.5446.762.787.728.553.263.9646.1517.816.658.9950.5548.042.517.508.283.465.6148.2117.46.248.8952.0350.11.937.298.003.666.9949.7417.255.858.7553.3751.381.997.097.713.868.0751.7616.315.518.5655.2653.132.136.947.434.069.1853.5415.645.288.3256.7654.582.186.837.174.269.9354.815.135.208.0357.8159.891.926.816.944.471.156.5714.535.327.6859.4557.272.186.896.774.671.957.9813.925.717.2760.358.471.837.096.664.873.5559.4614.096.436.8062.0859.972.117.456.655.074.8360.7614.077.566.2763.6461.292.358.016.75Figure 2
Traffic data of car 1 and car 2: (a) headway-time curves of car 1 and car 2; (b) velocity-time curves of car 1 and car 2.
(a)(b)Figures2(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 Figure1:(1)d2xjtdt2=aVΔxjt−vjt,(2)d2xjtdt2=aVΔxjt−vjt+λΔvjt.Equations (1) and (2) are, respectively, the expressions of OVM and FVDM, where a denotes sensitivity; xjt and vjt are, respectively, the position and velocity of the vehicle j at time t; Δxjt and Δvjt are the headway and the velocity difference between the preceding vehicle and the following vehicle, respectively; VΔxjt presents the optimal velocity function; and λ denotes the sensitivity of velocity difference Δvjt.Before the simulation, the parameters in OVM and FVDM need to be fitted. Parameters of OVM are the sensitivity coefficienta, maximum velocity of vehicle vmax, and safety distance hc. Parameters of FVDM are, respectively, a, hc, and the sensitivity of velocity difference λ. The 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. The results of fitting are a=0.0877, vmax=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.Figure 3
Comparison between simulation curves of classic car-following models and the curve of real velocity.From Figure3, 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. The simulation velocities of OVM continuously increase with the increase of headway. The 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. The 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 car-following models accounting for the lateral influence are proposed in the next section. In Section6, new improved models are validated by multiple sets of measured traffic data.


## 2.1. 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 Figure1 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. This phenomenon cannot be described by car-following 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
Real traffic.


## 2.2. Analysis of Driving Phenomenon

Figure1 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. The 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. The 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 Figure1 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.Table 1
Data of the vehicle in Figure1.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayVelocity of leading carVelocity of following carPosition of leading carPosition of following carHeadwayVelocity of leading carVelocity of following car037.0917.7519.3411.189.3622.5517.445.1110.3610.670.239.5519.8619.6910.189.1024.9119.895.029.9210.410.441.3221.6119.719.488.9126.7921.884.919.6010.220.642.7823.0519.739.028.7928.2123.544.679.3610.080.84525.1119.898.748.7230.5725.844.739.199.991.046.5426.5120.038.598.6932.1127.484.639.079.921.248.3728.5819.798.528.7034.1829.94.288.979.861.450.1430.4119.738.498.7435.8531.913.948.899.801.651.631.7719.838.478.7937.3333.3348.809.731.853.8333.9719.868.438.8539.4235.763.668.719.652.055.0235.6719.358.358.9241.1837.863.328.609.542.256.5137.0719.448.218.9842.6339.233.48.469.402.458.6239.2319.398.019.0344.5141.293.228.319.232.659.7340.6919.047.759.0646.0142.923.098.139.032.861.4842.9218.567.439.0647.9944.973.027.938.813.063.0144.6918.327.069.0449.5446.762.787.728.553.263.9646.1517.816.658.9950.5548.042.517.508.283.465.6148.2117.46.248.8952.0350.11.937.298.003.666.9949.7417.255.858.7553.3751.381.997.097.713.868.0751.7616.315.518.5655.2653.132.136.947.434.069.1853.5415.645.288.3256.7654.582.186.837.174.269.9354.815.135.208.0357.8159.891.926.816.944.471.156.5714.535.327.6859.4557.272.186.896.774.671.957.9813.925.717.2760.358.471.837.096.664.873.5559.4614.096.436.8062.0859.972.117.456.655.074.8360.7614.077.566.2763.6461.292.358.016.75Figure 2
Traffic data of car 1 and car 2: (a) headway-time curves of car 1 and car 2; (b) velocity-time curves of car 1 and car 2.
(a)(b)Figures2(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 Figure1:(1)d2xjtdt2=aVΔxjt−vjt,(2)d2xjtdt2=aVΔxjt−vjt+λΔvjt.Equations (1) and (2) are, respectively, the expressions of OVM and FVDM, where a denotes sensitivity; xjt and vjt are, respectively, the position and velocity of the vehicle j at time t; Δxjt and Δvjt are the headway and the velocity difference between the preceding vehicle and the following vehicle, respectively; VΔxjt presents the optimal velocity function; and λ denotes the sensitivity of velocity difference Δvjt.Before the simulation, the parameters in OVM and FVDM need to be fitted. Parameters of OVM are the sensitivity coefficienta, maximum velocity of vehicle vmax, and safety distance hc. Parameters of FVDM are, respectively, a, hc, and the sensitivity of velocity difference λ. The 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. The results of fitting are a=0.0877, vmax=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.Figure 3
Comparison between simulation curves of classic car-following models and the curve of real velocity.From Figure3, 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. The simulation velocities of OVM continuously increase with the increase of headway. The 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. The 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 car-following models accounting for the lateral influence are proposed in the next section. In Section6, new improved models are validated by multiple sets of measured traffic data.


## 3. Model

The second-order car-following model belongs to the stimulus response model, which can be expressed as follows:(3)d2xjdt2t=fstivjt,Δxjt,Δvjt.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. fsti 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.The new improved model with the general force framework can be expressed as follows:(4)d2xjdt2t=fj+t+fj−t,where fj+t and fj−t represent the influence of the leading car and the lateral vehicle’s effect, respectively. Because of mainly depending on Δxjt, fj+t can be expressed as OVM. fj−t is determined by the impact of lateral adjacent vehicle and plays a role only when headway Δxjt is in a certain range. When hc>Δxjt, 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<Δxjt, the driver will consider the lateral effect.fj−t can be defined as(5)fn−t=γ∗d2xj′tdt2,where xj′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. Therefore, we set the parameter γ=constΔv¯jt<ζvjtorΔx¯jt<Δxj′t,0,other, where Δxj′t=xj+1′t−xj′t is the headway between the jth vehicle’s lateral adjacent vehicle and its leading vehicle. Δx¯jt=xj′t−xjt and Δv¯jt=dxj′t/dt−dxjt/dt, respectively, represent the headway and velocity difference between the jth vehicle and its lateral adjacent vehicle. 0<ζ≤0.1 is a constant. Thus, the first improved model (model 1) is given as follows:(6)d2xjtdt2=aVΔxjt−dxjtdt+γ∗d2xj′tdt2.According to Reference [25], the second improved model (model 2) can be defined as(7)d2xjtdt2=1−P∗aVΔxjt−dxjtdt+p∗d2xj′tdt2,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 γ=1.0246. The 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
Comparison between simulation curves of improved models and the curve of real velocity: (a) velocity-time curves of improved model 1; (b) velocity-time curves of improved model 2.
(a)(b)Figure4 shows the comparison between simulation results of improved models and a set of traffic data for car 1 in Figure 1. The 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.


## 4. The Linear Stability Analysis

The linear stability theory is applied to derive linear stability condition of the improved model by the linear stability theory. The uniform traffic flow is defined by such a state that all vehicles on lanel move with the optimal velocity VlΔxl,nt and the identical headway hl. The solution xl,j0t is given by(8)xl,j0t=hl+Vhlt,where l=1,2, hl=D/Nl, D is length of road section, and Nl represents the total number of vehicles in lane l. Assuming xl,jt=xl,j0t+yl,jt, where xl,j0t=hl+Vhlt, and yl,jt be a small deviation from the steady state xl,j0t, then we can rewrite the two improved models (equations (6) and (7)) as(9)d2yl,jtdt2=alV′hlΔyl,jt−dyl,jtdt+γd2yl′,jtdt2,d2yl,jtdt2=1−palV′hlΔyl,jt−dyl,jtdt+pd2yl′,jtdt2,where Δyl,jt=yl,j+1t−yl,jt, yl′,jt=y3−l,jt. For convenience of calculation, we set y3−l,jt=μyl,j, μ=h3−l/hl.Expanding yl,jt in the Fourier-modes yl,jt∝expiklj+zlt, we have(10)zl21−μγ−alV′hlexpikl−1−zl=0,(11)zl21−μp−1−palV′hlexpikl−1−zl=0.Substitutingzl=z1lik+z2lik2+⋯ into equations (10) and (11), we obtain the first- and second-order terms of coefficients in the expression of zl as follows:(12)z1l=V′hl,z2l=V′hl2−1−μγV′2hlal,z1l=V′hl,z2l=V′hl2−1−μpV′2hlal1−p.According to the long wave theory, the uniform steady state will become unstable whenz2l is negative. Thus, the neutral stability curves for model 1 and model 2 are given by(13)als=21−μγV′hl,(14)als=21−μpV′hl1−p.For model 1, the uniform traffic flow will be stable if(15)al>21−μγV′hl.For model 2, the uniform traffic flow will be stable if(16)al>21−μpV′hl1−p.Whenγ=0 and P=0, the stability condition (equations (15) and (16)) is reduced to the stability condition in the literature [3].Figure5 shows the neutral stability line in the headway-sensitivity space Δx,al with μ=0.5 and different sets of parameters of lateral influence, where (a) and (b) are, respectively, for model 1 and model 2. The apex of each curve represents the critical point hc,alc. It can be seen clearly from Figure 5(a) that for model 1, the critical points decrease with increasing γ. 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.Figure 5
Critical stability curves: (a) for model 1; (b) for model 2.
(a)(b)


## 5. Nonlinear Analysis

In this section, nonlinear analysis is carried out to investigate the traffic performances near the critical pointhc,ac of improved car-following model 2. The 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:(17)X=j+bt,T=ε3t,where b is a constant to be determined. We set the headway as(18)Δxnt=hc+εRX,T.We rewrite equation (7) as follows:(19)d2Δxjtdt2=1−paVΔxj+1t−VΔxjt−−dΔxjtdt+pd2Δxj′tdt2.For convenience of calculation, we assume that the vehicles around the investigated vehicle are in normal driving process and have similar driving state. Therefore, the driving condition of the lateral vehicle can be approximately expressed by that of other vehicles around it. In view of above considerations, we setd2Δxj′t/dt2=λ1d2Δxj+1t/dt2+λ2d2Δxj−1t/dt2, 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:(20)ε2a1−pb−V′h∂XR+ε3b2−12a1−pV′h−pb2∂X2R+ε4a1−p∂TR−16V′h∂X3R−16V″′h∂XR3−pλ1−λ2b2∂X3R+ε52b1−p∂T∂XR−12b2∂X4R−a1−p124V′h∂X4R+112V″′h∂X2R3.Near the critical pointhc,ac, taking ac/a=1+ε2, b=V′hc, and eliminating both the second-order and the third-order term of ε, equation (14) can be simplified as(21)∂TR−g1∂X3R+g2∂X3R+εg3∂X2R+g4∂X4R+g5∂X2R3=0,where(22)g1=16+pλ1−λ2V′a1−pV′,g2=−16V′″,g3=−12V′,g4=4V′−a12aV′″,g5=V′3a+2pλ1−λ2a21−pV′2−V′2a1−p−124V′.Making such transformations asT=1/g1Tm and R=g1/g2Rm, we obtain the modified KDV equation with a oε correction term:(23)∂TRm−∂X3Rm+g2∂XRm3+εg1g3∂X2Rm+g4∂X4Rm+g1g5g2∂X2Rm3=0.If we ignoreoε term, this is just the modified KDV equation, and the solution is as follows:(24)RX,T=g1g2BtanhB2X−Bg1T,where B=5g2g3/2g2g4−3g1g5.Thus, we obtain the kink solution of the headway:(25)Δxnt=hc+g1Bg2aca−1×tanhB2aca−1×n+V′hct−Bg1aca−1t.


## 6. Numerical Simulations

### 6.1. Model Validation with Measured Data

In this subsection, the OVM and improved models are, respectively, verified by using real traffic data. The traffic data can be seen in Tables2–5. The 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.Table 2
The first sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayVelocity of following carPosition of leading carPosition of following carHeadwayVelocity of following car026.2218.937.296.8540.7118.2722.448.430.227.8520.57.356.8642.6520.1822.478.300.429.0221.567.466.8443.7621.4622.38.190.630.6523.177.486.7945.5423.4522.098.110.831.6124.247.376.7246.8124.6822.138.041.033.225.827.386.6348.6626.5522.117.991.234.527.157.356.5450.2928.1822.117.941.435.828.437.376.4551.8429.7322.117.911.636.9929.787.216.3653.2631.3821.887.881.837.9430.827.126.2854.3332.6421.697.862.039.2832.017.276.2155.934.1721.737.842.240.8333.557.286.1657.5636.221.367.822.441.8534.57.356.1258.8137.421.417.792.643.62367.626.1060.6539.2121.447.772.844.7837.17.686.0962.0340.7521.287.743.045.8538.177.686.1063.4841.9221.567.703.247.6339.538.16.1364.9943.8221.177.663.449.0240.748.286.1766.2645.3220.947.613.650.4842.028.466.2267.9146.7821.137.553.851.8443.398.456.2769.2948.3420.957.484.053.0544.48.656.3370.3749.520.877.404.254.7745.7796.3971.8551.2620.597.314.456.1346.99.236.4372.6252.5620.067.224.657.8148.449.376.4674.254.1420.067.124.859.3249.79.626.4775.4355.8119.627.015.060.3850.849.546.4476.756.6920.016.895.261.7552.239.526.3877.1358.318.836.77Table 3
The second sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following caHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car022.517.215.298.0935.8919.0716.829.350.224.2918.85.497.4638.0220.9617.068.890.426.0320.445.597.0240.1322.9717.168.600.627.2521.495.766.7341.5524.2317.328.440.829.0823.0866.5743.5526.2217.338.381.030.3324.116.226.5045.0527.5517.58.391.231.9325.456.486.4946.8229.3117.518.441.433.5926.756.846.5248.5431.0517.498.521.635.4228.357.076.5750.5833.0817.58.601.837.0929.677.426.6452.1234.8217.38.672.038.5330.727.816.7053.536.217.38.722.240.2232.078.156.7555.2137.9917.228.752.442.3333.698.646.7957.0239.9617.068.762.643.5534.738.826.8258.3941.3317.068.732.845.6736.389.296.8460.4143.5516.868.693.047.0737.59.576.8761.6344.916.738.633.248.7538.99.856.9163.2146.6316.588.583.450.5540.310.256.9864.8848.3416.548.543.652.4741.9210.557.1166.2750.3615.918.543.854.0743.4210.657.3268.152.1815.928.614.055.7544.6211.137.6369.2453.5415.78.774.257.5446.1911.358.0970.8355.1715.669.054.459.7348.0211.718.7272.4857.3715.119.49Table 4
The third sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car039.8619.1420.729.6224.1117.486.639.620.241.7221.2420.488.8826.2219.616.619.220.443.3422.9420.48.4027.9721.466.519.010.644.6324.3120.328.1229.3522.96.458.920.846.1625.9920.177.9831.0924.766.338.931.047.827.6520.157.9432.8726.556.328.981.249.1929.3719.827.9634.5928.46.199.051.451.0531.3219.738.0136.6930.66.099.111.652.4432.6419.88.0637.9432.075.879.151.854.1134.3519.768.0939.6333.855.789.142.055.4435.9319.518.0841.4735.675.89.092.257.5637.8919.678.0243.3437.695.658.972.458.8139.5319.287.9245.0839.535.558.802.660.1140.7519.367.7746.3240.865.468.582.861.7542.4119.347.5947.9742.535.448.333.063.4844.4918.997.3949.9244.495.438.063.264.6945.7518.947.1851.2545.895.367.803.466.5847.7118.877.0152.8547.715.147.583.667.9149.1618.756.9154.1149.164.957.423.868.9450.1918.756.9054.9950.364.637.384.070.7351.9918.747.0556.8451.994.857.494.271.8553.3418.517.4058.353.344.967.824.473.855.3918.418.0159.8455.174.678.41Table 5
The fourth sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car036.3816.4219.9610.5026.3516.539.8210.140.238.2618.5619.710.5228.2718.659.6210.210.439.9820.6419.3410.5330.1320.739.410.190.642.1222.7219.410.5431.9422.769.1810.110.844.0524.8519.210.5333.8624.79.169.971.045.8526.9118.9410.5135.6726.718.969.791.247.6329.0518.5810.4737.4128.568.859.601.449.3731.2118.1610.4139.1830.468.729.401.651.2533.1718.0810.3440.9732.328.659.211.852.8535.2317.6210.2442.7634.158.619.022.054.7737.3617.4110.1244.21368.218.852.256.1339.3316.89.9945.9237.638.298.712.457.5641.6515.919.8447.95407.958.582.658.8142.9215.899.6748.5440.997.558.482.860.9245.1215.89.4950.5842.977.618.403.062.3147.0815.239.3052.1244.47.728.343.264.0848.8815.29.1053.7446.257.498.293.465.6250.6414.988.9155.2147.917.38.243.667.2452.3414.98.7256.7649.517.258.183.868.9454.1614.788.5457.8451.036.818.094.070.7356.3414.398.3959.5352.646.897.974.271.4857.513.988.2560.7154.146.577.804.473.0158.9614.058.1662.2556.196.067.554.674.260.7613.448.1063.5457.46.147.214.875.8562.3913.468.1064.8858.676.216.755.077.1364.1113.028.1666.2759.996.286.15Table 6
Fitting parameters in every sets of data.
ahcγpThe 1st sets of data0.02338.44820.42590.9263The 2nd sets of data0.15449.00740.60040.8195The 3rd sets of data0.029816.68181.18861The 4th sets of data0.12939.42270.76890.6038The skewness test and kurtosis test of the fitted parameters are carried out.U1 and U2 are, respectively, denote sample skewness and sample kurtosis of the fitted parameters. Calculation results for skewness and kurtosis are shown in Table 7.Table 7
The skewness and kurtosis of the fitted parameters.
ahcγpU10.10911.90930.9596−0.9933U2−1.16300.83650.20850.1299We 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 U1 and U2 are all outside the reject region. Therefore, 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 Figure6.Figure 6
Comparison of simulation result: (a) the first sets of data; (b) the second sets of data; (c) the third sets of data; (d) the fourth sets of data.
(a)(b)(c)(d)In Figure6, 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. The 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.Table8 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 1st and 3rd 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. The results indicate that improved model 2 is more suitable than OVM and improved model 1 for describing the traffic phenomena mentioned in this article.Table 8
Comparison of evaluation indexes.
Mean square deviationMaximum absolute errorMinimum absolute errorOVMModel 1Model 2OVMModel 1Model 2OVMModel 1Model 2The 1st sets2.63282.44391.71593.14692.97311.76041.61161.62701.6094The 2nd sets5.26314.43451.64316.75896.09421.86302.11851.67171.3226The 3rd sets2.51560.85570.57223.70241.13960.81460.78030.18690.0462The 4th sets3.87832.61831.36966.94914.51162.47110.14390.23920.1219Average value3.57132.58811.32525.13933.67961.72731.16360.93120.7750


### 6.2. 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 ofp 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.0m. Under the periodic boundary condition, the following initial conditions in the two lanes are given as follows:(26)Δx1,j0=Δx1,j1=⋯=Δx1,jn=4.0,(27)Δx2,j0=Δx2,j1=4.0,forn≠50,51,Δx2,j1=4.0+0.1,forn=50,Δx2,j1=4.0−0.1,forn=51.Equations (26) and (27) are, respectively, the initial conditions in lane 1 and lane 2. The total number of vehicles in each lane is n=100.The safety distance is hc=4.0. Other input parameters for the simulation are a=2,vmax=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 car-following model. The vehicles in lane 1 are investigated. Then, 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
Space-time evolution of the headway aftert=10000s for small disturbance with different p: (a) p=0, (b) p=0.2, (c) p=0.4, and (d) p=0.6.
(a)(b)(c)(d)Figure 8
Hysteresis loops aftert=10000s with different sets of p.Figure7 shows the space-time evolution of the headway for different sets of p after t=10000s. The 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.Figure8 exhibits the hysteresis loops after t=10000s with different sets of p. The 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. The 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.


## 6.1. Model Validation with Measured Data

In this subsection, the OVM and improved models are, respectively, verified by using real traffic data. The traffic data can be seen in Tables2–5. The 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.Table 2
The first sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayVelocity of following carPosition of leading carPosition of following carHeadwayVelocity of following car026.2218.937.296.8540.7118.2722.448.430.227.8520.57.356.8642.6520.1822.478.300.429.0221.567.466.8443.7621.4622.38.190.630.6523.177.486.7945.5423.4522.098.110.831.6124.247.376.7246.8124.6822.138.041.033.225.827.386.6348.6626.5522.117.991.234.527.157.356.5450.2928.1822.117.941.435.828.437.376.4551.8429.7322.117.911.636.9929.787.216.3653.2631.3821.887.881.837.9430.827.126.2854.3332.6421.697.862.039.2832.017.276.2155.934.1721.737.842.240.8333.557.286.1657.5636.221.367.822.441.8534.57.356.1258.8137.421.417.792.643.62367.626.1060.6539.2121.447.772.844.7837.17.686.0962.0340.7521.287.743.045.8538.177.686.1063.4841.9221.567.703.247.6339.538.16.1364.9943.8221.177.663.449.0240.748.286.1766.2645.3220.947.613.650.4842.028.466.2267.9146.7821.137.553.851.8443.398.456.2769.2948.3420.957.484.053.0544.48.656.3370.3749.520.877.404.254.7745.7796.3971.8551.2620.597.314.456.1346.99.236.4372.6252.5620.067.224.657.8148.449.376.4674.254.1420.067.124.859.3249.79.626.4775.4355.8119.627.015.060.3850.849.546.4476.756.6920.016.895.261.7552.239.526.3877.1358.318.836.77Table 3
The second sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following caHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car022.517.215.298.0935.8919.0716.829.350.224.2918.85.497.4638.0220.9617.068.890.426.0320.445.597.0240.1322.9717.168.600.627.2521.495.766.7341.5524.2317.328.440.829.0823.0866.5743.5526.2217.338.381.030.3324.116.226.5045.0527.5517.58.391.231.9325.456.486.4946.8229.3117.518.441.433.5926.756.846.5248.5431.0517.498.521.635.4228.357.076.5750.5833.0817.58.601.837.0929.677.426.6452.1234.8217.38.672.038.5330.727.816.7053.536.217.38.722.240.2232.078.156.7555.2137.9917.228.752.442.3333.698.646.7957.0239.9617.068.762.643.5534.738.826.8258.3941.3317.068.732.845.6736.389.296.8460.4143.5516.868.693.047.0737.59.576.8761.6344.916.738.633.248.7538.99.856.9163.2146.6316.588.583.450.5540.310.256.9864.8848.3416.548.543.652.4741.9210.557.1166.2750.3615.918.543.854.0743.4210.657.3268.152.1815.928.614.055.7544.6211.137.6369.2453.5415.78.774.257.5446.1911.358.0970.8355.1715.669.054.459.7348.0211.718.7272.4857.3715.119.49Table 4
The third sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car039.8619.1420.729.6224.1117.486.639.620.241.7221.2420.488.8826.2219.616.619.220.443.3422.9420.48.4027.9721.466.519.010.644.6324.3120.328.1229.3522.96.458.920.846.1625.9920.177.9831.0924.766.338.931.047.827.6520.157.9432.8726.556.328.981.249.1929.3719.827.9634.5928.46.199.051.451.0531.3219.738.0136.6930.66.099.111.652.4432.6419.88.0637.9432.075.879.151.854.1134.3519.768.0939.6333.855.789.142.055.4435.9319.518.0841.4735.675.89.092.257.5637.8919.678.0243.3437.695.658.972.458.8139.5319.287.9245.0839.535.558.802.660.1140.7519.367.7746.3240.865.468.582.861.7542.4119.347.5947.9742.535.448.333.063.4844.4918.997.3949.9244.495.438.063.264.6945.7518.947.1851.2545.895.367.803.466.5847.7118.877.0152.8547.715.147.583.667.9149.1618.756.9154.1149.164.957.423.868.9450.1918.756.9054.9950.364.637.384.070.7351.9918.747.0556.8451.994.857.494.271.8553.3418.517.4058.353.344.967.824.473.855.3918.418.0159.8455.174.678.41Table 5
The fourth sets of vehicle data.
Time (s)The left laneThe right lanePosition of leading carPosition of following carHeadwayPosition of leading carPosition of following carHeadwayPosition of leading carPosition of following car036.3816.4219.9610.5026.3516.539.8210.140.238.2618.5619.710.5228.2718.659.6210.210.439.9820.6419.3410.5330.1320.739.410.190.642.1222.7219.410.5431.9422.769.1810.110.844.0524.8519.210.5333.8624.79.169.971.045.8526.9118.9410.5135.6726.718.969.791.247.6329.0518.5810.4737.4128.568.859.601.449.3731.2118.1610.4139.1830.468.729.401.651.2533.1718.0810.3440.9732.328.659.211.852.8535.2317.6210.2442.7634.158.619.022.054.7737.3617.4110.1244.21368.218.852.256.1339.3316.89.9945.9237.638.298.712.457.5641.6515.919.8447.95407.958.582.658.8142.9215.899.6748.5440.997.558.482.860.9245.1215.89.4950.5842.977.618.403.062.3147.0815.239.3052.1244.47.728.343.264.0848.8815.29.1053.7446.257.498.293.465.6250.6414.988.9155.2147.917.38.243.667.2452.3414.98.7256.7649.517.258.183.868.9454.1614.788.5457.8451.036.818.094.070.7356.3414.398.3959.5352.646.897.974.271.4857.513.988.2560.7154.146.577.804.473.0158.9614.058.1662.2556.196.067.554.674.260.7613.448.1063.5457.46.147.214.875.8562.3913.468.1064.8858.676.216.755.077.1364.1113.028.1666.2759.996.286.15Table 6
Fitting parameters in every sets of data.
ahcγpThe 1st sets of data0.02338.44820.42590.9263The 2nd sets of data0.15449.00740.60040.8195The 3rd sets of data0.029816.68181.18861The 4th sets of data0.12939.42270.76890.6038The skewness test and kurtosis test of the fitted parameters are carried out.U1 and U2 are, respectively, denote sample skewness and sample kurtosis of the fitted parameters. Calculation results for skewness and kurtosis are shown in Table 7.Table 7
The skewness and kurtosis of the fitted parameters.
ahcγpU10.10911.90930.9596−0.9933U2−1.16300.83650.20850.1299We 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 U1 and U2 are all outside the reject region. Therefore, 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 Figure6.Figure 6
Comparison of simulation result: (a) the first sets of data; (b) the second sets of data; (c) the third sets of data; (d) the fourth sets of data.
(a)(b)(c)(d)In Figure6, 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. The 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.Table8 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 1st and 3rd 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. The results indicate that improved model 2 is more suitable than OVM and improved model 1 for describing the traffic phenomena mentioned in this article.Table 8
Comparison of evaluation indexes.
Mean square deviationMaximum absolute errorMinimum absolute errorOVMModel 1Model 2OVMModel 1Model 2OVMModel 1Model 2The 1st sets2.63282.44391.71593.14692.97311.76041.61161.62701.6094The 2nd sets5.26314.43451.64316.75896.09421.86302.11851.67171.3226The 3rd sets2.51560.85570.57223.70241.13960.81460.78030.18690.0462The 4th sets3.87832.61831.36966.94914.51162.47110.14390.23920.1219Average value3.57132.58811.32525.13933.67961.72731.16360.93120.7750


## 6.2. 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 ofp 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.0m. Under the periodic boundary condition, the following initial conditions in the two lanes are given as follows:(26)Δx1,j0=Δx1,j1=⋯=Δx1,jn=4.0,(27)Δx2,j0=Δx2,j1=4.0,forn≠50,51,Δx2,j1=4.0+0.1,forn=50,Δx2,j1=4.0−0.1,forn=51.Equations (26) and (27) are, respectively, the initial conditions in lane 1 and lane 2. The total number of vehicles in each lane is n=100.The safety distance is hc=4.0. Other input parameters for the simulation are a=2,vmax=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 car-following model. The vehicles in lane 1 are investigated. Then, 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
Space-time evolution of the headway aftert=10000s for small disturbance with different p: (a) p=0, (b) p=0.2, (c) p=0.4, and (d) p=0.6.
(a)(b)(c)(d)Figure 8
Hysteresis loops aftert=10000s with different sets of p.Figure7 shows the space-time evolution of the headway for different sets of p after t=10000s. The 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.Figure8 exhibits the hysteresis loops after t=10000s with different sets of p. The 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. The 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.


## 7. 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. This 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. The 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. The 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.



---

*Source: 8340283-2020-09-19.xml*