Cooperative Adaptive Cruise Algorithm Based on Trajectory Prediction for Driverless Buses

: Cooperative adaptive cruise control (CACC) technology offers a proven solution to the current trafﬁc congestion problems caused by the yearly growth of car ownership. Coping with random lane changes of bypass vehicles under the condition of trafﬁc congestion is a challenge for urban driverless vehicles. In this paper, to meet the demand for high comfort driverless buses driving on urban roads, an active anti-disturbance following control method for driverless buses based on bystander vehicle intention recognition and trajectory prediction is proposed for the scenario of bystander vehicle cut-in during driving to alleviate the disturbance caused by bystander vehicles, improve passenger comfort, and suppress multi-vehicle oscillation. The simulation results show that the intelligent prediction system-based queue reduces the trafﬁc oscillation rate by an average of 9.8% and improves the comfort level by an average of 11% under side-car insertion conditions. The results of the real vehicle test show that the vehicles based on the intelligent prediction algorithm have a 25.5% reduction in maximum speed adjustment, 14.5 m average reduction in following distance, 6% improvement in comfort, and 27% improvement in rear vehicle comfort.


Introduction
Most CACC algorithms consider only the traffic information at the current moment as input, and their performance may be limited under some complex operating conditions, for example, the side lane vehicles cut into the working conditions. It is thus highly desirable to use predictive algorithms that predict future driving conditions to achieve smoother speed control, thereby reducing traffic oscillations and passenger line comfort. Therefore, on this basis, this paper addresses the cut-in working conditions of the by-pass vehicles, by predicting the lane change intentions and trajectories of the by-pass vehicles, and flexing the following distance and speed in advance in order to improve passenger comfort and reduce the distance.

Related Research
CACC is the use of V2V communication technology in combination with the ACC system to ensure stable following of the queueing vehicles. Compared to ACC, CACC has more powerful collision avoidance and maximizes traffic flow throughput. CACC can improve traffic flow by reducing headway time distance and reducing traffic disturbance.
In the queueing truck for the CACC development at Auburn University [1], when a cut-in is detected, the CACC switches the target from the lead truck to the cut-in vehicle in order to create a larger spacing between the following truck and the cut-in vehicle. This process results in reduced fuel savings because it requires more work from the longitudinal controller, which is dangerous if the clearance is not wide enough when the vehicle cuts in, especially considering that an 80,000 lb. trailer requires 400 to 500 feet of stopping distance idea of reinforcement learning (RL), which does not rely on the accurate modeling of bus vehicle dynamics. Considering a time-varying topology, each autonomous vehicle received information only from previous vehicles within its communication range and learned the distributed controller in real time from online time-distance, velocity, and acceleration data collected from the system trajectory. Results showed that dedicated bus lane travel times were close to current travel times even with a 30% increase in transit capacity. Jurj et al. [12] trained the reviewer (SAC) reinforcement learning (RL) algorithm, which used physical knowledge, such as distances, to avoid interference, to automatically adjust the ideal longitudinal distance between the self and the leading vehicle. The physically guided (PG) RL method was better to avoid collisions in any selective deceleration condition and any queueing condition, and the method was more reliable.
The boundaries of several types of the control methods analyzed above are located in the controlled vehicle. Furthermore, with the advancement of 5G communication technology in recent years, shop floor communication and vehicle infrastructure communication have attracted increasing attention from researchers as a way to address the limitation of too short ACC fields of view.
Ploeg et al. [13] funded by the Dutch Ministry of Economic Affairs through the High-Tech Automotive Systems (HTAS) project, from the California Advanced Transportation and Freeway Partnership Program at the University of California, Berkeley, Institute of Transportation Studies, addressed the situation arising from special working conditions through the popularization of devices and the application of 5G communication. Such methods can guarantee short-range following, as well as string stability. Bosch, together with Huawei and Vodafone, the U.K. mobile network operator, demonstrated that direct, delayfree data exchange between vehicles can also improve the functionality of driver assistance systems such as ACC [14][15][16].
The contributions of this paper are summarized as follows: (1) Based on the convolutional social pool trajectory prediction architecture, the driver personality factor algorithm and the model and sampling short-term optimization algorithm are used as a means to improve the distribution of future trajectory points during the time domain of autonomous driving prediction and to reduce the trajectory error of vehicle prediction.
(2) A single/multi-vehicle intelligent predictive adaptive cruise control system based on an immunity controller, a softening corrector, and a flexible fading corrector is used and outperforms the other method by simulation and experimental results in real vehicles.
The rest of this paper is structured as follows. In Section 2, a multi-vehicle collaborative following framework based on by-pass lane prediction is proposed and each module is explained in detail in different subsections. Test results and analysis of the proposed approach are discussed in Section 3. Finally, Section 4 concludes the paper. Figure 1 shows the framework of the cooperative adaptive cruise system with intelligent prediction proposed in this paper, including the prediction module, flexible following factor corrector, self-anti-disturbance controller, and flexible fading corrector.

