[go: up one dir, main page]

CN108805217A - A kind of health state of lithium ion battery method of estimation and system based on support vector machines - Google Patents

A kind of health state of lithium ion battery method of estimation and system based on support vector machines Download PDF

Info

Publication number
CN108805217A
CN108805217A CN201810636667.9A CN201810636667A CN108805217A CN 108805217 A CN108805217 A CN 108805217A CN 201810636667 A CN201810636667 A CN 201810636667A CN 108805217 A CN108805217 A CN 108805217A
Authority
CN
China
Prior art keywords
regression
set data
training
support vector
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810636667.9A
Other languages
Chinese (zh)
Other versions
CN108805217B (en
Inventor
崔纳新
方浩然
杨亚宁
王春雨
王光臣
张承慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN201810636667.9A priority Critical patent/CN108805217B/en
Publication of CN108805217A publication Critical patent/CN108805217A/en
Application granted granted Critical
Publication of CN108805217B publication Critical patent/CN108805217B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Molecular Biology (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

The invention discloses a kind of health state of lithium ion battery method of estimation and system based on support vector machines, including:Determine the input variable and output variable of Support vector regression prediction;Input variable and output variable are divided into training set data group and test set data group;Regression model foundation is carried out to the training set data after normalization, obtains regression function;Test set data are brought into regression model after training, to predict the electricity being filled with when constant-current charging of battery reaches blanking voltage;The electricity that constant-current charge in training set data to blanking voltage is filled with is fitted with the current test capacity obtained after volume test, the equation that the electricity that test set data prediction obtains is brought into after fitting is obtained into current predictive capacity, to estimate cell health state.The present invention can reach the electricity that blanking voltage is filled with to the lithium ion battery constant-current charge under various constant-current charge environment and predict there is wide applicability.

Description

一种基于支持向量机的锂离子电池健康状态估计方法及系统A method and system for estimating the state of health of a lithium-ion battery based on a support vector machine

技术领域technical field

本发明属于锂离子电池健康状态估计技术领域,尤其涉及一种基于支持向量机的锂离子电池健康状态估计方法及系统。The invention belongs to the technical field of estimation of the state of health of lithium-ion batteries, and in particular relates to a method and system for estimating the state of health of lithium-ion batteries based on a support vector machine.

背景技术Background technique

电动汽车在节能减排和环境保护方面与燃油车相比有较大的优势,因此各国争相出台各种激励措施来大力推动电动汽车的发展。而锂电池具有比能量高、使用寿命长、额定电压高、具备高功率承受力、自放电率低等优点,自上市以来就被电动汽车厂所青睐。电动汽车的快速发展,使得人们对锂电池的性能提出了更高的要求,从而锂电池的健康状态(State Of Health;SOH)估计备受关注。由于电池老化机理的复杂性,使得电池SOH的快速、准确估计具有难度。Compared with fuel vehicles, electric vehicles have greater advantages in terms of energy saving, emission reduction and environmental protection. Therefore, various countries are scrambling to introduce various incentives to vigorously promote the development of electric vehicles. Lithium batteries have the advantages of high specific energy, long service life, high rated voltage, high power tolerance, and low self-discharge rate, and have been favored by electric vehicle manufacturers since their listing. The rapid development of electric vehicles has made people put forward higher requirements on the performance of lithium batteries, so the state of health (State Of Health; SOH) estimation of lithium batteries has attracted much attention. Due to the complexity of the battery aging mechanism, it is difficult to quickly and accurately estimate the battery SOH.

判断SOH的典型方法是依据容量减少或者电阻增加来估计其值。其中,扩展卡尔曼滤波方法被认为是最可靠的在线估计SOH的方法。但该方法的准确度依赖于所建模型的参数,精度容易受到模型准确度的影响。因此,方便、简单准确的获取SOH值的方法备受关注。A typical method of judging SOH is to estimate its value based on a decrease in capacity or an increase in resistance. Among them, the extended Kalman filter method is considered to be the most reliable method for estimating SOH online. However, the accuracy of this method depends on the parameters of the built model, and the accuracy is easily affected by the accuracy of the model. Therefore, a convenient, simple and accurate method for obtaining the SOH value has attracted much attention.

由电池循环数据可知,在相同充电条件下,电池恒流充电到达截止电压所充入的电量Q随SOH的降低而减少且两者存在较高的关联度。由于实际充电情况,锂电池无法在每次恒流充电时都将端电压充电到截止电压。因此,准确预测相同充电条件下电池恒流充电到达截止电压所充入的电量Q成为预测锂电池SOH值的关键一步。同时,由于电池充电系统的迟滞性,锂电池在恒流充电过程中到达充电截止电压之后,各充电相关设备才会做出动作,必然会造成电池的过冲并对电池造成损害,从而缩短电池的使用寿命。为解决上述问题,需要提前预测出电池恒流充电到达截止电压所充入的电量Q以防止端电压超过截止电压并进一步做到对电池SOH的估计。It can be seen from the battery cycle data that under the same charging conditions, the amount of electricity Q charged by constant current charging to the cut-off voltage of the battery decreases with the decrease of SOH and there is a high degree of correlation between the two. Due to the actual charging situation, the lithium battery cannot charge the terminal voltage to the cut-off voltage every time the constant current is charged. Therefore, accurately predicting the charge Q of the battery constant current charging to the cut-off voltage under the same charging conditions becomes a key step in predicting the SOH value of lithium batteries. At the same time, due to the hysteresis of the battery charging system, after the lithium battery reaches the charging cut-off voltage in the constant current charging process, the charging related equipment will take action, which will inevitably cause the battery to overshoot and cause damage to the battery, thereby shortening the battery life. service life. In order to solve the above problems, it is necessary to predict in advance the charge Q of the battery constant current charging to the cut-off voltage to prevent the terminal voltage from exceeding the cut-off voltage and to further estimate the battery SOH.

SOH一般定义为其中cM为电池当前的测试容量,c0为电池额定容量。由于测试cM值时的试验周期长并且会造成浪费,因此准确、可靠、方便的估计电池的SOH是电池管理系统的重要任务。SOH is generally defined as Among them, c M is the current test capacity of the battery, and c 0 is the rated capacity of the battery. Since the test cycle is long and wasteful when testing the c M value, it is an important task of the battery management system to estimate the SOH of the battery accurately, reliably and conveniently.

发明内容Contents of the invention

本发明的目的就是为了解决上述问题,提出一种基于支持向量机算法的锂离子电池健康状态估计方法及系统。该方法根据锂电池工作的历史数据、工作状态及环境数据由支持向量机对锂电池恒流充电到达截止电压所充入的电量Q做出准确预测,同时使用网格搜寻法对支持向量机的参数选择进行了相应的优化,提高了参数选择的准确度。由实验可知,相同充电条件下,锂电池恒流充电到达截止电压所充入的电量随着SOH值的降低而减少且两者存在较高的关联度。因此根据之前测试的相关数据,可由最小二乘法拟合测试集恒流充电电量Q11与SOH值两者之间的关系,从而由上述支持向量机算法预测的恒流充电电量Q1对当前锂电池的SOH做出估计。The object of the present invention is to solve the above-mentioned problems, and propose a method and system for estimating the state of health of a lithium-ion battery based on a support vector machine algorithm. In this method, based on the historical data, working status and environmental data of the lithium battery, the support vector machine makes an accurate prediction of the amount of electricity Q charged by the constant current charging of the lithium battery to the cut-off voltage. The parameter selection is optimized accordingly, which improves the accuracy of parameter selection. It can be seen from the experiment that under the same charging conditions, the amount of electricity charged by constant current charging of lithium batteries to the cut-off voltage decreases as the SOH value decreases, and there is a high degree of correlation between the two. Therefore, according to the relevant data of the previous test, the relationship between the test set constant-current charging quantity Q11 and the SOH value can be fitted by the least square method, so that the constant-current charging quantity Q11 predicted by the above-mentioned support vector machine algorithm has a significant impact on the current lithium battery. The SOH of the battery is estimated.

为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:

本发明的第一目的是公开一种基于支持向量机的锂离子电池健康状态估计方法,包括:The first object of the present invention is to disclose a method for estimating the state of health of a lithium-ion battery based on a support vector machine, comprising:

对锂电池进行循环充放电实验,实时记录锂电池各项工作状态的历史数据,确定支持向量机回归预测的输入变量和输出变量;Carry out cyclic charging and discharging experiments on lithium batteries, record the historical data of various working states of lithium batteries in real time, and determine the input variables and output variables for support vector machine regression prediction;

将输入变量和输出变量分成两组:训练集数据组和测试集数据组;Divide the input variables and output variables into two groups: the training set data group and the test set data group;

对训练集数据组和测试集数据组进行归一化处理;Normalize the training set data set and the test set data set;

对归一化后的训练集数据进行回归模型建立,得到回归函数;Establish a regression model on the normalized training set data to obtain a regression function;

选取RBF作为核函数,运用网格搜寻法选择最优的RBF核函数参数组合:核函数的宽度参数,惩罚系数和损失函数;Select RBF as the kernel function, and use the grid search method to select the optimal combination of RBF kernel function parameters: the width parameter of the kernel function, the penalty coefficient and the loss function;

根据最优的RBF核函数参数组合确定最优参数回归模型;Determine the optimal parameter regression model according to the optimal RBF kernel function parameter combination;

将测试集数据带入训练后的回归模型,从而预测出电池恒流充电到达截止电压时所充入的电量Q1Bring the test set data into the trained regression model to predict the charge Q 1 when the constant current charging of the battery reaches the cut-off voltage;

在相同恒流充电条件下,将训练集数据中的恒流充电到截止电压所充入的电量Q11和进行容量测试后得到的当前测试容量Q0运用最小二乘法进行拟合,将测试集数据预测得到的Q1带入拟合后的方程得到当前预测容量cM,从而对电池健康状态进行估计。Under the same constant-current charging condition, use the least squares method to fit the electric quantity Q 11 charged from the constant-current charge in the training set data to the cut-off voltage and the current test capacity Q 0 obtained after the capacity test, and use the least squares method to fit the test set The Q 1 obtained from the data prediction is brought into the fitted equation to obtain the current predicted capacity c M , thereby estimating the state of health of the battery.

进一步地,所述确定支持向量机回归预测的输入变量和输出变量具体为:Further, the input variables and output variables for determining the regression prediction of the support vector machine are specifically:

实时记录锂电池各项工作状态的历史数据,包括循环次数、循环充放电电流、放电深度、温度和下次将要进行的恒流充电的电流,作为支持向量机回归预测的输入变量;Record the historical data of various working states of lithium batteries in real time, including the number of cycles, cycle charge and discharge current, discharge depth, temperature and the current of the next constant current charge, as input variables for support vector machine regression prediction;

在设定的循环周期后,对锂电池进行一次恒流充电和容量测试,统计电池恒流充电到达截止电压时所充入的电量Q作为支持向量机回归预测的输出变量。After the set cycle period, conduct a constant current charging and capacity test on the lithium battery, and count the charged power Q when the battery constant current charging reaches the cut-off voltage as the output variable of the support vector machine regression prediction.

进一步地,所述对训练集数据组和测试集数据组进行归一化处理,具体为:Further, the normalization processing of the training set data group and the test set data group is specifically:

归一化后的样本值为该样本值与样本的最小值之差与该样本值最大值与最小值之差的比值。The normalized sample value is the ratio of the difference between the sample value and the minimum value of the sample to the difference between the maximum value and the minimum value of the sample value.

进一步地,当参数训练预测模型的性能相同时,优先选择惩罚系数相对较小的参数组合。Furthermore, when the performance of the parameter training prediction model is the same, the parameter combination with a relatively small penalty coefficient is preferred.

进一步地,所述对归一化后的训练集数据进行回归模型建立,得到回归函数;具体为:Further, the regression model is established on the normalized training set data to obtain a regression function; specifically:

设含有m个训练样本的训练集样本对为其中,是第i个训练样本的输入列向量, 为对应的输出值;Let the training set sample pair containing m training samples be in, is the input column vector of the ith training sample, is the corresponding output value;

设在高维特征空间中建立的线性回归函数为:The linear regression function established in the high-dimensional feature space is:

定义ε线性不敏感损失函数为:Define the ε linear insensitive loss function as:

为了寻找到一个最优分类面使得所有训练样本离该最优分类面的误差最小,引入松弛变量ξii *,则需要满足如下约束关系式:In order to find an optimal classification surface so that all training samples have the smallest error from the optimal classification surface, the slack variables ξ i , ξ i * need to be satisfied as follows:

其中,为非线性映射函数,为回归系数向量,b为阈值;为回归函数返回的预测值,y为对应的真实值,ε为设定的大于0的数;C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大,ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小。in, is a nonlinear mapping function, is the regression coefficient vector, b is the threshold; is the predicted value returned by the regression function, y is the corresponding real value, ε is a set number greater than 0; C is the penalty factor, and the larger the C, the greater the penalty for the sample whose training error is greater than ε, and ε specifies the regression function The smaller the ε, the smaller the error of the regression function.

进一步地,对所述约束关系式进行求解,具体为:Further, the constraint relational expression is solved, specifically:

采用拉格朗日对偶理论将约束关系式转化为对偶问题;Using the Lagrangian dual theory to transform the constraint relation into a dual problem;

由拉格朗日对偶性,将原拉格朗日函数中最小值问题转化为最大化问题;By the Lagrangian duality, the minimum problem in the original Lagrangian function is transformed into a maximum problem;

采用SMO算法求解对偶问题中的拉格朗日乘子;进而得到回归系数向量和阈值b,确定回归函数。Use the SMO algorithm to solve the Lagrangian multipliers in the dual problem; and then obtain the regression coefficient vector and the threshold b to determine the regression function.

进一步地,所述回归函数具体为:Further, the regression function is specifically:

其中,分别为样本i对应的拉格朗日乘子的最优解,为核函数,b*为根据拉格朗日乘子的最优解求解得到的阈值。in, are the optimal solutions of the Lagrangian multipliers corresponding to sample i, respectively, is the kernel function, and b * is the threshold obtained by solving the optimal solution of Lagrangian multipliers.

进一步地,所述选取RBF作为核函数,运用网格搜寻法选择最优的RBF核函数参数组合;具体为:Further, the RBF is selected as the kernel function, and the grid search method is used to select the optimal RBF kernel function parameter combination; specifically:

步骤一:对各种可能的RBF核函数参数组合值进行交叉验证,找出使交叉验证精确度最高的组合对;将各个参数可能的取值进行排列组合,列出所有可能的组合结果生成“网格”;Step 1: Carry out cross-validation on various possible RBF kernel function parameter combination values, and find out the combination pair with the highest cross-validation accuracy; arrange and combine the possible values of each parameter, and list all possible combination results to generate " grid";

步骤二:将各RBF核函数参数组合用于SVM训练,并使用交叉验证对表现进行评估;Step 2: Combine the parameters of each RBF kernel function for SVM training, and use cross-validation to evaluate the performance;

步骤三:在拟合函数尝试了所有的参数组合后,返回一个最优分类面,自动调整至最佳参数组合,使得所有训练样本离该最优分类面的误差最小。Step 3: After the fitting function has tried all the parameter combinations, return an optimal classification surface, and automatically adjust to the optimal parameter combination, so that the error of all training samples from the optimal classification surface is minimized.

本发明的第二目的是公开一种基于支持向量机的锂离子电池健康状态估计系统,包括服务器,所述服务器包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:The second object of the present invention is to disclose a system for estimating the state of health of a lithium-ion battery based on a support vector machine, including a server, the server including a memory, a processor and a computer program stored on the memory and operable on the processor, When the processor executes the program, the following steps are implemented:

确定支持向量机回归预测的输入变量和输出变量;Determine input variables and output variables for SVM regression prediction;

将输入变量和输出变量分成两组:训练集数据组和测试集数据组;Divide the input variables and output variables into two groups: the training set data group and the test set data group;

对训练集数据组和测试集数据组进行归一化处理;Normalize the training set data set and the test set data set;

对归一化后的训练集数据进行回归模型建立,得到回归函数;Establish a regression model on the normalized training set data to obtain a regression function;

选取RBF作为核函数,运用网格搜寻法选择最优的RBF核函数参数组合:核函数的宽度参数,惩罚系数和损失函数;Select RBF as the kernel function, and use the grid search method to select the optimal combination of RBF kernel function parameters: the width parameter of the kernel function, the penalty coefficient and the loss function;

根据最优的RBF核函数参数组合确定最优参数回归模型;Determine the optimal parameter regression model according to the optimal RBF kernel function parameter combination;

将测试集数据带入训练后的回归模型,从而预测出电池恒流充电到达截止电压时所充入的电量Q1Bring the test set data into the trained regression model to predict the charge Q 1 when the constant current charging of the battery reaches the cut-off voltage;

在相同恒流充电条件下,将训练集数据中的恒流充电到截止电压所充入的电量Q11和进行容量测试后得到的当前测试容量Q0运用最小二乘法进行拟合,将测试集数据预测得到的Q1带入拟合后的方程得到当前预测容量cM,从而对电池健康状态进行估计。Under the same constant-current charging condition, use the least squares method to fit the electric quantity Q 11 charged from the constant-current charge in the training set data to the cut-off voltage and the current test capacity Q 0 obtained after the capacity test, and use the least squares method to fit the test set The Q 1 obtained from the data prediction is brought into the fitted equation to obtain the current predicted capacity c M , thereby estimating the state of health of the battery.

本发明的第三目的是公开一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时执行以下步骤:The third object of the present invention is to disclose a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the following steps are performed:

确定支持向量机回归预测的输入变量和输出变量;Determine input variables and output variables for SVM regression prediction;

将输入变量和输出变量分成两组:训练集数据组和测试集数据组;Divide the input variables and output variables into two groups: the training set data group and the test set data group;

对训练集数据组和测试集数据组进行归一化处理;Normalize the training set data set and the test set data set;

对归一化后的训练集数据进行回归模型建立,得到回归函数;Establish a regression model on the normalized training set data to obtain a regression function;

选取RBF作为核函数,运用网格搜寻法选择最优的RBF核函数参数组合:核函数的宽度参数,惩罚系数和损失函数;Select RBF as the kernel function, and use the grid search method to select the optimal combination of RBF kernel function parameters: the width parameter of the kernel function, the penalty coefficient and the loss function;

根据最优的RBF核函数参数组合确定最优参数回归模型;Determine the optimal parameter regression model according to the optimal RBF kernel function parameter combination;

将测试集数据带入训练后的回归模型,从而预测出电池恒流充电到达截止电压时所充入的电量Q1Bring the test set data into the trained regression model to predict the charge Q 1 when the constant current charging of the battery reaches the cut-off voltage;

在相同恒流充电条件下,将训练集数据中的恒流充电到截止电压所充入的电量Q11和进行容量测试后得到的当前测试容量Q0运用最小二乘法进行拟合,将测试集数据预测得到的Q1带入拟合后的方程得到当前预测容量cM,从而对电池健康状态进行估计。Under the same constant-current charging condition, use the least squares method to fit the electric quantity Q 11 charged from the constant-current charge in the training set data to the cut-off voltage and the current test capacity Q 0 obtained after the capacity test, and use the least squares method to fit the test set The Q 1 obtained from the data prediction is brought into the fitted equation to obtain the current predicted capacity c M , thereby estimating the state of health of the battery.

本发明有益效果:Beneficial effects of the present invention:

(1)不需要针对特定材料的锂电池,能够对各种恒流充电环境下的锂离子电池恒流充电到达截止电压所充入的电量Q1进行预测,具有广泛的适用性;(1) There is no need for a lithium battery of a specific material, and it is possible to predict the amount of electricity Q1 charged into the lithium-ion battery under various constant-current charging environments when the constant-current charging reaches the cut-off voltage, and has wide applicability;

(2)不需要深入研究锂电池复杂的电化学反应机理,简化了计算过程;(2) There is no need to study the complex electrochemical reaction mechanism of lithium batteries, which simplifies the calculation process;

(3)能够根据充电的数据进行锂电池SOH的估计,不需要单独进行容量测试,缩短了试验周期长并且减少了浪费;(3) It is possible to estimate the SOH of the lithium battery based on the charging data, without the need for a separate capacity test, shortening the long test cycle and reducing waste;

(4)考虑到实际应用,电池不必每次都充电到截止电压。(4) Considering the practical application, the battery does not have to be charged to the cut-off voltage every time.

(5)可确定电池的允许充电量,防止电池恒流充电过程的过充现象的出现,做到对充电设备和锂电池的有效保护。(5) It can determine the allowable charging capacity of the battery, prevent the occurrence of overcharging during the constant current charging process of the battery, and effectively protect the charging equipment and the lithium battery.

附图说明Description of drawings

图1为电池SOH估计的流程图;Figure 1 is a flowchart of battery SOH estimation;

图2为不同SOH值下的锂电池恒流恒压充电图;Figure 2 is a constant current and constant voltage charging diagram of a lithium battery under different SOH values;

图3(a)为训练集恒流充电到达截止电压所充入的电量Q11预测图;Fig. 3 (a) is the electric quantity Q 11 that is charged into the electric quantity Q that the constant current charging of training set reaches the cut-off voltage prediction figure;

图3(b)为测试集恒流充电到达截止电压所充入的电量Q1预测图;Fig. 3(b) is a prediction diagram of the electric quantity Q 1 charged into the test set when the constant current charging reaches the cut-off voltage;

图3(c)为训练集恒流充电的电量Q11和容量测试得到的当前测试容量Q0的拟合关系图。Fig. 3(c) is a fitting relation diagram of the electric quantity Q 11 charged by the constant current in the training set and the current test capacity Q 0 obtained from the capacity test.

具体实施方式Detailed ways

下面结合附图与具体实施方式对本发明作进一步的说明。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.

本发明公开了一种基于支持向量机的锂离子电池健康状态估计方法,如图1所示,包括以下步骤:The invention discloses a method for estimating the state of health of a lithium-ion battery based on a support vector machine, as shown in Figure 1, comprising the following steps:

步骤(1):在对锂电池进行循环充放电的过程中,实时记录锂电池各项工作状态的历史数据,这包括循环次数N、循环充放电电流I和I、放电深度D、温度Tem和下次将要进行的恒流充电的电流大小I,作为支持向量机的回归预测的输入变量。Step (1): During the cycle charge and discharge process of the lithium battery, record the historical data of the various working states of the lithium battery in real time, including the cycle number N, cycle charge and discharge current I charge and I discharge , discharge depth D, temperature T em and the current size I of the constant current charging to be carried out next time are used as input variables for the regression prediction of the support vector machine.

在经历一定的循环周期后,对锂电池进行一次恒流充电和容量测试,统计电池恒流充电到达截止电压所充入的电量Q,作为回归预测的输出变量。After going through a certain cycle period, carry out a constant current charging and capacity test on the lithium battery, and count the amount of electricity Q charged by the constant current charging of the battery to the cut-off voltage, which is used as the output variable of the regression prediction.

步骤(2):将输入和输出变量组分成训练集数据组和测试集数据组,其中训练集用来训练得到回归模型,测试集用来作为回归预测以及对模型的性能做出评价。由于输入和输出是成对存在的,可以将前半部分的输入和输出作为训练集,后半部分的输入和输出作为预测集。Step (2): The input and output variable groups are divided into a training set data set and a test set data set, wherein the training set is used to train the regression model, and the test set is used for regression prediction and evaluation of the performance of the model. Since the input and output exist in pairs, the input and output of the first half can be used as the training set, and the input and output of the second half can be used as the prediction set.

并将数据调整到符合Libsvm工具箱要求的格式。And adjust the data to a format that meets the requirements of the Libsvm toolbox.

步骤(3):对得到的输入和输出变量集进行归一化处理,以加快交叉验证方法求最优解的速度和提高预测的精度。Step (3): Normalize the obtained input and output variable sets to speed up the cross-validation method to find the optimal solution and improve the prediction accuracy.

在建立回归模型之前,先对数据进行归一化处理,即根据公式(1)将输入和输出变量统一归化到[0,1]的取值范围。Before building the regression model, the data is normalized first, that is, the input and output variables are uniformly normalized to the value range of [0,1] according to the formula (1).

其中,xi、yi分别为转换前后的值,MaxValue、MinValue分别为样本的最大值和最小值。Among them, x i and y i are the values before and after conversion respectively, and MaxValue and MinValue are the maximum value and minimum value of the sample respectively.

步骤(4):由上一步得到的归一化后的训练集数据进行回归模型建立,得到回归函数。Step (4): Establish a regression model from the normalized training set data obtained in the previous step to obtain a regression function.

不失一般性,设含有m个训练样本的训练集样本对为其中,Without loss of generality, let the training set sample pair containing m training samples be in,

是第i个训练样本的输入列向量, 为对应的输出值。在该预测的实际应用中,针对此类线性不可分问题,需要通过非线性映射Φ:Rd→H,将原输入空间的样本映射到高维的特征空间H中,再将高维特征空间H中构造最优分类面使得所有训练样本离该最优分类面的误差最小。设在高维特征空间中建立的线性回归函数为公式(2)所示。 is the input column vector of the ith training sample, is the corresponding output value. In the practical application of this prediction, for such linear inseparable problems, it is necessary to map the samples of the original input space to the high-dimensional feature space H through the nonlinear mapping Φ:R d →H, and then map the high-dimensional feature space H The optimal classification surface is constructed in order to minimize the error of all training samples from the optimal classification surface. Assume that the linear regression function established in the high-dimensional feature space is shown in formula (2).

其中,为非线性映射函数,为回归系数向量,为阈值。in, is a nonlinear mapping function, is the regression coefficient vector, is the threshold.

定义ε线性不敏感损失函数为:Define the ε linear insensitive loss function as:

其中,f(x)为回归函数返回的预测值,y为对应的真实值并且ε为事先取定的一个大于0的数。该函数表示若与y之间的差别小于等于ε,则损失等于0。Among them, f(x) is the predicted value returned by the regression function, y is the corresponding real value and ε is a predetermined number greater than 0. This function means that if The difference between y and y is less than or equal to ε, then the loss is equal to 0.

为了寻找到一个最优分类面使得所有训练样本离该最优分类面的误差最小,但考虑到少数样本满足不了约束条件,导致寻找不到最优分类面。针对此类情况,在满足上述约束的条件下使得最小,需要引入松弛变量ξii *,并将上述寻找b的问题用数学语言描述出来,即式(4)所示:In order to find an optimal classification surface, the error of all training samples from the optimal classification surface is minimized, but considering that a small number of samples cannot meet the constraints, the optimal classification surface cannot be found. For such cases, under the conditions of satisfying the above constraints, the minimum, need to introduce slack variables ξ i , ξ i * , and search for The problem of b is described in mathematical language, which is shown in formula (4):

其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大,ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小。Among them, C is the penalty factor. The larger the C, the greater the penalty for the sample whose training error is greater than ε. ε specifies the error requirement of the regression function. The smaller the ε, the smaller the error of the regression function.

对于上述问题的求解是一个凸二次规划问题,由于计算的复杂性,直接求解相对困难。虽然能直接用现成的优化计算包求解,但相比之下,将采用更高效的办法。即依据拉格朗日对偶理论将式(4)转化为对偶问题,这样做的优点在于:一者对偶问题往往更容易求解;二者可以自然的引入核函数,进而推广到非线性分类问题。其中构造的拉格朗日函数如式(5)所示。The solution to the above problem is a convex quadratic programming problem. Due to the complexity of the calculation, it is relatively difficult to solve it directly. Although it can be solved directly with the ready-made optimization calculation package, in comparison, a more efficient method will be adopted. That is, according to the Lagrangian dual theory, formula (4) is transformed into a dual problem. The advantage of this is that: one dual problem is often easier to solve; the two can naturally introduce kernel functions, and then be extended to nonlinear classification problems. The constructed Lagrangian function is shown in formula (5).

其中,a=(a1,a2,...,am)T,a*=(a1 *,a2 *,...,am *)T,η=(η12,...,ηm)T,η*=(η1 *2 *,...,ηm *)T是拉格朗日乘子向量。Among them, a=(a 1 ,a 2 ,...,a m ) T , a * =(a 1 * ,a 2 * ,...,a m * ) T , η=(η 12 ,...,η m ) T , η * = (η 1 *2 * ,...,η m * ) T is a Lagrange multiplier vector.

而求解这个对偶问题,分为3个步骤:首先要让L(ω,b,ξ,ξ*,a,a*,η,η*)关于ω,b和ξ,ξ*最小化,然后求对a,a*,η,η*的极大,最后SMO算法被用来求解对偶问题中的拉格朗日乘子。To solve this dual problem, it is divided into three steps: First, let L(ω,b,ξ,ξ*,a,a*,η,η*) be minimized with respect to ω,b and ξ,ξ*, and then find For the maxima of a, a*, η, η*, finally the SMO algorithm is used to solve the Lagrangian multipliers in the dual problem.

首先固定a,a*,η,η*,要让L关于ω,b和ξ,ξ*最小化,我们分别对ω,b和ξ,ξ*求偏导数,即令等于零,从而得到式(6)First fix a, a*, η, η*, to minimize L with respect to ω, b and ξ, ξ*, we take partial derivatives for ω, b and ξ, ξ* respectively, that is and is equal to zero, so that formula (6)

由拉格朗日对偶性,将式(6)带入到原拉格朗日函数(5)中可将最小值问题转化为最大化问题,如式(7)所示Due to the Lagrangian duality, bringing formula (6) into the original Lagrangian function (5) can transform the minimum value problem into a maximization problem, as shown in formula (7)

其中,为核函数。in, is the kernel function.

设求解式(7)时得到的最优解为则有Suppose the optimal solution obtained when solving formula (7) is then there is

其中,Nnsv为支持向量个数σ。Among them, N nsv is the number σ of support vectors.

于是,回归函数可化为Therefore, the regression function can be transformed into

其中,只有部分参数不等于0,其对应的样本xi即为问题中的支持向量。Among them, only some parameters is not equal to 0, and its corresponding sample xi is the support vector in the problem.

步骤(5):由式(9)可知,要想求解支持向量机算法的回归函数问题,必须要选择一个合适的核函数。其中常用的核函数包括:1.线性核函数;2.d阶多项式核函数;3.径向基核函数(RBF);4.Sigmoid核函数。RBF核函数相比于其它核函数应用相对较广,无论是小样本还是大样本,高维还是低维等情况,RBF核函数均适用。它相比其他的函数有以下优点:Step (5): It can be seen from formula (9) that in order to solve the regression function problem of the support vector machine algorithm, an appropriate kernel function must be selected. The commonly used kernel functions include: 1. Linear kernel function; 2. d-order polynomial kernel function; 3. Radial basis kernel function (RBF); 4. Sigmoid kernel function. Compared with other kernel functions, the RBF kernel function is relatively widely used, whether it is a small sample or a large sample, high-dimensional or low-dimensional, etc., the RBF kernel function is applicable. It has the following advantages over other functions:

1)RBF核函数可以将一个样本映射到一个更高维的空间,而且线性核函数是RBF的一个特例,即如果考虑使用RBF,那么就无需考虑线性核函数。1) The RBF kernel function can map a sample to a higher-dimensional space, and the linear kernel function is a special case of RBF, that is, if RBF is considered, then there is no need to consider the linear kernel function.

