CN114910796A - 基于miaukf算法的锂离子电池荷电状态估计方法 - Google Patents

基于miaukf算法的锂离子电池荷电状态估计方法 Download PDF

Info

Publication number
CN114910796A
CN114910796A CN202210607430.4A CN202210607430A CN114910796A CN 114910796 A CN114910796 A CN 114910796A CN 202210607430 A CN202210607430 A CN 202210607430A CN 114910796 A CN114910796 A CN 114910796A
Authority
CN
China
Prior art keywords
state
lithium ion
algorithm
ion battery
miaukf
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.)
Pending
Application number
CN202210607430.4A
Other languages
English (en)
Inventor
姚冲城
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Dianji University
Original Assignee
Shanghai Dianji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Dianji University filed Critical Shanghai Dianji University
Priority to CN202210607430.4A priority Critical patent/CN114910796A/zh
Publication of CN114910796A publication Critical patent/CN114910796A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/382Arrangements for monitoring battery or accumulator variables, e.g. SoC
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/385Arrangements for measuring battery or accumulator variables
    • G01R31/387Determining ampere-hour charge capacity or SoC
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/396Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tests Of Electric Status Of Batteries (AREA)
  • Secondary Cells (AREA)

Abstract

本发明公开了一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其包括:S1针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV‑SOC特性曲线;S2基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度;S3将多新息算法与AUKF算法结合,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,得到MIAUKF算法模型,以提高锂离子电池荷电状态的估计精度;S4基于MIAUKF算法模型对锂离子电池进行荷电状态估计。本方法的估计精度高,解决了传统卡尔曼滤波类算法中,噪声方差固定,估计精度不高,从而造成有偏估计的问题。

Description

