CN102831105B - 一种用excel和minitab15软件编制*-r控制图用系数表的方法 - Google Patents
一种用excel和minitab15软件编制*-r控制图用系数表的方法 Download PDFInfo
- Publication number
- CN102831105B CN102831105B CN201210250483.1A CN201210250483A CN102831105B CN 102831105 B CN102831105 B CN 102831105B CN 201210250483 A CN201210250483 A CN 201210250483A CN 102831105 B CN102831105 B CN 102831105B
- Authority
- CN
- China
- Prior art keywords
- control chart
- coefficient table
- data
- chart coefficient
- control
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种用EXCEL和MINITAB15软件编制-R控制图用系数表的方法,以
Description
技术领域
本发明涉及一种用EXCEL和MINITAB15软件编制控制图用系数表的方法。
背景技术
目前,控制图是对过程质量加以测定、记录从而进行控制管理的一种用科学方法设计的图。图上有中心线(CL)、上控制界限(UCL)和下控制界限(LCL),并有按时间顺序抽取的样本统计量数值的描点序列,参见控制图示例图。
如果数据越多,分组越密,则机螺丝直径直方图的直方图也越趋近一条光滑曲线,如直方图趋近光滑曲线图所示。在极限情况下得到的光滑曲线即为分布曲线,它反映了产品质量的统计规律,如分布曲线图所示。在质量特性值为连续值时,最常见的典型分布为正态分布。例如机螺丝直径直方图中机螺丝直径的分布就是如此,它的特点是中间高、两头低、左右对称并延伸至无限。正态分布可用两个参数即均值μ和标准差σ来决定。
正态分布有一个结论对质量管理很有用,即无论均值μ和标准差σ。取何值,产品质量特性值落在μ±3σ之间的概率为99.73,于是落在μ±3σ之外的概率为100%一99.73%=0.27%,而超过一侧,即大于μ-3σ或小于μ+3σ的概率为0.27%/2=0.135%≈1‰,如正态分布曲线图。这个结论十分重要。美国休哈特就根据这一事实提出了控制图。控制图的演变过程参见控制图的演变图。首先把正态分布曲线图按顺时针方向转90°成下图(控制图的演变a图),由于上下的数值大小不合常规,故再把控制图的演变图上下翻转180°而成下图(控制图的演变b图),这样就得到一张控制图,具体说是单值(χ)控制图。换个角度再来研究控制图的原理。根据来源的不同,质量因素可以分成4M1E五个方面。但从对质量的影响大小来看,质量因素可分成偶然因素(简称偶因)与异常因素(简称异因)两类。偶因是始终存在的,对质量的影响微小,但难以除去,例如机床开动时的轻微振动等。异因则有时存在,对质量影响大,但不难除去,例如车刀磨损、固定机床的螺母松动等。
偶因引起质量的偶然波动(简称偶波),异因引起质量的异常波动(简称异波)。偶波是不可避免的,但对质量的影响微小,故可把它看作背景噪声而听之任之。异波则不然,它对质量的影响大,且采取措施不难消除,故在过程中异波及造成异波的异因是我们注意的对象,一旦发生,就应该尽快找出,采取措施加以消除,并纳入标准化,保证它不再出现。
偶波与异波都是产品质量的波动,如何能发现异波的到来呢?经验与理论分析表明,当生产过程中只存在偶波时,产品质量将形成某种典型分布。例如,在车制螺丝的例子中形成正态分布。如果除去偶波外还有异波,则产品质量的分布必将偏离原来的典型分布。因此,根据典型分布是否偏离就能判断异波,即异因是否发生,而典型分布的偏离可由控制图检出。在上述车制螺丝的例子中,由于发生了车刀磨损的异因,螺丝直径的分布偏离了原来的正态分布而向上移动,于是点子超出上控制界的概率大为增加,从而点子频频出界,表明存在异波。控制图上的控制界限就是区分偶波与异波的科学界限。
根据上述,可以说休哈特控制图的实质是区分偶然因素与异常因素两类因素。
由于控制图是通过抽查来监控产品质量的,故两类错误是不可避免的。在控制图上,中心线一般是对称轴,所能变动的只是上下控制限的间距。若将间距增大,则α减小而β增大,反之,则α增大而β减小。因此,只能根据这两类错误造成的总损失最小来确定上下控制界限。
在控制图原理时曾经提到点子出界就判异,有的读者可能疑惑,如果是生产正常、点子偶然出界呢?现在,如果控制图是根据两类错误造成的总损失最小来确定的,那么根据“点子出界就判异”这样的准则来做,即使有时判断错误,但从长远看仍是合算的。
休哈特控制图(简称休图)的设计思想是先确定第I类错误的概率σ,然后再根据第Ⅱ类错误的概率β的大小来考虑是否需要采取必要的措施。通常σ取为1%,5%,10%。为了增加使用者的信心,休哈特将σ取得特别小,小到2.7‰~3‰这样,对于“点出界就判异”这条判异准则来讲,虽不百发百中,也是千发九九七中了。但σ小,β就大。为了减少第Ⅱ类错误,对于控制图中的界内点增添了第Ⅱ类判异准则,即“界内点排列不随机判异”。于是判断异常的准则就有两大类:
(1)点子出界就判断异常。
(2)界内点排列不随机判断异常。
其中,第(2)类准则是防止β大的。
休图的设计并未根据两类错误所造成的总损失最小这点来进行。从80年代以来,经济质量管理(EQC)兴起,学术代表人物是德国乌尔茨堡大学EQC中心的冯·考拉尼教授。EQC强调经济上最优,所以控制图设计的发展趋势之一就是根据两类错误所造成的总损失最小这点来确定控制界限。
4(均值-极差)控制图的原理及制作
对于计量值数据,(均值一极差)控制图是最常用、最重要的控制图,因为它具有下列优点:
a)适用范围广。对于图而言,计量值数据x服从正态分布是经常出现的。若x非正态分布,则当样本大小n≥4或5时,根据中心极限定理,知道近似正态分布。对于R图而言,通过在计算机上的统计模拟实验证实,只要总体分布不是太不对称的,R的分布没有大的变化。这就从理论上说明了 图适用的范围广泛。
b)灵敏度高。图的统计量为均值,反映在x上的偶然波动是随机的,通过均值的平均作用,这种偶然波动得到一定程度的抵消;而反映在x上的异常波动往往是在同一个方向的,它不会通过均值的平均作用抵消。因此,图检出异常的能力高。至于R图的灵敏度则不如图高。
现在说明一下图的统计基础。假定质量特性服从正态N(μ,σ2),且μ、σ均已知。若x1、x2、…xn是大小为n的样本,则样本均值为
由于服从正态分布N(μ,σ2/n),样本均值落入下列两个界限
间的概率为1-α。因此若μ与σ已知,则式(1)与式(2)可分别作为样本均值的控制图的上下控制界限。如前述,通常取Zα/2=3,即采用3σ控制界限。当然,即使x的分布是非正态的,但由于中心极限定理,上述结果也近似成立。
在实际工作中,μ与σ通常未知,这时就必须应用从稳态过程所取的预备样本的数据对它们进行估计。预备样本通常至少取25个(根据判稳准则,最好至少取35个预备样本)。设取m个样本,每个样本包含n个观测值。样本大小n主要取决于合理分组的结构,抽样与检查的费用,参数估计的效率等因素,n通常取为4、5或6。令所取的m个样本的均值分别为
则过程的μ的最佳估计量为总均值即
于是可作为图的中心线。
为了建立控制界限,需要估计过程的标准差σ可以根据m个样本的极差或标准差来进行估计。应用极差进行估计的优点是极差计算简单,所以至今R图的应用较s图为广。
现在讨论极差法。设为一大小为n的样本,则此样本的极差R为最大观测值xmax与最小观测值xmix之差,即
R=xmax-xmix(4)
若样本取自正态总体,可以证明样本极差R与总体标准差σ有下列关系:令W=R/σ,可以证明E(W)=d2,为一与样本大小n有关的常数,于是,σ的估计量为
令m个样本的极差为R1、R2、…、Rm,则样本平均极差为
故σ的估计量为
若样本大小n较小,则用极差法估计总体方差与用样本方差去估计总体方差的效果是一样的。但当n较大,如n>10,则由于极差没有考虑样本在xmax与xmix之间的观测值的信息,故极差法的效率迅速降低。但在实际工作中,R图一般取n=4、5或6,所以极差法是令人满意的。
若取μ的估计量为σ的估计量为则图的控制线为
式中
为一与样本大小n有关的常数,这可以通过查找计量值控制图系数表。
由上述,已知样本极差R与过程标准差σ有关,因此可以通过R来控制过程的变异度,这就是R图。R图的中心线即
为了确定R图的控制界限,需要对σR进行估计。若质量特性服从正态分布,令W=R/σ,可以证明σw=d3(d3为一与样本大小n有关的常数),于是从R=Wσ知σR=σwσ=d3σ。由于σ未知,故从式(6)得σR的估计量为
根据上述,得到R图的控制线如下:
令D3=1-3d3/d2,D4=1+3d3/d2,则代入上式后,得R图的控制线为
式中,系数D3、D4可以通过查找计量值控制图系数表。
随着现代科学技术的发展和计算机技术的广泛应用,质量控制领域对计算机的依赖程度也越来越高,广大工程技术人员可以充分利用EXCEL和MINITAB15强大的图表和数据处理功能为质量控制过程服务。
EXCEL和MINITAB15作为功能极强的数理统计计算软件将较复杂的编程简化为方便的菜单和工具栏操作,使数据处理及分析问题变得较为直观,EXCEL和MINITAB15软件具有丰富的库函数和强大的计算功能以及卓越的图表功能,灵活运用EXCEL和MINITAB软件可以极大地提高公式解决问题的能力。
质量控制过程计算涉及面广,内容庞杂,采集的数据多、处理复杂、计算量大,从现场采集的数据和模拟生成的数据找出相关量的关系以及服从的规律,需要对数据进行正确的整理分析和归纳计算。传统的控制图用系数表计算方法,手工计算量大、庞杂且不够精确,将EXCEL和MINITAB15软件应用于控制图用系数表编制,设计好数学模型,只要在工作表的单元格中输入公式,求解变得方便快捷,思路清晰、准确度高。
发明内容
本发明所解决的技术问题在于提供一种克服质量控制计算涉及面广、内容庞杂、采集的数据多、计算量大、且不够精确的特点,利用EXCEL和MINITAB15软件编制控制图用系数表的方法,在相应的工作表的单元格中输入相应数据,即可得到优化的计算结果。
本发明所解决的技术问题采用以下技术方案来实现:
一种利用EXCEL和MINITAB15软件编制控制图用系数表的方法,选择不同的样本,求取控制图用系数表中的各项系数,其特征在于:
(1)建立数学模型
d2=R均/бR
d3=d2*бR/R均
D4=(1+3*d3/d2)
D3=(1-3*d3/d2)
(2)依据(1)建立数学模型编制多种标准偏差区间范围的控制图用系数表的方法,其步骤如下:
1首选必须获得用于编制控制图用系数表中A2、d2、d3、D4、D3系数的必要的数据源,取得数据源的途径有以下两种:
A从制造过程现场获得;
B应用MINITAB15软件生成随机数据;
2采用1中的A或B骤获得的数据源,按照(1)所列的公式,分别以2,3,…,25,为样本数,分别计算A2、d2、d3、D4、D3;
3然后依次按照列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的表格依次填入;
4针对同一区间的观测值、标准偏差可以反复取样或反复生成模拟随机数据,按照1至3步骤,反复计算编制列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的控制图用系数表;
5将4步骤编制的不同的控制图用系数表中属于同一区间数据源的控制图用系数表数据进行数据叠加处理,编制最终的控制图用系数表。
本发明的有益效果是采用从制造现场取样或应用MINITAB软件生成随机数据,应用控制图用系数表原理,编制不同观测值和标准偏差区间范围的控制图用系数表;对应不同的观测值和标准偏差的数据特征,使用不同的控制图用系数表,从而精确估计控制图所对应的真实的标准偏差,达到改善工序控制质量的目的。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施例,进一步阐述本发明。
一种利用EXCEL和MINITAB15软件编制X-R控制图用系数表的方法,选择不同的样本,求取控制图用系数表中的各项系数,
(1)建立数学模型
d2=R均/бR
d3=d2*бR/R均
D4=(1+3*d3/d2)
D3=(1-3*d3/d2)
(2)依据(1)建立数学模型编制多种标准偏差区间范围的X-R控制图用系数表编制方法,其步骤如下:
1首选必须获得用于编制控制图用系数表中A2、d2、d3、D4、D3系数的必要的数据源,取得数据源的途径有以下两种:
A从制造过程现场获得;
B应用MINITAB15软件生成随机数据;
2采用1中的A或B骤获得的数据源,按照(1)所列的公式,分别以2,3,…,25,为样本数,分别计算A2、d2、d3、D4、D3;
3然后依次按照列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的表格依次填入;
4针对同一区间的观测值、标准偏差可以反复取样或反复生成模拟随机数据,按照1至3步骤,反复计算编制列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的控制图用系数表;
5将4步骤编制的不同的控制图用系数表中属于同一区间数据源的控制图用系数表数据进行数据叠加处理,编制最终的控制图用系数表。
(3)具体计算过程
1建立控制图用系数表中各项系数的数学模型,输入相应的计算公式:
在AB2-AB25单元格中输入A2的计算公式
在AC2-AC25单元格中输入d2的计算公式
在AD2-AD25单元格中输入d3的计算公式
在AE2-AE25单元格中输入D4的计算公式
在AF2-AF25单元格中输入D3的计算公式
2应用MINITAB软件中随机数据生成工具,产生符合正态分布特征的随机数据,粘贴入EXCEL工作表中B2-Z1000单元格内。
3依次在AG-DX列中分别输入б、R均、R和бR的计算公式。
4针对2,,25不同的样本,分别计算不同的б、R均、R和бR的数值。
5分别代入A2、d2、d3、D4、D3的计算公式进行计算;
6生成不同标准值和标准偏差区间的控制图用系数表。
项目达到的技术考指标:
通过将新编控制图用系数表与传统控制图用系数表系数进行对比,取样产品的观测值在0-1之间,标准偏差在0.1-0.01之间;随着观测值和标准偏差的增大,控制图用系数表的系数A2或d2逐渐增大或缩小;当观测值1-10,标准偏差在0.1-0.01之间;随着观测值和标准偏差的增大,控制图用系数表的系数A2或d2会继续增大或缩小,尤其是d2缩小幅度最大;当观测值10-100和100-1000,标准偏差在0-1和1-100之间,随着观测值和标准偏差的增大,控制图用系数表的系数A2或d2变化的速度会相应减缓。传统的控制图用系数表适用于观测值在0-1之间,标准偏差在0.1-0.01之间,而针对观测值1-10,标准偏差在0.1-0.01之间、观测值10-100之间,标准偏差在0-1之间、观测值100-1000之间、标准偏差在1-100之间,则需要重新编制适用于相应区间范围的控制图用系数表。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (1)
1.一种利用EXCEL和MINITAB15软件编制控制图用系数表的方法,选择不同的样本,求取控制图用系数表中的各项系数,其特征在于:
(a)建立数学模型
d2=R均/бR
d3=d2*бR/R均
D4=(1+3*d3/d2)
D3=(1-3*d3/d2)
(b)依据(a)建立数学模型编制多种标准偏差区间范围的控制图用系数表的方法,其步骤如下:
首先必须获得用于编制控制图用系数表中A2、d2、d3、D4、D3的必要的数据源,取得数据源的途径有以下两种:
A从制造过程现场获得;
B应用MINITAB15软件生成随机数据;
其次采用A或B获得的数据源,按照(a)所列的公式,分别以2,3,…,25,为样本数,分别计算A2、d2、d3、D4、D3;
然后依次按照列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的表格依次填入;
针对同一区间的观测值、标准偏差可以反复取样或反复生成模拟随机数据,按照(b)步骤,反复计算编制列为2,3,…,25样本数,行为A2、d2、d3、D4、D3的控制图用系数表;
将按照(b)步骤编制的不同的控制图用系数表中属于同一区间数据源的控制图用系数表数据进行数据叠加处理,编制最终的控制图用系数表。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210250483.1A CN102831105B (zh) | 2012-07-19 | 2012-07-19 | 一种用excel和minitab15软件编制*-r控制图用系数表的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210250483.1A CN102831105B (zh) | 2012-07-19 | 2012-07-19 | 一种用excel和minitab15软件编制*-r控制图用系数表的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102831105A CN102831105A (zh) | 2012-12-19 |
CN102831105B true CN102831105B (zh) | 2016-01-13 |
Family
ID=47334250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210250483.1A Expired - Fee Related CN102831105B (zh) | 2012-07-19 | 2012-07-19 | 一种用excel和minitab15软件编制*-r控制图用系数表的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102831105B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223909A (zh) * | 2014-06-12 | 2016-01-06 | 汕尾市智能岛电子科技有限公司 | 一种机床运动控制系统 |
CN105912817A (zh) * | 2016-05-06 | 2016-08-31 | 姜凌宇 | 一种基于Excel软件判定混凝土强度保证率的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6496347B1 (en) * | 2000-03-08 | 2002-12-17 | General Electric Company | System and method for optimization of a circuit breaker mechanism |
CN101459399A (zh) * | 2007-12-13 | 2009-06-17 | 现代自动车株式会社 | 混合动力电动车中控制电机转矩的方法 |
-
2012
- 2012-07-19 CN CN201210250483.1A patent/CN102831105B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6496347B1 (en) * | 2000-03-08 | 2002-12-17 | General Electric Company | System and method for optimization of a circuit breaker mechanism |
CN101459399A (zh) * | 2007-12-13 | 2009-06-17 | 现代自动车株式会社 | 混合动力电动车中控制电机转矩的方法 |
Non-Patent Citations (1)
Title |
---|
均值-极差控制图在阀门质量控制中的应用;郑明君 等;《价值工程》;20111108;第26-27页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102831105A (zh) | 2012-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104808587B (zh) | 一种基于机加工设备运行状态的稼动率统计方法 | |
Bhinge et al. | Toward a generalized energy prediction model for machine tools | |
Shu et al. | A new EWMA chart for monitoring process dispersion | |
CN104808585B (zh) | 一种机床健康状态快速检查方法 | |
CN101520652B (zh) | 一种数控装备服役可靠性的评估方法 | |
CN111046532B (zh) | 一种基于肘形判据的同调机群聚类识别方法 | |
Guh | Robustness of the neural network based control chart pattern recognition system to non‐normality | |
US11982996B2 (en) | Method and apparatus for configuring processing parameters of production equipment, and computer-readable medium | |
CN107944005B (zh) | 一种数据展示方法及装置 | |
CN103345200A (zh) | 一种基于广义区间的切削加工颤振辨识方法 | |
CN111724203B (zh) | 基于多重回归模型的产品加工价格预估方法及系统 | |
CN102831105B (zh) | 一种用excel和minitab15软件编制*-r控制图用系数表的方法 | |
Ghute et al. | A nonparametric signed-rank control chart for bivariate process location | |
Khaw et al. | A re-evaluation of the run rules control charts for monitoring the coefficient of variation | |
Lim et al. | AN ALTERNATIVE DESIGN FOR THE VARIABLE SAMPLE SIZE COEFFICIENT OF VARIATION CHART BASED ON THE MEDIAN RUN LENGTH AND EXPECTED MEDIAN RUN LENGTH. | |
Movaffagh et al. | Monotonic change point estimation in the mean vector of a multivariate normal process | |
CN109887253B (zh) | 石油化工装置报警的关联分析方法 | |
CN113836774B (zh) | 一种机理与数据融合的金属切削仿真过程不确定校准方法 | |
Filipovich et al. | Simulation model of selective assembly of two parts with sorting by the estimated values | |
Klyachkin et al. | Analysis of stable functioning of objects using machine learning | |
Nimbale et al. | A neural network based individual control chart | |
CN104050277A (zh) | 一种基于聚类分析的3mad-mmmd过失误差侦破方法 | |
CN110543908A (zh) | 一种基于动态观测窗口的控制图模式识别方法 | |
CN109284320A (zh) | 大数据平台上的自动回归诊断方法 | |
JPH04223560A (ja) | 品質改善方法及びそのシステム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190423 Address after: 230601 Kechuang Center 3F326, Jaicui Road 393, Peach Blossom Industrial Park, Jingkai District, Hefei, Anhui Province Patentee after: Hefei HEMPEL Mdt InfoTech Ltd Address before: 230000 Jade Road 367, Hefei Economic and Technological Development Zone, Anhui Province Patentee before: Jiu Jianjun |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160113 Termination date: 20210719 |