Open AccessArticle

Development Method for the Driving Cycle of Electric Vehicles

Zhecheng Jing

Tianxiao Wang

Shupei Zhang

and

Guolin Wang

School of Automotive and Traffic Engineering, Jiangsu University, Zhenjiang 212013, China

Author to whom correspondence should be addressed.

Energies 2022, 15(22), 8715; https://doi.org/10.3390/en15228715

Submission received: 25 October 2022 / Revised: 17 November 2022 / Accepted: 18 November 2022 / Published: 20 November 2022

Download

Browse Figures

Versions Notes

Abstract

With the development of electric vehicles, more attention has been paid to the role of the driving cycle in vehicle performance testing. At present, the K-means algorithm is often used in the development of driving cycles. However, it is sensitive to the outlier points and also difficult to determine the K value. To solve this problem, the hierarchical cluster method is applied in this study. First, the real-world driving data are collected and denoised through wavelet domain denoising. Then, the data are divided into micro-trips and the characteristic parameters are extracted. The hierarchical cluster method is adopted to classify the micro-trips into different categories. An appropriate number of micro-trips are selected from each group in proportion to each category to assemble the driving cycle. Finally, both the economic simulation and the statistical analysis prove the accuracy of the generated driving cycle and the feasibility of the development method proposed in this paper.

Keywords:

electric vehicle; driving cycle; hierarchical cluster method; principal component analysis

1. Introduction

The continuous and rapid increase in car ownership is an important cause of air pollution and the energy crisis. The emission pollution of traditional vehicles harms the environment and human health. Vehicle emissions have become the main source of air pollution in most countries worldwide. Twenty-five percent of the global CO₂ emissions come from automobiles, and the total emissions of greenhouse gases are proportional to the number of automobiles. Additionally, traditional cars are strongly dependent on oil, and volatile crude oil prices and the energy security strategies of countries have constrained the development of the auto industry. The development of electric vehicles cannot only effectively solve the emissions problem but also reduce the use of fuel, which is of great significance. The transformation of the automotive industry has already begun, and it is fast-tracked for further development, especially in China. The sales volume of new energy vehicles in China was less than 20,000 in 2013, but in 2021, 3.521 million units were sold, and 3.545 million units were produced, ranking first in the world in terms of production and sales for several consecutive years. This effectively reduces China’s pollution problems and energy problems.

The electric vehicle industry is transitioning from high-speed development to high-quality development, with higher performance testing requirements. Specifying a uniform condition for tests is necessary, which is why driving cycles need to be developed. Moreover, proper driving cycles are good for research on reducing energy consumption [1] A driving cycle is a series of data points representing the speed of a vehicle versus time [2], and it is the basis of testing. Europe uses the New European Driving Cycle (NEDC), while JC15 is used in Japan, and FTP75 is adopted in the USA. These are the world’s three most widely used driving cycles [3,4]. However, traffic varies significantly from country to country, and NEDC, JC15, and FTP75 cannot represent the local traffic conditions, which leads to inaccurate test results. Therefore, building a driving cycle suitable for the local area is of practical significance.

Constructing driving cycles has always been the frontier and focus of research in the automotive field. The micro-trips method is the most common approach [5]. The driving data are cut into micro-trips and then divided into several clusters. The K-means clustering method is the most widely used clustering method for generating a driving cycle [6,7,8]. Zhang proposed combining the distance optimization and density methods to construct a data set density metric method to construct a driving cycle as an improved K-means algorithm method [9]. A global K-means clustering algorithm was used to construct the Zhengzhou passenger car driving cycles [10]. Yuan proposed an improved K-means text clustering algorithm based on density peaks that avoided the pitfalls of randomly picking initial clusters [11]. However, K-means clustering has significant defects that will negatively affect the construction of a driving cycle. For instance, K-means is sensitive to the outlier points, which will introduce uncertainty in the clustering results.

Additionally, it is difficult to determine the K value when the amount of data is large and the K value is affected by the initial clustering center. To solve this problem, this paper proposes the methodology of generating a driving cycle based on a hierarchical cluster method. A hierarchical cluster is not sensitive to outlier points, and the clusters can be determined according to actual needs, which means the clustering result is controllable with a faster and more effective process.

