Abstract
The construction of typical driving condition of vehicles in line with the actual road traffic conditions in China requires the selection of the same vehicle for two months to collect driving data and the obtention of 496000 driving condition data of light vehicles. The sample data are preprocessed by using multivariate statistical theory and MATLAB. After the elimination of abnormal data, the effective data are extracted before being divided into 3020 kinematic segments. Then, it takes a principal component analysis to reduce the dimension of the characteristic parameter matrix. Through Kmeans++ clustering algorithm, the six principal components obtained by principal component analysis are clustered into two categories. Then the typical kinematic segments are selected from various fragment libraries by using correlation coefficient method, so as to construct a typical driving condition of the vehicles in a certain city. With the application of PCAKmeans and PCAKmeans++ clustering algorithm, a driving condition curve with a duration of 1200s is constructed before its effectiveness and accuracy being compared and analyzed. The results show that the error rate of driving condition between sample data and driving condition constructed by PCAKmean++ clustering algorithm is less than 6 % and the error rate of average speed and acceleration standard deviation is less than 1 %. The correlation degree between working condition curve constructed by PCAKmeans ++ clustering algorithm and sample data is increased by 4.08 %. The proportion of deceleration time and idle time in vehicle driving state is obviously different, which indicates that PCAKmeans++ is a better way to solve the problem and the clustering algorithm can effectively construct the driving condition curve of light vehicles suitable for local cities.
1. Introduction
Vehicle driving condition is the main calibration for optimization of vehicle performance indicators, which can be used to evaluate vehicle exhaust emissions, fuel consumption and other indicators. Since this century, China has mainly used NEDC driving condition of Europe and world light test cycle (WLTC) working conditions as references to optimize vehicle performance and carry out calibration, so as to give authoritative results on vehicle economy and emission.
With the acceleration of urban construction, the urban road network is becoming increasingly complex and the traffic flow situation has changed greatly. The two driving condition standards of NEDC and WLTC can no longer be adapted to the traffic situation under the background of intelligent transportation. According to statistics, the error between the driving condition of vehicles based on NEDC and WLTC and the actual operating conditions of automobiles in China is as high as 35 %. Therefore, in a certain area, the construction of accurate vehicle driving condition can provide reference for the test of vehicle emission and fuel consumption, as well as the technical development and evaluation of new models.
At present, researchers all over the world have carried out indepth study on the driving condition of vehicles, fitted the driving condition of vehicles in line with the traffic characteristics of various regions. In the process of constructing the vehicle driving condition, most researchers combine principal component analysis (PCA) with Kmeans clustering algorithm in the construction of vehicle driving condition, which brought the following results. Cai [1] et al. applied Kmeans clustering method to construct the driving condition of Xi’an City, whereas an optimization of the clustering centre has not been carried out. Hence, its application has limitations to some extent. Qin [2] et al. used Kmeans clustering algorithm to construct driving conditions, which did not optimize the initial value and affected the accuracy. Shi [3] et al. improved the selection of initial value in Kmeans clustering method by using neural network method, but the clustering quality was reduced due to the lack of clear definition of parameter value K. Gao [4] et al. used the global Kmeans clustering algorithm to construct the driving conditions to optimize the cluster centre, which improved the robustness of traditional Kmeans clustering. However, the classification results of the Kmeans algorithm still depend on the selection of the initial point. Therefore, on the basis of Kmeans, this paper uses Kmeans ++ clustering algorithm to construct vehicle driving condition, analyses the accuracy and effectiveness of vehicle driving condition, and finally constructs the driving condition of light vehicle in a city of China.
Through the investigation of vehicles and road conditions in a city, it is found that the car ownership is more than 2.7 million, while the construction of passenger car driving condition is less. Therefore, this paper uses the realtime vehicle driving condition data collected in a city as the sample, using MATLAB to preprocess the original data, dividing the kinematic segments and Kmeans++ clustering, and in the end obtaining the best among all kinds of cluster centres. At the same time, the abnormal data are processed to extract the kinematic segments, which accord with the actual motion characteristics of the vehicle. Then, the dimension of the collected data is reduced based on kinematic fragment analysis and principal component analysis, 3020 segments are calculated with Java language, the corresponding kinematic segment eigenvalues are obtained, the driving condition of a city road is constructed. Finally, the accuracy of the driving condition is verified by comparing with kmeans algorithm and the total sample data collected.
2. Model approach
2.1. Kinematic fragment
The process of frequent starting, accelerating and decelerating is actually a kinematic process. This process is affected by traffic conditions or traffic flow. Therefore, there will be continuous starting and stopping operations during the process from starting to stopping, so it can be regarded as composed of multiple “stop drive stop” segments. In order to simulate the frequent starting, accelerating, decelerating and other driving states of the vehicle in the actual driving process, the driving process from one idle state to the next is defined as a kinematic segment, which is usually composed of an idle part and a moving part [5], as shown in Fig. 1.
Fig. 1Schematic diagram of kinematics
In order to fully express the characteristics of each kinematic segment, according to research [6, 7], the division criteria of four operating conditions are defined in this article.
1) Idle speed: speed $V=$ 0, the engine is idling.
2) Acceleration: Driving state with speed $V\ne $ 0 and $a\ge \mathrm{}$0.1 m/s^{2}.
3) Deceleration: Driving state with speed $V\ne $ 0 and $a\le $ –0.1 m/s^{2}.
4) Constant speed: Driving state with speed $V\ne $ 0 and $\lefta\right\le $ 0.1 m/s^{2}.
Among them, v represents the speed, a represents the acceleration.
2.2. Principal component analysis
Principal component analysis (PCA) is a statistical method to reduce the dimension of highdimensional data. According to the projection principle, the highdimensional data is projected on the lowdimensional space with the minimum possibility, which can greatly simplify the data structure. At the same time, it is also a method to combine multiple related variables into a few unrelated variables based on the principle of minimizing data loss. Specific steps are shown in Fig. 2.
Fig. 2Calculation steps of principal component analysis
The calculation process is as follows:
Step 1: standardization of index parameters
Since the dimension of each parameter is different, the value of variables is dispersed greatly and the changing range of data affects the clustering effect. Therefore, it is necessary to standardize and dimensionless process each parameter. The Zscore method, which is a nondimensional method, is used to normalize the parameters of the raw matrix before obtaining the Eq. (1):
where $\overline{{X}_{j}}$ is the mean value of the $j$th object and ${S}_{ij}$ is the standard deviation of the $j$th object.
Step 2: Determination of the correlation coefficient matrix $R$ between parameters:
Step 3: Determination of the characteristic roots of $R$ under the form of ${\lambda}_{g}=(\mathrm{1,2},\dots ,m)$, which describes the importance of each component:
Step 4: Determination of the eigenvectors of the $R$ matrix, where $L$ is an mdimensional vector:
Step 5: Determination of the contribution rate of the $R$array:
Step 6: Determination of the number $K$ of principal components. Sort the components according to the contribution rate. If the original number is large, the first $K$ are taken to represent all the quantities:
According to determining the number of principal components, $m$ related indicators can be reduced to $K$ unrelated principal components:
2.3. Kmeans++ clustering method
2.3.1. Kmeans++ algorithm theory
Aiming at big data clustering, Kmeans++ algorithm is the optimization of random initialization centroid method of Kmeans. Comparing with other Kmeans clustering algorithms, the advantage of Kmeans++ clustering algorithm is that it takes the farthest possible distance as the data selection principle as the initial clustering centre. Besides, the selected initial cluster centre can be optimised. Therefore, the algorithm solves the problem of initial value selection, improves the stability of the algorithm and the robustness of clustering results that rely too much on the initial cluster centre.
The method of Kmeans++ clustering algorithm to select the initial center point is: first, considering the distribution of all samples in the data set, the distance between the selected initial center points should be as far as possible, so as to avoid the problem of large clusters being divided or several smaller clusters being wrongly merged.
Assuming that$D=\left\{{X}_{i}\left{X}_{i}\right.\in {R}^{m},i=\mathrm{1,2},\dots ,n\right\}$ represents the data set of the sample, $m$ represents the data dimension and n is the set size. Assuming that $C=\{{C}_{k}\left{C}_{k}\right.\in {R}^{m},i=\mathrm{1,2},\dots ,k\}$ represents the class to which the sample set belongs, then $k$ represents the number of category belonged, ${C}^{0}$ represents the initial clustering centre, the formula for calculating the distance between samples is as following:
Cluster centre:
The basic theory of Kmeans++ clustering [8] is that based on the set $K$ value and algorithm, samples are continuously selected from the class whose cluster centre is randomly chosen. When $E$ is the smallest (setting function is optimal), the process terminates. The result of the algorithm is to divide a large number of samples into different classes, showing different population aggregation distributions. The characteristic is that the distance is small between the internal sample elements of the class, while the distance is large between the classes:
The workflow of Kmeans++ algorithm is as follows:
1) A sample in dataset $D$ is randomly selected as the first initial cluster centre;
2) The shortest distance between each sample and the current cluster centre is calculated, which is represented by $D\left(x\right)$;
3) The probability that each sample is selected as the next cluster centre is calculated;
4) According to the principle that probability is proportional to $D\left(x\right)$, a new sample is selected as the cluster centre;
5) When the number of initial cluster centers $K$ is selected, the iteration process is terminated;
6) Re implement the traditional Kmeans algorithm.
2.4. Correlation analysis
Correlation analysis is an analysis of two or more related variables to describe the closeness between the variables. Correlation coefficient is a statistical analysis index for analysing the correlation degree of straight lines, and the value range is between –1 and 1. –1 indicates that the variables are completely uncorrelated, 0 indicates that they are not related and 1 indicates that they are completely related. The closer the data tends to 0, the weaker the correlation is. The specific calculation formula is shown in Eq. (13):
Among them, $Cov(X,Y)$ is the covariance of variables $X$ and $Y$, $D\left(X\right)$ and $D\left(Y\right)$ are variances of variables $X$ and $Y$, respectively. When ${\rho}_{xy}\ge \text{0.9}$, high correlation is generally believed between two groups of variables.
3. Data acquisition and processing
3.1. Data collection
The commonly used vehicle data acquisition methods include vehicle tracking method, average vehicle flow statistics method, autonomous driving method, etc. Among them, the autonomous driving method does not need to plan the specific test driving route in advance. Instead, the owner can drive autonomously according to the normal driving habits. In order to make the collected data more effective and fully express the traffic characteristics of a certain sound City, the selected vehicle should consider the location of the owner and the road condition of the owner every day. Therefore, the daily driving route of the selected vehicle should include places with frequent traffic interaction, such as residential areas, schools, stations, parks, etc.
In this paper, the autonomous driving method is used to collect the driving track data of different driving areas and driving periods and the vehicle terminal (GPS) is used to collect the driving data. The data acquisition frequency is set at 1Hz and the realtime vehicle driving data is recorded for about two months from November 1, 2017 to December 24, 2017, including vehicle position (longitude and latitude), speed, time, vehicle fuel consumption, etc. Finally, 496014 driving data were obtained.
3.2. Data preprocessing
In the construction process of vehicle driving condition, the quality of collected data will inevitably decline due to various reasons such as unstable transmission signal, electromagnetic interference, decoding error and so on. In order to improve the data quality and ensure the credibility of the research results, this paper uses MATLAB software to preprocess the original data.
In order to convert a given date string (date) into a date number (time), the datenum function in MATLAB software is used to process the time item according to the following format: Time=datenum(date,'yyyy/mm/dd HH:SS.000').
Since there are some errors in the original data collected directly, this paper analyzes and summarizes the error data according to the structure types, processes them in batches. For the discrete data, linear interpolation is carried out for the abnormal data before the abnormal data being eliminated.
3.2.1. Intermittent data processing
Discontinuities will occur in the data due to interference from weather changes, obstructions, poor signal reception, etc. The data in this case is defined as discontinuous data. The solution for this problem is demonstrated as follows. Clips are deleted that time intervals are greater than 10 s; For data with a time interval less than 10 s, linear interpolation is used to supplement the points. Before and after the discontinuous data is eliminated, the variance of the velocity is shown in Fig. 3 (a) and (b), respectively.
Fig. 3Intermittent data processing
a) Before processing
b) After processing
3.2.2. Acceleration and deceleration abnormal data processing
Due to the limitation of the road environment, the accuracy of GPS collected data will be affected, resulting in some abnormal signals. The abnormal data of car acceleration and deceleration is one of them. The acceleration time from 0 to 100 km/h is more than 7 seconds for an ordinary car, and the maximum deceleration is between 7.58 m/s^{2} when emergency brakes. According to related bibliography [9], it is defined as abnormal data for the tested vehicle when the acceleration exceeds 3.97 m/s^{2} or the deceleration exceeds 8 m/s^{2}, then such data will be deleted. Acceleration and deceleration abnormal data before and after removal are shown in Fig. 4(a) and (b), respectively.
Fig. 4Acceleration and deceleration abnormal data processing
a) Before elimination
b) After removal
3.2.3. Idle abnormal data processing
Idle is a continuous process in which the car stops moving but the engine is running at the lowest speed. When the vehicle stopped for a long time, abnormal idling data will appear in the collected valid data due to equipment accuracy. One is the case of stopping and not turning off, that is, the clip with zero speed is less than 20 s; the other is that the car is stopped and the device is to continue running. It is also in a state of zero speed. For this data, a rule aiming for idling abnormal data is formulated, which is filtered and eliminated according to the rules. Fig. 5(a) and 5(b) show the variation of speed before and after the idling abnormal data being eliminated.
Fig. 5Idle abnormal data processing
a) Before processing
b) After processing
3.2.4. Glitch data processing
During the driving process, the idling data collected has abnormal with nonzero speed due to factors such as poor signal. This data is defined as glitch data. This situation does not exist during the actual driving of the car. This paper uses such data, whose continuous driving time is no more than 3 seconds, at the same time its speed is less than 10 km/h, as glitch data. Reference [10] modified this data to zero to identify it as the idling point. Fig. 6(a) and (b) show the glitch data before and after removal.
Fig. 6Glitch data processing
a) Before processing
b) After processing
3.2.5. Long idle time abnormal data
The abnormal idling time refers to the data in which the segment of the holding speed is more than 180 seconds when the vehicle is stationary and the engine continues to work. If the idle is not in this range, one possibility is that the device is abnormal, and the car is not powered off after the car stops, and the device continues to collect information; another is that the car has been left in the same location for a long time without driving. These inconsistent data should be eliminated to reduce the error rate of the construction of the vehicle operating conditions. Fig. 7(a) and (b) before and after removal of the abnormal idling time abnormal data, respectively.
Valid 457550 driving data was obtained after discontinuous data being interpolated, acceleration and deceleration abnormal data being eliminated and idling being deleted, glitch being modified, and long idling time being proceeded.
Fig. 7Abnormal data processing at idle time
a) Before processing
b) After processing
4. Construction of driving condition model and result analysis
On the basis of data preprocessing, the driving condition data can be divided into kinematic segments and different number of characteristic parameters can be selected as principal component analysis.
4.1. Kinematic segmentation
According to the definition of kinematic fragments, the first data is divided into multiple motion segments. Then fragments are select based on valid segment filtering rules. Python programming is used to divide the microstroke due to the large amount of experimental data. Since the frequency is 1 Hz, when the fragments are performed, each line of data is reduced to one line, which includes the vehicle speed only [11, 12]. A total of 3020 kinematic segments are obtained by dividing the preprocessed data into kinematic segments. The vehicle operating conditions will be built based on this. The schematic diagram of kinematics is shown in Fig. 8.
Fig. 8Schematic diagram of kinematics
4.2. Construction of driving conditions model
This paper defines 16 feature parameters to describe the operating conditions of all kinematic segments, and 22 to describe the overall distribution characteristics of kinematic segments according to the importance of feature parameters to vehicle driving information. See Table 1 and Table 2 for information on each characteristic parameter. By combining the calculated feature parameter values of all the motion segments, a sample of each segment can be obtained before establishing a matrix of sample numbers (rows) × feature parameters (columns), which is used for the construction of driving conditions.
Table 1Kinematics segment running characteristic parameter information table
Number  Characteristic parameters  Unit  Number  Characteristic parameters  Unit  
1  Acceleration time ${T}_{a}$  %  9  Maximum speed ${V}_{max}$  km/h  
2  Deceleration time ${T}_{d}$  %  10  Minimum speed ${V}_{min}$  km/h  
3  Uniform time ${T}_{c}$  %  11  Acceleration ${A}_{max}$  m/s^{2}  
4  Idle time ${T}_{i}$  %  12  Deceleration ${D}_{max}$  m/s^{2}  
5  Operation hours ${T}_{r}$  S  13  Standard deviation of speed ${V}_{s}$  km/h  
6  Cumulative driving distance $S$  km  14  Standard deviation of acceleration ${A}_{s}$  m/s^{2}  
7  Average speed ${V}_{m}$  km/h  15  Average acceleration ${A}_{m}$  m/s^{2}  
8  Average running speed ${V}_{mr}$  km/h  16  Average deceleration ${D}_{m}$  m/s^{2} 
Table 2Feature parameter information table of the overall distribution of kinematics
Number  Characteristic parameters  unit  number  Characteristic Parameters  Unit 
1  Acceleration time ratio ${P}_{a}$  %  12  Speed ratio of 7080 km/h ${P}_{7080}$  % 
2  Slowdown ratio ${P}_{d}$  %  13  Ratio of speed greater than 80 km/h ${P}_{80}$  % 
4  Idle time ratio ${P}_{i}$  %  15  –2 m/s^{2} –1.5 m/s^{2} acceleration ratio ${A}_{2~1.5}$  % 
5  Speed ratio between 0~ 0 km/h ${P}_{010}$  %  16  –1.5 m/s^{2} –1 m/s^{2} acceleration ratio ${A}_{1.5~1}$  % 
6  Speed ratio of 1020 km/h ${P}_{1020}$  %  17  – l m/s^{2} –0.5 m/s^{2} acceleration ratio ${A}_{1~0.5}$  % 
7  Speed ratio of 2030 km/h ${P}_{2030}$  %  18  – 0.5 m/s^{2}0.5 m/s^{2} acceleration ratio ${A}_{0.5~0.5}$  % 
8  Speed ratio of 3040 km/h ${P}_{3040}$  %  19  0.5 m/s^{2}1 m/s^{2} acceleration ratio ${A}_{0.5~1}$  % 
9  Speed ratio of 4050 km/h ${P}_{4050}$  %  20  1m/s^{2}1.5 m/s^{2} acceleration ratio ${A}_{1~1.5}$  % 
10  Speed ratio of 5060 km/h ${P}_{5060}$  %  21  1.5 m/s^{2}2 m/s^{2} acceleration ratio ${A}_{1.5~2}$  % 
11  Speed ratio of 6070 km/h ${P}_{6070}$  %  22  Ratio of acceleration section above 2 m/s^{2}${A}_{2}$  % 
A program based on the Java language is written to serve for the calculation 3020 fragments and obtain the corresponding running feature values of the kinematic fragments. The characteristic values of all fragments are shown in Table 3.
4.2.1. PCA analysis results
According to the steps of principal component analysis, the eigenvalue matrix is analysed, and six principal components are obtained, which are labelled as ${M}_{1}$, ${M}_{2}$, ${M}_{3}$, ${M}_{4}$, ${M}_{5}$, ${M}_{6}$, and the contribution rate and cumulative contribution rate are calculated, the results are shown in Table 4.
Table 3Kinematic segment feature parameter values
Fragment number  ${T}_{d}$  ${T}_{c}$  ${V}_{m}$  ...  ${V}_{s}$ 
1  24  76  65.31  ...  17 
2  39  45  51.11  ...  16.78 
3  40  29  53.29  ...  23.44 
4  12  33  34.22  ...  0.64 
5  26  24  65.8  ...  0.53 
6  17  31  115.31  ...  0.66 
7  25  19  87.94  ...  19.84 
8  33  17  57.4  ...  17.78 
9  41  60  107.31  ...  23.61 
...  ...  ...  ...  ...  … 
3019  9  22  69.66  ...  24.21 
3020  36  34  78.17  ...  0.78 
Table 4Analysis results of the main components
Main ingredient  Eigenvalues  Contribution rate (%)  Cumulative contribution rate (%) 
${M}_{1}$  5.6174  45.8563  95.6584 
${M}_{2}$  3.1933  19.2535  92.2542 
${M}_{3}$  1.7534  16.2543  85.5463 
${M}_{4}$  0.8524  7.2514  63.2015 
${M}_{5}$  0.2563  3.2541  45.8563 
${M}_{6}$  0.0452  1.2542  39.7815 
By solving the characteristic equation $\left{\lambda}_{g}{I}_{m}R\right=0$, m nonnegative eigenvalues of the correlation coefficient matrix $R$ can be obtained and arranged in order of magnitude. The larger the eigenvalue of the correlation coefficient matrix is, the better the initial variable can be represented by the representative principal component to a greater extent. The contribution rate represents the degree of reflection on it. Generally, when the cumulative contribution rate reaches 85 %, the information of the initial variable can be reflected. It is concluded from Table 4 that the contribution rate of ${M}_{1}$, ${M}_{2}$, and ${M}_{3}$ reaches 85 %, which can represent the information of the initial variables, and their eigenvalues are greater than 1, so the initial characteristic parameters can be characterized by ${M}_{1}$, ${M}_{2}$, and ${M}_{3}$.
Table 5Calculation results of the principal component factor load matrix
Characteristic parameters  ${M}_{1}$  ${M}_{2}$  ${M}_{3}$ 
$T$  0.559  0.299  0.101 
$S$  –0.022  0.846  –0.136 
${T}_{i}$  0.697  –0.194  0.376 
${T}_{a}$  0.043  0.714  –0.083 
${T}_{d}$  0.091  0.615  –0.384 
${T}_{e}$  0.116  0.263  –0.217 
${V}_{max}$  –0.225  0.947  0.083 
${V}_{m}$  –0.442  0.866  –0.101 
${V}_{mr}$  –0.135  0.918  0.018 
${V}_{s}$  –0.316  0.941  0.290 
${A}_{max}$  0.643  0.535  –0.595 
${A}_{m}$  0.478  0.464  –0.430 
${D}_{max}$  0.064  0.689  0.533 
${D}_{m}$  0.064  0.680  0.532 
${A}_{s}$  0.386  0.544  –0.038 
According to the contribution table of the principal components, we choose to end the accumulation when the cumulative contribution rate exceeds 80 %. The factor load matrix is calculated by the characteristics of the first three components whose cumulative contribution rate is greater than 85 %. The calculation results are shown in Table 5.
As can be seen from Table 5, ${M}_{1}$ mainly reflects the idle time; ${M}_{2}$ mainly reflects the running distance, acceleration time, deceleration time, maximum speed, average speed, etc.; ${M}_{3}$ mainly reflects the maximum acceleration and the average deceleration zone. These three factors can fully reflect the characteristics of the kinematic segment.
4.2.2. Kmeans clustering results and analysis
According to the algorithm steps, ${M}_{1}$, ${M}_{2}$, and ${M}_{3}$ are divided into categories, and the data is analyzed through python. The original 3020 fragments are divided into two categories to ulterior analyses the driving characteristics of each kind of motion segment. According to the formula, the comprehensive eigenvalues of all kinds of motion segments and whole data are obtained, which have been shown in Table 6.
Table 6Comprehensive eigenvalues of different classes
Characteristic parameters  The first type  Second category 
$T\left(s\right)$  57339  73472 
$S\left(m\right)$  385658.046  481129.346 
${T}_{i}\left(s\right)$  15986  17870 
${T}_{a}\left(s\right)$  15124  20880 
${T}_{d}\left(s\right)$  16127  12177 
${T}_{e}\left(s\right)$  10102  22545 
${V}_{max}$_{}(km/h)  61  76 
${V}_{m}$ (km/h)  18.576  31.435 
${A}_{max}$ (km/h)  3.56  2.56 
${A}_{m}$ (km/h)  0.525  0.748 
${A}_{max}$ (km/h)  6.611  4.833 
${D}_{m}$ (km/h)  0.749  0.525 
According to the transition probability matrix of category 1 and category 2, the driving condition of 1200 second typical road vehicle is synthesized by principal component analysis and Kmeans algorithm, shown in Fig. 9.
Fig. 9Representative driving condition curve of Kmeans clustering
4.2.3. Kmeans++ clustering results and analysis
According to the transition probability matrices of category 1 and category 2, the principal component analysisKmeans++ method is used to synthesize a typical road driving condition map for 1200 seconds, as shown in Fig. 10. The traffic congestion in category 1 takes 824 seconds, while Category 2 is 381 seconds under relatively smooth conditions.
Fig. 10Representative driving conditions of Kmeans++ clustering
4.3. Comparative analysis of effectiveness
4.3.1. Comparative analysis of operating characteristics
The feature parameters such as ${V}_{m}$ (km/h), ${V}_{mr}$ (km/h), ${A}_{m}$ (m/s^{2}), ${D}_{m}$ (m/s^{2}), ${P}_{i}$ (%), ${P}_{a}$ (%), ${P}_{d}$ (%), ${P}_{c}$ (%), ${V}_{s}$ (km/h), ${A}_{s}$ (m/s^{2}) are selected to evaluate the effectiveness of automobile motion characteristics the same as the operating condition curve constructed by PCAKmeans. Then it is compared and analysed between the representative construction conditions of the final construction and the sample parameters, in the end the error rates are calculated respectively (see Table 7).
Table 7Error probability analysis of Kmeans++ construction conditions relative to sample data
Characteristic parameters  Test data  PCAKmeans combination model  PCAKmeans++ combination model  
PCAKmeans  Error  PCAKmeans++  Error  
${V}_{m}$  16.481  16.995  3.12%  16.515  0.21% 
${V}_{mr}$  22.497  24.585  9.28 %  23.620  4.99 % 
${V}_{s}$  11.462  12.237  6.76 %  11.884  3.68 % 
${A}_{m}$  0.545  0.590  8.34 %  0.572  5.02 % 
${D}_{m}$  0.543  0.593  9.12 %  0.580  5.75 % 
${A}_{s}$  0.749  0.773  3.25 %  0.755  0.74 % 
${P}_{i}$  0.224  0.236  5.43 %  0.231  3.26 % 
${P}_{a}$  0.285  0.303  6.21 %  0.296  3.70 % 
${P}_{d}$  0.292  0.310  6.17 %  0.302  3.45 % 
${P}_{e}$  0.199  0.210  5.59 %  0.208  4.37 % 
It can be seen from Table 7 that the error rate of representative working conditions and collected data constructed by PCAKmeans combination model is as high as 9 % and the minimum error rate is more than 3%; while the one constructed by PCAKmeans++ combination model is less than 6 %, and the error rate of some parameters (such as average speed and acceleration standard deviation) is less than 1 %. This fully shows that the statistical distribution characteristic parameters and driving condition proportion distribution law constructed by PCAKmeans++ combination model are consistent with the collected data, which are obviously higher than the working condition error rate constructed by PCAKmeans algorithm. It can be seen that the application of PCAKmeans++ combined model to build driving condition can better reflect the distribution of collected data in each speed segment, significantly improve the accuracy of driving condition construction and can more accurately reflect the driving state of a city vehicle, with stronger reliability.
Taking speed frequency as an example, the paper compares and analyzes the distribution of vehicle driving state (Fig. 11). Kmeans and Kmeans++ methods have less error with sample data in reflecting vehicle driving state. However, the error rate between the representative working conditions constructed by Kmeans++ algorithm and the overall collected data is smaller, especially the proportion of deceleration time and idle time. Therefore, it can be considered that the driving condition constructed by Kmeans++ method is more suitable for local needs and can reflect the actual driving characteristics of local vehicles. It has certain theoretical and practical value for performance calibration of fuel consumption and emission.
Fig. 11Proportion distribution of vehicle driving conditions
4.3.2. Correlation analysis
The correlation coefficient between the representative working condition curve and the sample data curve constructed by PCAKmeans and PCAKmeans++ combination model was calculated respectively, shown in Table 8.
Table 8Comparison of correlation coefficients
Model  Correlation coefficient 
PCAKmeans combination model  0.9276 
PCAKmeans++ combination model  0.9684 
The results show that the correlation coefficient between the representative working condition curve constructed by PCAKmeans combination model and sample data is 0.9276, while the one constructed by PCAKmeans++ combination model and sample data is 0.9684. Compared with PCAKmeans++ combination model, the correlation degree between driving condition curve and sample data is higher, which is 4.08 % higher than that of PCAKmeans. Therefore, PCAKmeans++ combination model can better reflect the actual driving level and road traffic characteristics of vehicles on the road.
5. Conclusions
In this paper, taking a city in China as an example, the realtime track of the same vehicle in two months is recorded, and a total of 496014 driving data are obtained.
1) In order to make the data accurate and reasonable, MATLAB software is used to preprocess the original data. Through interpolation of discontinuous data, screening and elimination of abnormal data, burr data processing and other operations, 457550 pieces of effective vehicle driving data are finally obtained. Python programming is used as well to divide the big data into micro stroke, the preprocessed data is divided into kinematic segments. A total of 3020 kinematic segments are obtained. After several selections of different number of feature parameters as principal component analysis, 16 feature parameters are selected to describe the operation conditions of all kinematic segments and 22 feature parameters are used to describe the overall distribution characteristics of kinematic segments. Besides, principal component analysis and Kmeans++ clustering algorithm are used to reduce dimension and classify the characteristic parameter matrix, so as to construct a vehicle driving condition of 1200 s in accordance with the traffic characteristics of a city.
2) Compared with the traditional Kmeans clustering algorithm, the vehicle driving condition model constructed by Kmeans++ algorithm can effectively solve the problem of unstable initial centre point and improve the accuracy of driving condition construction. It is found that the average relative error of driving condition generated by Kmeans++ algorithm is less than 6 %, the standard deviation error rate of average speed and acceleration is less than 1 %. The ratio of deceleration time to idle time is smaller in the running state of the vehicle. At the same time, the correlation between the vehicle driving condition curve and the sample data is higher than 4.08 %, which verifies the accuracy and effectiveness of the Kmeans++ clustering algorithm. Therefore, compared with NEDC, the driving condition constructed by Kmeans++ algorithm is more suitable for local needs and can reflect the actual driving characteristics of local vehicles. It has a certain theoretical and practical value for the calibration of vehicle fuel consumption, emissions and other performance indicators
3) Through the successful construction of the vehicle driving condition suitable for a certain city, it can provide certain reference for the construction of the actual driving characteristics of vehicles in other cities.
References

