未知協(xié)方差條件下的VSSI T2控制圖的經(jīng)濟性設(shè)計

顧麗君，唐慶國

（南京理工大學(xué) 經(jīng)濟管理學(xué)院，南京 210094）

0 引言

自從1947年，Hotelling提出了第一張均值向量T2控制圖之后，多元統(tǒng)計質(zhì)量控制便進入了人們的視野。現(xiàn)代產(chǎn)品指標(biāo)多樣，且關(guān)系復(fù)雜，多變量控制圖更加符合生產(chǎn)的需求。

與單變量控制圖相同，大多數(shù)多元均值向量T2控制圖是在協(xié)方差已知且保持不變的前提之下討論對均值的控制。但在實際生產(chǎn)中，協(xié)方差并不是一成不變的。一些學(xué)者進而研究了多變量協(xié)方差控制圖[1-3]，如 ||S圖、W圖、L圖等。但是這些控制圖存在明顯的缺陷，即要求“協(xié)方差已知”，這在實際生產(chǎn)控制中是十分苛刻的。

本文結(jié)合經(jīng)典經(jīng)濟模型研究了未知協(xié)方差矩陣的多元T2控制圖的經(jīng)濟性能。利用數(shù)值實例，對其靈敏度進行分析，為企業(yè)參數(shù)的選擇提供依據(jù)和建議。

1 總體協(xié)方差未知的VSSIT2控制圖

假設(shè)產(chǎn)品有p個質(zhì)量特征，產(chǎn)品質(zhì)量特征X服從均值向量為μ，協(xié)方差矩陣為Σ的p維正態(tài)分布。Xij'i=1'…'n(i)'j=1'...'m為來自X的第i個樣本子組，m表示子組容量，Xˉi表示子組 i的均值向量，即假設(shè)過程中只受到一個異常原因的影響，使得過程均值從μ0偏移到 μ1，μ0是目標(biāo)值，則T2控制圖統(tǒng)計量為：當(dāng)總體的協(xié)方差矩陣 Σ 未知時，需要利用已有的樣本數(shù)據(jù)估計Σ。式中的μ0和Σ需要用各子組均值向量的平均向量Xˉ和各子組協(xié)方差矩陣的平均值SP代替其中Si為子組i的協(xié)方差矩陣。假設(shè)用來建立控制圖的n個子組均為受控的子組，則和SP均為 μ0和Σ的無偏估計。構(gòu)造T2控制圖統(tǒng)計量為-Xˉ)。

此時服從第一自由度為p，第二自由度為(n-p)的F分布[4,5]。即

對于VSSIT2控制圖，假設(shè)抽樣間隔為(h1'h2)，h1＞h2；抽樣大小為 (n1'n2)，n2＞n1；警戒限為ω，大小不變，k為控制限。于是控制圖可分為以下三部分：I1=[0'ω]；I2=[ω'k]；I3=[k'+∞)。當(dāng)取樣落在受控區(qū)域時，則下一次取樣的參數(shù)選取如下：

當(dāng)則判定過程失控。

2 衡量指標(biāo)

通常衡量控制圖性能的指標(biāo)主要有三個：ATS（平均運行時間）；ANOS（整個周期的平均檢測數(shù)）；ANF（平均誤報次數(shù)）。未知協(xié)方差控制圖的研究目前僅局限于概率分布研究，并沒有考慮其經(jīng)濟性設(shè)計，也沒有將其與動態(tài)控制圖結(jié)合起來。本文用馬爾科夫鏈預(yù)測法研究VSSI控制圖的統(tǒng)計性指標(biāo)，同時對其經(jīng)濟性進行了設(shè)計和分析。

設(shè)受控狀態(tài)為I，失控狀態(tài)為N，本文一共構(gòu)建了6個馬爾科夫鏈狀態(tài)空間

假設(shè)異常原因發(fā)生的時間服從參數(shù)為λ的指數(shù)分布，那么異常原因發(fā)生的概率為1-e-λhi'i=1'2。當(dāng)過程處于狀態(tài){I 3 'I}時，屬于誤報警；當(dāng)過程處于狀態(tài){I3'N }時，即為吸附態(tài)。則該過程的轉(zhuǎn)移概率矩陣為：

其中Pij表示從狀態(tài)i轉(zhuǎn)移到狀態(tài)j的概率：

對于評判控制圖的指標(biāo)，計算公式如下：

其中，h為各個狀態(tài)對應(yīng)的樣本間隔，r為過程的初始概率，f為各個狀態(tài)對應(yīng)的誤報警概率，n為各個狀態(tài)對應(yīng)的樣本大小。為避免生產(chǎn)最初就存在問題，本文中取（r'為r的轉(zhuǎn)置，其他同理）：

3 數(shù)值實例

每個生產(chǎn)周期都包含了受控階段、失控階段、檢測階段和維修階段。關(guān)于費用函數(shù)，已有不少的經(jīng)典模型。本文中的費用函數(shù)引用Costa經(jīng)濟性設(shè)計模型[6]。

式（5）中V0表示受控狀態(tài)下每小時的利潤；V1表示失控狀態(tài)下每小時的利潤；C0表示單個抽樣的成本；C1表示誤報警時的平均搜索費用；A1表示搜索糾正異常原因的平均費用；T0表示搜索誤報警的平均時間；T1表示搜索糾正異常原因的平均時間。

結(jié)合數(shù)值實例研究VSSIT2控制圖的統(tǒng)計性和經(jīng)濟性，本文取p=2，式（5）中的參數(shù)選取為：λ=0.01,T0=1.5h,T1=1.5h,C0=3.6元/個 ，C1=115元/次 ，A1=465元/次 ，V0=135元/次，V1=60元/次。控制圖的設(shè)計參數(shù)如表1所示。

表1 控制圖的設(shè)計參數(shù)

當(dāng)n=4時，～3F(2'2)，即1/3服從第一自由度和第二自由度均為2的F分布；當(dāng)n=6時，0.4服從第一自由度和第二自由度分別為2和4的F分布。基于以上參數(shù)，可以推導(dǎo)出矩陣P中的元素。例如：

其中F(x'p'n-p'η)表示非中心F分布的累計概率分布函數(shù)，p為第一自由度，n-p為第二自由度，η為非中心參數(shù)。 η0=0，ηi=niδ2，δ為中心從 μ0偏移到 μ1的馬氏距離。a1=a4=1/3，a2=a3=a5=0.4。經(jīng)計算，VSSIT2控制圖的統(tǒng)計性和經(jīng)濟性指標(biāo)如表2所示。

表2 VSSIT2控制圖的統(tǒng)計性經(jīng)濟性指標(biāo)

由表2可以看出，其統(tǒng)計性和經(jīng)濟性指標(biāo)基本符合該控制圖的一般規(guī)律。隨著δ2的不斷增大，AATS、ANOS和EL的值都不斷減小，且減小的頻率越來越小。即過程偏移較大時，VSSIT2控制圖發(fā)現(xiàn)異常的速度較快，且成本較小。

假定過程的偏移為δ2=0.5，運用Matlab軟件中的optimtool的遺傳算法工具箱來搜索控制圖經(jīng)濟性設(shè)計的最優(yōu)解，以便相應(yīng)的單位時間損失EL達到最小值。考慮到實際情況和計算簡便，求解下面的最優(yōu)化問題：

minE(L)(n1'n2'h1'h2'w)

約束：0＜n1＜n2＜20'ni為整數(shù)

20＞h1＞h2＞0

0＜ω＜k

得到控制圖的最優(yōu)參數(shù)及其EL值如表3所示。

表3 VSSIT2控制圖的優(yōu)化參數(shù)

4 靈敏度分析

為了研究其他參數(shù)（如 δ、V0、V1）對EL值的影響，本文結(jié)合因子分析方法，三個參數(shù)各取兩水平，即作三因子兩水平23的全因子試驗設(shè)計。參數(shù)的兩水平選取如表4所示，該試驗設(shè)計的運行順序和該條件下優(yōu)化的結(jié)果如表5所示。

表4 三個參數(shù)的兩個水平

表5 全因子試驗設(shè)計及優(yōu)化結(jié)果

對于每個試驗子組，都用Matlab里的遺傳算法工具箱求得最優(yōu)解。根據(jù)上表數(shù)據(jù)，可研究表中的參數(shù)對VSSI T2控制圖經(jīng)濟性設(shè)計的影響作用。

本文采用Minitab軟件針對上述三個參數(shù)作方差分析，得到主效應(yīng)圖（見圖1）、效應(yīng)排列圖（見圖2）和方差分析結(jié)果（見表6、表7）。

圖1 三個參數(shù)對EL值的主效應(yīng)圖

由上述圖表可得以下結(jié)論：V1越大，V0越小；δ2越大，EL值越小，生產(chǎn)過程中造成的利潤所失越小。A、B、C三個因子（V1、V0、δ2）對EL均有顯著影響，另外 δ2和V0的交互作用對其也有顯著影響。從方差分析可以看出δ2、V1、V0的P-Value均小于0.05，可見其確實顯著影響了利潤所失的大小。結(jié)合實際生產(chǎn)情況，可以得出：在不影響正常生產(chǎn)的情況下，受控時利潤和失控時利潤接近，VSSIT2控制圖檢測大偏移時，利潤所失更加小。

5 總結(jié)

多元T2控制圖被廣泛應(yīng)用于監(jiān)測多元質(zhì)量問題，而VSSI控制圖已經(jīng)被證明在檢測微小偏移時性能明顯優(yōu)于一般控制圖。本文研究了未知協(xié)方差的VSSIT2控制圖，基于已有的研究成果，結(jié)合馬爾科夫鏈方法，推導(dǎo)出了轉(zhuǎn)移概率。通過數(shù)值實例，計算了在不同偏移下控制圖的統(tǒng)計性和經(jīng)濟性指標(biāo)。利用遺傳算法確定了利潤所失最小化時控制圖各參數(shù)的取值。最后，結(jié)合全因子試驗設(shè)計，分析了模型參數(shù)等對EL值的影響效應(yīng)。通過分析發(fā)現(xiàn)，各參數(shù)對利潤所失均具有顯著影響。受控時利潤與其正相關(guān)，失控時利潤和偏移大小均與其負(fù)相關(guān)。因此，在實際生產(chǎn)中，需要根據(jù)需求選擇合適的參數(shù)值。

表6 各因子的方差分析

表7 模型概述

[1]劉艷永,張公緒,孫靜.基于最大、最小特征根的多元協(xié)方差控制圖[J].方法與應(yīng)用,1998,(6).

[2]趙永滿.基于威沙特分布的二元過程均值與協(xié)方差監(jiān)控[J].統(tǒng)計與決策,2015,(5).

[3]Chou C Y,Chen C H,Chen C H.Economic Design of Variable Sampling Intervals T2Control Charts Using Genetic Algorithms[J].Expert Systems with Applications,2006,(30).

[4]Robert L,Mason,Nola D.Decomposition of T2for Multivariate Control Chart Interpretation[J].Journal of Quality Technology,1995,(27).

[5]Hayter A J,Tsui K L.Identification and Quantification in Multivariate Quality Control[J].Journal of Quality Technology,1994,(26).

[6]Maysa,Costa.Economic-Statistical Control Chart Design:A Sensitivity Study[J].Brazilian Journal of Operations&Production Management,2005,10(2).