CN104714188A - 量测噪声方差阵相匹配的电池电荷状态估计方法与*** - Google Patents
量测噪声方差阵相匹配的电池电荷状态估计方法与*** Download PDFInfo
- Publication number
- CN104714188A CN104714188A CN201510147926.8A CN201510147926A CN104714188A CN 104714188 A CN104714188 A CN 104714188A CN 201510147926 A CN201510147926 A CN 201510147926A CN 104714188 A CN104714188 A CN 104714188A
- Authority
- CN
- China
- Prior art keywords
- soc
- battle array
- moment
- chi
- battery
- 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
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000011159 matrix material Substances 0.000 title claims abstract description 30
- 230000003044 adaptive effect Effects 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 claims description 37
- 238000005259 measurement Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 17
- 108090000623 proteins and genes Proteins 0.000 claims description 17
- 230000006978 adaptation Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 4
- 230000009897 systematic effect Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000013011 mating Effects 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005530 etching Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Landscapes
- Tests Of Electric Status Of Batteries (AREA)
Abstract
本发明为量测噪声方差阵相匹配的电池电荷状态估计方法与***,本法步骤为:建立电池模型、电压电流采样,最小二乘法参数辨识模型参数;自适应无迹卡尔曼滤波估计SOC;调整滤波增益的方差阵相匹配的鲁棒自适应无迹卡尔曼滤波,计算***量测噪声实际残差方差阵值和估计的理论残差方差阵,得到εk。εk大于1自适应调整因子否则调整滤波增益,求得SOC估计值本***动力电池所接电压、电流传感器经模数转换模块连接微控制器。微控制器含模型参数辨识模块、调整因子计算模块、滤波增益自适应调整、鲁棒自适应无迹卡尔曼滤波SOC估计模块。SOC直接显示或传送到设备CAN网络。本发明动态调整滤波增益,提高了SOC估计的鲁棒性。
Description
技术领域
本发明涉及电动汽车动力电池电荷状态估计领域,具体为量测噪声方差阵相匹配的电池电荷状态估计方法与***。
背景技术
电池是电动汽车整车***中最重要组件之一,为了确保电动汽车的安全性,合理管理电池至关重要。近年来,锂离子电池以其能量密度高,循环寿命长,无记忆效应,环境友好等优点成为动力电池应用领域研究的热点,并逐渐成为电动车辆动力电池的最佳选择。
电荷状态(State of Charge,SOC)是电池管理***(Battery ManagementSystem,BMS)中的重要参数,然而它不能直接通过传感器测量所得,只能依据所建立的模型,运用相应的算法,间接估计得到。由于电池运行时受充、放电电流、温度、自放电和寿命等因素的影响,并且在使用过程中表现出的复杂非线性,对精确估计SOC造成一定难度。
目前SOC估计方法主要有:存在累积误差问题的安时法、需长时间静置的开路电压法、卡尔曼滤波法、需大量数据进行学习的模糊逻辑法和神经网络法等。
对于电池内部的化学特性表现出的复杂非线性特点,将经典卡尔曼滤波的状态估计扩展到非线性随机***的状态估计中,出现了扩展卡尔曼滤波和无迹卡尔曼滤波(Unscented Kalman Filter,UKF)。UKF克服了扩展卡尔曼滤波中需要计算雅克比矩阵,要求状态函数和量测函数必须是连续可微的缺点,但应用扩展卡尔曼滤波是以已知噪声统计特性为先决条件。当噪声统计特性随实际工况条件剧烈变化时,无迹卡尔曼滤波器不具有自适应能力,易出现滤波的状态估计精度下降甚至发散。为了克服这个缺点,又出现了自适应无迹卡尔曼滤波(Adaptive Unscented Kalman Filter,AUKF)。
AUKF虽然同时对量测噪声方差阵R和过程噪声方差阵Q进行调节时,但R与Q的调整未能保证实际残差方差阵与估计得到的理论残差方差阵完全相匹配,并未区别对待***过程噪声引起的误差和量测噪声引起的误差,其鲁棒性较差。
发明内容
本发明的目的是设计一种量测噪声方差阵相匹配的电池电荷状态估计方法,R本身自适应调整的方法结合实际残差方差阵与估计的理论残差方差阵相匹配的方法,以方差阵的完全匹配为目标,反向推算由虚拟的量测噪声方差阵与量测噪声方差阵R之比构成的调整因子εk,实现滤波增益的自适应调整,即方差阵相匹配的鲁棒自适应无迹卡尔曼滤波算法。对滤波增益进行自适应调节,解决过程噪声引起的误差、量测噪声引起的误差所造成的鲁棒性较差的问题,准确估计电池电荷状态SOC,即使在不同电池电荷状态SOC初值下,也能较快地收敛到真值。
本发明的另一目的是设计一种量测噪声方差阵相匹配的电池电荷状态估计***,其包括微处理器,其含有实现量测噪声方差阵相匹配的电池电荷状态估计方法的对应的软件模块,可嵌入到使用动力电池的装备中,实现动力电池的实时电荷状态在线估计与显示。
本发明设计的量测噪声方差阵相匹配的电池电荷状态估计方法包括如下步骤:
Ⅰ、建立电池模型及参数辨识
Ⅰ-1、建立电池模型
目前,电动汽车电池模型主要有:电化学模型、神经网络模型、阻抗模型、等效电路模型等。本发明根据电池的工作原理,采用联合模型,相对其他的模型,联合模型结构简单、辨识容易,仿真精度较好,具有通用性。其数学表达式:
其中Uk,xk和Ik分别是当前k时刻电池端电压,SOC和电流。电池放电过程Ik为正值,充电过程Ik为负值;a0,a1,a2,a3和a4是通过最小二乘法得到的相关参数;R0是欧姆内阻。ηk是k时刻库伦效率,Cn是电池的标称容量,△t表示采样间隔。ωk是k时刻***过程噪声,νk是k时刻量测噪声。采样得到Uk、Ik。
Ⅰ-2、模型参数辨识
模型参数辨识方法采用最小二乘法得:
c=[a0,a1,a2,a3,a4,R0]T=(MTM)-1MTU (3)其中,U=[U1,U2,...,UN]T是端电压序列,M=[M1,M2,...,MN]T是递归矩阵,
通过最小二乘法离线估计得到模型参数a0,a1,a2,a3,a4和R0,以此作为方差阵相匹配的鲁棒自适应无迹卡尔曼滤波(RAUKF)算法的电池模型基础。
Ⅱ、自适应无迹卡尔曼滤波(AUKF)
自适应无迹卡尔曼滤波估计SOC过程如下:
Ⅱ-1、初始状态统计特性为:
x0是SOC初始值,是SOC初始值的平均值,P0初始误差协方差阵。
Ⅱ-2、SOC的估计值的过程
SOC的估计值
参数更新
是k时刻自协方差阵,是k时刻互协方差阵。Kk是k时刻的滤波增益,zk是k时刻电池端电压的真实值,是k时刻的SOC值。
其中参数:
①
②
③
④
其中,λ=α2(n+κ)-n,α=1e-3,α为正值的比例缩放因子,通过调整α的取值来调节Sigma点与的距离,同时调整α使高阶项的影响达到最小;α的取值范围为0≤α≤1。当***为弱非线性时,α取较大的正值,α≥0.8当***非线性程度严重时,α取一个非常小的正值,α≤1×10-3,以避免采样点非局域效应的影响。
n为状态空间维数,n=1,
κ为比例参数,选择κ=0,
对于高斯函数分布的情况,最佳选择是β=2,
是k-1时刻SOC值,ξ0,k-1是k-1时刻比例对称修正采样的中心点,ξ1,k-1和ξ2,k-1分别是k-1时刻比例对称修正采样相对于中心点的对称点,Pk-1是k-1时刻的误差方差阵,W0 m是一阶统计特性的中心点权值,Wi m是一阶统计特性的对称点权值,W0 c是二阶统计特性的中心点权值,Wi c是二阶统计特性的对称点权值。
库伦效率ηk-1=1。
ξi,k-1是k-1时刻中心点或者对称点的SOC值,γi,k|k-1是基于前一时刻(k-1时刻)状态的当前k时刻的预测的中心点或者对称点的SOC值,ηk-1是k-1时刻的库伦效率,Cn是电池的标称容量,△t表示采样间隔,Ik-1是k-1时刻电流,是基于前一时刻(k-1时刻)状态下k时刻的预测的SOC值,Pk|k-1是基于前一时刻(k-1时刻)状态下k时刻预测的误差方差阵,Qk-1是k-1时刻过程噪声方差阵。ξ0,k|k-1是基于前一时刻(k-1时刻)状态的当前k时刻预测比例对称修正采样的中心点,ξ1,k|k-1和ξ2,k|k-1分别是基于前一时刻(k-1时刻)状态的当前k时刻的预测比例对称修正采样相对于中心点的对称点,χi,k|k-1是基于前一时刻(k-1时刻)状态的当前k时刻的预测中心或者是对称点的电压值,ξi,k|k-1是基于前一时刻(k-1时刻)状态的当前k时刻的预测中心点或者对称点的SOC值,Ik是k时刻电流值,R0是欧姆电阻,a0,a1,a2,a3,a4是模型参数,是一步预测的电压值,Rk是k时刻量测噪声方差阵,是k时刻自协方差阵,是k时刻互协方差阵。
uk是电池模型在k时刻的电压残差,Fk是k时刻电压残差方差阵的近似。Qk是k时刻过程噪声方差阵。
Ⅲ、方差阵相匹配的鲁棒自适应无迹卡尔曼滤波(Robust AdaptiveUnscented Kalman Filter,RAUKF)
卡尔曼滤波的状态估计结果受到许多参数的影响,例如量测噪声方差阵R、过程噪声方差阵Q及滤波增益K。主要调整量测噪声方差阵R,过程噪声方差阵Q的无迹卡尔曼滤波属于自适应无迹卡尔曼滤波(Adaptive UnscentedKalman Filter,AUKF);部分调整滤波增益的无迹自适应卡尔曼滤波属于鲁棒自适应无迹卡尔曼滤波(RAUKF)。
在公式(1)-(8)的自适应卡尔曼滤波算法估计SOC过程基础上,本发明结合鲁棒思想,提出了实际残差方差阵与估计的理论残差方差阵有效匹配方法。
鲁棒卡尔曼滤波中的***信息
定义***实际残差方差阵Ck:
其中L为所选取的平滑窗口,根据经验法,取L=16;***理论残差方差阵为两者的比值为:
其中tr(·)表示矩阵的迹。
当描述***噪声统计特性准确时,即***量测噪声的实际值与给定值一致,用于描述两者的匹配程度的残差方差实际值的迹与估计的理论值的迹之比值θ为1。
当***量测噪声的实际值与给定值不一致时,θ偏离1。为了使卡尔曼滤波稳定,基于由虚拟的量测噪声方差阵和量测噪声方差阵R之比构成的调整因子εk,结合Rk本身自适应调节,对估计的理论残差方差阵进行调整,从而达到调整增益K的目的。
令θ=1,即进行反推计算,过程如下:
得到调整因子εk
其中
为虚拟的量测噪声方差阵。
以θ=1为目标使得估计的理论残差方差阵能够反映实际残差方差阵,实现其方差的完全匹配。
Ⅳ、滤波增益自适应调整
AUKF算法中,量测噪声方差阵R是定值,认为主要受传感器偏差影响,量测噪声统计特性可根据传感器统计特性得知,但当传感器存在偏差时,量测噪声方差阵R将会比真实R值要小。AUKF算法认为R是固定,未考虑到R本身的变化,只是通过大于1的εk系数,以εkR形式,修正R。
本发明以量测噪声方差阵R本身自适应调整为前提,从新的角度,分析量测噪声方差阵R的所具有特殊性。量测噪声方差阵R的具体分析公式Rk的适应算法,
由此式得知,R不仅与传感器有关,还与状态变量有关。
通过公式(13)可知,虚拟的量测噪声方差阵可能为正值,或者是负值,对应的调整因子εk可能出现小于0、小于1或者是大于1的情况.所以当R本身自适应调整时,调整因子εk情况复杂,不局限于RUKF算法仅考虑的调整因子εk大于1的情况。
εk大于1时,把εk大于1带入公式(13)得知,***实际残差方差阵Ck大于***估计的理论残差方差阵为即θ大于1。也就是说,在R固定情况下,以εkR形式,修正对应量测噪声方差阵,只是考虑的θ的一种特殊情况,即对R只进行了部分修正。
θ反映的是实际残差方差阵与估计的理论残差方差阵匹配情况,完全匹配时,比值θ应为1,但实际中,***实际残差方差阵Ck与***估计的理论残差方差阵不相等,即比值θ可能出现θ>1或者0<θ<1。
εk<0时,算法初始阶段,初始值与真实值偏差大且噪声统计特性不准确,导致θ严重偏离1,为趋于0的小数值,此时εk可能为负值,εk不能作为调节因子直接调节Rk来对估计的理论残差方差阵进行调整,故使用自适应调整因子代替εk,取修正Rk,以避免对滤波增益过度调节所带来的协方差负定性问题,δ为衰减因子,0.3<δ<0.7。
当εk>0时,取自适应调整因子修正Rk,对估计的理论残差方差阵进行调整,使估计的理论的量测噪声统计特性动态跟踪实际的量测噪声特性。
本发明以实际的量测噪声统计特性为目标,通过反向推算,求取自适应调整因子并与Rk本身自适应调整相结合,以修正估计的理论计算得到量测噪声统计特性,实现实时跟踪实际的量测噪声特性。具体实现过程如下:
自适应调整因子带入公式(13)调整滤波增益K,得到新的滤波增益NK:
Pk是k时刻的误差协方差,NKk T是NKk的转置。
求得k时刻SOC估计值
当量测噪声水平发生变化时,比值θ为θ>1或者0<θ<1,与只调整比值θ大于1的情况不同,方差阵相匹配的RAUKF算法可以调整比值θ不为1的任何情况,提高了算法对随机量测噪声变化的鲁棒性,实现实际残差方差阵与估计的理论残差方差阵的完全匹配。
本发明量测噪声方差阵相匹配的电池电荷状态估计***,包括微控制器及其连接的显示器,所述微控制器是嵌入式微控制器,动力电池输出端接有电压传感器和电流传感器。电压传感器和电流传感器经模数转换模块连接嵌入微控制器,嵌入微控制器含有模型参数辨识模块、调整因子计算模块、滤波增益自适应调整、鲁棒自适应无迹卡尔曼滤波SOC估计模块。嵌入微控制器连接显示器,还接有CAN(控制器局域网络Controller Area Network)总线接口和/或RS232接口。本***与动力电池一起,嵌入在使用动力电池的设备中,在一个采样周期内完成电压、电流的采集、电池模型参数辨识及SOC在线估计,所得SOC结果在显示器上显示或直接传送到该设备的控制器局域网络。
与现有技术相比,本发明量测噪声方差阵相匹配的电池电荷状态估计方法与***的优点为:提高电池电荷状态估计的鲁棒性,将量测噪声方差阵R本身自适应调整、实际残差方差阵与估计的理论残差方差阵相匹配相结合,与过程噪声方差阵Q相比,观察到量测噪声方差阵R对电池电荷状态估计具有更大的影响,对R本身进行自适应调整,以方差阵的完全匹配--残差方差实际值的迹与估计的理论值的迹相等为目标,分析了基于虚拟的量测噪声方差阵和量测噪声方差阵R之比的调整因子εk出现的各种情况,进行分类,动态地调整滤波增益,保证估计的理论残差方差阵与实际残差方差阵的完全匹配。与无迹卡尔曼滤波(AUKF)、现有的鲁棒无迹卡尔曼滤波(RAUKF)比较,本发明的方差阵相匹配的鲁棒无迹卡尔曼滤波(RAUKF)方法能准确估计电池电荷状态,即使在不同电池电荷状态初值下,也能较快地收敛到真值,具有良好的鲁棒性和收敛性,对噪声有较强的抑制作用。本发明方法所得电池电荷状态估计值的精度提高达0.64%
附图说明
图1为本量测噪声方差阵相匹配的电池电荷状态估计方法实施例的流程示意图;
图2为本量测噪声方差阵相匹配的电池电荷状态估计***实施例的结构示意图;
图3为本量测噪声方差阵相匹配的电池电荷状态估计方法实施例所得SOC曲线、无迹卡尔曼滤波(AUKF)和现有的鲁棒无迹卡尔曼滤波(RAUKF)所得SOC曲线以及实际SOC曲线的比较图
图4为与实际SOC值比较,本量测噪声方差阵相匹配的电池电荷状态估计方法实施例所得SOC的误差曲线、无迹卡尔曼滤波(AUKF)和现有的鲁棒无迹卡尔曼滤波(RAUKF)所得SOC的误差曲线的比较图。
具体实施方式
量测噪声方差阵相匹配的电池电荷状态估计方法实施例
本量测噪声方差阵相匹配的电池电荷状态估计方法实施例的流程如图1所示,包括如下步骤:
Ⅰ、建立电池模型及参数辨识
Ⅰ-1、建立电池模型
本例采用联合模型,其数学表达式:
其中Uk,xk和Ik分别是当前k时刻电池端电压,SOC和电流。电池放电过程Ik为正值,充电过程Ik为负值;a0,a1,a2,a3和a4是通过最小二乘法得到的相关参数;R0是欧姆内阻。ηk是k时刻库伦效率,Cn是电池的标称容量,△t表示采样间隔。ωk是k时刻***过程噪声,νk是k时刻量测噪声。采样得到Uk、Ik。
Ⅰ-2、模型参数辨识
模型参数辨识方法采用最小二乘法得:
c=[a0,a1,a2,a3,a4,R0]T=(MTM)-1MTU (3)其中,U=[U1,U2,...,UN]T是端电压序列,M=[M1,M2,...,MN]T是递归矩阵,
Ⅱ、自适应无迹卡尔曼滤波(AUKF)
自适应无迹卡尔曼滤波估计SOC过程如下:
Ⅱ-1、初始状态统计特性为:
x0是SOC初始值,是SOC初始值的平均值,P0初始误差协方差阵。
Ⅱ-2、SOC的估计值的过程
SOC的估计值
参数更新
是k时刻自协方差阵,是k时刻互协方差阵。Kk是k时刻的滤波增益,zk是k时刻电池端电压的真实值,是k时刻的SOC值。
Ⅲ、方差阵相匹配的鲁棒自适应无迹卡尔曼滤波
鲁棒卡尔曼滤波中的***信息
定义***实际残差方差阵Ck:
其中L为所选取的平滑窗口,根据经验法,取L=16;***估计的理论残差方差阵为两者的比值为:
其中tr(·)表示矩阵的迹。
当描述***噪声统计特性准确时,即***量测噪声的实际值与给定值一致,用于描述两者的匹配程度的残差方差实际值的迹与估计的理论值的迹之比值θ为1。
当***量测噪声的实际值与给定值不一致时,θ偏离1。调整因子εk结合Rk本身自适应调节,对估计的理论残差方差阵进行调整。
得到调整因子εk
其中
为虚拟的量测噪声方差阵。
Ⅳ、滤波增益自适应调整
本例量测噪声方差阵R进行本身自适应调整,Rk的适应算法如下,
εk<0时,取自适应调整因子修正Rk、对估计的理论残差方差阵进行调整以避免对滤波增益Kk过度调节所带来的协方差负定性问题,δ为衰减因子,本例即自适应调整因子
当εk>0时,取自适应调整因子修正Rk、对估计的理论残差方差阵进行调整,使估计的理论量测噪声统计特性动态跟踪实际的量测噪声特性。
求取自适应调整因子并与R本身自适应调整相结合,以修正Rk理论计算得到量测噪声统计特性,具体如下:
调整因子带入公式(13)调整滤波增益K,得到新的滤波增益NKk:
Pk是k时刻的误差协方差,NKk T是NKk的转置。
求得k时刻SOC估计值
量测噪声方差阵相匹配的电池电荷状态估计***实施例
本量测噪声方差阵相匹配的电池电荷状态估计***实施例如图2所示,嵌入式微控制器连接显示器,动力电池输出端接有电压传感器和电流传感器。电压传感器和电流传感器经模数转换模块连接嵌入微控制器,嵌入微控制器含有模型参数辨识模块、调整因子计算模块、滤波增益自适应调整、鲁棒自适应无迹卡尔曼滤波SOC估计模块。嵌入微控制器还接有CAN总线接口和RS232接口。本***与动力电池一起,嵌入在使用动力电池的设备中,在一个采样周期内完成电压、电流的采集、电池模型参数辨识及SOC在线估计,所得SOC结果在显示器上显示或直接传送到该设备的控制器局域网络。
用上述本量测噪声方差阵相匹配的电池电荷状态估计***实施例、按本量测噪声方差阵相匹配的电池电荷状态估计方法对某型号动力电池的电荷状态SOC的全量程变化进行十小时以上的实验得到电池的电荷状态估计值,同时用无迹卡尔曼滤波(AUKF)方法和现有的鲁棒无迹卡尔曼滤波(RAUKF)方法对相同型号的动力电池的电荷状态进行估计,传统实验所得的该型号动力电池的SOC作为SOC实际值。不同方法所得SOC的估计值曲线如图3所示,其中,纵坐标为SOC值,横坐标为时间,大图时间单位为秒(s),局部放大图中时间单位为10-4秒。图中实线为实际SOC曲线,点虚线为AUKF估计SOC曲线,短横虚线为RAUKF估计SOC曲线,点划线为本实施例估计SOC曲线。图中可清楚看到本实施例所得SOC曲线点划线与实际SOC曲线实线最为接近。当本实施例电荷状态初值与实际值相差较大的情况下,也能较快地收敛到真值。
图4所示为与实际SOC值比较,本量测噪声方差阵相匹配的电池电荷状态估计方法实施例所得SOC的误差曲线、无迹卡尔曼滤波(AUKF)和现有的鲁棒无迹卡尔曼滤波(RAUKF)所得SOC的误差曲线,其中,纵坐标为误差值,横坐标为时间,大图时间单位为秒(s),局部放大图中时间单位为10-4秒。点虚线为AUKF估计SOC误差曲线,短横虚线为RAUKF估计SOC误差曲线,点划线为本实施例估计SOC误差曲线。从图中清楚看到本实施例的点划线在实验过程中平稳且基本处于最低位置。
上述实施例,仅为对本发明的目的、技术方案和有益效果进一步详细说明的具体个例,本发明并非限定于此。凡在本发明的公开的范围之内所做的任何修改、等同替换、改进等,均包含在本发明的保护范围之内。
Claims (3)
1.量测噪声方差阵相匹配的电池电荷状态估计方法,其特征在于包括如下步骤:
Ⅰ、建立电池模型及参数辨识
Ⅰ-1、建立电池模型
采用联合模型,其数学表达式:
其中Uk,xk和Ik分别是当前k时刻电池端电压,SOC和电流;电池放电过程Ik为正值,充电过程Ik为负值;a0,a1,a2,a3和a4是通过最小二乘法得到的相关参数;R0是欧姆内阻;ηk是k时刻库伦效率,Cn是电池的标称容量,△t表示采样间隔;ωk是k时刻***过程噪声,νk是k时刻量测噪声;采样得到Uk、Ik;
Ⅰ-2、模型参数辨识
模型参数辨识方法采用最小二乘法得:
c=[a0,a1,a2,a3,a4,R0]T=(MTM)-1MTU (3)
其中,U=[U1,U2,...,UN]T是端电压序列,M=[M1,M2,...,MN]T是递归矩阵,
Ⅱ、自适应无迹卡尔曼滤波
自适应无迹卡尔曼滤波估计SOC过程如下:
Ⅱ-1、初始状态统计特性为:
x0是SOC初始值,是SOC初始值的平均值,P0为初始误差协方差阵;
Ⅱ-2、SOC的估计值的过程
SOC的估计值
参数更新
是k时刻自协方差阵,是k时刻互协方差阵,Kk是k时刻的滤波增益,zk是k时刻电池端电压的真实值,是k时刻的SOC值;
Ⅲ、方差阵相匹配的鲁棒自适应无迹卡尔曼滤波
鲁棒卡尔曼滤波中的***信息
定义***实际残差方差阵Ck:
其中L为所选取的平滑窗口,根据经验法,取L=16;***估计的理论残差方差阵为两者的比值为:
其中tr(·)表示矩阵的迹;
当描述***噪声统计特性准确时,即***量测噪声的实际值与给定值一致,用于描述两者的匹配程度的残差方差实际值的迹与估计的理论值的迹之比值θ为1;
当***量测噪声的实际值与给定值不一致时,θ偏离1;调整因子εk结合Rk本身自适应调节,对估计的理论残差方差阵进行调整;
令θ=1,即Ck代入:
得到调整因子εk
其中
为虚拟的量测噪声方差阵;
Ⅳ、滤波增益自适应调整
本例量测噪声方差阵R进行本身自适应调整,Rk的适应算法如下,
εk<0时,取自适应调整因子修正Rk,δ为衰减因子;
当εk为正值时,取自适应调整因子修正Rk;
求取自适应调整因子并与R本身自适应调整相结合,以修正理论计算得到量测噪声统计特性,具体如下:
自适应调整因子带入公式(13)调整滤波增益,得到新的滤波增益NKk:
Pk是k时刻的误差协方差,NKk T是NKk的转置;
求得k时刻SOC估计值
2.根据权利要求1所述的量测噪声方差阵相匹配的电池电荷状态估计方法,其特征在于:
所述步骤Ⅳ中衰减因子0.3<δ<0.7。
3.根据权利要求1或2所述的量测噪声方差阵相匹配的电池电荷状态估计方法设计的量测噪声方差阵相匹配的电池电荷状态估计***,包括微控制器及其连接的显示器,其特征在于:
所述微控制器是嵌入式微控制器,动力电池输出端接有电压传感器和电流传感器;电压传感器和电流传感器经模数转换模块连接嵌入微控制器,嵌入微控制器含有模型参数辨识模块、调整因子计算模块、滤波增益自适应调整、鲁棒自适应无迹卡尔曼滤波SOC估计模块,还接有CAN总线接口和/或RS232接口;本***与动力电池一起,嵌入在使用动力电池的设备中,在一个采样周期内完成电压、电流的采集、电池模型参数辨识及SOC在线估计,所得SOC结果在显示器上显示或直接传送到该设备的控制器局域网络。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510147926.8A CN104714188B (zh) | 2015-03-31 | 2015-03-31 | 量测噪声方差阵相匹配的电池电荷状态估计方法与*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510147926.8A CN104714188B (zh) | 2015-03-31 | 2015-03-31 | 量测噪声方差阵相匹配的电池电荷状态估计方法与*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104714188A true CN104714188A (zh) | 2015-06-17 |
CN104714188B CN104714188B (zh) | 2017-05-24 |
Family
ID=53413679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510147926.8A Active CN104714188B (zh) | 2015-03-31 | 2015-03-31 | 量测噪声方差阵相匹配的电池电荷状态估计方法与*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104714188B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977545A (zh) * | 2015-08-03 | 2015-10-14 | 电子科技大学 | 一种动力电池的荷电状态估计方法及*** |
CN106646259A (zh) * | 2016-12-26 | 2017-05-10 | 上海思致汽车工程技术有限公司 | 一种电池功率衰减程度的检测方法 |
CN106786702A (zh) * | 2017-01-19 | 2017-05-31 | 上海电力学院 | 用于混合储能***的全状态模型预测能量调度方法 |
CN107505578A (zh) * | 2017-10-16 | 2017-12-22 | 西南交通大学 | 一种锂电池电荷状态估计的方法 |
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN109239604A (zh) * | 2018-10-24 | 2019-01-18 | 江苏理工学院 | 一种基于状态检测机制的无迹卡尔曼滤波车载电池soc估计方法 |
CN109459705A (zh) * | 2018-10-24 | 2019-03-12 | 江苏理工学院 | 一种抗野值鲁棒无迹卡尔曼滤波的动力电池soc估计方法 |
CN109884548A (zh) * | 2019-03-04 | 2019-06-14 | 武汉科技大学 | 一种基于gasvm-aukf算法的锂电池剩余寿命预测方法 |
CN111505506A (zh) * | 2020-05-15 | 2020-08-07 | 吉林大学 | 一种多尺度卡尔曼滤波与无迹卡尔曼滤波融合的电池soc估算方法 |
CN112470017A (zh) * | 2019-02-07 | 2021-03-09 | 株式会社Lg化学 | 电池管理装置、电池管理方法和电池组 |
CN112910430A (zh) * | 2021-01-15 | 2021-06-04 | 北京格润海泰科技有限公司 | 自动调节射频信号功率增益的控制方法及装置 |
CN113125969A (zh) * | 2020-01-14 | 2021-07-16 | 比亚迪股份有限公司 | 基于aukf的电池数据处理方法、设备和介质 |
CN113626950A (zh) * | 2021-08-11 | 2021-11-09 | 广东皓耘科技有限公司 | 一种打捆机设计方法、装置、存储介质及电子设备 |
CN114509947A (zh) * | 2022-02-18 | 2022-05-17 | 淮阴工学院 | 基于鲁棒卡尔曼滤波器的质子交换膜燃料电池故障检测方法 |
CN116449221A (zh) * | 2023-06-14 | 2023-07-18 | 浙江天能新材料有限公司 | 锂电池荷电状态预测方法、装置、设备及存储介质 |
CN117420447A (zh) * | 2023-12-18 | 2024-01-19 | 四川华泰电气股份有限公司 | 考虑噪声偏差补偿的锂电池soc估算方法、***和电子装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101604005A (zh) * | 2009-06-29 | 2009-12-16 | 杭州电子科技大学 | 一种基于组合采样点卡尔曼滤波的电池剩余电量估计方法 |
WO2012098770A1 (ja) * | 2011-01-18 | 2012-07-26 | カルソニックカンセイ株式会社 | バッテリの充電率推定装置 |
CN102645637A (zh) * | 2012-04-25 | 2012-08-22 | 中航锂电(洛阳)有限公司 | 一种均衡电池的soc估算方法 |
CN103472398A (zh) * | 2013-08-19 | 2013-12-25 | 南京航空航天大学 | 基于扩展卡尔曼粒子滤波算法的动力电池soc估计方法 |
CN103529398A (zh) * | 2013-10-28 | 2014-01-22 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的锂离子电池soc在线估计方法 |
CN103744026A (zh) * | 2013-12-19 | 2014-04-23 | 广西科技大学 | 基于自适应无迹卡尔曼滤波的蓄电池荷电状态估算方法 |
US20140257726A1 (en) * | 2011-10-07 | 2014-09-11 | Keio University | Apparatus and method for battery state of charge estimation |
-
2015
- 2015-03-31 CN CN201510147926.8A patent/CN104714188B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101604005A (zh) * | 2009-06-29 | 2009-12-16 | 杭州电子科技大学 | 一种基于组合采样点卡尔曼滤波的电池剩余电量估计方法 |
WO2012098770A1 (ja) * | 2011-01-18 | 2012-07-26 | カルソニックカンセイ株式会社 | バッテリの充電率推定装置 |
US20140257726A1 (en) * | 2011-10-07 | 2014-09-11 | Keio University | Apparatus and method for battery state of charge estimation |
CN102645637A (zh) * | 2012-04-25 | 2012-08-22 | 中航锂电(洛阳)有限公司 | 一种均衡电池的soc估算方法 |
CN103472398A (zh) * | 2013-08-19 | 2013-12-25 | 南京航空航天大学 | 基于扩展卡尔曼粒子滤波算法的动力电池soc估计方法 |
CN103529398A (zh) * | 2013-10-28 | 2014-01-22 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的锂离子电池soc在线估计方法 |
CN103744026A (zh) * | 2013-12-19 | 2014-04-23 | 广西科技大学 | 基于自适应无迹卡尔曼滤波的蓄电池荷电状态估算方法 |
Non-Patent Citations (1)
Title |
---|
魏克新等: "基于自适应无迹卡尔曼滤波算法的锂离子动力电池状态估计", 《中国电机工程学报》 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977545A (zh) * | 2015-08-03 | 2015-10-14 | 电子科技大学 | 一种动力电池的荷电状态估计方法及*** |
CN106646259A (zh) * | 2016-12-26 | 2017-05-10 | 上海思致汽车工程技术有限公司 | 一种电池功率衰减程度的检测方法 |
CN106646259B (zh) * | 2016-12-26 | 2017-11-17 | 上海思致汽车工程技术有限公司 | 一种电池功率衰减程度的检测方法 |
CN106786702B (zh) * | 2017-01-19 | 2019-10-18 | 上海电力学院 | 用于混合储能***的全状态模型预测能量调度方法 |
CN106786702A (zh) * | 2017-01-19 | 2017-05-31 | 上海电力学院 | 用于混合储能***的全状态模型预测能量调度方法 |
CN107505578A (zh) * | 2017-10-16 | 2017-12-22 | 西南交通大学 | 一种锂电池电荷状态估计的方法 |
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN108490472B (zh) * | 2018-01-29 | 2021-12-03 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN109239604A (zh) * | 2018-10-24 | 2019-01-18 | 江苏理工学院 | 一种基于状态检测机制的无迹卡尔曼滤波车载电池soc估计方法 |
CN109459705A (zh) * | 2018-10-24 | 2019-03-12 | 江苏理工学院 | 一种抗野值鲁棒无迹卡尔曼滤波的动力电池soc估计方法 |
CN112470017A (zh) * | 2019-02-07 | 2021-03-09 | 株式会社Lg化学 | 电池管理装置、电池管理方法和电池组 |
US11923710B2 (en) | 2019-02-07 | 2024-03-05 | Lg Energy Solution, Ltd. | Battery management apparatus, battery management method and battery pack |
CN112470017B (zh) * | 2019-02-07 | 2023-12-01 | 株式会社Lg新能源 | 电池管理装置、电池管理方法和电池组 |
CN109884548A (zh) * | 2019-03-04 | 2019-06-14 | 武汉科技大学 | 一种基于gasvm-aukf算法的锂电池剩余寿命预测方法 |
CN113125969B (zh) * | 2020-01-14 | 2022-07-15 | 比亚迪股份有限公司 | 基于aukf的电池数据处理方法、设备和介质 |
CN113125969A (zh) * | 2020-01-14 | 2021-07-16 | 比亚迪股份有限公司 | 基于aukf的电池数据处理方法、设备和介质 |
CN111505506A (zh) * | 2020-05-15 | 2020-08-07 | 吉林大学 | 一种多尺度卡尔曼滤波与无迹卡尔曼滤波融合的电池soc估算方法 |
CN112910430A (zh) * | 2021-01-15 | 2021-06-04 | 北京格润海泰科技有限公司 | 自动调节射频信号功率增益的控制方法及装置 |
CN112910430B (zh) * | 2021-01-15 | 2024-03-12 | 北京格润海泰科技有限公司 | 自动调节射频信号功率增益的控制方法及装置 |
CN113626950A (zh) * | 2021-08-11 | 2021-11-09 | 广东皓耘科技有限公司 | 一种打捆机设计方法、装置、存储介质及电子设备 |
CN113626950B (zh) * | 2021-08-11 | 2024-03-19 | 广东皓耘科技有限公司 | 一种打捆机设计方法、装置、存储介质及电子设备 |
CN114509947A (zh) * | 2022-02-18 | 2022-05-17 | 淮阴工学院 | 基于鲁棒卡尔曼滤波器的质子交换膜燃料电池故障检测方法 |
CN114509947B (zh) * | 2022-02-18 | 2024-07-16 | 淮阴工学院 | 基于鲁棒卡尔曼滤波器的质子交换膜燃料电池故障检测方法 |
CN116449221A (zh) * | 2023-06-14 | 2023-07-18 | 浙江天能新材料有限公司 | 锂电池荷电状态预测方法、装置、设备及存储介质 |
CN116449221B (zh) * | 2023-06-14 | 2023-09-29 | 浙江天能新材料有限公司 | 锂电池荷电状态预测方法、装置、设备及存储介质 |
CN117420447A (zh) * | 2023-12-18 | 2024-01-19 | 四川华泰电气股份有限公司 | 考虑噪声偏差补偿的锂电池soc估算方法、***和电子装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104714188B (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104714188A (zh) | 量测噪声方差阵相匹配的电池电荷状态估计方法与*** | |
Lin et al. | State of charge estimation with the adaptive unscented Kalman filter based on an accurate equivalent circuit model | |
Li et al. | A combination Kalman filter approach for State of Charge estimation of lithium-ion battery considering model uncertainty | |
CN109669134A (zh) | 一种基于卡尔曼滤波法的soc的估算方法 | |
Li et al. | A combination state of charge estimation method for ternary polymer lithium battery considering temperature influence | |
CN102998629B (zh) | 一种锂电池荷电状态的估计方法 | |
CN105717460B (zh) | 一种基于非线性观测器的动力电池soc估计方法和*** | |
Fleischer et al. | Adaptive on-line state-of-available-power prediction of lithium-ion batteries | |
Wang et al. | A novel energy management strategy for the ternary lithium batteries based on the dynamic equivalent circuit modeling and differential Kalman filtering under time-varying conditions | |
CN105301509A (zh) | 锂离子电池荷电状态、健康状态与功率状态的联合估计方法 | |
CN111060824B (zh) | 一种基于模型与算法分步融合的荷电状态估计方法 | |
CN110703114B (zh) | 一种基于电-热-神经网络耦合模型的动力电池soc和sot联合状态估计方法 | |
CN113625174B (zh) | 一种锂离子电池soc与容量联合估计方法 | |
CN101604005A (zh) | 一种基于组合采样点卡尔曼滤波的电池剩余电量估计方法 | |
CN111781503B (zh) | 一种锂离子储能电池soc在线估算方法 | |
Ouyang et al. | A novel state of charge estimation method for lithium-ion batteries based on bias compensation | |
CN111983457A (zh) | 一种基于lstm神经网络的电池组soh估算方法 | |
Li et al. | A novel state estimation approach based on adaptive unscented Kalman filter for electric vehicles | |
CN105223512A (zh) | 基于电池特性的动态矫正剩余电量的方法 | |
CN106154168A (zh) | 数据驱动的动力电池荷电状态估计方法 | |
CN112710955A (zh) | 一种用于提高电池容量估计精度的算法 | |
CN113805075A (zh) | 一种基于bcrls-ukf的锂电池荷电状态估计方法 | |
CN111060822B (zh) | 一种基于模型切换及融合的荷电状态估计方法 | |
CN116047308A (zh) | 一种基于pid控制与dekf的锂电池soc估算方法 | |
Liu et al. | Dynamic adaptive square-root unscented Kalman filter and rectangular window recursive least square method for the accurate state of charge estimation of lithium-ion batteries |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20150617 Assignee: NANNING YONGBAO SOLAR ENERGY CO.,LTD. Assignor: GUILIN University OF ELECTRONIC TECHNOLOGY Contract record no.: X2023980046705 Denomination of invention: A Method and System for Estimating Battery Charge State by Matching Measurement Noise Variance Matrix Granted publication date: 20170524 License type: Common License Record date: 20231109 |