Materials and Methods
corporates both a vehicle dynamics model and a motor-battery model. In the figure, a des is the expected acceleration, v act−bus1 is the actual velocity of vehicle 1, v act−busi is the actual velocity of vehicle i, and T des and T br_des are the desired driving force and desired braking force of the vehicle.
The flexible attenuation corrector is a longitudinal following pattern correction for vehicles in the rear in response to bystander vehicles cutting into the open. ρ flex−damp is the flexible attenuation following factor in the speed control framework.  The prediction module combines two algorithms: data-driven trajectory prediction and prototype trajectory prediction algorithms. The data-driven algorithm includes two modules: the first for driver unpredictable elements and the second for social long-term and short-term memory neural network prediction, as shown in the upper part of Figure 1. The trajectory prototype-driven method contains a module for kinematics, trajectory matching, and trajectory sampling. In the figure social-LSTM is a long and short-term neural network considering social interaction properties, CTRA-model is a constant-turn-rate and acceleration model, and V2V is the communication between vehicles.
The flexible following factor corrector is the distribution degree of the road centerline which is then employed as the lane change overlap degree, and the flexible following factor and flexible switching factor of the following target are produced by a double second-order time-varying filter for flexible speed tracking after target switching. ρ flex is the flexible following factor in the speed control framework, and τ flex flexibly switches the target switching factor _object at the time of target switching, preventing target switching when the switching factor is below a certain threshold. d obj and v obj are the speed and relative distance of the target vehicle.
The self-anti-disturbance controller, an instantaneous observation approach based on the expansion state observer, is proposed to correct for and suppress the internal disturbance generated by changes in the internal vehicle parameters and random disturbances induced by changes in the external road environment. The autonomous bus is finally directed to follow the target vehicle using an inverse longitudinal dynamics model that incorporates both a vehicle dynamics model and a motor-battery model. In the figure, a des is the expected acceleration, v act−bus1 is the actual velocity of vehicle 1, v act−busi is the actual velocity of vehicle i, and T des and T br_des are the desired driving force and desired braking force of the vehicle.
The flexible attenuation corrector is a longitudinal following pattern correction for vehicles in the rear in response to bystander vehicles cutting into the open. ρ flex−damp is the flexible attenuation following factor in the speed control framework.

Fusion Social-LSTM and Trajectory Prototype Prediction Algorithm
The prediction module is a vehicle trajectory dynamic prediction algorithm based on driving intention recognition, which provides the basis for the ACC and CACC of driverless buses during the process of vehicles cutting into the home lane from the side lane. First, a recurrent neural network model (Social-LSTM) based on social long-term and short-term memory is used as the base model of the trajectory dynamic prediction algorithm; then, a trajectory prototype is generated by setting a minimum value function to match the vehicle model with the planned sampled trajectory, and the trajectory prototype is used to modify the base prediction model to improve the correct rate of the short-term lane change trajectory prediction of the bypass vehicle. Finally, the driver adventurous factor is introduced into the driver's lane change behavior character to further modify the base model and improve the correct rate of long-term lane change trajectory prediction of bypass vehicles. x(t) ..
Traj sampling (i) = S(t) where τ character is the driver aggression factor, P lanechange− f usion is the social tensor input from the social pool layer, and P predicted−driver is the optimized social tensor. T is the headway time distance, and ω is the headway time distance as a percentage of the negative influence on the influence on the driver aggression factor. Lane ID is the lane type, β is the percentage of headway time distance to the positive influencing factor of driver aggressiveness. Agg is the aggressiveness due to acceleration magnitude, and α is the percentage of acceleration/deceleration to the positive influencing factor of driver aggressiveness.
Wid vehicle is the model size, and γ is the percentage of the model with a negative influence on the driver adventurous factor. The α, β, ω, and γ parameters are assigned according to the probability of lane change. And u x , u y , σ x , σ y , ρ are the five two-dimensional Gaussian distribution parameters based on the data-driven output. T intention is the trajectory of the data-driven output, Traj model is the kinematic model, Traj sampling is the trajectory sampling model, Traj matching is the matching function of the kinematic model and trajectory sampling, T f usion is the fused trajectory, and ω model and ω intention are the trajectory prediction weights of different methods.

Flexible following Factor Corrector
On the basis of the anti-turbulence controller for speed following control, the flexibility corrector ensures the comfort of the driverless bus by softening the speed of the side lane vehicle with lane change intention and inserting it into the process of this lane. The overlap degree under vehicle prediction, due to its additional spikes, cannot be directly used in the following mode; therefore, the flexible following factor is proposed. To calculate the flexible following factor, a dual second-order time-varying filter is selected. In the signal processing, the digital dual fourth-order filter is a second-order recursive linear filter containing two poles and two zeros.
where p lane (x|µ, Σ) is the calculated value of the two-dimensional Gaussian distribution points under different positions of the vehicle, Overlap is the overlap degree, ρ f lex and τ object are the flexible following factor and target switching factor, and a 0 , a 1 , a 2 , b 0 , b 1 , and b 2 are the curve change parameters, which are used as different parameters to ensure that the curve presents the change of the Gaussian curve, as well as the incremental change.

