[发明专利]一种基于两阶段随机退化建模的剩余寿命预测方法有效
申请号: | 201710658819.0 | 申请日: | 2017-08-04 |
公开(公告)号: | CN107480440B | 公开(公告)日: | 2020-01-21 |
发明(设计)人: | 周东华;张峻峰;何潇;张建勋;张海峰;卢晓 | 申请(专利权)人: | 山东科技大学 |
主分类号: | G06F30/20 | 分类号: | G06F30/20;G06F119/04 |
代理公司: | 37252 青岛智地领创专利代理有限公司 | 代理人: | 种艳丽 |
地址: | 266590 山东省青*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 阶段 随机 退化 建模 剩余 寿命 预测 方法 | ||
1.一种基于两阶段随机退化建模的剩余寿命预测方法,其特征在于:该方法包括离线建模以及在线参数更新与剩余寿命预测两部分,具体包括如下步骤:
步骤1:离线建模过程,具体包括如下步骤:
步骤1.1:收集n个电池历史退化数据,该退化数据指的是电池的电容量,建立训练数据集X={X1,X2,...,Xn},其中表示第i个电池在时刻处共有mi个退化数据;假设所有监测为等间隔监测,令Δt=ti,j-ti,j-1;
步骤1.2:定义两阶段退化模型如下所示:
其中,μ1与μ2为两阶段Wiener过程模型的漂移系数,σ1与σ2为扩散参数,为了描述不同样本间的差异性,令以及和为正态分布的期望与方差,但对于每个单独样本,μ1与μ2为固定常数;
步骤1.3:对每组退化数据X={X1,X2,...,Xn}分别构建如公式(2)所示的对数似然函数,并通过极大似然估计得到每一组退化数据的变点τi的估计值,如公式(3)所示:
其中,μi,σi表示第i组退化模型的漂移与扩散系数,Xi∈{X1,X2,...,Xn},xi,j表示第i组退化中第j个退化值,
其中,表示τi的估计值,并且通过统计分析得到变点的概率分布函数p(τ);
步骤1.4:基于EM算法离散估计参数,得到第k+1步参数迭代更新与第k步估计值关系如下:
其中,分别表示μ1p,σ1p,μ2p,σ2p,σ1,σ2在第k+1次迭代的估计值,
其中
其中,表示对中括号内取期望;
步骤1.5:对公式(4)与公式(5)进行不断迭代直至所有参数收敛,该收敛时的参数估计值即为离线得到的模型参数估计值,并将该估计值作为之后在线更新的先验信息;
步骤2:在线模型更新与剩余寿命预测,具体包括如下步骤:
步骤2.1:收集电池在线退化数据,假设当前时刻为tκ,其相应的退化数据为X0:κ={x0,x1,...,xκ},一共κ+1个数据;
步骤2.2:基于极大似然估计的变点检测,方法如下:
其中,xi∈{x0,x1,...,xκ},为变点时间τ的估计值;
若估计得到的那么说明到tκ时刻为止变点尚未出现,则执行步骤2.3;反之,则说明变点出现了且变点为则执行步骤2.5;
步骤2.3:基于贝叶斯理论利用先验信息更新第一阶段模型参数以及变点概率分布,如下:
其中,Pr(·)表示括号中事件的发生概率,μ1p,0表示μ1先验分布的期望,σ1p,0表示μ1先验分布的标准差;
步骤2.4:随机变点下基于两阶段退化模型的剩余寿命在线预测结果如下:
其中,表示在tκ剩余寿命为lκ的概率密度函数,p(τ)表示变点的概率分布函数;
其中,
μa4=μ2p(lκ-τ+tκ),μb4=ξ-xκ-μ1p(τ-tκ),
步骤2.5:基于贝叶斯理论利用先验信息更新第二阶段模型参数,如下:
其中,μ2p,0表示μ2先验分布的期望,σ2p,0表示μ2先验分布的标准差;
步骤2.6:根据单阶段退化建模方法得到剩余寿命预测结果,如下:
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于山东科技大学,未经山东科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201710658819.0/1.html,转载请声明来源钻瓜专利网。