Research Article  Open Access
Kai Zhang, Jian Ma, Xuan Zhao, Xiaodong Liu, Yixi Zhang, "Parameter Identification and State of Charge Estimation of NMC Cells Based on Improved Ant Lion Optimizer", Mathematical Problems in Engineering, vol. 2019, Article ID 4961045, 18 pages, 2019. https://doi.org/10.1155/2019/4961045
Parameter Identification and State of Charge Estimation of NMC Cells Based on Improved Ant Lion Optimizer
Abstract
For lithium battery, which is widely utilized as energy storage system in electric vehicles (EVs), accurate estimating of the battery parameters and state of charge (SOC) has a significant effect on the prediction of energy power, the estimation of remaining mileage, and the extension of usage life. This paper develops an improved ant lion optimizer (IALO) which introduces the chaotic mapping theory into the initialization and random walk processes to improve the population homogeneity and ergodicity. After the elite (best) individual is obtained, the individual mutant operator is conducted on the elite individual to further exploit the area around elite and avoid local optimum. Then the battery model parameters are optimized by IALO algorithm. As for the SOC estimation, unscented Kalman filter (UKF) is a common algorithm for SOC estimation. However, a disadvantage of UKF is that the noise information is always unknown, and it is usually tuned manually by “trialanderror” method which is irregular and timeconsuming. In this paper, noise information is optimized by IALO algorithm. The singular value decomposition (SVD) which is utilized in the process of unscented transformation to solve the problem of the covariance matrix may lose positive definiteness. The experiment results verify that the developed IALO algorithm has superior performance of battery model parameters estimation. After the noise information is optimized by IALO, the UKF can estimate the SOC accurately and the maximum errors rate is less than 1%.
1. Introduction
With the increasingly serious problems of global energy crisis and environmental pollution, electric vehicles (EVs) have experienced explosive development and drawn increasing attention from institute and industry in the last decade [1–3]. The energy storage system is the most important part in EVs which is often comprised by different kinds of batteries, including leadacid batteries, nickelhydrogen batteries, and secondary lithium batteries. Lithium battery is the most popular energy source nowadays and occupies more than 90% of the market share due to the characteristic of high energy density, high efficiency, no memory effect, long cycle life, and environmentfriendly performance [4–6]. Due to the flammability of liquid electrolyte and the lithium dendrite growth which may puncture battery separator and cause shortcircuit, lithium battery may degrade rapidly even thermal runaway under abuse condition, like overcharge or over discharge. Therefore, lithium battery should be operated stable and safe, which is achieved by advanced Battery Management System (BMS) [7]. One of the most important functions of BMS is state of charge (SOC) estimation, which can improve the control efficiency of energy and extend the usage life of the batteries [8]. However, SOC cannot be obtained directly by measuring and must be deduced from other measurable signals, like voltage and current. Also, several factors can influence the accuracy of SOC estimation and should be thoroughly investigated [9].
Many studies have focused on the research of SOC estimation strategies. The most widely used strategy is coulomb counting which is based on load current integration [10–13]. However, this strategy is openloop without feedback and may easily cause cumulative calculation error due to inaccurate initial SOC value and stochastic disturbances. The other one popular technology is intelligent algorithm control strategies including neural network and machine learning. This method is not interested in the battery internal state and only cares about the relationship between the output and input which can extract from experimental data set [14]. However, this method needs a larger amount of experimental data for algorithms training and always is timeconsuming [15–18]. Another dominating method is state estimate algorithms based on battery model, like equivalent circuit models (ECMs) and electrochemical models [19–26]. This method is often utilized in BMS due to the easy implementation, high accuracy, and high efficiency. As for the battery models, ECMs are more suitable for online deployment than electrochemical models since ECMs utilize electrical parameters like current, voltage, and temperature, which are easier to obtain online than the parameters used in electrochemical models (e.g., electrolyte concentration and diffusion coefficients). Among the ECMs modelbased methods, the two widely used SOC estimation techniques are extended Kalman filter (EKF) and unscented Kalman filter (UKF) [27–32]. However, the noise information (covariance matrixes) of the state space should be known to apply EKF and UKF algorithms, which is often obtained manually by “trialanderror” method and this may be difficult and timeconsuming. Some researches proposed adaptive technique to calculate the covariance matrixes based on the estimation error [28, 31]. However, this method may cause the covariance matrix becoming singular/negative and shut down the estimation process. Another important aspect of SOC estimation is the accuracy of battery model. The battery model parameters are always obtained by recursive least square (RLS) algorithm [7]. However, the results of RLS are always fluctuation and highly affected by system noise.
This paper introduces natureinspired metaheuristic algorithm to find the optimal value of battery model parameters and the covariance matrixes. The key features of metaheuristic algorithms are randomness and the utilization of stochastic operators when searching in the solution spaces [33]. Some of the most popular metaheuristic algorithms contain particle swarm optimization (PSO) [34], simulated annealing algorithm (SA) [35], ant colony optimization (ACO) [36], genetic algorithm (GA) [37], firefly algorithm [38], and BAT Search Algorithm [39]. Ant lion optimizer (ALO), proposed by Mirjalili, is imitating the unique predation behavior of ant lions to find the optimal value [40, 41]. The performance of ALO in solving optimization problem has been vindicated and only a few parameters are needed. E.S. Ali uses the ALO to find the optimal positions and determine the sizing of Distributed Generation (DG) in various distribution systems [42]. However, ALO has a drawback of early prematurity and local optimum, especially in solving complicated problems [43]. To solving this problem, different methods have been combined with ALO to enhance the performance. Tian applied PSO after the solution space has been searched by ALO, which can combine the search ability of ALO and PSO [43]. Zhongqiang Wu applied chaotic sequence for generating the initial position of population (ants and antlions) to increase the uniformity and ergodicity [44].
Several improvements are introduced into ALO algorithm to enhance the optimization performance. Chaotic mapping theory is introduced in the initialization and random walk processes. Then after the elite (best) individual is found by ALO in one iteration, individual mutant operator based on chaotic mapping is conducted on every elements of the elite individual, which can further exploit the region around elite. Based on the proposed improved ant lion optimization (IALO), the parameters of battery model are optimized by IALO. Also, the noise covariance matrixes in UKF are optimized by IALO. Another small improvement of this paper is that singular value decomposition (SVD) is utilized in UKF to solve the problem of the covariance matrix which may become singular and lose positive definiteness.
The remainder of the paper is arranged as follows: Section 2 introduces the battery model, the parameters estimation method, and the SOC estimation algorithm. Section 3 gives a brief review of the antlion algorithm. In Section 4, the improved ALO is proposed. Section 5 describes the experimental results and discussion. Section 6 provides the conclusion.
2. Battery Model and SOC Estimate Algorithm
2.1. Battery Model
Equivalent circuit model, electrochemical model, and black box model are three widely used models for modeling sophisticated battery behaviors. The three models have their own characteristic in complexity, accuracy, and computational cost [45]. Considering the hardware limits and realtime control of EVs, including enough accuracy and reasonable complexity, equivalent circuit model is widely used in EVs application. Many kinds of equivalent circuit models have been studied in the past. Twelve equivalent circuit models have been studied in [46] which include zerostate hysteresis, onestate hysteresis model, selfcorrecting (ESC) model, firstorder RC model, secondorder RC model, thirdorder RC model, and so on. Their experiments and comparison results demonstrated that, considering both model accuracy and model robustness, the firstorder RC model, which is also called Thevenin model, is preferred for lithium nickelmanganesecobalt oxide (LiNMC) cells. Therefore, the Thevenin battery model is adopted in this paper and demonstrated in Figure 1. The battery model consists of a voltage component which represents the open circuit voltage (OCV), a resistor component , which denote the ohmic resistance, and a RC network component, which describes the concentration polarization and diffusion processes during battery operation, the resistance element denotes the polarization resistance, and the capacitance denotes the polarization capacitance. The electrical characteristic of the Thevenin model is calculated aswhere is the load current (assume greater than 0 for discharging) and represents the terminal voltage. is the polarization voltage over the polarization resistance, and the represents the OCV and can be modeled with function which is a simple electrochemical model [7], and can be calculated aswhere are multinomial coefficients and they need to be fitted based on experimental data and z denotes the state of charge (SOC) of battery.
2.2. Parameters Identification Method
For iteration calculation purpose, the discrete form of state space equations can be written aswhere is fixed sample interval and assumed to be 1s in this study. is the time constant which equals . k is an integer variable which denotes the discrete index. is the kth sampling time of load current. and are the terminal voltage and open circuit voltage at kth sampling moment, respectively. is the polarization voltage at kth sampling moment.
The load current and the terminal voltage are directly measured. The open circuit voltage is calculated based on SOC which is integral of load current from a known initial SOC value. The model parameters, including , , and , are identified by the improved ant lion optimization (IALO) algorithm which will be discussed in Sections 3 and 4. The optimization objective function is to minimize the quadratic sum of the voltage estimation error.
2.3. SOC Estimation Based on Singular Value Decomposition (SVD) UKF
After the parameters of battery model are identified by IALO, SOC is estimated by UKF (unscented Kalman filter). The state space equations of battery model are nonlinear. According the widely used EKF SOC estimation algorithm, the nonlinear state equations must be linearization, and the oftenused linearization method is based on Taylor expansion which may lead into truncation error. The UKF algorithm does not need to linearize the model and apply nonlinear unscented transformation (UT) to generate sigma points from the original state. The sigma points are calculated based on known mean and covariance and then substituted into the state function to obtain the new sigma points, then updating mean and covariance matrixes by the statistical characteristics of new sigma points.
The SOC denotes the proportion of the residual charge in battery to the rated capacity. The discrete form of SOC function is calculated bywhere is the SOC at kth sampling moment, is the rated capacity, and is the coulomb efficiency and related to current rate, battery chemistry, and temperature. For simplicity, the coulomb efficiency is set as .
To utilize the UKF algorithm for SOC estimation, the nonlinear battery dynamic state equation and measurement equation are described aswhere denotes the state vector and the system input vector is the measured load current . The output vector is the terminal voltage , and represent process noise signal and measurement noise signal, and their covariance matrixes are and , respectively.
The very important part of UKF nonlinear unscented transformation (UT) is calculated as follows:where is the sigma points, and are the weight coefficients of mean and covariance, respectively, n is dimensional of the state vector , is a scale adjustment parameter which can adjust the distribution of sigma points and enhance accuracy, and are set to 1 and 2 in this paper, and is normally set to 0 or 3n. is the Cholesky decomposition of covariance matrix .
There is a problem existing in Cholesky decomposition, namely, if the covariance matrix is singular, the Cholesky decomposition cannot be carried out, and the iterative calculation would be stopped. Moreover, the Cholesky decomposition is very sensitive to the calculation error, and the finite word length or rounding error in the computational process will often make the covariance matrix lose its symmetry or positive definiteness, thus resulting in the problem of numerical instability of the filtering and reducing the accuracy of the filtering algorithm. In order to solve this problem, we introduce singular value decomposition (SVD) of matrix into the generating of sigma points, which can enhance the numerical stability of iterative updating state covariance matrix and decompose covariance matrix even morbidity occurs in matrix.
According to the singular value decomposition theory, a matrix can be decomposed by orthogonal matrix and aswhere is the eigenvalue of matrix and called singular value, . When matrix is symmetric positive definite matrix, the SVD of matrix is rewritten as
The SVD of covariance matrix is calculated aswhere is covariance matrix of last iteration k1th and then the sigma points at kth iteration are calculated by
The detailed SVDUKF algorithm process is shown in Table 1. In the SVDUKF algorithm process, the process noise covariance , and measurement noise covariance matrix are often considered as constant value. The and are always decided by trialanderror procedure which is arbitrary and timeconsuming. In this study, the improved ant lion optimization (IALO) algorithm is utilized to search the optimal value of and which will be discussed in Sections 3 and 4. The optimization objective function is to minimize the quadratic sum of the voltage estimation error and SOC estimation error with different weights.