In this study, the data collection and preprocessing are introduced first. Then, the real-world data are divided into micro-trips, and the characteristic parameters are extracted. The principal component analysis is applied for dimensionality reduction, and then the micro-trips are divided into three categories using the hierarchical cluster method: low-speed, medium-speed, and high-speed. Next, the driving cycle construction is completed. Finally, the validity of the constructed working conditions and the feasibility of the proposed theory is verified from the perspective of statistical analysis and economic simulation.

2. Data Collection and Preprocessing

Driving data is the basis of driving cycle development. To obtain the driving data reflecting the actual driving characteristics, this paper uses the autonomous driving method to collect the driving data of an electric vehicle. The driver drives the vehicle according to their own will without a designated driving route, which can reflect the real traffic situation to the greatest extent. The test system comprises an onboard data acquisition terminal (sampling frequency is 1 Hz) (Figure 1) and a data management platform. The onboard data acquisition terminal codes the acquired information according to the unified data protocol and sends it to the online data management platform in real-time through the GPRS network. The data source of the onboard data acquisition terminal includes two parts, including the GPS signal and the onboard diagnostic system (OBD) signal.

Due to the complex electromagnetic environment in electric vehicles, signal noise is inevitable, so it is necessary to filter the data. In this study, wavelet domain denoising is adopted, which can remove noise and preserve the information for the original data to the greatest extent. The db5 wavelet is selected to denoise and analyze the speed data; the level is 5 and the SURE threshold is used to eliminate the signal noise [12]. To make the comparison more intuitive, the denoised case is shown in Figure 2.

The micro-trip refers to the data period from the start of one idle to the start of the next idle, and it is the basic unit for constructing the driving cycle. In this study, a total of 1162 micro-trips are obtained after the division of the driving data.

To describe the micro-trips, 15 characteristic parameters are counted [13], including the average speed, maximum speed, maximum acceleration, average acceleration, maximum deceleration, average deceleration, acceleration time, deceleration time, idle time, uniform time, total time, acceleration time ratio, deceleration time ratio, idle time ratio, and uniform time ratio. The characteristic parameters and the symbols are shown in Table 1.

The 15 characteristic parameters are defined as follows:

Average speed:

V_{a v g} = \sum_{i = 1}^{T} v_{i} / T

(1)

v_{i}

is the vehicle speed at the

i_{t h}

second.

Maximum speed: the maximum speed in the micro-trip.

Maximum acceleration: the maximum acceleration in the micro-trip.

Average acceleration:

a_{a v g_a} = \sum_{i = 1}^{T} a_{i} / T

(2)

a_{i}

corresponds to the acceleration at the

i_{t h}

second and

n

is the total time of the micro-trip.

Maximum acceleration: the maximum acceleration in the micro-trip.

Average deceleration:

a_{a v g_d} = \sum_{i = 1}^{T_{d}} a_{d} / T_{d}

(3)

a_{d}

corresponds to the deceleration for each deceleration point.

Acceleration time: total time of the acceleration state.

Deceleration time: total time of the deceleration state.

Idle time: total time of the idle state.

Uniform time: total time of the uniform state.

Acceleration time ratio:

T_{r_a} = T_{a} / T

(4)

Deceleration time ratio:

T_{r_d} = T_{d} / T

(5)

Idle time ratio:

T_{r_i} = T_{i} / T

(6)

Uniform time ratio:

T_{r_u} = T_{u} / T

(7)

Considering the slight speed change during the driving process of the car, the definitions of the uniform speed, idle speed, acceleration, and deceleration are as follows [4,14]:

(1): Uniform: the continuous running state, in which | $a_{i}$ | ≤ 0.15 m/s² and $v_{i} \neq 0$ .
(2): Idle: the continuous state in which $v_{i} = 0$ and the engine is working.
(3): Acceleration: the continuous running state, with $a_{i}$ > 0.15 m/s².
(4): Deceleration: the continuous running state, with $a_{i}$ < −0.15 m/s².

3. Principal Component Analysis

There is a linear relationship between some parameters, so if all 15 characteristic parameters are used, it will cause overlap, redundancy of information, and wasted computing resources. Therefore, principal component analysis is used to reduce the dimensionality of the characteristic parameters, which reduces the complexity of data processing and improves the effectiveness of the computing.

The brief process of principal component analysis is as follows:

(1): Standardized processing.

The data matrix A of the characteristic parameters is

