常规控制图(休哈特控制图)· 互动讲解

依据 GB/T 17989.2-2020《控制图 第2部分:常规控制图》(ISO 7870-2:2013, MOD) | 面向统计学小白,由浅入深 + 可拖拽演示

目录你想先看哪块?

1缘起:休哈特与控制图想解决什么

一句话版:控制图是 1924 年休哈特发明的一张“过程心电图”,用来实时判断——过程到底是在正常“呼吸”,还是生病了。

传统质检是事后全检、挑出废品,浪费且不经济。休哈特的思路是:一开始就预防。他把过程波动分成两类:

偶然原因(常见原因 / 随机原因)

过程固有的、始终存在的很多小因素叠加而成。单个都很小、无法逐个区分,但总和可测量。
对策:要消除得从系统层面改(属于管理层责任)。

可查明原因(特殊原因 / 系统原因)

可识别的异常:坏料、破损工具、程序缺陷、设备异常、环境变化……
对策:当场找到并消除(操作员/工程师责任)。

受控(统计控制状态) = 过程变异来自偶然原因。一旦确定这个水平,再出现的偏差就该去抓“可查明原因”了。控制图就是帮你区分这两类变异的图形工具。

标准引言里还特别强调:控制图不只是“报警”,更能帮你理解过程、识别改进机会。这正是它作为 SPC 核心工具的价值。

2控制图长啥样 & 3σ 原理

控制图 = 一条中心线 CL + 上下两条控制限 UCL/LCL + 按时间顺序打的点。限设在中心线两侧 3σ 处。

原理示意:拖动 k 看控制限如何变化

为什么是 3σ?不是 2σ 也不是 4σ?

休哈特用 3σ 是经济平衡:限太窄→老把正常当异常去查(浪费);限太宽→真出问题时查不出来(损失大)。在“绘图统计量近似正态”假设下:

落在 3σ 限内的概率 ≈ 99.7% → 受控时约 0.3% 的点会偶然跑出去
即:1000 个点里平均约 3 个会“虚惊一场”(第一类错误 α≈0.3‰)

标准还提到:有人偏好用 3.09 对应精确的 0.2% 单边概率;也有用实际非正态概率——但休哈特坚持 ±3σ 是经验性的“行动限”。

2σ 警戒限 & 两类错误

第一类错误(α)
过程其实受控,点却随机跑出限外 → 误判失控、白查一顿。
3σ 图下 α≈0.3%。常规图的设计原则就是“控住第一类错误”(第二类错误很难预先估计)。
第二类错误(β)
过程已失控,点却还在限内 → 没报警,继续产废品。
风险取决于限宽、子组大小、失控幅度——很难事先算准
2σ 警戒限:在中心线两侧 2σ 处画线,点落在 2σ 外算“即将失控的警告”,不一定立刻行动(有人会再采一组同样的子组确认)。

阶段 1(建立) vs 阶段 2(监控)

阶段 1(回顾阶段):用历史/新采集数据估计 CL 与限,反复“剔除异常→重算”,直到无报警、过程稳定可预测。此时限是试验性的。
阶段 2(监控阶段):把阶段 1 定下的 CL/限用于持续监控,快速抓偶发异常。不消除异常原因会让波动被高估、限变宽、灵敏度下降。

3控制图类型全景(点选找图)

先回答三个问题:数据是计量的还是计数的?是正态/位置还是不合格/缺陷?子组大小固定吗?标准图 2 就是这套选法。

参数“给定”还是“未给定”?

未给定(估计):过程 μ、σ 未知,用本过程数据估计 CL 和限,目标是把过程搞到受控。最常见情形。
已给定:μ₀、σ₀ 是目标值/长期值/规格值。用来保持过程在期望水平。公式里用给定值。
计量控制图总是成对出现:一张管“位置”(X̄/中位数/单值),一张管“散布”(R/s/移动极差)。先看散布图,再位置图——因为散布图给出的 σ 估计才能用来定位置图的控制限。

4计量控制图原理与因子表

最经典的是 X̄-R 图(n<10 用 R,n≥10 用 s)。下面把标准表 1 的公式和表 2 的因子讲清楚。

4.1 控制限公式(标准表 1,大白话版)

X̄ 图(估计) CL=X̄̄ ; UCL/LCL = X̄̄ ± A2·R̄ (或 ± A3·s̄)
X̄ 图(给定) CL=μ₀ ; UCL/LCL = μ₀ ± A·σ₀
R 图 CL=R̄ ; UCL=D4·R̄ ; LCL=D3·R̄ (n<7 时 D3=0,下限不画)
s 图 CL=s̄ ; UCL=B4·s̄ ; LCL=B3·s̄ (n<6 时 B3=0)
σ 估计 σ̂ = R̄ / d2 (或 = s̄ / c4)
提示:A2 = 3/(d2√n),A = 3/√n,A3 = 3/(c4√n)。这三个“A 因子”本质都是把 3σ 折算到对应统计量上。

4.2 标准差 s 图什么时候用?

子组小(n<10)用 R 图;子组大(n≥10)用 s 图——因为 n 大时极差估计 σ 的效率下降。有电子设备算标准差时,更倾向 s 图(见附录 B.1.2 电池示例)。

4.3 因子查表(标准表 2,可查)

注:本表 n=2..25 的 A/A2/A3/c4/d2/D3/D4 与标准表 2 一致;B3=max(0,1−3/(c4√n))、B4=1+3/(c4√n)、B6=c4+3√(1−c4²) 由定义算出。

5互动①:X̄-R 图计算器(轴承车削示例)

标准附录 B.1.1:水泵轴承直径,每小时测 1 子组、子组大小 5,共 25 组。μ、σ 未知 → 估计控制限。点“载入轴承示例”即可看到完整结果。

当前:含全部 25 组
X̄ 控制图(子组均值)
R 控制图(子组极差)

6互动②:单值 X + 移动极差 Rm 图

当“合理子组”不现实(如每批只测一次),就用单值图。用相邻观测的移动极差 Rm 估计波动。标准附录 B.1.3:脱脂奶粉湿度百分比。

单值 X 控制图
移动极差 Rm 控制图
注意:单值图对过程变化不如子组图敏感;且要求数据是“时序列”、相邻观测间过程没大变。标准特别点名:别把两次不连续化学反应的首尾凑成一对算移动极差。

7计数控制图原理(p / np / c / u)

计量图要“位置+散布”两张;计数图只要一张——因为计数分布只有一个参数(均值水平)。p、np 基于二项分布;c、u 基于泊松分布。

p 图(不合格品率)

比例 = 不合格数 / 检查数。n 可变也可用。公式:p̄ ± 3√[p̄(1−p̄)/n]。

np 图(不合格品数)

n 固定时。中心线 np̄,限 np̄ ± 3√[np̄(1−p̄)]。

c 图(缺陷数)

固定机会数下的缺陷次数(泊松)。限 c̄ ± 3√c̄。

u 图(单位缺陷数)

机会数变化时,每单位缺陷数(泊松)。限 ū ± 3√(ū/n)。
p 图(估计) p̄ = Σd / Σn ; 每子组 UCL/LCL = p̄ ± 3·√[ p̄(1−p̄)/nᵢ ] (n 变化时每子组单独算)
标准化 z 图 当 n 变化明显:画 z = (p−p̄)/√[p̄(1−p̄)/n],则 CL=0、UCL=+3、LCL=−3(与 n 无关)
n 变化的简化:若各子组 n 与目标 n 相差在 ±25% 内,可用平均 n 算一套统一限(标准 B.2.1 晶体管示例就这么处理)。下控制限为负时不显示(不合格率/缺陷数不能是负数)。

8互动③:p 图计算器(晶体管示例)

标准附录 B.2.1:无线电晶体管,每天随机抽一批检不合格品,26 天。n 不完全相同(135~165)。先算 p̄,再逐子组算限。

当前:含全部 26 组
p 控制图(不合格品率,n 变化时逐子组算限)

9控制程序 7 步(先 R 后 X)

以 X̄-R 图为例(标准第 7 章 + 图 6)。核心纪律:先确认散布受控,再判断位置

  1. 收集初步数据:标准操作下采合理子组,最少约 25 组,算每组的 X̄、R。
  2. 检查 R 图(或 s 图):画试验 CL/限,看有无出限点或异常模式。每个报警去查可查明原因。
  3. 剔除可查明原因 + 修改图:删掉被识别原因影响的子组,重算 CL/限重画;必要时重复。保留 ≥2/3 子组,不够就补采。
  4. 检查 X̄ 图:R 图受控后,散布稳定,再分析位置。画 X̄ 图限,同样剔除失控点重算。(构造 X̄ 图时也要剔除 R 图已删的子组)
  5. 持续监测(阶段 2):无报警即受控,用修订后的限监控未来;过程变了就调限。
  6. 判过程能力:受控后,用 R̄/d2 得 σ 估计,算 Cp/Cpk(见第 11 节)。
  7. 过程改进:能力不足就改进过程/重定中心/改规格,循环往复。