Self-Tampering Controller
The core control of speed control is the target vehicle speed following, and this subsection mainly establishes the control-oriented power system model and the inverse model based on the whole vehicle power transmission system. Among them, the whole vehicle dynamics model mainly includes the power source model, the power transmission system model, and the whole vehicle longitudinal dynamics model. In the conventional ACC system, the output of the lower controller in the hierarchical control, i.e., the input of the actuator, is the drive pedal opening and brake pedal opening transformed by the motor torque and brake pressure; thus, the role of the lower controller is to transform the output of the upper controller, i.e., the desired acceleration, into the input parameters of the actuator. In order to clearly express the input and output variables of the observer, the improved observer equation is as follows: where u c = [u y] T is the input variable of the observer, including the control quantity output by the control algorithm and the state variable output by the system. y c is the output of the observer, including the vehicle speed v, the first-order derivative of the vehicle v, and the estimate of the disturbance f . After parameterization, the poles of the characteristic equation can be placed at the same position (−ω o ), and ω o is the observer bandwidth, i.e., the gain matrix L = [3ω 0 3ω 0 2 ω 0 3 ] T of the observer is used.

Flexible Decay Algorithm
As shown in Figure 2, according to the analysis of the actual multi-vehicle following scenario, the impact of the speed change of the lead vehicle on the rear vehicle during the queue following process presents a weaker and weaker situation; therefore, the flexible attenuation factor of the queue following in this paper also presents a trend of smaller and smaller changes, i.e., the attenuation becomes smaller on the basis of the flexible factor of the lead vehicle to avoid the impact of the cut-in condition on the rear vehicle passenger comfort.
Machines 2022, 10, x FOR PEER REVIEW 7 of 30 queue following process presents a weaker and weaker situation; therefore, the flexible attenuation factor of the queue following in this paper also presents a trend of smaller and smaller changes, i.e., the attenuation becomes smaller on the basis of the flexible factor of the lead vehicle to avoid the impact of the cut-in condition on the rear vehicle passenger comfort. The variation in the state equation transformation function of the dual second-order time-varying filter is used to calculate the following flexible attenuation factor, while multiple spikes still exist in the late stage of the constructed function in the channel change, which cannot be logically switched; the factor is constructed as follows: where and are the attenuation factors of the fleet; such parameters need to be constructed for different working conditions to verify their influence factors. To explore the laws of the two factors, different vehicle serial numbers and different speed ratios are constructed to explore the influence factors. The speed ratio is 0.83 for case 1 with a bypass vehicle cut-in speed of 10 km/h and a queue vehicle speed of 14 km/h, and 0.85 for case 2 with a bypass vehicle cut-in speed of 12 km/h and a queue vehicle speed of 14 km/h. The results of the three simulations in Figure 3 show that the value of the attenuation factor needs to progressively decrease as the vehicle serial number increases, i.e., the impact of the cut-in vehicle on the vehicle influence becomes smaller and smaller. In Figure  4, the simulation results of different ranges of attenuation factors at different vehicle speeds show that, as the speed of the cut-in vehicle in the side lane approaches the speed ratio of the driverless bus queue vehicles, the value of the attenuation factor progressively decreases, i.e., if the speed of the cut-in vehicle is close to the speed of this lane, the attenuation factor is smaller. The variation in the state equation transformation function of the dual second-order time-varying filter is used to calculate the following flexible attenuation factor, while multiple spikes still exist in the late stage of the constructed function in the channel change, which cannot be logically switched; the factor is constructed as follows: where α i and β i are the attenuation factors of the fleet; such parameters need to be constructed for different working conditions to verify their influence factors. To explore the laws of the two factors, different vehicle serial numbers and different speed ratios are constructed to explore the influence factors. The speed ratio is 0.83 for case 1 with a by-pass vehicle cut-in speed of 10 km/h and a queue vehicle speed of 14 km/h, and 0.85 for case 2 with a bypass vehicle cut-in speed of 12 km/h and a queue vehicle speed of 14 km/h. The results of the three simulations in Figure 3 show that the value of the attenuation factor needs to progressively decrease as the vehicle serial number increases, i.e., the impact of the cut-in vehicle on the vehicle influence becomes smaller and smaller. In Figure 4, the simulation results of different ranges of attenuation factors at different vehicle speeds show that, as the speed of the cut-in vehicle in the side lane approaches the speed ratio of the driverless bus queue vehicles, the value of the attenuation factor progressively decreases, i.e., if the speed of the cut-in vehicle is close to the speed of this lane, the attenuation factor is smaller.   The results of the above simulation show that the attenuation factor is related to the cut-in speed of the vehicle cutting into the side lane and the proximity to the speed of the train convoy, as well as the number of the team train (the head train number is 1, with numbers increasing in order); hence, the following function can be constructed: where is the queue train serial number, is the decay factor of queue train i, is the decay ratio (linearly increasing with the speed difference), is the speed of the sidecar, and is the speed of queue train i. Regarding safety, this paper defines a safety cost to evaluate the safety of the strategy, and the safety evaluation of vehicle at moment is shown in the following equation: where ( ) is the headway time distance of vehicle k at moment i. This represents the time interval between the head ends of two consecutive vehicles passing through a The results of the above simulation show that the attenuation factor is related to the cut-in speed of the vehicle cutting into the side lane and the proximity to the speed of the train convoy, as well as the number of the team train (the head train number is 1, with numbers increasing in order); hence, the following function can be constructed: where i is the queue train serial number, ε i is the decay factor of queue train i, µ is the decay ratio (linearly increasing with the speed difference), v side is the speed of the side-car, and v i is the speed of queue train i.
Regarding safety, this paper defines a safety cost to evaluate the safety of the strategy, and the safety evaluation of vehicle i at moment k is shown in the following equation: (15) where HWT i (k) is the headway time distance of vehicle k at moment i. This represents the time interval between the head ends of two consecutive vehicles passing through a certain section in a queue of vehicles traveling in the same lane. Regarding comfort, adopts jerk as the evaluation criterion of comfort, and its smaller value indicates the higher comfort level of the vehicle. However, for the queue vehicles, the following formula can be used to evaluate the comfort of the remaining vehicles in the queue in reverse.
where a i denotes the acceleration of the i vehicle of the fleet, and a 1 denotes the acceleration of the head vehicle of the fleet. Specifically, if d p,i is smaller, the value of a i is smaller than that of a 1 , i.e., its comfort level is higher than that of the head vehicle. When the head vehicle improves its comfort level through the prediction module, the value of d p,i of the rear queue vehicles becomes smaller and smaller, which means that the comfort level of the rear vehicles becomes higher and higher.