A = [\begin{matrix} a_{11} & a_{12} & \dots & a_{1 n} \\ a_{21} & a_{22} & \dots & a_{2 n} \\ \dots & \dots & a_{i j} & \dots \\ a_{m 1} & a_{m 2} & \dots & a_{m n} \end{matrix}]

(8)

The normalized matrix

X

X = [\begin{matrix} x_{11} & x_{12} & \dots & x_{1 n} \\ x_{21} & x_{22} & \dots & x_{2 n} \\ \dots & \dots & x_{i j} & \dots \\ x_{m 1} & x_{m 2} & \dots & x_{m n} \end{matrix}]

(9)

where

x_{i j} = (a_{i j} - {\bar{a}}_{j}) / \sqrt{Var (a_{j})}

(10)

{\bar{a}}_{j} = \sum_{i = 1}^{m} a_{i j} / m

(11)

V a r (a_{j}) = \sum_{i = 1}^{m} {(a_{i j} - {\bar{a}}_{j})}^{2} / (m - 1)

(12)

(2): The sample correlation coefficient matrix R is calculated.

R = [\begin{matrix} r_{11} & r_{12} & \dots & r_{1 n} \\ r_{21} & r_{22} & \dots & r_{2 n} \\ \dots & \dots & r_{i j} & \dots \\ r_{n 1} & r_{n 2} & \dots & r_{n n} \end{matrix}]

(13)

r_{i j} = \frac{C o v (x_{i}, x_{j})}{\sqrt[]{V a r (x_{1})} \sqrt[]{V a r (x_{2})}} = \frac{\sum_{k = 1}^{m} (x_{k i} - {\bar{x}}_{i}) (x_{i j} - {\bar{x}}_{j})}{\sqrt[]{\sum_{k = 1}^{m} {(x_{k i} - {\bar{x}}_{i})}^{2}} \sqrt[]{\sum_{k = 1}^{m} {(x_{i j} - {\bar{x}}_{j})}^{2}}} (n > 1)

(14)

{\bar{x}}_{i} = \sum_{i = 1}^{m} x_{i j} / m

(15)

{\bar{x}}_{j} = \sum_{i = 1}^{n} x_{i j} / n

(16)

(3): The characteristic eigenvalues and the characteristic eigenvectors are calculated.

The characteristic eigenvalues of the correlation coefficient matrix R

(λ_{1} \geq λ_{2} \geq \dots \geq λ_{p} \geq 0)

and the corresponding characteristic eigenvectors

e_{1}, e_{2}, \dots, e_{n}

, in which,

e_{j} = {(e_{1 j}, e_{2 j}, \dots, e_{n j})}^{T}

, are calculated, and the new parameters calculated with characteristic eigenvectors can be expressed as