标准强调:X̄ 图对组间均值变化敏感,R/s 图对组内波动敏感;不先稳住散布就去读位置图,σ 估计不可靠

10检验模式(第 8 章 · 可点看)

点出限外只是“检验 1”。很多小偏移/趋势不够大、点不出去,但会呈现非随机模式。标准图 3 给出 4 种,西方电气规则有 8 种(附录 A)。

检验模式示意
区域划分:CL 两侧各分 A、B、C 三区,每区宽 1σ(共 6σ)。稳定过程 >2/3 的点应在 C 区;远少于就提示非随机模式。
代价:只用法 1 虚发率 3‰;同时用法 1~3,X̄/单值图虚发率升到约 10‰。补充规则提高小偏移检出力,但以更高虚警为代价。

11过程能力 Cp / Cpk

过程受控、可预测后,才能谈“能力”。标准第 9 章。

Cp = (USL−LSL) / (6σ) ← 只比“散布”和规格带宽,不看中心偏不偏
Cpk = min[ (USL−μ)/(3σ) , (μ−LSL)/(3σ) ] ← 再考虑中心位置
Cp < 1
能力不足
Cp = 1
勉强可接受
Cp ≥ 1.33
通常作最低可接受(留了采样波动余量)
顺序不能乱:先确认受控,再评能力。Cp 高也可能因中心偏移而大量超规格——所以要看 Cpk。实务中常要求 ≥25 子组证明受控后再算。

12前期准备与 3σ 替代原则

11 章前期准备要点

13.3 3σ 替代原则(小偏移怎么办)

X̄ 图对大幅偏移敏感,但对 ≤1.5σ 的小偏移表现差。标准给两条出路:

① 用补充检验模式(第 8 章)或改用 EWMA / CUSUM 图——对小偏移更灵敏(你之前做过的 EWMA 工具正是为此)。
② “三点中两点”准则:连续 3 点有 2 点落在 2σ 外即报警;此时把常规 3σ 限改成 中心线 ±1.78σ,虚发率与常规“点出 3σ 限”相同,但对中小偏移检出力大增。
13.2 数据相关:若观测自相关,X̄ 图基本方程 σ²(X̄)=σ²(个体)/n 不成立,不宜用常规法——应建模取残差,或用 X̄ 的总变差算限,并请教专家。

13标准示例一览(附录 B)

标准给了 8 个完整示例,覆盖各类图。下面速览,前几个在本工具的互动计算器里可直接复现。

B.1.1 X̄-R 图(μσ 未知)· 轴承直径

n=5,25 组。R̄=0.0177,X̄̄=14.07317;R 图受控,X̄ 图第 12 组失控→剔除后 X̄̄=14.07385、R̄=0.01800,过程受控。→ 见第 5 节计算器

B.1.2 X̄-s 图(μσ 给定)· 电池质量

μ₀=29.87g、σ₀=0.062g、n=5。X̄ 限=μ₀±Aσ₀=29.87±1.342×0.062;s 图 CL=0.94σ₀,UCL=1.964σ₀。过程受控。

B.1.3 单值 X + Rm 图(μσ 未知)· 脱脂奶粉湿度

25 批单测。X̄=3.44%,Rm̄=0.33%;X 限=3.44±2.66×0.33,Rm 限=3.267×0.33。→ 见第 6 节计算器

B.1.4 中位数 + R 图(μσ 未知)· DVD 厚度

n=5,20 组。中位数 X̃=11.70、R̄=5.90;中位数图限=X̃±A4·R̄=11.70±0.691×5.90。

B.2.1 p 图(未给定)· 晶体管不合格品率

26 天,n=135~165。p̄=233/3893≈0.0598;第 17、26 组超上限→剔除后 p̄=195/3596≈0.054,用平均 n=150 算统一限。→ 见第 8 节计算器

B.2.2 np 图(未给定)· 开关 100% 检测

n=4000 固定。np̄=10.76;限=10.76±3√[10.76×0.9973]→UCL≈20.59、LCL≈0.93。

B.2.3 c 图(未给定)· 录像带疵点

20 卷,每卷 350m。c̄=3.4;限=3.4±3√3.4→UCL≈8.9,LCL 为负不画。

B.2.4 u 图(未给定)· 轮胎单位缺陷

n=50 固定。ū=0.077;限=0.077±3√(0.077/50)。n 固定时也可用 c 图。
沟通