CN107290688A - 一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 - Google Patents
一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 Download PDFInfo
- Publication number
- CN107290688A CN107290688A CN201710735042.3A CN201710735042A CN107290688A CN 107290688 A CN107290688 A CN 107290688A CN 201710735042 A CN201710735042 A CN 201710735042A CN 107290688 A CN107290688 A CN 107290688A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msubsup
- mtr
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法,对于不同的锂电池,确定卡尔曼滤波离散状态空间模型后,利用自适应模糊卡尔曼滤波的方法估计SOC,其中以锂电池等效模型中的端电压的残差方差匹配度和残差均值作为模糊控制***的输入,得到***噪声方差和量测噪声方差的调整系数,从而对这两个方差经行调整。再把调整后的***噪声方差和量测噪声方差代入到卡尔曼滤波算法中预估各时刻的SOC值。本发明可以精确的对动力锂电池荷点状态SOC经行估计,能解决现有估计方法存在的不能满足在线估计,累计误差大,发散,易受噪声影响等问题,估计精度高。
Description
技术领域
本发明涉及锂离子电池电荷预测技术领域,尤其涉及一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法。
背景技术
电动汽车作为未来主要的交通工具,对其启动、加速、爬坡等性能及续航里程具有一定的要求。而这些性能很大程度上取决于动力电池的性能。电池核电状态(SOC)是反映电池能量状态的重要参数。只有准确的估算电池SOC才能有效的提高动力汽车的利用率、优化驾驶,延长电池的使用寿命。而SOC是动力电池中隐含的状态量,难以直接测量和估算。影响SOC的因素主要有工作电压、电流、温度等,因而建立精确可靠的SOC估计方法才能获得准确的估计结果。
SOC的估计方法有很多种,有安时法、开路电压法、扩展卡尔曼滤波法等。安时法是一种评估电池剩余容量的经典方法,优点是简单易用,但缺点是必须知道电池荷电状态的初值。开路电压法是根据OCV-SOC曲线查表从而得到SOC的值。开路电压法虽然可以精确的估计SOC值,但是此方法不能实现实时在线估计。扩展卡尔曼滤波法是一种高效率的递归滤波器,它能够从一系列不完全包含噪声的测量中,估计动态***的状态。此方法的核心是由已建立的电池状态方程,根据电池工作时采集的电压值、电流值进行状态地推,得到SOC的实时估算值。然而,卡尔曼滤波需要预设噪声初值信息,不合适的噪声初值会是使估算结果不准确甚至发散。
发明内容
本发明目的就是为了弥补已有技术的缺陷,提供一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法。
本发明是通过以下技术方案实现的:
一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法,包括如下步骤:
(1)建立电池一阶RC模型确定状态方程和观测方程,确定电池模型的卡尔曼滤波离散状态空间模型及状态变量和观测变量;
(2)根据步骤(1)中确定的卡尔曼滤波离散状态空间模型计算该模型的状态估计更新值和噪声协方差矩阵的估计值;
(3)以某一时刻的端电压残差均值和端电压残差匹配度为Sugeno型模糊推理***的两个输入,输出***噪声方差阵调整系数和量测噪声方差阵调整系数;
(4)根据所得到的观测量端电压计算离散状态空间模型的离散状态协方差和卡尔曼增益矩阵,并提供状态变量估计值和滤波信息,计算状态变量的估计更新值;
(5)根据步骤(3)中的调整系数计算出当前的***噪声和量测噪声,并作为更新值进行状态误差先验值和卡尔曼增益矩阵的计算,从而实现模糊控制辅助卡尔曼滤波的地推计算,SOC估计值由每一步递推计算中的的状态估计更新值得出。
所述步骤(1)中的电池一阶RC模型为:
状态方程:
观测方程:Vk=F(SOCK)-R0Ik-uk+vk (17)
R0指一阶RC模型中电池欧姆内阻,R1指极化内阻。uk指电池极化电压,Ik指流过电池的电流,C值电池一阶RC模型中极化电容的容量,CN为电池额定容量,η为库伦效率,Vk为电池端电压,Ek为电池开路电动势,Ek与SOCk存在着非线性关系,wk、vk分别为***噪声和量测噪声;
将方程(16)(17)离散化可得:
其中的XK为k时刻的状态变量,YK为k时刻端电压的观测值,IK为k时刻的输入控制变量,ωk和Vk是互不相关的***噪声,式(18)中:
所述步骤(1)中,建立的电池一阶RC模型中需要辨识的参数为R0、R1、C、K0、K1,其中的R0、R1、C可以通过动力电池充放电采集电压、电流、温度数据,使用含有遗忘因子的地推最小二乘法在线辨识获得,其具体过程如下:
S1:对式(16)(17)经行离散化处理得到模型差分方程:
VK=m0+m1VK-1+m2IK (19)
式(19)中,m0,m1,m2为模型差分方程的待定系数,其值与模型中待辨识的参数成函数关系;
S2:将式(19)写成的形式,其中:
θ={m0,m1,m2} (21)
S3:基于遗忘因子的递推最小二乘参数辨识方法的具体估计过程:
确定最小二乘协方差P0和参数矩阵θ的初值;
确定最小二乘增益矩阵KK:
式中μ为最小二乘加权因子,hK为观测矩阵,根据随时间变化的最小二乘增益矩阵KK计算参数估计矩阵:
式中yK为K时刻端电压的测量值,θK为θK-1在k-1时刻对K时刻的参数估计值;根据KK计算协方差矩阵的更新为:
重复上述基于遗忘因子的递推最小二乘法的地推过程,直接辨识出m0,m1,m2的值,进而通过m0,m1,m2和R0、R1、C的函数关系得出R0、R1、C的值;
S4:K0、K1由电池充、放电的实验数据拟合得出,具体的过程为:先给电池充满电,然后静置5个小时;以C/3恒流放电,放出电池容量的10%后停止放电,搁置5小时,测量电池的开路电压;重复上一过程,直至放电至截止电压;在以C/3恒流充电,充到电池容量的10%后停止充电,搁置5小时后,测量电池的开路电压;重复上一过程,直到充电电流小于C/20;其中,充、放电相应开路电压的测量值的平均值作为电池开路电压EK,由10%间隔0.1s到100%对应的EK值和关系式:
EK=K0+K1ln(SOC),通过曲线拟合,求出K0,K1。
所述步骤(3)中的确定模糊推理***输入、输出值的方法,步骤如下:
步骤A,确定模糊推理***输入值的过程:
计算电压残差:
计算电压残差均值:
计算电压残差理论方差:
电压残差实际方差:
计算电压残差的匹配度:Hk=Mk-Nk (5)
其中,VK为K时刻的端电压值,为K时刻的端电压的估计值,DK为K时刻的端电压残差值,rk为K时刻端电压残差的平均值,MK为K时刻断电压残差的方差,NK为K端电压残差的实际方差值,CK,AK为K时刻的系数矩阵,QK为K时刻的***噪声方差,RK-1为K-1时刻的量测噪声方差,PK-1为K-1时刻的状态误差协方差,M,N为一段时间内统计次数,由式(1)(2)(3)(4)(5)可得到rK和HK,将所得到得rK和HK作为模糊推理***的两个输入;
步骤B,确定模糊推理***输出值的过程:
Sugeno模糊推理***可描述为:
其中是模糊集合,为常量。则***的总输出为:
wi为权重,每条规则Ri的重要性以权重的形式体现,为隶属度函数;
根据模糊推理***的模糊化方法对输入量Hk、rk进行模糊化,Hk和rk的隶属函数均采用三角形隶属函数,根据经验所得到的mf1、mf2、mf3分别代表三个模糊集合,Hk的论域为[-1,1],对应的mf1、mf2、mf3的参数值范围分别为[-1 -0.2],[-0.6 0 0.6],[1 1.7];rk的论域为[0,3],对应的mf1、mf2、mf3的参数值范围分别为[0 0.5],[0 0.5 3],[0.5 3];
设计两套模糊推理***分别输出调整系数α和β:
(1)输出为β的模糊推理***的第1、4、7条规则的输出参数
(2)输出为α的模糊推理***的第2、4、9条规则输出参数如下:
为了降低计算复杂度,将输出仅表示为残差方差匹配度的线性组合:
由式(8)(9)可得到调整系数α和β。
所述步骤(4)中,计算状态更新值,其步骤如下:
将模糊推理***得到的调整系数α、β代入(10),(11):
Qk=βQk-1 (10)
Rk=αRk-1 (11)
其中Qk为***噪声方差,Rk为量测噪声方差;
对***噪声方差QK和量测噪声方差RK经行调整;再将调整后的QK和RK代入(12),(13)
Pk/k-1=APk-1AT+Qk-1 (13)
将所得的结果进一步代入(14):
Pk=(I-KkCk)Pk/k-1 (14)
由和所得Kk可进一步得到:
所得到的即为当前SOC估计值;
其中,Kk为滤波增益,分别是当前状态的先验估计和现状态的最优估计值,Pk-1、Pk/k-1分别是状态误差协方差和状态误差协方差的先验值,VK为端电压,为VK的估计值。式中A,B,C,如下:
本发明的优点是:本发明可以精确的对动力锂电池荷点状态SOC经行估计,能解决现有估计方法存在的不能满足在线估计,累计误差大,发散,易受噪声影响等问题,估计精度高。
附图说明
图1为本发明实施的流程图。
图2为Sugeno型模糊推理***的流程图。
具体实施方式
如图1所示,一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法,包括如下步骤:
(1)建立电池一阶RC模型确定状态方程和观测方程,确定电池模型的卡尔曼滤波离散状态空间模型及状态变量和观测变量;
(2)根据步骤(1)中确定的卡尔曼滤波离散状态空间模型计算该模型的状态估计更新值和噪声协方差矩阵的估计值;
(3)以某一时刻的端电压残差均值和端电压残差匹配度为Sugeno型模糊推理***的两个输入,输出***噪声方差阵调整系数和量测噪声方差阵调整系数;
(4)根据所得到的观测量端电压计算离散状态空间模型的离散状态协方差和卡尔曼增益矩阵,并提供状态变量估计值和滤波信息,计算状态变量的估计更新值;
(5)根据步骤(3)中的调整系数计算出当前的***噪声和量测噪声,并作为更新值进行状态误差先验值和卡尔曼增益矩阵的计算,从而实现模糊控制辅助卡尔曼滤波的地推计算,SOC估计值由每一步递推计算中的的状态估计更新值得出。
所述步骤(1)中的电池一阶RC模型为:
状态方程:
观测方程:Vk=F(SOCK)-R0Ik-uk+vk (17)
R0指一阶RC模型中电池欧姆内阻,R1指极化内阻。uk指电池极化电压,Ik指流过电池的电流,C值电池一阶RC模型中极化电容的容量,CN为电池额定容量,η为库伦效率,Vk为电池端电压,Ek为电池开路电动势,Ek与SOCk存在着非线性关系,wk、vk分别为***噪声和量测噪声;
将方程(16)(17)离散化可得:
其中的XK为k时刻的状态变量,YK为k时刻端电压的观测值,IK为k时刻的输入控制变量,ωk和Vk是互不相关的***噪声,式(18)中:
所述步骤(1)中,建立的电池一阶RC模型中需要辨识的参数为R0、R1、C、K0、K1,其中的R0、R1、C可以通过动力电池充放电采集电压、电流、温度数据,使用含有遗忘因子的地推最小二乘法在线辨识获得,其具体过程如下:
S1:对式(16)(17)经行离散化处理得到模型差分方程:
VK=m0+m1VK-1+m2IK (19)
式(19)中,m0,m1,m2为模型差分方程的待定系数,其值与模型中待辨识的参数成函数关系;
S2:将式(19)写成的形式,其中:
θ={m0,m1,m2} (21)
S3:基于遗忘因子的递推最小二乘参数辨识方法的具体估计过程:
确定最小二乘协方差P0和参数矩阵θ的初值;
确定最小二乘增益矩阵KK:
式中μ为最小二乘加权因子,hK为观测矩阵,根据随时间变化的最小二乘增益矩阵KK计算参数估计矩阵:
式中yK为K时刻端电压的测量值,θK为θK-1在k-1时刻对K时刻的参数估计值;根据KK计算协方差矩阵的更新为:
重复上述基于遗忘因子的递推最小二乘法的地推过程,直接辨识出m0,m1,m2的值,进而通过m0,m1,m2和R0、R1、C的函数关系得出R0、R1、C的值;
S4:K0、K1由电池充、放电的实验数据拟合得出,具体的过程为:先给电池充满电,然后静置5个小时;以C/3恒流放电,放出电池容量的10%后停止放电,搁置5小时,测量电池的开路电压;重复上一过程,直至放电至截止电压;在以C/3恒流充电,充到电池容量的10%后停止充电,搁置5小时后,测量电池的开路电压;重复上一过程,直到充电电流小于C/20;其中,充、放电相应开路电压的测量值的平均值作为电池开路电压EK,由10%间隔0.1s到100%对应的EK值和关系式:
EK=K0+K1ln(SOC),通过曲线拟合,求出K0,K1。
所述步骤(3)中的确定模糊推理***输入、输出值的方法,步骤如下:
步骤A,确定模糊推理***输入值的过程:
计算电压残差:
计算电压残差均值:
计算电压残差理论方差:
电压残差实际方差:
计算电压残差的匹配度:Hk=Mk-Nk (5)
其中,VK为K时刻的端电压值,为K时刻的端电压的估计值,DK为K时刻的端电压残差值,rk为K时刻端电压残差的平均值,MK为K时刻断电压残差的方差,NK为K端电压残差的实际方差值,CK,AK为K时刻的系数矩阵,QK为K时刻的***噪声方差,RK-1为K-1时刻的量测噪声方差,PK-1为K-1时刻的状态误差协方差,M,N为一段时间内统计次数;
由式(1)(2)(3)(4)(5)可得到rK和HK,将所得到得rK和HK作为模糊推理***的两个输入;
步骤B,确定模糊推理***输出值的过程:
Sugeno模糊推理***可描述为:
其中是模糊集合,为常量。则***的总输出为:
wi为权重,每条规则Ri的重要性以权重的形式体现,为隶属度函数;
根据模糊推理***的模糊化方法对输入量Hk、rk进行模糊化,Hk和rk的隶属函数均采用三角形隶属函数,根据经验所得到的mf1、mf2、mf3分别代表三个模糊集合,Hk的论域为[-1,1],对应的mf1、mf2、mf3的参数值范围分别为[-1 -0.2],[-0.6 0 0.6],[1 1.7];rk的论域为[0,3],对应的mf1、mf2、mf3的参数值范围分别为[0 0.5],[0 0.5 3],[0.5 3];
设计两套模糊推理***分别输出调整系数α和β:
(1)输出为β的模糊推理***的第1、4、7条规则的输出参数
(2)输出为α的模糊推理***的第2、4、9条规则输出参数如下:
为了降低计算复杂度,将输出仅表示为残差方差匹配度的线性组合:
由式(8)(9)可得到调整系数α和β。
所述步骤(4)中,计算状态更新值,其步骤如下:
将模糊推理***得到的调整系数α、β代入(10),(11):
Qk=βQk-1 (10)
Rk=αRk-1 (11)
其中Qk为***噪声方差,Rk为量测噪声方差;
对***噪声方差QK和量测噪声方差RK经行调整;再将调整后的QK和RK代入(12),(13)
Pk/k-1=APk-1AT+Qk-1 (13)
将所得的结果进一步代入(14):
Pk=(I-KkCk)Pk/k-1 (14)
由和所得Kk可进一步得到:
所得到的即为当前SOC估计值;
其中,Kk为滤波增益,分别是当前状态的先验估计和现状态的最优估计值,Pk-1、Pk/k-1分别是状态误差协方差和状态误差协方差的先验值,VK为端电压,为VK的估计值。式中A,B,C,如下:
如图2所示,重复多次试验,将每组实验结果与先前结果进行纵向对比,从而得到效果更好的一组数据。最后经过100次左右的仿真实验,确定了模糊逻辑推理***的相应输出参数。
本发明提供了锂电池估计方法,基于自适应模糊推理辅助卡尔曼滤波滤波进行估计,根据电池电流、电压采集数据,计算端电压残差的均值和残差方差匹配度做为模糊推理***的两个输入,输出调整系数,从而对***噪声和量测噪声进行调整,进而对卡尔曼滤波增益和状态误差协方差经行调整,结合了模糊推理和卡尔曼滤波两种算法。该步骤中的卡尔曼滤波增益和状态误差协方差是根据模糊推理***在线辨识获得,能更好反应实际噪声的统计特性,因此其辅助卡尔曼滤波算法具有较高的精度和准确行。
以上所述仅是本发明优先实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可做出适当的改进和润饰,这些改进和润饰也是本发明的保护范围。
Claims (3)
1.一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法,其特征在于:包括如下步骤:
(1)建立电池一阶RC模型确定状态方程和观测方程,确定电池模型的卡尔曼滤波离散状态空间模型及状态变量和观测变量;
(2)根据步骤(1)中确定的卡尔曼滤波离散状态空间模型计算该模型的状态估计更新值和噪声协方差矩阵的估计值;
(3)以某一时刻的端电压残差均值和端电压残差匹配度为Sugeno型模糊推理***的两个输入,输出***噪声方差阵调整系数和量测噪声方差阵调整系数;
(4)根据所得到的观测量端电压计算离散状态空间模型的离散状态协方差和卡尔曼增益矩阵,并提供状态变量估计值和滤波信息,计算状态变量的估计更新值;
(5)根据步骤(3)中的调整系数计算出当前的***噪声和量测噪声,并作为更新值进行状态误差先验值和卡尔曼增益矩阵的计算,从而实现模糊控制辅助卡尔曼滤波的地推计算,SOC估计值由每一步递推计算中的的状态估计更新值得出。
2.根据权利要求1中所述的一种基于自适应模糊卡尔曼滤波的SOC估计的方法,其特征在于:所述步骤(3)中的确定模糊控制***输入、输出值的方法,步骤如下:
步骤A,确定模糊推理***输入值的过程:
计算电压残差:
计算电压残差均值:
计算电压残差理论方差:
电压残差实际方差:
计算电压残差方差匹配度:Hk=Mk-Nk (5)
其中,VK为K时刻的端电压值,为K时刻的端电压的估计值,DK为K时刻的端电压残差值,rk为K时刻端电压残差的平均值,MK为K时刻断电压残差的方差,NK为K端电压残差的实际方差值,CK,AK为K时刻的系数矩阵,QK为K时刻的***噪声方差,RK-1为K-1时刻的量测噪声方差,PK-1为K-1时刻的状态误差协方差,M,N为一段时间内的统计次数,由式(1)(2)(3)(4)(5)可得到rK和HK,将所得到得rK和HK作为模糊推理***的两个输入;
步骤B,确定模糊推理***输出值的过程:
Sugeno模糊推理***可描述为:
Ri:
其中是模糊集合,为常量。则***的总输出为:
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<msubsup>
<mi>A</mi>
<mn>1</mn>
<mi>i</mi>
</msubsup>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&Pi;</mo>
<msub>
<mi>u</mi>
<msubsup>
<mi>A</mi>
<mn>2</mn>
<mi>i</mi>
</msubsup>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
wi为权重,为隶属度函数;
根据模糊推理***的模糊化方法对输入量Hk、rk进行模糊化,Hk和rk的隶属函数均采用三角形隶属函数,根据得到的mf1、mf2、mf3分别代表三个模糊集合;
为降低计算复杂度,将输出仅表示为残差方差匹配度的线性组合,如下式所示:
<mrow>
<msup>
<mi>&beta;</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>c</mi>
<mn>1</mn>
<mi>i</mi>
</msubsup>
<mo>*</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>c</mi>
<mn>0</mn>
<mi>i</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&alpha;</mi>
<mi>i</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>d</mi>
<mn>1</mn>
<mi>i</mi>
</msubsup>
<mo>*</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>d</mi>
<mn>0</mn>
<mi>i</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
设计两套模糊推理***分别输出调整系数α和β,输出为α的模糊推理***输出参数输出为β的模糊推理***输出参数由式(8)(9)可得到调整系数α和β的值。
3.根据权利要求2中所述的一种基于自适应模糊卡尔曼滤波的锂电池SOC估计方法,其特征在于:所述的步骤(5)中计算SOC估计值,其步骤如下:
将模糊推理***得到的调整系数α、β代入(10),(11)式:
Qk=βQk-1 (10)
Rk=αRk-1 (11)
其中Qk为***噪声方差,Rk为量测噪声方差;
对***噪声方差QK和量测噪声方差RK经行调整;再将调整后的QK和RK代入(12),(13)式中:
<mrow>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>C</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>C</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
Pk/k-1=APk-1AT+Qk-1 (13)
将所得的结果进一步代入(14):
Pk=(I-KkCk)Pk/k-1 (14)
由和所得Kk可进一步得到:
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
所得到的即为当前SOC估计值;
其中,Kk为滤波增益,分别是当前状态的先验估计和现状态的最优估计值,Pk-1、Pk/k-1分别是状态误差协方差和状态误差协方差的先验值,VK为端电压,为VK的估计值。式中A,B,C,如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>C</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<msub>
<mi>C</mi>
<mi>N</mi>
</msub>
</mfrac>
<mi>&eta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>C</mi>
</mrow>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>SOC</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>SOC</mi>
<mi>k</mi>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>SOC</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710735042.3A CN107290688B (zh) | 2017-08-24 | 2017-08-24 | 一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710735042.3A CN107290688B (zh) | 2017-08-24 | 2017-08-24 | 一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107290688A true CN107290688A (zh) | 2017-10-24 |
CN107290688B CN107290688B (zh) | 2019-08-09 |
Family
ID=60107252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710735042.3A Active CN107290688B (zh) | 2017-08-24 | 2017-08-24 | 一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107290688B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108427079A (zh) * | 2018-06-11 | 2018-08-21 | 西南交通大学 | 一种动力电池剩余电量估计方法 |
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN109617839A (zh) * | 2018-11-21 | 2019-04-12 | 重庆邮电大学 | 一种基于卡尔曼滤波算法的莫尔斯信号检测方法 |
CN110048444A (zh) * | 2019-05-31 | 2019-07-23 | 闽江学院 | 一种基于soc状态估计的铅酸电池模糊控制方法 |
CN110133524A (zh) * | 2019-04-19 | 2019-08-16 | 顺丰科技有限公司 | 电池电荷状态计算方法、装置、服务器及介质 |
CN110187282A (zh) * | 2019-06-03 | 2019-08-30 | 珠海东帆科技有限公司 | 电池荷电状态估算方法和估算装置 |
CN110320473A (zh) * | 2019-07-26 | 2019-10-11 | 上海理工大学 | 一种基于卡尔曼滤波及模糊逻辑汽车锂电池容量估计方法 |
CN110395141A (zh) * | 2019-06-27 | 2019-11-01 | 武汉理工大学 | 基于自适应卡尔曼滤波法的动力锂电池soc估算方法 |
CN111751750A (zh) * | 2020-06-19 | 2020-10-09 | 杭州电子科技大学 | 基于模糊ekf的多阶段闭环锂电池soc估算方法 |
CN111781503A (zh) * | 2020-06-15 | 2020-10-16 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN111811503A (zh) * | 2020-07-15 | 2020-10-23 | 桂林电子科技大学 | 基于超宽带和二维码的无迹卡尔曼滤波融合定位方法 |
CN111983472A (zh) * | 2020-08-24 | 2020-11-24 | 哈尔滨理工大学 | 一种基于自适应卡尔曼滤波的锂离子动力电池安全度估算方法及估算装置 |
CN112130077A (zh) * | 2020-09-30 | 2020-12-25 | 东风汽车集团有限公司 | 一种不同工况下动力电池组的soc估算方法 |
CN112310985A (zh) * | 2020-07-06 | 2021-02-02 | 长沙理工大学 | 考虑储能***健康状态的风电功率平滑控制策略 |
CN112748665A (zh) * | 2020-12-22 | 2021-05-04 | 国网江苏省电力有限公司电力科学研究院 | 基于模糊卡尔曼滤波的氢燃料电池迭代控制方法及装置 |
CN112816879A (zh) * | 2021-01-05 | 2021-05-18 | 南京航空航天大学 | 一种用于变工况车用动力电池SoE在线估计方法 |
CN112858929A (zh) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | 一种基于模糊逻辑与扩展卡尔曼滤波的电池soc估计方法 |
CN113190969A (zh) * | 2021-04-07 | 2021-07-30 | 四川大学 | 一种基于信息评估机制的锂电池模型参数辨识方法 |
CN113758503A (zh) * | 2021-08-12 | 2021-12-07 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
CN114625013A (zh) * | 2022-03-31 | 2022-06-14 | 东风汽车集团股份有限公司 | 一种能量回馈的控制方法和三输入一输出模糊控制器 |
CN115598541A (zh) * | 2022-12-01 | 2023-01-13 | 西南科技大学(Cn) | 基于遗忘因子自适应反馈修正的电池能量状态评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102831100A (zh) * | 2012-07-18 | 2012-12-19 | 深圳职业技术学院 | 电池荷电状态估算方法及装置 |
US20150377974A1 (en) * | 2014-06-30 | 2015-12-31 | Foundation of Soongsil University-lndustry Cooperation | Battery state estimation method and system using dual extended kalman filter, and recording medium for performing the method |
CN105954679A (zh) * | 2016-04-29 | 2016-09-21 | 华中科技大学 | 一种锂电池荷电状态的在线估计方法 |
CN106054084A (zh) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | 一种动力电池soc估计方法 |
-
2017
- 2017-08-24 CN CN201710735042.3A patent/CN107290688B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102831100A (zh) * | 2012-07-18 | 2012-12-19 | 深圳职业技术学院 | 电池荷电状态估算方法及装置 |
US20150377974A1 (en) * | 2014-06-30 | 2015-12-31 | Foundation of Soongsil University-lndustry Cooperation | Battery state estimation method and system using dual extended kalman filter, and recording medium for performing the method |
CN105954679A (zh) * | 2016-04-29 | 2016-09-21 | 华中科技大学 | 一种锂电池荷电状态的在线估计方法 |
CN106054084A (zh) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | 一种动力电池soc估计方法 |
Non-Patent Citations (3)
Title |
---|
曲润峰等: "方差模糊自适应卡尔曼滤波及其应用", 《弹箭与制导学报》》 * |
朱浩等: "锂离子电池参数辨识与SOC估算研究", 《湖南大学学报》 * |
林程等: "基于模糊卡尔曼滤波算法的动力电池SOC估计", 《电源技术》 * |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN108427079A (zh) * | 2018-06-11 | 2018-08-21 | 西南交通大学 | 一种动力电池剩余电量估计方法 |
CN109617839B (zh) * | 2018-11-21 | 2021-03-02 | 重庆邮电大学 | 一种基于卡尔曼滤波算法的莫尔斯信号检测方法 |
CN109617839A (zh) * | 2018-11-21 | 2019-04-12 | 重庆邮电大学 | 一种基于卡尔曼滤波算法的莫尔斯信号检测方法 |
CN110133524A (zh) * | 2019-04-19 | 2019-08-16 | 顺丰科技有限公司 | 电池电荷状态计算方法、装置、服务器及介质 |
CN110133524B (zh) * | 2019-04-19 | 2022-07-08 | 丰翼科技(深圳)有限公司 | 电池电荷状态计算方法、装置、服务器及介质 |
CN110048444A (zh) * | 2019-05-31 | 2019-07-23 | 闽江学院 | 一种基于soc状态估计的铅酸电池模糊控制方法 |
CN110048444B (zh) * | 2019-05-31 | 2020-07-31 | 闽江学院 | 一种基于soc状态估计的铅酸电池模糊控制方法 |
CN110187282A (zh) * | 2019-06-03 | 2019-08-30 | 珠海东帆科技有限公司 | 电池荷电状态估算方法和估算装置 |
CN110187282B (zh) * | 2019-06-03 | 2021-09-10 | 珠海东帆科技有限公司 | 电池荷电状态估算方法和估算装置 |
CN110395141A (zh) * | 2019-06-27 | 2019-11-01 | 武汉理工大学 | 基于自适应卡尔曼滤波法的动力锂电池soc估算方法 |
CN110320473B (zh) * | 2019-07-26 | 2021-05-11 | 上海理工大学 | 一种基于卡尔曼滤波及模糊逻辑汽车锂电池容量估计方法 |
CN110320473A (zh) * | 2019-07-26 | 2019-10-11 | 上海理工大学 | 一种基于卡尔曼滤波及模糊逻辑汽车锂电池容量估计方法 |
CN111781503A (zh) * | 2020-06-15 | 2020-10-16 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN111781503B (zh) * | 2020-06-15 | 2023-09-22 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN111751750B (zh) * | 2020-06-19 | 2023-10-20 | 杭州电子科技大学 | 基于模糊ekf的多阶段闭环锂电池soc估算方法 |
CN111751750A (zh) * | 2020-06-19 | 2020-10-09 | 杭州电子科技大学 | 基于模糊ekf的多阶段闭环锂电池soc估算方法 |
CN112310985A (zh) * | 2020-07-06 | 2021-02-02 | 长沙理工大学 | 考虑储能***健康状态的风电功率平滑控制策略 |
CN112310985B (zh) * | 2020-07-06 | 2022-07-12 | 长沙理工大学 | 考虑储能***健康状态的风电功率平滑控制策略 |
CN111811503A (zh) * | 2020-07-15 | 2020-10-23 | 桂林电子科技大学 | 基于超宽带和二维码的无迹卡尔曼滤波融合定位方法 |
CN111983472B (zh) * | 2020-08-24 | 2022-11-25 | 哈尔滨理工大学 | 一种基于自适应卡尔曼滤波的锂离子动力电池安全度估算方法及估算装置 |
CN111983472A (zh) * | 2020-08-24 | 2020-11-24 | 哈尔滨理工大学 | 一种基于自适应卡尔曼滤波的锂离子动力电池安全度估算方法及估算装置 |
CN112130077A (zh) * | 2020-09-30 | 2020-12-25 | 东风汽车集团有限公司 | 一种不同工况下动力电池组的soc估算方法 |
CN112748665A (zh) * | 2020-12-22 | 2021-05-04 | 国网江苏省电力有限公司电力科学研究院 | 基于模糊卡尔曼滤波的氢燃料电池迭代控制方法及装置 |
CN112748665B (zh) * | 2020-12-22 | 2022-06-14 | 国网江苏省电力有限公司电力科学研究院 | 基于模糊卡尔曼滤波的氢燃料电池迭代控制方法及装置 |
CN112816879A (zh) * | 2021-01-05 | 2021-05-18 | 南京航空航天大学 | 一种用于变工况车用动力电池SoE在线估计方法 |
CN112858929A (zh) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | 一种基于模糊逻辑与扩展卡尔曼滤波的电池soc估计方法 |
CN113190969A (zh) * | 2021-04-07 | 2021-07-30 | 四川大学 | 一种基于信息评估机制的锂电池模型参数辨识方法 |
CN113190969B (zh) * | 2021-04-07 | 2022-02-11 | 四川大学 | 一种基于信息评估机制的锂电池模型参数辨识方法 |
CN113758503A (zh) * | 2021-08-12 | 2021-12-07 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
CN113758503B (zh) * | 2021-08-12 | 2022-03-18 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
CN114625013A (zh) * | 2022-03-31 | 2022-06-14 | 东风汽车集团股份有限公司 | 一种能量回馈的控制方法和三输入一输出模糊控制器 |
CN114625013B (zh) * | 2022-03-31 | 2024-04-16 | 东风汽车集团股份有限公司 | 一种能量回馈的控制方法和三输入一输出模糊控制器 |
CN115598541A (zh) * | 2022-12-01 | 2023-01-13 | 西南科技大学(Cn) | 基于遗忘因子自适应反馈修正的电池能量状态评估方法 |
CN115598541B (zh) * | 2022-12-01 | 2023-03-10 | 西南科技大学 | 基于遗忘因子自适应反馈修正的电池能量状态评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107290688B (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107290688B (zh) | 一种基于自适应模糊卡尔曼滤波的锂电池soc估计方法 | |
CN104535932B (zh) | 一种锂离子电池荷电状态估计方法 | |
CN106814329A (zh) | 一种基于双卡尔曼滤波算法的电池soc在线估计方法 | |
Song et al. | Improved SOC estimation of lithium-ion batteries with novel SOC-OCV curve estimation method using equivalent circuit model | |
CN105954679B (zh) | 一种锂电池荷电状态的在线估计方法 | |
CN103472398B (zh) | 基于扩展卡尔曼粒子滤波算法的动力电池soc估计方法 | |
CN102831100B (zh) | 电池荷电状态估算方法及装置 | |
CN103675706B (zh) | 一种动力电池电荷量估算方法 | |
CN106405433A (zh) | 一种基于扩展卡尔曼粒子滤波的soc估计方法及*** | |
CN107576919A (zh) | 基于armax模型的动力电池荷电状态估算***及方法 | |
CN101625397B (zh) | 一种电池剩余电量的混合快速估计方法 | |
CN102788957B (zh) | 一种动力电池荷电状态估算方法 | |
CN111323719A (zh) | 一种电动汽车动力电池组健康状态在线确定方法和*** | |
CN104267354B (zh) | 一种动力电池的峰值功率预测方法 | |
CN105425153B (zh) | 一种估计电动车辆的动力电池的荷电状态的方法 | |
CN107367699A (zh) | 一种基于分数阶模型的锂电池soc估算新方法 | |
CN106842056B (zh) | 基于双步在线智能优化算法的动力电池峰值功率估计方法 | |
CN105974327A (zh) | 一种基于神经网络和ukf的锂电池组soc预测方法 | |
CN105425154B (zh) | 一种估计电动汽车的动力电池组的荷电状态的方法 | |
CN109633479B (zh) | 基于嵌入式容积卡尔曼滤波的锂电池soc在线估算方法 | |
CN103020445A (zh) | 一种电动车车载磷酸铁锂电池的soc与soh预测方法 | |
CN102169168B (zh) | 一种基于粒子滤波的电池剩余电量估计方法 | |
CN108872869B (zh) | 一种基于bp神经网络的锂离子电池劣化分类方法 | |
CN104122504A (zh) | 一种电池的soc估算方法 | |
CN107589379A (zh) | 一种在线估计锂电池soc和阻抗的方法 |
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 |