小樣本條件下基于Bayes方法的裝備狀態(tài)評(píng)估

谷廣宇 劉建敏 喬新勇

摘要：??針對(duì)小樣本條件下單純采用使用期無(wú)法準(zhǔn)確衡量裝甲裝備技術(shù)狀況的問(wèn)題，本文以裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估為例，分析了小樣本條件下，裝備各技術(shù)狀況等級(jí)特征參數(shù)分布及其有效性和適用性。同時(shí)，采用Bayes決策法，從裝備屬于各技術(shù)狀況等級(jí)的概率和錯(cuò)判風(fēng)險(xiǎn)兩方面出發(fā)，對(duì)裝備狀態(tài)進(jìn)行評(píng)估。評(píng)估結(jié)果表明，利用兩種Bayes方法均能準(zhǔn)確評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況，評(píng)估準(zhǔn)確率達(dá)90%以上，并且基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法具有更好評(píng)估效果，評(píng)估準(zhǔn)確率較高，符合裝甲車輛技術(shù)狀況隨使用時(shí)間增加而逐漸變差的一般規(guī)律。該方法在裝備狀態(tài)評(píng)估中具有實(shí)用性及可行性。

關(guān)鍵詞：??狀態(tài)評(píng)估;?小樣本;?Bayes決策法;?發(fā)動(dòng)機(jī)

中圖分類號(hào)：?TK428;?E923.1?文獻(xiàn)標(biāo)識(shí)碼：?A

收稿日期：?2019-03-25;?修回日期：?2019-06-27

基金項(xiàng)目：??國(guó)家自然科學(xué)基金資助項(xiàng)目（U1836101）

作者簡(jiǎn)介：??谷廣宇（1989-），男，漢族，博士研究生，主要研究方向?yàn)榘l(fā)動(dòng)機(jī)狀態(tài)評(píng)估與預(yù)測(cè)。

裝甲裝備在使用過(guò)程中，其技術(shù)狀況將隨著使用時(shí)間的增長(zhǎng)發(fā)生變化，且呈現(xiàn)劣化趨勢(shì)。雖然其技術(shù)狀況在一定程度上可以采用運(yùn)行時(shí)間來(lái)衡量[1]，但由于使用環(huán)境、使用強(qiáng)度及工作地點(diǎn)等因素不盡相同，該方法有時(shí)不能完全真實(shí)地反映裝備的實(shí)際技術(shù)狀況。因此，需要運(yùn)用科學(xué)的方法對(duì)裝甲裝備進(jìn)行技術(shù)狀況評(píng)估[2]。通常情況下，在裝備的狀態(tài)檢測(cè)評(píng)估過(guò)程中，首先要采集裝備狀態(tài)信息，并對(duì)數(shù)據(jù)樣本進(jìn)行統(tǒng)計(jì)分析，才能得到相關(guān)結(jié)論。目前，基于數(shù)據(jù)驅(qū)動(dòng)的評(píng)估方法主要有模糊聚類[3]、支持向量機(jī)[4]、神經(jīng)網(wǎng)絡(luò)[5]、Bayes決策法[6]等。這些方法就統(tǒng)計(jì)分析而言，均需要盡可能大的樣本量，然而在實(shí)際評(píng)估過(guò)程中，受試驗(yàn)成本、試驗(yàn)周期等條件影響，無(wú)法進(jìn)行大量試驗(yàn)來(lái)獲取大樣本。在此情況下，仍采用傳統(tǒng)大樣本條件下的統(tǒng)計(jì)方法，難以保證其最終得到評(píng)估結(jié)果的可靠性。因此，研究在小樣本條件下的狀態(tài)檢測(cè)評(píng)估方法具有實(shí)際意義[7]。另外由于在狀態(tài)信息采集測(cè)量過(guò)程中，用以描述各個(gè)技術(shù)狀況的特征參數(shù)都是隨機(jī)變量，同一技術(shù)狀況下所有采集測(cè)量的樣本只可能是相似的，不可能完全一樣，而且某一特征值可能不只出現(xiàn)在某一特定狀況下，在其他狀況中也可能出現(xiàn)。鑒于Bayes決策理論既能參考各類狀態(tài)的分類錯(cuò)誤率，又能考慮分類錯(cuò)誤的風(fēng)險(xiǎn)程度，具有較強(qiáng)的分類決策能力。因此，本文采用Bayes決策法對(duì)裝備技術(shù)狀況進(jìn)行評(píng)估，并利用小樣本統(tǒng)計(jì)方法構(gòu)造各技術(shù)狀況下特征參數(shù)的分布，以獲取Bayes決策中的類條件概率密度。該研究可以對(duì)即將變“差”的車輛提前進(jìn)行預(yù)防性維修，降低故障風(fēng)險(xiǎn)。

