2. 统计学基础回顾

可靠性试验的样本量设计与置信度计算,本质上是基于概率统计的推断过程。在进入具体公式之前,必须牢固掌握几个核心统计学概念:概率分布、均值与方差、以及置信区间。本章将围绕可靠性工程中最常用的三种分布展开,并建立从分布参数到置信推断的桥梁。

2.1 三种核心概率分布

在可靠性工程中,产品的寿命数据通常服从以下三种分布之一。选择正确的分布是后续计算样本量和置信区间的前提。

分布类型 适用场景 关键参数 失效率特征
指数分布 电子元器件、随机失效期(浴盆曲线底部) 失效率 λ(或平均寿命 MTBF = 1/λ) 恒定失效率
威布尔分布 机械零件、磨损/疲劳失效、早期失效 形状参数 β、尺度参数 η β<1 递减,β=1 恒定,β>1 递增
正态分布 耗损失效、产品性能退化、测量误差 均值 μ、标准差 σ 随时间急剧上升

2.1.1 指数分布

指数分布是可靠性工程中最基础的分布,其概率密度函数(PDF)为:

f(t) = λ e^{-λ t},  t ≥ 0

其中 λ 为失效率,MTBF(平均无故障时间)= 1/λ。指数分布具有“无记忆性”,即产品在 t 时刻正常工作的条件下,剩余寿命的分布与 t 无关。这一特性使得指数分布下的样本量计算相对简洁。

2.1.2 威布尔分布

威布尔分布通过形状参数 β 灵活描述不同失效模式,其 PDF 为:

f(t) = (β/η) (t/η)^{β-1} e^{-(t/η)^β},  t ≥ 0
  • β < 1:失效率递减,对应早期失效(磨合期)。
  • β = 1:退化为指数分布,对应随机失效期。
  • β > 1:失效率递增,对应磨损或疲劳失效。

威布尔分布是可靠性试验中最灵活的模型,样本量设计通常需要基于 β 的估计值。

2.1.3 正态分布

正态分布在耗损失效阶段常见,其 PDF 为:

f(t) = (1/(σ√(2π))) e^{-(t-μ)^2/(2σ^2)}

在可靠性试验中,正态分布常用于分析产品性能退化数据或测量误差。当样本量较大时,样本均值的分布近似正态分布(中心极限定理),这是置信区间计算的重要依据。

分布选择提示: 如果产品失效机理明确为随机失效(如电子元件),优先选择指数分布;如果存在磨损或疲劳(如机械轴承),优先选择威布尔分布;如果关注性能退化或寿命分布对称(如灯泡寿命),可考虑正态分布。
三种分布失效率 λ(t) 曲线对比 时间 t 失效率 λ(t) 指数分布(恒定) 威布尔 β<1(递减) 威布尔 β>1(递增) 正态(急剧上升) 注:威布尔分布通过 β 值可模拟指数(β=1)、早期失效(β<1)和磨损失效(β>1),是可靠性试验中最灵活的模型。

2.2 均值与方差

均值和方差是描述分布集中趋势和离散程度的两个核心数字特征。在可靠性试验中,它们直接关联到产品寿命的估计精度。

分布 均值 E(T) 方差 Var(T)
指数分布 1/λ 1/λ²
威布尔分布 η · Γ(1 + 1/β) η² · [Γ(1 + 2/β) - Γ²(1 + 1/β)]
正态分布 μ σ²

其中 Γ(·) 为伽马函数。从表中可以看出:

  • 指数分布的方差等于均值的平方,离散程度完全由均值决定。
  • 威布尔分布的均值和方差均受形状参数 β 和尺度参数 η 共同影响,β 越小,数据越分散。
  • 正态分布的均值和方差相互独立,σ 越大,寿命数据越分散,需要的样本量通常也越大。
工程直觉: 方差越大,意味着产品寿命的个体差异越大。在样本量设计时,需要更大的样本量才能以相同的置信度估计出均值的真实范围。反之,如果产品一致性很好(方差小),少量样本即可获得可靠结论。

2.3 置信区间概念

置信区间是统计学中用于量化估计不确定性的核心工具。在可靠性试验中,我们通常用样本数据来估计总体参数(如 MTBF、可靠度),并给出一个区间,该区间以一定的概率包含真实参数值。

2.3.1 置信区间的定义

一个 100(1-α)% 的置信区间表示:如果重复进行大量相同试验,大约有 100(1-α)% 的区间会包含真实的总体参数。α 是显著性水平,通常取 0.05 或 0.10,对应 95% 或 90% 的置信水平。

2.3.2 置信区间的一般形式

置信区间 = 点估计 ± 临界值 × 标准误
  • 点估计:由样本数据计算得到的参数估计值(如样本均值)。
  • 临界值:取决于置信水平和分布类型(如正态分布的 z 值或 t 分布的值)。
  • 标准误:点估计的标准差,反映估计量的抽样波动。

2.3.3 可靠性试验中的常见置信区间

参数 分布假设 置信区间形式(近似)
MTBF(指数分布) 卡方分布 [2T/χ²_{α/2, 2r+2}, 2T/χ²_{1-α/2, 2r}]
均值 μ(正态分布,σ 已知) 正态分布 [x̄ - z_{α/2}·σ/√n, x̄ + z_{α/2}·σ/√n]
均值 μ(正态分布,σ 未知) t 分布 [x̄ - t_{α/2, n-1}·s/√n, x̄ + t_{α/2, n-1}·s/√n]
可靠度 R(t)(非参数) 二项分布 基于 Clopper-Pearson 精确区间
常见误区: 置信区间不是“参数落在区间内的概率为 95%”。正确的解释是:在重复抽样中,95% 的置信区间会包含真实参数。对于一次具体的试验,参数要么在区间内,要么不在,不存在概率。
置信区间概念示意图(95% 置信水平) 真实参数 θ 包含 ✓ 包含 ✓ 不包含 ✗ 包含 ✓ 包含 ✓ 在 20 次重复试验中,约有 19 次(95%)的置信区间会包含真实参数 θ。

2.4 从分布到置信度:样本量设计的逻辑链

理解上述概念后,可以梳理出可靠性试验样本量设计的核心逻辑链:

  1. 确定目标分布:根据产品失效机理选择指数、威布尔或正态分布。
  2. 设定置信水平:通常为 90% 或 95%,对应 α = 0.10 或 0.05。
  3. 指定允许误差:即置信区间的半宽度(精度要求),例如 MTBF 估计误差不超过 ±20%。
  4. 估计方差:利用历史数据或工程经验估计分布参数的方差(如 σ² 或 β)。
  5. 计算样本量:基于置信区间公式反推所需样本量 n。