威布尔分布用于其他场景时的具体参数估计

回顾

威布尔累积分布为
F ( t ) = 1 − e x p [ − ( t η ) β ] F(t)=1-exp\left [ -\left ( \frac{t}{\eta } \right )^{\beta } \right ] F(t)=1exp[(ηt)β]
此分布可以用在很多领域(不仅仅是质量可靠性来反映时间长度与故障累积概率)
由前面我们了解到常用最小二乘法用于求得威布尔分布的参数,详情回顾见此
博客威布尔分布的参数估计提到,对于样本 ( x i , F ( x i ) ) (x_i,F(x_i)) (xi,F(xi))的获取,对于可靠性而言:

X是设备发生故障的时间点,这可以从故障数据中获得。假如取100次故障发生的时间点,则:

X = [ x 1 , x 2 , x 3 , … , x 100 ] X = [x_1,x_2,x_3,…,x_{100}] X=[x1,x2,x3,,x100]

其中 x 1 < x 2 < x 3 < … < x 100 x_1<x_2<x_3<…<x_{100} x1<x2<x3<<x100

这里故障发生的时间点是按照时间顺序排列的, F ( x i ) F(x_i) F(xi)是截止到时间xi,设备发生故障的概率

显然是没有问题的,因此在按照顺序求和的时候,作为累积概率 F ( x i ) F(x_i) F(xi)应当是跟着累积的个数均匀递增的,因此理所当然可以设为 F ( x i ) = i n + 1 F(x_i)=\frac{i}{n+1} F(xi)=n+1i
但当可靠性之外的适用场景时,自变量 x i x_i xi往往不是按照顺序排布的,因此在不进行排序之前,不能将 F ( x i ) F(x_i) F(xi)简单的设为 i n + 1 \frac{i}{n+1} n+1i

思路

应当保留 F ( x i ) = i n + 1 F(x_i)=\frac{i}{n+1} F(xi)=n+1i,因为此时 F ( x i ) F(x_i) F(xi)是最简单的。个人倾向于在参数估计和直线拟合的时候提前对 x i x_i xi进行排序,再按照前面所说的进行。
其中排序算法可以考虑使用快速排序算法,而且需要每次送来新的样本参加排序。
在初始化完成之后,每次排序只需要最新样本进行排序即可,然后重新迭代更新分布。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
估计威布尔分布的参数有多种方法,以下是几种常见的方法: 1. 最大似然估计(Maximum Likelihood Estimation, MLE):最大似然估计是一种常用的参数估计方法。它通过寻找能够最大化观测数据的似然函数的参数值来进行估计。对于威布尔分布,可以通过最大化威布尔分布的似然函数来估计形状参数和尺度参数。 2. 图形法(Graphical Method):图形法是一种简单直观的估计方法,它基于观察数据的累积分布函数(CDF)和经验分位数图。通过绘制经验分位数图并与威布尔分布的理论分位数进行比较,可以粗略地估计威布尔分布的参数。 3. 矩估计(Moment Estimation):矩估计是一种基于样本矩的参数估计方法。通过计算样本矩和理论矩之间的差异,可以得到威布尔分布参数估计值。对于威布尔分布,可以使用样本均值和样本方差来进行矩估计。 4. 最小二乘法(Least Squares Estimation):最小二乘法是一种常用的拟合方法,它通过最小化观测数据与拟合模型之间的平方差来进行参数估计。对于威布尔分布,可以将观测数据的对数转换为线性形式,然后使用最小二乘法来拟合直线,从而估计威布尔分布的参数。 值得注意的是,不同的参数估计方法可能会产生略微不同的结果。选择合适的估计方法取决于数据的特点、具体应用场景和估计结果的准确性要求。在实际应用中,可以尝试多种方法来进行参数估计,并进行对比和验证,以选择最合适的估计方法。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值