双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

根据故障数据求威布尔分布参数步骤

根据众多文献研究,设备故障数据大多服从指数分布、威布尔分布或对数正态分布。其中威布尔分布较为常用,本文主要介绍威布尔分布的简单应用,并根据某设备故障数据为例简单介绍威布尔分布参数的计算方法。
参考文献11参考文献22

1.双参数威布尔分布的概率密度函数、概率分布函数、可靠度函数:

双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

2.预处理故障数据

将故障数据按时间排序,假设设备分别在以下时间点失效:11.1 、23.5、58.4……则将此时间从小到大排序。并计算中位秩。中位秩的计算公式如下:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
其中i表示第几个数据,n表示一共多少数据,例如下图所示例子中共43个故障时间点,即n=43,下图所示例子中F(t1)=(1-0.3)/(43+0.4)=0.016129。示例:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
以ti为横坐标,中位秩为纵坐标,画图,从威布尔概率图上可以大致看出,绝大多数数据点分布在直线的周边,
可以认为该样本服从威布尔分布。示例如下:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

3.参数估计的简单原理

上文已经介绍过威布尔分布的概率分布函数为:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
则根据此公式可进行如下推导转化:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
因此,若看成线性函数y=Bx+A的话,A和B分别等于:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

因此可通过回归直线的,求出参数 A 与 B,从而求出两个参数。

3.参数估计的具体方法

将数据进行处理,根据上述推导, xi 和 yi 的公式如下,其中 ti 表示故障数据中的时间,F(ti)表示中位秩:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
A和B的估计求法公式为:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
其中 lxx 、lyy 和 lxy 公式为:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

根据方法,示例求解的各项数据结果如下:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
计算可得A和B结果为:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
根据A和B 求出结果为:
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)

4.相关性检验和模型检验

双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)
双参数威布尔分布(Weibull分布)参数计算(根据故障数据,求双参数威布尔分布的参数详细步骤)


  1. 吕言. 基于多目标优化的核电站系统维修决策研究[D].清华大学,2016. ↩︎

  2. 徐慧. 动态不确定环境下流程制造企业设备维修计划研究[D].上海交通大学,2009. ↩︎