\{\begin{matrix} y_{1} = e_{11} x_{1} + e_{21} x_{2} + e_{m 1} x_{m} \\ y_{2} = e_{12} x_{1} + e_{22} x_{2} + e_{m 2} x_{m} \\ \dots \\ y_{n} = e_{1 n} x_{1} + e_{2 n} x_{2} + e_{m n} x_{m} \end{matrix}

(17)

where

x_{m}

is the original component,

y_{n}

is called the

n_{t h}

principal component.

(4): The p principal components are selected and the contribution rate is calculated:

Contribution rate = \frac{λ_{i}}{\sum_{i = 1}^{p} λ_{i}}

(18)

The larger the contribution rate, the stronger the information about the original variable contained in the principal component is. The number of principal components is mainly determined by the cumulative contribution rate. In this study, the cumulative contribution rate is selected to exceed 80% to ensure that the comprehensive variables can include the greatest amount of information for the original variables. The cumulative contribution value of the 15 characteristic parameters extracted in this study is shown in Table 2.

After the dimensionality reduction process with principal component analysis, the cumulative contribution rate of the first three principal components reaches 82.96%, which means these three principal components can represent the original 15 characteristic parameters. The correlation coefficients between the first five principal components and the 15 feature parameters are shown in Table 3.

The greater the absolute value of the correlation coefficient is, the greater the representation of the principal component is to this parameter. Table 4 displays the parameters of each principal component. The first principal component mainly reflects the average speed, maximum speed, maximum acceleration, average acceleration, maximum deceleration, average deceleration, acceleration time ratio, and deceleration time ratio. The second principal component mainly reflects the acceleration time, deceleration time, uniform time, and total time. The third principal component mainly reflects the idle time, idle time ratio, and uniform time ratio.

4. Hierarchical Cluster Method

The hierarchical cluster method is not sensitive to noise in the data, and the number of clusters can be selected according to the actual need in the research, which is suitable for building various driving cycles for different areas. The principle of the hierarchical cluster method is to regard each sample as a separate class and select two classes with the smallest distance to merge into a new class with the condition of specifying the distance between classes. Then the distance between the new class and other classes is calculated, and then the merging process is repeated. One category at a time is decreased until the given termination condition is met.

The general steps for the hierarchical cluster method are listed below:

(1): Each sample is set as an initial class.
(2): The inter-class distance matrix is calculated, and the two classes with the closest distance are combined into a new class. The inter-class distance is calculated as follows:

D_{G} (p, q) = \frac{1}{l k} \sum d_{i j}

(19)

In the formula:

D_{G} (p, q)

is the inter-class distance.

l

and

k

are the number of samples in class

G_{P}

and class

G_{q}

, respectively, and

d_{i j}

is the distance between the

i

sample in class

G_{P}

and the

j

sample in class

G_{q}

(3): Step (2) is repeated.
(4): The process is ended when the termination condition is met.

According to the local traffic law, the traffic condition is generally divided into three categories: low speed (low average running speed, frequent acceleration, and deceleration), medium speed (average running speed is moderate, acceleration and deceleration are frequent, normal deceleration), and high speed (high average running speed, large acceleration ratio). Therefore, the termination condition of the hierarchical cluster method in this study is the acquisition of three final classes. Finally, three clusters are obtained, and the number of micro-trips in each category is 572, 956, and 59, respectively.

5. Construction and Verification of Vehicle Driving Cycle

5.1. Construction of Vehicle Driving Cycles

S is set as the average relative error of the characteristic parameters of the constructed driving cycles and the real-world data:

S = (S_{1} + S_{2} + \dots + S_{j}) / j

(20)

In the expression,

S_{j}

is the relative error of the

j_{t h}

characteristic parameter. The calculation formula is as follows:

S_{j} = \frac{|a_{a v g_c j} - a_{r j}|}{a_{r j}}

(21)

In the formula:

a_{a v g_c j}

is the average value of the

j_{t h}

characteristic parameter of all the micro-trips of the constructed driving cycle, and

a_{r j}

is the

j_{t h}

characteristic parameter of the real-world data.

The reasonable duration used for most urban driving cycles is about 1200 s [15], so according to the proportion of the number of samples in each category, 10 samples with the smallest S values are selected and combined into a driving cycle. The values of characteristic parameters of the 10 micro-trips are shown in Table 5, and the constructed driving cycle is shown in Figure 3. The acceleration time, deceleration time, idle time, uniform time, and the total time of real-world data are the summary values of all micro-trips, so these parameters are not included in the error calculation.

5.2. Construction of Vehicle Driving Cycles

The economy is one of the important indicators of electric vehicles, and the power consumption per 100 km is often used as an evaluation index. In addition, overall efficiency is a commonly used evaluation index. In this study, the constructed driving cycle is used as the test cycle condition for the economy simulation, and the NEDC is also used for comparison. The simulation results are shown in Figure 4:

It can be seen from Figure 4 that both the power consumption per 100 km and the overall efficiency are closer for the construction driving cycles and the actual driving data. The relative errors are 0.23% and 3.6%, while the NEDC relative errors are 7.47% and 41.34%.

Furthermore, the motor’s efficiency at each second is also calculated, and the frequency of the efficiency is added, as seen in Figure 5.

In Figure 5, the scale on the horizontal axis represents the range of efficiency; for instance, “0.1” represents the efficiency from 0.1 to 0.2. Figure 5 indicates that the efficiency distribution for the constructed driving cycle is closer to real-world data in most cases. This is why the electric vehicle has a more realistic economy test result in the constructed driving cycle.

5.3. Statistical Verification of Constructed Driving Cycle

The reason the electric vehicle performs more “accurately” in the constructed driving cycle than in NEDC is the difference in the statistical characteristics of the driving cycles. The comparison of the characteristic parameters of different conditions is shown in Table 6.

The calculation results indicate that the overall relative error of the characteristic parameters of the constructed driving cycle and the real-world data is 5.46% (<10%). In comparison, the NEDC is 36.8%, which reveals why the constructed driving cycle is more accurate.

The VA matrix reflects the statistical distribution law of velocity acceleration and is often used for building driving cycles. From the perspective of the VA matrix, the constructed driving cycle is also closer to the real-world data. The 3D plots of the VA probability distributions of the real-world data constructed driving cycle and NEDC are shown in Figure 6, Figure 7 and Figure 8.

Figure 8 indicates that the v-a distribution of the NEDC is simple and centralized, which is very different from the real-world data. However, the VA matrix of a constructed driving cycle covers most values for real-world data, which is because NEDC covers few speeds and acceleration. For the matrices of the constructed driving cycle and real-world data, most of the acceleration from −1.4 m/s² to +1.2 m/s² appears in the low-speed zone. Additionally, the peak area distributions of the two matrices are highly similar, which means that the constructed driving cycle is more representative of the real traffic conditions.

6. Summary

This study proposes a method based on principal component analysis and the hierarchical cluster method for constructing electric vehicle driving cycles. First, the data collected with the onboard data acquisition terminal are denoised using the wavelet data denoising principle and divided into micro-trips. The characteristic parameters of the micro-trips are extracted. The principal component analysis is then applied to reduce the dimensions. Next, the hierarchical cluster method is used to classify the micro-trips into three categories, and 10 micro-trips are selected to combine the driving cycle according to the proportion of each category. Finally, the generated driving cycle is verified based on the aspects of both the electric vehicle economy test and statistical analysis.

The economy test verification shows that the constructed driving cycle has a 0.23% error in power consumption per 100 km, a 3.6% error in overall efficiency, and a closer distribution in terms of the motor efficiency frequency, which indicates that the tested electric vehicle performs similarly in both driving cycle. The statistical analysis shows that the overall relative error of the characteristic parameters of the constructed driving cycle and the real-world data is 5.46%, and the built driving cycle has a similar VA matrix to the real-world data, which both indicates that the constructed driving cycle can represent real-world data for related research. The verification results prove the validity and feasibility of the methodology proposed in this study, and this study provides a new method for the development of driving cycles.

Author Contributions

Conceptualization, Z.J. and T.W.; methodology, Z.J. and T.W.; software, T.W.; validation, Z.J. and T.W.; formal analysis, Z.J. and S.Z. investigation, Z.J. and T.W.; resources, Z.J. and T.W.; data curation, Z.J. and S.Z.; writing—original draft preparation, Z.J.; writing—review and editing, Z.J.; visualization, Z.J. and T.W.; supervision, G.W.; project administration, S.Z. and G.W.; funding acquisition, G.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National High Technology Research and Development Program of China, grant number 2011AA11A286.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

Komorska, I.; Puchalski, A.; Niewczas, A.; Ślęzak, M.; Szczepański, T. Adaptive driving cycles of EVs for reducing energy consumption. Energies 2021, 14, 2592. [Google Scholar] [CrossRef]
Tong, H.Y.; Hung, W.T. A framework for developing driving cycles with on-road driving data. Transp. Rev. 2010, 30, 589–615. [Google Scholar] [CrossRef]
Huang, D.; Xie, H.; Ma, H.; Sun, Q. Driving cycle prediction model based on bus route features. Transp. Res. Part D Transp. Environ. 2017, 54, 99–113. [Google Scholar] [CrossRef]
Shi, S.; Lin, N.; Zhang, Y.; Cheng, J.; Huang, C.; Liu, L.; Lu, B. Research on Markov property analysis of driving cycles and its application. Transp. Res. Part D Transp. Environ. 2016, 47, 171–181. [Google Scholar] [CrossRef]
Lin, J.; Niemeier, D.A. An exploratory analysis comparing a stochastic driving cycle to California’s regulatory cycle. Atmos. Environ. 2002, 36, 5759–5770. [Google Scholar] [CrossRef]
Berzi, L.; Delogu, M.; Pierini, M. Development of driving cycles for electric vehicles in the context of the city of Florence. Transp. Res. Part D Transp. Environ. 2016, 47, 299–322. [Google Scholar] [CrossRef]
Brady, J.; O’Mahony, M. Development of a driving cycle to evaluate the energy economy of electric vehicles in urban areas. Appl. Energy 2016, 177, 165–178. [Google Scholar] [CrossRef]
Nouri, P.; Morency, C. Evaluating Microtrip Definitions for Developing Driving Cycles. Transp. Res. Rec. J. Transp. Res. Board 2017, 2627, 86–92. [Google Scholar] [CrossRef]
Zhang, Y.X.; Su, X.H.; Gao, G.K.; Shang, Y. Driving cycles with improved principal component and K-means clustering algorithm. Sci. Technol. Eng. 2021, 21, 3199–3250. [Google Scholar]
Gao, J.P.; Ren, D.X.; Xi, J.G. Construction of vehicle driving cycles based on global K-means clustering algorithm. J. Henan Univ. Technol. Nat. Sci. Ed. 2019, 38, 112–118. [Google Scholar] [CrossRef]
Yuan, Y.M.; Liu, H.Z.; Li, H.S. Improved K-Means text clustering algorithm based on density peaks and its parallelization. J. Wuhan Univ. Sci. Ed. 2019, 5, 457–464. [Google Scholar] [CrossRef]
Chen, J.Y. Research on Construction of Tracked Vehicle Driving Conditions and Performance Matching of Transmission. Zhejiang University. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbname=CMFD202002&filename=1020725586.nh (accessed on 1 March 2020).
Hung, W.T.; Tong, H.Y.; Lee, C.P.; Ha, K.; Pao, L.Y. Development of a practical driving cycle construction methodology: A case study in Hong Kong. Transp. Res. Part D Transp. Environ. 2007, 12, 115–128. [Google Scholar] [CrossRef]
Arun, N.H.; Mahesh, S.; Ramadurai, G.; Shiva Nagendra, S.M. Development of driving cycles for passenger cars and motorcycles in Chennai, India. Sustain. Cities Soc. 2017, 32, 508–512. [Google Scholar] [CrossRef]
Galgamuwa, U.; Perera, L.; Bandara, S. A representative driving cycle for the southern expressway compared to existing driving cycles. Transp. Dev. Econ. 2016, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]

Figure 1. On-board data acquisition terminal.

Figure 2. An example of a denoised data period.

Figure 3. Constructed driving cycle.

Figure 4. Comparison of power consumption per 100 km.

Figure 5. The frequency of the efficiency of the motor.

Figure 6. VA matrix of real-world data.

Figure 7. VA matrix of constructed driving cycle.

Figure 8. VA matrix of NEDC.

Table 1. Micro-trip characteristic parameters.

	Symbol	Characteristic Parameters
1	$V_{a v g}$	Average speed (km/h)
2	$V_{m a x}$	Maximum speed (km/h)
3	$a_{m a x_a}$	Maximum acceleration (m/s²)
4	$a_{a v g_a}$	Average acceleration (m/s²)
5	$a_{m a x_d}$	Maximum deceleration (m/s²)
6	$a_{a v g_d}$	Average deceleration (m/s²)
7	$T_{a}$	Acceleration time (s)
8	$T_{d}$	Deceleration time (s)
9	$T_{i}$	Idle time (s)
10	$T_{u}$	Uniform time (s)
11	$T$	Total time (s)
12	$T_{r_a}$	Acceleration time ratio (%)
13	$T_{r_d}$	Deceleration time ratio (%)
14	$T_{r_i}$	Idle time ratio (%)
15	$T_{r_u}$	Uniform time ratio (%)