Experiment and Analysis of Results
This section contains a detailed description of the simulation platform and the real vehicle platform used to test the above method as well as the test results and analysis by setting up different working conditions.

Experimental Platform
The experimental platform used in this paper includes an in-loop simulation platform and a real vehicle test platform, where the in-loop simulation platform uses scenario-building software and strategy-building software for joint simulation. The real vehicle platform uses multiple driverless buses driving together in a closed environment for testing and analysis.

Simulation Experiment Platform
To verify the effectiveness of the algorithm in this chapter, a joint simulation platform of prediction-CACC (intelligent prediction for queueing cruise) system based on SIMULINK/ Prescan is designed. The joint simulation platform consists of four main parts: the Prescanbased vehicle dynamics model, information parsing module, information output module, and prediction-CACC predictor and prediction-CACC upper controller module and prediction-CACC lower controller module.
In order to complete the cut-in simulation scenario, a straight experimental road is specially designed, with a total length of 1400 m and a flat road type with no undulations and slopes on the road surface. Both the driverless bus prototype and the team train are small buses with the same parameters, and the cut-in vehicle is a small vehicle whose performance is not considered.
In order to verify the following effect of the single vehicle, two scenarios are designed in the joint simulation, which are the side lane vehicle cut-in for straight road close following and the side lane vehicle cut-in for straight road medium following. In addition, in order to verify the control effect of the queue train, multiple scenarios of queue following conditions are designed, and the conditions are similar to the single vehicle conditions for verifying the following safety and control ability of the flexible decay following factor.

Realistic Autonomous Bus Experiment Platform
The main body of this study relies on a purely electric new energy minibus, which has a complete vehicle wiring control system, radar and vision sensor configuration, and intelligent computing platform, providing a solid foundation and verification platform for the study, as shown in Figure 5. In order to complete the CACC system and the working condition design for detecting vehicles in the side lane in this paper, two additional vehicles are configured on the basis of the experimental prototype. The team train relies on a purely electric new energy bus, which is also equipped with a complete vehicle wire control system, vision sensor system, intelligent technology platform, and workshop communication equipment. To meet the needs of the research scenario, the cut-in vehicle uses the BYD SUV, which is equipped with a 2.0 L direct injection four-cylinder gasoline engine. The main body of this study relies on a purely electric new energy minibus, which has a complete vehicle wiring control system, radar and vision sensor configuration, and intelligent computing platform, providing a solid foundation and verification platform for the study, as shown in Figure 5. In order to complete the CACC system and the working condition design for detecting vehicles in the side lane in this paper, two additional vehicles are configured on the basis of the experimental prototype. The team train relies on a purely electric new energy bus, which is also equipped with a complete vehicle wire control system, vision sensor system, intelligent technology platform, and workshop communication equipment. To meet the needs of the research scenario, the cut-in vehicle uses the BYD SUV, which is equipped with a 2.0 L direct injection four-cylinder gasoline engine.

Performance Verification and Results Analysis
In order to verify the effectiveness of the proposed method, the proposed method is validated by designing a case study for different working conditions and following the vehicle. The validation scheme includes the determination of parameters in the analysis of the flexible decay function, the effect of multi-vehicle following in different operating conditions in the ring simulation, and the effectiveness of stable multi-vehicle following in the actual multi-vehicle following process.