基于MIAUKF算法的锂离子电池荷电状态估计方法
技术领域
本发明涉及电池技术领域,尤其涉及一种基于MIAUKF算法的锂离子电池荷电状态估计方法。
背景技术
目前针对电池的SOC估计算法主要有:开路电压法、安时积分法、神经网络法、卡尔曼滤波法等。其中开路电压和安时积分法属于开环算法,会导致SOC测量误差逐渐累积,降低准确性;而神经网络法在理论上具有高度准确性,但对数据要求较高,需要大量训练集才能保证估计值的精度,不利于实际应用;现今,卡尔曼滤波法(Kalman filter,KF)在线性系统中得到了广泛的使用,且经过学者们的不断改进,这类算法正在逐渐完善。有学者提出了一种扩展卡尔曼滤波算法(Extended Kalman filter,EKF),将其推广到非线性系统中,但EKF算法无法顾及到实际过程中的噪声,会对估计结果带来一定误差;有学者利用无迹卡尔曼滤波算法(Unscented Kalman filter,UKF)对电池SOC进行估计,精度较EKF有所提高,并采用在线参数辨识方法,与离线UKF相比,估算精度更高,但在估算过程中不能消除噪声协方差以及状态协方差造成的影响;有学者利用权重定义自适应无迹卡尔曼滤波算法(AUKF),使SOC的估计精度有所提高,但由于电池在放电过程中,内部特性会发生变化,会直接影响SOC的估计结果的准确性;有学者在二阶RC等效电路模型的基础上,采用脉冲放电实验进行参数辨识,最后用分数阶卡尔曼滤波法对电池SOC进行估计,但该方法估计的模型精度较低,导致SOC估计精度降低。
开路电压法与安时积分法属于开环算法,导致SOC测量误差累积,降低准确性;神经网络法对数据要求较高,需大量训练数据才能保证估计值的精度;扩展卡尔曼滤波无法顾及实际过程中的噪声,会对估计结果带来一定误差;无迹卡尔曼滤波在估算SOC过程中不能消除噪声协方差以及状态协方差的影响;分数阶卡尔曼滤波法的模型精度较低,导致SOC估计精度降低。
[1]戴尉阳,唐忠,史晨豪,等.基于QGA-QBP的动力电池SOC估计[J].电源技术,2021,45(07):893-896.
[2]陈德海,华铭,徐王娟,等.优化安时积分法实时预测电动自行车电池SOC[J].电池,2019,49(1):41-45.
[3]Tong S,Lacap J H,Park J W.Battery state of charge estimation usinga load-classifying neural network[J].Journal of Energy Storage,2016,7:236-243.
[4]郭向伟,康龙云,李文彪,等.基于SOH及离线数据分段矫正的锂电池SOC估计[J].电源技术,2017,41(11):1530-1532.
[5]郑力得,董建园,马强.基于参数辨识与AEKF的锂电池SOC估计[J].电源技术,2020.
[6]Xile D,Caiping Z,Jiuchun J.Evaluation of SOC estimation methodbased on EKF/AEKF under noise interference[J].Energy Procedia,2018,152:520-525.
[7]徐劲力,马国庆.基于UKF的在线锂离子电池SOC估算研究[J].电源技术,2019,43(10):1615-1618.
[8]黄鹏超,鄂加强.基于双自适应卡尔曼滤波的锂电池状态估算[J/OL].储能科学与技术[2021-12-14].https://doi.org/10.19799/j.cnki.2095-4239.
[9]鲁伟,续丹,杨晴霞,等.锂电池分数阶建模与荷电状态研究[J].西安交通大学学报,2017(7):124-129.
发明内容
有鉴于现有技术的上述缺陷,本发明提供一种基于MIAUKF算法的锂离子电池荷电状态估计方法,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,以提高锂离子电池荷电状态的估计精度。
为实现上述目的本发明提供了一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其包括
基于MIAUKF算法的锂离子电池荷电状态估计方法,其包括:
S1针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV-SOC特性曲线;
S2基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度;
S3将多新息算法与AUKF算法结合,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,得到MIAUKF算法模型,以提高锂离子电池荷电状态的估计精度;
S4基于MIAUKF算法模型对锂离子电池进行荷电状态估计。
本发明的进一步改进在于:所述二阶RC等效电路模型包括依次串联的锂离子电池的电压源UOC、极化电阻R1、极化电阻R2、电阻R0、与极化电阻R1并联的极化电容C1、与极化电阻R2并联的极化电容C2
本发明的进一步改进在于:开路电压特性实验包括以下步骤:
(1)设置环境温度为25℃,按CCCV方式将电池充满电,恒流0.5C,截至电流设置为0.05C;
(2)静置一段时间,约2h;
(3)以1C恒流放电3min;
(4)静置一段时间,约2h;
(5)重复步骤3和步骤4,直到端电压达到截至放电电压为止;
(6)静置一段时间,约2h。经过脉冲放电后得到的端电压曲线。
本发明的进一步改进在于:利用最小二乘法进行R0,R1,R2,C1,C2参数辨识具体包括:
(1)辨识参数R0:根据二阶RC等效电路模型,当电池在放电/充电电流被执行或停止时,其终端电压将立即下降或上升,因此利用电池开始放电/充电时端子电压的变化来确定欧姆电阻R0,在参数辨识过程中,取SOC>0.2的R0,并对辨识结果取平均值,其计算式如下式所示:
Figure BDA0003671958080000031
式中,Uk(ti)代表ti时刻的端电压;
(2)辨识参数R1,R2,C1,C2
根据二阶RC等效电路模型,并联系统端电压方程如下所示:
Figure BDA0003671958080000032
式中,t0为初始时间;
参数的具体辨识过程如下:
步骤一:τ1、τ2、U1(t)、U2(t)参数的确定:因电池在脉冲放电前需静置一段时间,故静置过程端电压可视为零状态响应,针对电池的端电压响应方程,利用拟合参数箱可确定静置过程的时间常数τ1,τ2,并将时间常数代入上式得U1(t),U2(t),由于在静置过程中电流为0,故并联系统在tc时刻的端电压方程为:
Figure BDA0003671958080000041
步骤二:R1,R2,C1,C2参数的确定:设a点是上一静置过程的结束点,故有U1(ta)=0,U2(ta)=0,由于在电流通断瞬间,电容两端的电压保持不变,因此电池在tc和ta时刻并联系统的端电压相同:
Figure BDA0003671958080000042
由上式确定参数R1和R2,将其代入τ1=R1C12=R2C2,可得参数C1和C2
本发明的进一步改进在于:步骤S3具体包括:
S31引入多新息矩阵Ep,k
Figure BDA0003671958080000043
式中,P为新息长度;
S32将卡尔曼增益Kk扩展为如下式所示的增益矩阵Kp,k
Kp,k=[K1,k,K2,k,…,Kp,k]∈Rn×p
S33利用上式对状态矩阵进行更新修正;
Figure BDA0003671958080000044
本发明的进一步改进在于:步骤S4具体包括:
S41确定状态值的初始值
Figure BDA0003671958080000045
和后验状态误差协方差初始值P0
Figure BDA0003671958080000051
S42在k-1时刻的采样点xk-1为:
Figure BDA0003671958080000052
式中,λ=α(L+b)-L;L为状态变量的维数,状态向量长度为3;α为采样间隔,b为可调参数;取α=0.01,b=0;
平均权重ωm和协方差权重ωc为:
Figure BDA0003671958080000053
其中,β为高斯随机变量,取β=2;
S43更新先验状态值
Figure BDA0003671958080000054
和系统方差预测值Pxx
Figure BDA0003671958080000055
式中,Qk是系统噪声wk协方差矩阵;
S44更新观测值
Figure BDA0003671958080000056
和观测方差预测值Pyy
Figure BDA0003671958080000057
式中,Rk是观测噪声vk协方差矩阵;
S45更新协方差Pxy,卡尔曼增益Kk,后验状态值
Figure BDA0003671958080000058
和后验状态误差协方差Pk
Figure BDA0003671958080000061
S46基于新息协方差矩阵Fk,实现自适应过程;
Figure BDA0003671958080000062
式中,M为协方差匹配值,ei为新息,ei为电压估计值和测量值的差值;
S47引入多新息矩阵,对状态矩阵进行更新修正后,将过去的多个状态值合并,以估计电池在k时刻SOC的状态值;
Figure BDA0003671958080000063
由于在状态测量更新中引入了一系列旧测量量作为校正项,导致输出电压和电流中存在一定的测量误差;另外,若旧测量值与新测量值的影响相同时,旧数据均会造成累积干扰现象,因此在不同新息中引入不同权重因子σ以降低旧测量值的校正效果;
S48引入不同权重因子对状态矩阵进一步修正;
Figure BDA0003671958080000064
为确保当前时间增益的最大权重,权重因子应满足;
σ1≥(σ23+…+σp)
将不同的权重因子做如下定义:
Figure BDA0003671958080000065
式中,a为可调系数,为保证当前增益权重最大,σ1=1。
本发明提供的装置及方法具有以下技术效果:
本发明引入一种新的锂离子电池荷电状态估计方法。本发明首先,针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV-SOC特性曲线;其次,基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度;然后,引入多新息思想改进AUKF算法,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,得到改进后的MIAUKF算法;最后,基于MIAUKF算法对锂离子电池进行荷电状态估计以提高锂离子电池荷电状态的估计精度,解决AUKF算法造成的后验校正信息丢失,以及在传统卡尔曼滤波类算法中,噪声方差固定,估计精度不高,从而造成有偏估计的问题。
附图说明
图1是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的锂离子电池二阶RC等效电路模型图;
图2是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的OCV-SOC拟合关系曲线图;
图3是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的脉冲放电电压曲线图;
图4是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的二阶RC等效电路端电压对比曲线图;
图5是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的MIAUKF算法流程图;
图6是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的在HWFET工况下SOC仿真结果图;
图7是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的在HWFET工况下端电压仿真结果图;
图8是本发明一种基于MIAUKF算法的锂离子电池荷电状态估计方法的在HWFET工况下各算法仿真结果对比曲线图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
为了阐释的目的而描述了本发明的一些示例性实施例,需要理解的是,本发明可通过附图中没有具体示出的其他方式来实现。
如图1至8所示,本发明的实施例包括基于MIAUKF算法的锂离子电池荷电状态估计方法,其包括:
S1针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV-SOC特性曲线;特性曲线如图2所示;
S2基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度;
S3将多新息算法与AUKF算法结合,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,得到MIAUKF算法模型,以提高锂离子电池荷电状态的估计精度;
S4基于MIAUKF算法模型对锂离子电池进行荷电状态估计。
如图1至4所示,本实施例中,为准确估计电动汽车电池组的SOC,需建立电池精确的等效电路模型。由于二阶RC等效电路模型相较于一阶RC模型精度较高,而三阶及以上RC等效电路模型在精度提高方面不明显,并且二阶RC模型的性能优于PNGV模型。因此综合考虑模型的精度和性能,本发明选择二阶RC等效电路模型建立电池的状态空间模型。其中,UOC为锂离子电池的电压源(open circuit voltage,OCV),且OCV随电池SOC的变化而变化,R1和R2为极化电阻,C1和C2为极化电容,Uk为端电压,R0为欧姆电阻,Ik为端电流。
利用二阶RC等效电路,结合基尔霍夫定律,得到电池状态空间模型。由于UOC、C1和C2均是储能元件,故该模型中含有三个状态变量:SOC,UR1C1,UR2C2
Figure BDA0003671958080000081
Figure BDA0003671958080000082
式中,锂电池的库伦效率η通常为1,时间常数τ1=R1C1,τ2=R2C2,URxCx为极化电阻Rx与相应极化电容Cx对应的电压(x取1、2)。
所述二阶RC等效电路模型包括依次串联的锂离子电池的电压源UOC、极化电阻R1、极化电阻R2、电阻R0、与极化电阻R1并联的极化电容C1、与极化电阻R2并联的极化电容C2
本实施例中,开路电压特性实验包括以下步骤:
(1)设置环境温度为25℃,按CCCV方式将电池充满电,恒流0.5C,截至电流设置为0.05C;
(2)静置一段时间,约2h;
(3)以1C恒流放电3min;
(4)静置一段时间,约2h;
(5)重复步骤3和步骤4,直到端电压达到截至放电电压为止;
(6)静置一段时间,约2h。经过脉冲放电后得到的端电压曲线(截取前50万组数据);脉冲放电的端电压波形图如图3所示。
由于电池在放电过程中,电压存在迟滞,故需要在电池放电后静置1~2h,使电压达到稳定值。
为使拟合曲线更加精确,本发明将电池在经放电并且静置2h后测得的OCV和SOC值进行八项式拟合,拟合后的曲线为标准曲线。
图2所示的OCV-SOC曲线是描述电池特性的重要参数,它将为后续的参数辨识与荷电状态估计算法提供准确的开路电压值。
本实施例中,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识具体包括:
(1)辨识参数R0:根据二阶RC等效电路模型,当电池在放电/充电电流被执行或停止时,其终端电压将立即下降或上升,因此利用电池开始放电/充电时端子电压的变化来确定欧姆电阻R0,在参数辨识过程中,取SOC>0.2的R0,并对辨识结果取平均值,其计算式如下式所示:
Figure BDA0003671958080000091
式中,Uk(ti)代表ti时刻的端电压;
(2)辨识参数R1,R2,C1,C2
根据二阶RC等效电路模型,并联系统端电压方程如下所示:
Figure BDA0003671958080000101
式中,t0为初始时间;
参数的具体辨识过程如下:
步骤一:τ1、τ2、U1(t)、U2(t)参数的确定:因电池在脉冲放电前需静置一段时间,故静置过程端电压可视为零状态响应,针对电池的端电压响应方程,利用拟合参数箱可确定静置过程的时间常数τ1,τ2,并将时间常数代入上式得U1(t),U2(t),由于在静置过程中电流为0,故并联系统在tc时刻的端电压方程为:
Figure BDA0003671958080000102
步骤二:R1,R2,C1,C2参数的确定:设a点是上一静置过程的结束点,故有U1(ta)=0,U2(ta)=0,由于在电流通断瞬间,电容两端的电压保持不变,因此电池在tc和ta时刻并联系统的端电压相同:
Figure BDA0003671958080000103
由上式确定参数R1和R2,将其代入τ1=R1C12=R2C2,可得参数C1和C2
本实施例中,步骤S3具体包括:
S31引入多新息矩阵Ep,k
Figure BDA0003671958080000104
式中,P为新息长度;
S32将卡尔曼增益Kk扩展为如下式所示的增益矩阵Kp,k
Kp,k=[K1,k,K2,k,…,Kp,k]∈Rn×p
S33利用上式对状态矩阵进行更新修正;
Figure BDA0003671958080000105
如图5所示,本实施例中,步骤S4具体包括:
S41确定状态值的初始值
Figure BDA0003671958080000111
和后验状态误差协方差初始值P0
Figure BDA0003671958080000112
S42在k-1时刻的采样点xk-1为:
Figure BDA0003671958080000113
式中,λ=α(L+b)-L;L为状态变量的维数,状态向量长度为3;α为采样间隔,b为可调参数;取α=0.01,b=0;
平均权重ωm和协方差权重ωc为:
Figure BDA0003671958080000114
其中,β为高斯随机变量,取β=2;
S43更新先验状态值
Figure BDA0003671958080000115
和系统方差预测值Pxx
Figure BDA0003671958080000116
式中,Qk是系统噪声wk协方差矩阵;
S44更新观测值
Figure BDA0003671958080000117
和观测方差预测值Pyy
Figure BDA0003671958080000118
式中,Rk是观测噪声vk协方差矩阵;
S45更新协方差Pxy,卡尔曼增益Kk,后验状态值
Figure BDA0003671958080000119
和后验状态误差协方差Pk
Figure BDA0003671958080000121
S46基于新息协方差矩阵Fk,实现自适应过程;
Figure BDA0003671958080000122
式中,M为协方差匹配值,ei为新息,ei为电压估计值和测量值的差值;
S47引入多新息矩阵,对状态矩阵进行更新修正后,将过去的多个状态值合并,以估计电池在k时刻SOC的状态值;
Figure BDA0003671958080000123
由于在状态测量更新中引入了一系列旧测量量作为校正项,导致输出电压和电流中存在一定的测量误差;另外,若旧测量值与新测量值的影响相同时,旧数据均会造成累积干扰现象,因此在不同新息中引入不同权重因子σ以降低旧测量值的校正效果;
S48引入不同权重因子对状态矩阵进一步修正;
Figure BDA0003671958080000124
为确保当前时间增益的最大权重,权重因子应满足;
σ1≥(σ23+…+σp)
将不同的权重因子做如下定义:
Figure BDA0003671958080000125
式中,a为可调系数,为保证当前增益权重最大,σ1=1。
本发明提出一种基于MIAUKF算法的锂离子电池荷电状态估计方法。针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV-SOC特性曲线。基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度。引入多新息思想改进AUKF算法,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,以提高锂离子电池荷电状态的估计精度。基于MIAUKF算法对锂离子电池进行荷电状态估计。如图6、图7、图8所示,本申请的算法具有较高的估计精度。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。

