指数加权滑动平均控制图与自回归模型的实例应用.doc

1、1指数加权滑动平均控制图与自回归模型的实例应用摘要：统计过程中的控制图如今已经成为质量管理中一种不可缺少的工具，已经成为一种控制的潮流和方法，它的中心思想就是预防为主，第一次就把事情做对，消除一些变异的因素和区分两种不同性质的波动。经过对大量文献的认真研究和学习，关于统计过程中的控制图，所提出的是假定过程均值不变等条件下，通过轴承的半径分析，用指数加权平均的控制方法来处理相关问题。目前为止还没有将其与常规控制图整合在一起，因此通过其与回归模型的方法相建立，来解决实际中的问题。 关键词：指数加权滑动平均（EWMA） ；控制图；自回归模型 中图分类号：F0 文献标志码：A 文章编号：1673-29

2、1X（2014）16-0010-03 1924 年，休哈特在美国贝尔实验室第一次发明质量控制图的概念，经过数年的发展已经逐渐成为质量工作管理中不可缺少的工具之一，过程控制是为了贯彻预防为主，在统计过程控制中，利用休哈特常规控制图时，我们总是假定过程均值不变，同时随机误差独立同分布，服从均值为 0、方差为 2 的正态分布。但在实际应用中，这些前提假设是很难实现的，无法获得相应的数据，许多过程中都存在自相关的现象，独立性也很难成立。这就使得常规控制图方法难以继续使用。这时可以用指数加权滑动平均（EWMA）控制图方法处理相关的问题。可以用时间序列2对于自相关过程的观测值进行适当的拟合，得出模型，再用

3、其预测值、残差来构造控制图。本文通过具体数据分别做指数加权滑动平均、自回归模型拟合及常规控制图，通过实际得到的图形简单地做了比较。 一、指数加权滑动平均（EWMA）控制图 EWMA 控制图上的点包含了所有前面子组的信息，它能探测过程的微小漂移。EWMA 控制图同样适用于单个观测值的情况。由前所诉，EWMA 控制图不受过程均值不变，同时随机误差独立同分布，服从均值为 0、方差为 2 的正态分布等条件限制，常用于处理序列数据。在 EWMA 控制图中绘制的统计量为当前值与历史数据的加权平均，即 Zt=xt+（1-）Zt-1 对于比较大的 t，当xt独立同分布，服从均值为 、方差为 2的正态分布时，样

4、本统计 Zt 量近似服从正态分布，且 E（Zt）=，Var（Zt）=2 1-（1-）2t 随着 t 的增大，1-（1-）2t趋近于 1。 绘制 EWMA 控制图时，要求给定权重 。 （1-）表示 EWMA 对历史测量值的权重，若（1-）越大（ 越小）表示历史对于现在的影响越重。权重 决定了 EWMA 对于当前数据突然发生变化的能力，具有时效性。 越小，时效性不强，则 Zt 的时间序列图更加平稳。EWMA 控制图的优点就是，当较小值或较大值进入计算时不会严重影响到这些控制图。对于权重 的选取，Montgomery（1991）推荐使用 0.050.25。一般 常取 0.08，0.10，0.15，更

5、要根据实际情况而定。 3下面给出具体例子。数据来源（见参考文献1中例 3.6.4） 。此例是单个观测值的情况，给出了不同时间测得的 50 个凸轮轴的轴承直径。给定权重 值为 0.2，用指数加权滑动平均的公式，具体得到 Zt 的值（初始值 Z0=x0） 。 Zt=0.2xt+（1-0.2）Zt-1 将目标值 50 作为 的估计（此处是以起点时间初始值 50 作为 的估计） ，对上述数据做 EWMA 控制图（如下页图 1 所示） 。 从该图可以看到控制限在第 11 个样本开始稳定。从 EWMA 控制图上可看出，样本点对目标值 50 的漂移。 二、自回归模型（auto regression mode

6、l，AR） 自回归（AR（n） ）模型是时间序列 ARMA 中的一种特例，即ARMA（n，0）只有自回归部分，没有移动平均部分。其实 ARMA 模型才是更一般的平稳的时间序列模型，只是由于文中例子的数据模拟后得到模型，便以此说明比较。对于后面所说得到 AR（n）模型，求出预测值和残差均是基于 AR（n）模型。对于一般的 ARMA 模型，同理亦可操作。 对于时间序列 AR（n）模型具有如下一般结构： xt=0+1xt-1+2xt-2+nxt-n+t E（t）=0，Var（t）=2 ，E（ts）=0，st E（txs）=0，st 对于本文数据，作出轴承直径数据的自相关图与偏自相关图分别（见下页图

7、2、图 3） 。 4可以看出大致的自相关图是拖尾的，偏自相关图时 1 阶截尾的，由ARMA 模型的判别条件符合模型 AR（1） 。拟合 AR（1）模型： xt=0+1xt-1+t 通过显著性检验，延迟各阶的 LB 统计量的 P 值均显著大于 0.05，残差通过白噪声检验，所以该拟合模型显著成立。得到 AR（1）模型为： （xt-48.55）=0.74（xt-1-48.55）+t xt=12.62+0.74xt-1+t 以此模型做预测值，用测量值减去预测值得到残差 r。下图给出观测值（） 、预测值（*） 、95%上下限（）： 三、基于 AR 模型所得残差做指数加权滑动平均（EWMA）控制图 由于

8、时间序列的预测值和观测值之间的残差可以极大地弱化过程的自相关性，并可以近似认为各个时刻的残差都是相互独立的，因此可以采用一般的各种控制图方法进行监控。在此，用上述所得的残差做 EWMA控制图（见图 5） ， （取目标值为 1.444 是对应于初始值 50 的残差值） 四、轴承直径的 EWMA 控制图对比其残差的 EWMA 控制图 对比轴承直径的 EWMA 控制图（见上页图 1）和残差 EWMA 控制图（见图 4） ，两图的曲线波动轮廓一致，但是轴承直径的 EWMA 控制图漂移幅度幅度远大于残差 EWMA 控制图。残差的 EWMA 控制图对于目标值的漂移超出控制限的值的个数少于轴承直径的 EWM

9、A 控制图。可见轴承直径的 EWMA控制图要比残差 EWMA 控制图更灵敏。但是轴承直径在本例中是自相关的，而残差是独立的。对于存在自相关，超出控制限的值对于其附近时间点值的的影响较大，所以连带一片的值都会超出控制限。 5五、总结 EWMA 控制图，不受正态假定的限定，图上的每个点包含着前面所有子组的信息，具有检出过程目标值小漂移的敏感性。对于时间序列来说，除了 EWMA 控制图，本文还用的拟合出时间序列模型，求出预测值、残差，再用残差作 EWMA 控制图。 参考文献： 1周纪芗，茆诗松.质量管理统计方法(第 2 版)M.北京：中国统计出版社，2008. 2王燕.应用时间序列分析(第 3 版)M.北京：中国人民大学出版社，2012. 3王海宇，等.平稳自相关过程的 EWMA 控制图J.工业工程，2006，(3). 4崔敬巍，谢里阳.自相关过程的统计控制方法研究J.机械制造，2006，(9). 责任编辑 吴高君 收稿日期：2014-03-08 作者简介：魏关军（1976-） ，男，河南潢川人，质量经理，硕士研究生，从事质量管理研究。