Simulation and Verification of Flexible Fading Algorithm
To verify the rationality and effectiveness of the flexible fading algorithm, in the simulation platform, the cut-in speeds of the cut-in vehicles are set to 10 km/h and 12 km/h, and the queue speed is set to 14 km/h; meanwhile, µ is set to 1-5. The headway time distance of the second vehicle and the jerk are used to evaluate the value range of u. The second following vehicle is used as the evaluation vehicle, whose headway time distance is used to assess its safety, while the jerk is used to assess its comfort, as shown in Figures 6 and 7.

Performance Verification and Results Analysis
In order to verify the effectiveness of the proposed method, the proposed method is validated by designing a case study for different working conditions and following the vehicle. The validation scheme includes the determination of parameters in the analysis of the flexible decay function, the effect of multi-vehicle following in different operating conditions in the ring simulation, and the effectiveness of stable multi-vehicle following in the actual multi-vehicle following process.

Simulation and Verification of Flexible Fading Algorithm
To verify the rationality and effectiveness of the flexible fading algorithm, in the simulation platform, the cut-in speeds of the cut-in vehicles are set to 10 km/h and 12 km/h, and the queue speed is set to 14 km/h; meanwhile, μ is set to 1-5. The headway time distance of the second vehicle and the jerk are used to evaluate the value range of u. The second following vehicle is used as the evaluation vehicle, whose headway time distance is used to assess its safety, while the jerk is used to assess its comfort, as shown in Figures  6 and 7.

Simulation Verification of Different Working Conditions
This subsection sets up the CACC algorithm with or without the flexible following factor and flexible attenuation factor optimization under four kinds of bypass lane insertion conditions, and the guiding vehicle in the four cases is set to run at a fixed speed. The initial speed of the rear vehicle is lower than the initial speed of the head vehicle; thus, there is more obvious acceleration following behavior in the initial stage, and this setting also corresponds to the actual real vehicle test. For the above four scenarios, this subsection configures five queue vehicles with CACC control for the overall test. This subsection mainly analyzes and compares the queues of the two types of fleets during the bypass vehicle cutover, and uses the following distance to assess the compactness of the two types of fleets, the headway time distance to assess the following safety, and the jerk to assess the comfort. For notational convenience, the fleet without flexible following factor and flexible decay factor is referred to as CACCNF in this paper (no factor CACC), and fleets with flexible following factor and flexible attenuation factor are CAVF (factor connected and automated vehicle).In the test effect picture, black, blue, coffee, dark green, purple

Simulation Verification of Different Working Conditions
This subsection sets up the CACC algorithm with or without the flexible following factor and flexible attenuation factor optimization under four kinds of bypass lane insertion conditions, and the guiding vehicle in the four cases is set to run at a fixed speed. The initial speed of the rear vehicle is lower than the initial speed of the head vehicle; thus, there is more obvious acceleration following behavior in the initial stage, and this setting also corresponds to the actual real vehicle test. For the above four scenarios, this subsection configures five queue vehicles with CACC control for the overall test. This subsection mainly analyzes and compares the queues of the two types of fleets during the bypass vehicle cutover, and uses the following distance to assess the compactness of the two types of fleets, the headway time distance to assess the following safety, and the jerk to assess the comfort. For notational convenience, the fleet without flexible following factor and flexible decay factor is referred to as CACCNF in this paper (no factor CACC), and fleets with flexible following factor and flexible attenuation factor are CAVF (factor connected and automated vehicle).In the test effect picture, black, blue, coffee, dark green, purple and light green indicate the first car to the sixth car respectively. Also, in the results of the comparison between the two teams, the dashed line indicates the CACCNF team and the solid line indicates the CAVF team.
Scenario 1: The vehicle in the side lane cuts into this lane with a constant speed lower than the fleet travel speed of 2 m/s; as the vehicle speed in the side lane is lower than 2 m/s in this lane, i.e., the average speed is 7.2 km/h lower than that of the team train, this speed difference is large, and the limit following the working condition may occur for the train. Therefore, the goal of scenario 1 is to verify whether the algorithm proposed in this paper is safe. The simulation results of the two fleets are shown in Figure 8. Due to the trajectory of the vehicles in the bypass lane being predicted in advance, the CAVF fleet adjusts its speed 3.5 s earlier. The statistical analysis based on the simulation results of Figure 8 is shown in Table  1. The overall HWT regulation of the convoy is reduced by 65% on average, while the average reduction in the convoy following distance is 49.6%. Therefore, in scenario 1, the CAVF fleet has a smaller and more compact spacing, thereby ensuring safety. By plotting the velocity and acceleration heatmaps for each vehicle in the CAVF fleet and CACCNF fleet, with time on the horizontal axis and the position of the vehicle on the vertical axis, and by calculating the jerk, the comfort level of the vehicle can be evaluated. As shown in Figure 9. The statistical analysis based on the simulation results of Figure 8 is shown in Table 1. The overall HWT regulation of the convoy is reduced by 65% on average, while the average reduction in the convoy following distance is 49.6%. Therefore, in scenario 1, the CAVF fleet has a smaller and more compact spacing, thereby ensuring safety. By plotting the velocity and acceleration heatmaps for each vehicle in the CAVF fleet and CACCNF fleet, with time on the horizontal axis and the position of the vehicle on the vertical axis, and by calculating the jerk, the comfort level of the vehicle can be evaluated. As shown in Figure 9. Based on the results shown in the figure, the following conclusions can be drawn: the acceleration changes in the CACCNF fleet are concentrated from 10 s to 35 s and the acceleration changes are more dispersed, while the acceleration changes of the CAVF fleet are concentrated from 5 to 20 s; after 25 s, the acceleration changes are more concentrated with a smaller change rate and higher comfort. As shown in Figure 10. The average CACCNF fleet jerk is 54%, and the average CAVF fleet jerk is 48%, representing a 6% increase in comfort. The comfort increase is smaller due to the safety constraints, as shown in Table 2. Based on the results shown in the figure, the following conclusions can be drawn: the acceleration changes in the CACCNF fleet are concentrated from 10 s to 35 s and the acceleration changes are more dispersed, while the acceleration changes of the CAVF fleet are concentrated from 5 to 20 s; after 25 s, the acceleration changes are more concentrated with a smaller change rate and higher comfort. As shown in Figure 10. Based on the results shown in the figure, the following conclusions can be drawn: the acceleration changes in the CACCNF fleet are concentrated from 10 s to 35 s and the acceleration changes are more dispersed, while the acceleration changes of the CAVF fleet are concentrated from 5 to 20 s; after 25 s, the acceleration changes are more concentrated with a smaller change rate and higher comfort. As shown in Figure 10. The average CACCNF fleet jerk is 54%, and the average CAVF fleet jerk is 48%, representing a 6% increase in comfort. The comfort increase is smaller due to the safety constraints, as shown in Table 2. The average CACCNF fleet jerk is 54%, and the average CAVF fleet jerk is 48%, representing a 6% increase in comfort. The comfort increase is smaller due to the safety constraints, as shown in Table 2. Scenario 2: The side lane vehicles cut into the lane at a higher speed than the fleet travel speed (constant at 2 m/s); because the side car speed insertion is higher than the team train speed (constant at 2 m/s), in the process of cutting in, due to the safety and comfort requirements, the train in the early stage suffers a certain degree of impact, before later converging to the team cruise speed, as shown in Figure 11.  Scenario 2: The side lane vehicles cut into the lane at a higher speed than the fleet travel speed (constant at 2 m/s); because the side car speed insertion is higher than the team train speed (constant at 2 m/s), in the process of cutting in, due to the safety and comfort requirements, the train in the early stage suffers a certain degree of impact, before later converging to the team cruise speed, as shown in Figure 11. The statistical analysis based on the simulation results of Figure 11 is shown in Table  3. The overall HWT regulation of the convoy is reduced by 65% on average, Meanwhile, according to Table 4, the average reduction in convoy following distance is 3.3%. Therefore, in scenario 1, the CAVF fleet has a smaller and more compact spacing, thereby ensuring safety.  The statistical analysis based on the simulation results of Figure 11 is shown in Table 3. The overall HWT regulation of the convoy is reduced by 65% on average, Meanwhile, according to Table 4, the average reduction in convoy following distance is 3.3%. Therefore, in scenario 1, the CAVF fleet has a smaller and more compact spacing, thereby ensuring safety. By plotting the velocity and acceleration heatmaps for each vehicle of the CAVF fleet and CACCNF fleet, with time as the horizontal axis and the position of the vehicle as the vertical axis, and by calculating the jerk, the comfort level of the vehicle can be evaluated. From the results displayed in Figure 12, it can be concluded that the speed of the CAVF fleet and CACCNF fleet varies more in the time domain of 5-15 s for CAVF, but the results also guarantee a tighter fleet in this scenario. As shown in Figure 13, the acceleration distribution of this condition is not obvious; thus, this condition is not compared. The side lane vehicle cuts into this lane, driving with a positive variable speed and an acceleration higher than the initial speed of the team train. As shown in Figure 14, in the CACCNF fleet, because the side lane vehicle accelerates into this lane and the speed is higher than the speed of the team train, the traffic oscillation changes gradually from the head vehicle to the tail vehicle, and the deceleration process of the head vehicle, as well as the speed change process, is longer, which leads to the back vehicles undergoing real-time speed changes to protect the stability of the vehicle. The time domain of speed change of the CAVF fleet is reduced by 2 s, and the fleet stability is reached more quickly. The side lane vehicle cuts into this lane, driving with a positive variable speed and an acceleration higher than the initial speed of the team train. As shown in Figure 14, in the CACCNF fleet, because the side lane vehicle accelerates into this lane and the speed is higher than the speed of the team train, the traffic oscillation changes gradually from the head vehicle to the tail vehicle, and the deceleration process of the head vehicle, as well as the speed change process, is longer, which leads to the back vehicles undergoing realtime speed changes to protect the stability of the vehicle. The time domain of speed change of the CAVF fleet is reduced by 2 s, and the fleet stability is reached more quickly. The statistical analysis based on the simulation results of Figure 14 is shown in Table  5. The HWT adjustment time of the two types of queues reveals that the adjustment time of the CAVF fleet is shorter compared to that of the CACCNF fleet by about 4 s, and the results of the following distance shown in the above figure indicate that the following distance of the CAVF fleet is lower and the spacing is smaller. As shown below, the heatmaps of velocity and acceleration for each vehicle of the CAVF fleet and CACCNF fleet are plotted separately, with time as the horizontal axis and the position of the vehicle as the vertical axis. From the results shown in Figures 15 and  16, it can be concluded that the velocity density of the CAVF fleet is concentrated in the range of 5-16 s, and the velocity density of the CACCNF fleet is concentrated in the range of 10-26 s. The velocity oscillation time of the CACCNF fleet is longer and the velocity fluctuation is higher compared with the CAVF fleet. Under this condition, the average fluctuation of both fleets is basically the same. The statistical analysis based on the simulation results of Figure 14 is shown in Table 5. The HWT adjustment time of the two types of queues reveals that the adjustment time of the CAVF fleet is shorter compared to that of the CACCNF fleet by about 4 s, and the results of the following distance shown in the above figure indicate that the following distance of the CAVF fleet is lower and the spacing is smaller. As shown below, the heatmaps of velocity and acceleration for each vehicle of the CAVF fleet and CACCNF fleet are plotted separately, with time as the horizontal axis and the position of the vehicle as the vertical axis. From the results shown in Figures 15 and 16, it can be concluded that the velocity density of the CAVF fleet is concentrated in the range of 5-16 s, and the velocity density of the CACCNF fleet is concentrated in the range of 10-26 s. The velocity oscillation time of the CACCNF fleet is longer and the velocity fluctuation is higher compared with the CAVF fleet. Under this condition, the average fluctuation of both fleets is basically the same. From the results shown in Figure 17, it can be concluded that the acceleration changes in the CACCNF team are concentrated between 10 and 35 s, whereas the acceleration changes in the CAVF team are concentrated between 5 and 20 s; after 25 s, the acceleration changes become more concentrated and the change rate is smaller. As shown in Table 6, the mean jerk value of the CACCNF team is 62.8%, while that of the CAVF team is 53.9%, representing an 8.9% improvement in comfort level. The vehicle in the side lane cuts into the lane with a negative acceleration higher than the initial speed of the team train, as shown in Figure 17. From the results shown in Figure 17, it can be concluded that the acceleration changes in the CACCNF team are concentrated between 10 and 35 s, whereas the acceleration changes in the CAVF team are concentrated between 5 and 20 s; after 25 s, the acceleration changes become more concentrated and the change rate is smaller. As shown in Table 6, the mean jerk value of the CACCNF team is 62.8%, while that of the CAVF team is 53.9%, representing an 8.9% improvement in comfort level. The vehicle in the side lane cuts into the lane with a negative acceleration higher than the initial speed of the team train, as shown in Figure 17. From the results shown in Figure 17, it can be concluded that the acceleration changes in the CACCNF team are concentrated between 10 and 35 s, whereas the acceleration changes in the CAVF team are concentrated between 5 and 20 s; after 25 s, the acceleration changes become more concentrated and the change rate is smaller. As shown in Table 6, the mean jerk value of the CACCNF team is 62.8%, while that of the CAVF team is 53.9%, representing an 8.9% improvement in comfort level. Scenario 4: The vehicle in the side lane cuts into the lane with a negative acceleration higher than the initial speed of the team train, as shown in Figure 17. The statistical analysis based on the simulation results of Figure 17 is shown in Table  7. The regulation of the overall HWT of the fleet is reduced by 13% on average. The regulation rate of the fleet following distance is reduced by 6% on average. The results according to the figure show that the following distance of the CAVF fleet is lower, and the spacing is smaller. The statistical analysis based on the simulation results of Figure 17 is shown in Table 7. The regulation of the overall HWT of the fleet is reduced by 13% on average. The regulation rate of the fleet following distance is reduced by 6% on average. The results according to the figure show that the following distance of the CAVF fleet is lower, and the spacing is smaller. As shown in the figures below, with time as the horizontal axis and the position of the vehicles as the vertical axis, the heatmaps of the velocity and acceleration of each vehicle of the CAVF fleet and the CACCNF fleet are plotted separately. Figure 18 shows that the high-speed density of the CAVF fleet is concentrated within 0-16 s, and the high-speed density of the CACCNF fleet is concentrated within 0-20 s. The velocity fluctuation time of the CACCNF fleet is higher compared to the CAVF fleet, while the speed fluctuation time of the CAVF fleet is longer and the speed fluctuation rate is higher. From the results shown in Figure 19, it can be concluded that the acceleration changes in the CACCNF fleet are concentrated between 10 and 35 s, and the acceleration changes are more dispersed, while the acceleration changes in the CAVF fleet are concentrated between 5 and 20 s; after 25 s, the acceleration changes are more concentrated.   As shown in Table 8, the average CACCNF team jerk is 88% and the average CAVF team jerk is 74.5%, representing a 13.5% increase in comfort.