3. Introduction of Ant Lion Optimization Algorithm
Ant lion optimizer (ALO), a recently proposed heuristics natureinspired algorithm, imitates the foraging behavior of ant lion’s larvae [40]. The ant lions belong to the genus of Myrmeleontidae, carnivorous, and feed on other insects. The larvae of ant lions live under dry ground and create funnelshaped traps in sandy soil to trap prey. After the trap is built, the ant lion larvae wait the prey (mainly ants) by hiding underneath the bottom of the trap. In natural environment, ants move randomly and casually; once ants move into pits, the ant lion will try to catch ants. However, ants usually may try to creep out the trap. In this situation, ant lions will cast sands towards the edge of the trap to make the trap side steeper and slide the ant down to the bottom of the trap and grab it. After an ant is seized, ant lion will pull the prey under the sand and eat, then throw out the leftovers, and rebuilt the pit for the next prey. The ant lion denotes the optimized solution of the optimization problem and the ant stand for the feasible solution. By catching the ant which has higher fitness, the ant lion can update its position and fitness by the information of ant.
In the search space, there are many ant lions and ants, whose positions are located randomly. To imitate the interaction between ant lions and ants, ants will move over the search space with random walk, which may go through or around the traps. The superior individual/ant lion will build a larger and sharper pit and catch more ants. The random walk steps of ants are influenced by ant lions and elite (best) ant lion. Ants will update their position based on ant lions and elite ant lion to ensure the diversity and optimization capability of the algorithm.
During the optimization, the algorithm abides the following conditions: ants are wandering over the solution space with random walks, and the random walks are influenced by the position of traps of ant lions; pits built by ant lions are proportional to the fitness of themselves, and the better fitness (the bigger pit) has the high probability to catch ants; the step of random walk is accordingly reduced when ant is sliding to the ant lion which is hiding under the bottom of the pit; if the fitness of an ant is better than that of an ant lion, this signifies that the ant is seized and eaten by the ant lion, and the ant lion changes itself with the position of the ant and growing fitter; after each hunting, an ant lion changes its position near the position where ant is caught and set up a new pit to improve its possibility of catching prey.
The random walk, which is chosen for modeling ants’ stochastically moving behavior, is obtained aswhere is the function of cumulative sum, is the biggest number of iterations, denotes the current iteration, and represents a stochastic function calculated bywhere is a random number between (0,1), which subject to uniform distribution.
The number of both ants and ant lions is N, and the dimension of position is D. The positions of ants are recorded in matrix ; represents the value of dimension of position value of the ant; A fitness function is utilized to evaluate each ant and matrix is utilized to record the fitness value of all ants; stores the position of every antlion and represents the value of dimension of position of the antlion; is utilized to evaluate each antlion and matrix is utilized to store the fitness value of all ant lions.
Random walk of ants is based on (12). Ant changes its position by random walk at each iteration during optimization. Considering the search space has its own boundary, in order to keep the position of ant after random walk remain in the search space, the following equation is utilized to normalize and calculate the new position:where is random walk step of ith ant; is the normalized random walk displacement; and are the lowest and highest of random walk displacement of ith ant; and indicate the lowest and highest of random walk displacement of ith ant at current iteration.
As mentioned before, ant lions’ position can influence the random walks of ants. This influence can be mathematically modeled by following equations:where denotes the position of selected jth ant lion at tth iteration and and indicate the lowest and highest of random walk displacement of all ants in the tth iteration. The selected antlion is determined by roulette wheel method. Ant lions throw sands out the trap after they perceive that an ant fall in the pit and slide ant towards the pit center. To model this mechanism, the region of ants’ random walks is shrunk adaptively as the increasing of iteration and calculated by following equations:where I is a ration, and shown as , t denotes the current iteration, T denotes the biggest number of iteration, and w is an index varying with the current iteration (w =2 if , w=3 if , w=4 if , w=5 if , and w=6 if ).
During optimization, the elite is set as the best ant lion and updated every iteration if better ant lion shows up. And the random walk of ants is also influenced by the elite ant lion. Therefore, it is assumed that the random walks of every ant are over the selected ant lion and the elite antlion and calculated bywhere denotes iteration, denotes the random walk over the selected ant lion, denotes the random walk over the elite ant lion, and denotes updated position of ith ant. After random walk, the fitness is calculated. Comparing the fitness of ant lion and ant, assign better one to ant lion, which can also be regarded as the ant which is hunted. The whole procedure of ALO algorithm is demonstrated in Figure 2.
4. Improvements on Ant Lion Optimizer
In the ALO algorithms optimization procedure, the updating of ants’ positions is based on the random walk over the chosen ant lion and the elite ant lion. Save the greatest agent (ant) as elite in each iteration. These can make ALO have low calculating cost, high efficiency, and high convergence rate. However, the ALO has the problem of premature convergence and local optimum when it is used for solving complex optimization problems. Several improvements are inserted into ALO to deal with these problems and further enhance the global optimization ability and accuracy.
4.1. Chaotic Map for Position Initialization and Random Walk
Random function has been used in many optimization algorithms. Chaos function has the same character of longscale indeterminism and irregular property with random function. Chaos can help orders arise from disorder, which is like many biologistic systems inspired optimization algorithms like GA, PSO, and ALO. Considering these common features between chaos function and optimization algorithms, it is natural to think about combining these theories to enhance the capability of the optimizer.
In this paper, to improve the population homogeneity and ergodicity and skip local optima, chaos function is utilized to generate random walk sequences. The random walk, which models the ants’ stochastic movement, is calculated by iterating the steps of chaotic map from a random initial value.
Logistic map, a wellknown chaotic map, has been used in many papers. Logistic map is based on polynomial model, which is a classical evidence of demonstrating how a very simple nonlinear dynamic system can generate complicated behavior [47, 48]. Logistic map is defined as follows:where is the nth chaotic value and n denotes the iteration number. is the control parameters of chaotic map. When , the system is in chaotic status. However, the logistic map demonstrates some shortages such as blank windows, stable windows, narrow key space, and uneven distribution of sequences. To obtain better chaotic behavior and larger key space and overcome this weakness, a 3D intertwining logistic map is adopted in this paper [49, 50] and can be defined by where and are control parameters and , , , and . The intertwining logistic map has better chaotic behavior and more even distribution, without the shortages of blank windows, stable windows. To more clearly demonstrate the uniformity and ergodicity of intertwining logistic map, the chaos sequence is generated by parameter set , , , , , , and , as shown in Figure 3, and the first 100 iterations are removed to avoid the transient effect. From this figure, the chaotic function is efficient and random and chaotic search is easier to get rid of local minimum comparing with other stochastic optimization algorithms [51].
The chaotic value belongs to the interval (0, 1); it needs scale to the solution space by the following equation:where is the updated position after random walk; and are the lowest and biggest value of search space, respectively.
4.2. Individual Mutation Based Chaotic Mapping
After the elite individual is obtained, to further explore of the solution space and exploit of the small region around elite, chaotic mapping is used to mutate every dimension of the elite ant lion with a new value. This approach can enhance the stochastic property and increase the variety of the individual. The process is summarized below.
Initialization. The obtained elite individual is assumed as ; denotes the dimension and is set as 3 in this paper. Set the maximum number of chaotic mutations as .
Start Iteration. Convert every element of elite individual to a normalized value in the range (0, 1).where and are maximum and minimum of the jth element, respectively.
Outer Loop . Generate new chaotic value through 3D intertwining logistic map theory as follows:
Inner Loop . Invert the chaotic value into the actual position as follows:
Replacing with , the new individual is .
Calculate the fitness and replace with when is better than .
End Inner Loop. Form the new individual . Then comparing the fitness of and , assign the bigger one to .
End Outer Loop. By means of individual mutant with chaotic mapping theory, there is better probability for the elite individual to skip over the local optimum and find a greater solution.
According to the above analysis of improvements on the ant lion optimization algorithm, the improved ant lion optimizer (IALO) procedure is demonstrated in Figure 4. In the following section, IALO algorithm will be used to optimize the value of battery model parameters. The state noise covariance matrix Q and measurement noise covariance matrix R values are also optimized by IALO algorithm.
5. Experiment and Results Discussion
In this section, to validate the performance of the developed IALOSVDUKF SOC estimation method, experiments on LiNiMnCoO2/Graphite lithiumion cells were conducted under DST, FUDS, US06, and BJDST, respectively. Section 5.1 presents a brief introduction of the experiments. The results are discussed in Section 5.2.
5.1. Experiment Detail
The test bench includes an Arbin BT2000 battery tester, a thermal chamber with temperature controllable, a PC with Arbin Software to control the test system, and monitor and storage test data. 18,650 LiNiMnCoO2/Graphite (NMC) lithiumion cells are tested in the experiment and the basic information is shown in Table 2. The experiment ambient temperature is set to 25°C by using the thermal chamber. Sampling interval is set to 1s during the test. The test results are available at https://web.calce.umd.edu/batteries/data.htm.

