软件验证数据集（11462-3）：软件算对了吗
一、先讲个真事
那家厂的喷油嘴针阀偶件，配合间隙是核心特性，客户要求 Cpk ≥ 1.33。系统里长期显示 1.6 上下，看着挺舒服。
争议起于一次客户抱怨——装配线反馈卡滞率偏高。按 Cpk 1.63 推算，超差比例应该是百万分之零点几，跟实际对不上。
我们做了件很简单的事：把 ISO/TR 11462-3 里的 DS1 数据集塞进他们的系统。这套数据是标准公开的、带"标准答案"的考题，正态分布、无异常值，125 个数、子组 n = 5、共 25 组。答案是白纸黑字印在标准里的：Cp = 2.07，Cpk = 1.95。
系统跑出来：Cp = 2.31，Cpk = 2.18。
差 11.6%。 一查配置就明白了：那套软件的能力模块，默认用整体标准差（s，即 22514 口径里的性能算法）去算 Cp/Cpk，而不是用组内标准差。放在他们那批离散度大的实际数据上，这个错误的方向反过来——把 1.42 算成了 1.63。
说白了，软件没坏，是算法口径和你以为的不一样。而你之所以三年没发现，是因为从来没拿标准答案考过它。
二、核心知识点拆解
2.1 这份 TR 到底是什么
ISO/TR 11462-3 是 SPC 软件的"带答案考卷"。它给出 11 套参考数据集和公布结果，用来验证实现了以下标准的软件：
| 覆盖标准 | 验证内容 |
|---|---|
| ISO 7870-2 / GB/T 17989.2 | 控制限计算、失控（判异）检测 |
| ISO 22514-2 | 位置估计量、离散度估计量、能力指数 |
| ISO 22514-8 | 多状态过程能力（对应 DS11） |
11 套数据覆盖不同分布模型（ISO 22514-2 规定的 A/B/C/D 时变分布族，A1/A2 等后缀区分具体设定），里面有正态的，也有威布尔、瑞利这类右偏、厚尾的。标准里明说：选哪个分布模型由统计专家决定，标准不规定选模流程——它只负责给答案，不替你做判断。
这些数据集的"标准答案"是用不同程序多次交叉验证过的，电子版可从 ISO 官网获取。
2.2 DS1：基准中的基准
DS1 是正态、无异常值的过程数据，最适合当第一道考题。标准公布的关键值我整理成一张校验表：
| 项目 | 公布值 | 说明 |
|---|---|---|
| 样本构成 | 125 个，子组 n = 5（25 组） | 分布模型 A1（正态） |
| 规格上限 USL | 14.075 | — |
| 规格下限 LSL | 14.060 | — |
| 过程均值 μ̂ | 14.067958 | 位置估计（估计量编号 1） |
| 组内标准差 σ | 0.001207 | 离散度估计（编号 d = 5），能力计算用 |
| 子组内标准差均值 s̄ | 0.001128 | — |
| 子组极差均值 R̄ | 0.005900 | — |
| Cp | 2.07 | 过程能力指数 |
| Cpk | 1.95 | 最小过程能力指数 |
| CpkL / CpkU | 2.20 / 1.95 | 下／上单侧能力指数 |
| Xbar 图控制限（用 s） | LCL 14.066349 / UCL 14.069568 | 依 ISO 7870-2 |
| I-MR 个体图控制限 | LCL 14.06429 / UCL 14.07163 | 依 ISO 7870-2 |
| 中位数图控制限 | LCL 14.06593 / UCL 14.06972 | 依 ISO 7870-2 |
| s 图控制限 | LCL 0 / UCL 0.002356 | — |
能力指数按 ISO 22514-2 的方法计算，位置估计量取编号 1、离散度估计量取编号 5；控制限用 7870-2 的校正因子，精确分位数与查表法的偏差小于 0.1%。
2.3 三条容易被写错的算法分叉
跑过十几家厂的软件，出问题的地方翻来覆去就这几处：
分叉一：σ 到底用哪个。 组内标准差（s̄/c₄ 或 R̄/d₂）算的是 Cp/Cpk，整体标准差算的是 Pp/Ppk。软件如果偷懒统一用整体 s，过程一稳定两者接近，看不出来；过程一有组间漂移，Cp 立刻虚高或虚低。
分叉二：c₄ 无偏修正做没做。 σ̂ = s̄ / c₄，n = 5 时 c₄ = 0.9400。不做修正，σ 会小 6%，Cp 直接虚高 6%。
分叉三：R 图 / s 图的下限截零。 n ≤ 6 时 D₃ = 0、B₃ = 0，下限就是 0。有的软件老老实实算出个负数还画出来，图上就多一条谁也不理解的线。
三、计算过程：拿 DS1 手算一遍
咱不依赖任何软件，纯手算，看看 2.07 和 1.95 是怎么来的。已知：USL = 14.075，LSL = 14.060，μ̂ = 14.067958，σ = 0.001207，s̄ = 0.001128，n = 5。
第 1 步，公差带宽度 T = USL − LSL = 14.075 − 14.060 = 0.015
第 2 步，Cp Cp = T / (6σ) = 0.015 / (6 × 0.001207) = 0.015 / 0.007242 = 2.0713 → 2.07 ✅
第 3 步，上侧能力 CpkU CpkU = (USL − μ̂) / (3σ) = (14.075 − 14.067958) / (3 × 0.001207) = 0.007042 / 0.003621 = 1.9448 → 1.95 ✅
第 4 步，下侧能力 CpkL CpkL = (μ̂ − LSL) / (3σ) = (14.067958 − 14.060) / 0.003621 = 0.007958 / 0.003621 = 2.1977 → 2.20 ✅
第 5 步，Cpk Cpk = min(CpkU, CpkL) = min(1.95, 2.20) = 1.95 ✅
第 6 步，Xbar 图控制限（用 s 路径） UCL/LCL = μ̂ ± A₃·s̄，n = 5 时 A₃ = 1.427 A₃·s̄ = 1.427 × 0.001128 = 0.0016096 UCL = 14.067958 + 0.0016096 = 14.069568 ✅ LCL = 14.067958 − 0.0016096 = 14.066348（公布 14.066349，末位舍入差）✅
第 7 步，s 图上控制限 UCL_s = B₄·s̄，n = 5 时 B₄ = 2.089 = 2.089 × 0.001128 = 0.0023563 → 0.002356 ✅ LCL_s = B₃·s̄ = 0 × 0.001128 = 0 ✅
七步全对上。这就是一份完整的软件验证记录——打印出来签字，就是审核时最硬的证据。
反过来看开头那家厂：他们的软件跑出 Cp = 2.31。咱倒推一下它用的 σ：σ = T / (6 × 2.31) = 0.015 / 13.86 = 0.001082。这个 0.001082 既不是组内 σ 的 0.001207，也不是 s̄ 的 0.001128——它比 s̄ 还小。再算：0.001128 × 0.9400 = 0.001060，接近。也就是说，那套软件把 c₄ 修正乘反了（应该除以 c₄，它去乘了）。一个符号级的 bug，藏了三年。
四、常见误区提醒
- "进口软件不用验"。算法 bug 跟国籍没关系。开头那个 c₄ 乘除方向搞反的，就是套价格不菲的进口系统。验证成本半天，不验的代价是三年的假数据。
- "验一套 DS1 就够了"。DS1 是正态无异常值的温柔题。真正暴露问题的是 DS2、DS5、DS9 那些威布尔、瑞利分布的数据集——很多软件遇上非正态数据一律按正态硬算，Cpk 差得离谱。别只考语文，理化生都得考。
- "结果差一点点无所谓"。Cpk 从 1.42 变成 1.63，看着差 0.21，换算成不合格率是从每百万 20 件掉到每百万 0.7 件，差了近 30 倍。能力指数是指数尺度的，小数点后第二位就是钱。
五、小结 & 互动提问
11462-3 干的事很朴素：给你一套考题和答案，让你有能力对软件说"不"。咱做质量的天天要求供应商提供检验报告，凭什么自己的 SPC 软件就能免检？跑一遍 DS1，半天时间，换一个心里踏实。
互动提问：你厂的 SPC 软件，算 Cp 用的是组内 σ 还是整体 σ？现在打开设置页看一眼——如果找不到这个选项，那才是真该担心的事。
如需配套互动练习、培训课件或软件试用，可搜索「东方大易」或访问大易官网 www.eastdayi.com（苏州东方大易管理咨询有限公司，电话 0512-65561989）联系客服获取。