王雪崢

(92493 部隊 98 分隊 葫蘆島 125000)

引 言

在裝備試驗鑒定中，經(jīng)常會涉及到服從二項分布的指標參數(shù)，例如彈炮命中概率的鑒定問題、復(fù)雜系統(tǒng)可靠性鑒定問題等。關(guān)于二項分布參數(shù)的檢驗方法，常用的有經(jīng)典假設(shè)檢驗［1］、貝葉斯假設(shè)檢驗［2］以及序貫檢驗［3］。不管是哪一種假設(shè)檢驗方法，都離不開統(tǒng)計假設(shè)。以導(dǎo)彈命中概率p(0＜p＜1)為例，給出命中概率指標為某一個數(shù)值p0和最低可接受值p1，則相應(yīng)的簡單假設(shè)為:H0:p=p0，H1:p=p1。試驗方案應(yīng)具有如下特性:當實際p≥p0時，應(yīng)以不低于(1-α)的概率通過驗證試驗;當p≤p1時，應(yīng)以不高于β的概率通過驗證試驗。α為生產(chǎn)方風險，β為使用方風險。在該假設(shè)條件下，不管是采用經(jīng)典假設(shè)檢驗、貝葉斯假設(shè)檢驗，還是序貫檢驗，最終得出的結(jié)論是接受H0:p=p0，認為實際導(dǎo)彈命中概率p≥p0;或者拒絕H0:p=p0，而接受H1:p=p1，認為實際導(dǎo)彈命中概率p≤p1。當導(dǎo)彈實際命中概率p處在［p1，p0］區(qū)間上時，這種簡單假設(shè)方案無法做出判斷，只會判給p≥p0或p≤p1，這對于子樣數(shù)較少的導(dǎo)彈類武器試驗來說，無疑減少了試驗信息，甚至是提供了錯誤的試驗信息。本文提出的多假設(shè)序貫驗后加權(quán)檢驗MSPOT(Multi-hypothesis Sequential Posterior Odd Test)方法很好地解決了上述問題。

所謂多假設(shè)，是指統(tǒng)計假設(shè)打破傳統(tǒng)兩個假設(shè)(即原假設(shè)和備擇假設(shè))的情況，根據(jù)實際需求，將參數(shù)空間劃分成多個，對應(yīng)出現(xiàn)多個統(tǒng)計假設(shè)。例如前文導(dǎo)彈命中概率問題，參數(shù)空間劃分為0＜p＜p1，p1≤p＜p0，p0≤p＜1三個區(qū)域，則對應(yīng)三個假設(shè)為H0:0＜p＜p1，H1:p1≤p＜p0，H2:p0≤p＜1。這里不再分原假設(shè)和備擇假設(shè)，而是根據(jù)具體參數(shù)空間的意義來判斷，比如H0為真時，說明0＜p＜p1，即命中概率沒有達到最低可接受值;當H1為真時，說明p1≤p＜p0，即命中概率達到最低可接受值，但尚未達到設(shè)計指標要求;當H2為真時，說明p0≤p＜1，即命中概率達到設(shè)計指標要求。

序貫驗后加權(quán)檢驗SPOT(Sequential Posterior Odd Test)方法通常是將A.Wald提出的序貫概率比檢驗SPRT(Sequential Probability Ratio Test)方法的似然比換作似然函數(shù)在Θ0、Θ1上的驗后加權(quán)比，Θ0、Θ1是原假設(shè)和備擇假設(shè)所對應(yīng)的參數(shù)空間。SPOT方法是一種在現(xiàn)場試驗中根據(jù)試驗樣本的實際情況進行實時統(tǒng)計決策的試驗方法，這種方法一方面可以充分利用驗前信息，另一方面又可視試驗參數(shù)樣本的實際情況決定是否停止試驗新的樣本，它在工程上，特別是小子樣問題的解決上，具有重要的應(yīng)用價值［4～6］。

MSPOT方法是在多假設(shè)序貫概率比檢驗MSPRT(Multi-hypothesis Sequential Probability Ratio Test)［7，8］方法基礎(chǔ)上，將多簡單假設(shè)推廣到更具一般性的多參數(shù)空間假設(shè)，即Θ0，Θ1，...，ΘM-1，M為假設(shè)的個數(shù)，Θi(i=0，...，M-1)可以是單點值，也可以是連續(xù)的參數(shù)空間。當為單點值時，便是MSPRT方法，即多簡單假設(shè)的MSPOT方法。當為連續(xù)的參數(shù)空間時，便是多復(fù)雜假設(shè)的MSPOT方法?？梢姡琈SPRT方法是MSPOT方法的一種特例，MSPOT方法更具一般性。

圖1給出了MSPOT方法與MSPRT、SPRT、SPOT方法的關(guān)系。由圖可見其他幾種方法都是MSPOT方法的特例，MSPOT方法的優(yōu)越性在于其既利用了驗前信息，又解決了多參數(shù)空間假設(shè)檢驗的問題，包括多簡單假設(shè)和多復(fù)雜假設(shè)，而MSPRT方法只考慮了多假設(shè)中的簡單假設(shè)問題，SPOT方法雖然包括簡單假設(shè)和連續(xù)參數(shù)空間假設(shè)，但只考慮了兩假設(shè)。由此可見，MSPRT、SPRT和SPOT都具有自身的局限性。

圖1 MSPOT方法與MSPRT、SPRT、SPOT方法的關(guān)系Fig.1 Relationship of the MSPOT，MSPRT，SPRT and SPOT methods

本文基于多簡單假設(shè)SPOT方法［9］和多參數(shù)空間復(fù)雜假設(shè)SPOT方法［10］，研究提出了MSPOT方法，并針對常見的成敗型參數(shù)試驗方案設(shè)計問題，具體推導(dǎo)了二項分布參數(shù)下MSPOT方法的計算公式，對解決成敗型參數(shù)檢驗問題具有一定的指導(dǎo)意義。

1 MSPOT方法

1.1 MSPOT方法的序貫決策過程

設(shè)X1，X2，…是獨立同分布隨機變量，概率密度為f。如果密度函數(shù)f的參數(shù)為θ，設(shè)有參數(shù)空間Θ，Θ=Θ0∪Θ1∪…∪ΘM-1，則多假設(shè)為

如果各參數(shù)空間為單點集，則假設(shè)為多簡單假設(shè);如果參數(shù)空間為連續(xù)參數(shù)空間集，則假設(shè)為多復(fù)雜假設(shè)。

已知θ的先驗分布為π(θ)。根據(jù)MSPRT方法的思想［7］，真實假設(shè)具有相對較大的后驗概率，則有如下的序貫決策規(guī)則。

①序貫停止規(guī)則

序貫停止時的樣本數(shù)N取第一個大于0且使得對于至少一個k成立的n值。Ak為閾值，為后驗概率。

當假設(shè)為多簡單假設(shè)時有

當假設(shè)為多復(fù)雜假設(shè)時有

其中fk(X|θ)表示概率密度函數(shù)f在假設(shè)Hk成立參數(shù)空間為Θk時的概率密度形式。

②序貫決策規(guī)則

上述序貫決策規(guī)則中，公式形式比較復(fù)雜，因此，將式(1)、式(2)代入序貫停止規(guī)則不等式進行化簡整理。序貫決策過程轉(zhuǎn)化成:N取第一個大于0且使得Onk＜Ak至少對一個k成立的n值，如果Onk=min{Onj:j=0，···，M-1}，則Hk為真實假設(shè)。

當假設(shè)為多簡單假設(shè)時有

當假設(shè)為多復(fù)雜假設(shè)時有

1.2 MSPOT方法中閾值A(chǔ)k的計算過程

對于上節(jié)的序貫決策過程，當樣本分布已知、先驗概率密度已知時，Onk的計算容易得出。至于閾值A(chǔ)k(k=0，1，…，M-1)的計算，雖然多簡單假設(shè)對應(yīng)的MSPRT方法中有明確的數(shù)值表達式［7，8］，但計算起來很復(fù)雜，而多復(fù)雜假設(shè)的MSPOT方法中Ak的計算還很難找到確切的數(shù)值表達式。因此，本文采用仿真搜索法計算Ak。從序貫停止規(guī)則可以看出，Ak的大小直接影響停止時的樣本數(shù)，進而影響決策的準確度。經(jīng)分析，Ak越小，相應(yīng)的N越大，決策結(jié)果越準確。因此Ak值并不是唯一確定的，在不影響決策判斷的前提下，為了方便計算，這里不妨設(shè)Ak均相等，采用仿真搜索方法進行求解，搜索策略如下:

設(shè)搜索步長為ε，這里取ε=0.05。令初值A(chǔ)k=A=1(k=0，1，…，M-1)，誤判率允許值α0=0.3。所謂誤判率是指Hk為真實假設(shè)時，決策判斷Hk為假的概率，這里用誤判的頻率代替誤判概率，即誤判的次數(shù)/序貫試驗總次數(shù)。

①產(chǎn)生分布密度為π(θ)的隨機參數(shù)值θi(i=1，2，…，N)，N是事先給定的，這里令N=10。

②生成服從f(Xi|θ)的隨機樣本Xi(i=1，2，…，n)，計算

按照1.1節(jié)序貫決策過程進行序貫決策判斷。

③對步驟②操作重復(fù)進行K次，K是事先給定的，這里令K=100。計算出θi情況下序貫試驗的誤判率

④計算平均誤判率

通過上述過程便可求出待檢驗樣本的閾值A(chǔ)，即Ak(k=0，1，…，M-1)。

1.3 MSPOT方法截尾序貫決策過程

給定試驗樣本數(shù)N，當根據(jù)1.1節(jié)的序貫決策過程無法停止試驗時，強制停止試驗。截尾序貫決策規(guī)則為:如果則Hk為真實假設(shè)，其中ONj的計算公式同式(3)、式(4)。

2 二項分布參數(shù)檢驗的MSPOT方法

假設(shè)參數(shù)p的最低可接受值為p1，指標值為p0，則對應(yīng)的多復(fù)雜假設(shè)為

其參數(shù)空間為

根據(jù)式(4)計算得到多復(fù)雜假設(shè)情況下

若采用多簡單假設(shè)，則可用參數(shù)區(qū)間中點代表參數(shù)空間，即取

對應(yīng)的多簡單假設(shè)為

H0:p=μ0H1:p=μ1H2:p=μ2

根據(jù)式(3)計算得到多簡單假設(shè)情況下

當參數(shù)p無先驗信息時，各假設(shè)的先驗概率相同，則

3 仿真算例及分析

在導(dǎo)彈試驗過程中，對導(dǎo)彈的命中概率進行檢驗。當導(dǎo)彈在規(guī)定區(qū)域爆炸，認為導(dǎo)彈命中目標，記為1;否則認為失敗，記為0。可見，導(dǎo)彈的命中概率試驗是一個獨立的n重伯努利試驗，因此在n次試驗中導(dǎo)彈命中目標的次數(shù)X服從二項分布，即X～b(n，p)。假設(shè)某型導(dǎo)彈命中概率p1=0.75為最低可接受值，設(shè)計指標p0=0.9，且沒有任何驗前信息。當p＜p1時，認為該指標沒有達到最低可接受值;當p1≤p＜p0時，認為該指標在最低可接受值和設(shè)計指標之間;當p≥p0時，認為該指標達到設(shè)計指標要求。由于沒有先驗信息可以利用，采用無驗前信息的多簡單假設(shè)MSPOT方法。給定允許誤判率α0=0.3，根據(jù)式(9)計算出μ0、μ1、μ2值，再根據(jù)1.2節(jié)步驟①～⑤計算出A=0.7，此時的平均誤判率α*=0.279。下面借助計算機模擬，應(yīng)用MSPOT方法進行導(dǎo)彈命中概率仿真模擬試驗。Matlab［12］仿真模擬流程如圖2所示。

分別在三個參數(shù)空間取出一個值作為指標真值代表，例如取p=0.5＜p1，按照圖2流程進行序貫仿真試驗，所得結(jié)果如表1所示。當n=3時，有O30=0.1887＜A，停止試驗，此時O30=min{O30，O31，O32} ，判斷假設(shè)H0為真，這與p=0.5＜p1相符。p=0.8和p=0.95的仿真結(jié)果同樣證明了本文方法的正確性和有效性。

圖2 一次MSPOT試驗(p給定)模擬計算流程Fig.2 Flow chart of a MSPOT test calculation with given p

表1 多簡單假設(shè)MSPOT方法仿真試驗結(jié)果Table 1 The result of MSPOTmethod with multiple simple hypothesis

若事先已規(guī)定試驗次數(shù)，則采用截尾的MSPOT方法。例如，試驗次數(shù)為7次，如果有5次命中目標，記為(7，5)，則7次試驗中可能出現(xiàn)的試驗結(jié)果見表2。

表2 試驗結(jié)果Table 2 Test results

從表2中可以看出，前五種試驗結(jié)果，根據(jù)截尾的MSPOT方法判斷O70最小，H0是真實假設(shè)，說明參數(shù)p滿足p＜0.75;試驗結(jié)果為(7，5)和(7，6)時，判斷O71最小，H1為真實假設(shè)，說明參數(shù)p滿足0.75≤p＜0.9;試驗結(jié)果為(7，7)時，判斷O72最小，H2為真實假設(shè)，說明參數(shù)p滿足p≥0.9。

4 結(jié)束語

從本文的計算結(jié)果可以看出，應(yīng)用MSPOT方法進行試驗，試驗結(jié)果并不是簡單拒絕或接受，而是能夠提供試驗參數(shù)所在的區(qū)間。它有效解決了傳統(tǒng)簡單假設(shè)中，最低可接受值和設(shè)計指標值之間的模糊參數(shù)區(qū)域問題，為試驗雙方提供更多的試驗信息。本文例子沒有驗前信息，如果在實際試驗中能夠得到相對準確的驗前信息，則應(yīng)用MSPOT方法進行試驗，結(jié)果將會更貼近實際。

［1］閆章更，濮曉龍.現(xiàn)代軍事抽樣檢驗方法及應(yīng)用［M］.北京:國防工業(yè)出版社，2008.

［2］武小悅，劉 琦.裝備試驗與評價［M］.北京:國防工業(yè)出版社，2008.

［3］楊萬軍，張士峰，蔡 洪.序貫檢驗方案分析——二項分布情形［J］.飛行器測控學(xué)報，2002，21(2):61～63.

［4］王國玉，申緒澗，汪連棟，等.電子系統(tǒng)小子樣試驗理論方法［M］.北京:國防工業(yè)出版社，2003.

［5］張金槐，唐雪梅.Bayes方法(修訂版)［M］.長沙:國防科技大學(xué)出版社，1993.

［6］陳孝強，武小悅.Bayes序貫驗后加權(quán)檢驗方法的仿真分析［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2004，9(5):60～62.Chen Xiaoqiang， Wu Xiaoyue.Simulation Analysis of Sequential Posterior Odd Test Methods［J］.Tactical Missile Technology 2004，9(5):60～62.

［7］Baum W， Vecravalli V.A Sequential Procedure for Multihypothesis Testing［J］.IEEE Transactions on Information Theory，1994，40(6).

［8］Woodroof M.Nonlinear Renewal Theory in Sequential Analysis［M］.SLAM，1982.

［9］王雪崢，武小悅.多簡單假設(shè)的SPOT方法［J］.系統(tǒng)工程與電子技術(shù)，2007，29(9):1581～1584.Wang Xuezheng， Wu Xiaoyue.Multi-simple-hypothesis SPOT Method［J］.Systems Engineering and Electronics， 2007，29(9):1581～1584.

［10］王雪崢，韓成哲.可靠性試驗中多參數(shù)空間復(fù)雜假設(shè)的SPOT方法研究［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2010，(1):48～52.Wang Xuezheng，Han Chengzhe.Development of SPOT Method for Multiple Parameter Spaces with Complex Hypothesis in Reliability Test［J］.Tactical Missile Technology， 2010， (1):48 ～52.

［11］曲寶忠，孫曉峰.海軍戰(zhàn)術(shù)導(dǎo)彈試驗與鑒定［M］.北京:國防工業(yè)出版社，2005.

［12］張志涌，等.精通MATLAB 6.5版［M］.北京:北京航空航天大學(xué)出版社，2003.