The pulsecurrent tests are conducted to obtain the OCVSOC function. The cell is first fully charged to 100% SOC by C/2 constant current charge phase and 4.2V constant voltage charge phase with a C/20 cutoff current. Second, the cell is discharged an amount of charge (10% SOC) with C/2, followed by 2 hours relaxation period to achieve equilibrium state inside the cell. The pulsecurrent, relaxation period is repeated every 10% SOC until reaching the lower cutoff voltage. Then, the cell is charged by following the similar steps until reaching the upper cutoff voltage: a constant voltage charge with a C/20 cutoff current is followed to make sure the cell is fully charged. Finally, the data from charge phase and discharge phase are averaged and nonlinear least square (LS) is applied to identify parameters of OCVSOC function which is calculated by (2) and the results are given in Table 3 and Figure 5. As can be seen from the figure, the OCV function can follow the experiment data accurately.

Several widely used cycle conditions are conducted in the experiment, including DST (dynamic stress test) [52], FUDS (Federal Urban Driving Schedule) [53], BJDST (Beijing Dynamic Stress Test) [54], and US06 (US6 Highway Driving Schedule) [55]. DST is the simplest condition and is utilized for identifying the battery parameters and the noise information. To validate the battery model and SOC estimation algorithm, other more complicated cycle conditions are utilized for SOC estimation.
5.2. Results and Discussion
(1) Certification of the Performance of IALO. Six widely used functions are selected to certify the optimization ability of IALO algorithm and compared with Genetic Algorithm (GA), particle swarm algorithm (PSO), and ALO algorithms. The test functions are shown in Table 4. Griewank function has many local extremum and multiple peaks. Ackley function has many local extremums around global extremum points in a narrow area. The global optimum of Rosenbrock function is located in a smooth, long, and narrow parabolic valley. Rastrigin function has many widespread local extremums. Branin function has 3 global minimum points and surrounded by many local extremums. SixHump Camelback function has six local extremums, and two of them are global optima. The parameters of GA algorithm are as follows: the population size is 100, the crossover probability is 0.7, and the mutation probability is 0.3. The parameters of PSO algorithm are as follows: the learning coefficients c1 and c2 are 1.5, and the population size is 200. The parameters of ALO algorithm are as follows: the number of search agents (ants) is 40, which is also the number of ant lions. The parameters of IALO algorithm are as follows: the number of ants or ant lions is the same as that of ALO, and the parameters of intertwining logistic map are , , , and .
 