2)与多项式核函数相比,RBF需要确定的参数要少,核函数参数的多少直接影响函数的复杂程度。另外,当多项式的阶数比较高时,核矩阵的元素值将趋于无穷大或无穷小。而RBF则在上述问题中,会减少数值的计算困难。2) Compared with the polynomial kernel function, RBF needs to determine fewer parameters, and the number of kernel function parameters directly affects the complexity of the function. In addition, when the order of the polynomial is relatively high, the element values of the kernel matrix will tend to be infinitely large or infinitely small. The RBF will reduce the numerical calculation difficulty in the above problems.

3)对于某些参数,RBF和Sigmoid具有相似的性能。3) For some parameters, RBF and Sigmoid have similar performance.

因此在恒流充电到达截止电压所充入的电量的预测中,选用RBF核函数,其表达式如式(10)所示。Therefore, the RBF kernel function is selected in the prediction of the amount of electricity charged when the constant current charging reaches the cut-off voltage, and its expression is shown in formula (10).

其中,σ为函数的宽度参数,其控制了函数的径向作用范围。Among them, σ is the width parameter of the function, which controls the radial range of the function.

由于核函数模型参数对模型的性能影响较大,因此,在使用Libsvm工具箱时,需要选择较佳的RBF核函数参数组合:核函数的宽度参数σ,惩罚系数C和损失函数P。考虑到需要选择的参数相对较少,“网格搜寻法”的复杂度相比高级算法并无太大差异,而且又具有可并行性高且能够避免过拟合的优势。因此,我们在参数选择中应用“网格搜寻法”。具体实施步骤如下:Since the kernel function model parameters have a great influence on the performance of the model, when using the Libsvm toolbox, it is necessary to select a better combination of RBF kernel function parameters: the width parameter σ of the kernel function, the penalty coefficient C and the loss function P. Considering that there are relatively few parameters to be selected, the complexity of the "grid search method" is not much different from that of advanced algorithms, and it has the advantages of high parallelism and the ability to avoid overfitting. Therefore, we apply the "grid search method" in the parameter selection. The specific implementation steps are as follows:

