CN104181470A - 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 - Google Patents
一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 Download PDFInfo
- Publication number
- CN104181470A CN104181470A CN201410459523.2A CN201410459523A CN104181470A CN 104181470 A CN104181470 A CN 104181470A CN 201410459523 A CN201410459523 A CN 201410459523A CN 104181470 A CN104181470 A CN 104181470A
- Authority
- CN
- China
- Prior art keywords
- battery
- matrix
- soc
- partiald
- battery system
- 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
Landscapes
- Secondary Cells (AREA)
- Tests Of Electric Status Of Batteries (AREA)
Abstract
本发明公开了一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,它具有运算简单与精度高的优点;它包括步骤:1)将整个方法所需要的时间平均分成N个时间段,每一个时间段代表一步,即第k个时间段为第k步;其中k≤N,k与N均为正整数;2)根据基尔霍夫电压和电流定理,建立电池***的模型;得到电池***模型的误差矩阵d(k)与电池***的模型误差分配矩阵G(k);3)对步骤2)所得的先验状态估计方程进行补偿;4)求解第k+1步电池***的后验状态估计方程,得到SOC值;5)根据步骤4)的第k+1步电池***的后验状态估计结果,与电池的SOC真值进行比较,验证非线性预测卡尔曼滤波NPEKF算法的有效性,k=n+1,并转到步骤3),直至第N步执行完。
Description
技术领域
本发明涉及电池的电荷状态SOC估计方法,尤其涉及一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法。
背景技术
车载动力电池作为电动汽车的关键部件,其性能对整车的动力性、经济性和安全性至关重要,是制约电动汽车规模发展的关键因素。为最大限度地发挥动力电池的性能并延长电池的使用寿命,对电池进行有效管理至关重要,而准确获得电池的荷电状态(state of charge,SOC)是电池管理最核心的技术。电池SOC估计是判断电池是否过充过放,是否需要均衡或更换某一单体电池的重要依据。因此,提高电池SOC估计精度,对于提高电池使用效率、延长电池循环寿命以及保障电池安全可靠,有着重要的意义。
目前常用的电池SOC估计方法可以分为以下四类:①基于安时计量的估计方法,主要有放电实验法和Ah计量法;②基于电池表征参数测量值的估计方法,主要有开路电压法和内阻法;③基于经验方程和数学模型的估计方法,主要有线性模型法、神经网络法、模糊逻辑法、支持向量机法和相关向量机法;④基于电池等效电路模型的估计方法,主要有卡尔曼滤波(KF)法、观测器法和粒子滤波法。近几年新兴的基于电池等效电路模型的估计方法因其模型物理意义明确,易于参数辨识,已成为电池研究的热点。其中扩展卡尔曼滤波(EKF)算法备受关注,它适应于电流波动比较剧烈的SOC估计,非常适合应用在电动汽车上。但是EKF有效应用的前提条件是:
1)要求***模型准确已知;
2)过程噪声被假设为协方差已知的零均值高斯白噪声过程。
在基于EKF的电池SOC估计中,电池动态模型中存在的误差被作为过程噪声来处理,而且一般假定为零均值、协方差已知的高斯白噪声。然而,这种假设往往缺少理论根据,不能正确反映电池模型误差的真实特性,从而导致SOC估计精度下降,甚至使滤波器发散。
发明内容
本发明的目的就是为解决上述问题,提出了一种基于非线性预测扩展卡尔曼滤波(NPEKF)的电池SOC估计方法,能够在电池模型存在较大模型误差的情况下准确估计电池SOC,有效地解决存在显著模型误差的非线性滤波问题,与传统基于EKF的SOC估计算法相比较,NPEKF算法具有以下优点:1)过程噪声的统计特性无需事先准确知道(不必满足高斯分布),对于未知的电池模型误差,作为滤波器解的一部分可以实时地得到它的一步预测估计值;2)电池模型误差的结构形式不受限制,可以是线性的,也可以是非线性的;3)省去了由于协方差的求解而带来的运算负担,而且不必计算非线性状态函数与测量函数的雅可比矩阵。
为了实现上述目的,本发明所采用的技术方案为:
一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,包括步骤:
1)将整个方法所需要的时间平均分成N个时间段,每一个时间段代表一步,即第k个时间段为第k步;其中k≤N,k与N均为正整数;
2)根据基尔霍夫电压和电流定理,建立电池***的模型;以电池的电流为输入变量,以SOC值与支路电压为状态变量,得到电池***的先验状态估计方程;以电池两端的电压为输出变量,以电池的电流为输入变量,以SOC值与支路电压为状态变量,得到电池***的测量方程;根据电池***模型的测量方程,得到电池***模型的误差矩阵d(k)与电池***的模型误差分配矩阵G(k);
3)对步骤2)所得的先验状态估计方程进行补偿:利用电池***的模型误差矩阵d(k)与电池***的模型误差分配矩阵G(k)对步骤2)所得的先验状态估计方程进行补偿,得到精准的先验状态估计方程,k=n,n为正整数,n≤N;
4)求解第k+1步电池***的后验状态估计方程,得到SOC值;
5)根据步骤4)得到的第k+1步电池***的后验状态估计结果,与电池的SOC真值进行比较,验证非线性预测卡尔曼滤波NPEKF算法的有效性,k=n+1,并转到步骤3),直至第N步执行完。
所述步骤3)中对先验状态估计方程进行补偿的具体步骤为:
31)对电池***的模型误差分配矩阵G(k)进行初始化;
32)求解电池***的模型误差矩阵d(k);
33)对步骤2)所得的先验状态估计方程进行补偿:
其中为第k+1步的先验状态估计,为关于的函数,xk为状态变量;uk为输入变量。
所述步骤31)还包括对NPEKF滤波器的初始化:
给定状态初始估计值给定电池***的模型误差加权矩阵W、过程激励噪声协方差矩阵Q、测量噪声方差矩阵R和误差协方差矩阵P0的初始值。
所述步骤3)中电池***的模型误差矩阵d(k)的具体求解过程为:
321)求解参数Λ(Δt):
Λ(Δt)是1×1维的对角阵,其值为:
Λ(Δt)=Δt2/2;
322)求解参数
的计算公式如下:
式中,分别是关于的第1阶和第2阶李导数,为输出变量估计关于的函数,为第k+1步的状态估计关于的函数,输入变量uk=ik,ik为第k个时刻点的电池端电流;
323)求解参数
是1×2维的矩阵,其行元素如下:
式中,g1和个g2分别是G(k)中的第1和第2列元素;是先沿后沿g1的李导数;是先沿后沿g2的李导数。
根据上述公式可得到:
式中,Uoc,k为第k步电池开路电压,SOC为电池的电荷状态;
324)求解模型误差矩阵d(k):
其中,R为测量噪声协方差矩阵,W为加权矩阵,Δt=tk+1-tk为采样时间间隔,为第k步的状态估计,为第k步的yk的估计,yk+1为第k+1步的测量真值,yk+1=Ubat,k+1,Ubat,k+1第k+1步的电池端电压。
所述步骤4)求解电池***的后验状态估计方程的具体步骤为:
41)计算第k+1步的增益矩阵Kk+1;
42)求解后验估计
其中,为第k+1步测量值估计值,为,yk+1为第k+1步的测量真值,具体指代为第k+1步电池的端电压Ubat,k+1,uk+1=ik+1,ik+1为第k+1步的电池电流,为第k+1步状态的先验估计。
所述步骤41)中第k+1步的增益矩阵Kk+1的具体求解过程为:
411)计算第k+1步误差协方差矩阵的先验估计
其中,Ak+1为第k+1步的***矩阵,为第k步的误差协方差矩阵的后验估计,Qk为第k步的过程激励噪声协方差矩阵,wk+1表示过程激励噪声,为正态分布的白色噪声:wk+1~N(0,Rk+1);
412)计算第k+1步的增益矩阵Kk+1:
其中,为第k+1步的误差协方差矩阵的先验估计,Hk+1为第k+1步的测量矩阵,vk+1为过程观测噪声,为正态分布vk+1~N(0,Rk+1)。
所述步骤411)第k步的误差协方差矩阵的后验估计的具体求解为:
其中Kk为第k步的增益矩阵,Hk为第k步的测量矩阵,为第k步的误差协方差矩阵的先验估计。
本发明的有益效果为:
1、对于未知的模型误差,该算法能够将其作为滤波器解的一部分,可以实时地得到它的一步预测估计值,并用于修正电池的动态模型,提高了SOC估计精度;
2、模型误差的结构形式不受限制;
3、算法运算简单;
4、有效地解决存在显著动态模型误差情况下的非线性电池***SOC估计问题,其估计误差在4%以内;
5、适用于电池的各种复杂工况。
附图说明
图1为本发明的二阶RC电池等效电路模型图,其中下标d为放电方向的参数;下标c为充电方向的参数。
图2为本发明的***流程图;
图3为本发明具体流程图;
图4为本发明的二阶RC模型在脉冲放电下的端电压及其误差曲线图。
图5为本发明的二阶RC模型在脉冲充电下的端电压及其误差曲线图。
图6为本发明的NPEKF滤波算法的流程图。
具体实施方式
为了更好的了解本发明的技术方案,下面结合附图对本发明作进一步说明。
1、二阶RC模型
利用非线性预测扩展卡尔曼滤波器估计电池SOC,需要精确建立电池模型。搭建电池模型是指应用数学理论尽量全面地去描述实际电池的响应特性和内部特性。所谓响应特性是指电池的端电压与负载电流的对应关系;内部特性是指电池的内部变量欧姆内阻、极化内阻和极化电压与SOC、温度间的关系。
如图1所示为本发明的二阶RC等效电路模型,包括运行时间电路和I-V特性电路。
所述运行时间电路包括电池的自放电电阻Rd、电容CQ和电流控制电流源电路,电阻Rd与电容CQ并联在电流控制电流源的受控源两端,独立电源的一端接地。所述I-V特性电路包括两路支路,每个支路包括三组一个电容与一个电阻并联组成的RC回路,且每个支路的两组RC回路并联有切换开关。所述I-V特性电路包括欧姆内阻R0、电化学极化内阻R1、电化学极化电容C1、浓差极化内阻R2、浓差极化电容C2和电流控制电流源、电压控制电压源电路,其中:电压控制电压源电路的受控源的正极连接两路,一路连接二极管Dd后连接电阻R1d、电阻R2d和电阻Rod后连接电池的正极,一路反接二极管Dc后连接电阻R1c、电阻R2c和电阻Roc后连接电池的正极,电容C1d并联在电阻R1d两端,电容C2d并联在电阻R2d两端,电容C1c并联在电阻R1c两端,电容C2c并联在电阻R2c两端,电压控制电压源电路的受控源正、负极之间的电压为电池开路电压OCV。
所述运行时间电路和I-V特性电路通过一个流控电流源和一个压控电压源建立联系,当对电池进行充放电时,负载电流ibat通过流控电流源对电容CQ进行充放电,从而改变CQ存储的电量,从而表征电池SOC的变化,CQ两端电压OCV也随之变化,进而控制I-V特性电路的压控电压源OCV=f(SOC)随SOC的变化而变化。负载电流ibat对电容CQ进行放电时,当CQ两端电压OCV达到放电终止电压时即可获得电池总的运行时间。
所述电容CQ表示电池的可用容量,CQ=3600·CAh·f1·f2,CAh为用安时为单位的电池容量,f1和f2分别是电池循环寿命和温度的修正因子。
所述电流控制电流源的受控源的电流为电池的端电流ibat,当电池进行充放电时负载电流ibat通过电流控制电流源对电容CQ进行充放电,改变电容CQ中存储的电量,从而表征电池SOC的变化。
所述电流控制电流源的受控源两端的电压为电池开路电压OCV(或记为Uoc),所述开路电压OCV与SOC存在较强的非线性关系,采用一个压控电压源来模拟(如图1所示),具体关系式为:
式中,a1-a5为常数,由实验数据基于最小二乘法辨识得到。
图2-3为本发明的流程图。
如图4-5所示分别为本发明的二阶RC模型在脉冲放电和充电下的端电压及其误差曲线图。从图中可以看出,采用本文提出的变阶数RC等效电路模型能很好地反应电池的脉冲充放电过程,这说明该模型是准确的,可以利用的。在恒流充放电的过程中所产生的误差要比在静止的阶段产生的误差要大一点,这是因为拟合的是零输入零状态响应。总体来说,静置阶段电池的误差几乎为零,充放电阶段电池的误差小于0.02V。
2、二阶RC电池模型的数学方程
以电池放电为例,在二阶RC模型的基础上,结合安时法,根据基尔霍夫电压和电流原理,建立了如下数学模型:
Ubat(t)=UOC(SOC)-ibat(t)R0-U1(t)-U2(t) (2)
式中,Ubat为电池的端电压;R0为欧姆内阻;R1为电化学极化内阻,C1为电化学极化电容;R2为浓差极化电阻,C2为浓差极化电容;Uoc为电池的开路电压;ibat(t)为电池的端电流,以下简记为i。
要利用NPEK计算动力电池的SOC,需要将SOC作为NPEK的状态变量之一。利用Ah法对SOC进行估算的离散化表达式为:
式中,η为库仑系数,可以通过电池充放电试验得到,电池充电时,η=1,放电时,η<1;ik为采样时刻点k处的电流,具有符号性,充电时ik>0,放电时ik<0;C为电池标称容量,单位是Ah,Δt=tk+1-tk是采样间隔。
根据图1所示的动力电池模型,选取RC网络中电容的电压和SOC作为***状态变量,***的状态方程为:
上述方程可记为:
其中,状态变量为:
输入变量为:
uk=ik (8)
***矩阵为:
输入控制矩阵为:
根据图1所示的动力电池模型,可得到电池***的测量方程的离散形式为:
上述(11)式可记为:
式中,随机信号vk+1表示过程观测噪声,假设为正态分布的白色噪声,且vk+1~N(0,Rk+1),Rk+1为第k+1步的观测噪声协方差。
测量矩阵为:
输出变量为:
yk+1=Ubat,k+1 (15)
3、输出估计的泰勒级数近似
对方程(5)中的进行泰勒级数展开,可得到
整理上式,得到
式中,Λ(Δt)是1×1维的对角阵,其值为:
Λ(Δt)=Δt2/2 (18)
的计算公式如下:
式中,是关于的m阶李导数,定义为:
根据上述定义可以得到:
式中,是1×2维的矩阵,其行元素给出如下:
式中,g1和个g2分别是G(k)中的第1和第2列元素;是先沿后沿g1的李导数;是先沿后沿g2的李导数。
根据上述公式可得到:
5、基于NPEKF的SOC估计
一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,如图2所示,它的步骤总结为:
第一步,建立电池数学模型。
以二阶RC电池模型为例,选取RC网络中电容的电压和SOC作为***状态变量,***的状态方程为:
上述方程可记为:
其中,状态变量为:
输入变量为:
uk=ik (28)
***矩阵为:
输入控制矩阵为:
根据图1所示的动力电池模型,可得到电池***的测量方程的离散形式为:
上述(27)式可记为:
式中,随机信号vk+1表示过程观测噪声,假设为正态分布的白色噪声,且vk+1~N(0,Rk+1),Rk+1为第k+1步的观测噪声协方差。
测量矩阵为:
输出变量为:
yk+1=Ubat,k+1 (35)
上式中,Ubat为电池的端电压;R0为欧姆内阻;R1为电化学极化内阻,C1为电化学极化电容;R2为浓差极化电阻,C2为浓差极化电容;U1,U2分别是电容C1和C2的端电压;Uoc为电池的开路电压;ik为k时刻的电池端电流。
第二步,NPEKF滤波器的初始化。
给定状态初始估计值加权矩阵W以及测量噪声方差矩阵R和误差协方差矩阵P0。
状态估计值初始化:
式中,SOC0为估计初值。
误差协方差初始化:
过程激励噪声协方差:
测量噪声协方差:
选取G(k)为:
结合先验知识选取加权矩阵W为:
其中,W为模型误差加权矩阵,可根据经验事先给定。值得说明的是,模型误差加权矩阵W的取值对估计效果的影响很大,但没有固定的选取规则。通过仿真发现,若模型误差是加性的或者初始估计误差过大时,W的值不应过大;若模型误差为乘性的,W的值可以大一些;测量误差一般对W的取值不构成影响。
第三步,计算参量
设已获得了k时刻的状态估计值利用计算出中间参数矩阵 在接到k+1时刻的测量信息y(k+1)后,根据公式(26)预测[k,k+1]间隔内的模型误差项d(k),k=n,n为自然数。
Λ(Δt)是1×1维的对角阵,其值为:
Λ(Δt)=Δt2/2 (42)
的计算公式如下:
式中,是关于的m阶李导数,定义为:
根据上述定义可以得到:
式中,是1×2维的矩阵,其行元素给出如下:
式中,g1和个g2分别是G(k)中的第1和第2列元素;是先沿后沿g1的李导数;是先沿后沿g2的李导数。
根据上述公式可得到:
第四步,预测模型误差d(k)。
把d(k)的预测值代入状态估计方程,将状态估计值传播到k+1时刻,即得到
式中,d(k)为2×1维的模型误差矩阵;G(k)为3×2维的模型误差分配矩阵。-代表先验,^代表估计,为在已知第k+1步以前状态情况下第k+1步的先验状态估计;为已知测量真值yk+1时,第k+1步的后验状态估计。
第六步,误差协方差矩阵预测
式中,Qk+1为第k+1步的过程噪声协方差,随机信号wk+1表示过程激励噪声,假设为正态分布的白色噪声:wk~N(0,Qk)。
第七步,计算增益矩阵Kk+1:
式中,Rk+1为第k+1步的观测噪声协方差,随机信号vk+1表示过程观测噪声,假设为正态分布的白色噪声,且vk+1~N(0,Rk+1)。
第六步,状态修正
式中,yk+1为第k+1步的测量真值,本发明中具体指代电池的端电压Ubat,k+1;第k+1步的测量估计值,即电池的端电压的估计值;uk+1为第k+1步的***输入值,本发明中具体指代电池的端电流ibat,k+1。
第九步,误差协方差矩阵修正Pk+1:
第十步,根据方程(54)得到第k+1步的估计结果,并与SOC真值进行比较,验证NPEKF算法的有效性,k=n+1,并转到第三步。
图6给出了在UDDS循环工况下基于NPEKF的SOC估计结果。NPEKF的初始参数设置为:SOC真值初始值SOCt0=0.9,SOC估计初始值SOCe0=0.8。从实验结果可以看出,SOC估计在初始SOC有较大误差的情况下能迅速向实际值逼近,之后SOC估计值在SOC真值附近波动,最大估计误差小于4%,证明了本发明提出的NPEKF算法的有效性。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,包括步骤:
1)将整个方法所需要的时间平均分成N个时间段,每一个时间段代表一步,即第k个时间段为第k步;其中k≤N,k与N均为正整数;
2)根据基尔霍夫电压和电流定理,建立电池***的模型;以电池的电流为输入变量,以SOC值与支路电压为状态变量,得到电池***的先验状态估计方程;以电池两端的电压为输出变量,以电池的电流为输入变量,以SOC值与支路电压为状态变量,得到电池***的测量方程;根据电池***的模型的测量方程,得到电池***模型的误差矩阵d(k)与电池***的模型误差分配矩阵G(k);
3)对步骤2)所得的先验状态估计方程进行补偿:利用电池***的模型误差矩阵d(k)与电池***的模型误差分配矩阵G(k)对步骤2)所得的先验状态估计方程进行补偿,得到精准的先验状态估计方程,k=n,n为正整数,n≤N;
4)求解第k+1步电池***的后验状态估计方程,得到SOC值;
5)根据步骤4)得到的第k+1步电池***的后验状态估计结果,与电池的SOC真值进行比较,验证非线性预测卡尔曼滤波NPEKF算法的有效性,k=n+1,并转到步骤3),直至第N步执行完。
2.如权利要求1所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤3)中对先验状态估计方程进行补偿的具体步骤为:
31)对电池***的模型误差分配矩阵G(k)进行初始化;
32)求解电池***的模型误差矩阵d(k);
33)对步骤2)所得的先验状态估计方程进行补偿:
其中为第k+1步的先验状态估计,为关于的函数,xk为状态变量;uk为输入变量。
3.如权利要求2所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤31)还包括对NPEKF滤波器的初始化:
给定状态初始估计值给定电池***的模型误差加权矩阵W、过程激励噪声协方差矩阵Q、测量噪声方差矩阵R和误差协方差矩阵P0的初始值。
4.如权利要求书1-2任一所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤3)中电池***的模型误差矩阵d(k)的具体求解过程为:
321)求解参数Λ(Δt):
Λ(Δt)是1×1维的对角阵,其值为:
Λ(Δt)=Δt2/2;
322)求解参数
的计算公式如下:
式中,分别是关于的第1阶和第2阶李导数,为输出变量估计关于的函数,为第k+1步的状态估计关于的函数,输入变量uk=ik,ik为第k个时刻点的电池端电流;
323)求解参数
是1×2维的矩阵,其行元素如下:
式中,g1和个g2分别是G(k)中的第1和第2列元素;是先沿后沿g1的李导数;是先沿后沿g2的李导数;
根据上述公式可得到:
式中,Uoc,k为第k步电池开路电压,SOC为电池的电荷状态;
324)求解模型误差矩阵d(k):
其中,R为测量噪声协方差矩阵,W为加权矩阵,Δt=tk+1-tk为采样时间间隔,为第k步的状态估计,为第k步的yk的估计,yk+1为第k+1步的测量真值,yk+1=Ubat,k+1,Ubat,k+1第k+1步的电池端电压。
5.如权利要求1所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤4)求解电池***的后验状态估计方程的具体步骤为:
41)计算第k+1步的增益矩阵Kk+1;
42)求解后验估计
其中,为第k+1步测量值估计值,yk+1为第k+1步的测量真值,具体指代为第k+1步电池的端电压Ubat,k+1,uk+1=ik+1,ik+1为第k+1步的电池电流,为第k+1步状态的先验估计。
6.如权利要求5所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤41)中第k+1步的增益矩阵Kk+1的具体求解过程为:
411)计算第k+1步误差协方差矩阵的先验估计
其中,Ak+1为第k+1步的***矩阵,为第k步的误差协方差矩阵的后验估计,Qk为第k步的过程激励噪声协方差矩阵,wk+1表示过程激励噪声,为正态分布的白色噪声:wk+1~N(0,Rk+1);
412)计算第k+1步的增益矩阵Kk+1:
其中,为第k+1步的误差协方差矩阵的先验估计,Hk+1为第k+1步的测量矩阵,vk+1为过程观测噪声,为正态分布vk+1~N(0,Rk+1)。
7.如权利要求6所述的一种基于非线性预测扩展卡尔曼滤波的电池SOC估计方法,其特征是,所述步骤411)第k步的误差协方差矩阵的后验估计的具体求解为:
其中Kk为第k步的增益矩阵,Hk为第k步的测量矩阵,为第k步的误差协方差矩阵的先验估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410459523.2A CN104181470B (zh) | 2014-09-10 | 2014-09-10 | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410459523.2A CN104181470B (zh) | 2014-09-10 | 2014-09-10 | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104181470A true CN104181470A (zh) | 2014-12-03 |
CN104181470B CN104181470B (zh) | 2017-04-26 |
Family
ID=51962676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410459523.2A Active CN104181470B (zh) | 2014-09-10 | 2014-09-10 | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104181470B (zh) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535934A (zh) * | 2014-12-31 | 2015-04-22 | 桂林电子科技大学 | 在线前馈补偿的动力电池电荷状态估计方法与*** |
CN104569835A (zh) * | 2014-12-16 | 2015-04-29 | 北京理工大学 | 一种估计电动汽车的动力电池的荷电状态的方法 |
CN105116338A (zh) * | 2015-07-15 | 2015-12-02 | 盐城工学院 | 一种基于soc补偿器的并联型电池***建模方法 |
CN105425154A (zh) * | 2015-11-02 | 2016-03-23 | 北京理工大学 | 一种估计电动汽车的动力电池组的荷电状态的方法 |
CN105842633A (zh) * | 2016-05-30 | 2016-08-10 | 广西大学 | 一种基于灰色扩展卡尔曼的锂离子电池soc估算方法 |
CN105929340A (zh) * | 2016-06-30 | 2016-09-07 | 四川普力科技有限公司 | 一种基于arima估算电池soc的方法 |
CN105974329A (zh) * | 2016-07-22 | 2016-09-28 | 深圳市沃特玛电池有限公司 | 一种估算电池组soh的方法 |
CN106451643A (zh) * | 2016-10-28 | 2017-02-22 | 四川普力科技有限公司 | 电能源管理***及其方法 |
CN106840458A (zh) * | 2017-03-03 | 2017-06-13 | 镇江海姆霍兹传热传动***有限公司 | 基于扩展卡尔曼滤波的多温度传感器融合方法 |
CN106855612A (zh) * | 2017-02-21 | 2017-06-16 | 山东大学 | 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法 |
CN106872901A (zh) * | 2017-02-21 | 2017-06-20 | 山东大学 | KiBaM‑分数阶等效电路综合特征电池模型及参数辨识方法 |
CN106896327A (zh) * | 2017-03-10 | 2017-06-27 | 山东大学 | 分数阶KiBaM‑等效电路综合特征电池模型及其参数辨识方法 |
CN106918787A (zh) * | 2017-03-20 | 2017-07-04 | 国网重庆市电力公司电力科学研究院 | 一种电动汽车锂电池剩余电荷估算方法及装置 |
CN107422269A (zh) * | 2017-06-16 | 2017-12-01 | 上海交通大学 | 一种锂电池在线soc测量方法 |
CN107505574A (zh) * | 2017-08-08 | 2017-12-22 | 中华人民共和国成都出入境检验检疫局 | 基于极大似然带渐消因子的ukf锂电池估算方法 |
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN108646710A (zh) * | 2018-05-10 | 2018-10-12 | 中国民航大学 | 一种基于改进容积粒子滤波的电液联合舵机状态估计方法 |
WO2018188321A1 (zh) * | 2017-04-13 | 2018-10-18 | 绵阳世睿科技有限公司 | 一种增强电池状态估计鲁棒性的方法 |
CN109061505A (zh) * | 2018-08-28 | 2018-12-21 | 淮阴工学院 | 一种锂电池soh的检测方法 |
CN109188293A (zh) * | 2018-11-08 | 2019-01-11 | 武汉理工大学 | 基于新息协方差带渐消因子的ekf锂离子电池soc估算方法 |
CN109613437A (zh) * | 2018-11-27 | 2019-04-12 | 南京航空航天大学 | 一种基于在线模型参数识别的电池荷电状态估计方法 |
CN111474481A (zh) * | 2020-04-13 | 2020-07-31 | 深圳埃瑞斯瓦特新能源有限公司 | 基于扩展卡尔曼滤波算法的电池soc估算方法及装置 |
CN112327182A (zh) * | 2020-08-02 | 2021-02-05 | 西北工业大学 | 基于量测值残差序列的自适应h无穷滤波soc估计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6534954B1 (en) * | 2002-01-10 | 2003-03-18 | Compact Power Inc. | Method and apparatus for a battery state of charge estimator |
CN101629992A (zh) * | 2009-05-27 | 2010-01-20 | 重庆大学 | 磷酸铁锂动力电池剩余容量的估算方法 |
CN102088118A (zh) * | 2010-12-28 | 2011-06-08 | 深圳市航盛电子股份有限公司 | 一种电池管理***、电动车及荷电状态的估算方法 |
CN103620432A (zh) * | 2012-03-16 | 2014-03-05 | 株式会社Lg化学 | 用于估计电池状态的设备和方法 |
CN103675683A (zh) * | 2012-09-02 | 2014-03-26 | 东莞市振华新能源科技有限公司 | 一种锂电池荷电状态(soc)估算方法 |
CN103926538A (zh) * | 2014-05-05 | 2014-07-16 | 山东大学 | 基于aic准则的变阶数rc等效电路模型及实现方法 |
CN103941195A (zh) * | 2014-05-05 | 2014-07-23 | 山东大学 | 基于小模型误差准则扩展卡尔曼滤波的电池soc估计方法 |
-
2014
- 2014-09-10 CN CN201410459523.2A patent/CN104181470B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6534954B1 (en) * | 2002-01-10 | 2003-03-18 | Compact Power Inc. | Method and apparatus for a battery state of charge estimator |
CN101629992A (zh) * | 2009-05-27 | 2010-01-20 | 重庆大学 | 磷酸铁锂动力电池剩余容量的估算方法 |
CN102088118A (zh) * | 2010-12-28 | 2011-06-08 | 深圳市航盛电子股份有限公司 | 一种电池管理***、电动车及荷电状态的估算方法 |
CN103620432A (zh) * | 2012-03-16 | 2014-03-05 | 株式会社Lg化学 | 用于估计电池状态的设备和方法 |
CN103675683A (zh) * | 2012-09-02 | 2014-03-26 | 东莞市振华新能源科技有限公司 | 一种锂电池荷电状态(soc)估算方法 |
CN103926538A (zh) * | 2014-05-05 | 2014-07-16 | 山东大学 | 基于aic准则的变阶数rc等效电路模型及实现方法 |
CN103941195A (zh) * | 2014-05-05 | 2014-07-23 | 山东大学 | 基于小模型误差准则扩展卡尔曼滤波的电池soc估计方法 |
Non-Patent Citations (3)
Title |
---|
卢杰祥: "锂离子电池特性建模与SOC估算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
毛华夫等: "基于卡尔曼滤波修正算法的电池SOC估算", 《电源技术》 * |
袁闪闪: "磷酸铁锂动力电池荷电状态估算方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569835A (zh) * | 2014-12-16 | 2015-04-29 | 北京理工大学 | 一种估计电动汽车的动力电池的荷电状态的方法 |
CN104569835B (zh) * | 2014-12-16 | 2017-11-17 | 北京理工大学 | 一种估计电动汽车的动力电池的荷电状态的方法 |
CN104535934B (zh) * | 2014-12-31 | 2017-07-21 | 桂林电子科技大学 | 在线前馈补偿的动力电池电荷状态估计方法与*** |
CN104535934A (zh) * | 2014-12-31 | 2015-04-22 | 桂林电子科技大学 | 在线前馈补偿的动力电池电荷状态估计方法与*** |
CN105116338A (zh) * | 2015-07-15 | 2015-12-02 | 盐城工学院 | 一种基于soc补偿器的并联型电池***建模方法 |
CN105425154A (zh) * | 2015-11-02 | 2016-03-23 | 北京理工大学 | 一种估计电动汽车的动力电池组的荷电状态的方法 |
CN105425154B (zh) * | 2015-11-02 | 2018-02-06 | 北京理工大学 | 一种估计电动汽车的动力电池组的荷电状态的方法 |
CN105842633A (zh) * | 2016-05-30 | 2016-08-10 | 广西大学 | 一种基于灰色扩展卡尔曼的锂离子电池soc估算方法 |
CN105929340A (zh) * | 2016-06-30 | 2016-09-07 | 四川普力科技有限公司 | 一种基于arima估算电池soc的方法 |
CN105929340B (zh) * | 2016-06-30 | 2019-08-20 | 四川普力科技有限公司 | 一种基于arima估算电池soc的方法 |
CN105974329A (zh) * | 2016-07-22 | 2016-09-28 | 深圳市沃特玛电池有限公司 | 一种估算电池组soh的方法 |
CN106451643A (zh) * | 2016-10-28 | 2017-02-22 | 四川普力科技有限公司 | 电能源管理***及其方法 |
CN106872901A (zh) * | 2017-02-21 | 2017-06-20 | 山东大学 | KiBaM‑分数阶等效电路综合特征电池模型及参数辨识方法 |
CN106872901B (zh) * | 2017-02-21 | 2019-09-24 | 山东大学 | KiBaM-分数阶等效电路综合特征电池模型及参数辨识方法 |
CN106855612A (zh) * | 2017-02-21 | 2017-06-16 | 山东大学 | 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法 |
CN106855612B (zh) * | 2017-02-21 | 2019-09-24 | 山东大学 | 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法 |
CN106840458B (zh) * | 2017-03-03 | 2019-04-05 | 镇江海姆霍兹传热传动***有限公司 | 基于扩展卡尔曼滤波的多温度传感器融合方法 |
CN106840458A (zh) * | 2017-03-03 | 2017-06-13 | 镇江海姆霍兹传热传动***有限公司 | 基于扩展卡尔曼滤波的多温度传感器融合方法 |
CN106896327B (zh) * | 2017-03-10 | 2019-04-16 | 山东大学 | 分数阶KiBaM-等效电路综合特征电池模型及其参数辨识方法 |
CN106896327A (zh) * | 2017-03-10 | 2017-06-27 | 山东大学 | 分数阶KiBaM‑等效电路综合特征电池模型及其参数辨识方法 |
CN106918787A (zh) * | 2017-03-20 | 2017-07-04 | 国网重庆市电力公司电力科学研究院 | 一种电动汽车锂电池剩余电荷估算方法及装置 |
WO2018188321A1 (zh) * | 2017-04-13 | 2018-10-18 | 绵阳世睿科技有限公司 | 一种增强电池状态估计鲁棒性的方法 |
CN107422269B (zh) * | 2017-06-16 | 2020-02-07 | 上海交通大学 | 一种锂电池在线soc测量方法 |
CN107422269A (zh) * | 2017-06-16 | 2017-12-01 | 上海交通大学 | 一种锂电池在线soc测量方法 |
CN107505574A (zh) * | 2017-08-08 | 2017-12-22 | 中华人民共和国成都出入境检验检疫局 | 基于极大似然带渐消因子的ukf锂电池估算方法 |
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN108646710A (zh) * | 2018-05-10 | 2018-10-12 | 中国民航大学 | 一种基于改进容积粒子滤波的电液联合舵机状态估计方法 |
CN109061505A (zh) * | 2018-08-28 | 2018-12-21 | 淮阴工学院 | 一种锂电池soh的检测方法 |
CN109061505B (zh) * | 2018-08-28 | 2021-02-19 | 淮阴工学院 | 一种锂电池soh的检测方法 |
CN109188293A (zh) * | 2018-11-08 | 2019-01-11 | 武汉理工大学 | 基于新息协方差带渐消因子的ekf锂离子电池soc估算方法 |
CN109613437A (zh) * | 2018-11-27 | 2019-04-12 | 南京航空航天大学 | 一种基于在线模型参数识别的电池荷电状态估计方法 |
CN111474481A (zh) * | 2020-04-13 | 2020-07-31 | 深圳埃瑞斯瓦特新能源有限公司 | 基于扩展卡尔曼滤波算法的电池soc估算方法及装置 |
CN111474481B (zh) * | 2020-04-13 | 2022-08-09 | 深圳埃瑞斯瓦特新能源有限公司 | 基于扩展卡尔曼滤波算法的电池soc估算方法及装置 |
CN112327182A (zh) * | 2020-08-02 | 2021-02-05 | 西北工业大学 | 基于量测值残差序列的自适应h无穷滤波soc估计方法 |
CN112327182B (zh) * | 2020-08-02 | 2021-11-16 | 西北工业大学 | 基于量测值残差序列的自适应h无穷滤波soc估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104181470B (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104181470A (zh) | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 | |
Ren et al. | Design and implementation of a battery management system with active charge balance based on the SOC and SOH online estimation | |
CN103926538B (zh) | 基于aic准则的变阶数rc等效电路模型及实现方法 | |
CN103941195B (zh) | 基于小模型误差准则扩展卡尔曼滤波的电池soc估计方法 | |
Aung et al. | State-of-charge estimation of lithium-ion battery using square root spherical unscented Kalman filter (Sqrt-UKFST) in nanosatellite | |
CN103439668B (zh) | 动力锂离子电池的电荷状态估算方法与*** | |
CN103616647B (zh) | 一种用于电动汽车电池管理***的电池剩余电量估计方法 | |
CN102680795B (zh) | 一种二次电池内阻的实时在线估计方法 | |
CN107741569A (zh) | 一种基于分段扩展卡尔曼滤波的锂电池荷电状态的估算方法 | |
CN103632018B (zh) | 一种基于Simscape平台的电池建模方法 | |
CN102645637A (zh) | 一种均衡电池的soc估算方法 | |
CN108008320A (zh) | 一种锂离子电池荷电状态与模型参数自适应联合估计方法 | |
CN103675683A (zh) | 一种锂电池荷电状态(soc)估算方法 | |
CN104392080A (zh) | 一种锂电池分数阶变阶等效电路模型及其辨识方法 | |
CN103529398A (zh) | 基于扩展卡尔曼滤波的锂离子电池soc在线估计方法 | |
CN104122504A (zh) | 一种电池的soc估算方法 | |
Wehbe et al. | Battery equivalent circuits and brief summary of components value determination of lithium ion: A review | |
CN110231567A (zh) | 一种电动汽车soc估算算法 | |
CN106872899A (zh) | 一种基于降维观测器的动力电池soc估计方法 | |
CN109085505A (zh) | 一种动力电池充放电状态估算方法 | |
CN105116346A (zh) | 一种串联型锂离子电池***及其荷电状态估计方法 | |
CN104573294A (zh) | 动力电池的自适应卡尔曼滤波器估计算法 | |
CN105044606A (zh) | 一种基于参数自适应电池模型的soc估计方法 | |
CN105445665A (zh) | 一种利用卡尔曼滤波预估电池荷电状态的方法 | |
CN107783057A (zh) | 一种动力电池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 |