D is set to 2 in this paper. 
The performance of optimal value search of GA, PSO, ALO, and IALO algorithm is demonstrated in Table 5. In the table, optimum value is the best optimal value that the algorithm can find. Convergence iteration represents the iteration when the optimum value keeps stable and there is no longer change. Convergence time is the time cost for the algorithm to find the optimal value. As shown in Table 5, four algorithms show different performance on different functions. GA algorithm and PSO algorithm have the similar performance and can only find the global value of Rastrigin function. When handling other functions, GA and PSO algorithms are stuck in local extremums with early iterations and cannot jump out local extremums. Due to the structural property, the convergence time of PSO is a little bit of higher than that of GA algorithms. ALO algorithm finds the global optimal value of all the test functions but Branin function with only a little deviation. IALO algorithm explores the solution spaces thoroughly and finds the global optimum values of all of the test function. Comparing with ALO algorithm, IALO algorithm can find the global value by less iteration and less computation time and show superior search ability and fast convergence property. In comparison, IALO algorithm has the best optimization performance among all the algorithms.

(2) Parameter Identification. The model parameters are optimized by the proposed IALO algorithm. The initial population is . The fitness of the individual is defined as follows:where N is the time length of the experiment data, denotes the measured terminal voltage at hth time, and is the estimated terminal voltage at hth time calculated by the corresponding individual (battery model parameters). The maximum iteration of IALO is set to 500, and iteration stops when the fitness of elite individual does not change in 10 continuous iterations.
The results of parameter identification are shown in Figure 6. As can be seen from the figure, the fitness convergence curves show that the ALO algorithm has early prematurity phenomenon and may fall in local optimal. On the contrary, the proposed improved ALO algorithm can further explore the solution space and find better solution. Comparing the conventions used battery parameters identification methods, recursive least square (RLS) algorithm, ALO and IALO can find the optimal solution in 50 iterations, while RLS needs more than 500 iterations to reach the parameter values. Also, there are fluctuations in the results of RLS, which is adverse to the optimizing process. The results of ALO and IALO are more stable and beneficial for optimization. The final identified parameters are shown in Table 6. The estimated terminal voltage is shown in Figure 7. The root mean square error (RMSE) of voltage estimation error by IALO algorithm is the smallest and equal to 0.0079 which is much smaller than the normal working range of voltage (2.5V4.2V). These verify that the parameters obtained by IALO optimization algorithm can fit the measurement data quite well.