1?小樣本統(tǒng)計(jì)原理與方法

為提高統(tǒng)計(jì)分析精度，在小樣本條件下，最有效的解決方法是設(shè)法增加信息量。在此類方法中，以美國(guó)斯坦福大學(xué)B.Efron教授提出的Bootstrap法[8]及B.R.Donald提出的Bayes?Bootstrap法[9]最為簡(jiǎn)單有效，而且在工程應(yīng)用中取得了豐碩成果[10-13]。其基本思路都是依據(jù)現(xiàn)有小樣本模仿未知分布，進(jìn)而獲取再生樣本，將小樣本問(wèn)題轉(zhuǎn)化成大樣本。因此，該方法適用于本文對(duì)不同技術(shù)狀況下小樣本狀態(tài)信息的統(tǒng)計(jì)分析。

1.1?Bootstrap方法

Bootstrap方法又稱自助估計(jì)法，其原理主要根據(jù)觀測(cè)到來(lái)自于未知總體分布F的隨機(jī)子樣X(jué)=（X1，…，Xn），估計(jì)總體分布F的某一分布特征R（X，F(xiàn)），如均值、方差等，從而推測(cè)總體分布F。具體方法如下：

設(shè)總體分布F的某個(gè)分布特征θ=θ（F）（如均值，方差等），由觀測(cè)子樣X(jué)=（X1，…，Xn）構(gòu)造經(jīng)驗(yàn)分布Fn，則有對(duì)θ的估計(jì)=（Fn），估計(jì)誤差為

Tn=（Fn）-（F）（1）

根據(jù)經(jīng)驗(yàn)分布Fn，重新抽取再生子樣X(jué)（1）=（X（1）1，…，X（1）n），進(jìn)而構(gòu)造經(jīng)驗(yàn)分布函數(shù)F（1）n。于是由X（1）又可得到θ的估計(jì)（F（1）n）。此時(shí)，可得到估計(jì)誤差Tn的Bootstrap統(tǒng)計(jì)量R（1）n，即

R（1）n=（F（1）n）-（Fn）（2）

重復(fù)抽取多組再生子樣X(jué)（i），i=1，2，…，m，可計(jì)算相應(yīng)R（i）n，i=1，2，…，m，并利用R（i）n的分布去逼近Tn的分布，并根據(jù)式（1）得到θ（F）的樣本，即

θ（i）（F）=（Fn）-Tn≌（Fn）-R（i）n（3）

該方法與傳統(tǒng)統(tǒng)計(jì)方法相比，在小樣本估計(jì)中具有較高精度。

1.2?Bayes?Bootstrap方法

Bayes?Bootstrap方法，即隨機(jī)加權(quán)自助法，也是將子樣信息“提攜”的方法，是對(duì)Bootstrap方法的改進(jìn)。其原理主要利用Dirichilet隨機(jī)向量對(duì)子樣X(jué)=（X1，…，Xn）進(jìn)行加權(quán)，從而計(jì)算隨機(jī)加權(quán)統(tǒng)計(jì)量Dn，進(jìn)而利用Dn模仿估計(jì)誤差Tn的分布。具體方法如下：

設(shè)子樣X(jué)=（X1，…，Xn）來(lái)自于未知總體分布F，μ和σ2為F的未知期望和方差，和S2為X的統(tǒng)計(jì)期望和方差，則估計(jì)誤差為

T1=-μ（4）

T2=nn-1S2-σ2（5）

式中，T1和T2為期望的估計(jì)誤差和方差的估計(jì)誤差。生成N個(gè)參數(shù)為（1，…，1）的Dirichilet隨機(jī)向量Vi=（vi1，vi2，…，vin），i=1，2，…，N，可獲取隨機(jī)加權(quán)統(tǒng)計(jì)量為

