首页 / 专利库 / 人工智能 / 状态空间 / 基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法

基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法

阅读:241发布:2020-05-12

专利汇可以提供基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于自适应双 扩展卡 尔曼滤波法的SOC估算方法,首先建立锂 电池 的二阶RC等效 电路 模型;然后通过脉冲充放电实验确定锂电池不同 荷电状态 SOC处的开路 电压 以及电池等效模型参数,得到开路电压与SOC的函数关系,及其它模型参数与不同SOC处的关系,其他模型参数包括欧姆内阻、电化学极化 电阻 、电化学极化电容、浓度差极化电阻、浓度差极化电容值;建立以SOC、极化电压为状态变量的 状态空间 方程和以欧姆内阻为状态变量的状态空间方程;最终 迭代 计算实时获得锂电池的荷电状态SOC值。本发明解决了 现有技术 中存在的噪声统计特性未知的问题,同时利用卡尔曼滤波 算法 估算电池的欧姆内阻,提高了模型 精度 。,下面是基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法专利的具体信息内容。

1.一种基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法,其特征在于,具体按照以下步骤实施:
步骤1、建立锂电池的二阶RC等效电路模型;
步骤2、通过脉冲充放电实验确定锂电池不同荷电状态SOC处的开路电压以及电池等效模型参数,然后通过函数拟合得到开路电压与SOC的具体函数关系,构建出其它模型参数与不同SOC处的关系,其它模型参数包括欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值;
步骤3、根据锂电池的等效电路模型,分别建立以SOC、极化电压为状态变量的状态空间方程和以欧姆内阻为状态变量的状态空间方程;
步骤4、运用自适应扩展卡尔曼滤波算法估算电池的SOC,同时运用卡尔曼滤波算法估算电池的欧姆内阻,最终迭代计算实时获得锂电池的荷电状态SOC值。
2.根据权利要求1所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤1中锂电池的二阶RC等效电路模型包括两个串联在一起的RC并联电路,其中一个RC并联电路的电化学极化电阻RL与电化学极化电容CL之间与充放电方向的开路电压源连接后接入电压,另一个RC并联电路的浓度差极化电阻RS与浓度差极化电容CS共同串联欧姆内阻Ro后接入电压。
3.根据权利要求2所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤2中将开路电压按照充、放电方向分为两个部分,当处于放电状态时,放电开路电压UOC及其所串联的二极管工作,充电开路电压U’OC不工作;当处于充电状态时,充电开路电压U’OC及其所串联的二极管工作,放电开路电压UOC不工作,在锂电池充电或放电两个过程中,每个过程因充电或放电电流方向不同而选择各自开路电压源UOC或U’OC,由此分别得到充电和放电过程的开路电压与SOC的关系,以及欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值在放电过程中与SOC的对应关系。
4.根据权利要求3所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述建立的以SOC、极化电压为状态变量的锂电池SOC、电化学极化电压UL与浓度差极化电压US的状态空间方程为:
锂电池端电压方程如下:
U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+vk
其中,ΔT为采样时间,Ro为电池的欧姆内阻,RL、CL分别为电池的电化学极化电阻和极化电容,RS、CS分别为电池的浓度差极化电阻和极化电容,τL、τS分别代表电化学极化时间常数与浓度差极化时间常数,其中τL=RLCL,τS=RSCS,US为RS两端的电压,UL为RL两端的电压,UOC为电池的开路电压,i为电池的工作电流,U为电池的工作电压,wk为SOC估算过程噪声,vk为SOC估算测量噪声,CN代表锂电池当前容量,k代表当前时刻迭代计算步数,SOC(k)代表第k次计算中的荷电状态SOC值;
所述步骤3中建立的以欧姆内阻为状态变量的状态空间方程为:
Ro,k=Ro,k-1+rk
U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+qk
其中,rk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算过程噪声,qk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算测量噪声,Ro,k为第k次迭代计算中的欧姆内阻估计值,Ro,k-1为前一次即第k-1次的欧姆内阻估计值。
5.根据权利要求4所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4具体按照以下步骤实施:
步骤4.1、设定锂电池荷电状态SOC和卡尔曼滤波算法中的状态误差协方差P的初值,选择系统噪声值为10-4,由SOC初值以及步骤2中的其他模型参数与SOC对应关系得到欧姆内阻的初值;
步骤4.2、根据此时荷电状态SOC的估算结果,结合步骤2中的其它模型参数与SOC对应关系,得到等效电路模型中的电化学极化内阻RL与电容CL、浓度差极化电阻RS与电容CS的值;对于不在对应关系表内的荷电状态SOC值,采用左侧临近SOC对应的参数值代替;
步骤4.3、由荷电状态SOC与欧姆内阻Ro各自的状态方程,分别迭代计算荷电状态SOC和欧姆内阻的状态预测值和误差协方差预测值;
步骤4.4、计算荷电状态SOC和欧姆内阻的卡尔曼滤波增益P;
步骤4.5、将荷电状态SOC和欧姆内阻的状态预测值代入观测方程得到观测量的预测值;
步骤4.6、计算系统量测新息;
步骤4.7、由量测新息、荷电状态SOC和欧姆内阻各自的卡尔曼滤波增益得到当前时刻的状态估计量并更新误差协方差;
步骤4.8、调整荷电状态SOC估算系统的过程噪声协方差;
步骤4.9、将步骤4.3得到的状态量、误差协方差和步骤4.4得到的增益P、步骤4.5得到的观测量预测值以及步骤4.6得到的新息代入步骤4.2中,得到锂电池荷电状态SOC的值,并开始新一轮的循环迭代。
6.根据权利要求5所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4.3具体如下:
计算状态预测值式子为:
其中, 为当前时刻的状态预测值, 为前一时刻的状态量,uk为当前时刻的输入变量,f为状态预测函数,k表示当前时刻的迭代计算步数;
计算误差协方差预测式子为:
Pk/k-1=Fk-1Pk-1Fk-1T+Qk
其中,Pk/k-1为当前时刻的状态预测值,Pk-1为前一时刻的状态量,Qk为当前时刻的过程噪声协方差,Fk-1由状态预测函数 对前一时刻的状态预测值求雅克比矩阵得到。
7.根据权利要求6所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4.4具体如下:
计算卡尔曼滤波增益为:
Kk=Pk/k-1GkT(GkPk/k-1GkT+R)-1
其中,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,R为测量噪声协方差,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
8.根据权利要求7所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4.6具体如下:
计算系统量测新息:
υj=yj-yj/j-1
其中,υj表示系统量测新息,yj表示前面某一时刻的系统量测值,yj/j-1表示前面某一时刻的系统量测的预测值,j表示某一历史时刻迭代计算步数。
9.根据权利要求8所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4.7具体如下:
计算当前时刻的状态估计量:
其中,yk为当前时刻观测值,g为观测函数;
更新误差协方差:
Pk=(E-KkGk)Pk/k-1
其中,E代表单位矩阵,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
10.根据权利要求9所述的一种基于自适应双扩展卡尔曼滤波法的SOC估算方法,其特征在于,所述步骤4.8具体如下:
调整荷电状态SOC的过程噪声协方差:
Qk=KkSkKkT
其中,m表示统计量测新息历史信息的窗口宽度,Sk表示前m个量测新息的均值,υj表示系统量测新息,Kk为当前时刻的卡尔曼滤波增益,Qk为当前时刻的过程噪声协方差,k表示当前时刻迭代计算步数,j表示某一历史时刻计算步数。

