International Scholarly Research Notices

International Scholarly Research Notices / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 120351 | 17 pages |

Estimation Strategies for the Condition Monitoring of a Battery System in a Hybrid Electric Vehicle

Academic Editor: F. Piazza
Received04 Jan 2011
Accepted06 Feb 2011
Published13 Apr 2011


This paper discusses the application of condition monitoring to a battery system used in a hybrid electric vehicle (HEV). Battery condition management systems (BCMSs) are employed to ensure the safe, efficient, and reliable operation of a battery, ultimately to guarantee the availability of electric power. This is critical for the case of the HEV to ensure greater overall energy efficiency and the availability of reliable electrical supply. This paper considers the use of state and parameter estimation techniques for the condition monitoring of batteries. A comparative study is presented in which the Kalman and the extended Kalman filters (KF/EKF), the particle filter (PF), the quadrature Kalman filter (QKF), and the smooth variable structure filter (SVSF) are used for battery condition monitoring. These comparisons are made based on estimation error, robustness, sensitivity to noise, and computational time.

1. Introduction

Condition monitoring is an essential process for fault detection and diagnosis. It involves monitoring system states or parameters over an operational period, where abnormal values or significant changes would indicate a fault. Quite often direct measurements of the states are not available due to limitations in design or cost. In these cases, state and parameter estimation techniques can be used for information extraction. Condition monitoring of systems allows proper maintenance to be scheduled, which helps reduce unscheduled downtime of manufacturing equipment, as well as the cost to repair damaged systems [1, 2]. An important area for condition monitoring is energy management for hybrid electric (HEVs) and battery electric vehicles (BEVs).

In general, HEVs have two power sources: a gasoline engine and an electric motor. In full hybrid vehicles, the engine and the motor can operate separately or simultaneously. The motor is used mainly during acceleration, startup, reverse mode, and in regenerative braking. A traction battery pack is used to provide power to the motor. It is recharged by a generator or during regenerative braking. The performance of an HEV is largely dependent on a balance between the gasoline engine and the electric motor, optimized with respect to fuel consumption based on vehicle conditions [3]. Many different types of control methods have been applied to balance the power and energy requirements of HEVs, including fuzzy logic [4ā€“6], genetic algorithms [7], dynamic programming [8, 9], Pareto optimization [10], and intelligent mechanism designs [11]. These control strategies rely heavily on the availability of battery power to balance the operation of the gasoline engine versus the electric motor. The available battery power may be obtained from the state-of-charge (SOC) information [3, 12]. Further to the SOC, the battery state-of-health (SOH) is required in order to help determine whether a battery would fail subject to a certain load [13]. The SOC cannot be measured directly with electric signals, and as such it often, needs to be estimated through current and voltage relationships [3]. Poor estimation or control of the SOC may lead to improper charging conditions and can degrade the efficiency and reliability of the batteries [13]. Hence, proper condition monitoring of batteries plays a pivotal role in the optimization of HEV performance, as well as extending the lifetime and increasing the reliability of the batteries [14].

State and parameter estimation techniques are an integral part of condition monitoring and are used when direct measurements of the states are not available. One of the most commonly studied methods for estimation is the Kalman filter (KF) for linear systems and its extended form (EKF) for nonlinear systems [15ā€“19]. Operating conditions such as battery SOC, power fade, capacity fade, resistance, and instantaneous available power have been estimated well using the EKF [17]. Other methods such as the sigma-point Kalman filtering (SPKF) and support vector machine (SVM) have also been used for condition monitoring and fault diagnosis [20ā€“22].

In critical applications, such as automotive that require added safety and reliability, the choice of the estimation method is very important and should be selected based on the linearity of the system, performance, robustness or sensitivity to noise and computational difficulty and time. In relation to battery condition monitoring, this paper presents a quantitative and qualitative comparison of the following methods: the Kalman and extended Kalman filters (KF/EKF), the particle filter (PF), the quadrature Kalman filter (QKF), and the smooth variable structure filter (SVSF).

2. State and Parameter Estimation Techniques

State and parameter estimation is essential for sensing and information processing in model-based condition monitoring. Estimation theory involves information extraction by tracking changes in physical parameters or operational states of the system. This paper studies four popular strategies: the Kalman and the extended Kalman filters (KF/EKF), the particle filter (PF), the quadrature Kalman filter (QKF), and the smooth variable structure filter (SVSF).

2.1. Kalman and Extended Kalman Filters

Even after 50 years, the Kalman filter (KF) remains the most studied and one of the most popular tools used in state estimation [23ā€“26]. It may be applied to linear dynamic systems in the presence of Gaussian white noise, and provides an elegant and statistically optimal solution by minimizing the mean-squared error. The impact that the KF has had on estimation and control problems is considered by some scientists and engineers to be one of the greatest achievement in engineering and signal processing [26]. It is a method that utilizes measurements linearly related to the states, and error covariance matrices, to generate a gain referred to as the Kalman gain. This gain is applied to the a priori state estimate, thus creating an a posteriori estimate. The estimation process continues in a predictor-corrector fashion while maintaining a statistically minimal state error covariance matrix for linear systems.

The following two equations describe the system dynamic model and the measurement model used in general for state estimation. Refer to Appendix A for a description of the nomenclature š‘„š‘˜+1=š“š‘˜š‘„š‘˜+šµš‘˜š‘¢š‘˜+š‘¤š‘˜,š‘§š‘˜+1=š¶š‘˜+1š‘„š‘˜+1+š‘£š‘˜+1.(1)

The next five equations form the KF algorithm and are used in an iterative fashion. Equation (2) extrapolates the a priori state estimate, and (3) is the corresponding a priori error covariance. The Kalman gain may be calculated by (4), and is used to update the state estimate and error covariance, described by (5) and (6), respectively. Ģ‚š‘„š‘˜+1|š‘˜=š“š‘˜Ģ‚š‘„š‘˜|š‘˜+šµš‘˜š‘¢š‘˜,š‘ƒ(2)š‘˜+1|š‘˜=š“š‘˜š‘ƒš‘˜|š‘˜š“š‘‡š‘˜+š‘„š‘˜š¾,(3)š‘˜=š‘ƒš‘˜+1|š‘˜š¶š‘‡š‘˜ī€ŗš¶š‘˜š‘ƒš‘˜+1|š‘˜š¶š‘‡š‘˜+š‘…š‘˜ī€»āˆ’1,(4)Ģ‚š‘„š‘˜+1|š‘˜+1=Ģ‚š‘„š‘˜+1|š‘˜+š¾š‘˜ī€ŗš‘§š‘˜āˆ’š¶š‘˜š‘„š‘˜+1|š‘˜ī€»š‘ƒ,(5)š‘˜+1|š‘˜+1=ī€ŗš¼āˆ’š¾š‘˜š¶š‘˜ī€»š‘ƒš‘˜+1|š‘˜.(6)

