Stability Analysis of Mixed Trafﬁc Flow Considering Personal Space under the Connected and Automated Environment

: A trafﬁc survey using an unmanned aerial vehicle on the Haixia road in Chongqing, China found a signiﬁcant correlation between the velocity of the vehicle and the distance of the vehicle when moving forward and laterally. To mitigate driving disruptions owing to vehicle intrusion at a desired distance, personal space (PS) is introduced to analyze the car-following behavior of continuous mixed trafﬁc ﬂow formed by human-driven (HD) vehicles and cooperative adaptive cruise control (CACC) vehicles. PS is a virtual boundary that refers to the space where psychological tension is caused by the invasion of others. Further, an intelligent driver model (IDM) was used to establish a mixed trafﬁc car-following PS-IDM. The stability of the disturbance transfer function analytical model in homogeneous and heterogeneous trafﬁc ﬂows was used to calculate the trafﬁc ﬂow stability region under different CACC vehicle permutations. The results show that the PS-IDM-based car-following model effectively improves the stability ﬁtting effect of mixed trafﬁc ﬂow, and the driving comfort is increased by up to 20.7% when compared with that of the single car-following model. In addition, there is a negative correlation between the PS and the unstable velocity range of the trafﬁc ﬂow. Compared with the homogeneous HD trafﬁc ﬂow, the average intrusion rate of homogeneous CACC trafﬁc ﬂow is reduced by 4.6%, and the driving comfort of vehicles is improved by approximately 65%.


Introduction
With the development of internet of vehicles technology, cooperative adaptive cruise control (CACC) vehicles have become the direction of advancement for future vehicles [1]. However, owing to technical and market limitations, human-driven (HD) vehicles will not be completely replaced by CACC in the short term. In the future, the urban road network will be a new type of mixed traffic flow composed of CACC and HD, which will subvert the current traffic organization scenario and bring unprecedented challenges to traffic managers [2].
In the context of mixed traffic flow, the car-following model remains a key factor in the study of traffic dynamics. The car-following model is an important representative of the microscopic traffic flow model. It identifies the microscopic state and internal mechanism of the traffic flow by quantifying the relationship between the longitudinal distances of vehicles in a single lane.
The car-following model was first proposed by Reuschel [3]. It uses dynamics to analyze and study train-like vehicles, marking the initial formation of the car-following theory. Hoogendoorn and Bovy [4] summarized the research results of traffic flow modeling from 1955 to 2000. They divided the traffic flow model into (sub-)microscopic, mesoscopic, and macroscopic groups according to the level of detail with which they represented the traffic systems. The development process of the car-following model was introduced in the microscopic model. On this basis, Yang et al. [5] reviewed 70 years of research on the car-following model. An overview of the main car-following model is shown in Table 1. The model calculated the safe following distance based on the velocity of the preceding and following vehicles.

Newell model 1961
The trajectory of the rear car was only a copy of the trajectory of the preceding car after a certain time interval and distance interval, and the trajectory trend in the front and rear cars was always consistent.
Optimal Velocity model 1995 The model gave a functional relationship between head spacing and optimal velocity.

IDM model 2000
The same model was used to describe the driver's following behavior in the traffic flow state from free driving state to traffic congestion state.

Full Velocity
Difference model

2001
The model considered the impact on vehicle acceleration when the velocity difference between the front and rear vehicles is positive or negative.

Safe-distance models
Kometani model 1959 The model calculated the safe following distance based on the velocity of the preceding and following vehicles.
Gipps model 1981 The model took into account the safety distance and acceleration constraints of the vehicle.

