CN1319212A - 用于统计过程控制的多维方法和系统 - Google Patents
用于统计过程控制的多维方法和系统 Download PDFInfo
- Publication number
- CN1319212A CN1319212A CN99811236A CN99811236A CN1319212A CN 1319212 A CN1319212 A CN 1319212A CN 99811236 A CN99811236 A CN 99811236A CN 99811236 A CN99811236 A CN 99811236A CN 1319212 A CN1319212 A CN 1319212A
- Authority
- CN
- China
- Prior art keywords
- control
- observation
- designator
- relevant
- point
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Feedback Control In General (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
本发明涉及一种过程统计方法和使用该方法的任何系统,其基于对输入、输出和控制的指示和测量数据的读取以及对运行参数的读取,并且能够在和其取样本相关的参考框架中用观察点表示出来。在本发明中,观察值被转化,以使结果值和多维的高斯分布模型一致;这些被转化的观察点在多维空间中,每一个维度与一个测量量值相关;失控观察点用同样方法计算,在某一特定方向上被识别;该方向和上述过程的偏移的一般原因相关,并且,每个“观察和异常对”和指示符相关,为了给出零个,一个或更多的可能和观察相关的原因。
Description
本发明涉及使用基于多维数据处理的统计过程控制的方法和系统。
其目的首先是当过程离开保证其生产质量要求的“正常”运行状态时引发警告,其次是为识别异常情况的可能原因提供建议。
统计过程控制(SPC)目前在世界各国(主要是工业化国家)众多领域的工业生产中正被广泛使用,例如工程学,电子学,化学,医药,农业食品,塑料材料,等等。
其目的是通过检查制造过程本身,而不仅仅检查产品的特征来保证产品质量。统计过程控制在实施“零缺陷运动”时和当企业寻求符合国际质量保证标准(国际标准化组织9000)时已是必不可少的阶段。
其技术目标是在不合格产品被制造之前检测出制造过程中可能的偏移并采取补救措施。
这种方法的使用现在已经从产品制造的领域扩展了,并包含了服务业领域(银行业务,保险,咨询,…)。
当过程启动时(参照图1),跟踪与同一过程相关的各种监测手段(指示符):输入特点(原材料);输出特点(产品);过程操作参数。每个观察单元(检测瞬间和生产元素)和多个测量获得的数值相关联,于是,使这组数据用所量度的多维的空间中的一点来表示。
统计过程控制的一般操作是通过追踪多个控制图来监控过程。控制图是观察的量值的变化路径的图形表示。同时也代表在一次测量中预先定义的控制极限值(见图2)。然后每个控制图都独立地被解释,且独立地引发警告。
目前使用的有各种控制图(象“Shewart,CuSum,EWMA,Me”)。上述后三种图被认为在检测小量的“偏移”方面优于第一种。
一般地,控制图被用于成组的数据:通过标记多次成组测量获得的平均值,可以较好地探测小量的偏移,另外,值的分布也与基于这种方法的正常情况下的假定更相吻合。通过标记每组数据的变化和延展,有可能探测到由于某种原因,测量可变性的有所提高。
通常的方法包括同时和独立地监视多个控制图,这种方法比较笨拙,而且对多维的统计数据控制中不太有效:
●这种方法有时发出太很多的假警告导致不必要的调整增多;这时,就需要很快地重新设定,从而导致整个过程控制需要多次调整,处于一种既混乱又昂贵的方式下运行状态。
●它检测到真实的异常情况太晚,并且
●该方法在异常原因不和测量直接相关联时,很难检测到它。于是,该方法鼓励了采用多样的测量,这样的方法很昂贵,并引入了多种控制图。
本发明的方法和系统寻求解决那些缺陷:该方法是一种统计过程控制方法,其基于读取输入、输出和控制端上的指示符和测量数据以及所述过程中的运行参数,可以使用基准结构中的观察点来表示它,其中基准结构将它们的值和它们的样本指标联系起来。根据本发明:
a)对观察值进行转化以使得结果值与多维的高斯分布模型一致,构成与本方法其余部分的观察点相对应的数据;
b)上述观察点位于多位空间,每个维度都和测量量值相关;
c)观察点中,那些与过程正常运行对应的所谓“控制下的”点和那些所谓“失控的”点有显著不同;
d)计算出控制下的点的分布中心作为控制下的观察点的点集之重心;
e)辨别出失控的观察点集中在“控制下的”点的分布中心的某一特定方向上;
f)这一方向与上述过程的通常的偏移原因相关;
g)每一个观察点和异常方向对与指示符相关,以便报告零个、一个或更多的所观察到的异常原因;并且
h)当以这种方式对一种异常情形进行分析后,将激发警告信号,并且,在工业生产过程中用这种办式检测到的偏移将被补救。
上述被检测的观察点重心对应于一个特定点,该点的分量是被检测的观察点的分量的平均值。
统计过程控制检测用传统方式,用统计规则或者用取样方法,有规律地观察p个连续的量值y1,y2,…yp。这些量值能同等地很好地反映原材料的特点,制造的产品的特点或制造过程的运行参数的特点。在某一给定瞬间,将这p个测量数据所组成的p-维向量表示为y,并作为过程观察向量。这个向量的终点是过程观察点,而起点是所用的基准结构的原点。
显然,在这段文章里“瞬间”的概念超过严格地时间概念:与同一个“瞬间”相关联的测量数据在可能的情况下应是同一批产品或同一制造单位的生产过程的参数测量数据。为了定义哪些测量数据和同一个相关瞬间相关,制造过程的完备的可跟踪性是必要的。
当过程是“控制下的”,y的值在不同的连续瞬间t0,t0+1,t0+2,…有y0(其是期望的目标)的微小变动,以保证生产达到令人满意的质量。这个微小的变化是由于原材料(原材料的硬度,原材料组份的化学成分,供应商,…)特性,或是环境因素(温度,湿度,…),或是过程因素(机器的设定,操作者的专注程度…)的随机变化而产生的。这些特性对y的一项或多项分量产生影响,并且被表示为z1,z2,…zm,这些变量的总体被表示为向量z。向量z在此被称为过程解释向量。
由于把过程特性看作yj过程的观察变量,必须在每个“瞬间”估值该变量。
由于把过程特性或输入特性看作过程的原因变量zk,必须由系统部的因素对其修改,如:故意或非故意的人为行动,环境变化,磨损或老化。一般的,由于费用或适用性的原因,这些变量不是在每个瞬间都进行测量(不然也将表示为变量yj),在这种意义上,它们组成了隐藏变量,影响着过程的行为。测量这些隐藏变量通常既昂贵、费时又不精确,仅仅在异常出现时才进行。
如果其可能值是数字表示的,而且是属于一个已知的数值范围(温度,压力,…)该变量可能是定量的。如果可能值是数字的或只有有限的个数(提供者,操作员,机器,…),改变量可能是定性的。本发明中的方法和模型假定y的分量全部是定量的。
过程的相同特征(例如炉子的控制温度)既可以作为z的分量又可以作为y的分量出现。
y和z之间的依赖关系设立如下模型:
y=f(z,t)+ε
其中t表示为观察瞬间,ε表示p维度上的随机向量,其平均值设定为0,并且有其协方差矩阵为∑ε。f是一向量函数,含有p个分量f1,f2,…fp,并有公式yi=fj(z,t)+ε
y的分量之间象z的分量之间一样相互关联。
在稳定条件下完美而稳定的过程的应有下列关系:
f(z,t)实际上与t无关,
每个“原因变量”zk都保持在固定值z0 k上,并且
y可用如下稳定过程的模型表示:
y=f(z0)+ε
实际中,定量的原因变量的不可能完全决定下来:变量z0已经加上了0期望值的随机误差和协方差矩阵∑e的随机误差。所以,模型变为
zk=z0 k+ek
y=f(z0+e)+ε
本发明建议的方法应用于以下统计过程控制通常所应用的环境中:
·函数f未知;
·解释变量zk无法全部识别;并且,
·n次观察由“瞬间t=1,…,t=n时的变量y1…yp产生。
这些观察结果以一个n行p列的矩阵y的形式表示。Yj代表Y的第j列;这一列的第i个元素表示为yi j,代表在t=i的瞬间变量yj的观察值。在t=i瞬间变量y1,…,yp的观察向量表示为yi。
观察向量yi的权重为pi,一般地等于1/n。以权重pi作为对角项的对角线矩阵(n,n)表示为Dp。
●当过程在适当的控制下时:
y0等于y的数学期望E[y];y0非常接近观察均值my(my代表由矩阵Y中P列Y1…Yp的均值y1…yp构成的向量)。
随机向量y的协方差矩阵不随时间变化。
●当过程发生偏移时,观察值yi偏离目标值y0很远。
这样的行为导致:
一个或多个原因变量zk的中心值随时间变化;或者
如果偏移是定性的,涉及原因变量zk,导致从z0k变化到z1k,则分布中心将从y0移动到y1。观察点沿y1-y0方向移动。
如果偏移是定量的,涉及zk,并且过程不是在y0周围关于z1k不稳定,则可以假定每个函数fj关于zk都有偏导数fjk。这种变化的演算结果表示了第一个顺序:
●均值中的偏移,z1k=z0k+d意味着观察点中心从y0沿着的偏导数向量定义的方向移动,该量为
(flk(y0),…,fjk(y0),…,fpk(y0));
●随机量ek的可变性的增加将增加观察点yi的分布在同一方向上延展:
(flk(y0),…,fjk(y0),…,fpk(y0))
●随机量ej的可变性的增加将使观察点yi的分布在第j个基本向量方向上延展:
(0,…,0,1,0,…,0)
l j p
本发明寻求管理上述过程的方法,以达到下述的技术目标:
识别出与已经在历史纪录识别出的任何偏移相关的方向,为每次观察定义参数,以便计算出这些方向的靠近指示符。
在过程控制的运行阶段:
探测是否在最后的观察中找到了运行过程的任何飘移,然后通过检查靠近指示符,识别出已发现与观察点最接近的原因的方向,从而可以给出偏移的大概原因;
在两个阶段中,给出图形表示,以便可以更快地对情况作出整体的估价。
为了使系统适应各种过程的特殊特性,为系统中的每种处理过程推荐了多个版本。
本发明提供一种统计过程控制的方法和使用该方法的任何系统,该方法基于在输入、输出和控制上检测到的指示和测量数据以及所述过程的操作参数,并且包括对按这种方式获得的数值进行的各种处理:该方法最好用计算机来进行计算机处理,完全自动地帮助使用者完成处理步骤。
这种操作包括两个阶段:
◆学习阶段:包括进行过程的历史数据分析;历史数据由一套与过程的运行有关联的量值构成;这些量值在连续的瞬时或者以取样的方式被测定。
这分析的目的是如下:
●为了识别
过程中偏移成为主流的瞬间;并且
●为了估计定义每个识别的特殊原因相关方向的参数的值
◆过程-追踪阶段
在这个阶段,系统从过程接受测量及数据(经过直接相联的传感器或通过人工输入)。当过程发生偏移的时候,它能激发警告。当为此次观察每个有关原因的量值都被估计,偏移的原因也被给出。这些原因是从以前的方法分析得出的原因中挑选出来的。如果没有建议的原因,将由人来检查实际过程,识别原因。新识别出的原因将被集成到系统中,以便以后发生时自动识别。
本发明系统的适当性和独创性来自于:
提出的模型和其应用的范围的适当;
采用了适当的图形表示方法:控制图用概率阈值作为刻度;并且
A1-数据的预先变换;
A2-识别控制下的测量数据和失控的测量数据;
A3-识别特殊原因类型和相关方向参数;
A4-在特殊原因指示符上建立和检查控制图;并且
A5-解释特殊的原因
这些步骤将和在实际的过程中那些不断向系统中引入新的信息的专家们一起互动的实施。(通过观察认定失控状态,联系多种原因,辨别可疑的测量数据,…)利用系统提供的部分数据来提取历史数据的解释。从某一步骤获得的结果经常可以重复一次或多次以前的步骤,只是修改可选步骤或某些参数。
过程的历史数据的特性由n行p列的矩阵Y描述。在第i行第j列交叉点的元素Yi j代表变量yj在观察瞬间i的值。
步骤A1:数据的预先变换
这涉及对原始观察数据yi j;变换因而采用的处理过程的值更加接近基于前述的方法的模型,所以提高了结果的适用性。
最为广泛和有用的变换是对每一个观察变量分开进行的,是将观察原始值yi j变换为新值
。以下提到一些变换方法,但不是全部。
其中w是每个数据组的大小。
计算连续的量度之间的差,以去掉量度之间的相关性:
变换给出离高斯分布“更近”的分布。例如:
其中FG -1是高斯分布函数的反函数,
假定为变量yj的分布的函数。
集中随机变化的推定的中心,在这种情况下:
中心值Ci j的选择依赖于实际情况和有关过程的知识。于是有:
Ci j=控制下变量yj的平均值;
Ci j=变量yj的目标点y0j,对时间是一常量;或者
Ci j=变量yj的目标点y0j,是时间变量(小的连已知的和可接收的变化…)。
其中,si j代表这次测量的挑选单元,例如:
si j=控制下的观察点的标准差yj。
步骤A2:分辨控制下的测量和失控的测量
这包括获取历史测量数据和分辨那些不是控制下的过程正常运行对应的点。
开始时,所有被过程控制负责人认为是失控的点都被作上标签,被认为是失控的。
然后,基于仍然被认为是在控制下的观察点建立一个HotellingT2表(见Douglas C.Montgomery的著作,题为“统计质量控制介绍”第二版,Wiley出版社1991年印刷)。在控制极限上限(极限值为α)以外的观察点被认为是失控点。α代表假预警的可接受概率。一般,α的值应在0.05-0.001之间。
这一步骤需要重复,每一次迭代,新的观察点将被标上失控的标签。
当在T2表的控制以外的观察点的数目达到假预警的可接受概率α的时候迭代停止。
这种一致性通过一种传统的假说试验来评价:给定超出极限值的观察频率,假说通过测量数据“出界”的概率小于等于α来检验。(见G.Saporta的著作,题为“概率分析和统计”,Technip出版社1990年出版)
A3-辨别特殊原因种类和相关方向参数。
这个想法是对于历史上在前续步骤中标记的失控的每个观察点都使其与一个很好标识异常原因(称为特殊原因)相关。
这种识别必须运用过程控制负责人的知识来进行。而且,一种很好的自动分类方法可以应用于那些被认为是失控的点,从而指导人们:这种分类方法可以汇总观察结果,把经常出现的原因和某种异常状况相对应。即:将在多维空间中离共同方向接近的观察点汇集。
适用这一目的的方法是一种层级提升的目的分类方法,该方法使用两个观察向量之间余弦函数的绝对值作为近似性指标,遵循最大约束量准则,即直径准则,以组成数据组(见以上引用的G.Saporta的著作)。分类适用于失控点的测量值,其中心在失控点数据的分布中心。
如果确定了与异常状况非常相关测量数据,异常方向则更容易解释。
为了达到这个目的,在没有“逐步”或“最佳子集”类型的常量项的情况进行线性回归,在变换和将控制分布中心作为解释变量之后把异常方向作为将要解释的变量和测量变量。
通过这样检查挑选的变量和所产生的回归系数的符号,过程负责人将被引导得出异常的根源。
A4-产生和检查特殊原因指示符的控制表
在前一步骤的结尾,得到了几组失控的观察点,每组都和特定的异常现象或特殊原因相关。
这一步首先把每种特殊原因和前面提到的首要的方向联系起来,然后,和两个指示符函数联系起来。这两个指示符的作用是:当应用到任何观察过程中,指示出观察是否能来自具有某种特殊原因的不正常运行过程。
图3和图4表示二维中与两个不同的异常现象相关的方向。图3表示用原始数值的数据观察空间的情况,图4是表示基于控制中的数据正常首要分量在基准结构中的情况。可以看出两种异常的分隔在图4中比在图3中更明显。
然后被考察的观察点位于基于控制下的数据的正常首要的分量基准结构中,它的原点和控制下数据的分布中心相吻合。
上面描述的理论模型引导我们去寻找代表某种异常的方向。它从失控观察的分布中心开始延展,靠向由于异常而产生的一组观察点。这个方向可以被定义为和上述异常相关的观察点的第一惯性轴,该轴通过原点(见前面所述的变换)。这个轴是由首要分量分析决定的第一主轴,其不集中,而且不被和异常相关的观察点团减弱。(见前面引用的G.Saporta的著作)
在评价了与特殊原因相关而偏离了受控点分布中心的观察点的方向之后,就有可能定义和异常相关的两个指示符了。
●第一个指示符测量观察点在上述方向上的距离;并且
●第二个指示符测量给定观察空间Rn中的一测定点到特定方向的靠近程度,用从失控点分布中心的角度临近方式。
在多维的观察空间中,第一个指示符沿着和特定原因相关的方向测定观察的距离。其将每个“观察和异常对”和观察和异常之间的关系的强度指示符相联系,称为“原因变量”,并作为观察点沿与原因相关的方向上的组成分量获得。此指示符作为观察向量和原因方向向量的标量乘积而计算出来。
通过在正规化的首要组份空间计算该标量乘积,其等于把最初的量度空间中的矩阵V-1作为标量乘积矩阵,在此当V为全秩时(fullrank),V-1代表受控数据的协差矩阵的逆,如果V不是全秩,V-1代表V的伪逆矩阵。
因而该指示符是初始测量值的线性组合,并且,假定这些测量值服从高斯定律(当过程在控制下时),这个指示符也服从高斯定律;其变化能用传统的控制图仿效。因而,在这样的强度图(intensity chart)上对失控的观察非常有可能和被记录在图上的指示符表示的异种类相对应。
但是,当观察点和控制下的点的分布中心相隔非常远的时候,它能在几个原因方向上显示出相隔很远,从而可能和几类异常情况相关。为了消除不确定性,考虑采用第二个指示符。
这是在代表异常的方向和以在控制下的数据的中心为原点并且指向观察点的向量之间的角度临近指示符。这个角度可以以其余弦求值,因为这样容易计算:通过用代表观察点的向量的范数除上述强度指示符来计算余弦。
因为这个计算是在赋范主分量空间执行的,它更适合,并且能被用来定义与这个新指示符相联的控制表上标出的控制限制。
在控制下运行的时候,该指示符具有决定分布函数F的概率分布。
函数F用来定义角度临近指示符Lα的控制极限。在此极限外对所述异常分配观察点,上述函数f代表在Rn上任意的给定方向和中心高斯随机向量Rn形成的角度(x)的余弦的分布函数,其中n≥2,并且随机高斯向量Rn具有单位矩阵作为协方差矩阵,有下列公式给出的随机变量的分布函数F:
如果x<-1 F(x)=0
如果x>1 F(x)=1
如果x≥-1且≤1
通过审查角度临近控制图,可能从被“强度”图指出为失控点的点中选出一些点来,这些点真正地和研究中的异常相依:这些点同时在两个图中被指出为异常点。
在实践中,这样的角度临近控制图易读性不强,因为在很少数量的维数中,控制极限非常接近1或-1。这就是为什么经过适当的变换后并且在有对数的标尺和用概率阀值表示的图表中,人们更愿意使用表示这些指示符的图表。
这种表示方法在所有检查到的量值的种类中都可以使用。它可以在展示和不同分布概率相关的图时更容易的解释和提供更大的一致性。下面定义了两种图:
●当控制用于单一的极限(象图T2的例子一样)时使用的单侧图Pu;以及
●当控制涉及下极限和上极限时使用的双侧图Pb。
和分布方程Fx的指标i相关的观察值x在单侧图中Pu中标出来。观察点分布在图中横坐标i和纵坐标y=min(log10(1-FX(x),4)的位置上。
纵坐标以与y的积分值相关的阀值做分度:1-Fx(x)=10(-y)。纵坐标上的4对应前面提到的“<0.0001”考虑到已经进行的截断。一条水平线标出选出的控制极限。图5表示常规的T2控制图和对应的Pu图。
和分布方程Fx的指标i相关的观察值x被标在双侧图Pb上。图上标出的观察值在横坐标i和纵坐标y上,定义如下:
如果x≤m y=min(-log10((1-FX(x)*2),4)
如果x>m y=min(-log10(FX(x)*2),4)
其中,m代表分布的平均值:FX(m)=0.5。
纵坐标以与整数值y相关的阀值做分度:10(-|y|)。水平线标出选定的控制阀值相关的控制极限。图6表示常规的余弦控制图和对应的Pb图。
A5-特殊原因的解释
前面定义的控制图提供给负责人一些附加信息,此信息使他们发现历史纪录中所有和正常过程的运行不符的点,并且对于每个点辨别异常的真正原因。上述工具使得每个识别出的异常的原因和观察空间中某一方向相关联,同时,两个指示符也使之可以识别。
●本发明方法的第二阶段涉及适当的控制,包括以下步骤:
S1-数据预先被变换
S2-通过对控制图的监视检测和识别异常
S3-必要时,组合原因或新的特别原因被识别并集成到本方法中。
在这第二步中,系统从过程接受测量和数据(直接从相连的传感器或手工输入得到),这些测量数据是基于分析过的历史测量数据的模型。当过程偏移时,能够引发报警。然后,从在学习过程中已经被辨别出来的原因中,系统指出偏移的可能原因。
从过程收到的每次观察按以下步骤来处理:
S1数据的预先变化
数据受到步骤A1定义的变换;
S2-检测和识别异常
与用在步骤A2中定义的计算模型做出的观察相关的T2的值对应T2图中的控制极限被计算并标出来。
如果控制图发现观察在控制之下,观察处理结束,否则,观察表示过程运行不正常。
在这种情况下,和学习阶段所识别出来的各种原因有关的指示符的取值被计算出来,代表新观察的点的位置将在和指示符有关的Pb型控制图中被检查。当和同一原因有关的两个指示符都指出观察失控,则这个原因被当作可能的过程误操作的起因。系统向过程控制者发出关于误操作的警报和误操作的起因的诊断。
观察过程将结束,除非系统又识别出来列入记入中的另一个的原因。在这种情况下,下面的附加步骤是必要的。
S3-必要时,原因的组合或一个新的特殊原因被识别和集成到系统中。
当在前一步骤中没有原因被识别出来,则要做检查看是否观察出的异常是源于两个已知原因同时出现。
为达到此目的,要不带常数项进行线性回归,其属于“逐步”(stepby step)类型或“最佳子集”(best subset)类型,把以受控分布中心为中心的观察点作为被解释变量,把原因变量作为解释变量。
当只使用两项解释变量时,如果获得了大量的相关系数R,并且,特别是R>0.95时,可以认为有两个和该两个变量相关的原因同时发生。也可以用类似方法考虑两个以上原因同时发生的可能性,尤其是三个原因或是四个原因。
当没有清单中的原因或原因组合被识别出来时,需要根据过程的负责人的知识和经验手工识别原因。
Claims (8)
1.一种用于控制工业过程的统计方法,该方法基于对输入、输出和控制的量值特征的指示符和测量数据的读取以及对所述过程的运行参数的读取,并且能够在将其值和其取样指标相关的基准框架中用观察点表示出来,所述指示符或传感器或测量数据可由手工输入给出,该方法的特征在于:
a)观察值被转化以使得结果值与多维的高斯分布模型一致,构成与本方法其余部分用的观察点对应的数据;
b)上述观察点位于多维空间,每个维度都和测量量值相关;
c)观察点中,那些与过程正常运行对应的所谓“控制下的”点和那些所谓“失控的”点有显著不同;
d)控制下的点的分布中心计算出来作为控制下的观察点的点集之重心;
e)可以识别出失控的观察点集中在“控制下的”点的分布中心的某一特定方向上;
f)这一方向与上述过程的通常的偏移原因相关;
g)每一个观察点和异常方向对与指示符相关,以便报告零个、一个或更多的所观察到的异常原因,其中所述异常原因可能与已进行的观察相关;并且
h)当一种异常情形用这种方式进行分析后,将触发警报,并且,在工业生产过程中用这种方式检测到的偏移将被补救。
2.根据权利要求1的统计过程控制方法,其特征在于:
每个“观察和异常对”和给出观察和异常之间的关系的强度的指示符相联系,被称为“原因变量”,指示符作为沿与原因相关的方向上的观察点的分量而获得,所述指示符作为观察向量和原因方向的方向向量相乘的标量积被计算,标量积矩阵被作为观察和异常之间的关系的强度指示符,即矩阵V-1,其中:当矩阵V是全秩时,V-1为受控数据矩阵V的协方差的逆矩阵,或者当矩阵V不是全秩,V-1为伪逆矩阵。
3.根据权利要求2的统计过程控制方法,其特征在于:在观察和异常之间使用了临近指示符,称作“角度临近指示符”,通过测量表示观察的向量和表示方向的向量之间的角度获得,角度临近指示符是所述角度的余弦,而此余弦值是通过用表示观察的向量的范数除密度指示符来计算的,所述范数通过所述矩阵V-1来定义。
5.根据权昨要求1到4中任何一个的统计过程控制方法,其特征在于将其组织成:
●在学习阶段,通过进行如下步骤,分析历史数据以分辨异常:
A1-数据的预先变换;
A2-划分控制下的测量数据和失控的测量数据;
A3-使用一种自动对观察数据分类的方法识别特殊原因类型和相关方向参数,该方法适用于识别和同一种异常情况相联系的观察数据,使得观察点靠近从控制下的点的分布中心引出的同一直线而成组;
A4-在特殊原因指示符上建立和检查控制图;并且
A5-解释特殊的原因
●控制阶段,该阶段对每个过程变量的新观察,判断观察是否与异常有关,并且如果这样,可能出现异常,该阶段包括如下步骤:
S1-数据的预先变换
S2-用检测控制图的方法探测和辨别异常;并且
S3-如果必要,将识别原因的组合或者新的原因并集成到方法中。
6.根据权利要求3到5的统计过程控制方法,其特征在于:
采用一种层级提升的自动分类方法,该方法使用两个观察向量之间余弦函数的绝对值作为近似性指标,遵循最大约束量准则,即直径准则,以组成数据组。
7.根据权利要求5的统计过程控制方法,其特征在于:
在学习阶段的步骤A3中用来代表异常的方向是和该异常相关的观察点的第一个惯性轴,该轴通过控制下的点的分布中心,该轴由不聚集且不减弱的和异常相关的观察点云的主要分量分析决定,即以前以控制下的观察点的分布中心为中心的点。
8.根据权利要求5到7中任何一个的统计过程控制方法,其特点在于:
控制图采用对数刻度表明极限阈值,在控制涉及单一的极限时(象图T2的例子一样)使用单侧图Pu,在控制涉及下限和上限时使用,双侧图Pb,该图提供了更容易解释的并且更统一的涉及不同分布概率的图形表示。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9812113A FR2783620B1 (fr) | 1998-09-22 | 1998-09-22 | Procede et systeme multidimensionnel de maitrise statistique des processus |
FR98/12113 | 1998-09-22 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1319212A true CN1319212A (zh) | 2001-10-24 |
CN100489870C CN100489870C (zh) | 2009-05-20 |
Family
ID=9530930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB998112364A Expired - Lifetime CN100489870C (zh) | 1998-09-22 | 1999-09-22 | 用于过程统计控制的方法和多维系统 |
Country Status (12)
Country | Link |
---|---|
US (1) | US6804563B1 (zh) |
EP (1) | EP1116150B1 (zh) |
JP (1) | JP4394286B2 (zh) |
CN (1) | CN100489870C (zh) |
AT (1) | ATE242513T1 (zh) |
AU (1) | AU5866999A (zh) |
CA (1) | CA2391716A1 (zh) |
DE (1) | DE69908620T2 (zh) |
ES (1) | ES2201771T3 (zh) |
FR (1) | FR2783620B1 (zh) |
HK (1) | HK1039668A1 (zh) |
WO (1) | WO2000017790A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251754B (zh) * | 2008-03-13 | 2010-06-02 | 西安交通大学 | 一种多工序加工过程误差流的处理及控制方法 |
CN102112933A (zh) * | 2008-10-09 | 2011-06-29 | 株式会社日立制作所 | 异常检测方法及系统 |
CN102246110A (zh) * | 2008-12-15 | 2011-11-16 | 斯奈克玛 | 用于监测航空发动机的数据标准化方法 |
CN102467699A (zh) * | 2010-10-29 | 2012-05-23 | 西门子公司 | 对缺陷管理建模及基于缺陷管理处理缺陷的方法 |
CN102542398A (zh) * | 2010-09-17 | 2012-07-04 | 波音公司 | 点云生成系统 |
CN103488135A (zh) * | 2013-08-14 | 2014-01-01 | 沈阳中科博微自动化技术有限公司 | 一种用于半导体生产加工过程监控的统计过程控制方法 |
CN103935120A (zh) * | 2014-03-25 | 2014-07-23 | 中山火炬职业技术学院 | 一种制定印刷企业质量控制参数标准的方法 |
CN112567302A (zh) * | 2018-08-13 | 2021-03-26 | 西门子股份公司 | 过程控制单元和用于在过程间交换过程变量的方法 |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL144358A (en) * | 2001-07-16 | 2006-10-31 | Oded Berkooz | Method for isolating sources of drifts in output properties for machines and processes |
WO2003009140A2 (en) * | 2001-07-20 | 2003-01-30 | Altaworks Corporation | System and method for adaptive threshold determination for performance metrics |
US7050933B2 (en) * | 2001-08-24 | 2006-05-23 | Bio-Rad Laboratories, Inc. | Biometric quality control process |
DE10151250A1 (de) * | 2001-10-17 | 2003-05-08 | Bayer Ag | Verfahren zur Ermittlung eines komplexen Korrelationsmusters aus Verfahrens- und Anlagendaten |
US7043403B1 (en) * | 2002-09-04 | 2006-05-09 | Advanced Micro Devices, Inc. | Fault detection and classification based on calculating distances between data points |
FI20030622A (fi) * | 2003-04-24 | 2004-10-25 | Tietoenator Oyj | Verkkopalveluihin liittyvien operaatioiden analysointi |
CN100419983C (zh) * | 2003-05-16 | 2008-09-17 | 东京毅力科创株式会社 | 处理系统健康指数及其使用方法 |
WO2005026893A2 (en) * | 2003-09-05 | 2005-03-24 | Sensitech Inc. | Automated generation of reports reflecting statistical analyses of supply chain processes |
US7369913B2 (en) * | 2004-04-02 | 2008-05-06 | Siemens Medical Solutions Usa, Inc. | Recipe editor and controller |
DE102004052302B4 (de) * | 2004-09-30 | 2013-05-23 | Immobiliengesellschaft Helmut Fischer Gmbh & Co. Kg | Verfahren zur Ausgabe von Messwerten und Anzeigeeinrichtung |
JP4399400B2 (ja) * | 2005-07-11 | 2010-01-13 | 株式会社日立ハイテクノロジーズ | 検査データ解析システムと検査データ解析プログラム |
KR100679721B1 (ko) * | 2005-11-01 | 2007-02-06 | (주)아이세미콘 | 반도체 공정장비의 변동 감지방법 |
ES2431863T5 (es) * | 2006-11-03 | 2017-07-27 | Air Products And Chemicals, Inc. | Sistema y método para la monitorización de procesos |
US7494893B1 (en) | 2007-01-17 | 2009-02-24 | Pdf Solutions, Inc. | Identifying yield-relevant process parameters in integrated circuit device fabrication processes |
US20080312875A1 (en) * | 2007-06-12 | 2008-12-18 | Yu Guanyuan M | Monitoring and control of integrated circuit device fabrication processes |
US8090676B2 (en) * | 2008-09-11 | 2012-01-03 | Honeywell International Inc. | Systems and methods for real time classification and performance monitoring of batch processes |
US8392009B2 (en) * | 2009-03-31 | 2013-03-05 | Taiwan Semiconductor Manufacturing Company, Ltd. | Advanced process control with novel sampling policy |
ES2424808B1 (es) * | 2011-01-25 | 2014-08-05 | Mondragón Goi Eskola Politeknikoa José María Arizmendiarrieta, S. Coop. | Método para detectar e identificar errores en procesos de fabricación |
CN103703083A (zh) | 2011-08-04 | 2014-04-02 | 亨斯迈先进材料(瑞士)有限公司 | 活性染料的混合物及其用途 |
US9070622B2 (en) * | 2013-09-13 | 2015-06-30 | Taiwan Semiconductor Manufacturing Company, Ltd. | Systems and methods for similarity-based semiconductor process control |
BR102014015604A2 (pt) | 2014-06-24 | 2015-03-31 | Gisele Castro Fontanella Pileggi | Máquina universal de ensaios reológicos e mecânicos |
US10089589B2 (en) * | 2015-01-30 | 2018-10-02 | Sap Se | Intelligent threshold editor |
ES2826348T3 (es) * | 2015-09-29 | 2021-05-18 | Siemens Ag | Procedimiento y sistema para el reconocimiento de fallos y la vigilancia de una parte de máquina regulada o controlada electrónicamente |
DE102016225081A1 (de) * | 2016-12-15 | 2018-06-21 | Robert Bosch Gmbh | Vorrichtung und Verfahren zum Bestimmen der Pinpoint-Fähigkeit möglicher Fehler einer oder mehrerer Komponenten |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0660826B2 (ja) * | 1989-02-07 | 1994-08-10 | 動力炉・核燃料開発事業団 | プラントの異常診断方法 |
US5150289A (en) * | 1990-07-30 | 1992-09-22 | The Foxboro Company | Method and apparatus for process control |
US5442562A (en) * | 1993-12-10 | 1995-08-15 | Eastman Kodak Company | Method of controlling a manufacturing process using multivariate analysis |
US6539267B1 (en) * | 1996-03-28 | 2003-03-25 | Rosemount Inc. | Device in a process system for determining statistical parameter |
US6442445B1 (en) * | 1999-03-19 | 2002-08-27 | International Business Machines Corporation, | User configurable multivariate time series reduction tool control method |
US6424876B1 (en) * | 1999-07-22 | 2002-07-23 | Advanced Micro Devices, Inc. | Statistical process control system with normalized control charting |
-
1998
- 1998-09-22 FR FR9812113A patent/FR2783620B1/fr not_active Expired - Fee Related
-
1999
- 1999-09-22 US US09/787,685 patent/US6804563B1/en not_active Expired - Fee Related
- 1999-09-22 WO PCT/FR1999/002254 patent/WO2000017790A1/fr active IP Right Grant
- 1999-09-22 CN CNB998112364A patent/CN100489870C/zh not_active Expired - Lifetime
- 1999-09-22 EP EP99946223A patent/EP1116150B1/fr not_active Expired - Lifetime
- 1999-09-22 JP JP2000571380A patent/JP4394286B2/ja not_active Expired - Fee Related
- 1999-09-22 AT AT99946223T patent/ATE242513T1/de not_active IP Right Cessation
- 1999-09-22 CA CA002391716A patent/CA2391716A1/fr not_active Abandoned
- 1999-09-22 ES ES99946223T patent/ES2201771T3/es not_active Expired - Lifetime
- 1999-09-22 DE DE69908620T patent/DE69908620T2/de not_active Expired - Lifetime
- 1999-09-22 AU AU58669/99A patent/AU5866999A/en not_active Abandoned
-
2002
- 2002-02-21 HK HK02101255.6A patent/HK1039668A1/zh unknown
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251754B (zh) * | 2008-03-13 | 2010-06-02 | 西安交通大学 | 一种多工序加工过程误差流的处理及控制方法 |
CN102112933A (zh) * | 2008-10-09 | 2011-06-29 | 株式会社日立制作所 | 异常检测方法及系统 |
CN102112933B (zh) * | 2008-10-09 | 2014-06-18 | 株式会社日立制作所 | 异常检测方法及系统 |
CN102246110B (zh) * | 2008-12-15 | 2014-12-03 | 斯奈克玛 | 用于监测航空发动机的数据标准化方法和系统 |
CN102246110A (zh) * | 2008-12-15 | 2011-11-16 | 斯奈克玛 | 用于监测航空发动机的数据标准化方法 |
CN102542398A (zh) * | 2010-09-17 | 2012-07-04 | 波音公司 | 点云生成系统 |
CN102542398B (zh) * | 2010-09-17 | 2016-04-13 | 波音公司 | 点云生成系统 |
CN102467699A (zh) * | 2010-10-29 | 2012-05-23 | 西门子公司 | 对缺陷管理建模及基于缺陷管理处理缺陷的方法 |
CN102467699B (zh) * | 2010-10-29 | 2016-08-03 | 西门子公司 | 对缺陷管理建模及基于缺陷管理处理缺陷的方法 |
CN103488135B (zh) * | 2013-08-14 | 2015-11-04 | 沈阳中科博微自动化技术有限公司 | 一种用于半导体生产加工过程监控的统计过程控制方法 |
CN103488135A (zh) * | 2013-08-14 | 2014-01-01 | 沈阳中科博微自动化技术有限公司 | 一种用于半导体生产加工过程监控的统计过程控制方法 |
CN103935120A (zh) * | 2014-03-25 | 2014-07-23 | 中山火炬职业技术学院 | 一种制定印刷企业质量控制参数标准的方法 |
CN112567302A (zh) * | 2018-08-13 | 2021-03-26 | 西门子股份公司 | 过程控制单元和用于在过程间交换过程变量的方法 |
Also Published As
Publication number | Publication date |
---|---|
AU5866999A (en) | 2000-04-10 |
WO2000017790A1 (fr) | 2000-03-30 |
CN100489870C (zh) | 2009-05-20 |
ES2201771T3 (es) | 2004-03-16 |
FR2783620B1 (fr) | 2002-03-29 |
JP2002525757A (ja) | 2002-08-13 |
EP1116150B1 (fr) | 2003-06-04 |
DE69908620T2 (de) | 2004-04-29 |
ATE242513T1 (de) | 2003-06-15 |
US6804563B1 (en) | 2004-10-12 |
EP1116150A1 (fr) | 2001-07-18 |
JP4394286B2 (ja) | 2010-01-06 |
DE69908620D1 (de) | 2003-07-10 |
HK1039668A1 (zh) | 2002-05-03 |
CA2391716A1 (fr) | 2000-03-30 |
FR2783620A1 (fr) | 2000-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1319212A (zh) | 用于统计过程控制的多维方法和系统 | |
CN1165869C (zh) | 通过均衡正交映射可视化和自组织多维数据 | |
US7079979B2 (en) | Inspection method, inspection apparatus, and facility diagnosis unit | |
CN1758034A (zh) | 知识作成装置及参数搜索方法 | |
JP4394728B2 (ja) | 影響要因特定装置 | |
CN1214340C (zh) | 多个神经网络的成像设备和方法 | |
CN1488933A (zh) | 良否判定装置、判定程序、方法以及多变量统计解析装置 | |
JP4911055B2 (ja) | バッチプロセスデータの解析装置およびそれを用いた異常検出/品質推定装置 | |
CN1945602A (zh) | 一种基于人工神经网络的特征选择方法 | |
CN101069184A (zh) | 信息处理装置和方法、及程序 | |
CN1929530A (zh) | 图像质量预测方法和设备以及故障诊断系统 | |
CN1256709C (zh) | 纸券类的真假判断方法 | |
CN1615434A (zh) | 适用于处理电化学信号的方法和装置 | |
CN1912891A (zh) | 人脸元数据生成设备和生成方法、人脸识别系统和方法 | |
CN1834607A (zh) | 检查装置及检查方法 | |
CN1924734A (zh) | 一种在线质量检测的控制方法 | |
CN100351610C (zh) | 状态判断方法 | |
JP4764490B2 (ja) | ハードウェア使用状況に応じたユーザー評価装置 | |
JP2008008708A (ja) | 判定知識作成装置、判定知識作成方法、プログラムおよび記録媒体 | |
CN1701430A (zh) | 用于数据分析的专家知识方法和系统 | |
Tello et al. | Automatic flower number evaluation in grapevine inflorescences using RGB images | |
CN1653492A (zh) | 纸片鉴别装置、纸片鉴别方法及程序 | |
JP2005092466A (ja) | 診断プロセス支援方法とそのためのプログラム | |
CN1756938A (zh) | 特征判定方法和特征预测方法以及装置 | |
CN1811802A (zh) | 质量变动显示装置、显示方法、显示程序及记录介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: WD Ref document number: 1039668 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CX01 | Expiry of patent term | ||
CX01 | Expiry of patent term |
Granted publication date: 20090520 |