The effects due to model uncertainties can have a large impact on the stability and performance of the KF [26, 27]. For nonlinear systems, the EKF may be used. It is conceptually similar to the KF process. The nonlinear system and measurement matrices are linearized according to their corresponding Jacobian, which is a first-order partial derivative. This linearization introduces uncertainties in the estimation process; such that overlooked nonlinearities in the system may cause the EKF to become unstable [26].

2.2. Particle Filter

The particle filter (PF) has many forms: Monte Carlo filters, interacting particle approximations [28], bootstrap filters [29], condensation algorithm [30], and survival of the fittest [31], to name a few. Compared to the KF, it is newer, being introduced in 1993. Since then, the PF has become a very popular method for solving nonlinear estimation problems, ranging from predicting chemical processes to target tracking. The PF takes the Bayesian approach to dynamic state estimation, in which one attempts to accurately represent the probability distribution function (PDF) of the values of interest [32]. The PDF contains all of the pertinent statistical information and may be considered as holding the solution to the estimation problem [32]. Essentially, the distribution holds a probability of values for the state being observed. The stronger or tighter the prediction PDF, the more accurate the state estimate.

The PF obtains its name from the use of weighted particles or ā€œpoint massesā€ that are distributed throughout the PDF to form an approximation. These particles are used in an iterative process to obtain new particles and associated importance weights, with the goal of creating a more accurate approximation of the PDF. In general, as the number of implemented particles becomes large, the PDF becomes more accurate [32]. An important step in the PF is that of resampling, which eliminates particles with low weights and multiplies those with high weights [32]. This helps to avoid the degeneracy problem with the PF, which refers to only one particle having a significant importance weight after a large number of recursions. Furthermore, it also increases the accuracy of the PDF approximation by replicating particles with high weights. The sequential importance resampling (SIR) algorithm is a very popular form of the PF and may be summarized by (7) to (10). The first equation draws samples or particles from the proposal distribution. š‘„š‘˜(š‘›)ī‚€š‘„āˆ¼šœ‹š‘˜āˆ£š‘„(š‘›)š‘˜āˆ’1,š‘¦š‘˜ī‚.(7) The next equation updates the importance weights up to a normalizing constant īšœ”š‘˜(š‘›)=šœ”(š‘›)š‘˜āˆ’1š‘ī‚€š‘¦š‘˜āˆ£š‘„š‘˜(š‘›)ī‚š‘ī‚€š‘„š‘˜(š‘›)āˆ£š‘„(š‘›)š‘˜āˆ’1ī‚šœ‹ī‚€š‘„š‘˜āˆ£š‘„(š‘›)š‘˜āˆ’1,š‘¦š‘˜ī‚.(8) Next, the normalized weights are calculated for each particle šœ”š‘˜(š‘›)=īšœ”š‘˜(š‘›)āˆ‘š‘›š‘–=1īšœ”š‘˜(š‘–).(9)

Finally, a constant known as the effective number of particles is calculated as shown in (10). Resampling is performed if the effective number of particles is lower than some design threshold īš‘eļ¬€=1āˆ‘š‘›š‘–=1ī‚€šœ”š‘˜(š‘–)ī‚2.(10)

2.3. Quadrature Kalman Filter

Similar to the PF, the quadrature Kalman filter (QKF) is a type of Bayesian filter that is able to model dynamic processes which are nonlinear and subject to non-Gaussian noise. In 2007, it was proposed that a set of Gauss-Hermite quadrature points could be used to parameterize the PDF [33, 34]. When compared with the EKF, it was found that the QKF method provides a more accurate least-squares solution [33]. Figure 1 shows the PDF of a nonlinearly transformed Gaussian random variable. The mean of the EKF appears to be biased, and the covariance is obviously far from the true covariance. However, the QKF mean and covariance match the true values quite well, with the 5 point QKF working the best. The main drawback to this method is the fact that the number of terms used in the Gaussian sum grows exponentially, which means that more memory will be used over a longer period of time [33]. The process of this filter is similar to the PF, in the sense that it may be solved recursively using two stages (time and measurement update). Please refer to Appendix B for the full QKF algorithm.

2.4. Smooth Variable Structure Filter

In 2002, the variable structure filter (VSF) was introduced as a new predictor-corrector method used for state and parameter estimation [27, 35]. It is a type of sliding mode estimator, where gain switching is used to ensure that the estimates converge to true state values. An internal model of the system, either linear or nonlinear, is used to predict an a priori state estimate. A corrective term is then applied to calculate the a posteriori state estimate, and the estimation process is repeated iteratively. The SVSF was later derived from the VSF and uses a simpler and less complex gain calculation [36]. In its present form, the SVSF is stable and robust to modeling uncertainties and noise, given an upperbound on uncertainties [36]. The basic concept of the SVSF is shown in Figure 2. Assume that the solid line in Figure 2 is a trajectory of some state (amplitude versus time). An initial value is selected for the state estimate. The estimated state is pushed towards the true value until it reaches a subspace around the actual state trajectory, referred to as the existence subspace. Once the value enters the existence subspace, the estimated state is forced to remain within it and into switching along the system state trajectory [36].

The SVSF method is model based and applies to smooth nonlinear dynamic equations. The estimation process may be summarized by (11) to (16) and is repeated iteratively. An a priori state estimate is calculated using an estimated model of the system. This value is then used to calculate an a priori estimate of the measurement defined by (12). A corrective term, referred to as the SVSF gain, is calculated as a function of the error in the predicted output, as well as a gain matrix and the smoothing boundary layer width. The corrective term calculated in (13) is then used in (14) to find the a posteriori state estimate Ģ‚š‘„š‘˜+1|š‘˜=īš¹ī€·Ģ‚š‘„š‘˜|š‘˜,š‘¢š‘˜ī€ø,(11)Ģ‚š‘§š‘˜+1|š‘˜=īš¶Ģ‚š‘„š‘˜+1|š‘˜š¾,(12)š‘˜+1=īš¶āˆ’1|||ī‚€|||š‘’š‘§š‘˜+1|š‘˜|||š“šµš‘†|||š‘’+š›¾š‘§š‘˜|š‘˜|||š“šµš‘†ī‚|||š“šµš‘†ī‚€š‘’āˆ˜satš‘§š‘˜+1|š‘˜ī‚,,ĪØ(13)Ģ‚š‘„š‘˜+1|š‘˜+1=Ģ‚š‘„š‘˜+1|š‘˜+š¾š‘˜+1š‘’,(14)š‘§š‘˜|š‘˜=š‘§š‘˜āˆ’Ģ‚š‘§š‘˜|š‘˜š‘’,(15)š‘§š‘˜+1|š‘˜=š‘§š‘˜+1āˆ’Ģ‚š‘§š‘˜+1|š‘˜.(16)

