当前位置：服务支持 > 软件文章 > 统计过程控制图SPC（2）进阶教程

统计过程控制图SPC（2）进阶教程

阅读数 9

文章提及的统计过程控制图是单一变量的控制图，它属于单值控制图。控制图的种类很多，我的老师在其书中如是分类：

（1）计量控制。对于产品的物理度量指标，如机器零件的尺寸，棉卷的重量，青砖的抗压强度等，在生产过程处于控制状态时，产品质量服从正态分布。根据这种指标进行抽样控制叫计量控制。常用的计量控制图有：单值控制图（即x控制图）、平均值控制图（即控制图）、极差控制图（即R控制图）、中位数极差控制图（即 -R控制图）等。

（2）计件控制。产品中具有某种标志的产品件数指标，如一批产品中的次品件数、废品件数、一等品件数等，服从二项分布。根据这种指标进行抽样控制叫计数控制。常用的计数控制图有不合格品数控制图（即Pn控制图）和不合格品率控制图（即P控制图）等。

（3）计点控制。对于产品上的疵点数，比如一件铸件表面的孔、砂眼等，一匹布上的挑纱、竹节、破眼等，均服从普哇松分布。根据这种指标进行的抽样控制叫计点控制图。常用的计点控制图有缺陷数控制图（即C控制图）和单位缺陷控制图（即控制图）。

这里涉及不同的分布函数，其中，正态分布中存在均值和方差，二项分布中的概率，普哇松分布（Poisson）的参数λ。三种分布对应存在着X-R,X-S,P,PN,C,U等控制图，控制图的控制线对应着相应的等式. 以二项分布为例，其控制线公式为：

minitab软件画统计过程控制图很方便，该软件有以上几种控制图菜单操作，简单易行。大伙可以在MINITAB软件的工具栏找到统计（S），进而进入控制图（C），控制图下拉菜单里存在具体的控制图。例如，X-R控制图的界面：

sas软件在统计过程控制运用方面专门有个模块：QC模块，该模块除了可以进行统计过程控制的功能外，还可以用来进行质量改进、实验设计、工艺改进、分析统计过程能力、开发和评估抽样计划等。就统计过程控制图来说，QC模块支持一下语句来完成相应图形。

BOXCHART box chart with optional trend chart

CCHART c chart

IRCHART individual and moving range charts

MCHART median chart with optional trend chart

MRCHART median and R charts

NPCHART np chart

PCHART p chart

RCHART R chart

SCHART s chart

UCHART u chart

XCHART X chart with optional trend chart

XRCHART X and R charts

XSCHART X and s charts

另外，SAS/QC User's Guede有详细信息，部分控制图案例摘录如下：

/*Example 13.8: Applying Tests for Special Causes*/data Trucks3;   input TruckID $ Defects @@;label TruckID='Truck Identification Number'       Defects='Number of Paint Defects';datalines;B1 12 B2 4 B3 4 B4 3B5 4 D1 2 D2 3 D3 3D4 2 D9 4 M2 9 M6 13L3 5 L4 4 L7 6 Z1 15Z2 8 Z3 9 Z7 6 Z9 8;run;ods graphics on;title1 'c Chart for Paint Defects in New Trucks';title2 'Tests=1 to 4';proc shewhart data=Trucks3;   cchart Defects*TruckID / tests = 1 to 4                            testlabel1 = 'Test=1 Signaled'                            testlabel2 = 'Test=2 signaled'                            odstitle = title                            zonelabels                            tabletests                            tablelegend;run;ods graphics off; /*Example 13.22: Applying Tests for Special Causes*/data Circuit3;   input Batch Fail @@;datalines;1 12 2 21 3 16 4 95 3 6 4 7 6 8 99 11 10 13 11 12 12 713 2 14 14 15 9 16 817 14 18 10 19 11 20 9;run;ods graphics on;title1'p Chart for the Proportion of Failing Circuits';title2 'Tests = 1 to 4';proc shewhart data=Circuit3;   pchart Fail*Batch / subgroupn = 500                       tests = 1 to 4                       zones                       zonelabels                       ltests = 20                       table                       tabletest                       tablelegend;run; /*Example 13.31: Applying Tests for Special Causes*/data Fabric3;   input Roll Defects @@;datalines;1 6 2 9 3 14 4 175 3 6 8 7 9 8 29 14 10 1 11 3 12 513 6 14 9 15 10 16 1217 11 18 4 19 9 20 4;run;ods graphics on;title1 'u Chart for Fabric Defects';title2 'Tests=1 to 4';proc shewhart data=Fabric3;   uchart Defects*Roll / subgroupn = 30                         tests = 1 to 4                         odstitle = title                         odstitle2 = title2                         tabletests                         zonelabels;run;ods graphics off; /*Example 13.38: Applying Tests for Special Causes*/data Tape;   input Sample $ WeightX WeightR;   WeightN=5;   label WeightX = 'Average Adhesive Amount'         Sample = 'Sample Code';datalines;C9 1270 35C4 1258 25A7 1248 24A1 1260 39A5 1273 29D3 1260 21D6 1259 37D1 1240 37R4 1260 28H7 1255 19H2 1268 36H6 1253 36P4 1273 29P9 1275 22J7 1257 24J2 1269 41J3 1249 36B2 1264 31G4 1258 25G6 1248 36G3 1248 30;run;title 'Tests for Special Causes Applied to Adhesive Tape Data';ods graphics on;proc shewhart history=Tape;   xrchart Weight*Sample / tests = 1 to 5                           odstitle = title                           tabletests                           zonelabels;run;ods graphics off;title 'Tests for Special Causes Applied to Adhesive Tape Data';proc shewhart history=Tape;   xrchart Weight*Sample / tests = 1 to 5                           odstitle = title                           tabletests                           nochart                           zonelabels;run;

结果如下: