现在的位置: 首页 > 综合 > 正文

统计过程控制图SPC(2)

2018年10月22日 ⁄ 综合 ⁄ 共 3228字 ⁄ 字号 评论关闭

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

(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 3
B5 4 D1 2 D2 3 D3 3
D4 2 D9 4 M2 9 M6 13
L3 5 L4 4 L7 6 Z1 15
Z2 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 9
5 3 6 4 7 6 8 9
9 11 10 13 11 12 12 7
13 2 14 14 15 9 16 8
17 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 17
5 3 6 8 7 9 8 2
9 14 10 1 11 3 12 5
13 6 14 9 15 10 16 12
17 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 35
C4 1258 25
A7 1248 24
A1 1260 39
A5 1273 29
D3 1260 21
D6 1259 37
D1 1240 37
R4 1260 28
H7 1255 19
H2 1268 36
H6 1253 36
P4 1273 29
P9 1275 22
J7 1257 24
J2 1269 41
J3 1249 36
B2 1264 31
G4 1258 25
G6 1248 36
G3 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;

结果如下:

minitab也可以完成上述操作,但具体细节上有所差别,下面是minitab的控制图结果:

 

 

抱歉!评论已关闭.