本发明涉及磷酸铁锂动力电池的技术领域，更具体地说，涉及一种基于电池滞回电压特性和回弹电压特性的扩展卡尔曼滤波soc估算方法。

在电动汽车用锂动力电池管理系统中，电池荷电状态(soc)的估算是电池管理系统(bms)的核心内容，且电池的均衡控制，健康状态(soh)控制均是以soc的估算为基础，soc的估算准确与否，将直接影响到电池管理系统的决策，因此精确的soc估计对改善电池性能，提高电池可靠性，延长电池使用寿命有重要意义。

电池soc的估算与电池的开路电压、充放电电流、蓄电池内阻。自放电及电池的循环寿命等参数有关。目前国内外常用的soc估算方法为开路电压法和安时积分法。开路电压法为电池停止工作测量电路开路电压值，通过ocv-soc查表估计此时的soc值，该方法需要静置足够的时间，无法实现在线估算；安时积分法通过电流的动态积分来估算soc，算法简单，但其为一种开环预测，无法实时修正，且无法确定初值，同时存在传感器的累积误差。

扩展卡尔曼滤波算法基于建立的精确的电池等效电路模型建立电池的状态空间方程，结合递推算法实现对soc的估计，可实现在线估计，同时可给出误差估计范围，满足电动汽车实际工况下的电流剧烈变化条件。

本发明的主要目的在于克服现有技术的缺点与不足，一种基于电池滞回电压和回弹电压特性的扩展卡尔曼滤波soc估算方法，克服了开路电压法无法在线估算，安时积分法电流误差累积的缺点，提高了电池soc的估算精度，实现了电池soc的在线估算。

为了达到上述目的，本发明采用以下技术方案：

本发明公开的一种基于电池滞回电压特性和回弹电压特性的扩展卡尔曼滤波soc估算方法，包括下述步骤：

(1)采集当前状态下的电压、电流、历史soc、累计充放电次数n这些参数信息；

(2)基于电池滞回电压特性和回弹电压特性建立电池的三阶rc等效电路模型；

(3)基于改进安时积分法建立soc状态方程；

(4)基于建立的状态方程和观测方程，根据扩展卡尔曼递推原理不断对当前时刻soc值进行修正，使估算值逐渐逼近真实值；

(5)soc误差小于2％时结束soc的扩展卡尔曼递推过程。

作为优选的技术方案，步骤(2)中，建立三阶rc等效电路模型的方法为：

建立三阶rc等效电路模型描述其响应特性和内部特性，响应特性是指电池的端电压与负载电流的对应关系；内部特性是指电池的内部变量欧姆内阻、极化内阻、极化电压与soc的关系，建立的观测方程如下：

uc(t)为电池的负载电压，eb为电池的等效电压，ur为电池欧姆内阻对应的电压，us、um、ul为三阶rc网络分别对应的电压，uemf为反电动势，与充放电soc有关，uvh为电池滞回电压，与电流i方向，电池soc有关；

由实验可知，电池内阻rr受soc影响较大，0≤soc≤50％时，rr随soc变化较大，50％＜soc≤100时，rr趋于稳定，因此，rr与soc关系描述如下：

由电化学原理可知，可采用rc网络描述电池的回弹特性，其精度与rc的阶数有关，采用三阶rc网络来表达电池的回弹电压特性，既满足精度要求，又不会增加计算的复杂度，表示如下：

作为优选的技术方案，所述步骤(3)中，基于改进安时积分法建立soc状态方程的方法为：

其中，α为电池的老化系数，与充放电次数n有关，soc(0)为电池的初始时刻soc，cn为电池的标称容量，单位为ah，表示电量完全充满的情况下电池以实际的工作温度和放电倍率下所放出的电量，η为电池效率系数，主要包括充放电系数ηc和温度系数ηt，它们之间满足η＝ηc·ηt的函数关系；

其中，soch为上次关机时系统记录的历史soc估计值，通过标准充放电实验分段拟合的平衡电势与soc一阶函数关系，sococv为通过采集电池开机电压、电流计算的soc值，β为与上次关机到此次开机的时间间隔t的相关系数。

作为优选的技术方案，在步骤3.1)中，α满足如下函数关系：

在步骤3.2)中，β满足如下函数关系：

作为优选的技术方案，步骤3.2)中在计算初始soc时，使用matlab进行一阶拟合，拟合方程为：

作为优选的技术方案，步骤(4)中，建立的状态方程和观测方程，根据扩展卡尔曼递推原理不断对当前时刻soc值进行修正的方法为：

4.1)获取扩展卡尔曼滤波状态估计向量初始值x0与误差协方差矩阵的初始值p0；

4.2)将k-1时刻的状态向量估计值与误差协方差矩阵pk-1/k-1进行时间更新进一步确定k时刻的状态向量先验估计值与误差协方差pk/k-1:

其中wk为均值为0的高斯白噪声，为wk的转置，e为数学期望

4.3)确定扩展卡尔曼增益矩阵kk

其中ck为雅可比矩阵，为ck的转置，rk为观测噪声协方差

vk为均值为0的高斯白噪声

4.4)获得卡尔曼滤波状态向量估计值及状态协方差

本发明与现有技术相比，具有如下优点和有益效果：