Real Vehicle Verification
The real vehicle test, conducted at the Tianjin University automated driving test site, features the BYD SUV parked in the left lane as the vehicle to be cut-in, and the automated driving fleet in the right lane, with the three vehicles separated by about 50 m behind the vehicle to be inserted, and the head vehicle starting as the signal. The vehicle to be inserted As shown in Table 8, the average CACCNF team jerk is 88% and the average CAVF team jerk is 74.5%, representing a 13.5% increase in comfort.

Real Vehicle Verification
The real vehicle test, conducted at the Tianjin University automated driving test site, features the BYD SUV parked in the left lane as the vehicle to be cut-in, and the automated driving fleet in the right lane, with the three vehicles separated by about 50 m behind the vehicle to be inserted, and the head vehicle starting as the signal. The vehicle to be inserted cuts into this lane after the speed of the convoy reaches the set speed of 10 km/h and drives at a uniform speed after cutting in. Figure 20 shows the road environment of the test site and the distribution of the vehicles.  During the actual car test, the test scenario and test environment of the two types of fleets are guaranteed to be the same as best as possible, such as the insertion distance and insertion speed of the vehicles in the side lane. The test results of the actual process are shown below (As shown in Figure 21-24, two CACCNF fleet tests and two CAVF fleet During the actual car test, the test scenario and test environment of the two types of fleets are guaranteed to be the same as best as possible, such as the insertion distance and insertion speed of the vehicles in the side lane. The test results of the actual process are shown below (As shown in Figures 21-24, two CACCNF fleet tests and two CAVF fleet tests at 10 km/h). Table 9 shows a comparison of the results. During the actual car test, the test scenario and test environment of the two types of fleets are guaranteed to be the same as best as possible, such as the insertion distance and insertion speed of the vehicles in the side lane. The test results of the actual process are shown below (As shown in Figure 21-24, two CACCNF fleet tests and two CAVF fleet tests at 10 km/h). Table 9 shows a comparison of the results.      As shown in the above table, the speed stabilization time domain of test 2 and test 4 is basically lower than the speed stabilization time domain of test 1 and test 3 by 1-2 s. The maximum speed difference between test 2 and test 4 is basically smaller than the maximum speed difference between test 1 and test 3 by 2.1-2.6 km/h. The maximum relative distance, i.e., the distance steadily followed by the queue vehicles, is also improved. The average speed of test 2 and test 4 is higher than the average speed of test 1 and test 3. The comfort level of the head vehicle is improved by 26%, and the comfort level of the rear vehicle is improved by 4.5%.

Discussion
The cooperative cruise control method described in this paper was implemented from the perspective of multi-vehicle flexible following. Therefore, the validation of the effectiveness of the flexible following algorithm was first performed on the joint simulation platform of Prescan 2019.3 and MATLAB2019/SIMULINK.
The multi-vehicle cooperative following algorithm was verified under four typical working conditions. The trajectory prediction algorithm could predict the future driving trajectory of the self-driving vehicle under four types of scenarios, as well as predict future driving maneuvers. Then, the prediction algorithm proposed in this paper and the cooperative algorithm of the CAVF fleet were comprehensively validated using the autonomous driving prototype platform. A multi-vehicle following test was conducted at the Tianjin University autonomous driving test site. The test results showed that the test prototype accurately identified the side lane vehicles and more accurately predicted the side lane change intention and lane change trajectory, while outputting the overlap degree of the side lane insertion into this lane. Lastly, the flexible following factor of the head vehicle and the flexible decay factor of the fleet train were output. In the cooperative control of multiple vehicles, the head vehicle speed stabilization time domain of the CAVF fleet was reduced by 1 s and 6 s, while the speed stabilization time domain of team train 2 was reduced by 5 s and 8 s compared with the fleet without the flexible following factor and flexible attenuation factor. The comfort level was improved by 26% and 25%, respectively, the fleet spacing was reduced by 5 m, and the average speed was more tangential to the bypass lane. The CAVF team is more comfortable than the CACCNF team, with a more compact workshop space. This paper develops a flexible following cooperative control strategy for queues based on vehicle prediction in the side lane. The strategy can control the longitudinal following distance as well as the speed of driverless networked vehicles in the queue by means of a flexible fading factor, flexibly adjusting the following distance, effectively suppressing traffic oscillations in the queue due to the insertion of manually driven vehicles in the side lane, and improving the overall stability of the hybrid fleet and reducing the fleet oscillation rate.
In future work, the approach developed in this paper could be applied not only to advanced driver assistance but also to autonomous driving systems. In the prediction layer, more advanced methods such as asymmetric game theory should be considered so that the prediction layer proposed in this paper can be extended to the remaining complex disclosures, such as intersections.

Conflicts of Interest:
The authors declare no conflict of interest.