说明书全文

基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法

技术领域

[0001] 本发明属于锂电池状态估算技术领域,具体涉及一种基于自适应双扩展卡尔曼滤波法的SOC估算方法。

背景技术

[0002] 随着清洁能源的发展,锂电池在、光储能、电动汽车等领域中获得了越来越多的应用。为了保证电池安全有效地运行,需要建立一套电池管理系统,实时监测电池的电压电流温度等参数,准确估算电池的荷电状态(State of Charge,SOC)、健康状态等状态信息。其中,准确估算SOC是其它状态估计的基础;可以避免过充过放,延长电池寿命;还可以帮助用户制定正确的产品使用计划,具有十分重要的意义。
[0003] 目前估算SOC的方法主要有安时积分法、开路电压法、神经网络法、扩展卡尔曼滤波法(Extend Kalman Filter,EKF)。安时积分法原理简单,易于实现,但是无法提供SOC的初值,而且存在积分运算使误差累积越来越大的问题。开路电压法(Open Circuit Voltage,OCV)依靠开路电压与SOC之间的SOC-OCV关系曲线,得到SOC的值,简单直接,但是由于需要测量开路电压,无法在电池充放电过程使用。神经网络法不需要考虑电池的内部结构,但前提是要用大量实验数据进行训练,估算精度依托于训练数据的针对性与全面性,而且此类算法普遍较为复杂,难以实现。EKF算法可以将安时积分法和开路电压法有效结合,对安时积分法的结果进行实时修正,消除累计误差,而且估算结果不依赖于SOC的初值,是目前研究的主流方法。但是EKF算法在估算过程中将系统噪声假定为确定的高斯白噪声,这与实际不符,可能会导致较大的估算误差。而且EKF算法是基于模型的算法,模型的精度很大程度上影响了估算结果的准确性,但是电池内部参数在不断变化,会使模型精度降低,这也会影响算法估算的准确性。

