首页 > 专利 > 杭州电子科技大学 > 一种基于NARX-UKF算法的电池剩余电量估计方法专利详情

一种基于NARX-UKF算法的电池剩余电量估计方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-04-23
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2019-09-10
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2021-09-14
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-04-23
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910329199.5 申请日 2019-04-23
公开/公告号 CN110133507B 公开/公告日 2021-09-14
授权日 2021-09-14 预估到期日 2039-04-23
申请年 2019年 公开/公告年 2021年
缴费截止日
分类号 G01R31/367G01R31/387G01R31/388 主分类号 G01R31/367
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 4
权利要求数量 5 非专利引证数量 1
引用专利数量 1 被引证专利数量 0
非专利引证 1、2006.08.16CN 105629175 A,2016.06.01Wei He等.State of charge estimationfor Li-ion batteries using neural networkmodeling and unscented Kalman filter-based error cancellation《.ElectricalPower and Energy Systems》.2014,付青等.基于NARX神经网络的光伏发电功率预测研究《.电气传动》.2016,第46卷(第4期),;
引用专利 EP1691209A 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 高明裕、秦潇涵、何志伟、朱晓帅、胡燕华 第一发明人 高明裕
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
朱月芬
摘要
本发明公开了一种基于NARX‑UKF算法的电池剩余电量估计方法。现有的各种方法存在各种问题。本发明方法首先测量电池在不同条件下的电压和电流,然后对测量数据进行预处理,再搭建NARX‑UKF网络,把处理后的电压电流训练数据输入NARX网络中来训练网络,达到训练目标后,输入测试数据,NARX网络的输入结果即为SOC的估计值同时也是UKF模型的输入量,该估计值经过UKF模型后,得到的状态更新值即为当前时刻所估计得到的电池剩余电量。本发明的不需要建立电池模型,只需要常见的可测量量便可以快速精确的进行电池剩余电量估计。有着模型训练速度快,所需参数少,估计精度高的优点。
  • 摘要附图
    一种基于NARX-UKF算法的电池剩余电量估计方法
  • 说明书附图:图1
    一种基于NARX-UKF算法的电池剩余电量估计方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2021-09-14 授权
2 2019-09-10 实质审查的生效 IPC(主分类): G01R 31/367 专利申请号: 201910329199.5 申请日: 2019.04.23
3 2019-08-16 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
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;
(1‑2).根据最小二乘方法拟合出Qi与Ci间的二次曲线关系,即在最小均方误差准则下求出同时满足 a、b、c为最优系数;
(1‑3).在放电电流为ik时,对应的库伦效率ηi为: 对于同一类型的电池最优系数a、b、c只确定一次,确定后作为已知常数直接用于所有同类型电池的剩余电量估计;
步骤(2).处理测量数据:
对电压电流进行标准化处理,处理后的数据x的范围在[‑1,1]之间,处理方法如下:
xmin和xmax分别为待标准化数据中的最小值和最大值;
步骤(3).搭建一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为R的NARX网络的数学模型,表示为: 其中,
fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值;
步骤(4).将步骤(2)得到的数据输入步骤(3)得到的网络中,对网络进行训练,训练的性能函数为: ym、tm分别是第m个测量点时的网络SOC预测值和SOC测
量值;
步骤(5).输入测试集,测试网络性能;
步骤(6).为NARX网络增加UKF模型,建立如下状态方程和测量方程:
Ek=SOCk+vk;其中,ik是k时刻的电流值,SOCk第k时刻的
SOC状态估计值,wk和vk是过程噪声和测量噪声,Ek是NARX网络在k时刻的估计值;UKF的输出即为SOC估计值;UKF修正估计值的方法如下:
(6‑1).系统初始化:初始化维度为K的状态 和协方差矩阵p0:
其中,x0为初始时刻的状态量;
(6‑2).计算k‑1时刻的2K+1个Sigma点集 对状态进行无迹变换并计算均值的权值和协方差的权值
其中, 和Pk‑1分别是K维状态变量的均值和方差, 表示矩阵
2
的第i列,λ=α(K+κ)‑K为比例缩放参数,α为控制采样点分布密集程度的参数,β为状态分布参数;
(6‑3).计算状态预测值,通过状态方程预测下一时刻的Sigma点: 过UT后的权重系数加权求得状态预测均值
计算状态预测协方差
(6‑4).测量更新:
通过测量方程预测下一时刻的Sigma点
过UT的权重系数加权求得测量预测均值
计算测量预测协方差Py,k:
预测和状态预测的交叉协方差矩阵Pxy,k:
‑1
卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k) ;
状态 及状态协方差矩阵Pk更新:

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算法的电池剩余电量估计方法,其特征在于:
‑4
步骤(6)中控制采样点分布密集程度的参数α取值范围为:1e ≤α≤1。