1、本发明采用改进的安时积分法和扩展卡尔曼滤波算法计算电池soc值，克服了安时积分在线估算时间长的问题。

2、本发明基于电池电压滞回特性和回弹特性建立较精确的三阶rc电池模型，提高了soc模型的估算精度。

3、本发明采用电池关机记录soc值与ocv-soc查表结合的方法计算电池soc初值，同时考虑到电池的老化效应，提高了soc初值的估算精度。

图2为本发明电池等效电路模型图；

图3为本发明内阻rs与soc关系图；

图4为本发明电池滞回电压特性图；

图5为本发明回弹电压特性图；

图6为本发明计算初始soc流程图。

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

如图1所示，本发明soc估算方法主要包括电池等效电路建模，构建扩展卡尔曼状态方程，扩展卡尔曼滤波递推算法状态变量更新计算几部分。

s1、采集当前状态下的电压、电流、历史soc、累计充放电次数n等参数信息；

s2、基于电池滞回电压特性和回弹电压特性建立电池的三阶rc等效电路模型；

图2所示等效电路模型，应用扩展卡尔曼滤波器估算电池系统soc，需要建立精确的电池模型。电池模型应尽量完善表达实际电池的响应特性和内部特性。响应特性是指电池的端电压与负载电流的对应关系；内部特性是指电池的内部变量欧姆内阻、极化内阻、极化电压与soc的关系，建立的观测方程如下：

图3为电池内阻rr与soc关系图，由图可知，电池欧姆内阻rr受soc影响较大，0≤soc≤50％时，rr随soc变化较大，50％＜soc≤100时，rr趋于稳定。因此，rr与soc关系描述如下：

图4为电池回弹电压特性图，由图可知，当回路电流为0时，电池内部由于其极化反应并未停止，因此其净反应速率不会立即降为0，电池端电压会继续上升，直至经历足够的静置时间，端电压才会恢复稳定，即与开路电压相同。由电化学原理可知，可采用rc网络描述电池的回弹特性，其精度与rc的阶数有关。本发明采用三阶rc网络来表达电池的回弹电压特性，既满足精度要求，又不会增加计算的复杂度。

由基尔霍夫电压电流定律可知三阶rc网络满足如下条件：

s3、基于改进安时积分法建立soc状态方程：

其中，α为电池的老化系数，与充放电次数n有关，soc(0)为电池的初始时刻soc，cn为电池的标称容量，单位为ah，表示电量完全充满的情况下电池以实际的工作温度和放电倍率下所放出的电量，η为电池效率系数，主要包括充放电系数ηc和温度系数ηt.，它们之间满足η＝ηc·ηt函数关系，α满足如下函数关系：

图5为电池滞回电压特性图，由图可知，相同的soc值在充放电条件下对应的电池平衡电压是不同的，相同条件下充电状态平衡电压较高，其主要由电池本身的电化学特性造成。同时由图5可知，磷酸铁锂电池在平台区(soc为20％～90％)电压变化不大，非平台区(soc为0～20％，90％～100％)电压变化较大。因此，根据磷酸铁锂电池充放电开路电压与soc关系,可使用matlab中将其一阶拟合，在精度要求前提下降低运算复杂度。本发明计算基于开路平衡电压的电池初始soc时，首先判断电流方向，即判断系统处于充电状态还是放电状态，然后根据开路平衡电压代入拟合方程式中计算soc值。拟合方程具体表示如下：

由于电池开端电压在较长静置时间才恢复稳定，为提高soc估算精度，本发明采用上次关机系统记录的历史soc值与电池开路电压一阶拟合相结合的方法计算初始soc，如图6所示。由实验数据及经验可知，在静置2小时以后电池的端电压趋于稳定。因此，本发明将开路平衡电压sococv与上次关机时系统记录的历史soc值soch结合计算初始时刻soc，具体表示如下：

其中，soch为系统记录的上次关机时的电池soc值，sococv为基于电池电压滞回特性计算的soc，β为时间t相关系数，0≤β≤1，其表示如下：

s4、基于电池滞回电压特性和回弹电压特性建立的三阶等效电池模型的的状态空间模型，选择电池的工作电流i作为系统激励、负载电压uc作为系统观测变量、电池的极化电压和soc值作为状态变量建立状态方程和观测方程，并将其离散化；

将其表示为矩阵形式如下：

s5、基于建立的状态方程和观测方程，根据扩展卡尔曼递推原理不断对当前时刻soc值进行修正，使估算值逐渐逼近真实值。基于扩展卡尔曼滤波的soc估算具体步骤：

s51)获取扩展卡尔曼滤波状态估计向量初始值x0与误差协方差矩阵的初始值p0；

s52)将k-1时刻的状态向量估计值与误差协方差矩阵pk-1/k-1进行时间更新进一步确定k时刻的状态向量估计值与误差协方差pk/k-1:

其中，f为状态转移函数，为k-1时刻状态向量估计值，ik-1为k-1时刻系统输入向量。

s53)确定系统噪声的方差qk为:

其中wk为均值为0的高斯白噪声，为wk的转置，e为数学期望

s54)确定扩展卡尔曼增益矩阵kk

其中ck为雅可比矩阵，为ck的转置，rk为观测噪声协方差

vk为均值为0的高斯白噪声

s55)获得卡尔曼滤波状态向量估计值及相应协方差矩阵

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。