Two critical variables in this process are the a priori and a posteriori output error estimates, defined by (15) and (16), respectively [36]. Note that (15) is the output error estimate from the previous time step and is used only in the gain calculation.

3. Condition Monitoring of a Battery System in a Hybrid Electric Vehicle

A variety of batteries have been studied in literature, most notably lead-acid and lithium-ion batteries [13, 15, 19, 20, 37]. Lead-acid batteries are the oldest type of rechargeable batteries, and are most commonly found in motor vehicles. Lithium-ion batteries are also a form of rechargeable battery, which contain lithium in its positive electrode (cathode). These batteries are usually found in portable consumer electronics (i.e., laptops or notebooks) due to particularly high energy-to-weight ratios, slow self-discharge, and a lack of memory effect (i.e., where a battery loses its maximum energy capacity over time) [16]. In recent years, lithium-ion batteries have slowly entered the hybrid electric vehicle market, due to the fact that they offer better energy density compared to standard batteries [38].

The operation of batteries may be studied by using the advanced vehicle simulator (ADVISOR), which was written in MATLAB and Simulink by the US Department of Energy and the National Renewable Energy Laboratory [39ā€“41]. ADVISOR is used for the analysis of performance and fuel economy of three vehicle types: conventional, electric, and hybrid vehicles [39]. In 2001, the resistance-capacitance (RC) battery model was first implemented in ADVISOR [42]. The electrical model consists of three resistors (š‘…š‘’, š‘…š‘, and š‘…š‘”) and two capacitors (š¶š‘ and š¶š‘). The first capacitor (š¶š‘) represents the capability of the battery to chemically store a charge, and the second capacitor (š¶š‘) represents the surface effects of a cell [41]. The resistances and capacitances vary with changing SOC and temperature (š‘‡) [41]. ADVISOR offers two different datasets for the RC battery model: lithium-ion and nickel-metal hydride chemistries. For the purposes of this study, the lithium-ion chemistry was used in conjunction with the RC battery model. Figure 3 illustrates the equivalent circuit diagram of the RC model.

A standard model of a parallel hybrid electric vehicle referred to within ADVISOR as the Annex VII PHEV was used for this study. This model has been developed by the International Energy Agency (IEA), which is an international research community for the development and commercialization of hybrid and electric vehicles [43]. The model is based on data obtained from published sources and national (U.S.) research test data [39]. The battery system of the HEV represents the battery pack which stores energy on board the HEV. The system accepts a power request and returns the available power from the battery, as well as the SOC, voltage and current [39].

The nonlinear equations that describe the system may be derived from the RC battery model of Figure 3, as shown in Appendix C. For the purposes of condition monitoring of the battery, two voltages š‘‰š¶š‘ and š‘‰š¶š‘ as well as two capacitance parameters š¶š‘ and š¶š‘ need to be estimated. Should a fault exist in either of the battery capacitors, one would be able to determine this given the corresponding change in the parametric value. Further to the equations found in Appendix C, a discrete-time state space model of the Capacitor voltages may be defined as follows: ī‚øš‘‰š¶š‘š‘˜+1š‘‰š¶š‘š‘˜+1ī‚¹=āŽ”āŽ¢āŽ¢āŽ¢āŽ£āˆ’š‘‡š‘ š¶š‘(š‘…š‘’+š‘…š‘)š‘‡+1š‘ š¶š‘(š‘…š‘’+š‘…š‘)š‘‡š‘ š¶š‘(š‘…š‘’+š‘…š‘)āˆ’š‘‡š‘ š¶š‘(š‘…š‘’+š‘…š‘)āŽ¤āŽ„āŽ„āŽ„āŽ¦+1š‘˜ī‚øš‘‰š¶š‘š‘˜š‘‰š¶š‘š‘˜ī‚¹+āŽ”āŽ¢āŽ¢āŽ¢āŽ£š‘‡š‘ š‘…š‘š¶š‘(š‘…š‘’+š‘…š‘)š‘‡š‘ š‘…š‘’š¶š‘(š‘…š‘’+š‘…š‘)āŽ¤āŽ„āŽ„āŽ„āŽ¦š‘˜š¼š‘†š‘˜.(17)

For parameter estimation, (17) is used to formulate a state vector that would include parameters such that ī€ŗš‘‰š‘„=š¶š‘š‘‰š¶š‘š¶š‘š¶š‘ī€»š‘‡.(18)

By rearranging (17), and further to Appendix C, the general model used for state and parameter estimation is obtained as š‘„š‘˜+1ī€·š‘„=š‘“š‘˜,š‘¢š‘˜ī€ø+š‘¤š‘˜,š‘§š‘˜+1ī€·š‘„=ā„Žš‘˜+1,š‘¢š‘˜+1ī€ø+š‘£š‘˜+1.(19)

To implement the EKF method, the nonlinear system had to be linearized using a first-order Taylor series approximation. The linearized form of the system equation is provided in Appendix D.

4. Estimation Results

This paper presents a comparative performance for the application of the extended Kalman filter, the particle filter, the quadrature Kalman filter, and the smooth variable structure filter for condition monitoring of RC batteries. The study is conducted by simulation using the ADVISOR battery model. The model parameters are documented in [39].

Figures 4 and 5 illustrate the input current (š¼š‘ ) and the output voltages š‘‰š¶š‘ and š‘‰š¶š‘ used in the simulation. The parameter values š¶š‘ and š¶š‘ are made to vary in order to simulate fault conditions as shown in Figures 6 and 7. Parameters š¶š‘ and š¶š‘ as well as the states š‘‰š¶š‘ and š‘‰š¶š‘ are estimated using the four methods described in Section 2.

A comparison and discussion of the results is provided using two simulated cases: one with noise and one with noise as well as modeling errors. Note that the sampling time used in the simulation was 0.01 seconds.

4.1. Extended Kalman Filter Results

For linear dynamic systems in the presence of Gaussian white noise, the KF provides an elegant and statistically optimal solution by minimizing the mean-squared estimation error. The EKF is used for nonlinear problems. The following covariance matrices (error, process, and measurement, resp.) were used for the EKF and were obtained by trial-and-error: āŽ”āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ¦,āŽ”āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ¦.š‘ƒ=1303030301303030301303030301š‘„=š‘…=1000010000100001(20)

The EKF worked very well in the absence of modeling errors, as shown in Figure 8. Its performance degraded considerably when modeling errors were added, as illustrated in Figure 9.

As shown in Figure 10, the EKF was very sensitive to the selection of the initial conditions; however, it often recovered and performed very well. This is important to note because without selecting good starting points, the filter performed the worst. However, when good initial conditions were selected, the EKF yields the best results (based on RMSE and the assumption that there were no uncertainties in the model). When uncertainties (besides white measurement and process noise) were present in the filter model, in this example, the EKF became unstable (as shown in Figure 9) and failed to yield reasonable results. The simulation time for the EKF was one of the fastest, and compared to the other methods, it was computationally easy.

4.2. Particle Filter Results

For the PF method, a large number of particles (500) were used. Its application results are shown in Figures 11 and 12. The PF was able to estimate fairly well in the presence of noise. However, when modeling uncertainties were introduced, the PF had difficulty tracking the true voltages. As shown in Figure 12(b), after about 20 seconds, the PF was able to recover and provide a good estimate.

The PF provided satisfactory results. However, the technique had the highest RMSE (for the first case). That being said, when compared to the EKF, no Jacobian matrix had to be calculated to linearize the system matrix, as weighted particles were used instead. This may be an attractive feature if the system is too difficult to linearize. A large number of particles (which were required to accurately estimate the PDF), and the resampling feature of the PF, contributed to a slower computational time. When uncertainties were added to the filter model (as shown in Figure 12), the PF was able to recover slowly from modeling errors and after some time, accurately represent the first two states. There were significant errors in the third and fourth states; however this, as previously mentioned, may be attributed to measurements being only available for the first two states.

4.3. Quadrature Kalman Filter Results

For the QKF, the initial values of the covariance matrices were set to the identity matrix. The following results shown in Figures 13 and 14 were generated using the QKF based on the simulation setup. The QKF performed very well for both cases. However, like the previous methods, when uncertainties were introduced, it was unable to accurately track the two parameters associated with the battery (š¶š‘ and š¶š‘) as shown in Figure 14.

The QKF provided good results for the system. Like the PF, there was no need to calculate the Jacobian, since weighted quadrature points were used instead. One of the drawbacks of the QKF is the computational demand when a large number of states and quadrature points are used. This leads to a slower calculation time. For example, since four states were required, an array of 81 quadrature points was required (34). If a more accurate model was required, five quadrature points per state could be used. However, this would further increase the computational time. For this system, it was found that increasing the number of quadrature points had a negligible effect on the accuracy. Note that when uncertainties were added to the filter model, the QKF worked extremely well. However, the slow computation time is its main hindrance.

4.4. Smooth Variable Structure Filter Results

The constant diagonal matrix (š›¾) was set to 0.4, and the smoothing boundary layer thicknesses (šœ“) were, respectively, set to 1Ɨ10āˆ’3 and 1Ɨ10āˆ’2 for the two states and the two parameters, respectively. The results are shown in Figures 15 and 16, and were generated using the SVSF based on the simulation setup. The SVSF method is very robust to noise and modeling uncertainties, as demonstrated by the accurate estimation in both cases.

The SVSF yielded very good results for condition monitoring of the battery system. Similar to the previous two methods, it did not require calculation of the Jacobian matrix, and also required a time delay (signal extraction for the last two states). The SVSF simulation time was just as fast as the EKF, and the results were similar for the case with noise. When errors were added to the filter model, the SVSF worked the best (in terms of RMSE). Although there were still large estimation errors for the third and fourth states, they were finite and stable. It is important to note that the SVSF was robust and not sensitive to changes in the filter parameters and the initial conditions.

4.5. Summary of Results

The results of the simulations were compared based on estimation error, robustness, sensitivity to noise, and computational difficulty and time. As shown in Table 1, the EKF performed best interms of estimation error, in the case of only noise (referred to as Caseā€‰ā€‰(1)). In the presence of modeling uncertainties (Caseā€‰ā€‰(2)), the EKFā€™s performance severely degraded. The SVSF yielded the most accurate estimation for the case involving modeling uncertainties. Both the PF and QKF performed moderately well, but at a cost of computational time, as shown in the last row of Table 1.

States and parametersCaseā€‰ā€‰(1)Caseā€‰ā€‰(2)Caseā€‰ā€‰(1)Caseā€‰ā€‰(2)Caseā€‰ā€‰(1)Caseā€‰ā€‰(2)Caseā€‰ā€‰(1)Caseā€‰ā€‰(2)

RMSE š‘‰ š¶ š‘ 2.25Eā€‰āˆ’ā€‰095.61Eā€‰āˆ’ā€‰012.34Eā€‰āˆ’ā€‰049.50Eā€‰āˆ’ā€‰032.07Eā€‰āˆ’ā€‰054.19Eā€‰āˆ’ā€‰052.82Eā€‰āˆ’ā€‰071.24Eā€‰āˆ’ā€‰08
š‘‰ š¶ š‘ 1.10Eā€‰āˆ’ā€‰089.23Eā€‰āˆ’ā€‰022.95Eā€‰āˆ’ā€‰041.30Eā€‰āˆ’ā€‰032.07Eā€‰āˆ’ā€‰054.19Eā€‰āˆ’ā€‰056.58Eā€‰āˆ’ā€‰079.62Eā€‰āˆ’ā€‰09
š¶ š‘ 12.0259741.19119628.07119618.481195
š¶ š‘ 1.40892.191771.201770.73177

Simulation time (sec)0.380.443.273.348.358.420.390.39

The error in Table 1 refers tothe root mean-squared error (RMSE), calculated as follows: 1RMSE=š‘ī‚™ī“š‘š‘–=1ī€·š‘„š‘–āˆ’Ģ‚š‘„š‘–ī€ø2.(21)

Table 2 gives a quantitative comparison of the four filters. The robustness was based on observations made on varying the noise levels and considering its impact on model errors. The sensitivities were determined by tuning the filters various parameters. Computation (simulation) time was determined for each method and may be compared on a relative basis.


Sensitivity to UncertaintiesHighHighLowLow
Sensitivity to Filter
Computation TimeFastSlowVery SlowFast

5. Conclusions

This paper discussed the application of condition monitoring to an RC battery system, typically found within a hybrid electric vehicle. State and parameter estimation techniques are important as they are responsible for providing accurate estimates of the states when reliable measurements are unavailable, and hence ensure successful condition monitoring. A comparative study was presented in which the Kalman and the extended Kalman filters (KF/EKF), the particle filter (PF), the quadrature Kalman filter (QKF), and the smooth variable structure filter (SVSF) were applied for condition monitoring. These estimation methods were compared based on estimation error, robustness, sensitivity to noise, and computation time. Table 3 summarizes the results of the comparison. For the case without modeling errors, the EKF worked the best in terms of RMSE and computational speed. When modeling errors were present, which is common in physical applications, the SVSF was shown to work significantly better than the other methods in terms of stability and RMSE.

Only noiseNoise and model errors



A. List of Nomenclature

š“:System matrixšµ: Input matrixš¶:Output matrixš‘’:State estimation errorš‘“:Nonlinear system or process equationā„Ž:Nonlinear measurement or output equationš‘˜:Time step indexš¾:Gain value (KF, QKF, or SVSF)š‘š:Number of measurementsš‘›:Number of statesš‘eļ¬€:Effective number of particlesš‘ƒ: Error covariance matrixš‘„:System noise covariance matrixš‘…:Measurement noise covariance matrixSat:Saturation functionš‘”:Simulation timeš‘‡š‘ :Sampling timeš‘¢: Inputš‘£:Measurement noiseš‘¤:System noiseš‘¤š‘™:Quadrature point weightš‘„:System statesš‘¦:Observed measurement or outputš‘§:Measurement outputš›¾:Constant diagonal gain matrix with elements having values between 0 and 1šœ”:Particle weightĪØ:Smoothing boundary layeršœ‹:Probability distributionšœ:Sampling timeĢ‚ā€Œ:Denotes an estimated value~:Denotes an error valueā‹…:On top of a parameter denotes a time derivative.

Furthermore, note that subscript š‘˜+1āˆ£š‘˜ refers to an a priori time step and that the subscript š‘˜+1āˆ£š‘˜+1 refers an a posteriori time step. A superscript of š‘‡ denotes a matrix transpose.

B. Quadrature Kalman Filter Algorithm

The following is a summary of the QKF algorithm, directly as presented in [33, 34]. This process was used to obtain the presented results. There are two updates at each step: the time update and the measurement update.

Time Update Step
(1)Assuming at time š‘˜ the posterior density function š‘(š‘„š‘˜āˆ’1āˆ£š‘§š‘˜āˆ’1)=š‘(Ģ‚š‘„š‘˜āˆ’1|š‘˜āˆ’1,š‘ƒš‘˜āˆ’1|š‘˜āˆ’1) is known, then we may factorize as follows: š‘ƒš‘˜āˆ’1|š‘˜āˆ’1=ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1ī‚€ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1ī‚š‘‡.(B.1)(2)Evaluate the quadrature points {š‘‹š‘™,š‘˜āˆ’1|š‘˜āˆ’1}š‘šš‘™=1 as š‘‹š‘™,š‘˜āˆ’1|š‘˜āˆ’1=ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1šœ‰š‘™+Ģ‚š‘„š‘˜āˆ’1|š‘˜āˆ’1.(B.2)(3)Evaluate the propagated quadrature points {š‘‹āˆ—š‘™,š‘˜|š‘˜āˆ’1}š‘šš‘™=1 as:š‘‹āˆ—š‘™,š‘˜|š‘˜āˆ’1ī€·š‘‹=š‘“š‘™,š‘˜āˆ’1|š‘˜āˆ’1,š‘¢š‘˜āˆ’1ī€ø,š‘˜āˆ’1.(B.3)(4)Estimate the predicted state: Ģ‚š‘„š‘˜|š‘˜āˆ’1=ī“š‘šš‘™=1š‘¤š‘™š‘‹āˆ—š‘™,š‘˜|š‘˜āˆ’1.(B.4)(5) Estimate the predicted error covariance: š‘ƒš‘˜|š‘˜āˆ’1=ī“š‘šš‘™=1š‘¤š‘™š‘‹āˆ—š‘™,š‘˜|š‘˜āˆ’1š‘‹āˆ—š‘‡š‘™,š‘˜|š‘˜āˆ’1āˆ’Ģ‚š‘„š‘˜|š‘˜āˆ’1Ģ‚š‘„š‘‡š‘˜|š‘˜āˆ’1+š‘„š‘˜.(B.5)

Measurement Update Step
(1)Factorize: š‘ƒš‘˜āˆ’1|š‘˜āˆ’1=ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1ī‚€ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1ī‚š‘‡.(B.6)(2)Evaluate the quadrature points {š‘‹š‘™,š‘˜āˆ’1|š‘˜āˆ’1}š‘šš‘™=1 as: š‘‹š‘™,š‘˜āˆ’1|š‘˜āˆ’1=ī”š‘ƒš‘˜āˆ’1|š‘˜āˆ’1šœ‰š‘™+Ģ‚š‘„š‘˜āˆ’1|š‘˜āˆ’1.(B.7)(3)Evaluate the propagated quadrature points {š‘š‘™,š‘˜|š‘˜āˆ’1}š‘šš‘™=1 as: š‘š‘™,š‘˜|š‘˜āˆ’1ī€·š‘‹=ā„Žš‘™,š‘˜|š‘˜āˆ’1,š‘¢š‘˜ī€ø,š‘˜.(B.8)(4)Estimate the predicted measurement: Ģ‚š‘§š‘˜|š‘˜āˆ’1=ī“š‘šš‘™=1š‘¤š‘™š‘š‘™,š‘˜|š‘˜āˆ’1.(B.9)(5)Estimate the innovation covariance matrix š‘ƒš‘§š‘§,š‘˜|š‘˜āˆ’1=š‘…š‘˜+ī“š‘šš‘™=1š‘¤š‘™š‘š‘™,š‘˜|š‘˜āˆ’1š‘š‘‡š‘™,š‘˜|š‘˜āˆ’1āˆ’Ģ‚š‘§š‘˜|š‘˜āˆ’1Ģ‚š‘§š‘‡š‘˜|š‘˜āˆ’1.(B.10)(6)Estimate the cross covariance matrix š‘ƒš‘„š‘§,š‘˜|š‘˜āˆ’1=ī“š‘šš‘™=1š‘¤š‘™š‘‹š‘™,š‘˜|š‘˜āˆ’1š‘š‘‡š‘™,š‘˜|š‘˜āˆ’1āˆ’Ģ‚š‘„š‘˜|š‘˜āˆ’1Ģ‚š‘§š‘‡š‘˜|š‘˜āˆ’1.(B.11)(7)Estimate the Kalman gain: š‘Šš‘˜=š‘ƒš‘„š‘§,š‘˜|š‘˜āˆ’1š‘ƒāˆ’1š‘§š‘§,š‘˜|š‘˜āˆ’1.(B.12)(8)Estimate the updated step: Ģ‚š‘„š‘˜|š‘˜=Ģ‚š‘„š‘˜|š‘˜āˆ’1+š‘Šš‘˜ī€·š‘§š‘˜āˆ’Ģ‚š‘§š‘˜|š‘˜āˆ’1ī€ø.(B.13)(9)Estimate the corresponding error covariance: š‘ƒš‘˜|š‘˜=š‘ƒš‘˜|š‘˜āˆ’1āˆ’š‘Šš‘˜š‘ƒš‘§š‘§,š‘˜|š‘˜āˆ’1š‘Šš‘‡š‘˜.(B.14)

