本发明涉及电动汽车电池参数估算技术领域,尤其涉及一种电池参数在线辨识方法及系统。
背景技术:
动力电池在新能源领域得到了广泛的应用,特别是在电动汽车上。为了实现对电池的精准管理,需对电池各种状态,例如,电池的剩余电量及电池健康状态等进行估算,而对电池状态的估算离不开电池参数,电池参数的准确性直接影响电池各种状态的估算精度。
电池参数通常包括电池开路电压、欧姆内阻、极化内阻及极化电容等参数目前,获取电池参数的方法主要有离线测试方法以及在线辨识方法两类。其中,离线测试方法通过专用的充放电测试设备,在特定工况下采集电池电压、电流及温度等数据,对测试数据进行处理后完成电池参数的标定。然而,该方法不能满足实时在线估算电池参数的需求,在实验室离线测试标定的电池参数初始值不能满足电池在使用过程中由于老化等原因引起的参数变化,亦即,该方法不能实现参数的实时在线跟踪。
现有的在线辨识方法主要为卡尔曼滤波法及最小二乘法。该类方法一般通过电池管理系统(Battery Management System,BMS)的采集模块对电池外特性参数进行采集,例如,电池实时电压、电流及温度等,然后通过BMS的运算单元在线估算电池模型参数。虽然,该类方法能够实现实时在线估算电池参数,但是选用的电池模型、参数方程以及计算复杂度差异较大,存在电池模型太复杂、运算量过大等缺点。
鉴于此,实有必要提供一种新的电池参数在线辨识方法及系统以克服上述缺陷。
技术实现要素:
本发明的目的是提供一种电池参数在线辨识方法及系统,所述电池参数在线辨识方法及系统能有效同时实现电池多个参数的在线辨识,且计算复杂度较低、精度较好。
为了实现上述目的,本发明提供一种电池参数在线辨识方法,其应用于电动汽车中,所述电动汽车包括为其提供动力的电池包,所述电池包包括若干单体电池;所述电池参数在线辨识方法包括如下步骤:
实时采集所述单体电池电压及电流;
建立电池电压及电流与各项待测电池参数的参数控制方程,其中,所述参数控制方程为
依据所述电压、电流并结合所述参数控制方程,利用遗忘递推最小二乘算法计算出θ1、θ2、θ3、θ4四个形式参数并输出;
依据所述θ1、θ2、θ3、θ4四个形式参数计算电池欧姆内阻Ro(k)、极化内阻Rc(k)、极化电容C(k)及开路电压Uocv(k)四个电池参数。
本发明还提供一种电池参数在线辨识系统,其应用于电动汽车中,所述电动汽车包括为其提供动力的电池包,所述电池包包括若干单体电池;所述电池参数在线辨识系统包括采集模块、参数控制方程模块、递推估算模块以及参数计算模块;所述采集模块用于实时采集所述单体电池电压及电流;所述参数控制方程模块用于建立所述电池电压及电流与各项待测电池参数的参数控制方程,其中,所述参数控制方程为所述递推估算模块用于依据所述电压、电流并结合所述参数控制方程,利用遗忘递推最小二乘算法计算出θ1、θ2、θ3、θ4四个形式参数并输出;所述参数计算模块用于依据所述递推估算模块输出的θ1、θ2、θ3、θ4四个形式参数计算电池欧姆内阻Ro(k)、极化内阻Rc(k)、极化电容C(k)及开路电压Uocv(k)四个电池参数。
本发明提供的电池参数在线辨识方法及系统,使用一种简单的电池模型以及参数控制方程,并结合遗忘递推最小二乘法算法实现电池参数的在线辨识。具有高度实时性、低计算复杂度及较高辨识精度等优点。
【附图说明】
图1为包含本发明提供的电池参数在线辨识系统的电动汽车的功能模块图。
图2为本发明提供的电池参数在线辨识方法的流程图。
图3为采集模块采集到的电压工况图。
图4为采集模块采集到的电流工况图。
图5为欧姆内阻Ro在线辨识结果图。
图6极化内阻Rc在线辨识结果图。
图7极化电容C在线辨识结果图。
图8开路电压Uocv在线辨识结果图。
【具体实施方式】
为了使本发明的目的、技术方案和有益技术效果更加清晰明白,以下结合附图和具体实施方式,对本发明进行进一步详细说明。应当理解的是,本说明书中描述的具体实施方式仅仅是为了解释本发明,并不是为了限定本发明。
请参阅图1,为包含本发明提供的电池参数在线辨识系统200的的电动汽车99的功能模块图。所述电动车99还包括电池包100、电池参数在线辨识系统200、存储器300以及处理器400。所述电池包100用于为所述电动汽车99提供动力。在本实施方式中,所述电池包100包括若干单体电池。所述电池参数在线辨识系统200包括采集模块10、参数控制方程模块20、递推估算模块30以及参数计算模块40。可以理解地,上述各功能模块可以软件程序的形式存储于存储器300中,并由处理器400执行。在替代实施例中,上述各功能模块也可为具有特定功能的硬件,例如,烧录有特定软件程序的芯片。
下面结合图2对上述各功能模块进行详细的介绍。
如图2所示,其为本发明实施例中电池参数在线辨识方法的流程图。所应说明的是,本发明的方法并不受限于下述步骤的顺序,且其他实施例中,本发明的方法可以只包括以下所述步骤的其中一部分,或者其中的部分步骤可以被删除。
S01,采集模块10用于实时采集单体电池电压U及电流Cur。
请再参阅图3及图4,为本实施方式中所述采集模块10采集到的电压、电流工况图。从电压工况图中可以看出,所述单体电池随着放电的进行电压逐渐减小。从电流工况图中可以看出,所述单体电池以7A的电流进行间歇性循环放电。具体地,实现电压、电流采集的方法有许多种,在本实施方式中,所述采集模块10以电池管理系统为例进行说明。电池管理系统的核心功能是对所述电池包100进行管理,而电池电压、电流采集是其基本功能之一,因此,任何一款电池管理系统均能实现电池电压、电流采集功能。在本实施方式中,单体电池电压记为U,电流记为Cur。可以理解,若所述采集模块10采集的电流为所述电池包100的总电流,则将所述总电流除以并联的单体电池数量即可得到单体电池的电流。
S02,所述参数控制方程模块20用于建立电池电压U及电流Cur与各项待测电池参数的参数控制方程。
具体地,所述参数控制方程如下:
其中,U(k)为第k时刻的电池端电压;Cur(k)为第k时刻的电池电流;为电流微分量,并且:为电压微分量,并且:θ1(k)为形式参数,并且θ1(k)=Ro(k)+Rc(k);θ2(k)为形式参数,并且θ2(k)=Ro(k)Rc(k)C(k);θ3(k)为形式参数,并且θ3(k)=-RC(k)C(k);θ4(k)为形式参数,并且θ4(k)=Uocv(k)。
进一步地,k表示当前时刻,k-1表示上一时刻,k的初始值为0。Ro(k)表示k时刻的电池欧姆内阻;Rc(k)表示k时刻的极化内阻;C(k)表示k时刻的极化电容,Uocv(k)表示k时刻的开路电压。
S03,所述递推估算模块30用于依据电压U、电流Cur并结合参数控制方程,利用遗忘递推最小二乘算法计算出θ1、θ2、θ3、θ4四个形式参数并输出。
具体的计算公式如下:
K(k)=p(k-1)ψ(k)/(α+ψ(k)'p(k-1)ψ(k));
p(k)=1/α[eye(4)-K(k)ψ(k)']p(k-1);
Upri(k)=ψ(k)'θ(k-1);
θ(k)=θ(k-1)+K(k)[U(k)-Upri(k)];
其中,k表示当前时刻,k-1表示上一时刻;ψ是由四个元素组成的向量,ψ(k)=[Cur(k);Cur(k)-Cur(k-1);U(k)-U(k-1);1];θ是由四个元素组成的向量,θ(k)=[θ1(k);θ2(k);θ3(k);θ4(k)];K、p、Upri是计算过程变量;U(k)、Cur(k)分别是所述采集模块10采集到的k时刻电压、电流值;α是遗忘因子,取值范围为0~1,在本实例中,α取值为0.995。
S04,所述参数计算模块40用于依据所述递推估算模块30输出的θ1、θ2、θ3、θ4四个形式参数计算电池欧姆内阻Ro(k)、极化内阻Rc(k)、极化电容C(k)及开路电压Uocv(k)四个电池参数。
请再一并参阅图5-8,分别为所述参数计算模块40依据所述递推估算模块30输出的θ1、θ2、θ3、θ4四个形式参数计算电池欧姆内阻Ro(k)、极化内阻Rc(k)、极化电容C(k)及开路电压Uocv(k)的在线辨识结果图。
具体的计算公式如下:
Ro(k)=-θ2(k)/θ3(k);
RC(k)=θ1(k)-Ro(k);
C(k)=-θ3(k)/RC(k);
Uocv(k)=θ4(k)。
本发明所提供的电池参数在线辨识方法及电池参数在线辨识系统200,使用了一种简单的电池模型以及参数控制方程,并结合遗忘递推最小二乘法算法,实现电池参数的在线辨识。具有高度实时性、低计算复杂度及较高辨识精度等优点。
本发明并不仅仅限于说明书和实施方式中所描述,因此对于熟悉领域的人员而言可容易地实现另外的优点和修改,故在不背离权利要求及等同范围所限定的一般概念的精神和范围的情况下,本发明并不限于特定的细节、代表性的设备和这里示出与描述的图示示例。