Claims (6)

1.一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其包括:
S1.针对锂离子电池的二阶RC等效电路模型进行锂离子电池开路电压特性实验,利用八项式拟合得到OCV-SOC特性曲线;
S2.基于所建二阶RC等效电路模型,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识并验证该模型的精度;
S3.将多新息算法与AUKF算法结合,通过利用多新息误差向量的新校正阶段和卡尔曼增益矩阵更新估计值,得到MIAUKF算法模型,以提高锂离子电池荷电状态的估计精度;
S4.基于MIAUKF算法模型对锂离子电池进行荷电状态估计。
2.根据权利1所述的一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其特征在于:所述二阶RC等效电路模型包括依次串联的锂离子电池的电压源UOC、极化电阻R1、极化电阻R2、电阻R0、与极化电阻R1并联的极化电容C1、与极化电阻R2并联的极化电容C2
3.根据权利1所述的一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其特征在于,开路电压特性实验包括以下步骤:
(1)设置环境温度为25℃,按CCCV方式将电池充满电,恒流0.5C,截至电流设置为0.05C;
(2)静置一段时间,约2h;
(3)以1C恒流放电3min;
(4)静置一段时间,约2h;
(5)重复步骤3和步骤4,直到端电压达到截至放电电压为止;
(6)静置一段时间,约2h。经过脉冲放电后得到的端电压曲线。
4.根据权利1所述的一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其特征在于,利用最小二乘法进行R0,R1,R2,C1,C2参数辨识具体包括:
(1)辨识参数R0:根据二阶RC等效电路模型,当电池在放电/充电电流被执行或停止时,其终端电压将立即下降或上升,因此利用电池开始放电/充电时端子电压的变化来确定欧姆电阻R0,在参数辨识过程中,取SOC>0.2的R0,并对辨识结果取平均值,其计算式如下式所示:
Figure FDA0003671958070000021
式中,Uk(ti)代表ti时刻的端电压;
(2)辨识参数R1,R2,C1,C2
根据二阶RC等效电路模型,并联系统端电压方程如下所示:
Figure FDA0003671958070000022
式中,t0为初始时间;
参数的具体辨识过程如下:
步骤一:τ1、τ2、U1(t)、U2(t)参数的确定:因电池在脉冲放电前需静置一段时间,故静置过程端电压可视为零状态响应,针对电池的端电压响应方程,利用拟合参数箱可确定静置过程的时间常数τ1,τ2,并将时间常数代入上式得U1(t),U2(t),由于在静置过程中电流为0,故并联系统在tc时刻的端电压方程为:
Figure FDA0003671958070000023
步骤二:R1,R2,C1,C2参数的确定:设a点是上一静置过程的结束点,故有U1(ta)=0,U2(ta)=0,由于在电流通断瞬间,电容两端的电压保持不变,因此电池在tc和ta时刻并联系统的端电压相同:
Figure FDA0003671958070000024
由上式确定参数R1和R2,将其代入τ1=R1C12=R2C2,可得参数C1和C2
5.根据权利4所述的一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其特征在于,步骤S3具体包括:
S31引入多新息矩阵Ep,k
Figure FDA0003671958070000031
式中,P为新息长度;
S32将卡尔曼增益Kk扩展为如下式所示的增益矩阵Kp,k
Kp,k=[K1,k,K2,k,…,Kp,k]∈Rn×p
S33利用上式对状态矩阵进行更新修正;
Figure FDA0003671958070000032
6.根据权利5所述的一种基于MIAUKF算法的锂离子电池荷电状态估计方法,其特征在于,步骤S4具体包括:
S41确定状态值的初始值
Figure FDA0003671958070000033
和后验状态误差协方差初始值P0
Figure FDA0003671958070000034
S42在k-1时刻的采样点xk-1为:
Figure FDA0003671958070000035
式中,λ=α(L+b)-L;L为状态变量的维数,状态向量长度为3;α为采样间隔,b为可调参数;取α=0.01,b=0;
平均权重ωm和协方差权重ωc为:
Figure FDA0003671958070000036
其中,β为高斯随机变量,取β=2;
S43更新先验状态值
Figure FDA0003671958070000041
和系统方差预测值Pxx
Figure FDA0003671958070000042
式中,Qk是系统噪声wk协方差矩阵;
S44更新观测值
Figure FDA0003671958070000043
和观测方差预测值Pyy
Figure FDA0003671958070000044
式中,Rk是观测噪声vk协方差矩阵;
S45更新协方差Pxy,卡尔曼增益Kk,后验状态值
Figure FDA0003671958070000045
和后验状态误差协方差Pk
Figure FDA0003671958070000046
S46基于新息协方差矩阵Fk,实现自适应过程;
Figure FDA0003671958070000047
式中,M为协方差匹配值,ei为新息,ei为电压估计值和测量值的差值;
S47引入多新息矩阵,对状态矩阵进行更新修正后,将过去的多个状态值合并,以估计电池在k时刻SOC的状态值;
Figure FDA0003671958070000048
由于在状态测量更新中引入了一系列旧测量量作为校正项,导致输出电压和电流中存在一定的测量误差;另外,若旧测量值与新测量值的影响相同时,旧数据均会造成累积干扰现象,因此在不同新息中引入不同权重因子σ以降低旧测量值的校正效果;
S48引入不同权重因子对状态矩阵进一步修正;
Figure FDA0003671958070000051
为确保当前时间增益的最大权重,权重因子应满足;
σ1≥(σ23+…+σp)
将不同的权重因子做如下定义:
Figure FDA0003671958070000052
式中,a为可调系数,为保证当前增益权重最大,σ1=1。
CN202210607430.4A 2022-05-31 2022-05-31 基于miaukf算法的锂离子电池荷电状态估计方法 Pending CN114910796A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210607430.4A CN114910796A (zh) 2022-05-31 2022-05-31 基于miaukf算法的锂离子电池荷电状态估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210607430.4A CN114910796A (zh) 2022-05-31 2022-05-31 基于miaukf算法的锂离子电池荷电状态估计方法

Publications (1)

Publication Number Publication Date
CN114910796A true CN114910796A (zh) 2022-08-16

Family

ID=82770412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210607430.4A Pending CN114910796A (zh) 2022-05-31 2022-05-31 基于miaukf算法的锂离子电池荷电状态估计方法

Country Status (1)

Country Link
CN (1) CN114910796A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116224099A (zh) * 2023-05-06 2023-06-06 力高(山东)新能源技术股份有限公司 一种动态自适应估算电池soc的方法
CN117074962A (zh) * 2023-08-22 2023-11-17 江南大学 锂离子电池状态联合估计方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116224099A (zh) * 2023-05-06 2023-06-06 力高(山东)新能源技术股份有限公司 一种动态自适应估算电池soc的方法
CN117074962A (zh) * 2023-08-22 2023-11-17 江南大学 锂离子电池状态联合估计方法及系统

Similar Documents

Publication Publication Date Title
CN110441694B (zh) 基于多重渐消因子卡尔曼滤波锂电池荷电状态估计方法
CN107422269B (zh) 一种锂电池在线soc测量方法
CN105301509B (zh) 锂离子电池荷电状态、健康状态与功率状态的联合估计方法
CN110398691B (zh) 基于改进自适应双无迹卡尔曼滤波器的锂离子动力电池SoC估计方法
CN110488194B (zh) 一种基于电化学阻抗模型的锂电池soc估算方法及其系统
CN111722118B (zh) 一种基于soc-ocv优化曲线的锂离子电池soc估计方法
CN114910796A (zh) 基于miaukf算法的锂离子电池荷电状态估计方法
CN111060834A (zh) 一种动力电池健康状态估算方法
CN109669132B (zh) 一种基于变分贝叶斯滤波的电池荷电状态估计方法
CN105093128A (zh) 基于扩展卡尔曼滤波的蓄电池荷电状态估算方法
CN111220920B (zh) 基于h∞无迹卡尔曼滤波算法的退役锂离子电池荷电状态计算方法
CN111060820A (zh) 一种基于二阶rc模型的锂电池soc、sop估计方法
CN106154168B (zh) 数据驱动的动力电池荷电状态估计方法
CN111044906B (zh) 一种基于极大似然准则的锂离子电池能量状态估算方法
CN112269133B (zh) 一种基于预充电路模型参数识别的soc估计方法
CN111198326B (zh) 一种具有抗扰动特性的电池单体短路电阻在线检测方法
CN111537903B (zh) 一种基于hckf的电池soc估计方法
CN107817448B (zh) 一种适用于复杂工况的在线实时监测电池电量的方法
CN110687462B (zh) 一种动力电池soc与容量全生命周期联合估计方法
CN112557906A (zh) 一种动力电池全生命周期内soc和容量在线联合估计方法
CN110716146A (zh) 一种动力电池开路电压的估计方法
CN112989690A (zh) 一种混合动力汽车锂电池多时间尺度荷电状态估计方法
CN112098849A (zh) 基于求积分卡尔曼滤波的锂电池剩余电量估计方法
CN110554320A (zh) 锂离子电池的soc估算方法
CN110412472B (zh) 一种基于正态伽马滤波的电池荷电状态估计方法

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