5.如权利要求1所述的一种基于NARX‑UKF算法的电池剩余电量估计方法,其特征在于:
步骤(6)中κ的取值需保证矩阵半正定性,通常取值为3‑K。
说明书

技术领域

[0001] 本发明属于电池技术领域,提出了一种基于NARX(基于非线性自适应回归神经网络)‑ UKF(无迹卡尔曼滤波器)算法的电池剩余电量估计的解决方案。

背景技术

[0002] 近些年来,随着经济的发展,汽车保有量的不断增加,环境问题也日益严重。使用高效环保的代替能源是解决污染问题和传统能源衰竭问题的有效途径之一,具有低噪声、低污染、低排放和能量利用率高等优点,因此新能源汽车得到了各国的大力发展。电池管理系统(BMS,Battery Management System)是电动汽电池的关键技术,电池管理系统的发展对电动汽车具有着相当重大的意义,BMS作为电池系统的大脑,控制着电池的充放电过程,避免电池出现过充过放的问题,实现对电池的有效的管理和监控,降低电池使用过程中存在的风险,最终达到延长电池使用寿命和对电池的最有效的利用的目的。
[0003] 荷电状态估计(SOC,State Of Charge)是电池管理系统需要计算的主要参数之一,SOC 可以直接反应电池剩余电量,定义为电池所剩电量和电池额定容量的比值,充当燃油汽车中燃油表。SOC估计对于电池系统的正常运转是必不可少的,是电动汽车续航行驶和防止电池过度充电过度放电的主要依据,对电动汽车的整体性能还有电池的寿命效率有重要的意义。电池的SOC参数具有十分重要的意义,但是不可以直接测量。电池是一个高度非线性的系统,影响SOC的因素很多,使得SOC估计具有一定的难度。
[0004] 目前,已经被提出的SOC估计方法主要有:安时积分法、开路电压法、卡尔曼、神经网络法等。安时积分法在实际应用条件下有噪声误差会造成计较大的累积误差,开路电压法需要电池长时间的静置,卡尔曼滤波法对电池建模的精确性要求较高,一般的神经网络法估计 SOC,估计结果的跳动较大。

发明内容

[0005] 本发明的目的就是克服已有技术的不足,为了满足较高精度的估计电池SOC值,提出了一种基于NARX‑UKF算法的电池剩余电量估计方法。
[0006] 本发明利用具有外部输入的非线性自回归网络(NARXNN,Nonlinear Autoregressive with Exogenous Inputs Neural Networks)来估计锂电池的SOC值,并且加入了无迹卡尔曼滤波器(UKF,Unscented Kalman Filter)来减小NARXNN模型的估计误差。具体步骤是:
[0007] 步骤(1).对特定的单节锂电池进行充放电实验,测量电池工作电流和工作电压,测量间隔为Δt为1s,并记录温度;使用安时积分法计算电池的SOC值,作为目标值:其中SOC(n)是电池在第n个测量点时的SOC值,η为库伦
效率,I为电流值的大小,充电时为负值,放电时为正值,QN是电池的额定容量。
[0008] 库伦电效率η按如下方法确定:
[0009] (1‑1).将以1/30倍额定电流放电时的额定容量为Qn的完全充满电的电池以不同放电速率Ci恒流放电N次,0<Ci≤C,1≤i≤N,N≥10,放空电池电量所需要的时间和放电电流值的乘积即为相应放电速率下的电池总电量Qi;
[0010] (1‑2).根据最小二乘方法拟合出Qi与Ci间的二次曲线关系,即在最小均方误差准则下求出同时满足 +bCi+c,a、b、c为最优系数;
[0011] (1‑3).在放电电流为ik时,对应的库伦效率ηi为: 对于同一类型的电池最优系数a、b、c只需确定一次,确定后作为已知常数直接用于所有同类型电池的剩余电量估计。
[0012] 步骤(2).处理测量数据:
[0013] 对电压电流进行标准化处理,处理后的数据x的范围在[‑1,1]之间,处理方法如下:
[0014] xmin和xmax分别为待标准化数据中的最小 值和最大 值。
[0015] 步骤(3).搭建NARXNN模型,输入延时步长设置为5~15,输出延时步长设置为2~5,神经元节点数设置为15;一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为 R的NARX网络的数学模型表示为:
[0016]
[0017] 其中,fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值。
[0018] 步骤(4).将步骤(2)得到的数据输入步骤(3)得到的网络中,对网络进行训练,训练的性能函数为: ym、tm分别是第m个测量点时的网络SOC预测值和 SOC测量值。网络的MSE训练目标在0℃、25℃和45℃下分别设置为8e‑06、1.1e‑05和 1.9e‑
05。
[0019] 步骤(5).输入测试集,测试网络性能。
[0020] 步骤(6).为NARX网络增加UKF模型,建立如下状态方程和测量方程:
[0021] Ek=SOCk+vk;其中,ik是k时刻的电流值,SOCk第k时刻的SOC状态估计值,wk和vk是过程噪声和测量噪声,Ek是NARX网络在k时刻的估计值。
[0022] UKF的输出即为SOC估计值。UKF修正估计值的方法如下:
[0023] (6‑1).系统初始化:初始化维度为K的状态 和协方差矩阵p0:
[0024]
[0025] 其中,x0为初始时刻的状态量。
[0026] (6‑2).计算k‑1时刻的2K+1个Sigma点集 对状态进行无迹变换并计算均值的权值 和协方差的权值
[0027]
[0028]
[0029] 其中, 和Pk‑1分别是K维状态变量的均值和方差, 表示矩阵2
的第i列,λ=α(K+κ)‑K为比例缩放参数,α为控制采样点分布密集程度的参‑4
数,取值范围为:1e ≤α≤1,κ的取值需保证矩阵半正定性,通常取值为3‑K,β为状态分布参数。
[0030] (6‑3).计算状态预测值,通过状态方程预测下一时刻的Sigma点: 过UT 后的权重系数加权求得状态预测均值
[0031] 计算状态预测协方差
[0032] (6‑4).测量更新:
[0033] 通过测量方程预测下一时刻的Sigma点
[0034] 过UT的权重系数加权求得测量预测均值
[0035] 计算测量预测协方差Py,k:
[0036] 预测和状态预测的交叉协方差矩阵Pxy,k:
[0037] 卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k)‑1;
[0038] 状态 及状态协方差矩阵Pk更新:
[0039] 本发明的训练集和测试集均采用工况下的放电数据,贴合电动汽车实际行驶时的电流电压变化,可以方便地进行电池SOC的快速估计,不需要对电池建模,只需要易于测量的量,该方法收敛速度快,模型训练时间只需要十几秒,估计精度高、结果误差小、计算量小、波动小,所需参数数目少,便于移植到硬件平台中,而且适用于各种电池SOC的快速估计。附图说明:
[0040] 图1为本发明系统的流程图。

实施方案

[0041] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
[0042] 参见图1,一种基于NARX‑UKF算法的电池剩余电量估计方法,具体步骤是:
[0043] 步骤(1).对特定的单节锂电池进行充放电实验,测量电池工作电流和工作电压,测量间隔为Δt为1s,并记录温度;使用安时积分法计算电池的SOC值,作为目标值:其中SOC(n)是电池在第n个测量点时的SOC值,η为库伦
效率,即充放电效率,I为电流值的大小,充电时为负值,放电时为正值,QN是电池的额定容量。
[0044] 充放电效率η按如下方法确定:
[0045] (1‑1).将以1/30倍额定电流放电时的额定容量为Qn的完全充满电的电池以不同放电速率Ci恒流放电N次,0<Ci≤C,1≤i≤N,N≥10,放空电池电量所需要的时间和放电电流值的乘积即为相应放电速率下的电池总电量Qi;
[0046] (1‑2).根据最小二乘方法拟合出Qi与Ci间的二次曲线关系,即在最小均方误差准则下求出同时满足 a、b、c为最优系数;
[0047] (1‑3).在放电电流为ik时,对应的库伦效率ηi为: 对于同一类型的电池最优系数a、b、c只需确定一次,确定后作为已知常数直接用于所有同类型电池的剩余电量估计。
[0048] 步骤(2).处理测量数据:
[0049] 对电压电流进行标准化处理,处理后的数据x的范围在[‑1,1]之间,处理方法如下:
[0050] xmin和xmax分别为待标准化数据中的最小 值和最大 值。
[0051] 步骤(3).搭建NARXNN模型,输入延时步长设置为5~15,输出延时步长设置为2~5,神经元节点数设置为15;一个具有H个隐藏层节点、输入延时步数为L、输出延时步数为 R,在本实施例中,0℃、25℃和45℃下设置输入延时步长分别为8、10和8,输出延时步长分别为3、4和3,隐藏层神经元个数设为15,输出层和隐藏层的激活函数分别为tansig 和purelin函数。搭建的NARX网络的数学模型表示为:
[0052]
[0053] 其中,fo(·)和fh(·)分别是输出层和隐藏层的非线性函数,bo和bh分别是输出层和隐藏层的阈值,who、wlh和wrh分别是第h个隐藏层节点到输出层节点、第l个输入层节点到第h个隐藏层节点和第r个输出反馈层节点到第h个隐藏层节点的权值。
[0054] 步骤(4).将步骤(2)得到的数据输入步骤(3)得到的网络中,对网络进行训练,利用 Levenberg‑Marquradt算法训练神经网络,保证目标函数均方误差(MSE,Mean Squared Error)最小,训练的性能函数为: ym、tm分别是第m个测量点时的网络SOC预测值和SOC测量值。为了防止过拟合,网络的MSE训练目标在0℃、25℃和 45℃下分别设置为8e‑06、1.1e‑05和1.9e‑05。
[0055] 步骤(5).输入测试集,测试网络性能。
[0056] 步骤(6) .为NARX网络增加UKF模型,建立如下状态方程和测量方程:Ek=SOCk+vk;其中,ik是k时刻的电流值,SOCk第k时刻的SOC
状态估计值,wk和vk是过程噪声和测量噪声,Ek是NARX网络在k时刻的估计值。
[0057] UKF的输出即为SOC估计值。UKF修正估计值的方法如下:
[0058] (6‑1).系统初始化:初始化维度为K的状态 和协方差矩阵p0:
[0059]
[0060] 其中,x0为初始时刻的状态量。
[0061] (6‑2).计算k‑1时刻的2K+1个Sigma点集 (又称之为采样点),对状态进行无迹变换并计算均值的权值 和协方差的权值
[0062]
[0063]
[0064] 其中, 和Pk‑1分别是K维状态变量的均值和方差, 表示矩阵2
的第i列,λ=α(K+κ)‑K为比例缩放参数,α为控制采样点分布密集程度的参‑4
数,取值范围为:1e ≤α≤1,κ的取值需保证矩阵半正定性,通常取值为3‑K。β为状态分布参数,本实施例取值为2。
[0065] (6‑3).计算状态预测值,通过状态方程预测下一时刻的Sigma点: 过UT 后的权重系数加权求得状态预测均值
[0066] 计算状态预测协方差
[0067] (6‑4).测量更新:
[0068] 通过测量方程预测下一时刻的Sigma点
[0069] 过UT的权重系数加权求得测量预测均值
[0070] 计算测量预测协方差Py,k:
[0071] 预测和状态预测的交叉协方差矩阵Pxy,k:
[0072] 卡尔曼增益矩阵Kk:Kk=Pxy,k(Py,k)‑1;
[0073] 状态 及状态协方差矩阵Pk更新:
[0074] 以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号