CN110133507B - 一种基于narx-ukf算法的电池剩余电量估计方法 - Google Patents
一种基于narx-ukf算法的电池剩余电量估计方法 Download PDFInfo
- Publication number
- CN110133507B CN110133507B CN201910329199.5A CN201910329199A CN110133507B CN 110133507 B CN110133507 B CN 110133507B CN 201910329199 A CN201910329199 A CN 201910329199A CN 110133507 B CN110133507 B CN 110133507B
- Authority
- CN
- China
- Prior art keywords
- value
- battery
- narx
- network
- ukf
- 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
- G01R31/387—Determining ampere-hour charge capacity or SoC
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种基于NARX‑UKF算法的电池剩余电量估计方法。现有的各种方法存在各种问题。本发明方法首先测量电池在不同条件下的电压和电流,然后对测量数据进行预处理,再搭建NARX‑UKF网络,把处理后的电压电流训练数据输入NARX网络中来训练网络,达到训练目标后,输入测试数据,NARX网络的输入结果即为SOC的估计值同时也是UKF模型的输入量,该估计值经过UKF模型后,得到的状态更新值即为当前时刻所估计得到的电池剩余电量。本发明的不需要建立电池模型,只需要常见的可测量量便可以快速精确的进行电池剩余电量估计。有着模型训练速度快,所需参数少,估计精度高的优点。
Description
技术领域
本发明属于电池技术领域,提出了一种基于NARX(基于非线性自适应回归神经网络)- UKF(无迹卡尔曼滤波器)算法的电池剩余电量估计的解决方案。
背景技术
近些年来,随着经济的发展,汽车保有量的不断增加,环境问题也日益严重。使用高效环保的代替能源是解决污染问题和传统能源衰竭问题的有效途径之一,具有低噪声、低污染、低排放和能量利用率高等优点,因此新能源汽车得到了各国的大力发展。电池管理系统(BMS,Battery Management System)是电动汽电池的关键技术,电池管理系统的发展对电动汽车具有着相当重大的意义,BMS作为电池系统的大脑,控制着电池的充放电过程,避免电池出现过充过放的问题,实现对电池的有效的管理和监控,降低电池使用过程中存在的风险,最终达到延长电池使用寿命和对电池的最有效的利用的目的。
荷电状态估计(SOC,State Of Charge)是电池管理系统需要计算的主要参数之一,SOC 可以直接反应电池剩余电量,定义为电池所剩电量和电池额定容量的比值,充当燃油汽车中燃油表。SOC估计对于电池系统的正常运转是必不可少的,是电动汽车续航行驶和防止电池过度充电过度放电的主要依据,对电动汽车的整体性能还有电池的寿命效率有重要的意义。电池的SOC参数具有十分重要的意义,但是不可以直接测量。电池是一个高度非线性的系统,影响SOC的因素很多,使得SOC估计具有一定的难度。
目前,已经被提出的SOC估计方法主要有:安时积分法、开路电压法、卡尔曼、神经网络法等。安时积分法在实际应用条件下有噪声误差会造成计较大的累积误差,开路电压法需要电池长时间的静置,卡尔曼滤波法对电池建模的精确性要求较高,一般的神经网络法估计 SOC,估计结果的跳动较大。
发明内容
本发明的目的就是克服已有技术的不足,为了满足较高精度的估计电池SOC值,提出了一种基于NARX-UKF算法的电池剩余电量估计方法。
本发明利用具有外部输入的非线性自回归网络(NARXNN,NonlinearAutoregressive with Exogenous Inputs Neural Networks)来估计锂电池的SOC值,并且加入了无迹卡尔曼滤波器(UKF,Unscented Kalman Filter)来减小NARXNN模型的估计误差。具体步骤是:
步骤(1).对特定的单节锂电池进行充放电实验,测量电池工作电流和工作电压,测量间隔为Δt为1s,并记录温度;使用安时积分法计算电池的SOC值,作为目标值:其中SOC(n)是电池在第n个测量点时的SOC值,η为库伦效率,I为电流值的大小,充电时为负值,放电时为正值,QN是电池的额定容量。
库伦电效率η按如下方法确定:
(1-1).将以1/30倍额定电流放电时的额定容量为Qn的完全充满电的电池以不同放电速率Ci恒流放电N次,0<Ci≤C,1≤i≤N,N≥10,放空电池电量所需要的时间和放电电流值的乘积即为相应放电速率下的电池总电量Qi;
步骤(2).处理测量数据:
对电压电流进行标准化处理,处理后的数据x的范围在[-1,1]之间,处理方法如下:
步骤(3).搭建NARXNN模型,输入延时步长设置为5~15,输出延时步长设置为2~5,神经元节点数设置为15;一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为 R的NARX网络的数学模型表示为:
其中,fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值。
步骤(4).将步骤(2)得到的数据输入步骤(3)得到的网络中,对网络进行训练,训练的性能函数为:ym、tm分别是第m个测量点时的网络SOC预测值和SOC测量值。网络的MSE训练目标在0℃、25℃和45℃下分别设置为8e-06、1.1e-05和 1.9e-05。
步骤(5).输入测试集,测试网络性能。
步骤(6).为NARX网络增加UKF模型,建立如下状态方程和测量方程:
UKF的输出即为SOC估计值。UKF修正估计值的方法如下:
其中,x0为初始时刻的状态量。
其中,和Pk-1分别是K维状态变量的均值和方差,表示矩阵的第i列,λ=α2(K+κ)-K为比例缩放参数,α为控制采样点分布密集程度的参数,取值范围为:1e-4≤α≤1,κ的取值需保证矩阵半正定性,通常取值为3-K,β为状态分布参数。
(6-4).测量更新:
卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k)-1;
本发明的训练集和测试集均采用工况下的放电数据,贴合电动汽车实际行驶时的电流电压变化,可以方便地进行电池SOC的快速估计,不需要对电池建模,只需要易于测量的量,该方法收敛速度快,模型训练时间只需要十几秒,估计精度高、结果误差小、计算量小、波动小,所需参数数目少,便于移植到硬件平台中,而且适用于各种电池SOC的快速估计。
附图说明:
图1为本发明系统的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
参见图1,一种基于NARX-UKF算法的电池剩余电量估计方法,具体步骤是:
步骤(1).对特定的单节锂电池进行充放电实验,测量电池工作电流和工作电压,测量间隔为Δt为1s,并记录温度;使用安时积分法计算电池的SOC值,作为目标值:其中SOC(n)是电池在第n个测量点时的SOC值,η为库伦效率,即充放电效率,I为电流值的大小,充电时为负值,放电时为正值,QN是电池的额定容量。
充放电效率η按如下方法确定:
(1-1).将以1/30倍额定电流放电时的额定容量为Qn的完全充满电的电池以不同放电速率Ci恒流放电N次,0<Ci≤C,1≤i≤N,N≥10,放空电池电量所需要的时间和放电电流值的乘积即为相应放电速率下的电池总电量Qi;
步骤(2).处理测量数据:
对电压电流进行标准化处理,处理后的数据x的范围在[-1,1]之间,处理方法如下:
步骤(3).搭建NARXNN模型,输入延时步长设置为5~15,输出延时步长设置为2~5,神经元节点数设置为15;一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为R,在本实施例中,0℃、25℃和45℃下设置输入延时步长分别为8、10和8,输出延时步长分别为3、4和3,隐藏层神经元个数设为15,输出层和隐藏层的激活函数分别为tansig 和purelin函数。搭建的NARX网络的数学模型表示为:
其中,fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值。
步骤(4).将步骤(2)得到的数据输入步骤(3)得到的网络中,对网络进行训练,利用 Levenberg-Marquradt算法训练神经网络,保证目标函数均方误差(MSE,Mean SquaredError)最小,训练的性能函数为:ym、tm分别是第m个测量点时的网络SOC预测值和SOC测量值。为了防止过拟合,网络的MSE训练目标在0℃、25℃和 45℃下分别设置为8e-06、1.1e-05和1.9e-05。
步骤(5).输入测试集,测试网络性能。
步骤(6).为NARX网络增加UKF模型,建立如下状态方程和测量方程:Ek=SOCk+vk;其中,ik是k时刻的电流值,SOCk第k时刻的SOC状态估计值,wk和vk是过程噪声和测量噪声,Ek是NARX网络在k时刻的估计值。
UKF的输出即为SOC估计值。UKF修正估计值的方法如下:
其中,x0为初始时刻的状态量。
其中,和Pk-1分别是K维状态变量的均值和方差,表示矩阵的第i列,λ=α2(K+κ)-K为比例缩放参数,α为控制采样点分布密集程度的参数,取值范围为:1e-4≤α≤1,κ的取值需保证矩阵半正定性,通常取值为3-K。β为状态分布参数,本实施例取值为2。
(6-4).测量更新:
卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k)-1;
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于NARX-UKF算法的电池剩余电量估计方法,其特征在于该方法的具体步骤是:
步骤(1).对单节锂电池进行充放电实验,测量电池工作电流和工作电压,测量间隔为Δt为1s,并记录温度;使用安时积分法计算电池的SOC值,作为目标值:
所述的库伦效率η即为充放电效率,按如下方法确定:
(1-1).将以1/30倍额定电流放电时的额定容量为Qn的完全充满电的电池以不同放电速率Ci恒流放电N次,0<Ci≤C,1≤i≤N,N≥10,放空电池电量所需要的时间和放电电流值的乘积即为相应放电速率下的电池总电量Qi;
步骤(2).处理测量数据:
对电压电流进行标准化处理,处理后的数据x的范围在[-1,1]之间,处理方法如下:
步骤(3).搭建一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为R的NARX网络的数学模型,表示为:其中,fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值;
步骤(5).输入测试集,测试网络性能;
步骤(6).为NARX网络增加UKF模型,建立如下状态方程和测量方程:
Ek=SOCk+vk;其中,ik是k时刻的电流值,SOCk第k时刻的SOC状态估计值,wk和vk是过程噪声和测量噪声,Ek是NARX网络在k时刻的估计值;UKF的输出即为SOC估计值;UKF修正估计值的方法如下:
(6-4).测量更新:
卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k)-1;
2.如权利要求1所述的一种基于NARX-UKF算法的电池剩余电量估计方法,其特征在于,步骤(3)中搭建的NARX网络的数学模型,输入延时步长设置为5~15,输出延时步长设置为2~5,神经元节点数设置为15。
3.如权利要求1所述的一种基于NARX-UKF算法的电池剩余电量估计方法,其特征在于:步骤(4)中MSE训练目标在0℃、25℃和45℃下分别设置为8e-06、1.1e-05和1.9e-05。
4.如权利要求1所述的一种基于NARX-UKF算法的电池剩余电量估计方法,其特征在于:步骤(6)中控制采样点分布密集程度的参数α取值范围为:1e-4≤α≤1。
5.如权利要求1所述的一种基于NARX-UKF算法的电池剩余电量估计方法,其特征在于:步骤(6)中κ的取值需保证矩阵半正定性,通常取值为3-K。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910329199.5A CN110133507B (zh) | 2019-04-23 | 2019-04-23 | 一种基于narx-ukf算法的电池剩余电量估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910329199.5A CN110133507B (zh) | 2019-04-23 | 2019-04-23 | 一种基于narx-ukf算法的电池剩余电量估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133507A CN110133507A (zh) | 2019-08-16 |
CN110133507B true CN110133507B (zh) | 2021-09-14 |
Family
ID=67570930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910329199.5A Active CN110133507B (zh) | 2019-04-23 | 2019-04-23 | 一种基于narx-ukf算法的电池剩余电量估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133507B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110888063A (zh) * | 2019-12-02 | 2020-03-17 | 上海国际港务(集团)股份有限公司 | 一种基于港机轮胎吊磷酸铁锂电池并联系统设计的soc估算法 |
CN112083331A (zh) * | 2020-08-09 | 2020-12-15 | 昆明理工大学 | 一种提升锂离子电池荷电状态估计精度的融合方法 |
CN112557907A (zh) * | 2020-12-17 | 2021-03-26 | 杭州六纪科技有限公司 | 一种基于gru-rnn的电动汽车锂离子电池的soc估计方法 |
CN112782591B (zh) * | 2021-03-22 | 2022-07-22 | 浙江大学 | 基于多电池数据融合的锂电池soh长期预测方法 |
CN113219347A (zh) * | 2021-04-27 | 2021-08-06 | 东软睿驰汽车技术(沈阳)有限公司 | 一种电池参数测量方法及装置 |
CN114050334A (zh) * | 2021-09-27 | 2022-02-15 | 深圳拓邦股份有限公司 | 电池主动放空修正方法及换电设备 |
CN117471327B (zh) * | 2023-12-25 | 2024-05-07 | 深圳市北测检测技术有限公司 | 一种汽车动力电池的安全性能检测方法 |
CN117686920B (zh) * | 2024-02-02 | 2024-04-12 | 河南科技学院 | 一种用于电池系统内单体电池的端电压预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1691209A1 (en) * | 2005-02-14 | 2006-08-16 | Denso Corporation | Method and apparatus for detecting charged state of secondary battery based on neural network calculation |
CN103675706A (zh) * | 2013-12-13 | 2014-03-26 | 桂林电子科技大学 | 一种动力电池电荷量估算方法 |
CN103744028A (zh) * | 2013-12-25 | 2014-04-23 | 广西科技大学 | 一种基于ukf的蓄电池荷电状态soc估算方法 |
CN105629175A (zh) * | 2015-12-29 | 2016-06-01 | 北京航天测控技术有限公司 | 一种基于无迹卡尔曼滤波的锂离子电池寿命预测方法 |
CN109444757A (zh) * | 2018-10-09 | 2019-03-08 | 杭州中恒云能源互联网技术有限公司 | 一种电动汽车动力电池剩余电量估算方法 |
-
2019
- 2019-04-23 CN CN201910329199.5A patent/CN110133507B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1691209A1 (en) * | 2005-02-14 | 2006-08-16 | Denso Corporation | Method and apparatus for detecting charged state of secondary battery based on neural network calculation |
CN103675706A (zh) * | 2013-12-13 | 2014-03-26 | 桂林电子科技大学 | 一种动力电池电荷量估算方法 |
CN103744028A (zh) * | 2013-12-25 | 2014-04-23 | 广西科技大学 | 一种基于ukf的蓄电池荷电状态soc估算方法 |
CN105629175A (zh) * | 2015-12-29 | 2016-06-01 | 北京航天测控技术有限公司 | 一种基于无迹卡尔曼滤波的锂离子电池寿命预测方法 |
CN109444757A (zh) * | 2018-10-09 | 2019-03-08 | 杭州中恒云能源互联网技术有限公司 | 一种电动汽车动力电池剩余电量估算方法 |
Non-Patent Citations (2)
Title |
---|
State of charge estimation for Li-ion batteries using neural network modeling and unscented Kalman filter-based error cancellation;Wei He等;《Electrical Power and Energy Systems》;20140618;第783-791页 * |
基于NARX神经网络的光伏发电功率预测研究;付青等;《电气传动》;20161231;第46卷(第4期);第42-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110133507A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110133507B (zh) | 一种基于narx-ukf算法的电池剩余电量估计方法 | |
Deng et al. | Data-driven state of charge estimation for lithium-ion battery packs based on Gaussian process regression | |
CN107957562B (zh) | 一种锂离子电池剩余寿命在线预测方法 | |
CN105183994B (zh) | 一种基于改进的i-elm的动力电池soc预测方法及装置 | |
CN109061506A (zh) | 基于神经网络优化ekf的锂离子动力电池soc估计方法 | |
CN112241608A (zh) | 一种基于lstm网络和迁移学习的锂电池寿命预测方法 | |
CN113064093B (zh) | 储能电池荷电状态与健康状态联合估算方法及系统 | |
CN111505506A (zh) | 一种多尺度卡尔曼滤波与无迹卡尔曼滤波融合的电池soc估算方法 | |
CN113253116A (zh) | 锂离子电池荷电状态估计方法、存储介质 | |
CN105974327A (zh) | 一种基于神经网络和ukf的锂电池组soc预测方法 | |
CN111208438B (zh) | 基于神经网络与无迹卡尔曼滤波器的锂电子电池剩余电量与传感器偏差协同估计的方法 | |
CN112434463B (zh) | 一种车辆复合电源能量管理系统 | |
CN112989690B (zh) | 一种混合动力汽车锂电池多时间尺度荷电状态估计方法 | |
CN112557907A (zh) | 一种基于gru-rnn的电动汽车锂离子电池的soc估计方法 | |
Li et al. | A novel state estimation approach based on adaptive unscented Kalman filter for electric vehicles | |
CN112630659A (zh) | 一种基于改进bp-ekf算法的锂电池soc估算方法 | |
CN112163372B (zh) | 一种动力电池的soc估算方法 | |
CN112269133B (zh) | 一种基于预充电路模型参数识别的soc估计方法 | |
CN112782594B (zh) | 考虑内阻的数据驱动算法估算锂电池soc的方法 | |
CN115308606B (zh) | 一种基于邻近特征的锂离子电池健康状态估计方法 | |
CN116298936A (zh) | 在不完全电压范围内的锂离子电池智能健康状态预测方法 | |
CN114861545A (zh) | 基于rnn神经网络与多参数约束的锂电池sop在线估算方法 | |
CN114814591A (zh) | 一种锂电池soe估算方法、装置及系统 | |
CN110232432B (zh) | 一种基于人工生命模型的锂电池组soc预测方法 | |
CN112946480B (zh) | 一种提高soc估计实时性的锂电池电路模型简化方法 |
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 |