[发明专利]一种基于灰色扩展卡尔曼的锂离子电池SOC估算方法在审

专利信息
申请号: 201610372662.0 申请日: 2016-05-30
公开(公告)号: CN105842633A 公开(公告)日: 2016-08-10
发明(设计)人: 潘海鸿;吕治强;陈琳 申请(专利权)人: 广西大学
主分类号: G01R31/36 分类号: G01R31/36
代理公司: 暂无信息 代理人: 暂无信息
地址: 530004 广西壮族自治*** 国省代码: 广西;45
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于灰色扩展卡尔曼滤波算法用于估算锂离子电池SOC的方法,包括步骤:首先通过灰色预测模型,预测当前时刻电池模型极化电压和SOC状态量的先验估算值,替代扩展卡尔曼滤波算法中Jacobian矩阵;然后利用扩展卡尔曼滤波算法通过观测值对先验估算值进行更新、修正得当前时刻锂离子电池SOC估算值。本发明为电动汽车电池管理系统提供了一种锂离子电池SOC估算方法,能提高锂离子电池SOC估算精度。
搜索关键词: 一种 基于 灰色 扩展 卡尔 锂离子电池 soc 估算 方法
【主权项】:
一种基于灰色扩展卡尔曼算法的锂离子电池SOC估算方法,将灰色预测模型和扩展卡尔曼滤波算法融合,通过灰色预测模型预测当前时刻锂离子电池极化电压和SOC状态先验估算值,再通过观测值对其进行更新、修正,进行锂离子电池SOC估算,其特征在于,基于灰色扩展卡尔曼算法估算锂离子电池SOC方法的包含以下步骤:(1)基于锂离子电池一阶RC模型,通过建立锂离子电池的数学模型,得到电池系统状态空间方程和系统观测方程,具体为:设电池荷电状态SOC和RC网络极化电压U1为系统状态量,电池端电压Ut为观测量,根据SOC定义,以安时法和锂离子电池一阶RC模型建立锂离子电池各时刻极化电压和SOC状态空间方程(1)和系统观测方程(2):<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>S</mi><mi>O</mi><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>e</mi><mfrac><mrow><mo>-</mo><mi>T</mi></mrow><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac></msup></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>S</mi><mi>O</mi><mi>C</mi><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msup><mi>e</mi><mfrac><mrow><mo>-</mo><mi>T</mi></mrow><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac></msup></mrow><mo>)</mo></mrow><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mfrac><mrow><mi>&eta;</mi><mi>T</mi></mrow><msub><mi>C</mi><mi>n</mi></msub></mfrac></mrow></mtd></mtr></mtable></mfenced><msub><mi>I</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>R1表示RC网络部分极化内阻,C1表示RC网络部分极化电容,η为容量系数,T为系统采样时间间隔,Cn为电池标称容量,IL表示电流,w为服从高斯分布的白噪声,均值为0,协方差为Qk,k为正整数,k=1,2,3···;<mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>dU</mi><mrow><mi>o</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mi>S</mi><mi>O</mi><mi>C</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>S</mi><mi>O</mi><mi>C</mi></mrow></mfrac></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>U</mi><mn>1</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><mi>S</mi><mi>O</mi><mi>C</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>R</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>I</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>Uoc表示锂离子电池开路电压,R0表示锂离子电池欧姆内阻部分,v为服从高斯分布的白噪声,均值为0,协方差为Rk,v和w相互独立;将方程(1)和方程(2)改写通用表达式为得k时刻电池系统系数矩阵:Dk=R0(k);系统状态量表达式电流u(k)=IL(k),模型估算端电压(2)进行灰色扩展卡尔曼滤波算法设计:首先,算法初始化;给定算法初始值:x(0),P0,Q0和初始时刻v的协方差R0,其中P0为极化电压和SOC状态估算误差协方差矩阵Pk初始时刻时值;其次,预测电池系统状态量的先验估计值;根据灰色预测模型预测k时刻系统状态量先验估算结果,如式(3)所示:<mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>式中是由k时刻的前n个时刻系统量后验估算值灰色预测得来,是当前时刻系统状态量的先验估算值,系统状态估算误差协方差由式(4)计算得来:<mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>A</mi><mi>k</mi></msub><msubsup><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>+</mo></msubsup><msubsup><mi>A</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>最后,更新估算结果和估算误差协方差;利用k时刻系统观测量Y(k)更新系统状态估算结果和状态估算误差协方差,首先通过式(5)计算系统测量余量:<mrow><msub><mi>e</mi><mi>k</mi></msub><mo>=</mo><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>C</mi><mi>k</mi></msub><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>D</mi><mi>k</mi></msub><msub><mi>I</mi><mi>L</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>其中Y(k)为电池实测端电压;ek为k时刻系统测量余量;最优卡尔曼增益Kk通过式(6)计算:<mrow><msub><mi>K</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>k</mi></msub><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>利用式(7)和式(8)更新预测结果:<mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi><mo>+</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mi>k</mi></msub><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>式中是k时刻电池系统状态量更新后估算结果,即后验估算值。同时通过式(8)对该值误差协方差矩阵进行更新:<mrow><msubsup><mi>P</mi><mi>k</mi><mo>+</mo></msubsup><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>K</mi><mi>k</mi></msub><msub><mi>C</mi><mi>k</mi></msub><mo>)</mo></mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>式中为误差协方差矩阵,I为单位矩阵。(3)锂离子电池SOC估算:在该算法的实际应用中,首先利用扩展卡尔曼算法使电池SOC估算收敛后再启用灰色预测模型估算电池SOC;以k时刻前n个系统状态量的后验估算值通过灰色预测模型预测k时刻系统状态先验估计值再依次以式(4)至式(8)分别计算ek,Kk根据中的SOC(k)即可得k时刻电池SOC值;在k+1时刻时,根据采集到的电池端电压Ut(k+1)和电流IL(k+1)辨识出电池模型参数R0(k+1),R1(k+1),C1(k+1)并更新状态空间方程(1)和观测方程(2)中的Ak+1,Bk+1,Ck+1,Dk+1,重复上述预测和更新的步骤,完成电池SOC估算。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于广西大学,未经广西大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610372662.0/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top