3. Weibull分布的统计特性:均值、方差、中位数、众数、失效率函数

Weibull分布之所以在可靠性工程与寿命数据分析中占据核心地位,很大程度上归功于其丰富的统计特性。通过调整形状参数 β 和尺度参数 η,Weibull分布能够灵活地模拟不同阶段的失效模式。本节将系统推导并解释其关键统计量:均值、方差、中位数、众数以及失效率函数。

3.1 概率密度函数与累积分布函数回顾

在深入统计特性之前,先明确两参数Weibull分布的标准数学表达。设随机变量 T 服从Weibull分布,其概率密度函数(PDF)和累积分布函数(CDF)分别为:

PDF:  f(t) = (β/η) * (t/η)^(β-1) * exp[-(t/η)^β],  t ≥ 0, β > 0, η > 0
CDF:  F(t) = 1 - exp[-(t/η)^β]

其中:

  • β(形状参数):决定分布曲线的形状,β < 1 时失效率递减(早期失效),β = 1 时退化为指数分布(恒定失效率),β > 1 时失效率递增(耗损失效)。
  • η(尺度参数):也称为特征寿命,表示63.2%的个体失效的时间点,即 F(η) = 1 - e⁻¹ ≈ 0.632。

3.2 均值(期望寿命)

Weibull分布的均值 E(T) 表示平均失效时间(MTTF),其计算公式涉及Gamma函数:

E(T) = η * Γ(1 + 1/β)

其中 Γ(·) 为Gamma函数,定义为 Γ(x) = ∫₀^∞ t^(x-1) e⁻ᵗ dt。Gamma函数在 β 取特定值时具有简洁形式:

形状参数 β Γ(1 + 1/β) 值 均值 E(T) 典型场景
0.5 Γ(3) = 2.000 2.000 η 早期失效(制造缺陷)
1.0 Γ(2) = 1.000 1.000 η 随机失效(指数分布)
2.0 Γ(1.5) ≈ 0.886 0.886 η 疲劳失效(瑞利分布)
3.5 Γ(1.286) ≈ 0.900 0.900 η 磨损失效
工程提示: 当 β 在 2.5~3.5 之间时,Weibull分布近似正态分布,此时均值接近 η 的 0.88~0.90 倍。若 β 趋近无穷大,分布退化为单点分布,均值趋近 η。

3.3 方差与标准差

方差衡量寿命数据的离散程度,其表达式为:

Var(T) = η² * [Γ(1 + 2/β) - Γ²(1 + 1/β)]
标准差 SD(T) = η * √[Γ(1 + 2/β) - Γ²(1 + 1/β)]

方差随 β 的变化呈现明显规律:当 β 较小时(β < 1),数据高度分散,方差较大;当 β 增大时,分布逐渐集中,方差减小。下表给出不同 β 下的标准化方差(Var/η²):

β Γ(1 + 2/β) Γ²(1 + 1/β) Var/η² 离散程度
0.5 Γ(5) = 24.000 4.000 20.000 极高
1.0 Γ(3) = 2.000 1.000 1.000 中等
2.0 Γ(2) = 1.000 0.785 0.215 较低
3.5 Γ(1.571) ≈ 0.890 0.810 0.080 低

3.4 中位数

中位数是寿命数据中50%个体失效的时间点,具有稳健性,不易受极端值影响。由CDF定义 F(t₅₀) = 0.5 可直接推导:

0.5 = 1 - exp[-(t₅₀/η)^β]
→ exp[-(t₅₀/η)^β] = 0.5
→ (t₅₀/η)^β = -ln(0.5) = ln(2)
→ t₅₀ = η * [ln(2)]^(1/β)

中位数与均值的比值 t₅₀ / E(T) 反映了分布的偏斜程度:

  • 当 β = 1 时,t₅₀ = η·ln(2) ≈ 0.693η,均值 = η,中位数小于均值,分布右偏。
  • 当 β = 2 时,t₅₀ = η·√(ln2) ≈ 0.833η,均值 ≈ 0.886η,两者接近。
  • 当 β = 3.5 时,t₅₀ ≈ η·(0.693)^(1/3.5) ≈ 0.901η,均值 ≈ 0.900η,几乎对称。

3.5 众数

众数是PDF取最大值时对应的时间点,表示最可能发生失效的时刻。通过对PDF求导并令导数为零,得到众数表达式:

当 β > 1 时:t_mode = η * [(β - 1)/β]^(1/β)
当 β ≤ 1 时:t_mode = 0(PDF在 t=0 处单调递减)

众数的存在性直接反映了失效模式:

  • β ≤ 1:失效率非增,失效概率在初始时刻最高,众数为0。适用于早期失效或随机失效。
  • β > 1:失效率递增,存在一个非零的众数,表明产品在经历一段“磨合期”后进入耗损阶段,失效概率先增后减。
关键对比: 均值、中位数、众数三者的关系可用于判断分布偏斜。对于Weibull分布,当 β > 1 时,通常有 众数 < 中位数 < 均值,呈现右偏特征;β 越大,偏斜程度越小。

3.6 失效率函数

失效率函数(Hazard Function)是可靠性工程中最具实际意义的统计量,定义为在时刻 t 尚未失效的产品,在下一瞬间发生失效的条件概率密度:

h(t) = f(t) / R(t) = f(t) / [1 - F(t)]

代入Weibull的PDF和可靠度函数 R(t) = exp[-(t/η)^β],得到简洁形式:

h(t) = (β/η) * (t/η)^(β-1)

失效率函数的三种基本形态直接对应产品的生命周期阶段:

β 范围 h(t) 趋势 生命周期阶段 工程实例
β < 1 递减(DFR) 早期失效期 电子元器件老化筛选、制造缺陷
β = 1 恒定(CFR) 偶然失效期 指数分布、随机外部冲击
β > 1 递增(IFR) 耗损失效期 机械磨损、材料疲劳、轴承失效
注意: 失效率函数 h(t) 与概率密度函数 f(t) 不同。f(t) 描述整体失效概率密度,而 h(t) 描述“幸存者”中的瞬时失效风险。在可靠性分析中,h(t) 比 f(t) 更能直观反映产品当前的健康状态。

3.7 统计特性关系图

下图展示了Weibull分布各统计量随形状参数 β 的变化趋势,以及失效率函数的三种基本形态:

Weibull分布统计特性与失效率函数形态 统计量随 β 变化趋势 均值 E(T) = η·Γ(1+1/β) β↑ → Γ(1+1/β)↓ → 均值趋近 0.886η~1.000η 中位数 t₅₀ = η·[ln(2)]^(1/β) β↑ → t₅₀ 从 0.693η 增至接近 η 众数 t_mode = η·[(β-1)/β]^(1/β) β≤1 时众数为0;β>1 时众数 > 0 失效率函数 h(t) 三种形态 β < 1:递减(DFR) 早期失效期:h(t) 随时间下降 β = 1:恒定(CFR) 偶然失效期:h(t) = 1/η 常数 β > 1:递增(IFR) 耗损失效期:h(t) 随时间上升 统计量关系总结(β > 1 时) 众数 < 中位数 中位数 < 均值 分布右偏 β 越大,三统计量越接近,分布越对称;β→∞ 时退化为单点分布

3.8 统计特性速查表

为便于工程应用,下表汇总了Weibull分布所有关键统计量的计算公式:

统计量 符号 计算公式 备注
均值 E(T) η · Γ(1 + 1/β) MTTF,需查Gamma函数表
方差 Var(T) η² · [Γ(1+2/β) - Γ²(1+1/β)] 离散程度度量
中位数 t₅₀ η · [ln(2)]^(1/β) 50%失效时间,稳健统计量
众数 t_mode η · [(β-1)/β]^(1/β), β>1 β≤1 时众数为0
失效率 h(t) (β/η) · (t/η)^(β-1) 三种形态对应生命周期
特征寿命 η F(η) = 0.632 63.2%失效时间
实用建议: 在实际工程中,建议优先使用中位数作为寿命数据的中心趋势度量,因为它对异常值不敏感。失效率函数 h(t) 的形态判断是选择维护策略的关键依据:β < 1 时宜做早期筛选,β > 1 时需规划预防性更换。