欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2017106588190
申请人: 山东科技大学
专利类型:发明专利
专利状态:已下证
专利领域: 计算;推算;计数
更新日期:2023-12-11
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.一种基于两阶段随机退化建模的剩余寿命预测方法,其特征在于:该方法包括离线建模以及在线参数更新与剩余寿命预测两部分,具体包括如下步骤:步骤1:离线建模过程,具体包括如下步骤:步骤1.1:收集n个退化设备历史退化数据,建立训练数据集X={X1,X2,...,Xn},其中表示第i个退化设备在时刻 处共有mi个退化数据;假设所有监测为等间隔监测,令Δt=ti,j-ti,j-1;

步骤1.2:定义两阶段退化模型如下所示:其中,μ1与μ2为两阶段Wiener过程模型的漂移系数,σ1与σ2为扩散参数,为了描述不同样本间的差异性,令 以及 μ1p, 和μ2p, 为正态分布的期望与方差,但对于每个单独样本,μ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:根据单阶段退化建模方法得到剩余寿命预测结果,如下: