CN111537894B - 一种用于估计锂电池soc和sop的方法 - Google Patents
一种用于估计锂电池soc和sop的方法 Download PDFInfo
- Publication number
- CN111537894B CN111537894B CN202010478880.9A CN202010478880A CN111537894B CN 111537894 B CN111537894 B CN 111537894B CN 202010478880 A CN202010478880 A CN 202010478880A CN 111537894 B CN111537894 B CN 111537894B
- Authority
- CN
- China
- Prior art keywords
- battery
- value
- equation
- current
- soc
- 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.)
- Active
Links
Images
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
-
- 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/385—Arrangements for measuring battery or accumulator variables
-
- 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/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
- G01R31/388—Determining ampere-hour charge capacity or SoC involving voltage measurements
Abstract
本发明公开了一种用于估计锂电池SOC和SOP的方法,通过在不同温度下进行充放电实验,测得不同温度下的电池可用容量,拟合出温度和电池实际放出的可用容量之间的关系;建立反映电池特性的二阶等效电路模型;通过最小二乘法实现对模型参数的辨识,对建立的等效电路模型方程进行离散化,得到离散的状态方程和观测方程;根据测得的数据和离散的状态方程和观测方程,使用基于Sage‑Husa滤波的平方根容积卡尔曼滤波算法对锂电池SOC和SOP进行估计。本发明能够减小噪声对SOC估计的影响;提高了数值计算的稳定性,并且考虑温度对电池实际放出容量的影响,对电池实际可用容量进行了修正,提高了估计结果的准确性。
Description
技术领域
本发明涉及锂电池技术领域,尤其涉及一种用于估计锂电池SOC和SOP的方法。
背景技术
由于传统石油资源的匮乏以及环境污染的加剧,电动汽车被引入了市场。随着电动汽车的迅速发展,锂电池作为其主流动力来源则备受关注。电池作为电动汽车的重要组成部分,如何对其进行有效的管理并且最大程度地发挥电池的性能提高整车的续航里程是一项关键技术。电池SOC是电池管理系统的一个重要参数,SOC反映的是电池的可用容量状态,对SOC的准确估计能够更有效地预测电动汽车的续航里程,电池SOP的准确预测对于实现整车功率分配、能量优化以及延长电池组的使用寿命具有重要意义。
目前,现有的电池SOC估计方法主要有安时积分法、神经网络法等,安时积分法由于在充放电运行过程中误差不断积累导致运行结果不准确,是一种开环的估计方法,神经网络法估计的准确性受到参数和训方法的影响,实际运行时需要获取大量的数据,适用性较差。SOP估计的方法主要有查表法、基于复合脉冲的方法,查表法原理简单、清晰,但是不适用于动态估计。基于复合脉冲法使用的模型比较简单,无法准确描述锂电池的动态特性,并且忽略了电池本身设计的电流限制,可能会导致电池出现过充和过放现象从而威胁电池本身安全。
发明内容
本发明目的就是为了弥补已有技术的缺陷,提供一种用于估计锂电池SOC和SOP的方法,采用了二阶等效电路模型,结合最小二乘法实现电池模型参数的辨识。针对卡尔曼滤波在状态估计时出现的噪声不确定性导致的精度降低问题,提出了基于Sage-Husa滤波的平方根容积卡尔曼滤波法来实现锂电池SOC的准确估计,通过电压、电流、SOC三者约束来完成电池SOP的预测。
本发明是通过以下技术方案实现的:
一种基于Sage-Husa滤波的平方根容积卡尔曼滤波锂电池SOC和SOP联合估计方法,对锂电池进行不同温度下的充放电实验获取锂电池的端电压和电流数据以及不同温度下锂电池实际放出的可用容量,将容量比定义为不同温度下锂电池实际放出的可用容量与电池额定容量的比值,拟合出温度与容量比的关系;
建立能反映电池特性的二阶等效电路模型,对等效电路模型方程进行离散化,得到离散化后的状态空间方程和观测方程;
根据不同温度下的充放电实验数据,确立开路电压、温度和SOC的关系,结合实验测得的电压和电流数据使用最小二乘法对电池模型参数进行辨识;
根据建立的状态方程和观测方程,利用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法对电池SOC和SOP进行估计。
具体如下:对电池进行不同温度下的充放电实验获取电池的端电压和电流值,以及不同温度下锂电池实际放出的可用容量值,拟合出温度与容量比的关系式,f为容量比,不同温度下电池可用容量可表示为:Q(T0)=fQ0,Q0为电池的额定容量;所建立的二阶等效电路模型数学方程如下:
式(1)中,UL表示电池模型输出端电压,Uoc表示电池的开路电压,在充放电实验停止时充分静置后测得,通过拟合的方式获得电池SOC和开路电压以及温度三者的关系;
Upol1为电化学极化内阻Rpol1端的电压,I为电流值,Upol2为浓度差极化内阻Rpol2端的电压,Cpol1为电化学极化电容,Cpol2为浓度差极化电容,Ri为电池的欧姆内阻。
对等效电路模型方程式(1)进行离散化得到如下方程:
式(2)中,UL,k表示k时刻的端电压,Ik为k时刻电流,Δt为采样周期,Upol1,k为k时刻Rpol1两端的电压,Upol2,k为k时刻Rpol2两端的电压。
在进行HPPC测试时,电池由静置状态变为开始放电状态时,端电压会产生一个突变值,利用该变化值除以当前的放电电流值可得到欧姆内阻Ri,在电流消失后,两个RC回路失去了外加激励,相当于零输入响应,此时电池端电压方程:
使用指数拟合方法,拟合出Upol1(0)、Upol2(0)、Rpol1Cpol1、Rpol2Cpol2的值,Upol1(0)、Upol2(0)分别为两个极化内阻的初始电压值。在进行放电之前,电池有一段时间处于静置状态,此时两个极化内阻的端电压值为0,两个RC回路会产生零状态响应,此时端电压方程为:
对此式进行拟合可得到IRpol1、IRpol2得值,进一步得到Rpol1、Rpol2的值,结合(3)式求出的Rpol1Cpol1、Rpol2Cpol2的值,可求得Cpol1、Cpol2的值。在不同的温度条件下利用上述方法对Ri、Rpol1、Rpol2、Cpol1、Cpol2进行辨识可以得到上述参数在不同温度下的值。
根据等效电路模型离散化方程,选取两个极化内阻端的电压Upol1、Upol2和SOC作为状态变量,电池端电压UL作为观测变量,建立下式所述的状态方程和观测方程:
观测方程:UL,k=Uoc(SOCk,T0)-Upol1,k-Upol2,k-IkRi+vk (6)
SOCk表示k时刻的SOC值,w1,k表示与其对应的噪声;Upol1,k表示k时刻电化学极化内阻两端的电压,w2,k表示与其对应的噪声;Upol2,k表示k时刻浓度差极化内阻两端的电压,w3,k表示与其对应的噪声。η为电池充放电效率,UL,k表示k时刻电池的端电压,vk是与其对应的噪声。
结合上述状态方程和观测方程,利用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法进行SOC和SOP估计,对应的状态方程和观测方程如下:
xk+1=F(xk,uk)+wk=Akxk+Bkuk+wk (7)
yk=H(xk,uk)+vk=Ckxk+Dkuk+vk (8)
上述式中,xk+1为状态变量表示下一时刻的状态值,xk表示当前时刻的状态值,uk表示当前时刻的输入变量,wk表示当前时刻的过程噪声,yk为观测变量表示当前时刻的观测值,vk表示当前时刻的测量噪声。Ak为状态转移矩阵,Bk为状态方程系数矩阵,Ck为观测矩阵,Dk为观测方程系数矩阵。F(xk,uk)表示系统状态模型函数,H(xk,uk)表示观测模型函数。
运用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法进行SOC和SOP估计步骤:
初始化,设置参数初始值,x0、S0、Q0、R0以及e的值,k=0, x0为状态变量初值,S0为系统误差协方差初始值的平方根,Q0为过程噪声协方差初值,R0为量测噪声协方差初值,e为遗忘因子值,通常取值为0.95~0.98。
时间更新,首先进行容积点计算,ξi为容积点集,为k时刻状态变量估计值,n为状态变量维数,然后进行容积点传播,即xi,k+1|k=F(xi,k,uk),i=1,2,...,2n,进一步得到k+1时刻的状态变量先验估计值:计算k+1时刻的系统状态变量先验估计值的误差协方差:
量测更新,先对容积点进行计算:然后进行容积点传播:Zi,k+1=H(xi,k+1,uk+1),i=1,2,...,2n,Zi,k+1为k+1时刻第i个容积点对应的观测值,观测值预测:观测误差协方差平方根计算:
卡尔曼增益计算、状态变量估计值及误差协方差平方根更新计算,卡尔曼增益计算:Gk+1为k+1时刻卡尔曼增益,状态变量估计值更新:zk+1为k+1时刻测量的端电压值,系统误差协方差平方根更新:过程噪声协方差和量测噪声协方差更新:
当k<μ,μ为观测样本数,令k=k+1,重复步骤(2)、(3)、(4)、(5),否则结束。
依据电压约束、电流约束、SOC约束来进行SOP预测:
(1)电压约束:利用端电压值可以求出当前时刻的峰值电流值:
电流值可用下式表示:
由于电池在实际进行充放电时电压会有一定的限制,将截止电压的上限值表示为UL,max,下限值表示为UL,min,可得到电压约束下的充放电峰值电流。
放电峰值电流:
充电峰值电流:
(2)电流约束:电池本身在设计就存在充电和放电电流的限制,通常是指电池放电的最大电流和充电的最大电流它们分别可表示为Idis和Icha。
(3)SOC约束:为了更好地对动力电池组进行保护,降低出现过充和过放的风险,需要对SOC设置一定的限制条件,其最小值可以表示为SOCmin,最大值可以表示为SOCmax,对应放电峰值电流和充电峰值电流可用下式表示:
基于上述三者约束条件下放电峰值电流和充电峰值电流可表示为:
本发明的优点是:本发明基于Sage-Husa滤波的平方根容积卡尔曼滤波算法对锂电池进行SOC估计,能够对系统噪声和过程噪声进行实时修正,减小噪声对SOC估计的影响;相比于传统的卡尔曼滤波算法,本发明所提算法能够改善传统卡尔曼滤波算法在状态估计时出现的数值计算发散现象从而导致估计精度降低的问题,通过矩阵分解形式直接传播和更新状态协方差矩阵的平方根,保证了协方差矩阵的非负定性,提高了数值计算的稳定性,并且考虑了温度对电池实际放出容量的影响,对电池实际可用容量进行了修正,提高了估计结果的准确性。
附图说明
图1本发明方法实现示意图。
图2锂电池等效电路模型图。
图3算法实现流程图。
图4温度与电池可用容量关系曲线图。
图5SOC估计结果图。
图6SOP估计结果图。
具体实施方式
如图1、3所示,一种用于估计锂电池SOC和SOP的方法,首先对锂电池进行不同温度下的充放电实验获取锂电池的端电压和电流数据以及不同温度下锂电池实际放出的可用容量,拟合出温度与容量比的关系式;具体如下:
对电池进行不同温度下的充放电实验获取电池的端电压和电流值,以及不同温度下锂电池实际放出的可用容量值,拟合出温度与容量比的关系式,如图4所示,不同温度下电池可用容量可表示为:Q(T0)=fQ0,Q0为电池的额定容量。
如图2所示,对电池建立二阶等效电路模型;然后通过最小二乘法实现对模型参数的辨识,对建立的等效电路模型的方程进行离散化,建立离散的状态方程和观测方程;具体如下:
所建立的二阶等效电路模型数学方程如下:
式(1)中,UL表示电池模型输出端电压,Uoc表示电池的开路电压,在充放电实验停止时充分静置后测得,通过拟合的方式获得电池SOC和开路电压以及温度三者的关系;
Upol1为电化学极化内阻Rpol1端的电压,I为电流值,Upol2为浓度差极化内阻Rpol2端的电压,Cpol1为电化学极化电容,Cpol2为浓度差极化电容,Ri为电池的欧姆内阻。
对等效电路模型方程式(1)进行离散化得到如下方程:
式(2)中,UL,k表示k时刻的端电压,Ik为k时刻电流,Δt为采样周期,Upol1,k为k时刻Rpol1两端的电压,Upol2,k为k时刻Rpol2两端的电压。
在进行HPPC测试时,电池由静置状态变为开始放电状态时,端电压会产生一个突变值,利用该变化值除以当前的放电电流值可得到欧姆内阻Ri,在电流消失后,两个RC回路失去了外加激励,相当于零输入响应,此时电池端电压方程:
使用指数拟合方法,拟合出Upol1(0)、Upol2(0)、Rpol1Cpol1、Rpol2Cpol2的值,Upol1(0)、Upol2(0)分别为两个极化内阻的初始电压值。在电池进行放电之前,电池有一段时间处于静置状态,此时两个极化内阻的端电压值为0,两个RC回路会产生零状态响应,此时端电压方程为:
对此式进行拟合可得到IRpol1、IRpol2得值,进一步得到Rpol1、Rpol2的值,结合(3)式求出的Rpol1Cpol1、Rpol2Cpol2的值,可求得Cpol1、Cpol2的值。在不同的温度条件下利用上述方法对Ri、Rpol1、Rpol2、Cpol1、Cpol2进行辨识可以得到上述参数在不同温度下的值。
根据等效电路模型离散化方程,选取两个极化内阻端的电压Upol1、Upol2和SOC作为状态变量,电池端电压UL作为观测变量,建立下式所述的状态方程和观测方程:
观测方程:UL,k=Uoc(SOCk,T0)-Upol1,k-Upol2,k-IkRi+vk (6)
SOCk表示k时刻的SOC值,w1,k表示与其对应的噪声;Upol1,k表示k时刻电化学极化内阻两端的电压,w2,k表示与其对应的噪声;Upol2,k表示k时刻浓度差极化内阻两端的电压,w3,k表示与其对应的噪声。η为电池充放电效率,UL,k表示k时刻电池的端电压,vk是与其对应的噪声。
最后根据测得的数据和离散的状态方程和观测方程,使用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法对锂电池SOC值进行估计,通过电压、电流以及SOC三者约束来完成SOP的预测;预测结果如图5、6所示,具体如下:对应的状态方程和观测方程如下:
xk+1=F(xk,uk)+wk=Akxk+Bkuk+wk (7)
yk=H(xk,uk)+vk=Ckxk+Dkuk+vk (8)
上述式中,xk+1为状态变量表示下一时刻的状态值,xk表示当前时刻的状态值,uk表示当前时刻的输入变量,wk表示当前时刻的过程噪声,yk为观测变量表示当前时刻的观测值,vk表示当前时刻的测量噪声。Ak为状态转移矩阵,Bk为状态方程系数矩阵,Ck为观测矩阵,Dk为观测方程系数矩阵。F(xk,uk)表示系统状态模型函数,H(xk,uk)表示观测模型函数。
运用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法进行SOC和SOP估计步骤:
初始化,设置参数初始值,x0、S0、Q0、R0以及e的值,k=0, x0为状态变量初值,S0为系统误差协方差初始值的平方根,Q0为过程噪声协方差初值,R0为量测噪声协方差初值,e为遗忘因子值,通常取值为0.95~0.98。
时间更新,首先进行容积点计算,ξi为容积点集,为k时刻状态变量估计值,n为状态变量维数,然后进行容积点传播,即xi,k+1|k=F(xi,k,uk),i=1,2,...,2n,进一步得到k+1时刻的状态变量先验估计值:计算k+1时刻的系统状态变量先验估计值的误差协方差:
量测更新,先对容积点进行计算:然后进行容积点传播:Zi,k+1=H(xi,k+1,uk+1),i=1,2,...,2n,Zi,k+1为k+1时刻第i个容积点对应的观测值,观测值预测:观测误差协方差平方根计算:
卡尔曼增益计算、状态变量估计值及误差协方差平方根更新计算,卡尔曼增益计算:Gk+1为k+1时刻卡尔曼增益,状态变量估计值更新:zk+1为k+1时刻测量的端电压值,系统误差协方差平方根更新:过程噪声协方差和量测噪声协方差更新:
当k<μ,μ为观测样本数,令k=k+1,重复步骤(2)、(3)、(4)、(5),否则结束。
依据电压约束、电流约束、SOC约束来进行SOP预测:
(4)电压约束:利用端电压值可以求出当前时刻的峰值电流值:
电流值可用下式表示:
由于电池在实际进行充放电时电压会有一定的限制,将截止电压的上限值表示为UL,max,下限值表示为UL,min,可得到电压约束下的充放电峰值电流。
放电峰值电流:
充电峰值电流:
(5)电流约束:电池本身在设计就存在充电和放电电流的限制,通常是指电池放电的最大电流和充电的最大电流它们分别可表示为Idis和Icha。
(6)SOC约束:为了更好地对动力电池组进行保护,降低出现过充和过放的风险,需要对SOC设置一定的限制条件,其最小值可以表示为SOCmin,最大值可以表示为SOCmax,对应放电峰值电流和充电峰值电流可用下式表示:
基于上述三者约束条件下放电峰值电流和充电峰值电流可表示为:
Claims (1)
1.一种用于估计锂电池SOC和SOP的方法,其特征在于:具体步骤如下:
对锂电池进行不同温度下的充放电实验获取锂电池的端电压和电流数据以及不同温度下锂电池实际放出的可用容量,将容量比定义为不同温度下锂电池实际放出的可用容量与电池额定容量的比值,拟合出温度与容量比的关系;
建立能反映电池特性的二阶等效电路模型,对等效电路模型方程进行离散化,得到离散的状态空间方程和观测方程;
根据不同温度下的充放电实验数据,确立开路电压、温度和SOC的关系,结合实验测得的电压和电流数据使用最小二乘法对电池模型参数进行辨识;
根据建立的状态方程和观测方程,利用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法对电池SOC和SOP进行估计;
所述的对锂电池进行不同温度下的充放电实验获取锂电池的端电压和电流数据以及不同温度下锂电池实际放出的可用容量,拟合出温度与可用容量的关系,具体如下:
对电池进行不同温度下的充放电实验获取电池的端电压和电流值,以及不同温度下锂电池实际放出的可用容量值,拟合出温度与容量比的关系式,f为容量比,a0、a1、a2为拟合系数值,T0表示温度值,不同温度下电池可用容量表示为:Q(T0)=fQ0,Q0为电池的额定容量;
所述的建立能反映电池特性的二阶等效电路模型,对等效电路模型方程进行离散化,得到离散的状态空间方程和观测方程;根据不同温度下的充放电实验数据,确立开路电压、温度和SOC的关系,结合实验测得的电压和电流数据使用最小二乘法对电池模型参数进行辨识,具体如下:
所建立的二阶等效电路模型数学方程如下:
式(1)中,UL表示电池模型输出端电压,Uoc表示电池的开路电压,在充放电实验停止时充分静置后测得,通过拟合的方式获得电池SOC和开路电压以及温度三者的关系;
Upol1为电化学极化内阻Rpol1端的电压,I为电流值,Upol2为浓度差极化内阻Rpol2端的电压,Cpol1为电化学极化电容,Cpol2为浓度差极化电容,Ri为电池的欧姆内阻;
对等效电路模型方程式(1)进行离散化得到如下方程:
式(2)中,UL,k表示k时刻的端电压,Ik为k时刻的电流,Δt为采样周期,Upol1,k为k时刻Rpol1两端的电压,Upol2,k为k时刻Rpol2两端的电压;
在进行HPPC测试时,电池由静置状态变为开始放电状态时,端电压产生一个突变值,利用该突变值除以当前的放电电流值得到欧姆内阻Ri,在电流消失后,两个RC回路失去了外加激励,相当于零输入响应,此时电池端电压方程:
使用指数拟合方法,拟合出Upol1(0)、Upol2(0)、Rpol1Cpol1、Rpol2Cpol2的值,Upol1(0)、Upol2(0)分别为两个极化内阻的初始电压值;在进行放电之前,电池有一段时间处于静置状态,此时两个极化内阻的端电压值为0,两个RC回路会产生零状态响应,此时端电压方程为:
对此式进行拟合得到IRpol1、IRpol2的值,进一步得到Rpol1、Rpol2的值,结合(3)式求出的Rpol1Cpol1、Rpol2Cpol2的值,再求得Cpol1、Cpol2的值;在不同的温度条件下利用上述方法对Ri、Rpol1、Rpol2、Cpol1、Cpol2进行辨识得到上述参数在不同温度下的值;
根据等效电路模型离散化后的方程,选取两个极化内阻端的电压Upol1、Upol2和SOC作为状态变量,UL作为观测变量,建立下式所述的状态方程和观测方程:
观测方程:UL,k=Uoc(SOCk,T0)-Upol1,k-Upol2,k-IkRi+vk (6)
SOCk表示k时刻的SOC值,w1,k表示与其对应的噪声;Upol1,k表示k时刻电化学极化内阻两端的电压,w2,k表示与其对应的噪声;Upol2,k表示k时刻浓度差极化内阻两端的电压,w3,k表示与其对应的噪声;η为电池充放电效率,UL,k表示k时刻模型输出的端电压,vk是与其对应的噪声;
根据建立的状态方程和观测方程,利用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法对电池SOC和SOP进行估计,具体如下:
系统对应的状态方程和观测方程如下:
xk+1=F(xk,uk)+wk=Akxk+Bkuk+wk (7)
yk=H(xk,uk)+vk=Ckxk+Dkuk+vk (8)
上述式中,xk+1为状态变量表示下一时刻的状态值,xk表示当前时刻的状态值,uk表示当前时刻的输入变量,wk表示当前时刻的过程噪声,yk为观测变量表示当前时刻的观测值,vk表示当前时刻的测量噪声;Ak为状态转移矩阵,Bk为状态方程系数矩阵,Ck为观测矩阵,Dk为观测方程系数矩阵;F(xk,uk)表示系统状态模型函数,H(xk,uk)表示观测模型函数;
运用基于Sage-Husa滤波的平方根容积卡尔曼滤波算法估计SOC和SOP,具体步骤如下:
(2)、时间更新,首先进行容积点计算,ξi为容积点集,为k时刻状态变量估计值,n为状态变量维数,然后进行容积点传播,即进一步得到k+1时刻的状态变量先验估计值:计算k+1时刻的系统状态变量先验估计值的误差协方差:式中,Tria()表示对()内的矩阵进行QR分解,为k时刻过程噪声协方差的平方根,
(3)、量测更新,先对容积点进行计算:然后进行容积点传播:Zi,k+1=H(xi,k+1,uk+1),i=1,2,...,2n,Zi,k+1为k+1时刻第i个容积点对应的观测值,观测值预测:观测误差协方差平方根计算:SRk=chol(Rk)为k时刻量测噪声协方差的平方根,互协方差矩阵平方根计算:
(5)、过程噪声协方差和量测噪声协方差更新:
(6)、当k<μ,μ为观测样本数,令k=k+1,重复步骤(2)、(3)、(4)、(5),否则结束;
(7)、依据电压约束、电流约束、SOC约束来进行SOP预测
电压约束:利用端电压值求出当前时刻的峰值电流值:
电流值用下式表示:
放电峰值电流:
充电峰值电流:
电流约束:电池本身就存在充电和放电电流的限制,指电池放电的最大电流和充电的最大电流,分别表示为Idis和Icha;
基于上述三者约束条件下放电峰值电流和充电峰值电流表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010478880.9A CN111537894B (zh) | 2020-05-29 | 2020-05-29 | 一种用于估计锂电池soc和sop的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010478880.9A CN111537894B (zh) | 2020-05-29 | 2020-05-29 | 一种用于估计锂电池soc和sop的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111537894A CN111537894A (zh) | 2020-08-14 |
CN111537894B true CN111537894B (zh) | 2022-11-25 |
Family
ID=71979549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010478880.9A Active CN111537894B (zh) | 2020-05-29 | 2020-05-29 | 一种用于估计锂电池soc和sop的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111537894B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111551869A (zh) * | 2020-05-15 | 2020-08-18 | 江苏科尚智能科技有限公司 | 锂电池低频参数测量方法、装置、计算机设备和存储介质 |
CN112068000B (zh) * | 2020-09-27 | 2022-03-11 | 哈尔滨工业大学(威海) | 一种考虑动力电池耐久性影响的峰值功率预测方法 |
CN112444749B (zh) * | 2020-11-06 | 2021-11-05 | 南京航空航天大学 | 一种基于温度修正模型的锂电池荷电状态联合估计方法 |
CN112379275A (zh) * | 2020-11-23 | 2021-02-19 | 中国电子科技集团公司第十八研究所 | 一种多参数修正的动力电池soc估计方法及估计系统 |
CN112816876B (zh) * | 2020-12-28 | 2021-12-07 | 湖南航天捷诚电子装备有限责任公司 | 一种用于可充电电池的低温电池剩余电量估算方法及装置 |
CN112800708B (zh) * | 2021-01-19 | 2022-08-09 | 武汉理工大学 | 一种基于滚动时域算法的全钒液流电池峰值功率估计方法 |
CN113189489A (zh) * | 2021-05-08 | 2021-07-30 | 东风汽车集团股份有限公司 | 基于电热耦合电池模型的sop预测系统及其预测方法、存储介质 |
CN113238150B (zh) * | 2021-05-24 | 2022-10-04 | 哈尔滨工业大学 | 基于状态估计算法的电池实时发热功率获取方法 |
CN114019397A (zh) * | 2021-11-04 | 2022-02-08 | 山东交通学院 | 一种梯次利用电池峰值功率的预测方法及系统 |
CN114217234B (zh) * | 2021-12-30 | 2023-08-22 | 南通大学 | 一种基于ide-asrckf的锂离子电池参数辨识与soc估计方法 |
CN114924190A (zh) * | 2022-04-29 | 2022-08-19 | 哈尔滨理工大学 | 一种在寒冷环境下对电动汽车锂电池荷电状态(soc)的估算方法 |
CN115144760B (zh) * | 2022-09-01 | 2022-12-13 | 中创新航科技股份有限公司 | 一种电池系统sop的估算方法及装置 |
CN115902667B (zh) * | 2023-02-15 | 2023-05-23 | 广东电网有限责任公司东莞供电局 | 基于权重和容积点自适应的锂电池soc估计方法 |
CN116500444B (zh) * | 2023-04-06 | 2024-03-26 | 重庆大学 | 一种面向运行安全的电动飞行汽车电池多状态联合估计方法 |
CN116930772B (zh) * | 2023-09-15 | 2023-11-24 | 东方电子股份有限公司 | 一种考虑边界约束的电池soc估计方法及装置 |
CN117388715B (zh) * | 2023-12-11 | 2024-02-27 | 西南交通大学 | 一种串联锂电池组的soc和sop联合估计方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301509B (zh) * | 2015-11-12 | 2019-03-29 | 清华大学 | 锂离子电池荷电状态、健康状态与功率状态的联合估计方法 |
CN107192961A (zh) * | 2017-07-12 | 2017-09-22 | 江苏维科新能源科技有限公司 | 新型动力电池soc估算方法 |
CN108508371B (zh) * | 2018-04-09 | 2019-10-22 | 重庆大学 | 一种基于等效电路模型的动力电池soc/soh/sop联合估计方法 |
CN109782181A (zh) * | 2018-12-20 | 2019-05-21 | 宁波飞拓电器有限公司 | 一种基于组合滤波的应急灯电池soc估计方法 |
CN110032812A (zh) * | 2019-04-18 | 2019-07-19 | 河海大学 | 一种基于自适应容积卡尔曼滤波的动态状态估计方法 |
CN110058160B (zh) * | 2019-04-29 | 2021-08-06 | 西安理工大学 | 基于平方根扩展卡尔曼滤波的锂电池健康状态预测方法 |
CN110395141B (zh) * | 2019-06-27 | 2020-10-30 | 武汉理工大学 | 基于自适应卡尔曼滤波法的动力锂电池soc估算方法 |
CN110361642A (zh) * | 2019-07-11 | 2019-10-22 | 中国科学院电工研究所 | 一种电容器荷电状态的预测方法、装置及电子设备 |
CN111060820B (zh) * | 2019-11-28 | 2022-06-17 | 上海理工大学 | 一种基于二阶rc模型的锂电池soc和sop估计方法 |
-
2020
- 2020-05-29 CN CN202010478880.9A patent/CN111537894B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111537894A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111537894B (zh) | 一种用于估计锂电池soc和sop的方法 | |
Xu et al. | State of charge estimation for lithium-ion batteries based on adaptive dual Kalman filter | |
Lai et al. | Co-estimation of state of charge and state of power for lithium-ion batteries based on fractional variable-order model | |
Liu et al. | A new method of modeling and state of charge estimation of the battery | |
CN108196200B (zh) | 一种锂电池健康和荷电状态的联合模拟评估方法 | |
Zhou et al. | A low-complexity state of charge estimation method for series-connected lithium-ion battery pack used in electric vehicles | |
Song et al. | Improved SOC estimation of lithium-ion batteries with novel SOC-OCV curve estimation method using equivalent circuit model | |
Sun et al. | Overview of the types of battery models | |
CN111722118B (zh) | 一种基于soc-ocv优化曲线的锂离子电池soc估计方法 | |
CN110824363B (zh) | 一种基于改进ckf的锂电池soc和soe联合估算方法 | |
CN105425154A (zh) | 一种估计电动汽车的动力电池组的荷电状态的方法 | |
CN110221221A (zh) | 锂离子电池荷电状态和健康状态联合估计方法 | |
Haoran et al. | Lithium battery soc estimation based on extended kalman filtering algorithm | |
CN111142025A (zh) | 一种电池soc估算方法、装置、存储介质及电动汽车 | |
Li et al. | An equivalent circuit model for state of energy estimation of lithium-ion battery | |
Yang et al. | State-of-health estimation for lithium-ion batteries based on decoupled dynamic characteristic of constant-voltage charging current | |
Yang et al. | SOC estimation of lithium-ion battery based on new adaptive fading extended Kalman filter | |
CN112858916A (zh) | 一种基于模型与数据驱动融合的电池包荷电状态估计方法 | |
CN112946481A (zh) | 基于联合h∞滤波的滑模观测器锂离子电池soc估计方法及电池管理系统 | |
Shen et al. | State of charge, state of health and state of function co-estimation of lithium-ion batteries for electric vehicles | |
CN115389936A (zh) | 一种数模混合驱动的锂电池持续峰值功率能力的在线预测方法 | |
Sarrafan et al. | Real-time estimation of model parameters and state-of-charge of lithiumion batteries in electric vehicles using recursive least-square with forgetting factor | |
Lee et al. | A new battery parameter identification considering current, SOC and Peukert's effect for hybrid electric vehicles | |
Sarrafan et al. | Real-time state-of-charge tracking system using mixed estimation algorithm for electric vehicle battery system | |
Zhu et al. | Online state of charge EKF estimation for LiFePO 4 battery management systems |
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 |