Table 2. Cumulative contribution value.

Component	Total	% of Variance	Cumulative (%)
1	7.1383	47.59	47.59
2	3.5963	23.98	71.56
3	1.7092	11.39	82.96
4	0.7362	4.91	87.87
5	0.6491	4.33	92.19
6	0.4340	2.89	95.09
7	0.3714	2.48	97.56
8	0.1551	1.03	98.60
9	0.0845	0.56	99.16
10	0.0605	0.40	99.56
11	0.0438	0.29	99.86
12	0.0184	0.12	99.98
13	0.0031	0.02	100.00
14	0.0000	0.00	100.00
15	0.0000	0.00	100.00

Table 3. Correlation coefficient between principal components and characteristic parameters.

Characteristic Parameters	Principal Component
Characteristic Parameters	1	2	3
Average speed	−0.32	0.10	0.11
Maximum speed	−0.35	−0.01	0.01
Maximum acceleration	−0.33	−0.09	−0.05
Average acceleration	−0.29	−0.23	−0.10
Maximum deceleration	0.33	0.12	0.05
Average deceleration	0.29	0.23	0.09
Acceleration time	−0.24	0.38	−0.04
Deceleration time	−0.23	0.39	−0.04
Idle time	0.01	0.03	−0.56
Uniform time	−0.15	0.45	−0.08
Total time	−0.18	0.44	−0.16
Acceleration time ratio	−0.30	−0.21	0.12
Deceleration time ratio	−0.29	−0.20	0.12
Idle time ratio	0.17	−0.04	−0.62
Uniform time ratio	0.16	0.26	0.45