Di1=∑nk=1vikxk-（6）

Di2=nn-1∑nk=1vik（xk-）2-nn-1S2（7）

利用Di1和Di2模仿估計(jì)誤差T1和T2的分布，根據(jù)式（4）和式（5），可得

μ=-T1≈-1（8）

S2=nn-1S2-T1≈nn-1S2-2（9）

已有成果的實(shí)際應(yīng)用及仿真計(jì)算表明，該方法分析精度較高于Bootstrap方法[14]。

2?Bayes決策理論

Bayes決策理論是模式識(shí)別問(wèn)題中經(jīng)典的分類決策理論，其基本思想是通過(guò)各狀態(tài)下的先驗(yàn)概率分布推定后驗(yàn)概率[15]。根據(jù)決策規(guī)則，Bayes決策又分為基于最小錯(cuò)誤率[16]和基于最小風(fēng)險(xiǎn)的Bayes決策[17]。

2.1?基于最小錯(cuò)誤率的貝葉斯決策

假設(shè)共有m個(gè)決策狀態(tài)總體W=（w1，w2，…，wm），已知其先驗(yàn)概率P（wi），i=1，2，…，m，依據(jù)前文小樣本統(tǒng)計(jì)理論可推定各狀態(tài)下特征參數(shù)x的分布模型，并獲得其類條件概率密度P（x|wi）。則根據(jù)Bayes公式，可得各狀態(tài)后驗(yàn)概率為

P（wi|x）=P（x|wi）P（wi）∑mi=1P（x|wi）P（wi）（10）

并有決策規(guī)則：若P（wi|x）=max＼[P（wi|x）＼]，則x∈wi。此類決策規(guī)則稱為基于最小錯(cuò)誤率的貝葉斯決策。

2.2?基于最小風(fēng)險(xiǎn)的貝葉斯決策

實(shí)際決策過(guò)程中除了錯(cuò)誤率外，有時(shí)也需要考慮決策錯(cuò)誤所帶來(lái)的風(fēng)險(xiǎn)，在這種情況下，通?？梢赃x擇基于最小風(fēng)險(xiǎn)的貝葉斯決策方法[18]。

基于最小風(fēng)險(xiǎn)的貝葉斯方法在基于最小錯(cuò)誤方法上引入了“損失”的概念，即把屬于wi的狀態(tài)錯(cuò)判為aj狀態(tài)的損失，可記為cji，一般決策損失表如表1所示。

根據(jù)表1及前文所得后驗(yàn)概率，可計(jì)算決策為wi的狀態(tài)條件風(fēng)險(xiǎn)為

R（wi|x）=E＼[c（i|j）＼]=∑mj=1c（j|i）P（wj|x）（11）

并有決策規(guī)則：若R（wi|x）=min＼[R（wi|x）＼]，則x∈wi。

3?應(yīng)用實(shí)例

基于貝葉斯方法的裝備狀態(tài)評(píng)估，實(shí)質(zhì)上就是將裝備的技術(shù)狀況劃分為從優(yōu)到劣幾個(gè)等級(jí)的問(wèn)題，本文以裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估為例，并依據(jù)摩托小時(shí)（使用時(shí)間）將其技術(shù)狀況分為“優(yōu)”、“良”、“一般”、“較差”、“差”這5個(gè)等級(jí)，分別以“0～200”，“200～300”，“300～400”，“400～500”，“500～600”這5個(gè)階段統(tǒng)計(jì)特征表示，即有W=（w1，w2，w3，w4，w5）={“優(yōu)”，“良”，“一般”，“較差”，“差”}。

3.1?特征提取

評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況的參數(shù)很多，但實(shí)踐證明，通過(guò)減速時(shí)間評(píng)估技術(shù)狀況是一種簡(jiǎn)單有效的辦法[19]，因此本文選取柴油機(jī)減速時(shí)間作為技術(shù)狀況評(píng)估的特征參數(shù)。減速時(shí)間是指發(fā)動(dòng)機(jī)在最高空轉(zhuǎn)轉(zhuǎn)速下運(yùn)轉(zhuǎn)時(shí)迅速停止供油，使其自行減速，計(jì)算從最高空轉(zhuǎn)轉(zhuǎn)速減到轉(zhuǎn)速為0的時(shí)間。隨著使用時(shí)間增長(zhǎng)，發(fā)動(dòng)機(jī)活塞缸套磨損增加，阻力減小，減速時(shí)間將變長(zhǎng)[20]。

