CN106383322A - 一种多时间尺度双ukf自适应估算soc和电池容量c的方法 - Google Patents
一种多时间尺度双ukf自适应估算soc和电池容量c的方法 Download PDFInfo
- Publication number
- CN106383322A CN106383322A CN201610917226.7A CN201610917226A CN106383322A CN 106383322 A CN106383322 A CN 106383322A CN 201610917226 A CN201610917226 A CN 201610917226A CN 106383322 A CN106383322 A CN 106383322A
- Authority
- CN
- China
- Prior art keywords
- ukf
- battery capacity
- micro
- soc
- 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.)
- Pending
Links
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种多时间尺度双UKF自适应估算SOC和电池容量C的方法,先建立电池状态空间方程;利用多时间尺度双UKF自适应方法估算SOC和电池容量C;对用于宏UKF和微UKF的***参数和状态变量以及***噪声进行初始化;更新宏UKF的预测方程;选择合适的sigma点更新微UKF的预测方程,并更新***的过程噪声;利用微UKF的校正方程,估计微Kalman滤波器的状态变量;当双时间尺度趋近时,将微Kalman滤波器的状态估计值传到宏Kalman滤波器中,作为测量值;对微Kalman滤波器的数据进行更新,宏Kalman滤波器的测量校正,更新电池参数。本发明能够提高BMS对动力电池SOC的估算精度及估算稳定性。
Description
技术领域
本发明涉及新能源汽车的动力电池管理***领域,具体是一种多时间尺度双UKF自适应估算SOC和电池容量C的方法。
背景技术
在节能环保的时代背景下,新能源汽车成为汽车行业发展的必然趋势。动力电池为新能源电动汽车提供动力,是电动汽车的核心组成部分之一,动力电池的性能完全决定了新能源电动车的性能。电池管理***(Battery Management System,简称BMS)主要就是为了能够提高电池的利用率,防止电池出现过充电与过放电,延长电池的使用寿命,监控电池的状态。其中,动力电池的荷电状态(State of Charge,简称SOC,用z表示)是BMS中重要参数,其估算方法是BMS中的关键技术。
基于电池模型估计动力电池SOC的方法是将电池的荷电状态SOC作为状态量建立起标准的状态空间模型从而实现动力电池SOC估算。该方法可以直接应用基于电池状态空间模型的Kalman滤波方法。相对于其他的SOC估算方法,Kalman滤波法应用的比较广泛,它是一种时域估计方法,将现代控制理论中的状态空间思想引入到最优滤波理论中,可处理时变***,并且Kalman滤波采用递归计算方法,计算量和数据存储量小,便于计算机在线计算,而且精确度、鲁棒性好。其中,无损Kalman滤波(Unscented Kalman Filter,简称UKF)方法是在动力电池SOC的估算中研究的比较多的一种方法,特别适用于像动力电池这种强非线性***,可以获得精确度较高的SOC值。
发明内容
本发明的目的在于提供一种提高BMS对动力电池SOC的估算精度及估算稳定性的多时间尺度双UKF自适应估算SOC和电池容量(用C表示)的方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种多时间尺度双UKF自适应估算SOC和电池容量C的方法,包括如下步骤:
步骤一,建立电池状态空间方程;电池状态空间方程如式(1)所示:
其中,xk,l为状态向量,uk,l为确定性的输入量,Yk,l为输出向量;非线性函数F(xk,l,θk,uk,l)和g(xk,uk)关于状态连续可微;wk,l和νk,l分别为***的扰动噪声和观测噪声,两者均为高斯白噪声;Qk为***扰动协方差矩阵,Rk为观测噪声协方差矩阵;
步骤二,利用多时间尺度双UKF自适应方法估算SOC和电池容量C;
(1)对用于宏UKF和微UKF的***参数和状态变量以及***噪声进行初始化;
(2)更新宏UKF的预测方程;
(3)选择合适的sigma点更新微UKF的预测方程,并更新***的过程噪声wk-1,l;
(4)利用微UKF的校正方程,估计微Kalman滤波器的状态变量;
(5)当双时间尺度趋近时,将微Kalman滤波器的状态估计值传到宏Kalman滤波器中,作为测量值;
(6)对微Kalman滤波器的数据进行更新,以进行下一个循环;
(7)宏Kalman滤波器的测量校正,更新电池参数。
作为本发明进一步的方案:在步骤一中,对于一阶RC电池模型,动力电池状态空间方程为:
其中,UOCV为电池开路电压;
通过实验辨识,确定UOCV:
作为本发明进一步的方案:在步骤二中,UKF自适应方法步骤描述如下:
(1)选择n+1个sigma点;
(2)预测状态值;
(3)预测滤波误差方差矩阵;
(4)量测更新方程;
(5)计算Kalman增益;
(6)状态和误差更新方程。
作为本发明进一步的方案:在步骤二中,在步骤二中,对过程噪声和观测噪声的处理使用UKF自适应方法,计算公式如下:
作为本发明进一步的方案:在步骤二中,对电池容量C的估算描述如下:
(1)对电池容量C进行更新,并将电池容量更新分为较大、准确、较小三种;
(2)将电池容量C代入微UKF中进行SOC的估算;
(3)将SOC估计值作为电池容量C估计的输入,从而获得准确的电池容量C的值。
考虑到电池容量C与可测量端电压Ui没有较强的关系,而负载电流i与C有直接关系,但是短期内电流的积累对C的影响很小,这使得准确估算出C比较困难,考虑到两个时间尺度双UKF自适应方法,可以通过较长时间的积累以放大这种效果,并且考虑SOC与电池容量C的关系有
这里运用两个滤波器的时间尺度不同,将i或者说SOC对电池容量C的作用进行放大,并将微Kalman滤波器的输出作为测量值对宏Kalman滤波器的状态变量进行修正,同时UKF的使用也将对微UKF估计的SOC值的精度要求有所降低。
与现有技术相比,本发明的有益效果是:
本发明提出的多时间尺度双UKF的自适应估算方法,通过使用两个UKF滤波器提高BMS对动力电池SOC的估算精度及估算稳定性。其中一个UKF滤波器为微Kalman滤波器,进行电池状态参数(如SOC)的估算;另一个UKF滤波器为宏Kalman滤波器,用来估算电池的参数(如电池容量C)。另外由于所估算量的变化程度快慢不同,两个Kalman滤波器之间以不同时间尺度周期进行数据的交换和更新,同时加入自适应环节来处理噪声的问题。
附图说明
图1为动力电池一阶RC模型图;
图2为多时间尺度双UKF自适应算法框架图;
图3为自适应环节流程图;
图4为电池容量的估算示意图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于电池模型的不确定性、动力电池周围环境的强干扰和强噪声、动力电池的老化导致电池参数具有易变性等因素,SOC和电池容量的估算精度会降低,因此需要对电池模型的参数进行更新。本发明提出的多时间尺度双UKF的自适应估算方法通过使用两个UKF滤波器来解决这些问题。其中一个UKF滤波器为微Kalman滤波器,进行电池状态参数(如SOC)的估算;另一个UKF滤波器为宏Kalman滤波器,用来估算电池的参数(如电池容量C)。另外由于所估算量的变化程度快慢不同,两个Kalman滤波器之间以不同时间尺度周期进行数据的交换和更新,同时加入自适应环节来处理噪声的问题。
实施例1
请参阅图1-4,本发明实施例中,一种多时间尺度双UKF自适应估算SOC与电池容量C的方法进行详细说明。
步骤一,建立电池状态空间方程;
步骤二,利用多时间尺度双UKF自适应方法估算SOC和电池容量C。
在步骤一中,建立的动力电池状态空间方程为:
其中,xk,l为状态向量,uk,l为确定性的输入量,Yk,l为输出向量;非线性函数F(xk,l,θk,uk,l)和g(xk,uk)关于状态连续可微;wk,l和νk,l分别为***的扰动噪声和观测噪声,两者均为高斯白噪声;Qk为***扰动协方差矩阵,Rk为观测噪声协方差矩阵。
对于一阶RC电池模型,可以表述成:
其中,UOCV为电池开路电压。
通过实验辨识,可以确定UOCV:
由于模型参数变化缓慢且采样时间足够小,并将UD解出同时表示为离散形式:
其中CDiffRDiff的值一般可以看作是常数τ,通过辨识τ=35.6,则式(1-4)可以写成:
SOC的理论计算公式为:
η(iL(τ))为不同放电倍率下修正因子;C为电池最大可用容量,是时间温度函数。则式(1-6)的离散形式为
将式(1-6)、(1-7)写成状态空间的形式有:
基于上述公式***的状态空间表达式如下:
其中,参数式如下:
在步骤二中,UKF自适应方法步骤描述如下:
(1)选择n+1个sigma点;
对sigma点使用对称采样策略,则对应的权重值为:
(2)预测状态值;
(3)预测滤波误差方差矩阵;
(4)量测更新方程;
(5)计算Kalman增益;
(6)状态和误差更新方程。
在步骤二中,自适应环节描述如下:
对过程噪声和观测噪声的处理使用自适应的算法,自适应算法如图3所示,计算公式如下:
在步骤二中,对电池容量C的估算如图4所示,描述如下:
(1)对电池容量C进行更新,并将容量更新分为较大、准确、较小三种;
(2)将电池容量C代入微UKF中进行SOC的估算,由于时间尺度的不同,会进行多次SOC的估算;
(3)将SOC估计值作为电池容量C估计的输入,从而获得准确的电池容量C的值。
考虑到电池容量C与可测量端电压Ui没有较强的关系,而负载电流i与电池容量C有直接关系,但是短期内电流的积累对电池容量C的影响很小,这使得准确估算出电池容量C比较困难,考虑到两个时间尺度双UKF自适应方法,可以通过较长时间的积累以放大这种效果,并且考虑SOC与电池容量C的关系有
这里运用两个滤波器的时间尺度不同,将i或者说SOC对电池容量C的作用进行放大,并将微Kalman滤波器的输出作为测量值对宏Kalman滤波器的状态变量进行修正,同时UKF的使用也将对微UKF估计的SOC值的精度要求有所降低。
步骤二,利用多时间尺度双UKF自适应方法估算SOC和电池容量C的步骤描述如下:
(1)对用于宏UKF和微UKF的***参数和状态变量以及***噪声进行初始化;
其中,k∈{1,…,∞}
(2)更新宏UKF的预测方程;
(3)根据公式(2-1)(2-2)选择合适的sigma点更新微UKF的预测方程,并根据公式(2-8)(2-9)(2-10)(2-11)更新***的过程噪声wk-1,l,第一次时使用初始化的噪声;
(4)利用微UKF的校正方程,估计微Kalman滤波器的状态变量。
(5)当双时间尺度趋近时,将微Kalman滤波器的状态估计值传到宏Kalman滤波器中,作为测量值。
(6)对微Kalman滤波器的数据进行更新,以进行下一个循环;
(7)宏Kalman滤波器的测量校正,更新电池参数。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (5)
1.一种多时间尺度双UKF自适应估算SOC和电池容量C的方法,其特征在于,包括如下步骤:
步骤一,建立电池状态空间方程;电池状态空间方程如式(1)所示:
其中,xk,l为状态向量,uk,l为确定性的输入量,Yk,l为输出向量;非线性函数F(xk,l,θk,uk,l)和g(xk,uk)关于状态连续可微;wk,l和νk,l分别为***的扰动噪声和观测噪声,两者均为高斯白噪声;Qk为***扰动协方差矩阵,Rk为观测噪声协方差矩阵;
步骤二,利用多时间尺度双UKF自适应方法估算SOC和电池容量C;
(1)对用于宏UKF和微UKF的***参数和状态变量以及***噪声进行初始化;
(2)更新宏UKF的预测方程;
(3)选择合适的sigma点更新微UKF的预测方程,并更新***的过程噪声wk-1,l;
(4)利用微UKF的校正方程,估计微Kalman滤波器的状态变量;
(5)当双时间尺度趋近时,将微Kalman滤波器的状态估计值传到宏Kalman滤波器中,作为测量值;
(6)对微Kalman滤波器的数据进行更新,以进行下一个循环;
(7)宏Kalman滤波器的测量校正,更新电池参数。
2.根据权利要求1所述的多时间尺度双UKF自适应估算SOC和电池容量C的方法,其特征在于,在步骤一中,对于一阶RC电池模型,动力电池状态空间方程为:
其中,UOCV为电池开路电压;
通过实验辨识,确定UOCV:
3.根据权利要求1所述的多时间尺度双UKF自适应估算SOC和电池容量C的方法,其特征在于,在步骤二中,UKF自适应方法步骤描述如下:
(1)选择n+1个sigma点;
(2)预测状态值;
(3)预测滤波误差方差矩阵;
(4)量测更新方程;
(5)计算Kalman增益;
(6)状态和误差更新方程。
4.根据权利要求1所述的多时间尺度双UKF自适应估算SOC和电池容量C的方法,其特征在于,在步骤二中,对过程噪声和观测噪声的处理使用UKF自适应方法,计算公式如下:
5.根据权利要求1所述的多时间尺度双UKF自适应估算SOC和电池容量C的方法,其特征在于,在步骤二中,对电池容量C的估算,具体步骤如下:
(1)对电池容量C进行更新,并将电池容量更新分为较大、准确、较小三种;
(2)将电池容量C代入微UKF中进行SOC的估算;
(3)将SOC估计值作为电池容量C估计的输入,从而获得准确的电池容量C的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610917226.7A CN106383322A (zh) | 2016-10-21 | 2016-10-21 | 一种多时间尺度双ukf自适应估算soc和电池容量c的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610917226.7A CN106383322A (zh) | 2016-10-21 | 2016-10-21 | 一种多时间尺度双ukf自适应估算soc和电池容量c的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106383322A true CN106383322A (zh) | 2017-02-08 |
Family
ID=57958735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610917226.7A Pending CN106383322A (zh) | 2016-10-21 | 2016-10-21 | 一种多时间尺度双ukf自适应估算soc和电池容量c的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106383322A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153163A (zh) * | 2017-06-20 | 2017-09-12 | 福建工程学院 | 一种基于自适应ukf的锂电池soc估算方法 |
CN109298340A (zh) * | 2018-10-19 | 2019-02-01 | 同济大学 | 一种基于可变时间尺度的电池容量在线估计方法 |
CN109814041A (zh) * | 2019-01-16 | 2019-05-28 | 上海理工大学 | 一种锂离子电池双卡尔曼滤波容量估计方法 |
CN112485672A (zh) * | 2020-11-13 | 2021-03-12 | 上海电气集团股份有限公司 | 一种电池状态确定方法及装置 |
CN112731157A (zh) * | 2020-12-16 | 2021-04-30 | 上海理工大学 | 一种基于数据驱动的锂离子电池容量估计方法 |
CN113030752A (zh) * | 2021-04-12 | 2021-06-25 | 安徽理工大学 | 一种基于变遗忘因子在线参数辨识和soc联合估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054084A (zh) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | 一种动力电池soc估计方法 |
-
2016
- 2016-10-21 CN CN201610917226.7A patent/CN106383322A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054084A (zh) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | 一种动力电池soc估计方法 |
Non-Patent Citations (1)
Title |
---|
熊刚: "车用锂离子电池SOC和电池容量估算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153163A (zh) * | 2017-06-20 | 2017-09-12 | 福建工程学院 | 一种基于自适应ukf的锂电池soc估算方法 |
CN109298340A (zh) * | 2018-10-19 | 2019-02-01 | 同济大学 | 一种基于可变时间尺度的电池容量在线估计方法 |
CN109814041A (zh) * | 2019-01-16 | 2019-05-28 | 上海理工大学 | 一种锂离子电池双卡尔曼滤波容量估计方法 |
CN112485672A (zh) * | 2020-11-13 | 2021-03-12 | 上海电气集团股份有限公司 | 一种电池状态确定方法及装置 |
CN112731157A (zh) * | 2020-12-16 | 2021-04-30 | 上海理工大学 | 一种基于数据驱动的锂离子电池容量估计方法 |
CN113030752A (zh) * | 2021-04-12 | 2021-06-25 | 安徽理工大学 | 一种基于变遗忘因子在线参数辨识和soc联合估计方法 |
CN113030752B (zh) * | 2021-04-12 | 2024-03-29 | 安徽理工大学 | 一种基于变遗忘因子在线参数辨识和soc联合估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106383322A (zh) | 一种多时间尺度双ukf自适应估算soc和电池容量c的方法 | |
Yu et al. | Online capacity estimation for lithium-ion batteries through joint estimation method | |
CN103995464B (zh) | 一种估计电动车辆的动力***的参数和状态的方法 | |
CN103941195B (zh) | 基于小模型误差准则扩展卡尔曼滤波的电池soc估计方法 | |
Farmann et al. | A comprehensive review of on-board State-of-Available-Power prediction techniques for lithium-ion batteries in electric vehicles | |
CN103616647B (zh) | 一种用于电动汽车电池管理***的电池剩余电量估计方法 | |
CN109188293B (zh) | 基于新息协方差带渐消因子的ekf锂离子电池soc估算方法 | |
CN101971043B (zh) | 在电池管理***中测量电池的充电状态的方法及其装置 | |
CN106054080A (zh) | 一种动力电池荷电状态与健康状态的联合估计方法 | |
CN104617623B (zh) | 一种电动汽车动力电池组均衡控制方法 | |
CN109061506A (zh) | 基于神经网络优化ekf的锂离子动力电池soc估计方法 | |
CN108896924B (zh) | 基于深度置信网络和卡尔曼滤波融合的锂离子电池荷电状态估计方法 | |
CN107576919A (zh) | 基于armax模型的动力电池荷电状态估算***及方法 | |
CN105974327A (zh) | 一种基于神经网络和ukf的锂电池组soc预测方法 | |
CN110346734A (zh) | 一种基于机器学习的锂离子动力电池健康状态估算方法 | |
CN106483470A (zh) | 一种基于未来运行工况预测的电池剩余放电能量预测方法 | |
CN107167743A (zh) | 基于电动车辆的荷电状态估算方法和装置 | |
CN104380128A (zh) | 蓄电池***和用于确定蓄电池***的蓄电池单元或蓄电池模块的内阻的所属的方法 | |
CN105929340A (zh) | 一种基于arima估算电池soc的方法 | |
CN109738810A (zh) | 一种电池剩余使用寿命的综合估算方法 | |
WO2010132303A1 (en) | Method for determining extent and type of capacity fade | |
CN109633472A (zh) | 一种锂电池单体的荷电状态估计算法 | |
Li et al. | Real-time estimation of lead-acid battery parameters: A dynamic data-driven approach | |
CN105223512A (zh) | 基于电池特性的动态矫正剩余电量的方法 | |
CN105044606A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170208 |
|
RJ01 | Rejection of invention patent application after publication |