一、尝试各种可能的组合值,然后进行交叉验证,找出使交叉验证精确度最高的组合对。将各个参数可能的取值进行排列组合,列出所有可能的组合结果生成“网格”。1. Try various possible combination values, and then perform cross-validation to find the combination pair that makes the cross-validation accuracy the highest. Arrange and combine the possible values of each parameter, and list all possible combinations to generate a "grid".

二、然后将各组合用于SVM训练,并使用交叉验证对表现进行评估。Second, each combination is then used for SVM training and performance is evaluated using cross-validation.

三、在拟合函数尝试了所有的参数组合后,返回一个合适的分类面,自动调整至最佳参数组合。3. After the fitting function tries all the parameter combinations, it returns a suitable classification surface and automatically adjusts to the best parameter combination.

步骤(6):利用“网格搜寻法”获得的最优参数训练预测模型。当模型的性能相同时,为了减少计算时间,优先选择惩罚因子C相对较小的参数组合。Step (6): Use the optimal parameters obtained by the "grid search method" to train the prediction model. When the performance of the models is the same, in order to reduce the calculation time, the parameter combination with a relatively small penalty factor C is preferred.

步骤(7):将测试集数据带入训练后的回归模型,从而预测出电池恒流充电到达截止电压所充入的电量Q1而进一步预测SOH值,并对最终的测试集预测结果的准确度做出评价。Step (7): Bring the test set data into the trained regression model, thereby predicting the charge Q1 of the battery constant current charging to the cut-off voltage to further predict the SOH value, and the accuracy of the final test set prediction results degree to make an evaluation.

步骤(8):在相同恒流充电条件下,将训练集数据中的恒流充电到达截止电压所充入的电量Q11和进行容量测试后得到的当前测试容量Q0运用最小二乘法进行拟合,然后将测试集数据预测得到的电量Q1带入拟合后的方程,从而得到预测的当前容量cM,根据做到对电池SOH的估计。Step (8): Under the same constant current charging conditions, use the least squares method to simulate the electric quantity Q 11 that is charged when the constant current charging reaches the cut-off voltage in the training set data and the current test capacity Q 0 obtained after the capacity test. Then put the power Q 1 predicted by the test set data into the fitted equation to get the predicted current capacity c M , according to Make an estimate of the battery SOH.

训练集和测试集恒流充电到达截止电压所充入的电量Q11与Q1预测结果如图3(a)和3(b)所示。训练集相同充电条件下的恒流充电的电量Q11和容量测试得到的当前测试容量Q0的拟合关系如图3(c)所示。Figure 3 (a) and 3 (b). The fitting relationship between the constant current charging power Q 11 and the current test capacity Q 0 obtained from the capacity test under the same charging conditions in the training set is shown in Figure 3(c).

训练集和测试集恒流充电到截止电压所充入电量Q11与Q1预测结果的精度通过均方误差E和决定系数R2评价,训练集和测试集对电池SOH估计的精度通过绝对误差率均值评价,具体结果如表1所示。从具体评价指标来说,基于支持向量机的恒流充电到截止电压所充入的电量Q1预测及SOH估计方法具有较高的预测精度。The accuracy of the prediction results of Q11 and Q1 from the constant current charging to the cut-off voltage of the training set and the test set is evaluated by the mean square error E and the coefficient of determination R2. The accuracy of the training set and the test set for the battery SOH estimation is evaluated by the absolute error The specific results are shown in Table 1. From the specific evaluation index, based on the support vector machine constant current charging to the cut-off voltage charging Q1 prediction and SOH estimation method has a higher prediction accuracy.