发明内容

[0004] 本发明的目的是提供一种基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法,解决了现有技术中存在的噪声统计特性未知的问题,同时利用卡尔曼滤波算法估算电池的欧姆内阻,提高了模型精度。
[0005] 本发明所采用的技术方案是,一种基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法,具体按照以下步骤实施:
[0006] 步骤1、建立锂电池的二阶RC等效电路模型;
[0007] 步骤2、通过脉冲充放电实验确定锂电池不同荷电状态SOC处的开路电压以及电池等效模型参数,然后通过函数拟合得到开路电压与SOC的具体函数关系,构建出其它模型参数与不同SOC处的关系,其它模型参数包括欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值;
[0008] 步骤3、根据锂电池的等效电路模型,分别建立以SOC、极化电压为状态变量的状态空间方程和以欧姆内阻为状态变量的状态空间方程;
[0009] 步骤4、运用自适应扩展卡尔曼滤波算法估算电池的SOC,同时运用卡尔曼滤波算法估算电池的欧姆内阻,最终迭代计算实时获得锂电池的荷电状态SOC值。
[0010] 本发明的特点还在于,
[0011] 步骤1中锂电池的二阶RC等效电路模型包括两个串联在一起的RC并联电路,其中一个RC并联电路的电化学极化电阻RL与电化学极化电容CL之间与充放电方向的开路电压源连接后接入电压,另一个RC并联电路的浓度差极化电阻RS与浓度差极化电容CS共同串联欧姆内阻Ro后接入电压。
[0012] 步骤2中将开路电压按照充、放电方向分为两个部分,当处于放电状态时,放电开路电压UOC及其所串联的二极管工作,充电开路电压U’OC不工作;当处于充电状态时,充电开路电压U’OC及其所串联的二极管工作,放电开路电压UOC不工作,在锂电池充电或放电两个过程中,每个过程因充电或放电电流方向不同而选择各自开路电压源UOC或U’OC,由此分别得到充电和放电过程的开路电压与SOC的关系,以及欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值在放电过程中与SOC的对应关系;
[0013] 步骤3中所述建立的以SOC、极化电压为状态变量的锂电池SOC、电化学极化电压UL与浓度差极化电压US的状态空间方程为:
[0014]
[0015] 锂电池端电压方程如下:
[0016] U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+vk
[0017] 其中,ΔT为采样时间,Ro为电池的欧姆内阻,RL、CL分别为电池的电化学极化电阻和极化电容,RS、CS分别为电池的浓度差极化电阻和极化电容,τL、τS分别代表电化学极化时间常数与浓度差极化时间常数,其中τL=RL CL,τS=RS CS,US为RS两端的电压,UL为RL两端的电压,UOC为电池的开路电压,i为电池的工作电流,U为电池的工作电压,wk为SOC估算过程噪声,vk为SOC估算测量噪声,CN代表锂电池当前容量,k代表当前时刻迭代计算步数,SOC(k)代表第k次计算中的荷电状态SOC值。
[0018] 所建立的以欧姆内阻为状态变量的状态空间方程为:
[0019] Ro,k=Ro,k-1+rk
[0020] U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+qk
[0021] 其中,rk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算过程噪声,qk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算测量噪声,Ro,k为第k次迭代计算中的欧姆内阻估计值,Ro,k-1为前一次即第k-1次的欧姆内阻估计值。
[0022] 步骤4具体按照以下步骤实施:
[0023] 步骤4.1、设定锂电池荷电状态SOC和卡尔曼滤波算法中的状态误差协方差P的初值,选择系统噪声值为10-4,由SOC初值以及步骤2中的其它模型参数与SOC对应关系得到欧姆内阻的初值;
[0024] 步骤4.2、根据此时荷电状态SOC的估算结果,结合步骤2中的其他模型参数与SOC对应关系,得到等效电路模型中的电化学极化内阻RL与电容CL、浓度差极化电阻RS与电容CS的值;对于不在对应关系表内的荷电状态SOC值,采用左侧临近SOC对应的参数值代替;
[0025] 步骤4.3、由荷电状态SOC与欧姆内阻Ro各自的状态方程,分别迭代计算SOC和欧姆内阻的状态预测值和误差协方差预测值;
[0026] 步骤4.4、计算荷电状态SOC和欧姆内阻的卡尔曼滤波增益P;
[0027] 步骤4.5、将荷电状态SOC和欧姆内阻的状态预测值代入观测方程得到观测量的预测值;
[0028] 步骤4.6、计算系统量测新息;
[0029] 步骤4.7、由量测新息,和荷电状态SOC与欧姆内阻各自的卡尔曼滤波增益得到当前时刻的状态估计量并更新误差协方差;
[0030] 步骤4.8、调整荷电状态SOC估算的过程噪声协方差;
[0031] 步骤4.9、将步骤4.3得到的状态量、误差协方差和步骤4.4得到的增益P、步骤4.5得到的观测量预测值以及步骤4.6得到的新息代入步骤4.2中,得到锂电池荷电状态SOC的值,并开始新一轮的循环迭代。
[0032] 步骤4.3具体如下:
[0033] 计算状态预测值式子为:
[0034]
[0035] 其中, 为当前时刻的状态预测值, 为前一时刻的状态量,uk为当前时刻的输入变量,f为状态预测函数,k表示当前时刻迭代计算步数;
[0036] 计算误差协方差预测式子为:
[0037] Pk/k-1=Fk-1Pk-1Fk-1T+Qk
[0038] 其中,Pk/k-1为当前时刻的状态预测值,Pk-1为前一时刻的状态量,Qk为当前时刻的过程噪声协方差,Fk-1由状态预测函数 对前一时刻的状态预测值求雅克比矩阵得到。
[0039] 步骤4.4具体如下:
[0040] 计算卡尔曼滤波增益为:
[0041] Kk=Pk/k-1GkT(GkPk/k-1GkT+R)-1
[0042] 其中,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,R为测量噪声协方差,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
[0043] 步骤4.6具体如下:
[0044] 计算系统量测新息:
[0045] υj=yj-yj/j-1
[0046] 其中,υj表示系统量测新息,yj表示前面某一时刻的系统量测值,yj/j-1表示前面某一时刻的系统量测的预测值,j表示某一历史时刻迭代计算步数。
[0047] 步骤4.7具体如下:
[0048] 计算当前时刻的状态估计量:
[0049]
[0050] 其中,yk为当前时刻观测值,g为观测函数;
[0051] 更新误差协方差:
[0052] Pk=(E-KkGk)Pk/k-1
[0053] 其中,E代表单位矩阵,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
[0054] 步骤4.8具体如下:
[0055] 调整荷电状态SOC估算的过程噪声协方差:
[0056]
[0057] Qk=KkSkKkT
[0058] 其中,m表示统计量测新息历史信息的窗口宽度,Sk表示前m个量测新息的均值,υj表示系统量测新息,Kk为当前时刻的卡尔曼滤波增益,Qk为当前时刻的过程噪声协方差,k表示当前时刻迭代计算步数,j表示某一历史时刻迭代计算步数。
[0059] 本发明的有益效果是:
[0060] (1)一种基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法,与传统安时积分法相比,本发明的估算结果不依赖于SOC初值,即使初值不准确,经过若干次迭代之后,也能快速收敛到真值。而且本发明具有抑制系统噪声的作用,可以消除积分累计误差的影响。
[0061] (2)与采用单一方向的SOC-OCV曲线相比,本发明在充、放电时分别采用充电和放电方向的SOC-OCV曲线,使电池的等效模型在充放电不断交替进行的工况条件下也能适用,算法的估算精度不会因电流方向的变化而下降。
[0062] (3)与一般的扩展卡尔曼滤波法相比,本发明在扩展卡尔曼滤波法的基础上采用自适应调整噪声协方差的方法,解决了实际中噪声统计特性未知的问题,能够适应不同工况,提高估算精度。除此以外,考虑到同时调整过程噪声协方差和测量噪声协方差容易造成滤波发散,本发明采用只调整过程噪声协方差的方法,降低了滤波发散的可能性。
[0063] (4)与对等效电路模型参数取定值的做法相比,本发明采用卡尔曼滤波算法对受电流、温度、老化情况等因素影响较大的欧姆内阻进行动态估算,能实时跟踪电池内部主要参数的变化,在电池内外条件改变的情况下也具有较高的估算精度。附图说明
[0064] 图1为锂电池二阶RC等效电路模型;
[0065] 图2为算法的原理图;
[0066] 图3为算法的流程图
[0067] 图4为动态应测试工况单周期电流;
[0068] 图5为采用本发明的算法估算SOC的估算结果和误差图。

具体实施方式

[0069] 下面结合附图和具体实施方式对本发明进行详细说明。
[0070] 本发明一种基于自适应双扩展卡尔曼滤波法的锂电池SOC估算方法,具体按照以下步骤实施:
[0071] 步骤1、如图1所示,建立锂电池的二阶RC等效电路模型;
[0072] 步骤1中锂电池的二阶RC等效电路模型包括两个串联在一起的RC并联电路,其中一个RC并联电路的电化学极化电阻RL与电化学极化电容CL之间与充放电方向的开路电压源连接后接入电压,另一个RC并联电路的浓度差极化电阻RS与浓度差极化电容CS共同串联欧姆内阻Ro后接入电压。
[0073] 步骤2、通过脉冲充放电实验确定锂电池不同荷电状态SOC处的开路电压以及电池等效模型参数,然后通过函数拟合得到开路电压与SOC的具体函数关系,构建出其它模型参数与不同SOC处的关系,其它模型参数包括欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值;
[0074] 电池的等效电路模型中,开路电压按照充、放电方向分为两个部分,每个部分因电流方向的不同与SOC具有不同的关系,需要分别得到充电和放电过程的开路电压与SOC的关系。
[0075] 步骤2中将开路电压按照充、放电方向分为两个部分,当处于放电状态时,放电开路电压UOC及其所串联的二极管工作,充电开路电压U’OC不工作;当处于充电状态时,充电开路电压U’OC及其所串联的二极管工作,放电开路电压UOC不工作,在锂电池充电或放电两个过程中,每个过程因充电或放电电流方向不同而选择各自开路电压源UOC或U’OC,由此分别得到充电和放电过程的开路电压与SOC的关系,以及欧姆内阻、电化学极化电阻、电化学极化电容、浓度差极化电阻、浓度差极化电容值在放电过程中与SOC的对应关系;
[0076] 步骤3、根据锂电池的等效电路模型,分别建立以SOC、极化电压为状态变量的状态空间方程和以欧姆内阻为状态变量的状态空间方程;
[0077] 步骤3中建立的以SOC、极化电压为状态变量的锂电池SOC、电化学极化电压UL与浓度差极化电压US的状态空间方程为:
[0078]
[0079] 锂电池端电压方程如下:
[0080] U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+vk
[0081] 其中,ΔT为采样时间,Ro为电池的欧姆内阻,RL、CL分别为电池的电化学极化电阻和极化电容,RS、CS分别为电池的浓度差极化电阻和极化电容,τL、τS分别代表电化学极化时间常数与浓度差极化时间常数,其中τL=RL CL,τS=RS CS,US为RS两端的电压,UL为RL两端的电压,UOC为电池的开路电压,i为电池的工作电流,U为电池的工作电压,wk为SOC过程噪声,vk为SOC测量噪声,CN代表锂电池当前容量,k代表当前时刻迭代计算步数,SOC(k)代表第k次计算中的荷电状态SOC值。
[0082] 所建立的以欧姆内阻为状态变量的状态空间方程为:
[0083] Ro,k=Ro,k-1+rk
[0084] U(k)=UOC(SOCk)-US(k)-UL(k)-Ro(k)i(k)+qk
[0085] 其中,rk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算的过程噪声,qk为以欧姆内阻Ro为状态变量的卡尔曼滤波法估算的测量噪声,Ro,k为第k次迭代计算中的欧姆内阻估计值,Ro,k-1为前一次即第k-1次的欧姆内阻估计值。
[0086] 步骤4、运用自适应扩展卡尔曼滤波算法估算电池的SOC,同时运用卡尔曼滤波算法估算电池的欧姆内阻,最终迭代计算实时获得锂电池的荷电状态SOC值。算法基本原理图如图2所示,按照该原理图,可以得到算法的执行流程,如图3所示。
[0087] 步骤4具体按照以下步骤实施:
[0088] 步骤4.1、设定锂电池荷电状态SOC和卡尔曼滤波算法中的状态误差协方差P的初-4值,选择系统噪声值为10 ,由SOC初值以及步骤2中的其他模型参数与SOC对应关系得到欧姆内阻的初值;
[0089] 步骤4.2、根据此时荷电状态SOC的估算结果,结合步骤2中的其它模型参数与SOC对应关系,得到等效电路模型中的电化学极化内阻RL与电容CL、浓度差极化电阻RS与电容CS的值;对于不在对应关系表内的荷电状态SOC值,采用左侧临近SOC对应的参数值代替;
[0090] 步骤4.3、由荷电状态SOC与欧姆内阻Ro各自的状态方程,分别迭代计算荷电状态SOC和欧姆内阻的状态预测值和误差协方差预测值;
[0091] 步骤4.4、计算荷电状态SOC和欧姆内阻的卡尔曼滤波增益P;
[0092] 步骤4.5、将荷电状态SOC和欧姆内阻的状态预测值代入观测方程得到观测量的预测值;
[0093] 步骤4.6、计算系统量测新息;
[0094] 步骤4.7、由量测新息、荷电状态SOC和欧姆内阻各自的卡尔曼滤波增益得到当前时刻的状态估计量并更新误差协方差;
[0095] 步骤4.8、调整荷电状态SOC估算的过程噪声协方差;
[0096] 步骤4.9、将步骤4.3得到的状态量、误差协方差和步骤4.4得到的增益P、步骤4.5得到的观测量预测值以及步骤4.6得到的新息代入步骤4.2中,得到锂电池荷电状态SOC的值,并开始新一轮的循环迭代。
[0097] 步骤4.3具体如下:
[0098] 计算状态预测值式子为:
[0099]
[0100] 其中, 为当前时刻的状态预测值, 为前一时刻的状态量,uk为当前时刻的输入变量,f为状态预测函数,k表示当前时刻的迭代计算步数;
[0101] 计算误差协方差预测式子为:
[0102] Pk/k-1=Fk-1Pk-1Fk-1T+Qk
[0103] 其中,Pk/k-1为当前时刻的状态预测值,Pk-1为前一时刻的状态量,Qk为当前时刻的过程噪声协方差,Fk-1由状态预测函数 对前一时刻的状态预测值求雅克比矩阵得到。
[0104] 步骤4.4具体如下:
[0105] 计算卡尔曼滤波增益为:
[0106] Kk=Pk/k-1GkT(GkPk/k-1GkT+R)-1
[0107] 其中,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,R为测量噪声协方差,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
[0108] 步骤4.6具体如下:
[0109] 计算系统量测新息:
[0110] υj=yj-yj/j-1
[0111] 其中,υj表示系统量测新息,yj表示前面某一时刻的系统量测值,yj/j-1表示前面某一时刻的系统量测的预测值,j表示某一历史时刻迭代计算步数。
[0112] 步骤4.7具体如下:
[0113] 计算当前时刻的状态估计量:
[0114]
[0115] 其中,yk为当前时刻观测值,g为观测函数;
[0116] 更新误差协方差:
[0117] Pk=(E-KkGk)Pk/k-1
[0118] 其中,E代表单位矩阵,Pk/k-1为当前时刻的状态预测值,Kk为当前时刻的卡尔曼滤波增益,Gk由观测函数对当前时刻的状态预测值求雅克比矩阵得到。
[0119] 步骤4.8具体如下:
[0120] 调整荷电状态SOC估算的过程噪声协方差:
[0121]
[0122] Qk=KkSkKkT
[0123] 其中,m表示统计量测新息历史信息的窗口宽度,Sk表示前m个量测新息的均值,υj表示系统量测新息,Kk为当前时刻的卡尔曼滤波增益,Qk为当前时刻的过程噪声协方差,k表示当前时刻迭代计算步数,j表示某一历史时刻计算步数。
[0124] 为了验证本发明估算SOC的准确度,以某种标称容量为1.7Ah的磷酸锂电池为研究对象,进行模拟工况实验。模拟的工况为动态应力测试工况DST(dynamic stress test),其单周期的充放电电流如图4所示。模拟工况实验完成之后,分别用扩展卡尔曼滤波法和本发明提出的自适应双扩展卡尔曼滤波算法对SOC进行估算,得到SOC估算结果和误差如图5所示,可以明显看出两种方法都能在SOC初值不准的情况下收敛到真值,但是自适应双扩展卡尔曼滤波算法的准确度更高。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