Fresim model 2003
The model takes the time headway into account from the driver's subjective awareness. The driver's visual information was the main factor affecting the acceleration The early car-following model considered the velocity difference between the front and rear cars and the safety distance separately. The optimal velocity model had more comprehensive considerations and was more widely used in simulations; however, it did not consider the acceleration of the leading car [6]. The full velocity difference model compensated for this shortcoming [7]. The intelligent driver model (IDM) is a car-following model proposed in the context of autonomous driving based on reasonable assumptions [8]. The model showed controllable stability in the simulation of actual highway data. It can produce real acceleration distribution and reasonable behavior in all single-lane traffic situations. Qin et al. [9] used IDM to analyze heterogeneous traffic flow mixed with CACC vehicles, and the verification analysis showed that the model can better describe the driving behavior in the intelligently networked environment. Subsequently, researchers mainly improved the IDM model from the perspective of traditional driving characteristics [10] and autonomous driving technologies [11].
However, these models assume that vehicles are driving on the centerline of the lane and do not fully consider the effect of lateral spacing [12]. Therefore, this paper introduces the personal space (PS) to improve the applicability of the car-following model. PS is a domain concept that refers to the space where psychological tension is caused by the invasion of others. In society, individuals have their own PS [13]. In transportation, Pham et al. [14] applied PS to personal mobility vehicles and established a PS-based avoidance pedestrian assistance system by quantifying the lateral demand distance and forward demand distance of traffic individuals. However, research on applying PS to the analysis of traffic flow is still lacking.
The stability analysis of mixed traffic flow considering PS under the connected and automated environment mainly focuses on two aspects: the improvement of the carfollowing model, and the analysis of the linear stability of the model.
(1) In the HD-CACC mixed traffic flow environment, based on the classical IDM, an IDM model considering PS is established. In the model, it is necessary to specifically consider how the PS is represented in the car-following model, and the difference between the CACC and HD vehicle car-following models.
(2) Traffic flow stability is an important indicator for judging traffic systems against disturbances [15]. Based on the above model, the traffic flow following stability under different parameter configurations is analyzed. Compared with the classic IDM, the effectiveness of PS for maintaining stability needs to be verified. This paper introduces PS into the study of the car-following behavior of individuals in traffic and constructs a car-following model by considering PS in the environment of HD-CACC mixed traffic flow. Further, the stability domain of the car-following model is analyzed through the linear control theory method, and the validity of the model is verified by comparing it with the classic IDM.
This paper is structured as follows. Section 1 is the literature review and problem statement. Section 2, Section 3, and Section 4 cover the test analysis model construction, and stability analysis, respectively. Section 5 discusses the model, and Section 6 verifies the stability and effectiveness through experimental simulation. Finally, Section 7 concludes this paper.

Test Analysis
In this test, Haixia Road in Chongqing, China was chosen as the test section. The road, an urban expressway, was officially opened to traffic in 2012, and the road surface is smooth and free of potholes. It is an important node of the main road in the urban area of Chongqing. The length of the road section is 1200 m, and the width of the road is 43.5 (2.5 + 18.75 + 1.0 + 18.75 + 2.5) m, which is a two-way 10-lane urban road, and a median strip separates the opposing lanes.
An unmanned aerial vehicle [16] was used to conduct a test survey on this expressway. The aerial vehicle model used was a DJI Air 2S. The vertical photography height was 269 m, and the video processing is shown in Figure 1. Based on the Matlab image recognition system detection, the vehicle position coordinates, velocity, acceleration, time headway, and space headway were collected. The test data considered the inner six lanes; lane numbers are shown on the right side of Figure 1, and the data without lane numbers were filtered. In addition, the research object of this paper is only light vehicles. Heavy vehicles such as buses and trucks were not selected in Figure 1. After filtering trajectories without lane numbers and trajectories with lane number changes, a total of 306 valid track-following straight trajectory data were obtained. The relationship between the continuous driving velocity of the vehicle and the distances to the (1) right and (2) front of the vehicle were obtained by fitting, as shown in Figure 2. From Figure 2, it can be observed that there is an evident correlation between the velocity and the distance of the vehicle in the forward and lateral distances. During the driving process of the vehicle, as the space headway to the leading vehicle decreases, the rear vehicle will be stimulated to slow down. Similarly, the closer the vehicle is to the sideways vehicle, the greater the incentive for the driver to slow down. Simultaneously, the desired following distance of the vehicle is intruded, prompting a larger deceleration. There are also differences in the PS of different drivers, resulting in the irregular curved shape observed in Figure 2.
Therefore, the forward and lateral distances of the vehicle are important factors affecting the vehicle's following behavior. To study the influence mechanism of this factor on traffic flow stability, the car-following model is established.

Car-Following IDM Considering PS
PS is variable and regional. Currently, the quantitative calculation method for the shape and distance of PS has not been determined. Hecht et al. [17] proved through the simple stopping distance method that the boundary of PS is approximately circular. Furthermore, the results of anisotropic testing showed that the human anterior thoracic area has a larger demand for space than the back area [18]. This paper focuses on the forward and lateral distances of the vehicle. An egg-shaped boundary is defined based on the above elliptical boundary. The variability is determined by the velocity of the vehicle and individual driver differences. The regionality is determined by the PS boundary formula expressed in Equation (1) below.
where S i (t) is the PS boundary of vehicle i at time t; x i (t) and y i (t) are the coordinates of vehicle i at time t; x and y are points on the PS boundary; A i is the forward distance of vehicle i; and B i is the lateral distance of vehicle i. The specific shapes are shown in Figure 3. Based on the intelligent driver model (IDM) [8], as shown in Equation (2), and by considering the influence of PS on car-following behavior, a PS-IDM analytical model of car-following behavior is established, as shown in Equation (3).
v is the vehicle acceleration; v 0 is the current velocity of the vehicle; a is the maximum acceleration; v f is the free-flow velocity; s 1 is the desired vehicle distance; h is the current space headway of the vehicle; l is the vehicle length; s 0 is the minimum parking distance; T s is the safe time headway; ∆v is the velocity difference from the vehicle ahead; and b is the comfortable deceleration.
Based on literature review and test analysis, PS has excellent coupling with carfollowing model. Considering the effects of A i and B i in Equation (2), the following equation for the car-following behavior was obtained.
In Equation (3), the time headway is an intuitive representation of the longitudinal dimension of the driver's difference, i.e., A i denotes the corrected time headway T = T s + ∆T, and ∆T is the time headway correction value. B i denotes the lateral distance factor term λc i , λ is the lateral distance utility parameter, and c i is the lateral intrusion distance of vehicle i. As shown in Figure 4, the intrusion distance is related to the vehicle position and the current velocity.

Conditions for Stability Judgment
The equilibrium velocity of the traffic flow is defined as v e , and the space headway of the traffic flow is defined as equilibrium h e . By considering only the velocity disturbance v n (t) − v e and the inter-vehicle distance disturbance h n (t) − h e in the traffic flow, the firstorder Taylor expansion of the model at f (h e , v e , 0) was obtained. Subsequently, the IDM traffic flow disturbance transfer function was derived from the Laplace transform of the velocity input signal and the Laplace transform of the velocity output signal.
The disturbance transfer function in the equilibrium state of traffic flow is expressed in Equation (4) as follows.
where S is the Laplace domain, and f ∆v , f h , and f v 0 are the partial differentials of ∆v, h and v 0 , respectively, in Equation (3). Subsequently, f ∆v , f h , and f v 0 are substituted to obtain the disturbance transfer function F(S 1 ) of the car-following model of the HD.
The feedback information of the preceding vehicle is added based on the CACC car-following model, and the acceleration of the preceding vehicle is often ignored in the car-following model. In this study, the acceleration feedback coefficient r(0 ≤ r ≤ 1) of the preceding vehicle was added to derive the disturbance transfer of the CACC car-following model. The function is given by Equation (6).
where r is acceleration feedback coefficient of CACC vehicles. In the mixed traffic flow, F(S 1 ) and F(S 2 ) act together; let S = jw according to the following system stability conditions [7]: where j and w are imaginary numbers and frequencies of frequency domains, |·| ∞ is the maximum amplitude of transfer function frequency domain, p is the penetration rate of the CACC vehicle.
The notations involved in the model are shown in Table 2.

Stability Analysis of Mixed Traffic Flow Based on PS-IDM
We performed a theoretical analysis on the stability of the car-following model and introduced PS-IDM with different values to analyze the stability of homogeneous and heterogeneous traffic flow and explored the influence of PS on the stability of traffic flow.

Model Parameter Calibration
In the model, the following parameter values were assigned according to Milanés and Shladover [11]: a = 1 m/s 2 , v f = 33.3 m/s, l = 5 m, and S 0 = 2 m. The value of b for the HD and CACC vehicles were 2 and 2.8 m/s, respectively. [19]. The time headway correction value ∆T and lateral distance utility parameters λ were calibrated using actual survey data.
When the time headway was greater than 3 s, the front and rear cars were less restrictive. After screening, 292 groups of valid data from the test analysis were extracted, with the time headway being between 0 and 3 s. To obtain the time headway correction value of different driver personalities, the data were clustered and analyzed based on the k-means algorithm. The driving styles of drivers are usually classified into three categories: aggressive, ordinary, and conservative [20]. Therefore, a clustering number of 3 was used. The results for clusters 1, 2, and 3 were 2.2, 1.6, and 1.1, respectively, and are shown in Figure 5. Based on the above analysis, the longitudinal stability study was conducted to explore individual differentiation at the values of 2.2, 1.6, and 1.1 s.
The difference between the maximum abscissas under free flow conditions was screened, and it was assumed that the current driver had obtained enough personal space, and the difference between the vehicle velocity and the abscissa was fitted to obtain the mathematical relationship between v 0 and 2B i . The lateral direction of PS was obtained through further processing. The relationship between distance B i and vehicle velocity v 0 is shown in Equation (8).
where B i is the lateral distance of the vehicle i, i is the vehicle number, and v 0 is the current velocity of the vehicle.

B i was substituted into the second item
obtain the partial PS boundary of vehicle i at x < x i (t), which was used to quantify c i , the lateral intrusion distance of vehicle i, as expressed in Equation (9).
where B i+1 is the PS lateral distance of the laterally adjacent vehicle, k is the width of the vehicle, y i is the ordinate of the vehicle, y i+1 is the ordinate of the laterally adjacent vehicle, and i + 1 is the number of the laterally adjacent vehicle. c i was substituted into Equation (3), and, using the nonlinear regression (nlinfit) and nonlinear curve fitting (lsqcurvefit) functions of Matlab to fit the two values of deceleration b, the HD vehicle lateral distance utility parameter was found to be 0.055. Similarly, the CACC vehicle lateral distance utility parameter was 0.048. The fitting residual error, the sums of squares, were all less than 0.001.
The results of parameter calibration are shown in Table 3.

Homogeneous Traffic Flow
In Equation (7), when p = 0, it degenerates to a homogeneous HDV traffic flow; when p = 1, it degenerates to a homogeneous CACC vehicle traffic flow stability condition.
Based on Equation (7) and ∀w > 0, the stability conditions of the homogeneous HDV traffic flow are given by Equation (10) below.
When the free flow velocity is 0-33.3 m/s, the stability of HDV at different equilibrium velocities at T = 1.1 s, T = 1.6 s, and T = 2.2 s is analyzed by using Equation (10), as shown in Figure 6. Similarly, the stability condition of the traffic flow of homogeneous CACC vehicles should satisfy Equation (11).
The sensitivity analysis of CACC vehicle acceleration feedback coefficients r = 0.1, r = 0.2, and r = 0.3 is shown in Figure 7.   Figure 7a shows that the unstable velocity range gradually decreases in the case of T = 1.1 s with an increase in r. When r = 0.1, the stability of the homogeneous CACC vehicle traffic flow increases with the velocity, reaching the zero point at 23 m/s. The unstable state turns to a stable state, and then it is stable within the free flow velocity. When r = 0.2, the stable zero point of the homogeneous CACC vehicle traffic flow shifts to the left, which is unstable at 0~21 m/s, and stabilizes after the traffic velocity exceeds 21 m/s. When r = 0.3, the unstable velocity boundary continues to shrink, and the zero-boundary point moves to the right to 17.9 m/s. Figure 7b shows that the traffic flow of homogeneous CACC vehicles is stable in the initial stage and then gradually decreases, when r = 0.1, in the case of T = 1.6 s. It reaches the first zero boundary point at 1.6 m/s and becomes unstable. This state continues to reach the second zero boundary point at 22 m/s, and then it is stable until 33.3 m/s; when r = 0.2, the first zero boundary point moves to the right, and the second zero boundary point moves to the left. The unstable velocity range of the CACC vehicle is 3~19.6 m/s; when r = 0.3, the traffic flow is only unstable at 7.5~14.5 m/s. Figure 7c shows that the stability of the homogeneous CACC vehicles' traffic flow gradually decreases within 0~4.6 m/s, and the traffic flow starts to be unstable at 4.6 m/s and stabilizes when the velocity reaches 20.8 m/s, when r = 0.1, in the case of T = 2.2 s; when r = 0.2, the unstable velocity range is reduced, the left boundary is shifted to 7.5 m/s, and the right boundary is shifted to 17.5 m/s; when r = 0.3, homogeneous CACC vehicles traffic flow is stable within the free flow velocity range.

Heterogeneous Traffic Flow
According to the discrimination conditions in Section 4, the discrimination conditions for the stability of heterogeneous traffic flow can be generalized using Equations (10) and (11) to obtain Equation (12) as follows.
where f is the stability of heterogeneous traffic flow, f 1 and f 2 are each the stability of homogeneous traffic flow, and p is the penetration rate of the CACC vehicle. Using a homogeneous traffic flow analysis method to adjust the CACC vehicle penetration rate, the mixed traffic flow stability region under different CACC vehicle penetration rates and different equilibrium velocities can be calculated. As shown in Figure 8, the stable region conditions for four types of r values and three types of T values are given. The black area in the figure is the unstable area, and the white area outside the black area is the stable area. In Figure 8, when the CACC vehicle penetration rate is 0, the boundary value of the black unstable area and the white stable area is the velocity boundary value of the homogeneous HDV traffic flow. It can be obtained from the analysis results in Figure 6. As shown in Figure 8, when the CACC vehicle penetration rate is 1, the velocity boundary and the values of r are negatively correlated, i.e., the larger the value of r, the smaller the unstable velocity range of the CACC vehicle. When r = 0.1, the velocity boundaries are 0~23 m/s, 1.6~22 m/s, and 4.6~20.8 m/s, corresponding to T = 1.1 s, 1.6 s, and 2.2 s, respectively. When r = 0.2, the velocity boundaries are 0~21 m/s, 3~19.6 m/s, and 7.5~17.5 m/s, corresponding to T = 1.1 s, 1.6 s, and 2.2 s, respectively. When r = 0.3, T = 1.1 s and 1.6 s correspond to the velocity boundaries 0~17.9 m/s and 7.5~14.5 m/s, respectively. When the value of T is 2.2 s, the traffic flow of homogeneous CACC vehicles is stable in the free-flow velocity range. When r = 1, the unstable region is the smallest.
Therefore, under the same value of T, as the feedback coefficient r increases, the unstable region of heterogeneous traffic flows shrinks. Furthermore, under the same value of r, with an increase in time headway T, the unstable region of heterogeneous traffic flow is significantly reduced.
For each individual unstable domain, when r ≥ 0.2, the traffic flow instability region of homogeneous CACC vehicles is within the homogeneous HDV traffic flow instability region, i.e., in the homogeneous HDV traffic flow stability region, and the heterogeneous traffic flow of any CACC vehicle permeability is stable. Especially, when r = 0.1, the unstable velocity range of the homogeneous CACC vehicle traffic flow is slightly larger than that of the homogeneous HD traffic flow.
Until r = 0.3, T = 2.2, p ≥ 0.72, a constant stability region appears. When r = 1, the CACC vehicle penetration rate reaches 0.395, 0.24, and 0.16, and the heterogeneous traffic flow is stable within the free flow velocity range of T = 1.1 s, 1.6 s, and 2.2 s, respectively.

Model Validation Analysis
We used a Matlab simulation experiment, compared with the IDM model, to analyze the effect of the PS-IDM model proposed in this paper on the stability of mixed traffic flow, to verify the correctness of the theoretical analysis.
Based on the international standard ISO2631-1 [21] related to human body vibration evaluation and driving comfort standards, the driving comfort during the disturbance propagation process of the car-following model in this study was evaluated, and the traffic flow stability was quantitatively evaluated. The driving comfort index can be calculated using Equation (13) as follows.
where CI is the driving comfort index, a i is the acceleration value of each vehicle obtained by statistics, and m is the total acceleration statistic. This index uses the root mean square value of the total vibration acceleration to determine the comfort of the occupants. The smaller the value of CI, the better the driving comfort, and vice versa, which verifies the consistency of the analytical results with the stability theory. Secondly, the intrusion rate index was used to evaluate the impact of lateral intrusion of adjacent vehicles into the PS, and to explore the internal relationship between the intrusion rate and longitudinal stability. Further, this index was used to study the lateral stability of the PS. The formula for calculating the intrusion rate is expressed in Equation (14).
where I is the percentage index of the invasion rate, c is the invasion length, and l a is the length of the PS in the invasion direction. l a is equal to the lateral distance of the vehicle PS B i . Therefore, Equation (14) can be rewritten as: The smaller the value of I, the smaller the PS intrusion rate, and the lower the effect on the vehicle PS. Conversely, the higher the intrusion rate, the greater the impact.

Experimental Design
A MATLAB simulation experiment was used to analyze the stability state of the mixed traffic flow of the disturbance transfer function under different CACC vehicle permeabilities. According to the analysis results of the stability of heterogeneous traffic flow in Section 5.2.2, the acceleration feedback coefficient of preceding vehicle r is set to 0.3, and the experimental results are satisfactory and without loss of generality. Based on MATLAB numerical simulation, the effectiveness of the PS-IDM proposed in this study is compared with that of a single IDM.
In the simulation experiment, the number of HD and CACC vehicles was determined using the set p value, and 50 vehicles were randomly placed on a 200 m single-lane road. Initially, the hybrid platoon drove at the equilibrium velocity v e , and then a disturbance was generated by randomly slowing down to break the equilibrium state. The disturbance propagates in pairs in the platoon, based on the disturbance transfer function. In the propagation process, the operating status of each vehicle was obtained, and the acceleration and distance between the heads of the vehicle were extracted to calculate the relevant indicators. According to the above stability analysis results, HD vehicles are always unstable at 4.4~20 m/s. To analyze the influence of CACC vehicles on the instability of HDV traffic flow, the initial equilibrium velocity of the fleet was selected as 10 m/s for the numerical simulation. The random slowdown probability was set to 0.3, the simulation step length was set to 0.1 s, and comfort was selected as the evaluation index of the simulation experiment.

Simulation Results
The results of the numerical simulation experiments are listed in Table 4. Taking the r value of 0.3, the comfort index was calculated, and the results are shown in the tables below. It can be noted from Table 4 and Figure 9 that the higher the penetration rate of CACC vehicles, the higher the driving comfort under the same T value.
Compared with the homogeneous CACC vehicle traffic flow, the CI value of the homogeneous HD traffic flow was reduced by 64.67, 66.06, and 72.16% when T = 1.1, 1.6, and 2.2 s, respectively. The results showed that an increase in PS improved driving comfort, which was consistent with the analytical results of the longitudinal stability theory. Figure 10 shows a comparison chart of the comfort index during the simulation process of the PS-IDM and single IDM models. The solid line is the PS-IDM comfort index fitting curve, and the dashed line is the single IDM comfort index fitting curve. Under the same circumstances, the CI value of a single IDM is greater than the CI value of the PS-IDM. Compared with a single IDM, the comfort of the PS-IDM can be increased by up to 20.7%.  The intrusion rate of each vehicle was discrete data, and the simulation experiment was performed using MATLAB. In addition, the average value of the intrusion rate of 50 for the simulated vehicles was taken as the simulation result of the intrusion rate index. Taking T = 1.1 s and r = 0.3 as an example, the result is shown in Figure 11.
In Figure 11, when the HD vehicle traffic flow is homogeneous, the average intrusion rate is higher than 24%. As the penetration rate of CACC vehicles increases, the intrusion rate decreases, when the traffic flow of the homogenous HD vehicle changes to the traffic flow of homogenous CACC vehicles. In particular, at the beginning of the transition period (p ∈ (0, 0.3)), the downward trend is evident. When the mixed traffic flow degenerates into a homogenous CACC vehicle traffic flow (p = 1), the average intrusion rate drops to 19.4%.
The simulation results show that there is a positive correlation between the intrusion rate index and the comfort index. The smaller the value of the comfort index, the smaller the intrusion rate index value, and the better the driving stability, which is consistent with the analytical results of the lateral stability theory.

Conclusions
(1) Aiming at the HD-CACC mixed traffic flow, this study establishes a car-following PS-IDM model based on an intelligent driver model by considering personal space. Using theoretical analysis and MATLAB simulation experiments, the stability regions of homogeneous traffic flows under different PS and heterogeneous traffic flows with different CACC vehicle permeabilities were analyzed. The results show that the PS-IDM model can improve the driving comfort by up to 20.7% compared with a single IDM model, and the average PS intrusion rate is reduced by 4.6%, which verifies the effectiveness of the PS-IDM model.
(2) The analytical results of the traffic flow stability of car-following behavior, considering PS, indicates that the unstable velocity range of homogeneous or heterogeneous traffic flows decreases with an increase in PS. In addition, the unstable velocity range of the CACC vehicle traffic flow has a negative correlation with the acceleration feedback coefficient. The increased penetration rate of CACC vehicles with the acceleration feedback coefficient below 0.2 could not turn an unstable traffic flow into a stable traffic flow; when the acceleration feedback coefficient exceeded 0.2, the increase in CACC vehicle permeability can effectively transform the stability of traffic flow. Regardless of the value of the acceleration feedback coefficient, increasing the time headway can make the traffic flow to a stable state. Therefore, an increase in the longitudinal distance of the vehicle PS is beneficial to the stability of the traffic flow.
(3) The results of the numerical simulation experiments were consistent with the analytical results of the stability theory. The model and experiment in this study made assumptions regarding the lane conditions. In the future, we will analyze the lateral factors influencing the multi-lane situation and conduct in-depth research on the lateral distance. In addition, because different values of the acceleration feedback coefficient can correspond to HD, multi-level individual intelligent vehicles, and multi-level connected automated vehicles, further analysis and research on the stability of complex heterogeneous traffic flows will be conducted.