表1锂电池恒流充电到截止电压所充入的电量Q预测及SOH估计方法精度评价Table 1 Prediction of the amount of electricity Q charged from constant current charging to the cut-off voltage of lithium batteries and the accuracy evaluation of SOH estimation methods

EE. R2 R 2 SOH估计绝对误差率均值SOH estimated absolute error rate mean 训练集Training set 0.000655660.00065566 0.992710.99271 0.48%0.48% 测试集test set 0.000860360.00086036 0.991120.99112 0.56%0.56%

由实验可知,相同充电条件下,锂电池恒流充电到达截止电压所充入的电量Q随着SOH值的减小而缩短且两者存在较高的关联度。因此根据之前测试的相关数据,可由最小二乘法拟合Q11与SOH值两者之间的关系,从而由上述支持向量机算法预测的恒流充电电量Q1对当前锂电池的SOH做出估计。It can be seen from the experiment that under the same charging conditions, the charge Q of lithium battery constant current charging to reach the cut-off voltage is shortened as the SOH value decreases, and there is a high degree of correlation between the two. Therefore, according to the relevant data of the previous test, the relationship between Q11 and SOH value can be fitted by the least square method, so that the constant current charging quantity Q1 predicted by the above-mentioned support vector machine algorithm can be used to estimate the SOH of the current lithium battery .

图2给出了不同SOH值下的锂电池恒流恒压充电图,由图2可知,在相同恒流充电条件下,恒流充电的电量随SOH值的降低呈下降趋势。Figure 2 shows the constant current and constant voltage charging diagram of lithium batteries under different SOH values. It can be seen from Figure 2 that under the same constant current charging conditions, the constant current charging power decreases with the decrease of SOH value.

上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。Although the specific implementation of the present invention has been described above in conjunction with the accompanying drawings, it does not limit the protection scope of the present invention. Those skilled in the art should understand that on the basis of the technical solution of the present invention, those skilled in the art do not need to pay creative work Various modifications or variations that can be made are still within the protection scope of the present invention.

Claims (10)

1. A lithium ion battery health state estimation method based on a support vector machine is characterized by comprising the following steps:
carrying out a cyclic charge and discharge experiment on the lithium battery, recording historical data of various working states of the lithium battery in real time, and determining an input variable and an output variable of regression prediction of a support vector machine;
the input variables and the output variables are divided into two groups: a training set data set and a test set data set;
carrying out normalization processing on the training set data group and the test set data group;
carrying out regression model establishment on the normalized training set data to obtain a regression function;
selecting RBF as kernel function, selecting optimal RBF kernel function parameter combination by using grid search method: width parameters of kernel functions, penalty coefficients and loss functions;
determining an optimal parameter regression model according to the optimal RBF kernel function parameter combination;
the data of the test set is substituted into the regression model after training, so that the electric quantity Q charged when the constant current charging of the battery reaches the cut-off voltage is predicted1
Under the same constant current charging condition, the constant current in the training set data is charged to the electric quantity Q charged by the cut-off voltage11And the current testing capacity Q obtained after the capacity testing0Performing fitting by using a least square method, and predicting the Q obtained by the data of the test set1The fitted equation is substituted to obtain the current predicted capacity cMAnd thus the state of health of the battery is estimated.
2. The method for estimating the state of health of a lithium ion battery based on a support vector machine according to claim 1, wherein the determining of the input variables and the output variables of the regression prediction of the support vector machine is specifically:
recording historical data of various working states of the lithium battery in real time, wherein the historical data comprises cycle times, cycle charging and discharging current, discharging depth, temperature and current of constant current charging to be carried out next time, and the historical data is used as an input variable for regression prediction of a support vector machine;
after a set cycle period, performing constant current charging and capacity testing on the lithium battery once, and counting the electric quantity Q charged when the constant current charging of the battery reaches a cut-off voltage, wherein the electric quantity Q is used as an output variable for regression prediction of the support vector machine.
3. The lithium ion battery health state estimation method based on the support vector machine according to claim 1, wherein the normalization processing is performed on the training set data group and the test set data group, specifically:
the normalized sample value is the ratio of the difference between the sample value and the minimum value of the sample to the difference between the maximum value and the minimum value of the sample value.
4. The lithium ion battery state of health estimation method based on support vector machine of claim 1, characterized in that, when the performance of the parameter training prediction model is the same, the parameter combination with relatively small penalty coefficient is selected preferentially.
5. The lithium ion battery state of health estimation method based on support vector machine of claim 1, characterized in that, the normalized training set data is subjected to regression model establishment to obtain a regression function; the method specifically comprises the following steps:
let the training set sample pair containing m training samples beWherein,is the input column vector for the ith training sample, is the corresponding output value;
the linear regression function established in the high-dimensional feature space is set as:
define the epsilon linear insensitive loss function as:
in order to find an optimal classification surface to minimize the error of all training samples from the optimal classification surface, a relaxation variable ξ is introducedii *Then, the following constraint relation needs to be satisfied:
wherein,in order to be a non-linear mapping function,is a regression coefficient vector, b is a threshold;a predicted value returned by the regression function is used, y is a corresponding true value, and epsilon is a set number larger than 0; and C is a penalty factor, wherein the larger C is the penalty larger for the sample with the training error larger than epsilon, epsilon specifies the error requirement of the regression function, and the smaller epsilon is the error of the regression function.
6. The lithium ion battery state of health estimation method based on a support vector machine according to claim 5, characterized in that the constraint relation is solved, specifically:
converting the constraint relation into a dual problem by adopting a Lagrange dual theory;
converting the minimum problem in the original Lagrange function into a maximization problem by Lagrange duality;
solving a Lagrange multiplier in the dual problem by adopting an SMO algorithm; further obtain the regression coefficient vectorAnd a threshold value b, determining a regression function.
7. The lithium ion battery state of health estimation method based on a support vector machine according to claim 6, characterized in that the regression function is specifically:
wherein,respectively the optimal solutions of the lagrange multipliers corresponding to sample i,is a kernel function, b*The obtained threshold value is solved according to the optimal solution of the Lagrange multiplier.
8. The method according to claim 1, wherein the RBFs are selected as kernel functions, and an optimal RBF kernel function parameter combination is selected by using a grid search method; the method specifically comprises the following steps:
the method comprises the following steps: performing cross validation on various possible RBF kernel function parameter combination values, and finding out a combination pair which enables the cross validation accuracy to be highest; arranging and combining the possible values of each parameter, and listing all possible combination results to generate a 'grid';
step two: using each RBF kernel function parameter combination for SVM training, and evaluating the performance by using cross validation;
step three: and after the fitting function tries all parameter combinations, returning to an optimal classification surface, and automatically adjusting to the optimal parameter combination to ensure that the error of all training samples from the optimal classification surface is minimum.
9. A system for estimating state of health of a lithium ion battery based on a support vector machine, comprising a server including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the program:
determining an input variable and an output variable of the regression prediction of the support vector machine;
the input variables and the output variables are divided into two groups: a training set data set and a test set data set;
carrying out normalization processing on the training set data group and the test set data group;
carrying out regression model establishment on the normalized training set data to obtain a regression function;
selecting RBF as kernel function, selecting optimal RBF kernel function parameter combination by using grid search method: width parameters of kernel functions, penalty coefficients and loss functions;
determining an optimal parameter regression model according to the optimal RBF kernel function parameter combination;
the data of the test set is substituted into the regression model after training, so that the electric quantity Q charged when the constant current charging of the battery reaches the cut-off voltage is predicted1
Under the same constant current charging condition, the constant current in the training set data is charged to the electric quantity Q charged by the cut-off voltage11And the current testing capacity Q obtained after the capacity testing0Performing fitting by using a least square method, and predicting the Q obtained by the data of the test set1The fitted equation is substituted to obtain the current predicted capacity cMAnd thus the state of health of the battery is estimated.
10. A computer-readable storage medium on which a computer program is stored, the program, when executed by a processor, performing the steps of:
determining an input variable and an output variable of the regression prediction of the support vector machine;
the input variables and the output variables are divided into two groups: a training set data set and a test set data set;
carrying out normalization processing on the training set data group and the test set data group;
carrying out regression model establishment on the normalized training set data to obtain a regression function;
selecting RBF as kernel function, selecting optimal RBF kernel function parameter combination by using grid search method: width parameters of kernel functions, penalty coefficients and loss functions;
determining an optimal parameter regression model according to the optimal RBF kernel function parameter combination;
the data of the test set is substituted into the regression model after training, so that the electric quantity Q charged when the constant current charging of the battery reaches the cut-off voltage is predicted1
Under the same constant current charging condition, the constant current in the training set data is charged to the electric quantity Q charged by the cut-off voltage11And the current testing capacity Q obtained after the capacity testing0Performing fitting by using a least square method, and predicting the Q obtained by the data of the test set1The fitted equation is substituted to obtain the current predicted capacity cMAnd thus the state of health of the battery is estimated.
CN201810636667.9A 2018-06-20 2018-06-20 A method and system for estimating state of health of lithium-ion battery based on support vector machine Active CN108805217B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810636667.9A CN108805217B (en) 2018-06-20 2018-06-20 A method and system for estimating state of health of lithium-ion battery based on support vector machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810636667.9A CN108805217B (en) 2018-06-20 2018-06-20 A method and system for estimating state of health of lithium-ion battery based on support vector machine

Publications (2)

Publication Number Publication Date
CN108805217A true CN108805217A (en) 2018-11-13
CN108805217B CN108805217B (en) 2020-10-23

Family

ID=64083726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810636667.9A Active CN108805217B (en) 2018-06-20 2018-06-20 A method and system for estimating state of health of lithium-ion battery based on support vector machine

Country Status (1)

Country Link
CN (1) CN108805217B (en)

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109543759A (en) * 2018-11-27 2019-03-29 北京石油化工学院 A kind of prediction technique of single flow gas-liquid cyclone separator separating property
CN109927575A (en) * 2019-02-28 2019-06-25 福建工程学院 A kind of battery performance detection method for direct-current charging post
CN110080882A (en) * 2019-04-16 2019-08-02 新奥能源动力科技(上海)有限公司 A kind of the starting method and starter of gas turbine
CN110275129A (en) * 2019-06-18 2019-09-24 中国电力科学研究院有限公司 A method and system for determining the synthesis error of a high-voltage electric energy metering device
CN110398697A (en) * 2019-07-23 2019-11-01 北京工业大学 A method for estimating the state of health of lithium ions based on the charging process
CN110501646A (en) * 2019-08-29 2019-11-26 中国人民解放军国防科技大学 Off-line lithium battery residual capacity estimation method
CN110568374A (en) * 2019-09-02 2019-12-13 东北电力大学 Prediction method of remaining service life of lithium-ion battery based on internal and external characteristics
CN110608660A (en) * 2019-08-28 2019-12-24 西安理工大学 A Displacement and Temperature Prediction Method of Eddy Current Sensor
CN110850315A (en) * 2019-11-29 2020-02-28 北京邮电大学 Method and device for estimating battery state of charge
CN111274539A (en) * 2020-02-18 2020-06-12 金陵科技学院 A Lithium Battery SOH Estimation Method Based on Alternating Least Squares
CN111308375A (en) * 2020-02-04 2020-06-19 浙江大学 LSTM-FFNN-based electric forklift lithium ion battery health state prediction method
CN111323705A (en) * 2020-03-19 2020-06-23 山东大学 Battery parameter identification method and system based on robust recursive least squares
CN111443293A (en) * 2020-03-30 2020-07-24 青岛大学 Lithium battery state of health (SOH) estimation method based on data driving
CN111999656A (en) * 2020-08-28 2020-11-27 广州小鹏汽车科技有限公司 Method and device for detecting short circuit in vehicle battery and electronic equipment
CN112213643A (en) * 2020-09-30 2021-01-12 蜂巢能源科技有限公司 Method, system and equipment for predicting initial capacity and health state of battery
CN112255559A (en) * 2020-10-12 2021-01-22 江苏慧智能源工程技术创新研究院有限公司 Method for predicting residual life of lithium battery energy storage power station based on multiple linear regression
CN112287597A (en) * 2020-09-22 2021-01-29 国网天津市电力公司电力科学研究院 Lead-acid storage battery SOH estimation method based on VPGA-GPR algorithm
CN112347692A (en) * 2020-09-21 2021-02-09 深圳前海有电物联科技有限公司 Method and device for realizing predictive maintenance of battery of uninterruptible power supply and electronic device
CN112540318A (en) * 2020-12-22 2021-03-23 武汉理工大学 Method for estimating health state of lead-acid storage battery for starting internal combustion engine
CN112580211A (en) * 2020-12-23 2021-03-30 天津大学 Lead-acid storage battery SOH estimation method based on SA and ANN algorithm
CN112684356A (en) * 2020-10-31 2021-04-20 浙江锋锂新能源科技有限公司 Cycle test method of lithium ion battery
CN112731184A (en) * 2020-12-28 2021-04-30 深圳供电局有限公司 Battery service life detection method and system
CN112834927A (en) * 2021-01-06 2021-05-25 合肥工业大学 Method, system, device and medium for predicting remaining life of lithium battery
CN112924886A (en) * 2021-01-23 2021-06-08 青岛大学 Battery state of health (SOH) prediction method and device
CN113011464A (en) * 2021-02-25 2021-06-22 沈阳工业大学 Comprehensive prediction method for running state of transformer based on multi-dimensional data evaluation
CN113219341A (en) * 2021-03-23 2021-08-06 陈九廷 Model generation and battery degradation estimation device, method, medium, and apparatus
CN113267733A (en) * 2021-04-13 2021-08-17 西安理工大学 Automatic configuration method for lithium battery health state estimation based on Gaussian process regression
CN113866644A (en) * 2021-09-30 2021-12-31 国网福建省电力有限公司龙岩供电公司 Method and device for predicting usable time and capacity of battery
CN113935225A (en) * 2020-06-29 2022-01-14 中国科学院大连化学物理研究所 A machine learning-based method for optimization and performance prediction of flow battery stacks
CN114355222A (en) * 2021-12-23 2022-04-15 厦门大学 Battery state of health estimation method, device and readable medium based on voltage curve
CN114895206A (en) * 2022-04-26 2022-08-12 合肥工业大学 SOH estimation method of lithium-ion battery based on RBF neural network based on improved gray wolf optimization algorithm
CN115267586A (en) * 2022-07-11 2022-11-01 国网综合能源服务集团有限公司 Lithium battery SOH evaluation method
CN115656855A (en) * 2022-03-15 2023-01-31 上海舞洋船舶科技有限公司 Lithium ion battery health assessment method and system
CN116413609A (en) * 2023-06-08 2023-07-11 江苏正力新能电池技术有限公司 Battery diving identification method and device, electronic equipment and storage medium
CN118275903A (en) * 2024-06-04 2024-07-02 河南师范大学 Battery performance test method based on data analysis
CN118348812A (en) * 2024-06-18 2024-07-16 浙江维度仪表有限公司 Intelligent control valve regulation and control method and system based on Internet of things
CN118625156A (en) * 2024-08-14 2024-09-10 长安绿电科技有限公司 A lithium battery status monitoring method, device, medium and equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102680907A (en) * 2012-05-31 2012-09-19 浙江大学 Battery charging stress optical coefficient (SOC) detection method in photovoltaic system
CN102749585A (en) * 2011-04-21 2012-10-24 李昌 Multi-layer SVM (support vector machine) based storage battery on-line monitoring method
CN104021238A (en) * 2014-03-25 2014-09-03 重庆邮电大学 Lead-acid power battery system fault diagnosis method
CN104391252A (en) * 2014-12-04 2015-03-04 上海理工大学 Automobile lead-acid battery health state detection method
EP2975421A1 (en) * 2014-07-18 2016-01-20 Samsung Electronics Co., Ltd Method and apparatus for estimating state of battery
CN105334465A (en) * 2015-09-15 2016-02-17 重庆长安汽车股份有限公司 Method for online evaluating state of health of lithium ion battery
CN107132490A (en) * 2017-07-05 2017-09-05 福州大学 A kind of method for realizing the estimation of lithium battery group state-of-charge
CN107290679A (en) * 2017-07-03 2017-10-24 南京能瑞电力科技有限公司 The Intelligentized battery method for detecting health status of charging pile is shared for electric automobile

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749585A (en) * 2011-04-21 2012-10-24 李昌 Multi-layer SVM (support vector machine) based storage battery on-line monitoring method
CN102680907A (en) * 2012-05-31 2012-09-19 浙江大学 Battery charging stress optical coefficient (SOC) detection method in photovoltaic system
CN104021238A (en) * 2014-03-25 2014-09-03 重庆邮电大学 Lead-acid power battery system fault diagnosis method
EP2975421A1 (en) * 2014-07-18 2016-01-20 Samsung Electronics Co., Ltd Method and apparatus for estimating state of battery
CN104391252A (en) * 2014-12-04 2015-03-04 上海理工大学 Automobile lead-acid battery health state detection method
CN105334465A (en) * 2015-09-15 2016-02-17 重庆长安汽车股份有限公司 Method for online evaluating state of health of lithium ion battery
CN107290679A (en) * 2017-07-03 2017-10-24 南京能瑞电力科技有限公司 The Intelligentized battery method for detecting health status of charging pile is shared for electric automobile
CN107132490A (en) * 2017-07-05 2017-09-05 福州大学 A kind of method for realizing the estimation of lithium battery group state-of-charge

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HUANG HAI ET AL: "Aging Performances and Cycle-life Predictions of Li-ion Battery", 《PROCEEDINGS OF THE 35TH CHINESE CONTROL CONFERENCE》 *
WENWEI HUANG ET AL: "SOC Prediction of Lithium Battery Based on Fuzzy Information Granulation and Support Vector Regression", 《 2017 4TH INTERNATIONAL CONFERENCE ON ELECTRICAL AND ELECTRONIC ENGINEERING (ICEEE)》 *
何旭等: "基于SVM的小样本条件下继电保护可靠性参数估计", 《电网技术》 *
李睿琪等: "一种基于支持向量机的锂电池健康状态评估方法", 《第17届中国系统仿真技术及其应用学术年会(17TH CCSSTA 2016)》 *

Cited By (56)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109543759A (en) * 2018-11-27 2019-03-29 北京石油化工学院 A kind of prediction technique of single flow gas-liquid cyclone separator separating property
CN109927575A (en) * 2019-02-28 2019-06-25 福建工程学院 A kind of battery performance detection method for direct-current charging post
CN110080882A (en) * 2019-04-16 2019-08-02 新奥能源动力科技(上海)有限公司 A kind of the starting method and starter of gas turbine
CN110275129B (en) * 2019-06-18 2022-10-04 中国电力科学研究院有限公司 Method and system for determining synthetic error of high-voltage electric energy metering device
CN110275129A (en) * 2019-06-18 2019-09-24 中国电力科学研究院有限公司 A method and system for determining the synthesis error of a high-voltage electric energy metering device
CN110398697A (en) * 2019-07-23 2019-11-01 北京工业大学 A method for estimating the state of health of lithium ions based on the charging process
CN110608660A (en) * 2019-08-28 2019-12-24 西安理工大学 A Displacement and Temperature Prediction Method of Eddy Current Sensor
CN110608660B (en) * 2019-08-28 2021-11-16 西安理工大学 Eddy current sensor displacement and temperature prediction method
CN110501646A (en) * 2019-08-29 2019-11-26 中国人民解放军国防科技大学 Off-line lithium battery residual capacity estimation method
CN110568374A (en) * 2019-09-02 2019-12-13 东北电力大学 Prediction method of remaining service life of lithium-ion battery based on internal and external characteristics
CN110568374B (en) * 2019-09-02 2021-04-27 东北电力大学 Prediction method of remaining service life of lithium-ion battery based on internal and external characteristics
CN110850315A (en) * 2019-11-29 2020-02-28 北京邮电大学 Method and device for estimating battery state of charge
CN111308375A (en) * 2020-02-04 2020-06-19 浙江大学 LSTM-FFNN-based electric forklift lithium ion battery health state prediction method
CN111274539A (en) * 2020-02-18 2020-06-12 金陵科技学院 A Lithium Battery SOH Estimation Method Based on Alternating Least Squares
CN111323705B (en) * 2020-03-19 2021-07-23 山东大学 Battery parameter identification method and system based on robust recursive least squares
CN111323705A (en) * 2020-03-19 2020-06-23 山东大学 Battery parameter identification method and system based on robust recursive least squares
CN111443293A (en) * 2020-03-30 2020-07-24 青岛大学 Lithium battery state of health (SOH) estimation method based on data driving
CN113935225A (en) * 2020-06-29 2022-01-14 中国科学院大连化学物理研究所 A machine learning-based method for optimization and performance prediction of flow battery stacks
CN113935225B (en) * 2020-06-29 2024-05-07 中国科学院大连化学物理研究所 Flow battery pile optimization and performance prediction method based on machine learning
CN111999656A (en) * 2020-08-28 2020-11-27 广州小鹏汽车科技有限公司 Method and device for detecting short circuit in vehicle battery and electronic equipment
CN112347692B (en) * 2020-09-21 2024-02-06 深圳有电物联科技有限公司 Method and device for realizing battery predictive maintenance of uninterruptible power supply and electronic device
CN112347692A (en) * 2020-09-21 2021-02-09 深圳前海有电物联科技有限公司 Method and device for realizing predictive maintenance of battery of uninterruptible power supply and electronic device
CN112287597B (en) * 2020-09-22 2023-10-03 国网天津市电力公司电力科学研究院 Lead-acid storage battery SOH estimation method based on VPGA-GPR algorithm
CN112287597A (en) * 2020-09-22 2021-01-29 国网天津市电力公司电力科学研究院 Lead-acid storage battery SOH estimation method based on VPGA-GPR algorithm
CN112213643B (en) * 2020-09-30 2023-06-23 蜂巢能源科技有限公司 Prediction method, system and equipment of battery initial capacity and battery health state
CN112213643A (en) * 2020-09-30 2021-01-12 蜂巢能源科技有限公司 Method, system and equipment for predicting initial capacity and health state of battery
CN112255559A (en) * 2020-10-12 2021-01-22 江苏慧智能源工程技术创新研究院有限公司 Method for predicting residual life of lithium battery energy storage power station based on multiple linear regression
CN112255559B (en) * 2020-10-12 2023-05-09 启东沃太新能源有限公司 Lithium battery energy storage power station residual life prediction method based on multiple linear regression
CN112684356B (en) * 2020-10-31 2024-06-04 浙江锋锂新能源科技有限公司 Circulation test method of lithium ion battery
CN112684356A (en) * 2020-10-31 2021-04-20 浙江锋锂新能源科技有限公司 Cycle test method of lithium ion battery
CN112540318A (en) * 2020-12-22 2021-03-23 武汉理工大学 Method for estimating health state of lead-acid storage battery for starting internal combustion engine
CN112540318B (en) * 2020-12-22 2022-01-04 武汉理工大学 A method for estimating the state of health of a lead-acid battery for starting an internal combustion engine
CN112580211A (en) * 2020-12-23 2021-03-30 天津大学 Lead-acid storage battery SOH estimation method based on SA and ANN algorithm
CN112731184A (en) * 2020-12-28 2021-04-30 深圳供电局有限公司 Battery service life detection method and system
CN112731184B (en) * 2020-12-28 2023-03-03 深圳供电局有限公司 Battery service life detection method and system
CN112834927A (en) * 2021-01-06 2021-05-25 合肥工业大学 Method, system, device and medium for predicting remaining life of lithium battery
CN112924886A (en) * 2021-01-23 2021-06-08 青岛大学 Battery state of health (SOH) prediction method and device
CN113011464A (en) * 2021-02-25 2021-06-22 沈阳工业大学 Comprehensive prediction method for running state of transformer based on multi-dimensional data evaluation
CN113219341A (en) * 2021-03-23 2021-08-06 陈九廷 Model generation and battery degradation estimation device, method, medium, and apparatus
CN113267733B (en) * 2021-04-13 2023-11-17 西安理工大学 Automatic configuration method for lithium battery health state estimation based on Gaussian process regression
CN113267733A (en) * 2021-04-13 2021-08-17 西安理工大学 Automatic configuration method for lithium battery health state estimation based on Gaussian process regression
CN113866644B (en) * 2021-09-30 2024-12-20 国网福建省电力有限公司龙岩供电公司 A method and device for predicting battery service life and capacity
CN113866644A (en) * 2021-09-30 2021-12-31 国网福建省电力有限公司龙岩供电公司 Method and device for predicting usable time and capacity of battery
CN114355222B (en) * 2021-12-23 2024-10-15 厦门大学 Method and device for estimating battery state of health based on voltage curve and readable medium
CN114355222A (en) * 2021-12-23 2022-04-15 厦门大学 Battery state of health estimation method, device and readable medium based on voltage curve
CN115656855A (en) * 2022-03-15 2023-01-31 上海舞洋船舶科技有限公司 Lithium ion battery health assessment method and system
CN114895206B (en) * 2022-04-26 2023-04-28 合肥工业大学 Lithium-ion battery SOH estimation method based on improved gray wolf optimization algorithm and RBF neural network
CN114895206A (en) * 2022-04-26 2022-08-12 合肥工业大学 SOH estimation method of lithium-ion battery based on RBF neural network based on improved gray wolf optimization algorithm
CN115267586A (en) * 2022-07-11 2022-11-01 国网综合能源服务集团有限公司 Lithium battery SOH evaluation method
CN116413609B (en) * 2023-06-08 2023-08-29 江苏正力新能电池技术有限公司 Battery diving identification method and device, electronic equipment and storage medium
CN116413609A (en) * 2023-06-08 2023-07-11 江苏正力新能电池技术有限公司 Battery diving identification method and device, electronic equipment and storage medium
CN118275903B (en) * 2024-06-04 2024-09-03 河南师范大学 Battery performance test method based on data analysis
CN118275903A (en) * 2024-06-04 2024-07-02 河南师范大学 Battery performance test method based on data analysis
CN118348812A (en) * 2024-06-18 2024-07-16 浙江维度仪表有限公司 Intelligent control valve regulation and control method and system based on Internet of things
CN118348812B (en) * 2024-06-18 2024-08-20 浙江维度仪表有限公司 Intelligent control valve regulation and control method and system based on Internet of things
CN118625156A (en) * 2024-08-14 2024-09-10 长安绿电科技有限公司 A lithium battery status monitoring method, device, medium and equipment

Also Published As

Publication number Publication date
CN108805217B (en) 2020-10-23

Similar Documents

Publication Publication Date Title
CN108805217B (en) A method and system for estimating state of health of lithium-ion battery based on support vector machine
Cui et al. A dynamic spatial-temporal attention-based GRU model with healthy features for state-of-health estimation of lithium-ion batteries
KR102650965B1 (en) Method of estimating battery states
WO2022198616A1 (en) Battery life prediction method and system, electronic device, and storage medium
Chen et al. State of health estimation for lithium-ion battery based on particle swarm optimization algorithm and extreme learning machine
CN114861527A (en) Lithium battery life prediction method based on time series characteristics
CN105911476B (en) A kind of battery energy storage system SOC prediction techniques based on data mining
CN108931729B (en) A dynamic identification method for capacity cycling degradation of lithium-ion batteries
CN114636932A (en) Method and system for predicting remaining service life of battery
CN111950205A (en) A Lithium Battery SOH Prediction Method Based on FWA Optimization Extreme Learning Machine
Li et al. A novel hybrid data-driven method based on uncertainty quantification to predict the remaining useful life of lithium battery
CN114781176B (en) An equivalent circuit parameter identification method for lumped parameters of lithium-ion battery energy storage system
CN115308608A (en) All-vanadium redox flow battery voltage prediction method, device and medium
Yüksek et al. A novel state of health estimation approach based on polynomial model for lithium-ion batteries
CN115270454A (en) Battery life prediction method and related equipment
Upashrutti et al. Estimation of State of Charge of EV Batteries-A Machine Learning Approach
CN114089204A (en) Battery capacity diving inflection point prediction method and device
CN118425815A (en) Energy storage aging state estimation method based on transposed attention model
Gao et al. A novel indirect health indicator extraction based on charging data for lithium-ion batteries remaining useful life prognostics
CN116774045A (en) A lithium battery health status prediction method based on HHO-SVR
CN116736171A (en) A data-driven method for estimating health status of lithium-ion batteries
Hua et al. Surrogate modelling for battery state-of-charge estimation in electric vehicles based on pseudo-2-dimensional model and gradient boosting machines
Lipu et al. Optimal neural network approach for estimating state of energy of lithium-ion battery using heuristic optimization techniques
Chen et al. Energy storage battery state of health estimation based on singular value decomposition for noise reduction and improved LSTM neural network
Zhu et al. State of Health Estimation of Lithium Ion Battery Based on CNN-LSTM Neural Network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant