逐步增加的Ⅱ型截尾壽命試驗下BS分布的統計分析

佘紀濤，徐曉嶺，顧蓓青

(上海對外經貿大學 統計與信息學院， 上海 201620)

【機械制造與檢測技術】

逐步增加的Ⅱ型截尾壽命試驗下BS分布的統計分析

佘紀濤，徐曉嶺，顧蓓青

(上海對外經貿大學 統計與信息學院， 上海 201620)

在逐步增加的Ⅱ型截尾壽命試驗場合下，研究兩參數Birnbaum-Saunders(簡稱BS)疲勞壽命分布產品的可靠性統計分析似然函數，運用泰勒展開得到了近似似然方程，通過求解僅含刻度參數的超越方程得到形狀參數的近似極大似然估計；最后通過Monte-Carlo模擬算例說明方法的可行性。

兩參數Birnbaum-Saunders疲勞壽命分布，逐步增加II型截尾，近似極大似然估計

在可靠性試驗中，逐步增加的Ⅱ型截尾壽命試驗在實際中經常碰到，具有理論與應用價值。譬如當組織一批價值昂貴的產品做壽命試驗，若已發生一部分產品失效時，考慮到試驗費用等因素，經常從未失效的產品中抽取一部分產品移離試驗現場，這樣可以大大節約成本。又例如在壽命試驗過程中為了解產品的失效機理和退化情況，以便更好設計產品，需要從尚未失效的產品中選取部分產品進行解剖分析。關于逐步增加的II型截尾壽命試驗的統計分析可見文獻[1-6]，其主要針對的是威布爾與對數正態等壽命分布。

Birnbaum-Saunders(簡稱BS)模型是概率物理方法的一個重要失效分布模型，這個模型是Birnbaum和Sauders在1969年研究主因裂紋擴展導致的材料失效過程中推導出來的，這一模型在機械產品可靠性研究中應用廣泛，主要應用于疲勞失效研究，對于電子產品性能退化失效分析也有重要應用。Birnbaum-Saunders疲勞壽命分布比常用壽命分布如威布爾分布,對數正態分布更適合描述某些由于疲勞而引起失效的產品壽命規律，因此該分布已經成為可靠性統計分析中的常用分布之一。關于Birnbaum-Saunders疲勞壽命的統計分析研究可見文獻[7-18]，其中關于統計分析的研究大多是全樣本場合，而針對Birnbaum-Saunder疲勞壽命分布逐步增加的Ⅱ型截尾壽命試驗的統計分析至今未見文獻報導。

逐步增加的Ⅱ型截尾試驗是Ⅱ型截尾試驗的一個推廣，它是在n個產品的樣本中先觀察到r1個失效，然后剩下的n-r1個未失效的產品中有n1個被移離試驗，留下n-r1-n1個產品繼續試驗；當再有r2-r1個產品失效后,將剩余仍未失效產品中的n2個被移離試驗。試驗如此進行下去，直到有一定數量的產品失效為止，終止試驗。試驗數據如表1所示，其中r0=0,n0=0,x(0)=0,ri≥1,ni≥0,i=1,2,…,k,k≥1，而nk為最后試驗結束時未失效的產品數。當ni=0,i=1,2,…,k-1時，便為通常的定數截尾壽命試驗；當r1=1,r2=2,…,rk=k時，通常稱為逐次截尾實驗。

表1 逐步增加的Ⅱ型截尾壽命試驗數據

本文給出兩參數Birnbaum-Saunders疲勞壽命分布在逐步增加Ⅱ型截尾壽命試驗下的近似極大似然估計，并通過實例分析說明該方法的應用可行。

1 兩參數Birnbaum-Saunders疲勞壽命分布參數的近似極大似然估計

將n個壽命服從兩參數Birnbaum-Saunders疲勞壽命分布(記為BS(α,β))的產品做逐步增加的Ⅱ截尾壽命試驗，試驗數據如表1所示。而兩參數BS(α,β)分布函數FX(x)與密度函數fX(x)分別為

x>0,α,β>0

易見

FX(x)=Φ(z)

而

似然函數：

注意到

(1)

(2)

(3)

將式(3)代入式(1)得：

化簡得：

即

α2-H(β)α-Q(β)≈0

(4)

將式(3)代入式(2)得：

于是得如下方程：

α2-G(β)α-L(β)≈0

(5)

由此可得：

(6)

進而可得僅含參數β的超越方程：

[Q(β)-L(β)]2-H(β)[Q(β)-L(β)][G(β)-H(β)]-

Q(β)[G(β)-H(β)]2≈0

(7)

(8)

再者

其中，

2 實例分析

例[19]：某種軸承的疲勞壽命，共10個數據(單位：小時)：

152.7， 172.0， 172.5， 173.5， 193.0， 204.7，

216.5， 234.9， 262.6， 422.6

先保留前5個數據，后按照隨機原則刪掉2個數據(刪掉234.9和262.6)后，留下3個數據(204.7，216.5，422.6)，利用本文方法可以得到參數的近似極大似然估計為

作為比較，文獻[19-20]利用全樣本數據得到如下6種點估計，見表2所示。

表2 全樣本場合下參數β,α的6種點估計

3 結論

介紹了在逐步增加的II型截尾壽命試驗條件下BS疲勞壽命分布產品的可靠性統計分析似然函數。及得到近似似然方程和形狀參數的近似極大似然估計的推導過程，通過實際算例說明了方法實用可行。

