本发明涉及电池领域,特别涉及一种带遗忘因子的最小二乘法锂电池模型参数辨识方法。
背景技术:
锂离子电池充放电过程的内部化学反应较为复杂,因此此过程是时变且非线性的,因而要通过理论分析获得模型的参数是很困难的。虽然目前有应用离线指数拟合方法辨识出了模型参数,但是由于电池系统的时变性,随着电池soc、外界温度、循环次数等因素的改变,其模型参数也会发生较大的变化,因此为了提高soc的估计精度,增强系统的适应能力,需要对电池模型参数在线辨识并做出实时修正。参数估计是在模型结构已知的情况下,通过采集的数据来确定模型参数的一种数学方法。通过构建锂离子电池的模型结构,结合锂离子电池ocv-soc关系曲线,目前最常用的方法是最小二乘法辨识理论,利用该理论可以对电池的二阶rc等效电路模型参数进行辨识。
递推最小二乘法是一种容易理解和掌握的辨识方法,实现起来也比较简单,大部分情况下,递推最小二乘法可以给出具有准确统计特性的参数辨识结果。只是在参数辨识中,递推最小二乘法是具有无限记忆长度的算法,对于电池系统,最小二乘法在递推运算过程中旧数据越来越多,会导致递推结果不能良好的反应新数据的特性。
技术实现要素:
本发明的目的在于提供一种带遗忘因子的最小二乘法锂电池模型参数辨识方法,大部分情况下,递推最小二乘法可以给出具有准确统计特性的参数辨识结果,只是在参数辨识中,递推最小二乘法是具有无限记忆长度的算法,对于电池系统,最小二乘法在递推运算过程中旧数据越来越多,会导致递推结果不能良好的反应新数据的特性,而引入遗忘因子能够有效的克服了“数据饱和”现象。
为实现上述目的,本发明提供的技术方案为:一种带遗忘因子的最小二乘法锂电池模型参数辨识方法,包括以下步骤:
步骤一,建立二阶rc等效电路模型,该等效电路模型的表达式为
通过对该表达式进行离散化后得到等效电路模型的拉普拉斯方程为
步骤二,对步骤一的等效电路模型拉普拉斯方程进行双线性变换,得出等效电路模型系统输入与系统输出的差分方程为
其中i(k)为系统输入,y(k)为系统输出,
θ=[a1a2a3a4a5],a2、a3、a4、a5为相应的常数系数;
步骤三,建立带遗忘因子λ的递推最小二乘法
其中
步骤四,建立系数方程
步骤五,采集电池的参数,以采样周期t来采集电池的端电压v(k)、端电流i(k),端电压v(k-1)、端电流i(k-1)、电池荷电状态soc(k-1)和端电压v(k-2)、端电流i(k-2);计算voc(k)-v(k),以获得辨识过程中的输入φ(k)、系统输出y(k)。
步骤六,对θ(0)、p(0)和遗忘因子λ进行初始化,根据步骤五采集的参数以及步骤三的递推最小二乘法求出步骤二的θ值,将θ值代入该系数方程中则得到系数方程中左边的系数,即辨识出等效电路模型的参数。
优选的,所述步骤四的采样周期t=1s,soc(0)=90%。
优选的,所述步骤六中的θ(0)为任意值,p(0)=αi,α=5000,λ=0.96。
本发明有益效果在于利用ocv-soc函数关系及带遗忘因子最小二乘法对二阶rc等效电路模型进行动态参数辨识,在常用的递推最小二乘法中引入遗忘因子,缓解在运算递推过程中旧数据叠加的问题,克服了“数据饱和”现象。
附图说明
图1是锂电池的二阶rc模型。
图2是不同倍率恒流间歇放电ocv-soc曲线。
具体实施方式
参照附图1至附图2介绍本发明的具体实施方式。
一种带遗忘因子的最小二乘法锂电池模型参数辨识方法,包括以下步骤:
步骤一,建立二阶rc等效电路模型,该等效电路模型的表达式为
通过对该表达式进行离散化后得到等效电路模型的拉普拉斯方程为对(1)式离散化,解得状态方程为:
其中:
由式(2)(3)可得电池模型拉普拉斯方程如下式。
步骤二,对步骤一的等效电路模型拉普拉斯方程(6)进行双线性变换进行离散化,令
得出等效电路模型系统输入与系统输出的差分方程为
其中i(k)为系统输入,y(k)为系统输出,
θ=[a1a2a3a4a5],a2、a3、a4、a5为相应的常数系数;
步骤三,建立带遗忘因子λ的递推最小二乘法,
设k时刻传感器采样误差为e(k),则:
将
取泛函数j(θ):
因为最小二乘法原理是使j(θ)取最小值,所以求j(θ)极值,令:
可得:
对上述(9)-(13)的过程进行递推运算,即得到式子(14),
递推最小二乘法是具有无限记忆长度的算法,对于电池系统,最小二乘法在递推运算过程中旧数据越来越多会导致递推结果不能良好的反应新数据的特性,为避免上述情况,引入遗忘因子λ,0<λ<1,即:
p-1(k+1)=λp-1(k)+φ(k+1)φt(k+1)(16);
所以,即使(n+1)很大,p(n+1)也不趋于0,有效的克服了“数据饱和”现象。带遗忘因子最小二乘算法的步骤为:
当λ=1时,为普通最小二乘法,λ越小跟踪能力越强,但波动也越大,一般取0.95<λ<1。
步骤四,建立系数方程,根据系统的采样测量信息,采样周期t=1s,soc(0)=90%。估计出等效电池模型的结构及未知参数,令:
由(6)(19)的系数对应相等可得:
其中t为采样周期,此式子右边的系数可以通过递推最小二乘法求出,左边即是模型的未知参数。
步骤五,采集电池的参数,以采样周期t来采集电池的端电压v(k)、端电流i(k),端电压v(k-1)、端电流i(k-1)、电池荷电状态soc(k-1)和端电压v(k-2)、端电流i(k-2);计算voc(k)-v(k),以获得辨识过程中的输入φ(k)、系统输出y(k)。
步骤六,对θ(0)、p(0)和遗忘因子λ进行初始化,θ(0)为任意值,p(0)=αi,α=5000,λ=0.96。根据步骤五采集的参数以及步骤三的递推最小二乘法求出步骤二的θ值,将θ值代入该系数方程中则得到系数方程中左边的系数,求得r,rs,rp,cp,cs,即辨识出等效电路模型的参数,进而可以按照自定义工况对电池模型参数进行动态、实时更新。
本发明以如图1的二阶rc等效电路模型为例,同时应用如图2的ocv-soc曲线,在常用的递推最小二乘法中引入遗忘因子,缓解在运算递推过程中旧数据叠加的问题,克服了“数据饱和”现象。
图2是对三元锂离子电池的充放电实验是在25摄氏度恒温条件下进行的,分别标定0.2c、0.3c、0.4c、0.5c、0.6c、0.75c、1c恒流间歇放电条件下的ocv-soc曲线。
下表为ocv-soc标定实验的设置和操作步骤。
每组标定步骤如下:
①采用先恒流(0.2c)后恒压(截止电压4.25v)的方式对电池进行充电;
②对电池进行恒流、恒容量(260mah)放电;
③放电结束,静置一小时以待极化效应消失;
④重复步骤②③,至电池放电结束。
以上所述并非对本发明的技术范围作任何限制,凡依据本发明技术实质,对以上的实施例所作的任何修改、等同变化与修饰,均仍属于本发明的技术方案的范围内。