(a)
(b)
(c)
(d)
(a)
(b)
(3) SOC Estimation. After the battery model parameters are obtained, the proposed SVDUKF is utilized for SOC estimation. However, the process noise covariance matrix Q and the measurement noise covariance matrix R are always unknown and usually tuned by “trialanderror” strategy manually by choosing random or experiential values, which may be very difficult, timeconsuming, and unreliable. To solving this problem, the proposed IALO algorithm is applied to search the optimal value for Q and R.
In UKF observer, Q and R matrixes are defined as diagonal and positive. In this paper, Q is matrix and R is matrix. Then the initial population is set as
The fitness of the individual is defined as follows:where is the weight of voltage and SOC estimation error, respectively, and and in this paper. N is the time length of the experiment data, denotes the measured terminal voltage at hth time, denotes the estimated terminal voltage at hth time, and denotes the true SOC which is obtained by calculating the integral of current over time via the high accuracy current sensors of the battery test system and the precision control initial SOC. denotes the estimated SOC based on the corresponding individual (diagonal element of Q and R matrixes). The maximum iteration of IALO is set to 500, and iteration stops when the fitness of elite individual does not change in 10 continuous iterations. Then the first two elements of the optimal results are assigned to the covariance matrixes Q and the third element of the optimal results is assigned to the covariance matrixes R. Based on the optimal Q and R and the former obtained battery model, the SVDUKF is conducted for the SOC estimation via DST. For comparison, the UKF and EKF with arbitrarily tuning parameters are also applied for SOC estimation. The estimation results are demonstrated in Figure 8 and Table 7. From the results, the proposed IALOSVDUKF algorithm can obtain precise SOC estimation and the RMSE is nearly 10 times smaller than that of UKF and 30 times smaller than that of EKF, and the maximum error rate is smaller than 1%, which can certify that the IALO algorithm can explore and exploit the solution space thoroughly and find the better solution quickly and effectively and is laborsaving. Also, the terminal voltage estimation of three methods is comparable with each other and the maximum of the error is 2 orders of magnitude smaller than the normal range of battery voltage (2.5V4.2V) and the maximum error rate is less than 0.5%.

To further validate the performance of IALO algorithm, the battery model parameters and the noise covariance matrixes are adopted to estimate the battery SOC under different dynamic conditions. The estimation results via FUDDS, US06, and BJDST are shown in Figures 9–11, similarly, comparing with UKF and EKF. The RMSEs of SOC and voltage estimation are shown in Table 8.