C. Derivation of the RC Battery Model Equations

The output voltage may be calculated by summing the voltages of each element in the circuit. Summation of the outer loop and inner loop voltages yield two equations for the output voltage, respectively: š‘‰š‘‚=š¼š‘†š‘…š‘”+š¼š‘š‘…š‘’+š‘‰š¶š‘,š‘‰(C.1)š‘‚=š¼š‘†š‘…š‘”+š¼š‘š‘…š‘+š‘‰š¶š‘.(C.2)

Recall Kirchhoffā€™s current law, which states that the total current entering a node must equal the total current leaving, such that the supply current is defined by š¼š‘†=š¼š‘+š¼š‘.(C.3) Rearranging for š¼š‘ in (C.3) and substituting into (C.2) yields š‘‰š‘‚=š¼š‘†š‘…š‘”+ī€·š¼š‘†āˆ’š¼š‘ī€øš‘…š‘+š‘‰š¶š‘.(C.4) Equating (C.4) with the first output voltage (C.1) and solving for š¼š‘ yields š¼š‘=š¼š‘†š‘…š‘š‘…š‘’+š‘…š‘+š‘‰š¶š‘š‘…š‘’+š‘…š‘āˆ’š‘‰š¶š‘š‘…š‘’+š‘…š‘.(C.5) Note that Ģ‡š‘‰š¶š‘=š¼š‘/š¶š‘ such that (C.5) becomes Ģ‡š‘‰š¶š‘=š¼š‘†š‘…š‘š¶š‘ī€·š‘…š‘’+š‘…š‘ī€ø+š‘‰š¶š‘š¶š‘ī€·š‘…š‘’+š‘…š‘ī€øāˆ’š‘‰š¶š‘š¶š‘ī€·š‘…š‘’+š‘…š‘ī€ø.(C.6) Similarly, the above approach may be used to solve for the second capacitor voltage rate of change: Ģ‡š‘‰š¶š‘=š¼š‘†š‘…š‘’š¶š‘ī€·š‘…š‘’+š‘…š‘ī€ø+š‘‰š¶š‘š¶š‘ī€·š‘…š‘’+š‘…š‘ī€øāˆ’š‘‰š¶š‘š¶š‘ī€·š‘…š‘’+š‘…š‘ī€ø.(C.7) The output voltage may be solved by manipulating (C.1) through (C.3). Setting š¼š‘=š¼š‘†āˆ’š¼š‘ in (C.2), and rearranging for š¼š‘, results in the following: š¼š‘=š‘…š‘”+š‘…š‘š‘…š‘š¼š‘†+š‘‰š¶š‘š‘…š‘āˆ’š‘‰š‘‚š‘…š‘.(C.8) Rearranging (C.1) for š¼š‘,and setting the result equal to (C.8) yields š‘‰š‘‚š‘…š‘’āˆ’š‘‰š¶š‘š‘…š‘’āˆ’š‘…š‘”š‘…š‘’š¼š‘†=š‘…š‘”+š‘…š‘š‘…š‘š¼š‘†+š‘‰š¶š‘š‘…š‘āˆ’š‘‰š‘‚š‘…š‘.(C.9) Manipulating (C.9) and simplifying for the output voltage results in š‘‰š‘‚=š‘…š‘š‘…š‘’+š‘…š‘š‘‰š¶š‘+š‘…š‘’š‘…š‘’+š‘…š‘š‘‰š¶š‘+ī‚øš‘…š‘”+š‘…š‘’š‘…š‘š‘…š‘’+š‘…š‘ī‚¹š¼š‘†.(C.10)

D. Linearization of the Battery System

The battery model described in Section 3 requires linearization for the EKF to be implemented. The linearization was based on the Jacobian matrix and the four states (where for neatness, š‘Šš¶š‘=1/š¶š‘ and š‘Šš¶š‘=1/š¶š‘) šœ™š‘˜=āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£šœ•š‘“1šœ•š‘‰š¶š‘šœ•š‘“1šœ•š‘‰š¶š‘šœ•š‘“1šœ•š‘Šš¶š‘šœ•š‘“1šœ•š‘Šš¶š‘šœ•š‘“2šœ•š‘‰š¶š‘šœ•š‘“2šœ•š‘‰š¶š‘šœ•š‘“2šœ•š‘Šš¶š‘šœ•š‘“2šœ•š‘Šš¶š‘šœ•š‘“3šœ•š‘‰š¶š‘šœ•š‘“3šœ•š‘‰š¶š‘šœ•š‘“3šœ•š‘Šš¶š‘šœ•š‘“3šœ•š‘Šš¶š‘šœ•š‘“4šœ•š‘‰š¶š‘šœ•š‘“4šœ•š‘‰š¶š‘šœ•š‘“4šœ•š‘Šš¶š‘šœ•š‘“4šœ•š‘Šš¶š‘āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š‘˜,(D.1) where the functions in (D.1) are described by š‘“1š‘‡=āˆ’š‘†š‘‰š¶š‘š‘Šš¶š‘š‘…š‘’+š‘…š‘+š‘‰š¶š‘+š‘‡š‘†š‘‰š¶š‘š‘Šš¶š‘š‘…š‘’+š‘…š‘+š‘‡š‘†š‘…š‘š‘Šš¶š‘š‘…š‘’+š‘…š‘š¼š‘†,š‘“2š‘‡=āˆ’š‘†š‘‰š¶š‘š‘Šš¶š‘š‘…š‘’+š‘…š‘āˆ’š‘‡š‘†š‘‰š¶š‘š‘Šš¶š‘š‘…š‘’+š‘…š‘+š‘‰š¶š‘+š‘‡š‘†š‘…š‘’š‘Šš¶š‘š‘…š‘’+š‘…š‘š¼š‘†,š‘“3=š‘Šš¶š‘,š‘“4=š‘Šš¶š‘.(D.2)

The linearized form (šœ™) of the system matrix, using the state vector described by (18), is as follows: šœ™š‘˜=āŽ”āŽ¢āŽ¢āŽ¢āŽ£šœ™11šœ™12šœ™130šœ™21šœ™220šœ™24āŽ¤āŽ„āŽ„āŽ„āŽ¦00100001š‘˜,(D.3) where the elements are described by: šœ™11=āˆ’š‘‡š‘†š‘Šš‘˜š¶š‘š‘…š‘’š‘˜+š‘…š‘š‘˜šœ™+1,12=š‘‡š‘†š‘Šš‘˜š¶š‘š‘…š‘’š‘˜+š‘…š‘š‘˜,šœ™13=š‘‡š‘†š‘…š‘’š‘˜+š‘…š‘š‘˜ī€·āˆ’š‘‰š‘˜š¶š‘+š‘‰š‘˜š¶š‘+š‘…š‘š‘˜š¼š‘†š‘˜ī€ø,šœ™21=š‘‡š‘†š‘Šš‘˜š¶š‘š‘…š‘’š‘˜+š‘…š‘š‘˜,šœ™22=āˆ’š‘‡š‘†š‘Šš‘˜š¶š‘š‘…š‘’š‘˜+š‘…š‘š‘˜šœ™+1,24=š‘‡š‘†š‘…š‘’š‘˜+š‘…š‘š‘˜ī€·š‘‰š‘˜š¶š‘āˆ’š‘‰š‘˜š¶š‘+š‘…š‘’š‘˜š¼š‘†š‘˜ī€ø.(D.4)


  1. T. L. Churchill, ā€œBattery condition management: an important way to protect a critical asset,ā€ IEEE Aerospace and Electronic Systems Magazine, vol. 15, no. 7, pp. 41ā€“47, 2000. View at: Publisher Site | Google Scholar
  2. F. V. Conte, ā€œBattery and battery management for hybrid electric vehicles: a review,ā€ Elektrotechnik und Informationstechnik, vol. 123, no. 10, pp. 424ā€“431, 2006. View at: Publisher Site | Google Scholar
  3. I. S. Kim, ā€œNonlinear state of charge estimator for hybrid electric vehicle battery,ā€ IEEE Transactions on Power Electronics, vol. 23, no. 4, pp. 2027ā€“2034, 2008. View at: Publisher Site | Google Scholar
  4. A. Piccolo, A. Vaccaro, and D. Villacci, ā€œFuzzy logic based optimal power flow management in parallel hybrid electric vehicles,ā€ Iranian Journal of Electrical and Computer Engineering, vol. 4, no. 2, p. 83, 2005. View at: Google Scholar
  5. X. Li, J. Li, L. Xu, and M. Ouyang, ā€œPower management and economic estimation of fuel cell hybrid vehicle using fuzzy logic,ā€ in Proceedings of the 5th IEEE Vehicle Power and Propulsion Conference (VPPC '09), pp. 1749ā€“1754, September 2009. View at: Publisher Site | Google Scholar
  6. P. Naderi, M. Mirsalim, M. T. Bathaee, and R. Chini, ā€œFuzzy controller design for parallel hybrid vehicle analysis using forward simulation,ā€ in Proceedings of the 5th IEEE Vehicle Power and Propulsion Conference (VPPC '09), pp. 234ā€“241, September 2009. View at: Publisher Site | Google Scholar
  7. M. Montazeri-Gh, A. Poursamad, and B. Ghalichi, ā€œApplication of genetic algorithm for optimization of control strategy in parallel hybrid electric vehicles,ā€ Journal of the Franklin Institute, vol. 343, no. 4-5, pp. 420ā€“435, 2006. View at: Publisher Site | Google Scholar
  8. L. V. PĆ©rez, G. R. Bossio, D. Moitre, and G. O. GarcĆ­a, ā€œOptimization of power management in an hybrid electric vehicle using dynamic programming,ā€ Mathematics and Computers in Simulation, vol. 73, no. 1ā€“4, pp. 244ā€“254, 2006. View at: Publisher Site | Google Scholar
  9. B. Sampathnarayanan, L. Serrao, S. Onori, G. Rizzoni, and S. Yurkovich, ā€œModel predictive control as an energy management strategy for hybrid electric vehicles,ā€ in Proceedings of the ASME Dynamic Systems and Control Conference (DSCC '09), pp. 1161ā€“1168, October 2009. View at: Google Scholar
  10. D. F. Opila, X. Wang, R. McGee, J. A. Cook, and J. W. Grizzle, ā€œFundamental structural limitations of an industrial energy management controller architecture for hybrid vehicles,ā€ in Proceedings of the ASME Dynamic Systems and Control Conference (DSCC '09), pp. 213ā€“221, October 2009. View at: Google Scholar
  11. K. D. Huang, S. C. Tzeng, T. M. Jeng, and C. C. Chen, ā€œIntegration mechanism for a parallel hybrid vehicle system,ā€ Applied Energy, vol. 82, no. 2, pp. 133ā€“147, 2005. View at: Publisher Site | Google Scholar
  12. F. CodecĆ , S. M. Savaresi, and V. Manzoni, ā€œThe mix estimation algorithm for battery state-of-charge estimator-analysis of the sensitivity to model errors,ā€ in Proceedings of the ASME Dynamic Systems and Control Conference (DSCC '09), pp. 97ā€“104, October 2009. View at: Google Scholar
  13. B. S. Bhangu, P. Bentley, D. A. Stone, and C. M. Bingham, ā€œNonlinear observers for predicting state-of-charge and state-of-health of lead-acid batteries for hybrid-electric vehicles,ā€ IEEE Transactions on Vehicular Technology, vol. 54, no. 3, pp. 783ā€“794, 2005. View at: Publisher Site | Google Scholar
  14. V. Marano, S. Onori, Y. Guezennec, G. Rizzoni, and N. Madella, ā€œLithium-ion batteries life estimation for plug-in hybrid electric vehicles,ā€ in Proceedings of the 5th IEEE Vehicle Power and Propulsion Conference (VPPC '09), pp. 536ā€“543, September 2009. View at: Publisher Site | Google Scholar
  15. G. L. Plett, ā€œExtended Kalman filtering for battery management systems of LiPB-based HEV battery packsā€”part 1. Background,ā€ Journal of Power Sources, vol. 134, no. 2, pp. 252ā€“261, 2004. View at: Publisher Site | Google Scholar
  16. G. L. Plett, ā€œExtended Kalman filtering for battery management systems of LiPB-based HEV battery packsā€”part 2. Modeling and identification,ā€ Journal of Power Sources, vol. 134, no. 2, pp. 262ā€“276, 2004. View at: Publisher Site | Google Scholar
  17. G. L. Plett, ā€œExtended Kalman filtering for battery management systems of LiPB-based HEV battery packsā€”part 3. State and parameter estimation,ā€ Journal of Power Sources, vol. 134, no. 2, pp. 277ā€“292, 2004. View at: Publisher Site | Google Scholar
  18. A. Vasebi, M. Partovibakhsh, and S. M. T. Bathaee, ā€œA novel combined battery model for state-of-charge estimation in lead-acid batteries based on extended Kalman filter for hybrid electric vehicle applications,ā€ Journal of Power Sources, vol. 174, no. 1, pp. 30ā€“40, 2007. View at: Publisher Site | Google Scholar
  19. A. Vasebi, S. M. T. Bathaee, and M. Partovibakhsh, ā€œPredicting state of charge of lead-acid batteries for hybrid electric vehicles by extended Kalman filter,ā€ Energy Conversion and Management, vol. 49, no. 1, pp. 75ā€“82, 2008. View at: Publisher Site | Google Scholar
  20. G. L. Plett, ā€œSigma-point Kalman filtering for battery management systems of LiPB-based HEV battery packsā€”part 1: introduction and state estimation,ā€ Journal of Power Sources, vol. 161, no. 2, pp. 1356ā€“1368, 2006. View at: Publisher Site | Google Scholar
  21. G. L. Plett, ā€œSigma-point Kalman filtering for battery management systems of LiPB-based HEV battery packsā€”part 2: simultaneous state and parameter estimation,ā€ Journal of Power Sources, vol. 161, no. 2, pp. 1369ā€“1384, 2006. View at: Publisher Site | Google Scholar
  22. A. Widodo and B. S. Yang, ā€œSupport vector machine in machine condition monitoring and fault diagnosis,ā€ Mechanical Systems and Signal Processing, vol. 21, no. 6, pp. 2560ā€“2574, 2007. View at: Publisher Site | Google Scholar
  23. Y. Bar-Shalom, X. Rong Li, and T. Kirubarajan, Estimation with Applications to Tracking and Navigation, John Wiley and Sons, New York, NY, USA, 2001.
  24. N. Nise, Control Systems Engineering, John Wiley and Sons, New York, NY, USA, 4th edition, 2004.
  25. R. E. Kalman, ā€œA new approach to linear filtering and prediction problems,ā€ Transactions of the ASME Journal of Basic Engineering, pp. 35ā€“45, 1960. View at: Google Scholar
  26. M. S. Grewal and A. P. Andrews, Kalman Filtering: Theory and Practice Using MATLAB, John Wiley and Sons, New York, NY, USA, 3rd edition, 2008.
  27. S. Habibi, R. Burton, and Y. Chinniah, ā€œEstimation using a new variable structure filter,ā€ in Proceedings of the American Control Conference, pp. 2937ā€“2942, May 2002. View at: Google Scholar
  28. P. Del Moral, ā€œMeasure-valued processes and interacting particle systems. Application to nonlinear filtering problems 1,ā€ Annals of Applied Probability, vol. 8, no. 2, pp. 438ā€“495, 1998. View at: Google Scholar
  29. N. J. Gordon, D. J. Salmond, and A. F. M. Smith, ā€œNovel approach to nonlinear/non-gaussian Bayesian state estimation,ā€ IEE Proceedings. Part F, vol. 140, no. 2, pp. 107ā€“113, 1993. View at: Google Scholar
  30. J. MacCormick and A. Blake, ā€œProbabilistic exclusion principle for tracking multiple objects,ā€ in Proceedings of the 1999 7th IEEE International Conference on Computer Vision (ICCV'99), pp. 572ā€“578, September 1999. View at: Google Scholar
  31. K. Kanazawa, D. Koller, and S. J. Russell, ā€œStochastic simulation algorithms for dynamic probabilistic networks,ā€ in Proceedings of the 11th Annual Conference on Uncertainty in Artificial Intelligence, pp. 346ā€“351, 1995. View at: Google Scholar
  32. B. Ristic, S. Arulampalam, and N. Gordon, Beyond the Kalman Filter: Particle Filters for Tracking Applications, Artech House, Boston, Mass, USA, 2004.
  33. I. Arasaratnam, S. Haykin, and R. J. Elliott, ā€œDiscrete-time nonlinear filtering algorithms using gauss-hermite quadrature,ā€ Proceedings of the IEEE, vol. 95, no. 5, pp. 953ā€“977, 2007. View at: Publisher Site | Google Scholar
  34. I. Arasaratnam and S. Haykin, ā€œSquare-root quadrature Kalman filtering,ā€ IEEE Transactions on Signal Processing, vol. 56, no. 6, pp. 2589ā€“2593, 2008. View at: Publisher Site | Google Scholar
  35. S. R. Habibi and R. Burton, ā€œThe variable structure filter,ā€ Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, vol. 125, no. 3, pp. 287ā€“293, 2003. View at: Publisher Site | Google Scholar
  36. S. Habibi, ā€œThe smooth variable structure filter,ā€ Proceedings of the IEEE, vol. 95, no. 5, pp. 1026ā€“1059, 2007. View at: Publisher Site | Google Scholar
  37. T. Okoshi, K. Yamada, T. Hirasawa, and A. Emori, ā€œBattery condition monitoring (BCM) technologies about lead-acid batteries,ā€ Journal of Power Sources, vol. 158, no. 2, pp. 874ā€“878, 2006. View at: Publisher Site | Google Scholar
  38. J. Voelcker, ā€œLithium batteries take to the road,ā€ IEEE Spectrum, vol. 44, no. 9, pp. 26ā€“31, 2007. View at: Publisher Site | Google Scholar
  39. ADVISOR 2004 User's Guide, Advanced Vehicle Labs, 2004.
  40. T. Markel, A. Brooker, T. Hendricks et al., ā€œADVISOR: a systems analysis tool for advanced vehicle modeling,ā€ Journal of Power Sources, vol. 110, no. 2, pp. 255ā€“266, 2002. View at: Publisher Site | Google Scholar
  41. V. H. Johnson, ā€œBattery performance models in ADVISOR,ā€ Journal of Power Sources, vol. 110, no. 1, pp. 321ā€“329, 2002. View at: Google Scholar
  42. V. H. Johnson, M. Zolot, and A. Pesaran, ā€œDevelopment and validation of a temperature-dependent resistance/capacitance battery model for ADVISOR,ā€ in Proceedings of the 18th Electric Vehicle Symposium, Berlin, Germany, 2001. View at: Google Scholar
  43. HEV Implementing Agreement, International Energy Agency, 2008.

Copyright Ā© 2011 S. A. Gadsden 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.

1090Ā Views | 642Ā Downloads | 10Ā Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder