基于APL的双抽样均值控制图的优化设计.doc

1、1基于 APL 的双抽样均值控制图的优化设计摘要：以平均产品长度作为度量标准，研究了双抽样均值控制图的优化设计问题。分析了双抽样均值控制图的平均产品长度的计算方法，在此基础上建立了控制图参数设计的优化模型，并通过与原有的双抽样均值控制图的性能对比，验证了该优化设计模型的有效性。 关键词：平均产品长度 双抽样控制图 优化设计 1 概述 控制图是企业中进行质量管理和控制的基本方法之一，对产品质量的好坏起着至关重要的作用。控制图首先由 Shewhart 于 1924 年提出，之后得到了广泛的研究与引用，取得了相当不错的经济效益与社会效益，尤其是对战后世界经济的发展做出了巨大的贡献。但随着制造业的发展

2、，生产的灵活性越来越高，原先的固定抽样区间和固定样本容量的静态控制图已不能满足质量监控的需要，于是 Reynolds 等提出了变化抽样区间的均值（Variable Sampling Intervals，VSI）控制图1，并由此形成了动态控制图这一新的研究领域。后来又有人提出变化样本容量（Variable Sampling Size，VSS）控制图、双抽样（Double Sampling，DS）控制图等，并取得了一系列关于这类动态控制图的研究成果2-4。Costa 曾比较了 VSS、VSI 和 DS 控制图，结果发现 DS 控制图的监测效率最高5。目前，DS 控制图正逐渐成为研究的热点。DS 控

3、制图最早是由 Daudin 提出的，他首先将二次抽样检测的概念结合到控制图2理论中，采用两个阶段的抽样来监控过程，这种方法能降低抽样数量并提升检测过程偏移能力6。He 和 Grigoryan 又提出了用于求解 DS 控制图各个参数的遗传算法7。上述 DS 控制图的研究，都是基于传统的效率评价方法，平均运行长度（Average Run Length，ARL） ，即平均需要多少个样本来发现过程异常，但由于 DS 控制图采用两次抽样，每次抽样的样本数量是不一样的，因此即使在 ARL 相同的情况下，其检测异常波动的效果仍有很大的差异。为此，本文将采用平均产品长度（Average Product Len

4、gth， APL）作为控制图的效率度量指标，来进行 DS 控制图的优化设计，并与其他的一些动态控制图方法进行性能比较。 2 DS 控制图 Daudin 提出的 DS 控制图将两个不同控制限的常规控制图结合在一起进行过程质量监控，并在双抽样的第一阶段的控制图中增加了警告界限，如图 1 所示。 第一阶段抽样及监控 第二阶段抽样及监控 图 1 DS 均值控制图 图 1 中，L1 和-L1 分别表示第一阶段的控制界限；w 和-w 分别表示第一阶段的警告界限；L2 和-L2 分别表示第二阶段的控制界限。在 DS 控制图的使用过程中，首先从第一阶段抽样和监控开始，先抽取样本容量为 n1 的一个样本，并计算

5、均值 X1，若 X1 落在 I1 区域，则判断过程处于受控状态；若 X1 落在 I2 区域，则判断过程处于失控状态；若 X1 落在 I3区域，则无法直接进行判断，需要进行第二阶段抽样及监控，再抽取样本容量为 n2 的一个样本，并计算均值 X2，在第二阶段监控时必须包含第3一和第二阶段的样本数据，因此第二阶段控制图的统计量为 X2=，若X2落在 I4 区域，则判断过程处于受控状态；若落在 I5 区域，则判断过程处于失控状态。 3 APL 的计算 平均产品长度的定义为从均值或方差的变化发生到这个变化在控制图中被发现之间生产工序平均生产的产品数量8。而受控情况下的APL，用 APL0 表示，指从方案

6、开始到产生一个错误报警信号之间平均生产的产品数。n 表示控制图上描一个点代表的样本容量的大小。S 表示在偏移发生到被发现之间抽样的样本的个数。Z 表示从偏移发生到之后第一个样本之间的产品数。一个常量 h 表示样本之间间隔的产品数。L，即产品运行长度，表示产品的均值或方差的变动发生到该变动被发现之间生产的总产品数。如图 2 所示，L 的数学表达可以用公式（1）表示： L=Z+h（S-1）+nS=Z+hS-h+nS； （1） 产品运行长度 L 的期望，即 APL，就可用样本容量 n 和抽样比例 r 表示出来： APL=E（L）=E（Z）+E（n）+hE（S）-h； （2） 对于 E（S） ，它就是

7、通常用于控制图性能度量的平均运行长度（ARL） ，若过程处于稳定状态，通常用 ARL0=1/a 表示，其中 a 为控制图第一类错误（即虚发报警）的概率；若过程出现异常，导致过程均值偏移了？啄倍标准差？啄 ，此时用 ARL？啄=，？啄0 表示，其中 ？啄为控制图第二类错误（漏发报警）的概率。DS 控制图的 ARL0 和 ARL？啄的计算见式（3）和（4） 。 4偏移产生 信号产生 图 2 运行长度的计算图解 ARL0=； = （3） ARL？啄=； 其中，Z1=；Z2=；Z=； C=；r=；I*=（-w+？啄，w+？啄） ； I*=（-L1+？啄，-w+？啄）（w+？啄，L1+？啄） 。P（?）

8、表示落在某区域的概率，？椎（?）表示标准正态分布的累计分布函数，？准（?）表示标准正态分布的概率密度函数。 对于 E（n） ，在常规控制图中，样本容量是一个常量，每次描点采用的样本容量都是一样的，但对于 DS 控制图来说，每次描点可能只需用进行一重抽样，也可能需用进行二重抽样，因此，每此描点需用的样本数是不固定的，在过程处于稳定状态下，需用的样本量是比较小的，而当过程出现一定程度的偏移时，就常常需用进行二重抽样，需用的样本量就会比较大，因此其数学期望可以分别考虑在过程稳定时的期望样本量E0（n）和过程偏移程度为？啄倍标准差？啄 的期望样本量 E？啄（n） ，见式（5）和（6） 。 E0（n）=

9、n1+n2P（Z1I3）=n1+n2?？椎（-w）-？椎（-L1）+？椎（L1）-？椎（W） E？啄（n）=n1+n2P（Z1I3）=n1+n2?？椎（-w+？啄）-？椎（-L1+？啄）+？椎（L1+？啄）-？椎（W+？啄） （6） 5对于 E（Z） ，当 n 的值为定值时，Z 服从均匀分布，E（Z）=，但由于 DS 控制图每次描点需用的样本数是不固定的，因此同样分别考虑过程稳定和过程出现偏移情况下两种状态的数学期望： E0（Z）= E？啄（Z）= 将式（3）（7）代入到式（2） ，就可以计算出 APL0 和 APL？啄：APL0=E0（Z）+E0（n）+h?ARL0-h APL？啄=E？啄（

10、Z）+E？啄（n）+h?ARL？啄-h 4 优化设计模型 DS 控制图的设计主要是确定 5 个参数：两个控制限参数（L1 和 L2） 、两个阶段的抽样数（n1 和 n2）以及第一阶段的警告限参数 w。控制图的一般设计原则是：当过程稳定，希望错误报警的概率要小，即 APL0 尽量的大；当过程出现异常，希望漏发报警的概率要小，即 APL？啄越小越好。但虚发和漏发报警的概率为反比关系，不能同时达到最优。因此，可以先令 APL0 取一个足够大的可以接受的固定值，在此基础上要求 APL？啄达到最小。于是，设 APL0 为一较大的常数 C，抽样间隔的产品数为常数h，则优化设计模型就可以表示为： Min A

11、PL？啄 约束：APL0=CAPL？啄 在上式中，两个阶段的抽样数 n1 和 n2 为整数，为了减少计算量，取小于 50 的数，事实上生产实践中也不可能取过大的样本容量。同样，两个控制限参数 L1 和 L2 以及第一阶段的警告6限参数 w 都取小于 6 的数。表 1 为分别设定 C=10000，h=100 和C=50000，h=100 两种情况下通过上式优化设计模型，运用 MATLAB 进行计算的结果。 表 1 在给定 APL0 和抽样间隔的产品数 h 时 DS 控制图的优化设计结果 5 性能对比 在我们的模型中，采用的控制参数都是根据 APL 选择最优值，那么这时的控制方案到底比以往的 DS

12、 控制图有多优越呢？现在我们就拿我们的模型与 Irianto 和 Shinozaki 提出的 DS 图设计方案9进行比较。分别用基于 APL 的优化设计模型（Optimal Design）和 Irianto 和Shinozaki 提出的设计方案（I&S）针对 1 倍和两倍标准差的过程偏移进行设计确定控制参数，再用两种方案对各种程度的过程偏移进行监控效率的对比，结果见表 2。 由于 Irianto 和 Shinozaki 提出的设计方案（I&S）是以 为监控效率的评价指标，一般都取 ARL0 在 370 左右，按照抽样样本间隔为 h=100的假设，换算出 APL0 约为 40000，为了便于对两

13、种设计方案进行比较，因此同样取 C=40000。表 2 中给出了各种程度过程异常偏移情况下的 APL、ARL和 E（n）的值，从表中的结果可以发现，不论是对较小的过程偏移？啄=1.0 还是对较大的过程偏移？啄=2.0，本文提出的优化设计方案发现异常偏移的效率都是高于 Irianto 和 Shinozaki 提出的设计方案的。 76 结论 本文主要讨论了双抽样均值控制图的平均产品长度的计算方法，并以此作为效率评价指标，建立了该控制图的优化设计模型。通过与原有的双抽样均值控制图设计方案的对比，说明新的优化设计模型能够有效提高双抽样均值控制图技术检测各种程度异常偏移的有效性。 参考文献： 1Reyn

14、olds M R，Amin R W，Arnold J C and Nachlas，J.A.X Charts with Variable Sampling IntervalsJ. Technometrics.1988，30：181-192. 2Costa A F.X Charts with Variable ParametersJ. Journal of Quality Technology.1999，31：408-416. 3Costa A F.Joint X and R Charts with Variable Sample Size and Sampling IntervalsJ.Jour

15、nal of Quality Technology，1999，33：387-397. 4王兆军.关于动态质量控制图的设计理论J.应用概率统计. 2002，18（3）：316-332. 5Costa A F.，Rahim M A.Joint X and R Charts with Two Stage Sampling C.Proceedings of IIE Research Conference.Orlando，Fla.，USA：IIE Press，2002：22-24. 6Daudin J J.Double Sampling X-bar Charts J.Journal of Quality

16、 Technology.1992，24（2）：78-87. 7He D.Grigoryan A.Construction of Double Sampling Control 8Charts for Agile ManufacturingJ.Quality and Reliability Engineering International，2002，18（4）：343-355. 8王海宇.基于平均产品长度的 EWMA 图性能评价与优化J. 管理工程学报，2007，21（3）：164-167. 9Irianto D.and Shinozaki N.An Optimal Double Sampling X-bar Control ChartJ.International Journal of Industrial Engineering.1998，5：226-243. 10Hsu，L.F.Note on“Design of Double and Triple Sampling X-bar Control Charts Using Genetic Algorithms”J.International Journal of Product Research.2004，42（5）：1043-1047. 作者简介： 娄琳华（1979-） ，女，河南新乡人，中级经济师，研究方向：质量管理。