第24卷 第20期 Vo1.24 No.20 电子设计工程 Electronic Desizn Engj 2016年1O月 Oct.2016 锂电池荷电状态估算Matlab仿真研究 涂涛,钟其水,李波 (电子科技大学航空航天学院,四川成都611731) 摘要:锂电池的剩余电量(SOC)估算不仅可以作为电动汽车续航里程的参考值,而且可以为电动汽车的能量管理策 略提供依据,具有重要意义。本文以法国SAFT公司生产的额定容量为6 AH,额定电压为l0.68 V的锂离子电池包为 研究对象,通过使用Matlab和Advisor等仿真软件,研究了福克斯电动汽车行驶在UDDS工况下,采用安时法和扩展 卡尔曼滤波算法结合估算锂电池的剩余电量。仿真结果证明,该方法有效提高了锂电池的SOC估算精度.电池SOC 估算误差在5%以内。 关键词:锂离子电池;SOC估算;安时法;扩展卡尔曼滤波算法;Matlab仿真 中图分类号:TN98 文献标识码:A 文章编号:1674—6236f2016)20—0129—04 The Matlab simulation research f0r the SOC estimation of lithium.ion battery 1’U Tao,ZHONG Qi—shui,LI Bo (School ofAeronautics and Astronautcis,University fEloectronic Science and Technology f oChina,Chengdu 61 1731,China) Abstract:The state of charge(SOC)of lithium-ion battery is the key pammeteI|0.r endurance mileage and energy management strategy of electric vehicles.In this paper,both AH algorithm and Extended Kalman Filtering algorithm have been used to estimate the SOC of he SAFY tlithium battery with capacity of 6AHand nominal voltage of 10.68V under the working condition of UDDS.Based on Matlab and Advisor,simulation results show that the proposed methods can effectively increase the SOC estimation accuracy and the error of SOC estimation is less than 5%. Key words:lithium battery;SOC estimation;AH algorithm;EKF algorithm;Matlab simulation 当前,汽车行业快速发展。初步判断2020年前,我国乘用 车市场仍将处于快速发展时期,增长率相当于GDP增长率的 1.5倍左右[11。然而传统汽车的发展面临着化石能源枯竭以及 多元线性回归算法实时在线辨识电池等效电路模型的参数 值,然后根据测得的电池端电压和充放电电流来计算得到电 池的开路电压值,从而获得电池的SOC值。该方法SOC估算 精度主要取决于电池等效模型的阶数和参数的辨识精度。卡 环境污染等问题,因此,发展使用清洁能源的电动汽车成了当 下的热点。电池管理系统是电动汽车最关键的部件之一.对于 尔曼滤波算法的核心思想是用观测值来进一步修正预测值, 其精度主要取决于电池等效电路模型的精度。文献[5】研究了 采用扩展卡尔曼滤波算法估算电池SOC值,取得了较高的 SOC估算精度。机器学习的方法主要是人工神经网络算法, BP网络模型是目前应用比较成熟的一种神经网络模型,它 能够很好的对复杂系统建模.并且具有多个输入和非线性特 点 。因此,BP神经网络法很适合用来估算电池的SOC值。该 优化电池组性能,延长电池组寿命,提高电池组能量使用效 率,增加续航里程,确保电动车的行驶安全等都具有十分重要 的意义[21,而电池的荷电状态(SOC)估算是电池管理系统的关 键技术之一,精确的SOC估算可有效防止电池过充过放,提高 电池使用寿命,提高能量利用率,具有重要意义。目前,国内外 在锂离子电池soc估算方面已取得大量研究成果.主要包括 安时积分法、开路电压法、扩展卡尔曼滤波法、机器学习等。 安时积分法是根据电池充入和放出的电量来计算电池 的SOC值。该方法没有考虑电池内阻、电池温度、电池老化等 因素的影响。文献【3】研究了采用安时法估算电池的soc值, 实验结果表明采用安时法估算电池SOC值在短时间内精度 较高,但是其估算误差会逐渐积累,以至到无法接受的地步。 开路电压法是依据电池SOC值和开路电压间的一一对应关 方法可以应用于各种锂电池.但是估计精度受采集数据精度 和训练方法影响较大 。 综合考虑安时积分法和扩展卡尔曼滤波算法的优缺点, 文中采用安时法和扩展卡尔曼滤波法结合估算电池的剩余 电量。通过在电池放电前期和后期采用安时法估算电池的 SOC值,在电池中间放电阶段采用扩展卡尔曼滤波算法估算 电池SOC值.并采用扩展卡尔曼滤波算法修正安时法SOC 估算的误差,有效提高了电池SOC估算精度,仿真结果表明, UDDS工况下。该算法SOC估算误差在5%以内。 系来估算电池SOC值。文献f4】通过Z域变换,采用最小二乘 收稿日期:2016—03—10 稿件编号:201603129 基金项目:四川省科技支撑计 ̄'l(2014GZ0079,2015GZ0129,2016GZ0020);成都市科技惠民项目(2014HM0100101SF) 作者简介:涂涛(199l一),男,江西上饶人,硕士研究生。研究方向:电池管理系统。 .129- 《电子设计工程)2016年第20期 1 安时法和扩展卡尔曼滤波法SoC估算原理 观测方程:Y(k+1)=C( +1)X(k+1)+D( +1 (Jj}+1) (k+1) 1.1 S0C的定义 电池的SOC是电池的一个非常重要的参数,目前SOC 的定义还没有统一。不同机构有不同的定义,如美国先进电 池联合会将SOC定义为 :电池在一定放电倍率下,剩余容量 与相同条件下额定容量的比值,即: s0c= ×100% (1) 式(1)中,剩余容量是电池在一定条件下从当前状态放 电至终止所放出的总电量,额定容量是指电池在一定条件下 充满电能够放出的总电量。 日本本田汽车公司将SOC修正定义为: .s0 硪 ×100%(2) 式(2)中的容量衰减因子是考虑到温度、充放电倍率、自 放电、电池老化等对电池总容量的影响。 起亚汽车公司将SOC定义如下式所示。 |soc_ ×100% (3) 文中采用新电池进行相应的实验研究,可以忽略容量衰 减因子的影响。因此采用式(1)作为SOC的定义方式。 1.2安时法SOC估算原理 安时法又称作库伦计数法,是当前比较常用的估计算法0o41, 安时法的基本原理是将电池看作一个黑箱子。然后计算电池充入 和放出的电量来获得电池的SOC值。在对温度、电池老化、充放电 效率、自放电有较好补偿时,短时间内,安时法有很好的估计效果, 安时法估算SOC的计算如式(4)所示。 r‘ S(£)=Is(如)+I pxl(t)dt/Q0 (4) 式(4)中,S(to)是电池在to时刻的剩余电量,I(t)是电池 充放电电流,充电电流为正,放电电流为负,()0是电池的额定 容量,t时刻的电量等于t。时刻的电量加上电池从to到t充放 电电量,P是充放电效率。安时法原理简单,易于实践,在保证 电流测量精度和初始电池容量已知的情况下.安时积分法具 有非常高的精度和计算简单的优点【I21。 1.3扩展卡尔曼滤波法SoC估算原理 卡尔曼滤波算法是一种数学方法,主要思想是用当前时 刻的观测值修正当前时刻状态估计值,使得状态估计值更接 近真实值。卡尔曼滤波的两个关键是状态一步预测值和新息. 状态一步预测可以对下一个状态的SOC进行估算.新息可以 对一步预测值进行校正【 q。但是卡尔曼滤波只适用于线性模 型,而电池本身是一个非常复杂的非线性系统,因此需要将非 线性系统线性化。将非线性系统线性化进行估算的方法即扩 展卡尔曼滤波算法[15-161。通过建立扩展卡尔曼滤波算法的状态 方程和观测方程,实现电池的SOC估算。由于实际处理器处 理的都是离散时间信号,因此文中研究离散型的卡尔曼滤波。 离散系统的状态方程和观测方程如下所示。 状态方程: (五十1);A( )x(k)柏(后)U(k) ( ) (5) 一】30- (6) 式(5)状态方程中, (|j}+1)是k+l时刻的状态向量,u(k) 是k时刻的输入向量。A(k)是k时刻的系统矩阵,B(k)是k 时刻的输入矩阵。w(k)是k时刻的系统噪声向量。式(6)观测 方程中,Y(k+1)是k+l时刻的观测向量,C(k+1)是k+l时刻 的输出矩阵,D(k+1)是k+l时刻的直接传输矩阵, (k+1)是 k+l时刻的测量噪声向量。这里w(k)和 (k)都是零均值的白 噪声,彼此之间相互独立,并与状态向量x(k)不相关。w(k)、 x(k)、 (k)的统计特性如下式所示。 『目t ( )=0,印( )】= ,Cov[w(k), ( ) (J})口(J}) {Co ̄[w(|i}), U)l= E (.i}) (,) }=Q( )6(J} ),Cov[X(|i}), U))=0(7) 【Cov[v(k),vq)l=E[v(k)v(i)q=R(k ( ),Co,IX(|j}),vq)t=o 文中选择SOC值作为状态量,电池端电压徘为观测值。 电流,作为输人值.建立电池的状态方程和观测方程如下所示: 状态方程:S(k+l/k)=S(|i})+,( )*t/Qo+w(k) (8) 观测方程:U(k+1) (.j}+l )】棚I,( +1) (k+1) (9) 式(9)中,(S(|j}+1 ))是开路电压和电池剩余电量SOC 的函数,因为是非线性函数,对其进行泰勒级数展开并忽略二 次以上的项,以实现方程的线性化,可得扩展卡尔曼滤波的观 测方程为: u(k+1);§(jc+1 )×. Is 撕n 1 +1)十,( +UR ̄+v(k+1)(1o) 根据卡尔曼滤波原理推导可得: 状态一步预测: S(k+l/k)=.s( ) ( )xAt/Q0 (11) 观测值一步预测: (k+1)=‘§( +1屈)×掣Is; —.R1,( +1)+,( +1)R (12) 状态一步预测误差自相关矩阵: P( +l )=P(J})+Q( ) (13) 卡尔曼增益: K( +1): ( +l/k)x ^)】×【立 xP(k+1 × ∞ 掣l “ 1】一 (14) 状态最优估计值: S(k+1)=S(k+l/k)+ ( +1)【u( +1)-U(k+1)】 (15) 状态最优估计误差自相关矩阵: P( +1):(1-K(k+1) } {5= ( +l ))P( +1, ) (16) 扩展卡尔曼滤波算法估算电池SOC的具体步骤为:根据 电池初始剩余电量估算值S(0)和充放电电流1(o),可以得到 状态一步预测值 (1,0),端电压一步预测值D(1);根据初始状 态估算误差自相关值P(O)计算得到一步预测误差自相关值 P(1/0),然后计算卡尔曼增益值 (1);在下一时刻测得电池 端电压U(1)和电流,(1)观测值,最终计算最优估计值S(1)以 涂涛。等锂电池荷电状态估算Matlab仿真研究 法SOC估算误差会逐渐积累变大,最终达到不可接受的程度。 在同样的工况下,采用扩展卡尔曼滤波算法估算电池的 及最优估计值的误差自相关值P(1)。以此类推可实时估算电 池SOC值。电池SOC值估算流程如图1所示。 ④ 面 t SOC。得到的Matlab仿真结果如图4所示。 UDDS ̄况下扩展卡尔曼滤波算法S0C估算结果 ● 时刻s1l修正值 l0c l 图1扩展卡尔曼滤波SOC估算流程 上 2 Maflab仿真及结果分析 在Advisor工具中,设定研究对象为Focus_in纯电动汽 f/s 图4扩展卡尔曼滤波SOC估算 由图4可知,采用扩展卡尔曼滤波法进行电池SOC估算 在放电前期和后期精度较低,在放电中间阶段精度较高。这主 要是由于扩展卡尔曼滤波算法的SOC估算精度主要由电池 车。Focus_in纯电动汽车的动力电池组是由25个额定容量为 6 AH、额定电压为1O.68 V的电池包串连而成,该电池包由法 国SAFT公司生产。该电池包的开路电压和SOC之间的关系 曲线拟合结果如图2所示。 模型的精度决定。而电池模型参数在放电前期和后期变化较 大,在放电中间阶段,电池模型参数基本恒定。因此,扩展卡尔 曼滤波算法估算电池SOC值在放电中间阶段精度较高。在放 电前期和后期精度较差。 根据安时法和扩展卡尔曼滤波法SOC估算的以上特性。 文中结合使用两种方法进行电池的SOC估算。具体方法为: 在电池放电前期,采用安时法估算电池的SOC值,经过200 秒后。采用扩展卡尔曼滤波算法估算电池的SOC值,当电池 的SOC值低于25%时,重新使用安时法进行SOC估算。在 UDDS工况下,两种方法结合使用的SOC估算MaⅡab仿真结 果如图5所示。 s0c 图2开路电压一SOC拟合结果 由图2可知。开路电压和SOC值的五次拟合精度最高。 因此.文中采用五次拟合的方法来获取电池充放电过程中 SOC对应的开路电压值。 在UDDS工况下。对安时法估算电池SOC进行Matlab仿 真实验.得到SOC估算结果如图3所示。 IJDDS3L况下安时法S0C估算结果 t/0 图5安时法一扩展卡尔曼滤波法结合SOC估算结果 估算误差如图6所示。 由图5和图6可知。前期采用安时法估算电池SOC,SOC 估算精度较高,但误差逐渐变大,中间阶段采用扩展卡尔曼滤 波算法进行SOC估算,可以修正安时法SOC估算的累积误 差.且扩展卡尔曼滤波算法进行SOC估算不存在累积误差, 后期采用安时法进行SOC估算,同样具有较高精度,但误差 图3安时法SOC估算 也会逐渐变大。最后再次用扩展卡尔曼滤波算法修正安时法 SOC估算的累积误差。整个放电过程中,电池SOC估算误差 一由图3可知,前期安时法SOC估算精度较高,但是该算 131一 《电子设计工程)2016年第20期 on Extended Kalman Filter[J].Applied Energy,2015,155 834-845. 【6]廖恩华.基于神经网络的电动汽车磷酸铁锂电池SOC估算 方法研究【D】.成都:电子科技大学,2011. g [7]Christian Fleischer,WladislawWaag,Hans—Martin Heyn et a1.On-line adaptive battery impedance parameterand state estimation considering physical principles in reduced order equivalent circuit battery models Part 1.Requirements, critical review of methods and modeling[J].Journal of Power £/ Sources,2014,260(2014):276—291. 图6安时法一扩展卡尔曼滤波法结合SOC估算误差 【8】周美兰,赵强,周永勤.改进g/PSO-BP神经网络估算磷酸 铁锂电池S0C 哈尔滨理工大学学报,2015,20(4):88-92. [9】Simon Schwunk,Nils Armbruster,Sebastian Straub.Particl— efilterforstateofchargeand stateofheal ̄estimationforlithiu— 始终保持在5%以内。 3结 论 安时法估算电池SOC在放电初始时刻精度较高,但是安 时法估算误差会随着时间逐渐累积,最终达到不可接受的程 度。扩展卡尔曼滤波算法SOC估算精度主要取决于所建立的 电池模型的精度,在电池高电量和低电量阶段,电池模型参数 变化较大,扩展卡尔曼滤波算法SOC估算精度较差。在电池 meiron phosphate batteries[J].Journal of Power Sources, 2013,239:705-710. 【101齐国光,李建民,郏航等.电动汽车电量计量技术的研究 【J】.清华大学学报:自然科学版,1997(3). [1 llSUN Feng—chun,RUI Xiong,HE Hong—wen.Estimation of state・・of-charge and state・-o ̄power capability of lithium--ion 放电的中间阶段,电池模型参数值变化小,扩展卡尔曼滤波算 法SOC估算精度较高。文K中采用安时法和扩展卡尔曼滤波 法结合估算电池SOC,在电池放电初期和后期采用安时法.在 电池放电中间阶段采用扩展卡尔曼滤波法,并采用扩展卡尔 battery considering varying health conditions[J】.Journal of Power Sources,2014,259(2014):166—176. [12]Piller S,Penn rM,Jossen A.Methods ofr state-of-charge deter- 曼滤波法的SOC估算值修正安时法的初值。Matlab仿真结果 表明,文中提出的方法有效提高了电池SOC估算精度。然而。 本文只研究了该算法在UDDS工况下的SOC估算效果.若要 将该算法运用于实践,将来仍需在电池充电过程中以及其它 工况下对该算法进行实验验证。 参考文献: [1】曹云平.电动车动力电源的发展现状[J].化工时刊,2001 (10):13—16. mination and their Applications[J].Journal of Power Source, 2oo1(96):l 13—120. 【13]Christian Fleischer,Christian Fleischeret,et a1.On—line ada— ptive battery impedance parameter and state estimation cons ̄ dering physical principles in reduced order equivalent Circuit 【JJ.Journal of Power Sources,2014,260:276—291. [14]KatlainGy rgy,AndrrsKelemen,Ldszl6Ddvid.Unscented kal_ man filtersandParticleFihermethodsfornonlinearstateestima- [2】Fengchun Sun,RuiXiong,Hongwen He.Estimation of State— of-Charge and State—of-Power Capability of Lithium.ion Bat— tion[J].Procedia Technology,2014,12:65-74. [15]RuiXiong,GONG Xian-zhi,Chunting Chris Mi.A robust state- o ̄eharge estimator for multiple types of lithium-・ion batteries—- tery Considering Varying Health Conditions[J].Journal of Power Sources,2014.259(2014):166—176 using adaptive extended kalmanfiher[J].Journal of Power Sources,2013,243:805—482. [3】王飞.动力电池运行状况测控技术的研究[D].武汉:武汉理 工大学自动化学院.2008. 【16]CHAO Hu,Byeng D.Youn,Jaesik Chung.A muhiscale  ̄amework wih exttended Kalman filter for lithium—ion 【4】张东华,马燕,陈思琪等.锂电池模型参数估计与荷电状态 估算研究[J],武汉理工大学学报,2015,37(2):179—182. [5】Perez G,Garmendia M,Reynaud JF,et a1.Enhanced clos— ed loop State of Charge estimator for lithium-ion batteries based battery SOC and capacity estimation[J】.Applied Enery,g 2012,92:694—704. 欢迎订阅2017年度《电子设计工程》(半月刊) 国内邮发代号:52—142 —.国际发行代号:M2996 订价:15.oo元/期360.O0元/年 132-