It is clearly shown that, under FUDS cycle condition, with the parameters and noise covariance matrixes identified by IALO algorithm, SVDUKF can obtain more than 30 times smaller SOC estimation RMSE comparing with UKF and EKF which are using arbitrary noise covariance matrixes. Under other more complicated and fluctuated cycle condition, like US06 and BJDST, IALOSVDUKF algorithm can still achieve nearly 4 times smaller SOC estimation RMSE. As for the terminal voltage estimation, the developed algorithm can acquire small RMSE which is less than that of UKF algorithm. EKF algorithm obtains the smallest voltage estimation RMSE but the largest SOC estimation RMSE. In the EV application, BMS is more concerned about SOC estimation. Therefore, the proposed algorithm can obtain the optimal value of the covariance matrix Q and R and then estimate SOC with highdegree precision, which is more preferable for EV and BMS application.
6. Conclusion
As a critical part of EV BMS algorithms, accurate SOC estimation and battery model can guarantee a precise control of battery inner state, lengthen the battery usage life, and reduce mileage anxiety. The widely used SOC estimation methods including EKF and UKF algorithms have the drawbacks of that the prior knowledge of noise covariance matrixes, which are often selected manually by “trialanderror” method and are arbitrary and timeconsuming.
This paper developed an improved ant lion optimization (IALO) algorithm to identify the battery model parameters and noise covariance matrixes. The IALO algorithm is based on a recently developed heuristics search method, namely, ant lion optimization (ALO) algorithm. To overcome the drawbacks of ALO, including early convergence and local optimum trapping, the chaotic mapping is utilized for generating the initial position of population and the random walk sequences to ensure the population homogeneity and ergodicity. Another improvement is the individual mutant which transfers each element of the elite ant individual according chaotic sequences generated by intertwining logistic map, to jump out local extremum and thoroughly explore the solution space.
To validate the performance of developed IALO algorithm for battery model parameters identification and SOC estimation, systematic experiments are conducted on NMC cells over different dynamic conditions. The results verify that the developed algorithm can obtain optimal battery model parameters in only 50 iterations. The accuracy of SOC estimation is less than 1% maximum error rate. The terminal voltage estimation is also desirable with less than 0.5% maximum error rate. Therefore, the proposed IALO algorithm has superior performance and can achieve fast and reliable results. Future work will be concentrated on online application of the IALO for SOC estimation.
Nomenclature
LiNMC:  Lithium nickelmanganesecobalt oxide 
EVs:  Electric vehicles 
SOC:  State of charge 
ALO:  Ant lion optimizer 
IALO:  Improved ant lion optimizer 
UKF:  Unscented Kalman filter 
EKF:  Extended Kalman filtering 
SVD:  Singular value decomposition 
BMS:  Battery Management System 
ECMs:  Equivalent circuit models 
RLS:  Recursive least square 
PSO:  Particle swarm optimization 
SA:  Simulated annealing 
GA:  Genetic algorithm 
DG:  Distributed Generation 
UT:  Unscented transformation 
DST:  Dynamic stress test 
FUDS:  Federal Urban Driving Schedule 
BJDST:  Beijing Dynamic Stress Test 
US06:  US6 Highway Driving Schedule 
:  Open circuit voltage 
:  Ohmic resistance 
:  Polarization resistance 
:  Polarization capacitance 
:  Load current 
:  Terminal voltage 
:  Polarization voltage over the polarization resistance 
:  Multinomial coefficients 
:  Sample interval 
:  Time constant 
:  The SOC at th sampling moment 
:  Rated capacity 
:  The coulomb efficiency 
:  Process noise signal 
:  Measurement noise signal 
:  Process noise covariance matrix 
:  Measurement noise covariance matrix 
:  Sigma points 
:  Random walk over the chosen ant lion 
:  Random walk over the elite ant lion 
:  Updated position of ith ant 
:  th chaotic value 
:  Weight coefficients of mean 
:  Weight coefficients of covariance 
:  Covariance matrix 
:  Scale adjustment parameter 
:  Scale adjustment parameter 
:  Scale adjustment parameter 
:  Output sigma points 
:  Kalman gain 
:  Random walk sequences 
:  Biggest number of iterations 
r(t):  A stochastic function 
N:  The number of ants and ant lions 
D:  The dimension of position 
:  The positions matrix of ants 
:  The value of th dimension of position value of the th ant 
:  Fitness function of ant 
:  Fitness value matrix of ants 
:  Position matrix of antlion 
:  The value of th dimension of position of the th antlion 
:  Fitness function of antlion 
:  Fitness value matrix of ant lions 
:  Normalized random walk displacement 
:  Lowest random walk displacement of ith ant 
:  Highest random walk displacement of ith ant 
:  The lowest of random walk displacement of th ant at current iteration 
:  The highest of random walk displacement of th ant at current iteration 
:  The lowest of random walk displacement of all ant at current iteration 
:  The highest of random walk displacement of all ant at current iteration 
:  Ration of the shrinking of random region 
:  Index varying with the current iteration 
:  Elite individual 
:  New individual 
:  New individual. 
Data Availability
The data used to support the findings of this study have been uploaded to website (https://web.calce.umd.edu/batteries/data.htm#type3) which can be downloaded freely.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was funded by the National Key R&D Program of China, 2017YFC0803904, China Postdoctoral Science Foundation, 2018T111006 and 2017M613034, Postdoctoral Science Foundation of Shaanxi Province, 2017BSHEDZZ36, Shaanxi Province Industrial Innovation Chain Project, 2018ZDCXLGY050301, Shaanxi Provincial Key Research and Development Plan Project, 2018ZDXMGY082, and Shaanxi Innovative Talents Promotion Plan Project, 2018KJXX005.
References
 X. Zhao, Q. Yu, J. Ma, Y. Wu, M. Yu, and Y. Ye, “Development of a representative EV urban driving cycle based on a kmeans and SVM hybrid clustering algorithm,” Journal of Advanced Transportation, vol. 2018, Article ID 1890753, 18 pages, 2018. View at: Publisher Site  Google Scholar
 H. Chen, B. Xie, J. Ma, and Y. Chen, “NOx emission of biodiesel compared to diesel: Higher or lower?” Applied Thermal Engineering, vol. 137, pp. 584–593, 2018. View at: Publisher Site  Google Scholar
 X. Zhao, J. Ma, S. Wang, Y. Ye, Y. Wu, and M. Yu, “Developing an electric vehicle urban driving cycle to study differences in energy consumption,” Environmental Science and Pollution Research, pp. 1–15, 2018. View at: Google Scholar
 Y. Zou, T. Liu, D. Liu, and F. Sun, “Reinforcement learningbased realtime energy management for a hybrid tracked vehicle,” Applied Energy, vol. 171, pp. 372–382, 2016. View at: Publisher Site  Google Scholar
 F. Zheng, Y. Xing, J. Jiang, B. Sun, J. Kim, and M. Pecht, “Influence of different open circuit voltage tests on state of charge online estimation for lithiumion batteries,” Applied Energy, vol. 183, pp. 513–525, 2016. View at: Publisher Site  Google Scholar
 X. Zhao, Q. Yu, M. Yu, and Z. Tang, “Research on an equal power allocation electronic differential system for electric vehicle with dualwheeledmotor front drive based on a wavelet controller,” Advances in Mechanical Engineering, vol. 10, no. 2, pp. 1–24, 2018. View at: Google Scholar
 R. Xiong, F. Sun, X. Gong, and C. Gao, “A datadriven based adaptive state of charge estimator of lithiumion polymer battery used in electric vehicles,” Applied Energy, vol. 113, pp. 1421–1433, 2014. View at: Publisher Site  Google Scholar
 Y. Xing, E. W. M. Ma, K. L. Tsui, and M. Pecht, “Battery management systems in electric and hybrid vehicles,” Energies, vol. 4, no. 11, pp. 1840–1857, 2011. View at: Publisher Site  Google Scholar
 Y. Xing, W. He, M. Pecht, and K. L. Tsui, “State of charge estimation of lithiumion batteries using the opencircuit voltage at various ambient temperatures,” Applied Energy, vol. 113, pp. 106–115, 2014. View at: Publisher Site  Google Scholar
 S. Piller, M. Perrin, and A. Jossen, “Methods for stateofcharge determination and their applications,” Journal of Power Sources, vol. 96, no. 1, pp. 113–120, 2001. View at: Publisher Site  Google Scholar
 J. H. Aylor, A. Thieme, and B. W. Johnson, “A battery stateofcharge indicator for electric wheelchairs,” IEEE Transactions on Industrial Electronics, vol. 39, no. 5, pp. 398–409, 1992. View at: Publisher Site  Google Scholar
 T. Liu, D. Chen, and C. Fang, “Design and implementation of a battery charger with a stateofcharge estimator,” International Journal of Electronics, vol. 87, no. 2, pp. 211–226, 2000. View at: Publisher Site  Google Scholar
 R. Xiong, H. He, F. Sun, X. Liu, and Z. Liu, “Modelbased state of charge and peak power capability joint estimation of lithiumion battery in plugin hybrid electric vehicles,” Journal of Power Sources, vol. 229, pp. 159–169, 2013. View at: Publisher Site  Google Scholar
 X. Zhao, S. Wang, J. Ma, and Q. Yu, “Identification of driver's braking intention based on a hybrid model of GHMM and GGAPRBFNN,” Neural Computing & Applications, pp. 1–24, 2018. View at: Google Scholar
 B. Cheng, Z. Bai, and B. Cao, “State of charge estimation based on evolutionary neural network,” Energy Conversion and Management, vol. 49, no. 10, pp. 2788–2794, 2008. View at: Publisher Site  Google Scholar
 X. Hu and F. Sun, “Fuzzy clustering based multimodel support vector regression state of charge estimator for lithiumion battery of electric vehicle,” in Proceedings of the 2009 International Conference on Intelligent HumanMachine Systems and Cybernetics, IHMSC 2009, pp. 392–396, China, August 2009. View at: Google Scholar
 P. Singh and D. Reisner, “Fuzzy logicbased stateofhealth determination of lead acid batteries,” in Proceedings of the 24th International Telecomunications Energy Conference (INTELEC 2002), pp. 583–590, IEEE, Canada, October 2002. View at: Google Scholar
 Y. Shen, “Adaptive online stateofcharge determination based on neurocontroller and neural network,” Energy Conversion and Management, vol. 51, no. 5, pp. 1093–1098, 2010. View at: Publisher Site  Google Scholar
 S. Lee and J. Kim, “Discrete wavelet transformbased denoising technique for advanced stateofcharge estimator of a lithiumion battery in electric vehicles,” Energy, vol. 83, pp. 462–473, 2015. View at: Publisher Site  Google Scholar
 Y. Li, C. Wang, and J. Gong, “A wavelet transformadaptive unscented Kalman filter approach for state of charge estimation of LiFePo4 battery,” International Journal of Energy Research, vol. 42, no. 2, pp. 587–600, 2018. View at: Publisher Site  Google Scholar
 H. Pan, Z. Lü, W. Lin, J. Li, and L. Chen, “State of charge estimation of lithiumion batteries using a grey extended Kalman filter and a novel opencircuit voltage model,” Energy, vol. 138, pp. 764–775, 2017. View at: Publisher Site  Google Scholar
 C. Weng, J. Sun, and H. Peng, “A unified opencircuitvoltage model of lithiumion batteries for stateofcharge estimation and stateofhealth monitoring,” Journal of Power Sources, vol. 258, pp. 228–237, 2014. View at: Publisher Site  Google Scholar
 B. Xia, C. Chen, Y. Tian, M. Wang, W. Sun, and Z. Xu, “State of charge estimation of lithiumion batteries based on an improved parameter identification method,” Energy, vol. 90, pp. 1426–1434, 2015. View at: Publisher Site  Google Scholar
 J. Li, L. Wang, C. Lyu, and M. Pecht, “State of charge estimation based on a simplified electrochemical model for a single LiCoO2 battery and battery pack,” Energy, vol. 133, pp. 572–583, 2017. View at: Publisher Site  Google Scholar
 J. Li, L. Wang, C. Lyu, H. Wang, and X. Liu, “New method for parameter estimation of an electrochemicalthermal coupling model for LiCoO 2 battery,” Journal of Power Sources, vol. 307, pp. 220–230, 2016. View at: Publisher Site  Google Scholar
 J. Li, Q. Lai, L. Wang, C. Lyu, and H. Wang, “A method for SOC estimation based on simplified mechanistic model for LiFePO4 battery,” Energy, vol. 114, pp. 1266–1276, 2016. View at: Publisher Site  Google Scholar
 Z. Chen, S. Qiu, M. A. Masrur, and Y. L. Murphey, “Battery state of charge estimation based on a combined model of extended kalman filter and neural networks,” in Proceedings of the International Joint Conference on Neural Network, IJCNN 2011, pp. 2156–2163, USA, August 2011. View at: Publisher Site  Google Scholar
 J. Han, D. Kim, and M. Sunwoo, “Stateofcharge estimation of leadacid batteries using an adaptive extended Kalman filter,” Journal of Power Sources, vol. 188, no. 2, pp. 606–612, 2009. View at: Publisher Site  Google Scholar
 R. Xiong, H. He, F. Sun, and K. Zhao, “Evaluation on state of charge estimation of batteries with adaptive extended kalman filter by experiment approach,” IEEE Transactions on Vehicular Technology, vol. 62, no. 1, pp. 108–117, 2013. View at: Publisher Site  Google Scholar
 Y. Li, C. Wang, and J. Gong, “A multimodel probability SOC fusion estimation approach using an improved adaptive unscented Kalman filter technique,” Energy, vol. 141, pp. 1402–1415, 2017. View at: Publisher Site  Google Scholar
 M. Partovibakhsh and G. Liu, “An adaptive unscented kalman filtering approach for online estimation of model parameters and stateofcharge of lithiumion batteries for autonomous mobile robots,” IEEE Transactions on Control Systems Technology, vol. 23, no. 1, pp. 357–363, 2015. View at: Publisher Site  Google Scholar
 W. He, N. Williard, C. Chen, and M. Pecht, “State of charge estimation for Liion batteries using neural network modeling and unscented Kalman filterbased error cancellation,” International Journal of Electrical Power & Energy Systems, vol. 62, pp. 783–791, 2014. View at: Publisher Site  Google Scholar
 S. Mirjalili and A. Lewis, “The whale optimization algorithm,” Advances in Engineering Software, vol. 95, pp. 51–67, 2016. View at: Publisher Site  Google Scholar
 M. A. Hossain, H. R. Pota, S. Squartini, and A. F. Abdou, “Modified PSO algorithm for realtime energy management in gridconnected microgrids,” Renew Energy, vol. 136, pp. 746–757, 2019. View at: Publisher Site  Google Scholar
 W. Zhang, A. Maleki, M. A. Rosen, and J. Liu, “Optimization with a simulated annealing algorithm of a hybrid system for renewable energy including battery and hydrogen storage,” Energy, vol. 163, pp. 191–207, 2018. View at: Publisher Site  Google Scholar
 Z. Chen, S. Zhou, and J. Luo, “A robust ant colony optimization for continuous functions,” Expert Systems with Applications, vol. 81, pp. 309–320, 2017. View at: Publisher Site  Google Scholar
 N. Ghorbani, A. Kasaeian, A. Toopshekan, L. Bahrami, and A. Maghami, “Optimizing a hybrid windPVbattery system using GAPSO and MOPSO for reducing cost and increasing reliability,” Energy, vol. 154, pp. 581–591, 2018. View at: Publisher Site  Google Scholar
 S. M. AbdElazim and E. S. Ali, “Load frequency controller design of a twoarea system composing of PV grid and thermal generator via firefly algorithm,” Neural Computing and Applications, vol. 30, no. 2, pp. 607–616, 2018. View at: Publisher Site  Google Scholar
 A. S. Oshaba, E. S. Ali, and S. M. Abd Elazim, “PI controller design for MPPT of photovoltaic system supplying SRM via BAT search algorithm,” Neural Computing and Applications, vol. 28, no. 4, pp. 651–667, 2017. View at: Publisher Site  Google Scholar
 S. Mirjalili, “The ant lion optimizer,” Advances in Engineering Software, vol. 83, pp. 80–98, 2015. View at: Publisher Site  Google Scholar
 S. Mirjalili, P. Jangir, and S. Saremi, “Multiobjective ant lion optimizer: a multiobjective optimization algorithm for solving engineering problems,” Applied Intelligence, vol. 46, no. 1, pp. 79–95, 2017. View at: Publisher Site  Google Scholar
 E. S. Ali, S. M. Abd Elazim, and A. Y. Abdelaziz, “Ant Lion Optimization Algorithm for optimal location and sizing of renewable distributed generations,” Journal of Renewable Energy, vol. 101, pp. 1311–1324, 2017. View at: Publisher Site  Google Scholar
 T. Tian, C. Liu, Q. Guo, Y. Yuan, W. Li, and Q. Yan, “An improved ant lion optimization algorithm and its application in hydraulic turbine governing system parameter identification,” Energies, vol. 11, no. 1, p. 95, 2018. View at: Google Scholar
 Z. Wu, D. Yu, and X. Kang, “Parameter identification of photovoltaic cell model based on improved ant lion optimizer,” Energy Conversion and Management, vol. 151, pp. 107–115, 2017. View at: Publisher Site  Google Scholar
 M. A. Hannan, M. S. H. Lipu, A. Hussain, and A. Mohamed, “A review of lithiumion battery state of charge estimation and management system in electric vehicle applications: Challenges and recommendations,” Renewable & Sustainable Energy Reviews, vol. 78, pp. 834–854, 2017. View at: Publisher Site  Google Scholar
 X. S. Hu, S. B. Li, and H. Peng, “A comparative study of equivalent circuit models for Liion batteries,” Journal of Power Sources, vol. 198, pp. 359–367, 2012. View at: Publisher Site  Google Scholar
 E. Ott, Chaos in Dynamical Systems, Cambridge University Press, Cambridge, UK, 2nd edition, 2002. View at: Publisher Site  MathSciNet
 A. Kaveh, R. Sheikholeslami, S. Talatahari, and M. KeshvariIlkhichi, “Chaotic swarming of particles: A new method for size optimization of truss structures,” Advances in Engineering Software, vol. 67, pp. 136–147, 2014. View at: Publisher Site  Google Scholar
 M. Ragulskis, Z. Navickas, R. Palivonaite, and M. Landauskas, “Algebraic approach for the exploration of the onset of chaos in discrete nonlinear dynamical systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 11, pp. 4304–4315, 2012. View at: Publisher Site  Google Scholar
 G. Ye and X. Huang, “An efficient symmetric image encryption algorithm based on an intertwining logistic map,” Neurocomputing, vol. 251, pp. 45–53, 2017. View at: Publisher Site  Google Scholar
 X. Wang and D. Luan, “A novel image encryption algorithm using chaos and reversible cellular automata,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 11, pp. 3075–3085, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 T. Q. Duong, “USABC and PNGV test procedures,” Journal of Power Sources, vol. 89, no. 2, pp. 244–248, 2000. View at: Publisher Site  Google Scholar
 W. He, N. Williard, C. C. Chen, and M. Pecht, “Robust and adaptive estimation of state of charge for lithiumion batteries,” International Journal of Power Electronics, vol. 62, pp. 783–791, 2014. View at: Google Scholar
 Z. Ma, J. Jiang, W. Shi, W. Zhang, and C. C. Mi, “Investigation of path dependence in commercial lithiumion cells for pure electric bus applications: aging mechanism identification,” Journal of Power Sources, vol. 274, pp. 29–40, 2015. View at: Publisher Site  Google Scholar
 M. M. Maricq, D. H. Podsiadlik, and R. E. Chase, “Gasoline vehicle particle size distributions: Comparison of steady state, FTP, and US06 measurements,” Environmental Science & Technology, vol. 33, no. 12, pp. 2007–2015, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Kai Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.