CN104569835B - 一种估计电动汽车的动力电池的荷电状态的方法 - Google Patents
一种估计电动汽车的动力电池的荷电状态的方法 Download PDFInfo
- Publication number
- CN104569835B CN104569835B CN201410784128.1A CN201410784128A CN104569835B CN 104569835 B CN104569835 B CN 104569835B CN 201410784128 A CN201410784128 A CN 201410784128A CN 104569835 B CN104569835 B CN 104569835B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- power battery
- state
- msubsup
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 92
- 239000011159 matrix material Substances 0.000 claims abstract description 85
- 238000007599 discharging Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 24
- 230000010287 polarization Effects 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000012546 transfer Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- QHGJSLXSVXVKHZ-UHFFFAOYSA-N dilithium;dioxido(dioxo)manganese Chemical compound [Li+].[Li+].[O-][Mn]([O-])(=O)=O QHGJSLXSVXVKHZ-UHFFFAOYSA-N 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007772 electrode material Substances 0.000 description 1
- 239000003792 electrolyte Substances 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Landscapes
- Secondary Cells (AREA)
Abstract
本发明涉及电动汽车的动力电池管理***。为提高估计精度,本发明提出一种估计电动汽车的动力电池的荷电状态的方法,采集动力电池的端电压、充放电电流及电池表面温度;根据其开路电压与荷电状态的对应关系解析得出开路电压模型,根据动力电池的动态电压***输入矩阵、参数矩阵以及输出矩阵,建立其动态电压模型;对状态观测器进行时间更新和测量更新,得到***状态的和***状态估计误差的协方差矩阵在k时刻的先验估计值和后验估计值,循环上述更新操作,直至估计完成,得到动力电池的荷电状态SoC。采用该估计方法可将荷电状态SoC的估计精度提高到5%以内,并在动力电池的荷电状态SoC的设定值偏离其参考值时,将荷电状态SoC迅速收敛至其参考值。
Description
技术领域
本发明涉及电动汽车的动力电池管理***,尤其涉及动力电池管理***估计动力电池的荷电状态的方法。
背景技术
在电动汽车中,其动力电池的荷电状态(State-of-Charge,简称SoC,用符号z表示)作为对电动汽车的能量进行管理的决策因素之一,对于优化电动汽车的能量管理、提高动力电池容量和能量的利用率、避免动力电池过充电和/或过放电、提高动力电池的使用安全性、延长动力电池的使用寿命均有重要作用。但是,SoC是动力电池的隐含状态量,难以直接通过测试或计算得出,只能通过估计方法估计得出。
目前,国内外使用的估计动力电池的荷电状态SoC的方法有如下三类:
一、利用动力电池表征参数的测量值估计其荷电状态SoC的方法
下面,以动力电池的当前剩余容量和开路电压Uoc作为动力电池的表征参数为例,说明该估计方法存在的缺陷。
在以动力电池的当前剩余容量为表征参数估计该动力电池的SoC时,需对动力电池进行恒流放电,从而推导计算得出该动力电池的SoC。由于恒流放电过程通常需要持续好几个小时,耗时较长,故该估计方法在电动汽车运行过程中难以实现。
在以动力电池的开路电压Uoc为表征参数估计该动力电池的SoC时,需先将该动力电池静置至少3个小时,才能够测量出其开路电压Uoc。由于在测量开路电压Uoc前所需的静置时间过长,故该估计方法在电动汽车的运行过程中也难以实现。
综上可见,在电动汽车运行过程中,难以利用其表征参数实时估计出其动力电池的SoC,且该类估计方法的估计精度在5%-8%之间,精度较低。
二、利用安时积分法估计动力电池的荷电状态SoC的方法
安时积分法又称为库仑计数法,在测出动力电池的荷电状态SoC的初始值SoC0和充放电电流iL后,根据式(1)积分得出该动力电池当前的SoC:
其中,Cn为该动力电池的额定容量,η为该动力电池的充放电效率,t为该动力电池的充放电时间。
由式(1)可知,在使用该估计方法估计动力电池的SoC时,必须精确地测量出该动力电池的SoC的初始值SoC0、充放电电流iL以及额定容量Cn。但是,在实际操作中,动力电池的SoC的初始值SoC0及额定容量Cn的测量精度会因动力电池的温度、充放电倍率的变化及电池老化等原因而降低;动力电池的充放电电流iL的测量精度易受测量用的电流传感器的测量精度的影响,而电流传感器的测量精度易因受到噪声、温度漂移及其他未知的随机扰动的干扰而降低。故在进行积分计算的过程中,计算结果容易因上述参数的测量精度低而形成误差,且该误差在长期积累后,越来越大,进而导致计算得到的SoC的估计值偏离其实际值5%-10%,估计精度较低。
三、基于动力电池模型估计动力电池的SoC的方法
该估计方法利用动力电池的表征参数比如开路电压、阻抗、极化特性等与其荷电状态SoC之间的映射关系建立动力电池模型,并根据相应的表征参数值推导计算得出动力电池的荷电状态SoC。但该估计方法采用单一的映射关系建立动力电池模型,所建立的动力电池模型固定,导致动力电池的荷电状态SoC的估计精度因动力电池的循环工况、所处环境温度及老化状态的变化而降低,通常在5%-10%之间,估计精度较低,
综上可见,现有的估计动力电池的荷电状态SoC的方法,估计精度低。
发明内容
为提高电动汽车的动力电池的荷电状态SoC的估计精度,本发明提出一种估计电动汽车的动力电池的荷电状态的方法,该方法包括如下步骤:
步骤1、在所述动力电池充放电的过程中,对该动力电池的端电压Ut、充放电电流iL及电池表面温度Tem进行采样,采样时间间隔为Δt;
步骤2、建立所述动力电池的开路电压模型和动态电压模型:
根据所述动力电池的开路电压Uoc与荷电状态SoC、电池表面温度Tem三者之间的对应关系解析得出该动力电池的开路电压模型Uoc(z),且
Uoc(z)=α0+α1z+α2z2+α3z3+α4/z+α5ln z+α6ln(1-z)
其中,
z为所述动力电池的荷电状态SoC的表示符号,
αj(j=0,1,2,……,6)为所述动力电池的开路电压模型的拟合系数,且该拟合系数与电池表面温度Tem相关;
根据所述动力电池的动态电压***输入矩阵参数矩阵θn,k以及输出矩阵yn,k,建立所述动力电池的动态电压模型,
其中,
EL,k为所述动力电池的动态电压EL在k时刻的值,EL,k-1为所述动力电池的动态电压EL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,且EL,k=Ut,k-Uoc(zk),
iL,k为所述动力电池的充放电电流iL在k时刻的值,iL,k-1为所述动力电池的充放电电流iL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,
为所述动态电压模型的模型参数,
n为所述动力电池的等效电路中RC网络的阶数,且n为正整数,
并采用最小二乘法辨识出所述动态电压模型的模型参数,得出所述动力电池的动态电压模型;
步骤3、估计所述动力电池的荷电状态SoC:
建立所述动力电池的状态方程和观测方程,
其中,
xk为所述动力电池在k时刻的***估计状态,且xk=[UD,k,zk]T,UD,k为所述动力电池在k时刻的极化电压,zk为所述动力电池在k时刻的荷电状态,
uk-1为所述动力电池在k-1时刻的***输入信息,
ωk为所述动力电池在k时刻的***状态白噪声,其均值为零,协方差为Qk,
Yk为所述动力电池在k时刻的***观测矩阵,
Uoc,k为所述动力电池在k时刻的开路电压,
ρk为所述动力电池在k时刻的***参数白噪声,其均值为零,协方差为Rk,
对状态观测器中的X,P,Q,R进行初始化设置得到X0、P0、Q0、R0,其中,
X0为所述动力电池的***状态X的初始值,
P0为所述动力电池的***状态估计误差的协方差矩阵P的初始值,
Q0为所述动力电池的过程噪声Q的初始值,
R0为所述动力电池的观测噪声的协方差R的初始值;
对所述状态观测器进行时间更新,且更新的时间长度为一个采样时间间隔Δt,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的先验估计值和
对所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和
循环上述更新操作,得到所述动力电池的荷电状态SoC。
采用该估计方法估计电动汽车的动力电池的荷电状态SoC时,对动力电池的开路电压和动态电压进行分别建模,增强了开路电压对动力电池的荷电状态SoC的估计误差的校正能力;根据动力电池管理***采集的数据实时在线地对其动态电压模型的模型参数进行修正,不仅提高了动力电池的荷电状态SoC的估计精度,还可以在电动汽车运行过程中对其动力电池的荷电状态SoC进行实时估计;根据需要融合不同类型的状态观测器对动力电池的荷电状态SoC进行估计,具有很强的扩展性和开放性,实用性强;估计过程中,考虑到噪声对估计结果的影响,并在对噪声进行抑制;从而在动力电池的荷电状态SoC的初值偏离其参考值时,将荷电状态SoC迅速收敛至其参考值,将荷电状态SoC的估计精度提高到5%以内。
优选地,在所述步骤1中,采集数据时,采样时间间隔Δt为均值,以提高动力电池的荷电状态SoC的估计精度。
优选地,在所述步骤1中,采样得到的端电压Ut、充放电电流iL和电池表面温度Tem分别存储在电压存储单元、电流存储单元和温度存储单元中。
优选地,在所述步骤2中,建立所述动力电池的动态电压模型时,所述动力电池的等效电路中的RC网络的阶数1≤n≤5,以减小运算量,降低估计成本。
优选地,在所述步骤2中,所述最小二乘法为基于遗忘因子的最小二乘法,以实现对动态电压模型的模型参数的实时在线辨识,提高模型参数的辨识精度,进而提高动力电池的荷电状态SoC的估计精度。
优选地,在所述步骤3中,所述状态观测器为基于卡尔曼滤波器或其衍生滤波器的状态观测器。
进一步地,
当所述卡尔曼滤波器的衍生滤波器为扩展卡尔曼滤波器时,
所述状态观测器进行时间更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的先验估计值和且
所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
卡尔曼增益矩阵:
***状态修正:
***状态估计误差的协方差矩阵更新:
其中,Ck为所述动力电池的观测方程在k时刻的雅可比矩阵;
当所述卡尔曼滤波器的衍生滤波器为自适应扩展卡尔曼滤波器时,
所述状态观测器进行时间更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的先验估计值和且
所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
新息矩阵:ek=Ut,k-(Uoc(z)-UD,k-iL,kRi)
卡尔曼增益矩阵:
自适应匹配:
噪声协方差:
***状态修正:
***状态估计误差的协方差更新:
其中,Ck为所述的动力电池的观测方程在k时刻的雅可比矩阵。
在使用本发明估计方法对动力电池的荷电状态SoC进行估计时,对动力电池的开路电压和动态电压进行分别建模,增强了开路电压对动力电池的荷电状态SoC的估计误差的校正能力;根据动力电池管理***采集的数据实时在线地对其动态电压模型的模型参数进行修正,不仅提高了动力电池的荷电状态SoC的估计精度,还可以在电动汽车运行过程中对其动力电池的荷电状态SoC进行实时估计;根据需要融合不同类型的状态观测器对动力电池的荷电状态SoC进行估计,具有很强的扩展性和开放性,实用性强;估计过程中,考虑到噪声对估计误差的影响,并在对噪声进行抑制,以提高估计精度。在采样时,采样时间间隔Δt为均值,以避免动力电池的荷电状态SoC的估计精度因采集到的数据分布不均匀而降低。在建立动力电池的动态电压模型时,该动力电池的等效电路中的RC网络的阶数不超过5,以减小建模的运算量,从而降低建模成本。采用基于遗忘因子的最小二乘法对动力电池的动态模型的模型参数进行辨识,可实现实时在线辨识,提高模型参数的辨识精度,进而提高动力电池的荷电状态SoC的估计精度。在经过上述一系列的改进后,本发明估计方法的估计精度提高到5%以内,甚至1%以内,同时实现对动力电池的荷电状态SoC的实时估计。
本发明还提出一种应用上述任意一种估计方法对电动汽车的动力电池的荷电状态进行估计的动力电池管理***,以对电动汽车的车载动力电池的荷电状态SoC进行实时估计。
附图说明
图1为本发明估计电动汽车的动力电池的荷电状态SoC的流程图;
图2为动力电池含有n阶RC网络时的等效电路图;
图3为基于最小二乘法在线辨识动力电池的动态电压模型的模型参数的流程图;
图4为电动汽车的动力电池进行循环工况试验得到的测试结果,其中,图4(a)为动力电池的充放电电流随时间变化的曲线,图4(b)为动力电池的端电压随时间变化的曲线,图4(c)为动力电池的荷电状态SoC随时间变化的曲线;
图5为试验用电动汽车的动力电池的开路电压与荷电状态SOC、电池表面温度三者之间的对应关系图;
图6为本发明采用基于扩展卡尔曼滤波器的状态观测器估计电动汽车的动力电池的荷电状态SoC的流程图;
图7为本发明采用基于扩展卡尔曼滤波器的状态观测器估计电动汽车的动力电池的荷电状态SoC的仿真结果,其中,图7(a)为动力电池的荷电状态SoC的估计结果与其参考值的对比曲线,图7(b)为动力电池的荷电状态SoC的估计结果与其参考值的差的变化曲线,即误差曲线,图7(c)为图7(a)中在0-300s时间内的局部放大图,图7(d)为图7(b)的放大图;
图8为本发明采用基于自适应扩展卡尔曼滤波器的状态观测器估计电动汽车的动力电池的荷电状态SoC的流程图;
图9为本发明采用自适应扩展卡尔曼滤波器的状态观测器估计电动汽车的动力电池的荷电状态SoC的仿真结果,其中,图9(a)为动力电池的荷电状态SoC的估计结果与其参考值的对比曲线,图9(b)为动力电池的荷电状态SoC的估计结果与其参考值的差的变化曲线,即误差曲线,图9(c)为图9(a)中在0-300s时间内的局部放大图,图9(d)为图9(b)的放大图。
具体实施方式
下面结合图1-9对本发明估计电动汽车的动力电池的荷电状态的方法进行详细说明。
如图1所示,从电动汽车的动力电池管理***(BATTERY MANAGEMENT SYSTEM,简称BMS)中采集实时数据,建立该动力电池的开路电压模型和动态电压模型,并通过对状态观测器进行时间更新和测量更新得出该动力电池的***估计状态,且在进行测量更新时对该动力电池的***估计状态进行修正,从而得出动力电池的荷电状态SoC。
步骤1:采集动力电池的端电压Ut、充放电电流iL及电池表面温度Tem。
从动力电池的管理***(BMS)检测到的关于动力电池的实时数据中采样得到动力电池的端电压Ut、充放电电流iL及电池表面温度Tem,并分别存储到动力电池管理***中的电压存储单元、电流存储单元和温度存储单元中,且采样时间间隔为Δt。优选地,在进行采样时,进行实时在线采样,这样,可实现对动力电池的荷电状态SoC进行实时在线估计。优选地,采样时间间隔Δt为均值,这样,可避免动力电池的荷电状态SoC的估计精度因采样数据分布不均而降低。
步骤2:建立动力电池的开路电压模型和动态电压模型。
由图2可知,该动力电池的等效电路由电压源-OCV、欧姆内阻-Ri和RC网络三部分组成,其中,电压源-OCV为动力电池的开路电压Uoc;欧姆内阻-Ri表示动力电池中电极材料、电解液、隔膜电阻及其他零件的接触电阻;RC网络使用极化内阻RDn和极化电容CDn(n为自然数,表示RC网络的阶数)来描述动力电池的动态特性,该动态特性包括动力电池的极化特性和扩散效应。由基尔霍夫定律可得,
Ut=Uoc-(UD1+UD2+……+UDn)-iLRi (2)
其中,UDn为动力电池的极化电压。
由此可见,可将该动力电池的端电压Ut分为开路电压Uoc和动态电压EL两部分,即Ut=Uoc+EL,故EL=Ut-Uoc。
由于动力电池的荷电状态SoC与其开路电压Uoc之间的映射关系呈单调递增,故选用开路电压Uoc作为估计动力电池SoC的反馈调节量,对动力电池的SoC的估计结果进行修正,以提高估计精度。
在建立动力电池的开路电压模型时,可根据动力电池出厂时配备的其开路电压Uoc与其荷电状态SoC、电池表面温度Tem三者之间的对应关系图表或通过充放电测试得到的三者之间的对应关系解析出该动力电池的开路电压模型Uoc(z),如式(3)所示,
Uoc(z)=α0+α1z+α2z2+α3z3+α4/z+α5ln z+α6ln(1-z) (3)
其中,αj(j=0,1,2,……,6)为动力电池的开路电压模型的拟合系数,且该拟合系数与电池表面温度Tem相关。
在建立动力电池的动态电压模型时,由拉普拉斯变换将式(2)变换为动力电池在频域下的状态方程,如式(4)所示:
其中,iL(s)为充放电电流iL的频域形式,Ut(s)为端电压Ut的频域形式,Uoc(s)为开路电压Uoc的频域形式,s为频域符号,n为动力电池的等效电路中RC网络的阶数。
动态电压模型的传递函数为:
即:
当n=1时,即动力电池的等效电路中RC网络的阶数为1,此时,其动态电压模型的传递函数如式(7)所示:
由于EL=Ut-Uoc,故:
采用如式(9)所示的双线性变换法则将式(7)所示的传递函数映射到Z平面上,可得到该传递函数基于Z平面的方程如式(10)所示:
并定义动态电压模型的模型参数与其系数af的关系为:
则,动态电压模型的模型参数为:
同时简化式(4)可得:
EL,k=a1EL,k-1+a2iL,k+a3iL,k-1 (13)
其中,
EL,k为动力电池的动态电压EL在k时刻的值,EL,k-1为动力电池的动态电压EL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,
iL,k为动力电池的充放电电流iL在k时刻的值,iL,k-1为动力电池的充放电电流iL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt。
由此可得动力电池的动态电压模型在k时刻的输入矩阵参数矩阵θ1,k和输出矩阵y1,k:
当n=2时,即动力电池的等效电路中RC网络的阶数为2,此时,式(4)可简化为:
其传递函数如式(16)所示:
基于式(9)所示的双线性变化法则,可将传递函数映射到Z平面并进行离散化操作得到:
其中,bf(f=1,2,……,5)为与动态电压模型参数相关的系数。
基于式(13)的推导结果,式(15)可转化为:
EL,k=b1EL,k-1+b2EL,k-2+b3iL,k+b4iL,k-1+b5iL,k-2 (18)
故可定义动力电池动态电压模型的输入矩阵为:
定义动力电池动态电压的参数矩阵θ2,k为:
θ2,k=[b1 b2 b3 b4 b5]T (20)
故式(18)可简化为:
相应的,动态电压模型的模型参数与其系数bf的关系式为:
以此类推,当n=n时,动力电池基于频域的输出电压与输入电流之间的数学关系如式(23)所示:
基于式(9)所示的双线性变化法则,可将式(23)化简,得到:
其中,为动态电压模型的模型参数相关的系数。
同理可得到动力电池的动态电压模型的输入矩阵和参数矩阵θn,k为:
式(24)简化为:
经试验计算可知,当n>60时,在单位采样时间间隔Δt内,动力电池的荷电状态SoC对时间的变化率逐渐接近2%,不再可以近似等于零,且随着等效电路中的RC网络的阶数n增大,该等效电路的复杂度也相应增加,不利于动力电池管理***的运算,故本发明优选RC网络阶数1≤n≤5的等效电路。
在确定动力电池的动态电压模型中RC网络的阶数后,采用最小二乘法辨识出动力电池的动态电压模型的模型参数,从而建立动力电池的动态电压模型。
下面,以动力电池含有一阶RC网络的等效电路的动态电压模型为例,具体说明如图3所示的采用该辨识方法辨识动态电压模型的模型参数的过程:
首先,定义辨识***输出量为yk,且
其中,δLs,k为均值为零的白噪声。
接着,利用采集到的输入矩阵和测量得到的输出矩阵yk,根据式(14)得出动态电压模型的在k时刻的模型参数矩阵θk的先验估计值。
然后,根据式(28)分别计算出***状态估计值的卡尔曼增益矩阵KLs,k及***状态估计误差的协方差矩阵PLs,k,
其中,
μ为遗忘因子,且0<μ<1;
I为单位矩阵。
接着,对该动力电池的动态电压模型的模型参数矩阵θk的先验估计值进行修正,得到模型参数矩阵θk的后验估计值且
最后,进行时间更新,并循环上述辨识操作,得出动力电池的动态电压模型的模型参数矩阵θk,进而得出动力电池的动态电压模型的模型参数Ri、RD和CD。
步骤3:估计动力电池的荷电状态SoC。
建立动力电池的状态方程和观测方程如式(30)所示,
其中,
xk为动力电池在k时刻的***估计状态,且xk=[UD,k,zk],UD,k为动力电池在k时刻的极化电压,zk为动力电池在k时刻的荷电状态,
uk-1为动力电池在k-1时刻的***输入矩阵,
Yk为动力电池在k时刻的***观测矩阵,
ωk为动力电池的***状态白噪声,其均值为零,
ρk为动力电池的***参数白噪声,其均值为零,
当动力电池的动态电压模型的等效电路中含有一阶RC网络时,
当动力电池的动态电压模型的等效电路中含有二阶RC网络时,
依次类推,当动力电池的动态电压模型的等效电路中含有n阶RC网络时,
对状态观测器中的X,P,Q,R进行初始化设置得到X0,P0,Q0,R0,其中,
X0为动力电池的***状态X的初始值,
P0为动力电池的***状态估计误差的协方差矩阵P的初始值,
Q0为动力电池的过程噪声Q的初始值,
R0为动力电池的观测噪声的协方差矩阵R的初始值。
对状态观测器进行时间更新,且更新的时间长度为一个采样时间间隔Δt,将动力电池的***状态X和***状态估计误差的协方差矩阵P从前一时刻(k-1)+推算到当前时刻k-,得到动力电池的***状态X和***状态估计误差的协方差矩阵在k时刻的先验估计值和Pk -;
对状态观测器进行测量更新,用动力电池的***观测矩阵Yk对该动力电池的***状态X的先验估计值和***状态估计误差的协方差矩阵P的先验估计值进行修正,得到该动力电池的***状态X和***状态估计误差的协方差矩阵P在时刻k的后验估计值和
循环上述更新操作,直至估计完成,估计得到动力电池的荷电状态SoC。
下面,采用容量为35Ah,电压为3.7V锰酸锂锂离子动力电池(LiMn2O4)作为实验对象,对该动力电池进行循环工况测试,且环境温度为40℃,并采集估计该动力电池的荷电状态SoC所需的端电压Ut,充放电电流iL以及参考用荷电状态SoC,如图4所示,以验证本发明在估计电动汽车上的动力电池的荷电状态SoC时相对于现有的估计方法存在的优势。
根据如图5所示的该动力电池的开路电压与荷电状态SoC、电池表面温度Tem三者之间的对应关系解析出该动力电池的开路电压模型,且该开路电压模型的拟合系数αj(j=0,1,2,……,6)如表1所示,其中,Tem∈[10,45]。
表1
下面,以基于卡尔曼滤波器或其衍生滤波器的状态观测器为例说明采用状态观测器融合动力电池的模型对动力电池的荷电状态SoC进行估计的过程。
Eg.1:采用基于扩展卡尔曼滤波器的状态观测器融合动力电池的模型对该动力电池的荷电状态SoC进行估计,其估计流程如图6所示。
首先,对基于扩展卡尔曼滤波器的状态观测器中的X,P,Q,R进行初始化设置得到X0,P0,Q0,R0,其中,
X0为动力电池的***状态X的初始值,
P0为动力电池的***状态估计误差的协方差矩阵P的初始值,
Q0为动力电池的过程噪声Q的初始值,
R0为动力电池的观测噪声的协方差矩阵R的初始值。
接着,对基于扩展卡尔曼滤波器的状态观测器进行时间更新即进行先验估计,且更新的时间长度为一个采样时间间隔Δt,将动力电池的***状态X和***状态估计误差的协方差矩阵P从前一时刻(k-1)+推算到当前时刻k-,得到动力电池的***状态X和***状态估计误差的协方差矩阵P在k时刻的先验估计值和且
***状态预估:
***状态估计误差的协方差矩阵预估:
然后,对基于该扩展卡尔曼滤波器的状态观测器进行测量更新即进行后验估计,用动力电池的***观测矩阵Yk对该动力电池的***状态X的先验估计值和***状态估计误差的协方差矩阵P的先验估计值进行修正,得到该动力电池的***状态X和***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
卡尔曼增益矩阵:
***状态修正:
***状态估计误差协方差矩阵更新:
其中,Ck为估计过程中电动汽车的动力电池的观测方程在k时刻的雅可比矩阵;
循环上述更新操作,直至估计完成。
采用该估计方法对动力电池的荷电状态进行仿真估计,估计仿真结果如图7所示。可见:该估计方法估计得到的动力电池的荷电状态SoC的最大估计误差在5%以内,相对于现有的估计方法,估计精度较高;当动力电池的荷电状态SoC的设定值偏离其参考值时,该估计方法能够将动力电池的荷电状态SoC迅速收敛至其参考值。
Eg.2:采用基于自适应扩展卡尔曼滤波器的状态观测器融合动力电池的模型对该动力电池的荷电状态SoC进行估计,其估计流程如图8所示。
首先,对基于自适应扩展卡尔曼滤波器的状态观测器中的X,P,Q,R进行初始化设置得到X0,P0,Q0,R0,其中,
X0为动力电池的***状态X的初始值,
P0为动力电池的***状态估计误差的协方差矩阵P的初始值,
Q0为动力电池的过程噪声Q的初始值,
R0为动力电池的观测噪声的协方差矩阵R的初始值。
接着,对基于扩展卡尔曼滤波器的状态观测器进行时间更新即进行先验估计,且更新的时间长度为一个采样时间间隔Δt,将动力电池的***状态X和***状态估计误差的协方差矩阵P从前一时刻(k-1)+推算到当前时刻k-,得到***状态X和***状态估计误差的协方差矩阵P在k时刻的先验估计值和且
***状态预估:
***状态估计误差的协方差矩阵预估:
然后,对基于该自适应扩展卡尔曼滤波器的状态观测器进行测量更新即进行后验估计,用动力电池的***观测矩阵Yk对该动力电池的***状态X的先验估计值和***状态估计误差的协方差矩阵P的先验估计值进行修正,得到该动力电池的***状态X和***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
新息矩阵:ek=Ut,k-(Uoc(z)-UD,k-iL,kRi)(38)
卡尔曼增益矩阵:
自适应匹配:
噪声协方差:
***状态修正:
***状态估计误差的协方差矩阵更新:
其中,Ck为估计过程中电动汽车的动力电池的观测方程在k时刻的雅可比矩阵;
循环上述更新操作,直至估计完成。
采用该估计方法对动力电池的荷电状态进行仿真估计,估计仿真结果如图9所示。可见:该估计方法估计动力电池的荷电状态SoC的最大估计误差在1%以内,其估计精度远高于现有的估计方法的估计精度;当动力电池的荷电状态SoC的设定值偏离其参考值时,该估计方法能够将动力电池的荷电状态SoC迅速收敛至其参考值。
综上可见,相对现有的估计动力电池的荷电状态SoC的方法,本发明所提出的估计动力电池的荷电状态SoC的方法具有如下优势:
(1)该估计方法使得动力电池的荷电状态SoC的估计值的精度可达到5%以内,甚至1%以内,极大地提高了动力电池的荷电状态SoC的估计精度;
(2)当动力电池的荷电状态SoC的设定值偏离其参考值时,该估计方法能够使动力电池的荷电状态SoC迅速收敛至其参考值,进而提高估计精度;
(3)该估计方法对动力电池的开路电压和动态电压进行分别建模,增强了开路电压对动力电池的荷电状态SoC的估计误差的校正能力,提高了估计精度;
(4)该估计方法能够根据动力电池管理***采集的数据实时在线地对其动态电压模型的模型参数进行修正,不仅提高了动力电池的荷电状态SoC的估计精度,还可以在电动汽车运行过程中对其动力电池的荷电状态SoC进行实时估计;
(5)该估计方法能够融合不同类型的状态观测器对动力电池的荷电状态SoC进行估计,比如基于扩展卡尔曼滤波器的状态观测器和基于自适应扩展卡尔曼滤波器的状态观测器,具有很强的扩展性和开放性,实用性强;
(6)该估计方法将动力电池的动态电压模型和状态观测器融合在一起,且在估计时,考虑到噪声对估计结果的影响,并在估计过程中对噪声进行抑制,以提高估计精度。
Claims (8)
1.一种估计电动汽车的动力电池的荷电状态的方法,其特征在于,该方法包括如下步骤:
步骤1、在所述动力电池充放电的过程中,对该动力电池的端电压Ut、充放电电流iL及电池表面温度Tem进行采样,采样时间间隔为Δt;
步骤2、建立所述动力电池的开路电压模型和动态电压模型:
根据所述动力电池的开路电压Uoc与荷电状态SoC、电池表面温度Tem三者之间的对应关系解析得出该动力电池的开路电压模型Uoc(z),且
Uoc(z)=α0+α1z+α2z2+α3z3+α4/z+α5ln z+α6ln(1-z)
其中,
z为所述动力电池的荷电状态SoC的表示符号,
αj(j=0,1,2,……,6)为所述动力电池的开路电压模型的拟合系数,且该拟合系数与电池表面温度Tem相关;
根据所述动力电池的动态电压***输入矩阵参数矩阵θn,k以及输出矩阵yn,k,建立所述动力电池的动态电压模型,
其中,
EL,k为所述动力电池的动态电压EL在k时刻的值,EL,k-1为所述动力电池的动态电压EL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,且EL,k=Ut,k-Uoc(zk),
iL,k为所述动力电池的充放电电流iL在k时刻的值,iL,k-1为所述动力电池的充放电电流iL在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,
为所述动态电压模型的模型参数,
n为所述动力电池的等效电路中RC网络的阶数,且n为正整数,
并采用最小二乘法辨识出所述动态电压模型的模型参数,得出所述动力电池的动态电压模型;
步骤3、估计所述动力电池的荷电状态SoC:
建立所述动力电池的状态方程和观测方程,
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>Ax</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Bu</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>oc</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>L</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,
xk为所述动力电池在k时刻的***估计状态,且xk=[UD,k,zk]T,UD,k为所述动力电池在k时刻的极化电压,zk为所述动力电池在k时刻的荷电状态,
uk-1为所述动力电池在k-1时刻的***输入信息,
ωk为所述动力电池在k时刻的***状态白噪声,其均值为零,协方差为Qk,
Yk为所述动力电池在k时刻的***观测矩阵,
Uoc,k为所述动力电池在k时刻的开路电压,
ρk为所述动力电池在k时刻的***参数白噪声,其均值为零,协方差为Rk,
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open='[' close=']'>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mi>&Delta;t</mi>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>D</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>D</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mi>&Delta;t</mi>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>D</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>D</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mi>&Delta;t</mi>
<mrow>
<msub>
<mi>R</mi>
<mi>Dn</mi>
</msub>
<msub>
<mi>C</mi>
<mi>Dn</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>Dn</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mfrac>
<mi>&Delta;t</mi>
<msub>
<mi>C</mi>
<mi>n</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
对状态观测器中的X,P,Q,R进行初始化设置得到X0、P0、Q0、R0,其中,
X0为所述动力电池的***状态X的初始值,
P0为所述动力电池的***状态估计误差的协方差矩阵P的初始值,
Q0为所述动力电池的过程噪声Q的初始值,
R0为所述动力电池的观测噪声的协方差R的初始值;
对所述状态观测器进行时间更新,且更新的时间长度为一个采样时间间隔Δt,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的先验估计值和
对所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和
循环上述更新操作,得到所述动力电池的荷电状态SoC。
2.根据权利要求1所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,在所述步骤1中,采集数据时,采样时间间隔Δt为均值。
3.根据权利要求1所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,在所述步骤1中,采样得到的端电压Ut、充放电电流iL和电池表面温度Tem分别存储在电压存储单元、电流存储单元和温度存储单元中。
4.根据权利要求1-3中任意一项所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,在所述步骤2中,建立所述动力电池的动态电压模型时,所述动力电池的等效电路中的RC网络的阶数1≤n≤5。
5.根据权利要求4所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,在所述步骤2中,所述最小二乘法为基于遗忘因子的最小二乘法。
6.根据权利要求5所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,在所述步骤3中,所述状态观测器为基于卡尔曼滤波器或其衍生滤波器的状态观测器。
7.根据权利要求6所述的估计电动汽车的动力电池的荷电状态的方法,其特征在于,
当所述卡尔曼滤波器的衍生滤波器为扩展卡尔曼滤波器时,
所述状态观测器进行时间更新,得到所述***状态X和所述***状态估计误差的协方
差矩阵P在k时刻的先验估计值和且
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>Ax</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Bu</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>AP</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
卡尔曼增益矩阵:
<mrow>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<msubsup>
<mi>C</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<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>
</mrow>
***状态修正:
<mrow>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
<mo>+</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>oc</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>L</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
***状态估计误差的协方差矩阵更新:
其中,Ck为所述动力电池的观测方程在k时刻的雅可比矩阵;
当所述卡尔曼滤波器的衍生滤波器为自适应扩展卡尔曼滤波器时,
所述状态观测器进行时间更新,得到所述***状态X和所述***状态估计误差的协方
差矩阵P在k时刻的先验估计值和且
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>Ax</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Bu</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>AP</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
所述状态观测器进行测量更新,得到所述***状态X和所述***状态估计误差的协方差矩阵P在k时刻的后验估计值和且
新息矩阵:ek=Ut,k-(Uoc(z)-UD,k-iL,kRi)
卡尔曼增益矩阵:
<mrow>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<msubsup>
<mi>C</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<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>
</mrow>
自适应匹配:
<mrow>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mi>M</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>e</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
</mrow>
噪声协方差:
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>P</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<msubsup>
<mi>C</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>K</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
***状态修正:
***状态估计误差的协方差更新:
其中,Ck为所述的动力电池的观测方程在k时刻的雅可比矩阵。
8.一种应用权利要求1-7中任意一项所述的估计方法对电动汽车的动力电池的荷电状态进行估计的动力电池管理***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410784128.1A CN104569835B (zh) | 2014-12-16 | 2014-12-16 | 一种估计电动汽车的动力电池的荷电状态的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410784128.1A CN104569835B (zh) | 2014-12-16 | 2014-12-16 | 一种估计电动汽车的动力电池的荷电状态的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104569835A CN104569835A (zh) | 2015-04-29 |
CN104569835B true CN104569835B (zh) | 2017-11-17 |
Family
ID=53086358
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410784128.1A Expired - Fee Related CN104569835B (zh) | 2014-12-16 | 2014-12-16 | 一种估计电动汽车的动力电池的荷电状态的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104569835B (zh) |
Families Citing this family (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044606B (zh) * | 2015-07-01 | 2018-03-02 | 西安交通大学 | 一种基于参数自适应电池模型的soc估计方法 |
CN105449739B (zh) * | 2015-09-30 | 2019-08-23 | 上海凌翼动力科技有限公司 | 基于单体电池的串联电池组荷电状态soc在线估计方法 |
CN105425153B (zh) * | 2015-11-02 | 2017-12-05 | 北京理工大学 | 一种估计电动车辆的动力电池的荷电状态的方法 |
CN105425154B (zh) * | 2015-11-02 | 2018-02-06 | 北京理工大学 | 一种估计电动汽车的动力电池组的荷电状态的方法 |
CN105467328B (zh) * | 2015-12-29 | 2018-07-03 | 哈尔滨工业大学 | 一种锂离子电池荷电状态估计方法 |
CN105842627B (zh) * | 2016-02-01 | 2018-06-01 | 北京理工大学 | 基于数据模型融合的动力电池容量和荷电状态的估计方法 |
CN105607010A (zh) * | 2016-02-02 | 2016-05-25 | 北京理工大学 | 一种估计电动车辆的动力电池的健康状态的方法 |
US10023064B2 (en) * | 2016-03-10 | 2018-07-17 | Ford Global Technologies, Llc | Power capability estimation for vehicle battery systems |
CN106324521B (zh) * | 2016-09-05 | 2018-09-11 | 北京理工大学 | 一种联合估计动力电池***参数与荷电状态的方法 |
CN106249171B (zh) * | 2016-09-18 | 2018-10-26 | 北京理工大学 | 一种用于宽采样间隔的动力电池***辨识和状态估计方法 |
CN110167783B (zh) | 2017-01-09 | 2022-09-20 | 沃尔沃卡车集团 | 一种用于确定电池组的充电状态的方法和装置 |
CN106802394B (zh) * | 2017-02-06 | 2020-05-22 | 清华大学深圳研究生院 | 汽车电池荷电状态的修正方法及装置 |
CN106842060A (zh) * | 2017-03-08 | 2017-06-13 | 深圳市海云图新能源有限公司 | 一种基于动态参数的动力电池soc估算方法及*** |
CN107870305B (zh) * | 2017-12-04 | 2019-10-18 | 浙江大学城市学院 | 基于温度参数的锂离子电池在线参数辨识和soh估计方法 |
CN108445408A (zh) * | 2018-03-20 | 2018-08-24 | 重庆大学 | 一种基于参数估计ocv的全温度soc估计方法 |
CN108318733B (zh) * | 2018-03-26 | 2020-05-12 | 深圳市天鹏建设工程有限公司 | 智慧楼宇物业控制*** |
CN108872865B (zh) * | 2018-05-29 | 2020-12-29 | 太原理工大学 | 一种防滤波发散的锂电池soc估算方法 |
CN111262896A (zh) * | 2018-11-30 | 2020-06-09 | 联合汽车电子有限公司 | 一种网联汽车电池管理*** |
US10950903B2 (en) * | 2018-12-07 | 2021-03-16 | GM Global Technology Operations LLC | Battery state estimation using electrode transient model |
US11034257B2 (en) * | 2019-01-11 | 2021-06-15 | GM Global Technology Operations LLC | Method and system for estimating remaining battery pack energy using cell-group state of charge spread |
CN109839596B (zh) * | 2019-03-25 | 2021-04-16 | 重庆邮电大学 | 基于ud分解的自适应扩展卡尔曼滤波的soc估算方法 |
CN110018419B (zh) * | 2019-03-29 | 2021-09-10 | 深圳猛犸电动科技有限公司 | 锂离子电池的soc与soh联合估计方法、装置及终端设备 |
CN109901083B (zh) * | 2019-04-01 | 2020-05-22 | 北京理工大学 | 一种动力电池ocv曲线在线重构方法 |
CN110261778B (zh) * | 2019-05-27 | 2020-08-21 | 南京理工自动化研究院有限公司 | 一种锂离子电池soc估测算法 |
CN110286332B (zh) * | 2019-06-17 | 2021-06-22 | 南通大学 | 一种基于多新息理论的电动汽车动力电池soc估计方法 |
CN110308394B (zh) * | 2019-06-25 | 2020-11-06 | 南京理工大学 | 一种新型状态观测器的锂离子电池荷电状态估测方法 |
CN110308396B (zh) * | 2019-07-03 | 2020-09-25 | 华人运通(江苏)技术有限公司 | 电池状态监测方法、边缘处理器、***及存储介质 |
CN110361642A (zh) * | 2019-07-11 | 2019-10-22 | 中国科学院电工研究所 | 一种电容器荷电状态的预测方法、装置及电子设备 |
CN110456279A (zh) * | 2019-08-15 | 2019-11-15 | 长安大学 | 一种基于数据驱动模型的动力电池云管理*** |
CN112433154A (zh) * | 2019-08-25 | 2021-03-02 | 南京理工大学 | 基于ffrls和ekf的锂离子电池soc估测算法 |
CN112601968A (zh) * | 2019-08-27 | 2021-04-02 | 淄博火炬能源有限责任公司 | 用于48v轻混汽车锂离子电池的荷电状态在线估算*** |
CN111308351B (zh) * | 2019-10-18 | 2021-07-23 | 南京航空航天大学 | 低温环境动力电池soc估计方法、存储介质和设备 |
CN110716141B (zh) * | 2019-11-27 | 2022-02-08 | 四川长虹电器股份有限公司 | 一种基于扩展卡尔曼滤波的电池荷电状态估计方法 |
CN111146514B (zh) * | 2019-12-19 | 2021-05-07 | 上海派能能源科技股份有限公司 | 锂离子电池模组运行安全性评估预测方法、***及电子设备 |
CN112649745A (zh) * | 2020-12-16 | 2021-04-13 | 北京理工大学 | 一种协同重组ocv的线性化类状态观测器的电池soc估算方法 |
CN114264966A (zh) * | 2021-12-06 | 2022-04-01 | 阳光储能技术有限公司 | 电池的荷电状态的评估方法、装置、终端设备以及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001343437A (ja) * | 2000-06-05 | 2001-12-14 | Matsushita Electric Ind Co Ltd | 蓄電池の寿命判定方法 |
CN101098029A (zh) * | 2006-06-26 | 2008-01-02 | 三星Sdi株式会社 | 估计电池充电状态的方法、电池管理***及其驱动方法 |
CN103439668A (zh) * | 2013-09-05 | 2013-12-11 | 桂林电子科技大学 | 动力锂离子电池的电荷状态估算方法与*** |
CN103901351A (zh) * | 2014-03-18 | 2014-07-02 | 浙江大学城市学院 | 一种基于滑动窗滤波的单体锂离子电池soc估计方法 |
CN104181470A (zh) * | 2014-09-10 | 2014-12-03 | 山东大学 | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 |
-
2014
- 2014-12-16 CN CN201410784128.1A patent/CN104569835B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001343437A (ja) * | 2000-06-05 | 2001-12-14 | Matsushita Electric Ind Co Ltd | 蓄電池の寿命判定方法 |
CN101098029A (zh) * | 2006-06-26 | 2008-01-02 | 三星Sdi株式会社 | 估计电池充电状态的方法、电池管理***及其驱动方法 |
CN103439668A (zh) * | 2013-09-05 | 2013-12-11 | 桂林电子科技大学 | 动力锂离子电池的电荷状态估算方法与*** |
CN103901351A (zh) * | 2014-03-18 | 2014-07-02 | 浙江大学城市学院 | 一种基于滑动窗滤波的单体锂离子电池soc估计方法 |
CN104181470A (zh) * | 2014-09-10 | 2014-12-03 | 山东大学 | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 |
Non-Patent Citations (2)
Title |
---|
A data-driven based adaptive state of charge estimator of lithium-ion polymer battery used in electric vehicles;Rui Xiong et.;《Applied Energy》;20140131;第113卷;1421-1433 * |
Evaluation on State of Charge Estimation of Batteries With Adaptive Extended Kalman Filter by Experiment Approach;Rui Xiong et.;《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》;20130131;第62卷(第1期);108-117 * |
Also Published As
Publication number | Publication date |
---|---|
CN104569835A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104569835B (zh) | 一种估计电动汽车的动力电池的荷电状态的方法 | |
CN105425153B (zh) | 一种估计电动车辆的动力电池的荷电状态的方法 | |
CN108196200B (zh) | 一种锂电池健康和荷电状态的联合模拟评估方法 | |
CN105301509B (zh) | 锂离子电池荷电状态、健康状态与功率状态的联合估计方法 | |
CN111060820B (zh) | 一种基于二阶rc模型的锂电池soc和sop估计方法 | |
CN110068772B (zh) | 基于改进分数阶模型的锂离子电池荷电状态估计方法 | |
CN106716158B (zh) | 电池荷电状态估算方法和装置 | |
CN106054085B (zh) | 一种基于温度用于估计电池soc的方法 | |
CN103020445B (zh) | 一种电动车车载磷酸铁锂电池的soc与soh预测方法 | |
CN104122504B (zh) | 一种电池的soc估算方法 | |
CN106324523B (zh) | 基于离散变结构观测器的锂电池soc估计方法 | |
CN110795851B (zh) | 一种考虑环境温度影响的锂离子电池建模方法 | |
CN101604005B (zh) | 一种基于组合采样点卡尔曼滤波的电池剩余电量估计方法 | |
CN105425154B (zh) | 一种估计电动汽车的动力电池组的荷电状态的方法 | |
CN109633479B (zh) | 基于嵌入式容积卡尔曼滤波的锂电池soc在线估算方法 | |
CN109839599B (zh) | 基于二阶ekf算法的锂离子电池soc估计方法 | |
CN105699910A (zh) | 一种锂电池剩余电量在线估计方法 | |
CN111060824A (zh) | 一种基于模型与算法分步融合的荷电状态估计方法 | |
CN111044906A (zh) | 一种基于极大似然准则的锂离子电池能量状态估算方法 | |
CN112269133B (zh) | 一种基于预充电路模型参数识别的soc估计方法 | |
CN111060823A (zh) | 一种基于dp模型的低温环境下电池sop在线估计方法 | |
CN115656848A (zh) | 一种基于容量修正的锂电池soc估算方法 | |
CN114720881A (zh) | 一种基于改进初值带遗忘因子递推最小二乘法的锂电池参数辨识方法 | |
CN114910796A (zh) | 基于miaukf算法的锂离子电池荷电状态估计方法 | |
CN112816873B (zh) | 一种基于改进电池模型的iekf锂电池soc估算方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171117 Termination date: 20211216 |