Cai Yan, Li Yangyang, Li Chunming, Tan Xiaowei, Liu Dongmin Research on synthetic technology of vehicle driving conditions in Xi’an based on KMeans clustering algorithm. Automobile Technology, Vol. 8, 2015, p. 3336.

Qin Datong, Zhan Sen, Qi Zhenggang, Chen Shujiang Construction method of driving conditions based on Kmeans clustering algorithm. Journal of Jilin University (Engineering and Technology Edition), Vol. 46, Issue 2, 2016, p. 383389.

Shi Qin, Qiu Duoyang, Zhou Jieyu Driving condition construction and accuracy analysis based on combined clustering method. Automotive Engineering, Vol. 34, Issue 2, 2012, p. 164169+158.

Gao Jianping, Wu Jianguo Construction of vehicle driving conditions based on global Kmeans clustering algorithm. Journal of Henan University of Science and Technology (Natural Science), Vol. 38, Issue 1, 2019, p. 112118.

Shi Qin, Zheng Yubo, Jiang Ping Research on driving conditions of urban roads based on kinematics segments. Automotive Engineering, Vol. 33, Issue 3, 2011, p. 256261.

Knez M., Muneer T., Jereb B., et al. The estimation of a driving cycle for Celjeanda comparison to other European cities. Sustainable Cities and Society, Vol. 11, 2014, p. 5660.