Table 4. Feature parameters represented by five comprehensive feature parameters.

Principal Component	Parameters
1	Average speed, Maximum speed, Maximum acceleration, Average acceleration, Maximum deceleration, Average deceleration, Acceleration time ratio, Deceleration time ratio
2	Acceleration time, Deceleration time, Uniform time, Total time
3	Idle time, Idle time ratio, Uniform time ratio

Table 5. Characteristic parameters of micro-trips.

	V_avg	V_max	a_{max_a}	a_{avg_a}	a_{max_d}	a_{avg_d}	T_{r_a}	T_{r_d}	T_{r_i}	T_{r_u}
1	11.57	44.74	0.7	0.34	−0.71	−0.48	20.13	15.72	50.94	13.21
2	0.18	0.53	0	0	0	0	0	0	53.85	46.15
3	26.34	60.02	0.91	0.38	−0.81	−0.37	23.66	25.27	29.03	22.04
4	18.89	45.29	0.83	0.42	−0.6	−0.32	26.12	34.33	22.39	17.16
5	4.21	17.78	0.56	0.35	−0.72	−0.31	13.82	15.13	48.03	23.03
6	16.47	49.81	0.96	0.52	−0.76	−0.38	13.33	17.78	18.33	50.56
7	33.82	60.21	1.15	0.58	−0.9	−0.33	20.86	33.09	7.19	38.85
8	30.16	58.10	0.8	0.44	−1.36	−0.79	50.7	28.17	7.04	14.08
9	0.21	0.46	0	0	0	0	0	0	45	55
10	36.46	58.73	1.29	0.47	−1.22	−0.49	28.8	28	8	35.2

Table 6. Comparison of characteristic parameters.

Characteristic Parameters	Real-World Data	Constructed Driving Cycle		NEDC
Characteristic Parameters	Real-World Data	Characteristic Parameters	Relative Error (%)	Characteristic Parameters	Relative Error (%)
V_avg	16.54811	17.831	7.75	19.35	16.93
V_max	36.09808	39.56778	9.61	39.08	8.26
a_{max_a}	0.698444	0.72	3.09	0.81	15.97
a_avg	0.353232	0.35	0.92	0.74	109.49
a_{max_d}	−0.74004	−0.708	4.33	−0.87	17.56
a_{avg_d}	−0.37009	−0.347	6.24	−0.79	113.46
T_{r_a}	21.16277	21.242	0.37	19.79	6.49
T_{r_d}	20.17054	21.696	7.56	17.03	15.57
T_{r_i}	22.17403	23.294	5.05	32.44	46.30
T_{r_u}	37.43831	33.8	9.72	30.73	17.92

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Jing, Z.; Wang, T.; Zhang, S.; Wang, G. Development Method for the Driving Cycle of Electric Vehicles. Energies 2022, 15, 8715. https://doi.org/10.3390/en15228715

AMA Style

Jing Z, Wang T, Zhang S, Wang G. Development Method for the Driving Cycle of Electric Vehicles. Energies. 2022; 15(22):8715. https://doi.org/10.3390/en15228715

Chicago/Turabian Style

Jing, Zhecheng, Tianxiao Wang, Shupei Zhang, and Guolin Wang. 2022. "Development Method for the Driving Cycle of Electric Vehicles" Energies 15, no. 22: 8715. https://doi.org/10.3390/en15228715

APA Style

Jing, Z., Wang, T., Zhang, S., & Wang, G. (2022). Development Method for the Driving Cycle of Electric Vehicles. Energies, 15(22), 8715. https://doi.org/10.3390/en15228715

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu