CN103389472A - 一种基于nd-ar模型的锂离子电池循环寿命的预测方法 - Google Patents

一种基于nd-ar模型的锂离子电池循环寿命的预测方法 Download PDF

Info

Publication number
CN103389472A
CN103389472A CN2013103319989A CN201310331998A CN103389472A CN 103389472 A CN103389472 A CN 103389472A CN 2013103319989 A CN2013103319989 A CN 2013103319989A CN 201310331998 A CN201310331998 A CN 201310331998A CN 103389472 A CN103389472 A CN 103389472A
Authority
CN
China
Prior art keywords
partiald
capacity
centerdot
parameter
battery
Prior art date
Application number
CN2013103319989A
Other languages
English (en)
Other versions
CN103389472B (zh
Inventor
刘大同
周建宝
郭力萌
彭宇
彭喜元
Original Assignee
哈尔滨工业大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 哈尔滨工业大学 filed Critical 哈尔滨工业大学
Priority to CN201310331998.9A priority Critical patent/CN103389472B/zh
Publication of CN103389472A publication Critical patent/CN103389472A/zh
Application granted granted Critical
Publication of CN103389472B publication Critical patent/CN103389472B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

一种基于ND-AR模型的锂离子电池循环寿命的预测方法,涉及锂离子电池循环寿命的预测方法,它为解决AR模型没有包含容量退化数据的非线性特征,对非线性数据预测的失配问题。一、选取预测起始点,得到容量数据并进行建模,二、建立AR模型,并对AR模型预测容量;三、获取非线性退化因子的真实值;四、获得离线ND-AR模型:五、进行容量前期退化特征相似度分析:六、在线预测电池ND-AR模型参数加权估计:七、对AR模型预测结果的非线性修正:八、完成锂离子电池容量长期退化趋势的预测。它可广泛适用于对锂离子电池循环寿命进行预测。

Description

—种基于ND-AR模型的锂离子电池循环寿命的预测方法
技术领域
[0001] 本发明涉及锂离子电池循环寿命的预测方法,属于电池寿命预测领域。
背景技术
[0002] 现有电池循环寿命预测一般均采用线性AR模型,但是现有的线性AR模型无法直接对随时间呈现非线性退化特征的电池容量退化曲线进行准确预测。由于AR模型没有包含容量退化数据的非线性特征,也就是说对于加速退化趋势或未来将要出现的退化趋势变化没有预测能力导致个别样本的预测效果可能较差的现象,为解决AR线性模型对非线性数据预测的失配问题,提出本申请改进AR模型下的预测框架的研究。
发明内容
[0003] 本发明目的是为了解决AR模型没有包含容量退化数据的非线性特征,对非线性数据预测的失配问题,提供了一种基于ND-AR模型的锂离子电池循环寿命的预测方法。
[0004] 一种基于ND-AR模型的锂离子电池循环寿命的预测方法,实现该方法的步骤为:
[0005] 步骤一、通过锂离子电池测试平台对与待预测电池样本同型号的离线测试电池样本进行充放电循环实验,获取与待预测电池样本同型号的离线测试电池样本的离线测试容量数据Capacity, Capacity为非负实数,并从获得的离线测试容量数据Capacity中选取预测起始点T,并将预测起始点T之前的电池容量数据作为容量数据F进行建模,F为有理数集合;将预测起始点T之后的电池容量数据定义为Capreal ;
[0006] 所述预测起始点T的选 取是根据待测电池样本的容量数据F’占整体寿命的容量数据的百分比来确定的;;
[0007] 步骤二、根据步骤一所获得的容量数据F建立AR模型,并通过建立的AR模型获取容量数据F对应于电池容量数据Capreal的AR模型预测容量ARpredict ;
[0008] 步骤三、将步骤二所获得的容量数据F对应于电池容量数据Capreal的AR模型预测容量ARpredict与步骤一得到的电池容量数据Capreal比对,获取该AR模型预测容量退化特征与实际电池容量退化特征的比例误差,即非线性退化因子的真实值KT,Mal ;
[0009] 步骤四、将步骤三所得到的非线性退化因子的真实值KT@al以具体的函数表达式进行表征,引入非线性退化因子Kt的公式(10)获得离线ND-AR模型:
[0010]
Figure CN103389472AD00101
[0011] 其中参数&为包含电池退化特征信息的非线性退化因子;φ身=IX…,P)为
AR模型的自回归系数,Xt为当前时刻系统状态,即离线电池样本当前时刻的容量值,xt_i(i=l,2,...,p)为系统t-Ι时刻至t-p时刻的状态,即离线电池样本在相应时刻的容量值,at为噪声,服从均值为0,方差为W的正态分布,即at〜N(0,W),其中W为实数;
[0012] 通过采用指数型因子来描述电池容量的非线性退化信息,将公式(11)所示因子与公式(12)所示因子同时进行对比实验,比较不同表示形式的因子对预测效果的影响,获得一种包含更多退化信息的非线性退化因子:
[0013] [0014]
Figure CN103389472AD00111
[0015] 式中参数k表示的是预测步长,k的取值范围从I到n,其中η为离线测试容量数据Capacity的总长度,参数a、b、C、d代表待确定参数;通过上述两种形式,从不同的角度出发,由于非线性退化因子的真实值KT,Mal是一个在I附近微小变化的值,而公式(12)表示的因子数值在I附近变化,且为随着预测步长增加而呈现不同退化速率的因子形式;因此公式12所示非线性退化因子Kt更接近非线性退化因子的真实值KT,Mal ;
[0016] 建立所述离线ND-AR模型的具体方法如下:
[0017] 步骤四一、提取预测起始点T之后的电池容量数据Capreal ;所述电池容量数据Capreal是指真实容量退化特征信息;
[0018] 步骤四二、采用下式计算非线性退化因子Kt的真实值KT@al:
[001 9]
Figure CN103389472AD00112
[0020] 步骤四三、基于EKF算法对得到的非线性退化因子Kt的真实值KT,Mal的未知参数进行状态跟踪,同时进行步骤四三一和步骤四三二 ;以获取每一次放电循环所对应的非线性退化因子Kt的具体因子参数:
[0021] 步骤四三一:将公式(11)的因子参数a、b、c和d进行参数估计后代入式(11)中,获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;
[0022] 步骤四三二、将公式(12)的因子参数a和b,进行参数估计后,代入式(12)中获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;
[0023] 步骤五、进行离线电池样本与在线待预测电池容量前期退化特征相似度分析:
[0024] 首先对离线基于真实退化信息进行ND-AR模型建模的电池样本和在线寿命待预测电池样本的容量退化趋势进行关联性分析:采用灰色关联分析方法得到离线建模容量序列与在线待预测容量序列的前期历史容量数据变化趋势之间的关联度,关联度越大说明退化趋势越相近,非线性退化因子Kt的参数越接近,对应的加权权重越大;
[0025] 使用基于关联度的加权方法计算并获取待测电池样本的非线性退化因子Kt参数的估计值:
[0026] 所述灰色关联分析为:
[0027] 首先,确定分析数列:
[0028] 确定反映系统行为特征的参考数列和影响系统行为的比较数列;其中,反映系统行为特征的数据序列,称为参考数列;影响系统行为的因素组成的数据序列,称为比较数列;设参考数列为 Y= {Y (k) k=l, 2,...,η};比较数列为 Xi=IXi (k) k=l, 2,...,η},i=l, 2,...,m ;[0029] 为消除观测数据中包含的噪声,对数据样本进行多项式拟合,拟合次数根据数据特征来确定,并且从连续的拟合曲线中等间隔采集相同个数的数据点;获取参考数列Y为对在线待预测电池在线采集容量数据拟合并采样后的数据集,比较数列为离线测试电池样本经相同处理后获取的数据集;
[0030] 计算关联系数:
[0031] X0 (k)与XiGO的关联系数为:
[0032] g
Figure CN103389472AD00121
[0033] 式中y(k)为经确定分析数列处理后的待预测电池样本数据,Xi (k)为经确定分析数列处理后的第i个离线测试电池样本数据;P e(o,⑴),称为分辨系数;P越小,分辨力越大,
[0034] 将各个时刻关联系数集中为一个值,求其平均值,并作为比较数列与参考数列间关联程度的数量表示,第i个离线测试电池样本Xi (k)与在线待预测电池样本y(k)的关联度&公式如下:
[0035]
Figure CN103389472AD00122
[0036] 由此获得离线ND-AR建模电池样本容量数据与在线待预测电池容量数据退化趋势间的关联度用于后期非线性退化因子参数的估计;
[0037] 步骤六、在线预测电池ND-AR模型参数加权估计:
[0038] 选用两个电池样本模拟离线建模电池样本,两组拟合参数分别记为m1和m2,m代表了参数,对于式(11)的因子,m可以为a、b、c或者d,对于式(12)的因子,m可以为a或者b,1、2作为组别的区分而不是指数;通过获得的两组离线建模电池与在线待预测电池容量退化趋势的关联度A和r2,通过式(45)即可得到待预测电池ND-AR模型参数估计结果:
Figure CN103389472AD00123
[0040] 步骤七、对AR模型预测结果的非线性修正:
[0041] 通过步骤六获得的模型参数加权估计结果后,直接利用AR模型对在线待预测电池的容量退化数据进行预测,预测完成后按公式(10)进行容量预测结果的非线性校正,完成基于ND-AR模型的锂离子电池循环寿命的预测;
[0042] 步骤八、通过上述过程完成锂离子电池容量长期退化趋势的预测。
[0043] 本发明的优点:通过选取某类电池中的一部分电池样本在锂离子电池测试平台上进行充放电循环实验,获取其真实容量信息,并且利用AR模型进行建模获取相应AR模型预测容量,通过预测容量与真实容量比对获取AR模型预测容量退化特征与实际电池容量退化特征的比例误差,将此误差以某一种具体的函数表达式进行表征,且利用EKF算法获取表达式中的具体参数,也就是说构建一个包含此误差信息或者说包含电池非线性退化特征信息的非线性退化因子Kt,通过多组重复试验获取所有测试电池样本的相应退化因子表达式,完成基于真实退化信息的各自独立ND-AR模型构建。随后,在线对待预测电池样本进行容量信息采集,通过对比待预测电池数据样本与离线测试电池相对应数据样本的相似程度(比如这类电池的寿命大约也就是100个循环,在线已经采集了 50个循环的数据,大概也就是占整个寿命的50%,那么也对离线测试电池样本提取各自前50%的数据,分析这些数据与在线测得数据间的相似度),根据此相似程度估计在线待测电池可能的非线性退化特性,即估计其非线性退化因子的参数,并将此估计的非线性退化因子应用于AR直接预测输出的校正,即进行在线ND-AR模型推广预测,最终达到非线性退化趋势预测效果的提升。
具体实施方式
[0044] 具体实施方式一:本实施方式所述一种基于ND-AR模型的锂离子电池循环寿命的预测方法,实现该方法的步骤为:
[0045] 步骤一、通过锂离子电池测试平台对与待预测电池样本同型号的离线电池样本进行充放电循环实验,获取与待预测电池样本同型号的离线电池样本的离线测试容量数据Capacity, Capacity为非负实数,并从获得的离线测试容量数据Capacity中选取预测起始点T,并将预测起始点T之前的电池容量数据作为容量数据F进行建模,F为有理数;将预测起始点T之后的电池容量数据定义为Capreal ;
[0046] 所述预测起始点T的选取是根据在线待预测电池样本已采集的容量数据F’占整体寿命的容量数据的百分比来确定的;;
[0047] 步骤二、根据步骤一所获得的容量数据F建立离线电池样本的AR模型,并通过建立的AR模型获取容量数据Capacity中对应于电池容量数据Capreal的AR模型预测容量ARpredict ;
[0048] 步骤三、将步骤二所获得的容量数据Capacity中对应于电池容量数据Capreal的AR模型预测容量ARpredict与步骤一得到的电池容量数据Capreal比对,获取该AR模型预测容量退化特征与实际电池容量退化特征的比例误差,即非线性退化因子的真实值KT,Mal ;
[0049] 步骤四、将步骤三所得到的非线性退化因子的真实值KT@al以具体的函数表达式进行表征,引入非线性退化因子Kt的公式(10)获得离线ND-AR模型:
[0050]
Figure CN103389472AD00131
[0051] 其中参数&为包含电池退化特征信息的非线性退化因子;《(/ = 1,2,…为
AR模型的自回归系数,Xt为当前时刻系统状态,即离线电池样本当前时刻的容量值,xt_i(i=l,2,...,p)为系统t-Ι时刻至t-p时刻的状态,即离线电池样本在相应时刻的容量值,at为噪声,服从均值为0,方差为W的正态分布,即at〜N(0,W),其中W为实数;
[0052] 通过采用指数型因子来描述电池容量的非线性退化信息,将公式(11)所示因子与公式(12)所示因子同时进行对比实验,比较不同表示形式的因子对预测效果的影响,获得一种包含更多退化信息的非线性退化因子:
Figure CN103389472AD00132
[0055] 式中参数k表示的是预测步长,k的取值范围从I到n,其中η为离线测试容量数据Capacity的总长度,参数a、b、C、d代表待确定参数;通过上述两种形式,从不同的角度出发,由于非线性退化因子的真实值KTpal是一个在I附近微小变化的值,而公式(12)表示的因子数值在I附近变化,且为随着预测步长增加而呈现不同退化速率的因子形式;因此公式12所示非线性退化因子Kt更接近非线性退化因子的真实值KT,Mal ;
[0056] 建立所述离线ND-AR模型的具体方法如下:
[0057] 步骤四一、提取预测起始点T之后的电池容量数据Capreal ;所述电池容量数据Capreal是指真实容量退化特征信息;
[0058] 步骤四二、采用下式计算非线性退化因子Kt的真实值KT@al:
Figure CN103389472AD00141
[0060] 步骤四三、基于EKF算法对得到的非线性退化因子Kt的真实值KT,Mal的未知参数进行状态跟踪,同时进行步骤四三一和步骤四三二 ;以获取每一次放电循环所对应的非线性退化因子Kt的具体因子参数:
[0061] 步骤四三一:将公式(11)的因子参数a、b、c和d进行参数估计后代入式(11)中,获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;
[0062] 步骤四三二、将公式(12)的因子参数a和b,进行参数估计后,代入式(12)中获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;
[0063] 步骤五、进行离线电池样本与在线待预测电池容量前期退化特征相似度分析:
[0064] 首先对离线基于真实退化信息进行ND-AR模型建模的电池样本和在线寿命待预测电池样本的容量退化趋势进行关联性分析:采用灰色关联分析方法得到离线建模容量序列与在线待预测容量序列的前期历史容量数据变化趋势之间的关联度,关联度越大说明退化趋势越相近,非线性退化因子Kt的参数越接近,对应的加权权重越大;
[0065] 使用基于关联度的加权方法计算并获取待测电池样本的非线性退化因子Kt参数的估计值:
[0066] 所述灰色关联分析为:
[0067] 首先,确定分析数列:
[0068] 确定反映系统行为特征的参考数列和影响系统行为的比较数列;其中,反映系统行为特征的数据序列,称为参考数列;影响系统行为的因素组成的数据序列,称为比较数列;设参考数列为 Y= {Y (k) k=l, 2,...,η};比较数列为 Xi=IXi (k) k=l, 2,...,η},i=l, 2,...,m ;
[0069] 为消除观测数据中包含的噪声,对数据样本进行多项式拟合,拟合次数根据数据特征来确定,并且从连续的拟合曲线中等间隔采集相同个数的数据点;获取参考数列Y为对在线待预测电池在线采集容量数据拟合并采样后的数据集,比较数列为离线测试电池样本经相同处理后获取的数据集;
[0070] 计算关联系数:
[0071] X0 (k)与XiGO的关联系数为:
Figure CN103389472AD00151
[0073] 式中y(k)为经确定分析数列处理后的待预测电池样本数据,Xi (k)为经确定分析数列处理后的第i个离线测试电池样本数据;P e(0, m),称为分辨系数;P越小,分辨力越大,
[0074] 将各个时刻关联系数集中为一个值,求其平均值,并作为比较数列与参考数列间关联程度的数量表示,第i个离线测试电池样本Xi (k)与在线待预测电池样本y(k)的关联度&公式如下:
Figure CN103389472AD00152
[0076] 由此获得离线ND-AR建模电池样本容量数据与在线待预测电池容量数据退化趋势间的关联度用于后期非线性退化因子参数的估计;
[0077] 步骤六、在线预测电池ND-AR模型参数加权估计:
[0078] 选用两个电池样本模拟离线建模电池样本,两组拟合参数分别记为m1和m2,m代表了参数,对于式(11)的因子,m可以为a、b、c或者d,对于式(12)的因子,m可以为a或者b,1、2作为组别的区分而不是指数;通过获得的两组离线建模电池与在线待预测电池容量退化趋势的关联度A和r2,通过式(45)即可得到待预测电池ND-AR模型参数估计结果:
Figure CN103389472AD00153
[0080] 步骤七、对AR模型预测结果的非线性修正:
[0081] 通过步骤六获得的模型参数加权估计结果后,直接利用AR模型对在线待预测电池的容量退化数据进行预测,预测完成后按公式(10)进行容量预测结果的非线性校正,完成基于ND-AR模型的锂离子电池循环寿命的预测;
[0082] 步骤八、通过上述过程完成锂离子电池容量长期退化趋势的预测。
[0083] 具体实施方式二:本实施方式对实施方式一作进一步说明,所述步骤一中利用AR模型进行建模获取相应AR模型预测容量的具体方法为:首先,是AR模型的建模过程,模型阶次P的求取和自回归系数识O' = 1,2,..., P)的求取:
[0084] 步骤A:首先将容量数据F作为阶次判断的原始数据输入,对容量数据F进行零均值化和方差标准化处理;
[0085] 所述零均值化处理是指求取容量数据F的均值Fmean,从而得到零均值化的序列
Figure CN103389472AD00154
[0086] 所述方差标准化处理是指求取零均值化序列f的标准差σ f,得到标准化的建模数据Y=f/ O f ;Y为非负有理数集合;
[0087] 步骤B:判断标准化后的建模数据Y是否适合建立AR模型:
Figure CN103389472AD00155
[0089] 其中,R0为O步自协方差;L1为容量数据F数据集的长度;[0090] I 〜20 步自协方差
Figure CN103389472AD00161
[0091] 其中,R (k)为k步自协方差;
[0092]自相关系数:
Figure CN103389472AD00162
[0093] 根据计算结果,绘制自相关系数曲线,判断截尾特性即随着k的增加自相关系数是否逐步趋向于O,若表现截尾特征则适合建立时间序列分析模型MA模型,由于MA模型可由高阶AR模型进行近似,因此,若呈现截尾特性则表明同样适合AR模型建模;
[0094] 偏相关系数:求解Yule-Wallker方程,根据求解结果绘制偏相关系数曲线,判断截
[0095] 尾特性,若截尾则适合AR建模;
[0096] 步骤C:AIC计算:
[0097]通过自协方差计算得到:S= [R。,R(I),R(2),R(3) ] (4);
[0098] 其中,S为由O至3步自协方差构成的向量;
[0099]计算托普利茨矩阵 Toeplitz 矩阵:G=toeplitz (S) (5);
[0100] 其中,G为通过Matlab自带函数计算获得的向量S的托普利茨矩阵;
[0101]计算参数:W=G-1.[R(I),R(2),R(3),R(4)]T (6);
[0102] 其中W为计算过程中的中间向量
[0103]模型残差方差计算
Figure CN103389472AD00163
[0104] 其中,Y (t)是指t时刻对应的标准化后的建模数据即第t个建模数据;P为模型阶次;
[0105] 由此,AIC计算公式为:
Figure CN103389472AD00164
[0107] 其中,N为序列元素个数,σΐ为P阶预报误差方差,P为模型阶次;
[0108] 步骤D:将AIC最小值所对应的模型阶次P作为最优阶次;
[0109] 步骤E:对离线建模的各个电池样本分别进行上述步骤所述的建模数据的提取和AIC准则下的最佳模型阶次的求取,用于后续建模;
[0110] 步骤F:分别使用Burg法和Yule-Wallker法,利用相同的历史建模数据Y计算模型自回归系数狖(/=〗,2,...,尸),得到独立的系数求取结果氣i和ft,i ;其中,%和分别为
Burg法获取的自回归系数结果和Yule-Wallker法获取的自回归系数结果;
[0111] 步骤H:设置初始融合系数P1和P2 ;
[0112] 步骤1:随着预测步长的增加,动态调整融合系数=P1 = Prf (i), P2 = P2+f(i),其中i为预测步长;f (i)为融合系数动态调整因子,用于随预测步长的增加动态调整融合系数结果;
[0113] 步骤J:融合系数计算:m =P1-^,+,将此系数作为最终用以容量长期退化趋势预测的AR模型的系数;
[0114] 步骤K:利用上述步骤建立获得AR模型,如公式(9)所示:
[0115]
Figure CN103389472AD00171
[0116] 其中识(/ = 1,2,…,/0是待确定的自回归系数,P为模型的阶次;at, t=0, ±1,…为相互独立的白噪声序列,且服从均值为0,方差为CTe2的正态分布;
[0117] 获取相应AR模型预测容量,得到每一时刻的电池容量预测结果,所述每一时刻的电池容量预测结果组成容量退化长期预测输出数据集ARpredict。
[0118] 具体实施方式三:本实施方式对实施方式一作进一步说明,步骤四基于EKF算法对得到的非线性退化因子Kt的真实值KT,al的未知参数进行状态跟踪还包括下述步骤:
[0119] 步骤a、利用EKF算法对未知参数进行状态跟踪,首先要建立相应状态空间模型,寻求状态转移方程和观测方程,将待确定参数作为系统状态向量,构造的状态空间模型如公式(14)所示:
Figure CN103389472AD00172
[0121]其中方程伐:.:::::¾)为参数估计的状态转移方程,用于描述
上一时刻和下一时刻之间状态关系,a、b、c和d是式(11)中的因子参数,构成系统状态向
Figure CN103389472AD00173
为高斯白噪声,描述系统过程噪声,分别服从均值为0,方差为Qa、Qb、Qc和Qd的高斯分布;
[0122] 方程
Figure CN103389472AD00174
为系统观测方程,将估计获得的参数带入此方程获取一个非线性退化因子的估计值;
[0123] 对状态空间模型公式(14)进行线性化处理,由于状态转移方程是典型的线性方程,因此不需要进行线性化展开,只需要直接输入状态转移矩阵(15)即可:
Figure CN103389472AD00175
[0125] 观测方程为指数形式KT,k=f (ak,bk, ck, dk),对其进行Taylor展开并利用其一阶部分进行非线性方程的线性化近似如式(16)至(19):
[0126]
Figure CN103389472AD00176
[0127]
Figure CN103389472AD00181
[0130] 得到线性化后的观测矩阵Hk如式(20)所示:
[0131]
Figure CN103389472AD00182
[0132] 系统过程噪声及观测噪声均为线性叠加噪声,故有线性化过程噪声及观测噪声系数矩阵如式(21)、(22)所示:
[0133]
Figure CN103389472AD00183
[0134]
[0135] 由于各噪声互相独立,则有系统过程噪声协方差矩阵Q如式(23):
Figure CN103389472AD00184
[0137] 系统状态空间模型进行线性化处理后,进行状态的估计及更新过程,即对每一时刻模型的参数进行估计与更新:
[0138]所述参数估计是通过状态转移方程位
Figure CN103389472AD00185
激对模型参数进
行估计:
Figure CN103389472AD00186
[0141] 式中<、K、6和4分别代表k时刻状态的估计值
Figure CN103389472AD00191
和分
别代表k-ι时刻状态更新值,/T为k时刻系统状态协方差矩阵的估计值,AUk-1时刻系
统状态协方差矩阵的更新值,Fk为系统状态转移矩阵,Wk为线性化系统过程噪声系数矩阵,Qk为过程噪声方差;
[0142] 所述参数更新是通过状态估计后,得到当前时刻参数的先验估计值,将先验估计值带入系统观测方程
Figure CN103389472AD00192
,得到观测值的估计值;
[0143] 步骤b:将观测值的估计值与观测值真值进行比较得到测量余差,对状态估计值进行基于最小方差原则下的状态更新,得到最终的状态预测结果;具体更新步骤如下:
[0144]因子估计值:
Figure CN103389472AD00193
[0145]测量余差协方差:
Figure CN103389472AD00194
[0146]卡尔曼增益
Figure CN103389472AD00195
0147]状态更新
Figure CN103389472AD00196
[0149] 以上各式中,是基于估计参数所计算得到的第k个周期的非线性退化因子估计结果,KT,k是k个周期的真实非线性退化因子值KT,Mal(k),Sk是测量余差的协方差矩阵,Pt是状态协方差矩阵的估计值,Hk为观测矩阵,Vk为观测噪声系数矩阵,Rk为观测噪声方
差,Kk即为当前最优卡尔曼增益,apbpCk+和^更新后的状态值,P:为更新后的协方差矩阵;
[0150] 由此得到每一个时刻的参数估计值;
[0151] 得到每一时刻的参数估计值后,需要通过这些估计值综合出一组统一的参数a_s、b_s、c_s和d_s,以明确非线性退化因子最终的表达式;通过多元高斯分布概率密度计算获
得在当前观测值估计值€.,及相应测量余差协方差Sk的条件下得到测量真值的概率P,即
当前的参数估计值为参数真值的概率为P ;基于此概率进行加权平均,P值越大,说明相应的参数预测结果更为接近真实的参数值,因此应具有更高的权重,即其可信度更高;模型参数的确定按照式(31)进行;
Figure CN103389472AD00197
[0153] 其中,N即为参数估计序列的长度,P(k)为估计获取的参数即为真实参数的概率,m代表模型中的参数a、b、c或d。
[0154] 具体实施方式四:本实施方式对实施方式一作进一步说明,基于公式(12)的因子参数a和b,进行参数估计,获得建模电池样本最终的参数a和b之后,代入式(12)中获取相应的非线性退化因子Kt的表达式,并将此因子代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;对不同离线测试电池个体,重复上述步骤,获得每个电池样本各自基于真实退化信息的ND-AR模型;
[0155] 对应式(14)的位置,换为
Figure CN103389472AD00201
[0157] 对应式(15)的位置,换为
Figure CN103389472AD00202
[0159] 式(16)〜(19)的非线性模型线性化过程,替换为:
Figure CN103389472AD00203
[0162] 对应式(20)位置,换为
[0163]
Figure CN103389472AD00204
[0164] 对应式(21)位置,换为
Figure CN103389472AD00205
[0166] 对应式(23)位置,换为
[0167]
Figure CN103389472AD00206
[0168] 对应式(24)位置,换为
[0169]
Figure CN103389472AD00207
[0170] 对应式(26)位置,换为
Figure CN103389472AD00208
[0172] 对应式(29)位置,换为
[0173]
Figure CN103389472AD00209
[0174] 对应式(31)位置,换为
Figure CN103389472AD00211
[0176] 未提及的公式保持不变,通过相同的流程获取式(12)形式下的因子具体表达形式。

Claims (4)

1.一种基于ND-AR模型的锂离子电池循环寿命的预测方法,其特征在于,实现该方法的步骤为: 步骤一、通过锂离子电池测试平台对与在线待预测电池样本同型号的若干离线电池样本进行充放电循环实验,获取与在线待预测电池样本同型号的若干离线电池样本的离线测试容量数据Capacity, Capacity为非负实数,并从获得的离线测试容量数据Capacity中选取预测起始点T,并将预测起始点T之前的电池容量数据作为容量数据F进行建模,F为有理数集合;将预测起始点T之后的电池容量数据定义为Capreal ; 所述预测起始点T的选取是根据在线待预测电池样本已采集的容量数据F’占整体寿命的容量数据的百分比来确定的;; 步骤二、根据步骤一所获得的容量数据F建立离线电池样本的AR模型,并通过建立的AR模型获取容量数据Capacity对应于电池容量数据Capreal的AR模型预测容量ARpredict ; 步骤三、将步骤二所获得的容量数据Capacity对应于电池容量数据Capreal的AR模型预测容量ARpredict与步骤一得到的电池容量数据Capreal比对,获取该AR模型预测容量退化特征与实际电池容量退化特征的比例误差,即非线性退化因子的真实值KT,Mal ;步骤四、将步骤三所得到的非线性退化因子的真实值KT,Mal以具体的函数表达式进行表征,引入非线性退化因子Kt的公式(10)获得离线ND-AR模型:Xt = (P1XU +: + …+ φρχ,^ρ +α,).Κτ (10);其中参数Kt为包含电池退化特征信息的非线性退化因子;伊,(/=1,2为AR模型的自回归系数,Xt为当前时刻系统状态,即离线电池样本当前时刻的容量值,xt_i(i=l,2,...,p)为系统t-Ι时刻至t-p时刻的状态,即离线电池样本在相应时刻的容量值,at为噪声,服从均值为0,方差为W`的正态分布,即at〜N(0,W),其中W为实数;通过采用指数型因子来描述电池容量的非线性退化信息,将公式(11)所示因子与公式(12)所示因子同时进行对比实验,比较不同表示形式的因子对预测效果的影响,获得一种包含更多退化信息的非线性退化因子: KT=a.eb.k+c.ed.k (11); Kr =............................1............................................................(12); \-¥ ύ ■ ik -\-b) 式中参数k表示的是预测步长,k的取值范围从I到n,其中η为离线测试容量数据Capacity的总长度,参数a、b、C、d代表待确定参数;通过上述两种形式,从不同的角度出发,由于非线性退化因子的真实值KT,Mal是一个在I附近微小变化的值,而公式(12)表示的因子数值在I附近变化,且为随着预测步长增加而呈现不同退化速率的因子形式;因此公式12所示非线性退化因子Kt更接近非线性退化因子的真实值KT,Mal ; 建立所述离线ND-AR模型的具体方法如下: 步骤四一、提取预测起始点T之后的电池容量数据Capreal ;所述电池容量数据Capreal是指真实容量退化特征信息; 步骤四二、采用下式计算非线性退化因子Kt的真实值Klreal:
Figure CN103389472AC00031
步骤四三、基于EKF算法对得到的非线性退化因子Kt的真实值KT,Mal的未知参数进行状态跟踪,同时进行步骤四三一和步骤四三二 ;以获取每一次放电循环所对应的非线性退化因子Kt的具体因子参数: 步骤四三一:将公式(11)的因子参数a、b、c和d进行参数估计后代入式(11)中,获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模; 步骤四三二、将公式(12)的因子参数a和b,进行参数估计后,代入式(12)中获取相应的非线性退化因子Kt的表达式,并将此非线性退化因子Kt代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模; 步骤五、进行离线电池样本与在线待预测电池容量前期退化特征相似度分析: 首先对离线基于真实退化信息进行ND-AR模型建模的电池样本和在线寿命待预测电池样本的容量退化趋势进行关联性分析:采用灰色关联分析方法得到离线建模容量序列与在线待预测容量序列的前期历史容量数据变化趋势之间的关联度,关联度越大说明退化趋势越相近,非线性退 化因子Kt的参数越接近,对应的加权权重越大; 使用基于关联度的加权方法计算并获取待测电池样本的非线性退化因子Kt参数的估计值: 所述灰色关联分析为: 首先,确定分析数列: 确定反映系统行为特征的参考数列和影响系统行为的比较数列;其中,反映系统行为特征的数据序列,称为参考数列;影响系统行为的因素组成的数据序列,称为比较数列;设参考数列为 Y= {Y (k) k=l, 2,..., η};比较数列为 Xi=(XiQi) k=l, 2,..., n}, i=l, 2,…,m ; 为消除观测数据中包含的噪声,对数据样本进行多项式拟合,拟合次数根据数据特征来确定,并且从连续的拟合曲线中等间隔采集相同个数的数据点;获取参考数列Y为对在线待预测电池在线采集容量数据拟合并采样后的数据集,比较数列为离线测试电池样本经相同处理后获取的数据集; 计算关联系数: X0 (k)与XiGO的关联系数为: min min max 刚-调I gik) = ^^^-1——ί-— (43); Iy(k)-φ)\ +Pmax max Iy(k)-x,(k)\ i k 式中y(k)为经确定分析数列处理后的待预测电池样本数据,Xi(k)为经确定分析数列处理后的第i个离线测试电池样本数据;P e(o,⑴),称为分辨系数;P越小,分辨力越大, 将各个时刻关联系数集中为一个值,求其平均值,并作为比较数列与参考数列间关联程度的数量表示,第i个离线测试电池样本Xi(k)与在线待预测电池样本y(k)的关联度&公式如下:ri =-έ《(々) (k=l, 2,..., η) (44); η Ti 由此获得离线ND-AR建模电池样本容量数据与在线待预测电池容量数据退化趋势间的关联度用于后期非线性退化因子参数的估计; 步骤六、在线预测电池ND-AR模型参数加权估计: 选用两个电池样本模拟离线建模电池样本,两组拟合参数分别记为m1和m2,m代表了参数,对于式(11)的因子,m可以为a、b、c或者d,对于式(12)的因子,m可以为a或者b,1、2作为组别的区分而不是指数;通过获得的两组离线建模电池与在线待预测电池容量退化趋势的关联度A和r2,通过式(45)即可得到待预测电池ND-AR模型参数估计结果: Τλ \ Τ.Ί ^ i AC\ m= 1 m + - wr (45); r,+r2 #i+r2 步骤七、对AR模型预测结果的非线性修正: 通过步骤六获得的模型参数加权估计结果后,直接利用AR模型对在线待预测电池的容量退化数据进行预测,预测完成后按公式(10)进行容量预测结果的非线性校正,完成基于ND-AR模型的锂离子电池循环寿命的预测; 步骤八、通过上述过程完成锂离子电池容量长期退化趋势的预测。
2.根据权利要求1 所述一种基于ND-AR模型的锂离子电池循环寿命的预测方法,其特征在于所述步骤一中利用AR模型进行建模获取相应AR模型预测容量的具体方法为:首先,是AR模型的建模过程,模型阶次P的求取和自回归系数识(/ = 1,2,…,P)的求取: 步骤A:首先将容量数据F作为阶次判断的原始数据输入,对容量数据F进行零均值化和方差标准化处理; 所述零均值化处理是指求取容量数据F的均值Fmean,从而得到零均值化的序列f=F-Fmean ; 所述方差标准化处理是指求取零均值化序列f的标准差a,得到标准化的建模数据Y=f/O f ;γ为非负有理数集合; 步骤B:判断标准化后的建模数据Y是否适合建立AR模型: 其 F2CAο步自协方差CD 其中,RO为O步自协方差;L1为容量数据F数据集的长度; I 〜20 步自协方差:轉、=Σ —:.....(i = 1,2,...,20) (2) 其中,R (k)为k步自协方差; 自相关系数:_) = ^^ (3) 根据计算结果,绘制自相关系数曲线,判断截尾特性即随着k的增加自相关系数是否逐步趋向于0,若表现截尾特征则适合建立时间序列分析模型MA模型,由于MA模型可由高阶AR模型进行近似,因此,若呈现截尾特性则表明同样适合AR模型建模;偏相关系数:求解Yule-Wallker方程,根据求解结果绘制偏相关系数曲线,判断截尾 特性,若截尾则适合AR建模; 步骤C =AIC计算: 通过自协方差计算得到:S= [R0, R(I),R(2),R(3) ] (4); 其中,S为由O至3步自协方差构成的向量; 计算托普利茨矩阵Toeplitz矩阵:G=toeplitz (S) (5); 其中,G为通过Matlab自带函数计算获得的向量S的托普利茨矩阵; 计算参数:W=G-1* [R(1),R(2),R(3),R(4)]T (6); 其中W为计算过程中的中间向量 模型残差方差计算:
Figure CN103389472AC00051
其中,Y (t)是指t时刻对应的标准化后的建模数据即第t个建模数据;P为模型阶次; 由此,AIC计算公式为:
Figure CN103389472AC00052
其中,N为序列元素个数,《为P阶预报误差方差,P为模型阶次; 步骤D:将AIC最小值所对应的模型阶次P作为最优阶次; 步骤E:对离线建模的各个电池样本分别进行上述步骤所述的建模数据的提取和AIC准则下的最佳模型阶次的求取,用于后续建模; 步骤F:分别使用Burg法和Yule-Wallker法,利用相同的历史建模数据Y计算模型自回归系数:识(/ = 1,2-%ρ) »得到独立的系数求取结果Ψη和Ψχχ ;其中,I和分别为Burg法获取的自回归系数结果和Yule-Wallker法获取的自回归系数结果; 步骤H:设置初始融合系数P1和P2 ; 步骤1:随着预测步长的增加,动态调整融合系数=P1 = P「f(i),P2 = P2+f (i),其中i为预测步长;f(i)为融合系数动态调整因子,用于随预测步长的增加动态调整融合系数结果; 步骤J:融合系数计算:= H +H,将此系数作为最终用以容量长期退化趋势预测的AR模型的系数; 步骤K:利用上述步骤建立获得AR模型,如公式(9)所示:
Figure CN103389472AC00053
其中识(/ = 1,2,...,/?)是待确定的自回归系数,P为模型的阶次;at, t=0, ±1,…为相互独立的白噪声序列,且服从均值为0,方差为(Tj的正态分布; 获取相应AR模型预测容量,得到每一时刻的电池容量预测结果,所述每一时刻的电池容量预测结果组成容量退化长期预测输出数据集ARpredict。
3.根据权利要求1所述一种基于ND-AR模型的锂离子电池循环寿命的预测方法,其特征在于步骤四基于EKF算法对得到的非线性退化因子Kt的真实值Klreal的未知参数进行状态跟踪还包括下述步骤:步骤a、利用EKF算法对未知参数进行状态跟踪,首先要建立相应状态空间模型,寻求状态转移方程和观测方程,将待确定参数作为系统状态向量,构造的状态空间模型如公式(14)所示:
Figure CN103389472AC00061
其中方程I
Figure CN103389472AC00062
为参数估计的状态转移方程,用于描述上一 时刻和下一时刻之间状态关系,a、b、c和d是式(11)中的因子参数,构成系统状态向量[a;b ;c ;d],wa、wb、w。和wd为高斯白噪声,描述系统过程噪声,分别服从均值为O,方差为Qa、Qb、Qc和Qd的高斯分布; 方程=¾ # +¾ ;为系统观测方程,将估计获得的参数带入此方程获取一个非线性退化因子的估计值; 对状态空间模型公式(14)进行线性化处理,由于状态转移方程是典型的线性方程,因此不需要进行线性化展开,只需要直接输入状态转移矩阵(15)即可:
Figure CN103389472AC00063
观测方程为指数形式KT,k=f(ak,bk, ck, dk),对其进行Taylor展开并利用其一阶部分进行非线性方程的线性化近似如式(16)至(19):
Figure CN103389472AC00064
得到线性化后的观测矩阵Hk如式(20)所示:
Figure CN103389472AC00065
系统过程噪声及观测噪声均为线性叠加噪声,故有线性化过程噪声及观测噪声系数矩阵如式(21)、(22)所示:
Figure CN103389472AC00071
由于各噪声互相独立,则有系统过程噪声协方差矩阵Q如式(23):
Figure CN103389472AC00072
系统状态空间模型进行线性化处理后,进行状态的估计及更新过程,即对每一时刻模型的参数进行估计与更新: 所述参数估计是通过状态转移方程
Figure CN103389472AC00073
对模型参数进行估 计:
Figure CN103389472AC00074
式中%_、K Λ 和4分别代表k时刻状态的估计值,C1、4和C1分别代表k-ι时刻状态更新值,!:ζ.为k时刻系统状态协方差矩阵的估计值,P“为k-ι时刻系统状态协方差矩阵的更新值,Fk为系统状态转移矩阵,Wk为线性化系统过程噪声系数矩阵,Qk为过程噪声方差; 所述参数更新是通过状态估计后,得到当前时刻参数的先验估计值,将先验估计值带入系统观测方程=Ch +¾ -e4* ,得到观测值的估计值; 步骤b:将观测值的估计值与观测值真值进行比较得到测量余差,对状态估计值进行基于最小方差原则下的状态更新,得到最终的状态预测结果;具体更新步骤如下: γ Sr »« Ί »«.^« *義 闵子估i+侑-Λ =a -e 1 +c -e 1 /Λ,、 凶丁WTT1 且.ri t * (26);测量余差协方差:Sk =HkPk Hl +VkRkV- (27); 卡尔曼增益:Kt =iTifIC1 (28); 状态更新:]+&((,(29).1'....=(/-ΛΓ"/.」/Γ 嶋 以上各式中,是基于估计参数所计算得到的第k个周期的非线性退化因子估计结果,KT,k是k个周期的真实非线性退化因子值KT,Mal(k),Sk是测量余差的协方差矩阵,P;是状态协方差矩阵的估计值,Hk为观测矩阵,Vk为观测噪声系数矩阵,Rk为观测噪声方差,Kk即为当前最优卡尔曼增益,%+、h\ %+和(<4+更新后的状态值,if为更新后的协方差矩阵; 由此得到每一个时刻的参数估计值; 得到每一时刻的参数估计值后,需要通过这些估计值综合出一组统一的参数a_s、b_s、c_s和d_s,以明确非线性退化因子最终的表达式;通过多元高斯分布概率密度计算获得在当前观测值估计值及相应测量余差协方差Sk的条件下得到测量真值的概率P,即当前的参数估计值为参数真值的概率为P ;基于此概率进行加权平均,P值越大,说明相应的参数预测结果更为接近真实的参数值,因此应具有更高的权重,即其可信度更高;模型参数的确定按照式(31)进行; 尽,,、P(k) m s=y m(k)...—S1 j^p(k) m=a,b,c或d (31); 其中,N即为参 数估计序列的长度,P(k)为估计获取的参数即为真实参数的概率,m代表模型中的参数a、b、c或d。
4.根据权利要求1所述一种基于ND-AR模型的锂离子电池循环寿命的预测方法,其特征在于,基于公式(12)的因子参数a和b,进行参数估计,获得建模电池样本最终的参数a和b之后,代入式(12)中获取相应的非线性退化因子Kt的表达式,并将此因子代入式(10),完成该电池样本基于真实退化信息的ND-AR模型建模;对不同离线测试电池个体,重复上述步骤,获得每个电池样本各自基于真实退化信息的ND-AR模型; 对应式(14)的位置,换为 Cft = Cik^ + Wa Wa ~ N(0, Oa) « hk = hk_x + wb wh - N(0,0^) (32); KTA=ak-eb-k+ck-ed-k 对应式(15)的位置,换为 Γ ί οΊ ,、Fk= (33); —0 K 式(16)〜(19)的非线性模型线性化过程,替换为:
Figure CN103389472AC00091
对应式(20)位置,换为
Figure CN103389472AC00092
对应式(21)位置,换为
Figure CN103389472AC00093
对应式(23)位置,换为
Figure CN103389472AC00094
对应式(24)位置,换为
Figure CN103389472AC00095
对应式(26)位置,换为
Figure CN103389472AC00096
对应式(29)位置,换为
Figure CN103389472AC00097
对应式(31)位置,换为
Figure CN103389472AC00098
未提及的公式保持不变,通过相同的流程获取式(12)形式下的因子具体表达形式。
CN201310331998.9A 2013-08-01 2013-08-01 一种基于nd-ar模型的锂离子电池循环寿命的预测方法 Active CN103389472B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310331998.9A CN103389472B (zh) 2013-08-01 2013-08-01 一种基于nd-ar模型的锂离子电池循环寿命的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310331998.9A CN103389472B (zh) 2013-08-01 2013-08-01 一种基于nd-ar模型的锂离子电池循环寿命的预测方法

Publications (2)

Publication Number Publication Date
CN103389472A true CN103389472A (zh) 2013-11-13
CN103389472B CN103389472B (zh) 2015-09-16

Family

ID=49533799

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310331998.9A Active CN103389472B (zh) 2013-08-01 2013-08-01 一种基于nd-ar模型的锂离子电池循环寿命的预测方法

Country Status (1)

Country Link
CN (1) CN103389472B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926536A (zh) * 2014-03-07 2014-07-16 南京航空航天大学 基于dst和bmc技术的锂离子电池剩余使用寿命预测方法
CN104090238A (zh) * 2014-06-19 2014-10-08 哈尔滨工业大学 基于ar集成学习模型的锂离子电池剩余寿命预测方法
CN105068009A (zh) * 2015-07-14 2015-11-18 盐城工学院 电池循环寿命预测方法
CN106226699A (zh) * 2016-07-11 2016-12-14 北京航空航天大学 一种基于时变权重最优匹配相似性的锂离子电池寿命预测方法
CN106965685A (zh) * 2015-10-20 2017-07-21 福特全球技术公司 用于指示电池老化的系统和方法
CN107015157A (zh) * 2017-04-01 2017-08-04 湖南银杏数据科技有限公司 基于恒流等压升片段的锂电池剩余循环寿命在线快速测试法
CN107192954A (zh) * 2017-05-02 2017-09-22 杭州科工电子科技有限公司 一种锂离子电池性能在线诊断方法
CN107783050A (zh) * 2016-08-31 2018-03-09 北京海博思创科技有限公司 电池寿命预测方法及装置
CN112257348A (zh) * 2020-12-22 2021-01-22 北京航空航天大学 一种锂电池长期退化趋势预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030034779A1 (en) * 2001-03-21 2003-02-20 Carsten Juncker Battery life estimation
CN102213747A (zh) * 2010-04-09 2011-10-12 深圳市比克电池有限公司 一种预测锂离子电池循环寿命的方法
JP2012122817A (ja) * 2010-12-07 2012-06-28 Gs Yuasa Corp 非水電解質二次電池の可逆容量推定方法、寿命予測方法、可逆容量推定装置、寿命予測装置及び蓄電システム
CN102778653A (zh) * 2012-06-20 2012-11-14 哈尔滨工业大学 基于ar模型和rpf算法的数据驱动的锂离子电池循环寿命预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030034779A1 (en) * 2001-03-21 2003-02-20 Carsten Juncker Battery life estimation
CN102213747A (zh) * 2010-04-09 2011-10-12 深圳市比克电池有限公司 一种预测锂离子电池循环寿命的方法
JP2012122817A (ja) * 2010-12-07 2012-06-28 Gs Yuasa Corp 非水電解質二次電池の可逆容量推定方法、寿命予測方法、可逆容量推定装置、寿命予測装置及び蓄電システム
CN102778653A (zh) * 2012-06-20 2012-11-14 哈尔滨工业大学 基于ar模型和rpf算法的数据驱动的锂离子电池循环寿命预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JINGLIANG ZHANG ET AL.: "A review on prognostics and health monitoring of Li-ion battery", 《JOURNAL OF POWER SOURCES》 *
彭宇等: "故障预测与健康管理技术综述", 《电子测量与仪器学报》 *
许参等: "一种锂离子蓄电池寿命的预测模型", 《应用科学学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926536A (zh) * 2014-03-07 2014-07-16 南京航空航天大学 基于dst和bmc技术的锂离子电池剩余使用寿命预测方法
CN104090238A (zh) * 2014-06-19 2014-10-08 哈尔滨工业大学 基于ar集成学习模型的锂离子电池剩余寿命预测方法
CN104090238B (zh) * 2014-06-19 2016-09-28 哈尔滨工业大学 基于ar集成学习模型的锂离子电池剩余寿命预测方法
CN105068009A (zh) * 2015-07-14 2015-11-18 盐城工学院 电池循环寿命预测方法
CN105068009B (zh) * 2015-07-14 2017-10-20 盐城工学院 电池循环寿命预测方法
CN106965685A (zh) * 2015-10-20 2017-07-21 福特全球技术公司 用于指示电池老化的系统和方法
CN106226699A (zh) * 2016-07-11 2016-12-14 北京航空航天大学 一种基于时变权重最优匹配相似性的锂离子电池寿命预测方法
CN106226699B (zh) * 2016-07-11 2019-01-22 北京航空航天大学 一种基于时变权重最优匹配相似性的锂离子电池寿命预测方法
CN107783050A (zh) * 2016-08-31 2018-03-09 北京海博思创科技有限公司 电池寿命预测方法及装置
CN107015157A (zh) * 2017-04-01 2017-08-04 湖南银杏数据科技有限公司 基于恒流等压升片段的锂电池剩余循环寿命在线快速测试法
CN107015157B (zh) * 2017-04-01 2019-04-16 湖南银杏数据科技有限公司 基于恒流等压升片段的锂电池剩余循环寿命在线快速测试法
CN107192954A (zh) * 2017-05-02 2017-09-22 杭州科工电子科技有限公司 一种锂离子电池性能在线诊断方法
CN112257348A (zh) * 2020-12-22 2021-01-22 北京航空航天大学 一种锂电池长期退化趋势预测方法

Also Published As

Publication number Publication date
CN103389472B (zh) 2015-09-16

Similar Documents

Publication Publication Date Title
WO2018161486A1 (zh) 一种基于动态参数的动力电池soc估算方法及系统
Shao et al. Adaptive soft sensor for quality prediction of chemical processes based on selective ensemble of local partial least squares models
Hung et al. An artificial neural network model for rainfall forecasting in Bangkok, Thailand
CN103730006B (zh) 一种短时交通流量的组合预测方法
JP6312630B2 (ja) 時系列データ内の異常を検出する方法
Morando et al. Fuel cells prognostics using echo state network
CN103336248B (zh) 基于电池退化状态模型的锂离子电池循环寿命预测方法
CN101899563B (zh) 基于pca模型的连续退火机组炉内温度、张力监测及故障追溯方法
CN104166787B (zh) 一种基于多阶段信息融合的航空发动机剩余寿命预测方法
Gholami et al. Uncertainty analysis of intelligent model of hybrid genetic algorithm and particle swarm optimization with ANFIS to predict threshold bank profile shape based on digital laser approach sensing
Saha et al. Comparison of prognostic algorithms for estimating remaining useful life of batteries
CN104778837A (zh) 一种道路交通运行态势多时间尺度预测方法
CN103941196A (zh) 锂离子电池荷电状态估计方法
CN106405433A (zh) 一种基于扩展卡尔曼粒子滤波的soc估计方法及系统
CN104408913B (zh) 一种考虑时空相关性的交通流三参数实时预测方法
CN102608542A (zh) 动力电池荷电状态估计方法
Zhao et al. On-line least squares support vector machine algorithm in gas prediction
Singh et al. Estimation of traffic densities for multilane roadways using a markov model approach
CN102789545B (zh) 基于退化模型匹配的涡轮发动机剩余寿命的预测方法
Su et al. Time-varying identification model for dam behavior considering structural reinforcement
Muralidharan et al. Imputation of ramp flow data for freeway traffic simulation
CN103268519B (zh) 基于改进Lyapunov指数的电力系统短期负荷预测方法及装置
WO2013090026A2 (en) Automated monitoring for changes in energy consumption patterns
CN103472398B (zh) 基于扩展卡尔曼粒子滤波算法的动力电池soc估计方法
CN103777627A (zh) 一种基于少量批次的批次过程在线监测方法

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model