Lin J., Niemeier D. A. Exploratory analysis comparing a stochastic driving cycle to California’s regulatory cycle. Atmospheric Environment, Vol. 36, Issue 38, 2002, p. 57595770.

Liu Ye, Wu Sheng, Zhou Haihe, Wu Xingyi, Han Linyi Research on Kmeans clustering algorithm optimization method. Information Technology, Vol. 43, Issue 1, 2019, p. 6670.

Li Yang Research on Car Driving Conditions Based on Clustering Algorithm. Beijing Institute of Technology, 2016.

Tian Yu Construction and Analysis of Road Driving Conditions for Light Vehicles in Taiyuan City. Taiyuan University of Technology, 2018.

Zhang Zhiyong, et al. Master of MATALB. Beijing University of Aeronautics and Astronautics Press, Beijing, 2011.

Sun Qiuyue, Wang Lizhen, Wu Fengping Introduction and application of data compression method. Journal of Yunnan University (Natural Science Edition), Vol. 29, Issue 1, 2007, p. 115118.

Ma Zhixiong, Li Mengliang, Zhang Fuxing, et al. Application of principal component analysis to the development of vehicle actual driving conditions. Journal of Wuhan University of Technology, Vol. 26, Issue 4, 2004, p. 3235.

Ai Guohe, Qiao Weigao, Li Mengliang, et al. Analysis of the kinematic parameters of vehicle driving. Highway Communications Science and Technology, Vol. 23, Issue 2, 2006, p. 154157.

Shi Shuming, Wei Shuying, Hai Linkui, Liu Liu Improvement soft the design method for transient driving cycle for passenger car. IEEE Vehicle Power and Propulsion Conference, 2009, p. 15811586.