發(fā)動(dòng)機(jī)瞬時(shí)轉(zhuǎn)速如圖1所示，將實(shí)際計(jì)算瞬時(shí)轉(zhuǎn)速?gòu)膎3減速到n4的時(shí)間作為發(fā)動(dòng)機(jī)減速時(shí)間，此時(shí)取n3=2?200?r/min，n4=800?r/min。

依此方法在5種技術(shù)狀況等級(jí)下，分別對(duì)同一型號(hào)6輛不同摩托小時(shí)的裝甲車輛進(jìn)行測(cè)試，提取減速時(shí)間，并作為統(tǒng)計(jì)樣本，用以計(jì)算該等級(jí)下分布特征，以檢驗(yàn)本文評(píng)估方法的準(zhǔn)確性。另外，隨機(jī)抽取5輛未知狀況的裝甲車，提取修正減速時(shí)間作為測(cè)試樣本，利用本文方法對(duì)其進(jìn)行評(píng)估。統(tǒng)計(jì)樣本減速時(shí)間如表2所示，測(cè)試樣本減速時(shí)間如表3所示。

3.2?小樣本統(tǒng)計(jì)

假設(shè)發(fā)動(dòng)機(jī)減速時(shí)間在各個(gè)使用階段均服從正態(tài)分布，利用Bayes?Bootstrap小樣本統(tǒng)計(jì)方法，對(duì)各個(gè)階段減速時(shí)間進(jìn)行統(tǒng)計(jì)，根據(jù)式（4）～式（9），得各階段分布的均值和標(biāo)準(zhǔn)差，Bayes?Bootstrap法下各技術(shù)狀況分布特征如表4所示，各技術(shù)狀況等級(jí)分布特征如圖2所示。由表4和圖2可以看出，減速時(shí)間隨技術(shù)狀況等級(jí)變化較為明顯，因此能夠較好地反映發(fā)動(dòng)機(jī)技術(shù)狀況的變化。

3.3?基于Bayes決策的狀態(tài)評(píng)估

由于測(cè)試車輛摩托小時(shí)在0～600?h內(nèi)呈均勻分布，因此各狀態(tài)的先驗(yàn)概率為

P（W）=＼[1/3，1/6，1/6，1/6，1/6＼]

考慮評(píng)估過(guò)程中錯(cuò)誤率的存在，根據(jù)經(jīng)驗(yàn)，引入發(fā)動(dòng)機(jī)狀態(tài)評(píng)估決策損失C，發(fā)動(dòng)機(jī)狀態(tài)評(píng)估決策損失如表5所示。

3.4?評(píng)估準(zhǔn)確性對(duì)比

根據(jù)表4及式（10）和式（11），分別利用兩種Bayes方法對(duì)統(tǒng)計(jì)樣本進(jìn)行評(píng)估，統(tǒng)計(jì)樣本評(píng)估結(jié)果如表6所示。

由表6可以看出，利用兩種Bayes方法均能夠準(zhǔn)確評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況，評(píng)估準(zhǔn)確率達(dá)90%以上，并且基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法具有更好評(píng)估效果。

3.5?測(cè)試樣本狀態(tài)評(píng)估

采用基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法，對(duì)表3測(cè)試樣本進(jìn)行評(píng)估，在5種狀態(tài)下，各測(cè)試樣本的條件風(fēng)險(xiǎn)及評(píng)估結(jié)果如表7所示。

由表7可以看出，評(píng)估結(jié)果符合裝甲車輛技術(shù)狀況隨使用時(shí)間增加而逐漸變差的一般規(guī)律，證明該方法在裝備狀態(tài)評(píng)估中的實(shí)用性及可行性。

該方法還能判斷測(cè)試樣本的變化趨勢(shì)，例如測(cè)試樣本1根據(jù)表7被評(píng)估為“優(yōu)”狀態(tài)，但其“良”狀態(tài)的風(fēng)險(xiǎn)與“優(yōu)”狀態(tài)只有很小差別。因此，可判斷該樣本接近“優(yōu)”狀態(tài)與“良”狀態(tài)的邊緣，即將退化為“良”狀態(tài)。依據(jù)該方法可以對(duì)即將變“差”的車輛提前進(jìn)行預(yù)防性維修，以降低故障風(fēng)險(xiǎn)。

4?結(jié)束語(yǔ)

本文通過(guò)結(jié)合Bootstrap小子樣統(tǒng)計(jì)方法及Bayes決策法，有效地解決了小樣本條件下的裝甲裝備狀態(tài)評(píng)估的問(wèn)題，并以裝甲車輛發(fā)動(dòng)機(jī)為例，驗(yàn)證了該方法的有效性。采用Bootstrap小子樣統(tǒng)計(jì)方法，能夠有效擴(kuò)充發(fā)動(dòng)機(jī)狀態(tài)參數(shù)樣本，實(shí)現(xiàn)小樣本條件下的裝備狀態(tài)評(píng)估的先驗(yàn)統(tǒng)計(jì);Bayes決策法在裝備狀態(tài)評(píng)估過(guò)程中，不但考慮了各技術(shù)狀況出現(xiàn)的概率，還考慮了錯(cuò)誤評(píng)估的風(fēng)險(xiǎn)，比其他技術(shù)狀況評(píng)估方法更具合理性;采用本文方法進(jìn)行裝甲車輛發(fā)動(dòng)機(jī)的技術(shù)狀況等級(jí)評(píng)估，評(píng)估結(jié)果較符合實(shí)際規(guī)律，正確率可達(dá)90%以上。該研究具有一定工程實(shí)用性及可行性。

參考文獻(xiàn)：

[1]?李國(guó)錦.?基于層次分析及模糊綜合評(píng)判的裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估[J].?車輛與動(dòng)力技術(shù)，?2017（4）：?40-44.

[2]?張會(huì)奇，?陳春良，?劉峻巖，?等.?基于BP神經(jīng)網(wǎng)絡(luò)的裝甲車輛發(fā)動(dòng)機(jī)使用狀態(tài)評(píng)價(jià)[J].?兵工自動(dòng)化，?2014，?33（1）：?32-34.

[3]?辛江慧，?鄒政耀，?任成龍.??基于灰色定權(quán)聚類的混合動(dòng)力汽車發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)評(píng)估[J].??重慶理工大學(xué)學(xué)報(bào)：自然科學(xué)，?2014，?28（7）：?26-29.

[4]?汪瑾，?陳果，?王洪偉，?等.??基于后驗(yàn)概率SVM的航空發(fā)動(dòng)機(jī)滾動(dòng)軸承狀態(tài)評(píng)估[J].??軸承，?2015（5）：?55-59.

[5]?劉伯運(yùn)，?歐陽(yáng)光耀，?常漢寶.??基于神經(jīng)網(wǎng)絡(luò)和DS證據(jù)理論的柴油機(jī)狀態(tài)評(píng)估[J].??車用發(fā)動(dòng)機(jī)，?2005（5）：?14-17.

[6]?朱偉明.??基于貝葉斯融合的柴油機(jī)性能狀態(tài)評(píng)估方法研究[J].??內(nèi)燃機(jī)，?2010（2）：?18-19，?36.

[7]?高攀東，?陳曉陽(yáng)，?沈雪瑾，?等.??小樣本下軸承可靠性評(píng)估方法對(duì)比分析[J].??軸承，?2014（12）：??44-48.

[8]?Efron?B，?Tibshirani?R?J.??An?introduction?to?the?bootstrap[M].??London：?Chapman?and?Hall，?1993.

[9]?Donald?B?R.??The?bayesian?bootstrap[J].??The?Annals?of?Statistics，?1981，?9（1）：?130-134.

[10]?李會(huì)鵬，?明廷鋒，?賀國(guó).??小樣本條件下一種旋轉(zhuǎn)機(jī)械振動(dòng)狀態(tài)閾值的確定方法及應(yīng)用[J].??武漢理工大學(xué)學(xué)報(bào)：?交通科學(xué)與工程版，?2012，?36（5）：?1059-1063.

[11]?喬新勇，?劉建敏，?張小明.??小子樣統(tǒng)計(jì)理論在裝備狀態(tài)檢測(cè)中的應(yīng)用[J].??裝甲兵工程學(xué)院學(xué)報(bào)，?2009，?23（2）：?41-43，?87.

[12]?唐雪梅，?周伯昭，?李榮.??武器裝備小子樣綜合試驗(yàn)設(shè)計(jì)與鑒定技術(shù)[J].??戰(zhàn)術(shù)導(dǎo)彈技術(shù)，?2007（2）：?51-56，?79.

[13]?袁?鵬，?張衛(wèi)兵，?韓自剛.?銀川地區(qū)建筑地基沉降計(jì)算驗(yàn)系數(shù)取值[J].?科學(xué)技術(shù)與工程，?2017，?17（20）：??236-240.

[14]?樊鶴紅，?劉盼，?趙興群，?等.??小樣本條件下行波管可靠性評(píng)估方法的研究[J].??電子學(xué)報(bào)，?2010，?38（6）：?1394-1398.

[15]?張湘平，?張金槐，?謝紅衛(wèi)，?等.??基于隨機(jī)加權(quán)法的BAYES精度評(píng)定[J].??國(guó)防科技大學(xué)學(xué)報(bào)，?2001，?23（3）：?98-102.

[16]?全星慧，?牟海維，?呂秀麗，?等.??基于PCA與貝葉斯決策的人臉識(shí)別算法[J].?光學(xué)儀器，?2014，?36（2）：?122-125.

[17]?徐智康，?李旸，?李德玉.?基于可變最小貝葉斯風(fēng)險(xiǎn)的層次多標(biāo)簽分類方法[J].?南京大學(xué)學(xué)報(bào)：?自然科學(xué)版，?2017，?53（6）：?1023-1032.

[18]?毛思敏.??基于最小風(fēng)險(xiǎn)貝葉斯準(zhǔn)則的工業(yè)過(guò)程故障診斷[D].?武漢：華中科技大學(xué)，?2016.

[19]?張更云，?何來(lái)龍.??一種評(píng)價(jià)發(fā)動(dòng)機(jī)實(shí)際技術(shù)狀況的方法[J].??裝甲兵工程學(xué)院學(xué)報(bào)，?2005，?19（4）：?67-70.

[20]?張小明，?劉建敏，?喬新勇，?等.?柴油機(jī)健康狀況評(píng)估及其方法[J].?裝甲兵工程學(xué)院學(xué)報(bào)，?2008，?22（2）：?30-34.

Equipment?States?Evaluation?Based?on?Bayes?Method?with?Small?Sample

GU?Guangyu，?LIU?Jianmin，?QIAO?Xinyong

（Vehicle?Engineering?Department，?Army?Academy?of?Armored?Forces，?Beijing?100072，?China）

Abstract：??In?order?to?solve?the?problem?that?the?technical?condition?of?armored?equipment?can?not?be?accurately?measured?by?using?the?service?life?alone?under?the?condition?of?small?sample，?this?paper?takes?the?engine?state?evaluation?of?armored?vehicle?as?an?example，?and?analyzes?the?effectiveness?and?applicability?of?the?characteristic?parameters?distribution?in?each?technical?condition?of?the?equipment?under?the?condition?of?small?sample.??At?the?same?time，?the?Bayes?decision?method?is?used?to?evaluate?the?equipment?status?from?two?aspects：??the?probability?of?the?equipment?belonging?to?each?technical?condition?grade?and?the?misjudgment?risk.??The?evaluation?results?show?that?both?Bayes?methods?can?accurately?evaluate?the?engine?technical?status，?and?the?evaluation?accuracy?is?larger?than?90%.??The?Bayes?evaluation?method?based?on?minimum?risk?has?better?evaluation?effect?and?higher?evaluation?accuracy，?which?accords?with?the?general?law?that?the?technical?condition?of?armored?vehicles?becomes?worse?with?the?increase?of?service?time.??This?method?is?practical?and?feasible?in?equipment?condition?evaluation.

Key?words：??states?evaluation;?small?sample;?Bayes?decision;?engine