[1] 徐曉嶺，王蓉華.Weibull分布逐步增加的Ⅱ型截尾試驗的統計分析[J].強度與環境，2003，30(2)：31-37.

[2] 王炳興.Weibull分布基于定數逐次截尾壽命數據的統計分析[J].科技通報，2004，20(6)：488-496.

[3] 李鳳，師義民，田亞愛.逐步增加II型截尾下Weibull分布的Bayes估計[J].工程數學學報，2008，25(4)：641-650.

[4] 衛超，師義民.逐步Ⅱ型混合截尾下Pareto分布的統計分析[J].火力與指揮控制，2014，39(4)：27-33.

[5] 楊君慧，師義民，曹弘毅.逐步增加Ⅱ型截尾試驗下廣義指數分布的統計分析[J].統計與決策，2014(16)：28-30.

[6] BALAKRISHNAN N A.Clifford Cohen.Order Statistics and Inference Estimation Methods[M].Boston:Academic Press,1991.

[7] KUNDU D,KANNAN N,BALAKRISHNAN N.On the hazard function of birnbaum-saunders distribution and associated inference[J].Computational Statistics & Data Analysis,2008,52:2692-2702.

[8] WANG Ronghua,FEI Heliang.Statistical analysis for the birnbaum-saunders fatigue life distribution under multiply type II censoring[J].Chin Quart of Math,2006,21(1):15-27.

[9] WANG Ronghua,FEI Heliang.Statistical analysis for the birnbaum-saunders fatigue life distribution under type II bilateral censoring and multiply type II censoring,chin.quart[J].J of Math,2004,19(2):126-132.

[10] 孫祝嶺.Birnbaum-Saunders 疲勞壽命分布尺度參數的區間估計[J].兵工學報，2009，30(11)：1558-1561.

[11] 孫祝嶺.Birnbaum-Saunders 疲勞壽命分布參數的回歸估計方法[J].兵工學報，2010，31(9)：1260-1262.

[12] 孫祝嶺.疲勞壽命分布變異系數的統計推斷[J].質量與可靠性，2013(1)：13-15.

[13] 孫祝嶺.Birnbaum-Saunders分布環境因子的置信限[J].強度與環境，2012，39(4)：51-55.

[14] WANG Bingxing.Generalized interval estimation for the birnbaum-saunders distribution[J].Computational Statistics and Data Analysis,2012,56:4320-4326.

[15] NIU Cuizhen,GUO Xu,XU Wangli,et al.Comparison of several birnbaum-saunders distributions[J].Journal of Statistical Computation and Simulation,2014,84(12):2721-2733.

[16] 趙建印,孫權,彭寶華,等.基于加速退化數據的BS分布的統計推斷，電子產品可靠性與環境試驗，2006(1)：11-14.

[17] BALAKRISHNAN N,ZHU Xiaojun Zhu,On the existence and uniqueness of the maximum likelihood estimates of the parameters of birnbaum-saunders distribution based on type-I,type-II and hybrid censored samples[J].Statistics,2014，48(5)：1013-1032.

[18] 徐曉嶺,王蓉華,顧蓓青.關于兩參數Birnbaum-Saunders疲勞壽命分布統計分析的2個注記[J].浙江大學學報(理學版)，2016，43(5)：539-544.

[19] NG H K T,KUNDU D,BALAKRISHNAN N.Modified moment estimation for the two-parameter birnbaum-saunders distribution[J].Computational Statistics and Data Analysis,2003,43:283-298.

[20] SYED E.Ahmed,kamon budsaba,supranee lisawadi,andvei volodin,parametric estimation for the birnbaum-saunders lifetime distribution based on a new parametrization[J].Thailand Statistician July,2008,6(2):213-240.

(責任編輯唐定國)

StatisticalAnalysisofBSDistributionUnderProgressiveType-ⅡCensoredLifeTest

SHE Jitao, XU Xiaoling, GU Beiqing

(School of Statistics and Information, Shanghai University of International Business and Economics, Shanghai 201620, China)

The reliability statistical analysis method of two parameters Birnbaum-Saunders (BS) fatigue life distribution products is studied under progressive type-II censored life test. The likelihood function is given, and the approximate likelihood equation is obtained by Taylor expansion. By solving the transcendental equation with only scale parameter, we get approximate maximum likelihood estimation of shape parameter. Finally, the examples illustrated the feasibility of the method by Monte-Carlo simulations.

two-parameter Birnbaum-Saunders fatigue life distribution; incremental type II truncated censoring; approximate maximum likelihood estimation

2017-05-12；

：2017-05-30

：2017年上海市大學生創新創業訓練計劃項目(102732017093)；國家自然科學基金資助項目(11671264)

佘紀濤 (1995—)，女，主要從事應用統計研究。

10.11809/scbgxb2017.09.033

format:SHE Jitao, XU Xiaoling, GU Beiqing.Statistical Analysis of BS Distribution Under Progressive Type-Ⅱ Censored Life Test[J].Journal of Ordnance Equipment Engineering,2017(9):156-160.

O251

：A

2096-2304(2017)09-0156-05

本文引用格式：佘紀濤，徐曉嶺，顧蓓青.逐步增加的Ⅱ型截尾壽命試驗下BS分布的統計分析[J].兵器裝備工程學報，2017(9):156-160.